The tomographic reconstruction of the state of a quantum-mechanical system is an essential component in the development of quantum technologies. We present an overview of different tomographic methods for determining the quantum-mechanical density matrix of a single qubit: (scaled) direct inversion, maximum likelihood estimation (MLE), minimum Fisher information distance and Bayesian mean estimation (BME). We discuss the different prior densities in the space of density matrices, on which both MLE and BME depend, as well as ways of including experimental errors and of estimating tomography errors. As a measure of the accuracy of these methods, we average the trace distance between a given density matrix and the tomographic density matrices it can give rise to through experimental measurements. We find that the BME provides the most accurate estimate of the density matrix, and suggest using either the pure-state prior, if the system is known to be in a rather pure state, or the Bures prior if any state is possible. The MLE is found to be slightly less accurate. We comment on the extrapolation of these results to larger systems.
The tomographic reconstruction of the state of a quantum-mechanical system is an essential component in the development of quantum technologies. We present an overview of different tomographic methods for determining the quantum-mechanical density matrix of a single qubit: (scaled) direct inversion, maximum likelihood estimation (MLE), minimum Fisher information distance and Bayesian mean estimation (BME). We discuss the different prior densities in the space of density matrices, on which both MLE and BME depend, as well as ways of including experimental errors and of estimating tomography errors. As a measure of the accuracy of these methods, we average the trace distance between a given density matrix and the tomographic density matrices it can give rise to through experimental measurements. We find that the BME provides the most accurate estimate of the density matrix, and suggest using either the pure-state prior, if the system is known to be in a rather pure state, or the Bures prior if any state is possible. The MLE is found to be slightly less accurate. We comment on the extrapolation of these results to larger systems.
Summary form only given. Experimental techniques for measurement of the quantum state of light have been the subject of intensive investigation for some time. Tomographic techniques have been applied to experiments such as the homodyne measurement of the Wigner function of a single mode of light and of the density matrix of the polarization degrees of freedom of a pair of entangled photons. In this technique, the density matrix (or Wigner function) which characterizes the quantum state of the system being measured is found from a linear transformation of experimental data. There are a number of drawbacks to the method, principally in that the recovered state might not, because of experimental noise, correspond to a physical state. For example, density matrices for any quantum state must be Hermitian, positive semi-definite matrices with unit trace. The tomographically measured matrices often fail to be positive semi-definite. To avoid this problem, the maximum likelihood approach to the estimation of quantum states has been developed. We evaluate two different approaches to the maximum likelihood technique, namely defining the likelihood function f as a weighted sum of squared variances from measured data, and defining it in terms of information content. These are compared with results of standard tomographical schemes.
In standard optical tomographic methods, the off-diagonal elements of a density matrix ρ are measured indirectly. Thus, the reconstruction of ρ, even if it is based on linear inversion, typically magnifies small errors in the experimental data. Recently, an optimal tomography solution measuring all the elements of ρ one-by-one without error magnification has been theoretically proposed. We implemented this method for two-qubit polarization states. For comparison, we also experimentally implemented other well-known tomographic protocols, either based solely on local measurements (of, e.g., the Pauli operators and James-Kwiat-Munro-White projectors) or with mutually unbiased bases requiring both local and global measurements. We reconstructed seventeen separable, partially and maximally entangled two-qubit polarization states. Our experiments show that our method has the highest stability against errors in comparison to other quantum tomographies. In particular, we demonstrate that each optimally-reconstructed state is embedded in an uncertainty circle of the smallest radius, both in terms of trace distance and disturbance. We explain how to experimentally estimate uncertainty radii for all the implemented tomographies and show that, for each reconstructed state, the relevant uncertainty circles intersect indicating the approximate location of the corresponding physical density matrix.
In this paper, we derive analytic expressions for the starting (initial) values of the parameters of the T-matrix that is frequently employed in the construction of a theoretical density matrix in a Maximum Likelihood Estimate (MLE) procedure; this optimization procedure is used to fit the experimental data pertaining to single qubit state tomography. Appropriate starting values are critical in achieving a global minimum in the fitting process. We also indicate an analytic way of making the experimentally determined density matrix physical without resorting to the MLE process, if optimality in quantum state tomography can be disregarded, since the ultimate goal is optimality in quantum process tomography.
We present a Bayesian approach to the problem of estimating density matrices in quantum state tomography. A general framework is presented based on a suitable mathematical formulation, where a study of the convergence of the Monte Carlo Markov Chain algorithm is given, including a comparison with other estimation methods, such as maximum likelihood estimation and linear inversion. This analysis indicates that our approach not only recovers the underlying parameters quite properly, but also produces physically acceptable punctual and interval estimates. A prior sensitive study was conducted indicating that when useful prior information is available and incorporated, more accurate results are obtained. This general framework, which is based on a reparameterization of the model, allows an easier choice of the prior and proposal distributions for the Metropolis–Hastings algorithm.
In quantum information transformation and quantum computation, the most critical issues are security and accuracy. These features, therefore, stimulate research on quantum state characterization. A characterization tool, Quantum state tomography, reconstructs the density matrix of an unknown quantum state. Theoretically, reconstructing an unknown state using this method can be arbitrarily accurate. However, this is less practical owing to the huge burden of measurements and data processing for large numbers of qubits. Even comprising an efficient estimator and a precise algorithm, an optimal tomographic framework can also be overburdened owing to the exponential growth of the measurements. Moreover, the consequential postprocessing of huge amounts of data challenges the capacity of computers. Thus, it is crucial to build an efficient framework that requires fewer measurements but yields an expected accuracy. To this end, we built a tomography schema by which only a few Pauli measurements enable an accurate tomographic reconstruction. Subsequently, this schema was verified as efficient and accurate through numerical simulations on the tomography of multi-qubit quantum states. Furthermore, this schema was proven to be robust through numerical simulations on a noisy superconducting qubit system. Therefore, the tomography schema paves an alternatively effective way to reconstruct the density matrix of a quantum state owing to its efficiency and accuracy, which are essential for quantum state tomography.
Recent experimental progress in the preparation and control of quantum systems has brought to light the importance of Quantum State Tomography (QST) in validating the results. In this thesis we investigate several aspects of QST, whose central problem is to devise estimation schemes for the recovery of an unknown state, given an ensemble of n independent identically prepared systems. The key issues in tackling QST for large dimensional systems is the construction of physically relevant low dimensional state models, and the design of appropriate measurements. Inspired by compressed sensing tomography, in chapters 4, 5 we consider the statistical problem of estimating low rank states (r ≪ d) in the set-up of Multiple Ions Tomography (MIT), where r and d are the rank and the dimension of the state respectively. We investigate how the estimation error behaves with a reduction in the number of measurement settings, compared to ‘full’ QST in two setups - Pauli and random bases measurement designs. We study the estimation errors in this ‘incomplete’ measurement setup in terms of a concentration of the Fisher information matrix. For the random bases design we demonstrate that O(r logd) settings suffice for the mean square error w.r.t the Frobenius norm to achieve the optimal O(1/n) rate of estimation. When the error functions are locally quadratic, like the Frobenius norm, then the expected error (or risk) of standard procedures achieves this optimal rate. However, for fidelity based errors such as the Bures distance we show that no ‘compressive’ recovery exists for states close to the boundary, and it is known that even with conventional ‘full’ tomography schemes the risk scales as O(1/√n) for such states and error functions. For qubit states this boundary is the surface of the Bloch sphere. Several estimators have been proposed to improve this scaling with ‘adaptive’ tomography. In chapter 6 we analyse this problem from the perspective of the maximum Bures risk over all qubit states. We propose two adaptive estimation strategies, one based on local measurements and another based on collective measurements utilising the results of quantum local asymptotic normality. We demonstrate a scaling of O(1/n) for the maximum Bures risk with both estimation strategies, and also discuss the construction of a minimax optimal estimator. In chapter 7 we return to the MIT setup and systematically compare several tomographic estimators in an extensive simulation study. We present and analyse results from this study, investigating the performance of the estimators across various states, measurement designs and error functions. Along with commonly used estimators like maximum likelihood, we propose and evaluate a few new ones. We finally introduce two web-based applications designed as tools for performing QST simulations online.