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:1205.2300 · New Journal of Physics

Quantum tomography via compressed sensing: error bounds, sample complexity and efficient estimators

Steven T. Flammia, David Gross, Yi-Kai Liu, Jens Eisert

Published Sep 27, 2012• 356 citations
arXivHTMLPDF

Abstract

Intuitively, if a density operator has small rank, then it should be easier to estimate from experimental data, since in this case only a few eigenvectors need to be learned. We prove two complementary results that confirm this intuition. Firstly, we show that a low-rank density matrix can be estimated using fewer copies of the state, i.e. the sample complexity of tomography decreases with the rank. Secondly, we show that unknown low-rank states can be reconstructed from an incomplete set of measurements, using techniques from compressed sensing and matrix completion. These techniques use simple Pauli measurements, and their output can be certified without making any assumptions about the unknown state. In this paper, we present a new theoretical analysis of compressed tomography, based on the restricted isometry property for low-rank matrices. Using these tools, we obtain near-optimal error bounds for the realistic situation where the data contain noise due to finite statistics, and the density matrix is full-rank with decaying eigenvalues. We also obtain upper bounds on the sample complexity of compressed tomography, and almost-matching lower bounds on the sample complexity of any procedure using adaptive sequences of Pauli measurements. Using numerical simulations, we compare the performance of two compressed sensing estimators—the matrix Dantzig selector and the matrix Lasso—with standard maximum-likelihood estimation (MLE). We find that, given comparable experimental resources, the compressed sensing estimators consistently produce higher fidelity state reconstructions than MLE. In addition, the use of an incomplete set of measurements leads to faster classical processing with no loss of accuracy. Finally, we show how to certify the accuracy of a low-rank estimate using direct fidelity estimation, and describe a method for compressed quantum process tomography that works for processes with small Kraus rank and requires only Pauli eigenstate preparations and Pauli measurements.

DOI
10.1088/1367-2630/14/9/095022

Related papers

Research on Technology of Quantum Tomography Based on Compressed Sensing

Wenying Li

This paper establishes the methods for quantum state tomography based on compressed sensing.These methods are specialized for quantum states that are fairly pure,and they offer a significant performance improvement on large quantum systems.In particular,they are able to reconstruct an unknown density matrix of dimension d and rank r using O(rdlog2d) measurement settings,compared to standard methods that require d2 settings.The methods have several features that make them amenable to experimental implementation: they require only simple Pauli measurements,use fast convex optimization,are stable against noise,and can be applied to states that are only approximately low rank.The acquired data can be used to certify that the state is indeed close to pure,so no a priori assumptions are needed.

Computer and Modernization•Jan 1, 2011•0 citations

Quantum State Tomography via Compressed Sensing

David Gross, Yi-Kai Liu, Steven T. Flammia +2 more

We establish methods for quantum state tomography based on compressed sensing. These methods are specialized for quantum states that are fairly pure, and they offer a significant performance improvement on large quantum systems. In particular, they are able to reconstruct an unknown density matrix of dimension d and rank r using O(rdlog²d) measurement settings, compared to standard methods that require d² settings. Our methods have several features that make them amenable to experimental implementation: they require only simple Pauli measurements, use fast convex optimization, are stable against noise, and can be applied to states that are only approximately low rank. The acquired data can be used to certify that the state is indeed close to pure, so no a priori assumptions are needed.

arXiv:0909.3304•Physical Review Letters•Oct 4, 2010•1,097 citations

Estimation of low rank density matrices by Pauli measurements

Dong Xia

Density matrices are positively semi-definite Hermitian matrices with unit trace that describe the states of quantum systems. Many quantum systems of physical interest can be represented as high-dimensional low rank density matrices. A popular problem in quantum state tomography (QST) is to estimate the unknown low rank density matrix of a quantum system by conducting Pauli measurements. Our main contribution is twofold. First, we establish the minimax lower bounds in Schatten $p$-norms with $1\leq p\leq+\infty$ for low rank density matrices estimation by Pauli measurements. In our previous paper [14], these minimax lower bounds are proved under the trace regression model with Gaussian noise and the noise is assumed to have common variance. In this paper, we prove these bounds under the Binomial observation model which meets the actual model in QST. Second, we study the Dantzig estimator (DE) for estimating the unknown low rank density matrix under the Binomial observation model by using Pauli measurements. In our previous papers [14] and [25], we studied the least squares estimator and the projection estimator, where we proved the optimal convergence rates for the least squares estimator in Schatten $p$-norms with $1\leq p\leq2$ and, under a stronger condition, the optimal convergence rates for the projection estimator in Schatten $p$-norms with $1\leq p\leq+\infty$. In this paper, we show that the results of these two distinct estimators can be simultaneously obtained by the Dantzig estimator. Moreover, better convergence rates in Schatten norm distances can be proved for Dantzig estimator under conditions weaker than those needed in [14] and [25]. When the objective function of DE is replaced by the negative von Neumann entropy, we obtain sharp convergence rate in Kullback-Leibler divergence.

