© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2603.21380

Experimental Quantum State Tomography of Multimode Gaussian States

Chan Roh, Geunhee Gwak, Young-Do Yoon, Young-Sik Ra

quant-phSubmitted Mar 22, 2026• Announced Mar 24, 2026• 0 citations

v1 · Mar 22, 2026 v2 · Apr 10, 2026

Abstract

Multimode Gaussian states are a versatile resource for quantum information technologies and have been realized across a wide range of physical platforms. Recent progress in the large-scale generation of such states provides a key ingredient for scalable quantum technologies. Despite the importance of accurately characterizing these states, conventional tomography methods are often impractical because they require large sample sizes and can yield unphysical states. Here we present a reliable and efficient tomography method for multimode Gaussian states based on maximum-likelihood estimation. By directly operating on covariance matrices, the method avoids the exponential overhead associated with density-matrix reconstruction. We consider two commonly used detection schemes--single and joint homodyne detection--and systematically analyze the reconstruction performance. Our method outperforms conventional approaches by ensuring physical covariance matrices and achieving better agreement with the true states. To demonstrate the experimental applicability of the method, we experimentally generate various multipartite entangled states--six-mode graph states with different connectivity, a six-mode GHZ state, and a fully connected ten-mode graph state--and reconstruct their covariance matrices. Using the reconstructed covariance matrices, we quantify fidelities, detect entanglement, and reveal the multimode structure of squeezing and noise. Our technique offers a practical diagnostic tool for developing scalable quantum technologies.

Related papers

Energy-independent tomography of Gaussian states

Lennart Bittel, Mele, Francesco Anna, Jens Eisert +1 more

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.

arXiv:2508.14979•Aug 20, 2025•0 citations

Machine Learning Enhanced Quantum State Tomography: A Covariance Matrix Approach

J.C. Rodríguez, Hsien-Yi Hsieh, Hua Li Chen +2 more

We present a machine learning-based approach for efficient quantum state tomography, reconstructing covariance matrices of single- and two-mode squeezed states from homodyne data. This method ensures precise, real-time quantum property estimation and noise robustness.

•Jan 1, 2025•0 citations

Efficient convex optimization for bosonic state tomography

Shengyong Li, Yanjin Yue, Ying Hu +7 more

Quantum states encoded in electromagnetic fields, also known as bosonic states, have been widely applied in quantum sensing, quantum communication, and quantum error correction. Accurate characterization is therefore essential yet difficult when states cannot be reconstructed with sparse Pauli measurements. Tomography must work with dense measurement bases, high-dimensional Hilbert spaces, and often sample-based data. However, existing convex optimization-based techniques are not efficient enough and scale poorly when extended to large and multimode systems. In this work, we explore convex optimization as an effective framework to address problems in bosonic state tomography, introducing three techniques to enhance efficiency and scalability: efficient displacement operator computation, Hilbert-space truncation, and stochastic convex optimization, which mitigate common limitations of existing approaches. Then we propose a sample-based, convex maximum-likelihood estimation method specifically designed for flying-mode tomography. Numerical simulations of flying four-mode and nine-mode problems demonstrate the accuracy and practicality of our methods. This method provides practical tools for reliable bosonic mode quantum state reconstruction in high-dimensional and multimode systems.

arXiv:2509.06305•Physical Review Applied•Dec 4, 2025•0 citations

Machine learning for quantum state tomography: robust covariance matrix estimation for squeezed vacuum states with thermal noise

Juan Camilo Rodríguez, Hsien-Yi Hsieh, Hua Li Chen +3 more

We present a supervised machine learning-based method using convolutional neural networks to estimate the covariance matrix of Gaussian quantum states in the presence of thermal noise. Unlike computationally intensive density matrix reconstructions, our machine learning-based method allows for the reconstruction of impure squeezed vacuum states using sparse measurements of quadrature sequences based on a model employing a two-component state mixed together from thermal and squeezed thermal states. The method achieves high fidelity and precision, notably also at high squeezing levels, while offering an effective characterization of physical quantities and accurately estimating the covariance matrix. We benchmark our machine against experimental data of single-mode squeezed vacuum states, demonstrating its accuracy and capability to quantify experimental degradation to squeezing and purity. We experimentally verify that our covariance matrix estimation exhibits robustness to state degradation induced by thermal state admixtures. We provide a method for lightweight, compact, and complete representation of lab-generated Gaussian states and lay the foundation for extending real-time quantum state tomography for thermal multi-component Gaussian states to multi-mode systems.

Optics Express•Oct 31, 2025•0 citations

Efficient Tomography of Microwave Photonic Cluster States

Yoshiki Sunada, Shingo Kono, Jesper Ilves +6 more

Entanglement among a large number of qubits is a crucial resource for many quantum algorithms. Such many-body states have been efficiently generated by entangling a chain of itinerant photonic qubits in the optical or microwave domain. However, it has remained challenging to fully characterize the generated many-body states by experimentally reconstructing their exponentially large density matrices. Here, we develop an efficient tomography method based on the matrix-product-operator formalism and demonstrate it on a cluster state of up to 35 microwave photonic qubits by reconstructing its 2 35 × 2 35 density matrix. The full characterization enables us to detect the performance degradation of our photon source which occurs only when generating a large cluster state. This tomography method is generally applicable to various physical realizations of entangled qubits and provides an efficient benchmarking method for guiding the development of high-fidelity sources of entangled photons.

arXiv:2410.03345•PRX Quantum•Oct 10, 2025•2 citations

Wigner function tomography via optical parametric amplification

Mahmoud Kalash, Maria V. Chekhova

