Continuous-variable (CV) quantum states are naturally characterized by their moments, defined as expectation values of products of single- or multimode ladder operators. Many CV Hamiltonians and quantum algorithms are formulated directly in terms of these moments, and therefore an efficient procedure to estimate moments with limited state measurements is necessary. In this paper, we present a protocol for shadow tomography of moment-generating functions (MGFs) of CV states based on randomized homodyne and heterodyne measurements. The resulting shadows enable an efficient and concurrent estimation of many multimode moments. Our complexity analysis shows that the number of measurements required to estimate a certain moment to a given precision grows exponentially in the order of the moments. Finally, we assess the precision of these moment estimators by applying them to two tasks: detecting entanglement in both Gaussian and non-Gaussian states via the Shchukin-Vogel protocol, and characterizing optical loss in a photonic chip. We demonstrate that both tasks can be accomplished with only a few thousand randomized measurements. The small number of required measurements, combined with efficient sample processing, makes our protocol applicable to a wide range of CV simulation tasks.
Continuous-variable (CV) quantum states are naturally characterized by their moments, defined as expectation values of products of single- or multimode ladder operators. Many CV Hamiltonians and quantum algorithms are formulated directly in terms of these moments, and therefore an efficient procedure to estimate moments with limited state measurements is necessary. In this paper, we present a protocol for shadow tomography of moment-generating functions (MGFs) of CV states based on randomized homodyne and heterodyne measurements. The resulting shadows enable an efficient and concurrent estimation of many multimode moments. Our complexity analysis shows that the number of measurements required to estimate a certain moment to a given precision grows exponentially in the order of the moments. Finally, we assess the precision of these moment estimators by applying them to two tasks: detecting entanglement in both Gaussian and non-Gaussian states via the Shchukin-Vogel protocol, and characterizing optical loss in a photonic chip. We demonstrate that both tasks can be accomplished with only a few thousand randomized measurements. The small number of required measurements, combined with efficient sample processing, makes our protocol applicable to a wide range of CV simulation tasks.
Shadow estimation provides an efficient framework for estimating observable expectation values using randomized measurements. While originally developed for discrete-variable systems, its recent extensions to continuous-variable (CV) quantum systems face practical limitations due to idealized assumptions of continuous phase modulation and infinite measurement resolution. In this work, we develop a practical shadow estimation protocol for CV systems using discretized homodyne detection with a finite number of phase settings and quadrature bins. We construct an unbiased estimator for the quantum state and establish both sufficient conditions and necessary conditions for informational completeness within a truncated Fock space up to nmax photons. We further provide a comprehensive variance analysis, showing that the shadow norm scales as O(nmax4), improving upon previous O(nmax13/3) bounds. Our work bridges the gap between theoretical shadow estimation and experimental implementations, enabling robust and scalable quantum state characterization in realistic CV systems.
Shadow tomography is a framework for constructing succinct descriptions of quantum states using randomized measurement bases, called “classical shadows,” with powerful methods to bound the estimators used. We recast existing experimental protocols for continuous-variable quantum state tomography in the classical-shadow framework, obtaining rigorous bounds on the number of independent measurements needed for estimating density matrices from these protocols. We analyze the efficiency of homodyne, heterodyne, photon-number-resolving, and photon-parity protocols. To reach a desired precision on the classical shadow of an N -photon density matrix with high probability, we show that homodyne detection requires order O(N4+1/3) measurements in the worst case, whereas photon-number-resolving and photon-parity detection require O(N4) measurements in the worst case (both up to logarithmic corrections). We benchmark these results against numerical simulation as well as experimental data from optical homodyne experiments. We find that numerical and experimental analyses of homodyne tomography match closely with our theoretical predictions. We extend our single-mode results to an efficient construction of multimode shadows based on local measurements. Published by the American Physical Society 2024
Non-Gaussian quantum states are critical resources in photonic quantum information processing, rendering their generation and characterization of increasing importance in quantum optics. In this work, we theoretically and numerically analyze the relative efficiency of homodyne versus heterodyne measurements for reconstructing non-Gaussian states, a major outstanding question in continuous-variable tomography. Combining a Fisher information-based formalism with simulated experiments, we find homodyne tomography to outperform heterodyne measurements for all non-Gaussian states tested, although the separation between the two modalities proves significantly narrower than suggested by the asymptotic Cramer-Rao lower bound. Our results should find use for optimizing measurement strategies in practical continuous-variable quantum systems.
In this article we develop a continuous variable (CV) shadow tomography scheme with wide ranging applications in quantum optics. Our work is motivated by the increasing experimental and technological relevance of CV systems in quantum information, quantum communication, quantum sensing, quantum simulations, quantum computing and error correction. We introduce two experimentally realisable schemes for obtaining classical shadows of CV (possibly non-Gaussian) quantum states using only randomised Gaussian unitaries and easily implementable Gaussian measurements such as homodyne and heterodyne detection. For both schemes, we show thatN=O(poly (1 /ϵ, log (1/δ),Mr+αn, log(m) )) samples of an unknownm-mode state ρ suffice to learn the expected value of anyr-local polynomial in the canonical observables of degree α, both with high probability 1 - δ and accuracy ϵ, as long as the state ρ has moments of ordern> α bounded byMn. By simultaneously truncating states and operators in energy and phase space, we are able to overcome new mathematical challenges that arise due to the infinite-dimensionality of CV systems. We also provide a scheme to learn nonlinear functionals of the state, such as entropies over any small number of modes, by leveraging recent energy-constrained entropic continuity bounds. Finally, we provide numerical evidence of the efficiency of our protocols in the case of CV states of relevance in quantum information theory, including ground states of quadratic Hamiltonians of many-body systems and cat qubit states. We expect our scheme to provide good recovery in learning relevant states of 2D materials and photonic crystals.
The entanglement of Gaussian continuous-variable (CV) states is fully determined by the state's second moments. In contrast, some entangled non-Gaussian states evade every second-moment criterion, and non-Gaussian entanglement detection remains an experimental challenge. The p3-PPT criterion detects entanglement using moments of the partial transpose of the density matrix. This criterion was recently extended to CV systems using photon-number-resolving detectors and multi-copy interferometry; here we introduce a single-copy homodyne protocol that detects bipartite CV entanglement via the same criterion. Unbiased U-statistic estimators for the partial-transpose moments p2 and p3 are constructed directly from randomized homodyne data and used to evaluate the p3-PPT entanglement witnesses: a linear one for detection, and a quadratic one whose violation yields a dimension-free lower bound on the entanglement negativity. The protocol estimates p2 and p3 up to additive error ε at Fock cutoff N from O((N+1)14/3/ε2) measurements at fixed confidence. We demonstrate the protocol on six families of Gaussian and non-Gaussian states, reaching 95% empirical one-sided detection probability from ∼103 to 104 homodyne measurements for states with nˉ≈2, placing non-Gaussian entanglement detection within reach of current homodyne experiments.
Cross-platform verification compares output states of independently operated bosonic and quantum optical systems when classical predictions and joint measurements are unavailable. Continuous-variable distributed inner product estimation (CV-DIPE) provides a local measurement approach, but infinite-dimensional state spaces require rigorous control of truncation error. We develop three CV-DIPE protocols for the natural class of energy-constrained continuous-variable quantum states. One uses coordinated displaced parity, another coordinated Haar measurements, and the third independent homodyne shadows. For arbitrary multimode inputs in this class, including entangled and non-Gaussian states, we establish truncation error bounds and worst-case sample complexity guarantees for the exact inner product. The leading terms of these bounds favor coordinated Haar measurements, whose global implementation is experimentally challenging. These protocols provide certified strategies for cross-platform verification of continuous-variable systems using only local measurement data.