Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:quant-ph/0311097 · Journal of Optics B Quantum and Semiclassical Optics

Iterative maximum-likelihood reconstruction in quantum homodyne tomography

Alexander I. Lvovsky

Published May 28, 2004• 304 citations
arXivHTMLPDF

Abstract

I propose an iterative expectation maximization algorithm for reconstructing the density matrix of an optical ensemble from a set of balanced homodyne measurements. The algorithm applies directly to the acquired data, bypassing the intermediate step of calculating marginal distributions. The advantages of the new method are made manifest by comparing it with the traditional inverse Radon transformation technique.

DOI
10.1088/1464-4266/6/6/014

Related papers

Diluted maximum-likelihood algorithm for quantum tomography

Jaroslav Řeháček, Z. Hradil, Emanuel H. Knill +1 more

We propose a refined iterative likelihood-maximization algorithm for reconstructing a quantum state from a set of tomographic measurements. The algorithm is characterized by a very high convergence rate and features a simple adaptive procedure that ensures likelihood increase in every iteration and convergence to the maximum-likelihood state. We apply the algorithm to homodyne tomography of optical states and quantum tomography of entangled spin states of trapped ions and investigate its convergence properties.

arXiv:quant-ph/0611244•Physical Review A•Apr 18, 2007•187 citations

RECONSTRUCTION OF THE DENSITY MATRIX AS A CONSTRAINED OPTIMIZATION PROBLEM

Konrad Banaszek, Matteo G. A. Paris, Massimiliano Federico Sacchi +1 more

We present a numerical algorithm for the maximum-likelihood estimation of the density matrix, and apply it to the homodyne tomography of a single-mode radiation field. The algorithm is based on a specific form of the Gauss decomposition for positive definite Hermitian matrices. Results from Monte Carlo simulated experiments are presented. PACS: 42.50.Dv, 03.65.Bz Over past several years a great deal of interest has been paid to the reconstruction of the complete quantum state from experimental data [1]. Theoretical research in this field, apart from designing new measurement schemes, explores also the problem of extracting the optimal amount of information from data available in a realistic experimental setup. In this communication, we present a numerical algorithm for reconstructing the density matrix using the maximum-likelihood method [2, 3]. We discuss this algorithm in the context of quantum homodyne tomography [4, 5], which is actually the most successful method for quantum state reconstruction, and, in fact, the unique technique that has been so far experimentally implemented to measure nonclassical states of light [4, 6]. The realistic homodyne measurement of the quadrature is described by the positive operator-valued measure

•Jan 1, 1999•0 citations

An inverse problem approach to optical homodyne tomography

Sze Meng Tan

The problem of optical homodyne tomography is considered in the context of a Bayesian model-fitting or inverse problem approach. An algorithm is formulated, based on matrix computation rather than the numerical approximation of an analytic inverse transform. This automatically takes into account the effects of noise, detector inefficiencies and incomplete sampling of the data. The relationships with conventional reconstruction schemes, methods for including various forms of prior information and for calculating error estimates are discussed. The process of reconstructing the photon number distribution and the density matrix are illustrated using both simulated and experimental data.

Journal of Modern Optics•Nov 1, 1997•40 citations

Quantum homodyne tomography with a priori constraints

Konrad Banaszek

I present an iterative algorithm for reconstructing the Wigner function from homodyne statistics. The proposed method, based on maximum-likelihood estimation, is capable of compensating for detection losses in a numerically stable way.

arXiv:quant-ph/9901064•Physical Review A•Jun 1, 1999•18 citations

Optical homodyne tomography with polynomial series expansion

Hugo Benichi, Akira Furusawa

We present and demonstrate a method for optical homodyne tomography based on the inverse Radon transform. Different from the usual filtered back-projection algorithm, this method uses an appropriate polynomial series to expand the Wigner function and the marginal distribution, and discretize Fourier space. We show that this technique solves most technical difficulties encountered with kernel deconvolution-based methods and reconstructs overall better and smoother Wigner functions. We also give estimators of the reconstruction errors for both methods and show improvement in noise handling properties and resilience to statistical errors.

arXiv:1107.0526•Physical Review A•Sep 6, 2011•9 citations

Maximum-likelihood algorithm for quantum tomography

Konrad Banaszek

Optical homodyne tomography is discussed in the context of classical image processing. Analogies between these two fields are traced and used to formulate an iterative numerical algorithm for reconstructing the Wigner function from homodyne statistics.

arXiv:quant-ph/9904063•Apr 15, 1999•1 citation
I propose an iterative expectation maximization algorithm for reconstructing the density matrix of an optical ensemble from a set of balanced homodyne measurements. The algorithm applies directly to the acquired data, bypassing the intermediate step of calculating marginal distributions. The advantages of the new method are made manifest by comparing it with the traditional inverse Radon transformation technique.
We propose a refined iterative likelihood-maximization algorithm for reconstructing a quantum state from a set of tomographic measurements. The algorithm is characterized by a very high convergence rate and features a simple adaptive procedure that ensures likelihood increase in every iteration and convergence to the maximum-likelihood state. We apply the algorithm to homodyne tomography of optical states and quantum tomography of entangled spin states of trapped ions and investigate its convergence properties.
We present a numerical algorithm for the maximum-likelihood estimation of the density matrix, and apply it to the homodyne tomography of a single-mode radiation field. The algorithm is based on a specific form of the Gauss decomposition for positive definite Hermitian matrices. Results from Monte Carlo simulated experiments are presented. PACS: 42.50.Dv, 03.65.Bz Over past several years a great deal of interest has been paid to the reconstruction of the complete quantum state from experimental data [1]. Theoretical research in this field, apart from designing new measurement schemes, explores also the problem of extracting the optimal amount of information from data available in a realistic experimental setup. In this communication, we present a numerical algorithm for reconstructing the density matrix using the maximum-likelihood method [2, 3]. We discuss this algorithm in the context of quantum homodyne tomography [4, 5], which is actually the most successful method for quantum state reconstruction, and, in fact, the unique technique that has been so far experimentally implemented to measure nonclassical states of light [4, 6]. The realistic homodyne measurement of the quadrature is described by the positive operator-valued measure
The problem of optical homodyne tomography is considered in the context of a Bayesian model-fitting or inverse problem approach. An algorithm is formulated, based on matrix computation rather than the numerical approximation of an analytic inverse transform. This automatically takes into account the effects of noise, detector inefficiencies and incomplete sampling of the data. The relationships with conventional reconstruction schemes, methods for including various forms of prior information and for calculating error estimates are discussed. The process of reconstructing the photon number distribution and the density matrix are illustrated using both simulated and experimental data.
We present and demonstrate a method for optical homodyne tomography based on the inverse Radon transform. Different from the usual filtered back-projection algorithm, this method uses an appropriate polynomial series to expand the Wigner function and the marginal distribution, and discretize Fourier space. We show that this technique solves most technical difficulties encountered with kernel deconvolution-based methods and reconstructs overall better and smoother Wigner functions. We also give estimators of the reconstruction errors for both methods and show improvement in noise handling properties and resilience to statistical errors.