arXiv:1610.04811•Electronic Journal of Statistics•Jan 1, 2017•4 citations

Compressed-sensing tomography for qudits in Hilbert spaces of non-power-of-two dimensions

Revanth Badveli, Vinayak Jagadish, R. Srikanth +1 more

The techniques of low-rank matrix recovery were adapted for quantum state tomography (QST) previously by Gross et al. [Phys. Rev. Lett. 105, 150401 (2010)] where they consider the tomography of $n$ spin-1/2 systems. For the density matrix of dimension $d={2}^{n}$ and rank $r$ with $r\ensuremath{\ll}{2}^{n}$, it was shown that randomly chosen Pauli measurements of the order $O[dr\phantom{\rule{0.16em}{0ex}}log{(d)}^{2}]$ are enough to fully reconstruct the density matrix by running a specific convex optimization algorithm. The result utilized the low operator norm of the Pauli operator basis, which makes it ``incoherent'' to low-rank matrices. For quantum systems of dimension $d$ not a power of two, Pauli measurements are not available, and one may consider using $\mathrm{SU}(d)$ measurements. Here, we point out that the $\mathrm{SU}(d)$ operators, owing to their high operator norm, do not provide a significant savings in the number of measurement settings required for successful recovery of all rank-$r$ states. We propose an alternative strategy in which the quantum information is swapped into the subspace of a power-two system using only $\mathrm{poly}[log{(d)}^{2}]$ gates at most with QST being implemented, subsequently, by performing $O[dr\phantom{\rule{0.16em}{0ex}}log{(d)}^{2}]$ Pauli measurements. We show that, despite the increased dimensionality, this method is more efficient than the one using $\mathrm{SU}(d)$ measurements.

arXiv:2006.01803•Physical Review A•Jun 22, 2020•3 citations

Optimal Low-Rank Quantum State Tomography with Bounded-Sample Joint Measurements

Ashwin Nayak, Xingyu Zhou

We determine the optimal sample complexity of low-rank quantum state tomography when each measurement may act jointly on at most $t$ samples. For sufficiently small $\varepsilon$, estimating an unknown state on $\mathbb{C}^d$ of rank at most $r$ to trace norm error $\varepsilon$ with constant success probability requires, and is achievable with, $$Θ\left(\frac{dr}{\varepsilon^2}\mathop{\mathrm{max}}\left\{1,\frac{r}{\sqrt{t}}\right\}\right)$$ samples. The lower bound allows the protocol to choose each joint measurement adaptively using all previous classical outcomes; the matching upper bound is nonadaptive. Thus joint measurements on at most $t$ samples improve the complexity of algorithms making single-sample measurements by at most a factor $\sqrt{t}$. Further, measuring order $r^2$ samples jointly is necessary and sufficient to attain the unrestricted collective rate. For the lower bound, we vary the support of a state with fixed uniform spectrum and bound the Fisher information trace of every joint measurement on $t$ samples. The adaptive Fisher chain rule and the van Trees inequality then give the trace norm lower bound. For the upper bound, we construct and analyze a nonadaptive tomography protocol based on a Gaussian joint measurement. An explicit second moment identity and a conditional Gaussian law outside the state's support give a rank-dependent error analysis, yielding the matching rate.

quant-pharXiv:2609.10514•Sep 10, 2026•0 citations

Experimental quantum compressed sensing for a seven-qubit system

Carlos A. Riofrío, D. Gross, Steven T. Flammia +4 more

