Journal of Modern Optics

An inverse problem approach to optical homodyne tomography

Sze Meng Tan

Published Nov 1, 1997• 40 citations

Abstract

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.

DOI: 10.1080/09500349708231881

Related papers

Iterative maximum-likelihood reconstruction in quantum homodyne tomography

Alexander I. Lvovsky

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.

arXiv:quant-ph/0311097Journal of Optics B Quantum and Semiclassical OpticsMay 28, 2004304 citations

Tomographic reconstruction of the density matrix via pattern functions

Ulf Leonhardt, Henryk Paul, Giacomo Mauro D’Ariano

We propose a general method for reconstructing directly the density matrix of a single light mode in optical homodyne tomography. In our scheme the density matrix 〈a\ensuremath{\Vert}\ensuremath{\rho}^\ensuremath{\Vert}${\mathit{a}}^{\ensuremath{'}}$〉 is obtained by averaging a set of pattern functions ${\mathit{F}}_{\mathit{aa}}^{\ensuremath{'}}$(${\mathit{x}}_{\mathrm{\ensuremath{\theta}}}$,\ensuremath{\theta}) with respect to the homodyne data ${\mathit{x}}_{\mathrm{\ensuremath{\theta}}}$. The functions show the typical features of the quadrature distributions for the corresponding density-matrix elements. It is also possible to compensate the effect of detection losses which requires, however, extra effort in both experimental and numerical precision. We calculate the pattern functions for the coherent-state and Fock representations and study their properties. We believe that our method is the most efficient way for reconstructing the density matrix from homodyne measurements.

Physical Review ADec 1, 1995130 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.0526Physical Review ASep 6, 20119 citations

Quantum Tomography: Renormalization of Incompatible Observations

Z. Hradil, Robert Myška, Johann Summhammer +1 more

Standard deterministic techniques for quantum state reconstruction, as for example optical homodyne tomography, photon chopping, unbalanced homodyning etc., are based on the deterministic inversion of measured data. Since the frequencies obtained in realistic experiments always differ from probabilities predicted by quantum theory due to fluctuations, imperfections and realistic restrictions, the algorithm of inversion cannot guarantee the positive definiteness of the reconstructed density matrix. Hence the estimation of the noises may appear as doubtful. Quantum states may be successfully reconstructed within quantum and information theories using the maximum likelihood estimation. The question of deterministic schemes: “What quantum state is determined by that measurement?” is replaced by the formulation consistent with quantum theory: “What quantum state(s) seems to be most likely for that measurement?” Nonlinear equation for reconstructed state is formulated. An exact solution may be approached by subsequent iterations. Reconstruction is formulated as a problem of proper normalization of incompatible (nonorthogonal) measurements. The results obtained by this novel method may differ significantly from the standard predictions. Data are fitted better keeping the constraint of positive definiteness of reconstructed density matrix. However, this interpretation may enlarge uncertainty in prediction of quantum state in comparison with deterministic schemes, since, in general, there is a whole family of states which fit the measured data equally well. The novel technique is nonlinear and the reinterpretation of existing reconstruction schemes represents an advanced program. These keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

Kluwer Academic Publishers eBooksDec 19, 20050 citations

Detection of the density matrix through optical homodyne tomography without filtered back projection

Giacomo Mauro D’Ariano, Chiara Macchiavello, Matteo G. A. Paris

We propose an alternative method for detecting the density matrix of a radiation field via optical homodyne tomography. The method needs no assumption on the state, and the density matrix is recovered directly from averages on data. The tomographic reconstruction is very fast and provides very reliable statistics on line with detection. As a test of the method we present numericla results for a set of computer-simulated experiments, corresponding to different quantum states of the radiation source. Then a tomographic reconstruction is given for experimental data provided by Smithey etal. [Phys. Rev. Lett. 70, 1244 (1993)].

Physical Review ANov 1, 1994218 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, 19990 citations