We study quantum process tomography given the prior information that the map is a unitary or close to a unitary process. We show that a unitary map on a $d$-level system is completely characterized by a minimal set of ${d}^{2}+d$ elements associated with a collection of POVMs, in contrast to the ${d}^{4}\ensuremath{-}{d}^{2}$ elements required for a general completely positive trace-preserving map. To achieve this lower bound, one must probe the map with particular sets of $d$ pure states. We further compare the performance of different estimators inspired by compressed sensing, to reconstruct a near-unitary process from such data. We find that when we have accurate prior information, an appropriate compressed sensing method reduces the required data needed for high-fidelity estimation, and different estimators applied to the same data are sensitive to different types of noise. Convex optimization inspired by the techniques of compressed sensing can therefore be used both as indicators of error models and to validate the use of the prior assumptions.
We study quantum process tomography given the prior information that the map is a unitary or close to a unitary process. We show that a unitary map on a d-level system is completely characterized by a minimal set of d2+d elements associated with a collection of POVMs, in contrast to the d4\ensuremath−d2 elements required for a general completely positive trace-preserving map. To achieve this lower bound, one must probe the map with particular sets of d pure states. We further compare the performance of different estimators inspired by compressed sensing, to reconstruct a near-unitary process from such data. We find that when we have accurate prior information, an appropriate compressed sensing method reduces the required data needed for high-fidelity estimation, and different estimators applied to the same data are sensitive to different types of noise. Convex optimization inspired by the techniques of compressed sensing can therefore be used both as indicators of error models and to validate the use of the prior assumptions.
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.
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.
The development of large-scale platforms for quantum information requires new methods for verification and validation of quantum behavior. Quantum tomography (QT) is the standard tool for diagnosing quantum states, process, and readout devices by providing complete information. However, QT is limited since it is expensive to not only implement experimentally, but also requires heavy classical post-processing of data. In this dissertation, we introduce new methods for QT that are more efficient to implement and robust to noise and errors, thereby making QT a more practical tool for current quantum information experiments. The crucial detail that makes these new, efficient, and robust methods possible is prior information about the quantum system. This prior information is prompted by the goals of most experiments in quantum information, which require pure states, unitary processes, and rank-1 POVM operators. Therefore, most experiments are designed to operate near this ideal regime. We show that when this is the case, QT can be accomplished with significantly fewer resources, and produce a robust estimate in the presence of noise and errors. Moreover, the estimate is also robust if the state is not exactly pure, the process is not exactly unitary, or the POVM is not exactly rank-1. Such compelling methods are only made possible by the positivity constraint on quantum states, processes, and POVMs. This requirement is an inherent feature of quantum mechanics, but has powerful consequences to QT. Since QT is necessarily an experimental tool, we discuss a test of these new methods in an experimental setting. The physical system is an ensemble of laser-cooled cesium atoms in the laboratory of Prof. Poul Jessen. Experiments were conducted by Hector Sosa-Martinez et al. to demonstrate different QT protocols. We compare the results, and conclude that the new methods are effective.
Quantum process tomography is the task of reconstructing unknown quantum channels from measured data. In this work, we introduce compressed sensing-based methods that facilitate the reconstruction of quantum channels of low Kraus rank. Our main contribution is the analysis of a natural measurement model for this task: We assume that data is obtained by sending pure states into the channel and measuring expectation values on the output. Neither ancillary systems nor coherent operations across multiple channel uses are required. Most previous results on compressed process reconstruction reduce the problem to quantum state tomography on the channel's Choi matrix. While this ansatz yields recovery guarantees from an essentially minimal number of measurements, physical implementations of such schemes would typically involve ancillary systems. A priori, it is unclear whether a measurement model tailored directly to quantum process tomography might require more measurements. We establish that this is not the case.Technically, we prove recovery guarantees for three different reconstruction algorithms. The reconstructions are based on a trace, diamond, and ℓ2 -norm minimization, respectively. Our recovery guarantees are uniform in the sense that with one random choice of measurement settings all quantum channels can be recovered equally well. Moreover, stability against arbitrary measurement noise and robustness against violations of the low-rank assumption is guaranteed. Numerical studies demonstrate the feasibility of the approach.
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.
Drawing inspiration from gradient-descent methods developed for data processing in quantum state tomography [\href{https://iopscience.iop.org/article/10.1088/2058-9565/ae0baa}{Quantum Sci.~Technol.~\textbf{10} 045055 (2025)}] and quantum process tomography [\href{https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.130.150402}{Phys.~Rev.~Lett.~\textbf{130}, 150402 (2023)}], we introduce stochastic gradient descent (SGD) algorithms for fast quantum measurement tomography (QMT), applicable to both discrete- and continuous-variable quantum systems -- thus completing the tomography trio. A measurement device or detector in a quantum experiment is characterized by a set of positive operator-valued measure (POVM) elements; the goal of QMT is to estimate these operators from experimental data. To ensure physically valid (positive and complete) POVM reconstructions, we propose two distinct parameterization schemes within the SGD framework: one leveraging optimization on a Stiefel manifold and one based on Hermitian operator normalization via eigenvalue scaling. Within the SGD-QMT framework, we further investigate two loss functions: mean squared error, equivalent to L2 or Euclidean norm, and average negative log-likelihood, inspired by maximum likelihood estimation. We benchmark performance against state-of-the-art constrained convex optimization methods. Numerical simulations demonstrate that, compared to standard methods, our SGD-QMT algorithms offer significantly lower computational cost, superior reconstruction fidelity, and enhanced robustness to noise. We make a Python implementation of the SGD-QMT algorithms publicly available at \href{https://github.com/agtomo/SGD-QMT}{github.com/agtomo/SGD-QMT}.