Well-controlled quantum devices with their increasing system size face a new roadblock hindering further development of quantum technologies. The effort of quantum tomography-the reconstruction of states and processes of a quantum device-scales unfavourably: state-of-the-art systems can no longer be characterized. Quantum compressed sensing mitigates this problem by reconstructing states from incomplete data. Here we present an experimental implementation of compressed tomography of a seven-qubit system-a topological colour code prepared in a trapped ion architecture. We are in the highly incomplete-127 Pauli basis measurement settings-and highly noisy-100 repetitions each-regime. Originally, compressed sensing was advocated for states with few non-zero eigenvalues. We argue that low-rank estimates are appropriate in general since statistical noise enables reliable reconstruction of only the leading eigenvectors. The remaining eigenvectors behave consistently with a random-matrix model that carries no information about the true state.

arXiv:1608.02263•Nature Communications•May 17, 2017•151 citations
Intuitively, if a density operator has small rank, then it should be easier to estimate from experimental data, since in this case only a few eigenvectors need to be learned. We prove two complementary results that confirm this intuition. Firstly, we show that a low-rank density matrix can be estimated using fewer copies of the state, i.e. the sample complexity of tomography decreases with the rank. Secondly, we show that unknown low-rank states can be reconstructed from an incomplete set of measurements, using techniques from compressed sensing and matrix completion. These techniques use simple Pauli measurements, and their output can be certified without making any assumptions about the unknown state. In this paper, we present a new theoretical analysis of compressed tomography, based on the restricted isometry property for low-rank matrices. Using these tools, we obtain near-optimal error bounds for the realistic situation where the data contain noise due to finite statistics, and the density matrix is full-rank with decaying eigenvalues. We also obtain upper bounds on the sample complexity of compressed tomography, and almost-matching lower bounds on the sample complexity of any procedure using adaptive sequences of Pauli measurements. Using numerical simulations, we compare the performance of two compressed sensing estimators—the matrix Dantzig selector and the matrix Lasso—with standard maximum-likelihood estimation (MLE). We find that, given comparable experimental resources, the compressed sensing estimators consistently produce higher fidelity state reconstructions than MLE. In addition, the use of an incomplete set of measurements leads to faster classical processing with no loss of accuracy. Finally, we show how to certify the accuracy of a low-rank estimate using direct fidelity estimation, and describe a method for compressed quantum process tomography that works for processes with small Kraus rank and requires only Pauli eigenstate preparations and Pauli measurements.
This paper establishes the methods for quantum state tomography based on compressed sensing.These methods are specialized for quantum states that are fairly pure,and they offer a significant performance improvement on large quantum systems.In particular,they are able to reconstruct an unknown density matrix of dimension d and rank r using O(rdlog2d) measurement settings,compared to standard methods that require d2 settings.The methods have several features that make them amenable to experimental implementation: they require only simple Pauli measurements,use fast convex optimization,are stable against noise,and can be applied to states that are only approximately low rank.The acquired data can be used to certify that the state is indeed close to pure,so no a priori assumptions are needed.
We establish methods for quantum state tomography based on compressed sensing. These methods are specialized for quantum states that are fairly pure, and they offer a significant performance improvement on large quantum systems. In particular, they are able to reconstruct an unknown density matrix of dimension d and rank r using O(rdlog²d) measurement settings, compared to standard methods that require d² settings. Our methods have several features that make them amenable to experimental implementation: they require only simple Pauli measurements, use fast convex optimization, are stable against noise, and can be applied to states that are only approximately low rank. The acquired data can be used to certify that the state is indeed close to pure, so no a priori assumptions are needed.
Density matrices are positively semi-definite Hermitian matrices with unit trace that describe the states of quantum systems. Many quantum systems of physical interest can be represented as high-dimensional low rank density matrices. A popular problem in quantum state tomography (QST) is to estimate the unknown low rank density matrix of a quantum system by conducting Pauli measurements. Our main contribution is twofold. First, we establish the minimax lower bounds in Schatten ppp-norms with 1≤p≤+∞1\leq p\leq+\infty1≤p≤+∞ for low rank density matrices estimation by Pauli measurements. In our previous paper [14], these minimax lower bounds are proved under the trace regression model with Gaussian noise and the noise is assumed to have common variance. In this paper, we prove these bounds under the Binomial observation model which meets the actual model in QST. Second, we study the Dantzig estimator (DE) for estimating the unknown low rank density matrix under the Binomial observation model by using Pauli measurements. In our previous papers [14] and [25], we studied the least squares estimator and the projection estimator, where we proved the optimal convergence rates for the least squares estimator in Schatten ppp-norms with 1≤p≤21\leq p\leq21≤p≤2 and, under a stronger condition, the optimal convergence rates for the projection estimator in Schatten ppp-norms with 1≤p≤+∞1\leq p\leq+\infty1≤p≤+∞. In this paper, we show that the results of these two distinct estimators can be simultaneously obtained by the Dantzig estimator. Moreover, better convergence rates in Schatten norm distances can be proved for Dantzig estimator under conditions weaker than those needed in [14] and [25]. When the objective function of DE is replaced by the negative von Neumann entropy, we obtain sharp convergence rate in Kullback-Leibler divergence.
The techniques of low-rank matrix recovery were adapted for quantum state tomography (QST) previously by Gross et al. [Phys. Rev. Lett. 105, 150401 (2010)] where they consider the tomography of nnn spin-1/2 systems. For the density matrix of dimension d=2nd={2}^{n}d=2n and rank rrr with $r\ensuremath{\ll}{2}^{n}$, it was shown that randomly chosen Pauli measurements of the order O[drlog(d)2]O[dr\phantom{\rule{0.16em}{0ex}}log{(d)}^{2}]O[drlog(d)2] are enough to fully reconstruct the density matrix by running a specific convex optimization algorithm. The result utilized the low operator norm of the Pauli operator basis, which makes it ``incoherent'' to low-rank matrices. For quantum systems of dimension ddd not a power of two, Pauli measurements are not available, and one may consider using SU(d)\mathrm{SU}(d)SU(d) measurements. Here, we point out that the SU(d)\mathrm{SU}(d)SU(d) operators, owing to their high operator norm, do not provide a significant savings in the number of measurement settings required for successful recovery of all rank-rrr states. We propose an alternative strategy in which the quantum information is swapped into the subspace of a power-two system using only poly[log(d)2]\mathrm{poly}[log{(d)}^{2}]poly[log(d)2] gates at most with QST being implemented, subsequently, by performing O[drlog(d)2]O[dr\phantom{\rule{0.16em}{0ex}}log{(d)}^{2}]O[drlog(d)2] Pauli measurements. We show that, despite the increased dimensionality, this method is more efficient than the one using SU(d)\mathrm{SU}(d)SU(d) measurements.
We determine the optimal sample complexity of low-rank quantum state tomography when each measurement may act jointly on at most ttt samples. For sufficiently small ε\varepsilonε, estimating an unknown state on Cd\mathbb{C}^dCd of rank at most rrr to trace norm error ε\varepsilonε with constant success probability requires, and is achievable with, Θ(drε2max{1,rt})Θ\left(\frac{dr}{\varepsilon^2}\mathop{\mathrm{max}}\left\{1,\frac{r}{\sqrt{t}}\right\}\right)Θ(ε2dr​max{1,t​r​}) samples. The lower bound allows the protocol to choose each joint measurement adaptively using all previous classical outcomes; the matching upper bound is nonadaptive. Thus joint measurements on at most ttt samples improve the complexity of algorithms making single-sample measurements by at most a factor t\sqrt{t}t​. Further, measuring order r2r^2r2 samples jointly is necessary and sufficient to attain the unrestricted collective rate. For the lower bound, we vary the support of a state with fixed uniform spectrum and bound the Fisher information trace of every joint measurement on ttt samples. The adaptive Fisher chain rule and the van Trees inequality then give the trace norm lower bound. For the upper bound, we construct and analyze a nonadaptive tomography protocol based on a Gaussian joint measurement. An explicit second moment identity and a conditional Gaussian law outside the state's support give a rank-dependent error analysis, yielding the matching rate.
Well-controlled quantum devices with their increasing system size face a new roadblock hindering further development of quantum technologies. The effort of quantum tomography-the reconstruction of states and processes of a quantum device-scales unfavourably: state-of-the-art systems can no longer be characterized. Quantum compressed sensing mitigates this problem by reconstructing states from incomplete data. Here we present an experimental implementation of compressed tomography of a seven-qubit system-a topological colour code prepared in a trapped ion architecture. We are in the highly incomplete-127 Pauli basis measurement settings-and highly noisy-100 repetitions each-regime. Originally, compressed sensing was advocated for states with few non-zero eigenvalues. We argue that low-rank estimates are appropriate in general since statistical noise enables reliable reconstruction of only the leading eigenvectors. The remaining eigenvectors behave consistently with a random-matrix model that carries no information about the true state.