Direct estimation of arbitrary observables of an oscillator

Quantum harmonic oscillators serve as fundamental building blocks for quantum information processing, particularly in the context of the bosonic circuit quantum electrodynamics (cQED) platform. Conventional methods for extracting oscillator properties rely on predefined analytical gate sequences to access a restricted set of observables or resource-intensive tomography processes. Here, we introduce the Optimized Routine for Estimation of any Observable (OREO), a numerically optimized protocol that maps the expectation value of arbitrary oscillator observables onto that of an ancillary qubit. We demonstrate OREO in a bosonic cQED system as a means to efficiently measure phase-space quadratures and their higher moments, directly obtain faithful non-Gaussianity ranks, and effectively achieve state preparation independent of initial conditions in the oscillator. These results position OREO as a valuable tool for direct and efficient information extraction from bosonic quantum states, unlocking possibilities for measurement, control, and state preparation in continuous-variable quantum information processing.