When learning descriptions of many-body quantum systems, assuming locality of interactions can dramatically reduce the required number of observations and computational time compared to unstructured approaches. In this work, we uncover an infinite-dimensional instance of this phenomenon, initiating the study of Hamiltonian learning for positive-temperature bosonic Gaussian states. We present sample- and computationally efficient protocols for inferring the parameters and graph structure of underlying quadratic Hamiltonians, provided that temperature, squeezing, displacement, and the interaction graph's maximal degree are bounded. Notably, our protocol relies solely on heterodyne measurements, which are often experimentally feasible, and achieves a sample complexity that scales only logarithmically with the number of modes. Furthermore, our techniques allow for the learning of positive-temperature Gaussian states in trace distance with quadratic scaling in precision and polynomial scaling in the number of modes.
When learning descriptions of many-body quantum systems, assuming locality of interactions can dramatically reduce the required number of observations and computational time compared to unstructured approaches. In this work, we uncover an infinite-dimensional instance of this phenomenon, initiating the study of Hamiltonian learning for positive-temperature bosonic Gaussian states. We present sample- and computationally efficient protocols for inferring the parameters and graph structure of underlying quadratic Hamiltonians, provided that temperature, squeezing, displacement, and the interaction graph's maximal degree are bounded. Notably, our protocol relies solely on heterodyne measurements, which are often experimentally feasible, and achieves a sample complexity that scales only logarithmically with the number of modes. Furthermore, our techniques allow for the learning of positive-temperature Gaussian states in trace distance with quadratic scaling in precision and polynomial scaling in the number of modes.
Continuous-variable systems enable key quantum technologies in computation, communication, and sensing. Bosonic Gaussian states emerge naturally in various such applications, including gravitational-wave and dark-matter detection. A fundamental question is how to characterize an unknown bosonic Gaussian state from as few samples as possible. Despite decades-long exploration, the ultimate efficiency limit remains unclear. In this work, we study the necessary and sufficient number of copies to learn an n-mode Gaussian state, with energy less than E, to ε trace distance with high probability. We prove a lower bound of Ω(n3/ε2) for Gaussian measurements, matching the best known upper bound up to doubly-log energy dependence, and Ω(n2/ε2) for arbitrary measurements. We further show an upper bound of O(n2/ε2) given that the Gaussian state is promised to be either pure or passive. Interestingly, while Gaussian measurements suffice for nearly optimal learning of pure Gaussian states, non-Gaussian measurements are provably required for optimal learning of passive Gaussian states. Finally, focusing on learning single-mode Gaussian states via non-entangling Gaussian measurements, we provide a nearly tight bound of Θ(E/ε2) for any non-adaptive schemes, showing adaptivity is indispensable for nearly energy-independent scaling. As a byproduct, we establish sharp bounds on the trace distance between Gaussian states in terms of the total variation distance between their Wigner distributions, and obtain a nearly tight sample complexity bound for learning the Wigner distribution of any Gaussian state to ε total variation distance. Our results significantly advance quantum learning theory in the bosonic regimes and have practical impact in quantum sensing and benchmarking applications.
The exploration of tomography of bosonic Gaussian states is presumably as old as quantum optics, but only recently, their precise and rigorous study have been moving into the focus of attention, motivated by technological developments. In this work, we present an efficient and experimentally feasible Gaussian state tomography algorithm with provable recovery trace-distance guarantees, whose sample complexity depends only on the number of modes, and - remarkably - is independent of the state's photon number or energy, up to doubly logarithmic factors. Our algorithm yields a doubly-exponential improvement over existing methods, and it employs operations that are readily accessible in experimental settings: the preparation of an auxiliary squeezed vacuum, passive Gaussian unitaries, and homodyne detection. At its core lies an adaptive strategy that systematically reduces the total squeezing of the system, enabling efficient tomography. Quite surprisingly, this proves that estimating a Gaussian state in trace distance is generally more efficient than directly estimating its covariance matrix. Our algorithm is particularly well-suited for applications in quantum metrology and sensing, where highly squeezed - and hence high-energy - states are commonly employed. As a further contribution, we establish improved sample complexity bounds for standard heterodyne tomography, equipping this widely used protocol with rigorous trace-norm guarantees.
We develop a protocol for learning a class of interacting bosonic Hamiltonians from dynamics with Heisenberg-limited scaling. For Hamiltonians with an underlying bounded-degree graph structure, we can learn all parameters with root mean square error ϵ using O(1/ϵ) total evolution time, which is independent of the system size, in a way that is robust against state-preparation and measurement error. In the protocol, we only use bosonic coherent states, beam splitters, phase shifters, and homodyne measurements, which are easy to implement on many experimental platforms. A key technique we develop is to apply random unitaries to enforce symmetry in the effective Hamiltonian, which may be of independent interest.
Learning the Hamiltonian underlying a quantum many-body system in thermal equilibrium is a fundamental task in quantum learning theory and experimental sciences. To learn the Gibbs state of local Hamiltonians at any constant inverse temperature, the state-of-the-art provable algorithms fall short of the optimal sample and computational complexity, in sharp contrast with the locality and simplicity in the classical cases. In this work, we present a learning algorithm that learns each local term of a n-qubit Hamiltonian on any bounded-degree graph to a constant additive error with the optimal sample complexity O(logn). The protocol uses parallelizable local quantum measurements that act within bounded neighborhoods of the graph and near-linear-time classical post-processing. We also give a learning algorithm for lattice Hamiltonians with near-optimal scaling on the learning precision and the inverse temperature. At the heart of our algorithm is the interplay between locality, the Kubo-MartinSchwinger condition, and the operator Fourier transform at arbitrary temperatures.
Bosonic Gaussian unitaries are fundamental building blocks of central continuous-variable quantum technologies such as quantum-optic interferometry and bosonic error-correction schemes. In this work, we present the first time-efficient algorithm for learning bosonic Gaussian unitaries with a rigorous analysis. Our algorithm produces an estimate of the unknown unitary that is accurate to small worst-case error, measured by the physically motivated energy-constrained diamond distance. Its runtime and query complexity scale polynomially with the number of modes, the inverse target accuracy, and natural energy parameters quantifying the allowed input energy and the unitary's output-energy growth. The protocol uses only experimentally friendly photonic resources: coherent and squeezed probes, passive linear optics, and heterodyne/homodyne detection. We then employ an efficient classical post-processing routine that leverages a symplectic regularization step to project matrix estimates onto the symplectic group. In the limit of unbounded input energy, our procedure attains arbitrarily high precision using only 2m+2 queries, where m is the number of modes. To our knowledge, this is the first provably efficient learning algorithm for a multiparameter family of continuous-variable unitaries.
A fundamental question in physics is: When does classical behavior emerge from quantum systems? Bosonic Gaussian states provide a natural setting to explore this quantum-classical boundary, as they capture both the classical field behavior and the intrinsic quantum nature of light. Here, we address this problem from a learning-theoretic perspective by asking: When are bosonic Gaussian states classical to learn? That is, under what conditions (if any) can an n-mode bosonic Gaussian state be learned with as few samples, and with operations as simple, as are needed to learn a classical 2n-variate Gaussian distribution? We establish a smooth crossover in learnability governed by the state's thermal fluctuations:
- Cold Gaussian states are non-classical to learn: When the covariance matrix satisfies Σ≤(21+O(n1))I, i.e. close to the vacuum covariance, tomography under single-copy (i.e., non-entangled) measurements fundamentally requires Ω(n3) copies, strictly exceeding the sample complexity Θ(n2) of learning classical Gaussian distributions. We show that this hardness persists even when few-copy entangled measurements are allowed.
- Warm Gaussian states are classical to learn: When thermal fluctuations exceed the vacuum noise, parameterized by Σ≥(21+ν)I for any parameter ν>0, we prove that single-copy tomography requires N=Θ(n2min(n,1+ν−1)) copies. This bound is tight and is achieved by simple, non-adaptive, unentangled heterodyne measurements. Crucially, for ν=Ω(1), the sample complexity drops to Θ(n2), matching the classical case.
Our results tightly characterize a quantum-to-classical crossover in the learnability of bosonic Gaussian states, reveal a novel connection between fundamental physics and statistical learning theory, and have implications for real-world sensing experiments.