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/0611244 · Physical Review A

Diluted maximum-likelihood algorithm for quantum tomography

Jaroslav Řeháček, Z. Hradil, Emanuel H. Knill, Alexander I. Lvovsky

Published Apr 18, 2007• 187 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.75.042108

Related papers

Maximum-Likelihood Quantum State Tomography by Cover's Method with Non-Asymptotic Analysis

Chien-Ming Lin, Hao–Chung Cheng, Yen-Huan Li

We propose an iterative algorithm that computes the maximum-likelihood estimate in quantum state tomography. The optimization error of the algorithm converges to zero at an $O ( ( 1 / k ) \log D )$ rate, where $k$ denotes the number of iterations and $D$ denotes the dimension of the quantum state. The per-iteration computational complexity of the algorithm is $O ( D ^ 3 + N D ^2 )$, where $N$ denotes the number of measurement outcomes. The algorithm can be considered as a parameter-free correction of the $R ρR$ method [A. I. Lvovsky. Iterative maximum-likelihood reconstruction in quantum homodyne tomography. \textit{J. Opt. B: Quantum Semiclass. Opt.} 2004] [G. Molina-Terriza et al. Triggered qutrits for quantum communication protocols. \textit{Phys. Rev. Lett.} 2004.].

arXiv:2110.00747•Oct 2, 2021•0 citations

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/0311097•Journal of Optics B Quantum and Semiclassical Optics•May 28, 2004•304 citations

Global convergence of diluted iterations in maximum-likelihood quantum tomography

Doglas S. Gonçalves, Márcia A. Gomes-Ruggiero, Carlile Lavor

In this paper we address convergence issues of the Diluted $R \rho R$ algorithm \cite{rehacek2007}, used to obtain the maximum likelihood estimate for the density matrix in quantum state tomography. We give a new interpretation to the diluted $R \rho R$ iterations that allows us to prove the global convergence under weaker assumptions. Thus, we propose a new algorithm which is globally convergent and suitable for practical implementation.

Quantum Information and Computation•Sep 1, 2014•4 citations

Global convergence of diluted iterations in maximum-likelihood quantum tomography

Gonçalves, D. S., Márcia A. Gomes-Ruggiero, C. Lavor

In this paper we present an inexact stepsize selection for the Diluted RρR algorithm, used to obtain the maximum likelihood estimate to the density matrix in quantum state tomography. We give a new interpretation for the diluted RρR iterations that allows us to prove the global convergence under weaker assumptions. Thus, we propose a new algorithm which is globally convergent and suitable for practical implementation.

arXiv:1306.3057•Jun 13, 2013•0 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

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
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 propose an iterative algorithm that computes the maximum-likelihood estimate in quantum state tomography. The optimization error of the algorithm converges to zero at an O((1/k)log⁡D)O ( ( 1 / k ) \log D )O((1/k)logD) rate, where kkk denotes the number of iterations and DDD denotes the dimension of the quantum state. The per-iteration computational complexity of the algorithm is O(D3+ND2)O ( D ^ 3 + N D ^2 )O(D3+ND2), where NNN denotes the number of measurement outcomes. The algorithm can be considered as a parameter-free correction of the RρRR ρRRρR method [A. I. Lvovsky. Iterative maximum-likelihood reconstruction in quantum homodyne tomography. \textit{J. Opt. B: Quantum Semiclass. Opt.} 2004] [G. Molina-Terriza et al. Triggered qutrits for quantum communication protocols. \textit{Phys. Rev. Lett.} 2004.].
In this paper we address convergence issues of the Diluted RρRR \rho RRρR algorithm \cite{rehacek2007}, used to obtain the maximum likelihood estimate for the density matrix in quantum state tomography. We give a new interpretation to the diluted RρRR \rho RRρR iterations that allows us to prove the global convergence under weaker assumptions. Thus, we propose a new algorithm which is globally convergent and suitable for practical implementation.