Wigner function tomography is indispensable for characterizing quantum states, but its commonly used version, balanced homodyne detection, suffers from several weaknesses. First, it requires efficient detection, which is critical for measuring fragile non-Gaussian states, especially bright ones. Second, it needs a local oscillator, tailored to match the spatiotemporal properties of the state under test, and fails for multimode and broadband states. Here we propose Wigner function tomography based on optical parametric amplification followed by direct detection. The method is immune to detection inefficiency and loss, and suitable for broadband, spatially and temporally multimode quantum states. To prove the principle, we experimentally reconstruct the Wigner function of squeezed vacuum occupying a single mode of a strongly multimode state. We obtain a squeezing of −7.5±0.4dB and purity of 0.91 −0.08 +0.09 despite more than 97% loss caused mainly by filtering. Theoretically, we also consider the reconstruction of a squeezed single photon—a bright non-Gaussian state. Due to multimode parametric amplification, the method allows for simultaneous tomography of multiple modes. This makes it a powerful tool for optical quantum information processing.

arXiv:2207.10030•Optica•Jul 26, 2023•33 citations

Multimode Gaussian states are a versatile resource for quantum information technologies and have been realized across a wide range of physical platforms. Recent progress in the large-scale generation of such states provides a key ingredient for scalable quantum technologies. Despite the importance of accurately characterizing these states, conventional tomography methods are often impractical because they require large sample sizes and can yield unphysical states. Here we present a reliable and efficient tomography method for multimode Gaussian states based on maximum-likelihood estimation. By directly operating on covariance matrices, the method avoids the exponential overhead associated with density-matrix reconstruction. We consider two commonly used detection schemes--single and joint homodyne detection--and systematically analyze the reconstruction performance. Our method outperforms conventional approaches by ensuring physical covariance matrices and achieving better agreement with the true states. To demonstrate the experimental applicability of the method, we experimentally generate various multipartite entangled states--six-mode graph states with different connectivity, a six-mode GHZ state, and a fully connected ten-mode graph state--and reconstruct their covariance matrices. Using the reconstructed covariance matrices, we quantify fidelities, detect entanglement, and reveal the multimode structure of squeezing and noise. Our technique offers a practical diagnostic tool for developing scalable quantum technologies.

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.

Quantum states encoded in electromagnetic fields, also known as bosonic states, have been widely applied in quantum sensing, quantum communication, and quantum error correction. Accurate characterization is therefore essential yet difficult when states cannot be reconstructed with sparse Pauli measurements. Tomography must work with dense measurement bases, high-dimensional Hilbert spaces, and often sample-based data. However, existing convex optimization-based techniques are not efficient enough and scale poorly when extended to large and multimode systems. In this work, we explore convex optimization as an effective framework to address problems in bosonic state tomography, introducing three techniques to enhance efficiency and scalability: efficient displacement operator computation, Hilbert-space truncation, and stochastic convex optimization, which mitigate common limitations of existing approaches. Then we propose a sample-based, convex maximum-likelihood estimation method specifically designed for flying-mode tomography. Numerical simulations of flying four-mode and nine-mode problems demonstrate the accuracy and practicality of our methods. This method provides practical tools for reliable bosonic mode quantum state reconstruction in high-dimensional and multimode systems.

We present a supervised machine learning-based method using convolutional neural networks to estimate the covariance matrix of Gaussian quantum states in the presence of thermal noise. Unlike computationally intensive density matrix reconstructions, our machine learning-based method allows for the reconstruction of impure squeezed vacuum states using sparse measurements of quadrature sequences based on a model employing a two-component state mixed together from thermal and squeezed thermal states. The method achieves high fidelity and precision, notably also at high squeezing levels, while offering an effective characterization of physical quantities and accurately estimating the covariance matrix. We benchmark our machine against experimental data of single-mode squeezed vacuum states, demonstrating its accuracy and capability to quantify experimental degradation to squeezing and purity. We experimentally verify that our covariance matrix estimation exhibits robustness to state degradation induced by thermal state admixtures. We provide a method for lightweight, compact, and complete representation of lab-generated Gaussian states and lay the foundation for extending real-time quantum state tomography for thermal multi-component Gaussian states to multi-mode systems.

Entanglement among a large number of qubits is a crucial resource for many quantum algorithms. Such many-body states have been efficiently generated by entangling a chain of itinerant photonic qubits in the optical or microwave domain. However, it has remained challenging to fully characterize the generated many-body states by experimentally reconstructing their exponentially large density matrices. Here, we develop an efficient tomography method based on the matrix-product-operator formalism and demonstrate it on a cluster state of up to 35 microwave photonic qubits by reconstructing its 2 35 × 2 35 density matrix. The full characterization enables us to detect the performance degradation of our photon source which occurs only when generating a large cluster state. This tomography method is generally applicable to various physical realizations of entangled qubits and provides an efficient benchmarking method for guiding the development of high-fidelity sources of entangled photons.

Wigner function tomography is indispensable for characterizing quantum states, but its commonly used version, balanced homodyne detection, suffers from several weaknesses. First, it requires efficient detection, which is critical for measuring fragile non-Gaussian states, especially bright ones. Second, it needs a local oscillator, tailored to match the spatiotemporal properties of the state under test, and fails for multimode and broadband states. Here we propose Wigner function tomography based on optical parametric amplification followed by direct detection. The method is immune to detection inefficiency and loss, and suitable for broadband, spatially and temporally multimode quantum states. To prove the principle, we experimentally reconstruct the Wigner function of squeezed vacuum occupying a single mode of a strongly multimode state. We obtain a squeezing of −7.5±0.4dB and purity of 0.91 −0.08 +0.09 despite more than 97% loss caused mainly by filtering. Theoretically, we also consider the reconstruction of a squeezed single photon—a bright non-Gaussian state. Due to multimode parametric amplification, the method allows for simultaneous tomography of multiple modes. This makes it a powerful tool for optical quantum information processing.