Results concerning the construction of quantum Bayesian error regions as a means to certify the quality of parameter point estimators have been reported in recent years. This task remains numerically formidable in practice for large dimensions and so far, no analytical expressions of the region size and credibility (probability of any given true parameter residing in the region) are known, which form the two principal region properties to be reported alongside a point estimator obtained from collected data. We first establish analytical formulas for the size and credibility that are valid for a uniform prior distribution over parameters, sufficiently large data samples and general constrained convex parameter estimation settings. These formulas provide a means to an efficient asymptotic error certification for parameters of arbitrary dimensions. Next, we demonstrate the accuracies of these analytical formulas as compared to numerically computed region quantities with simulated examples in qubit and qutrit quantum-state tomography where computations of the latter are feasible.
Results concerning the construction of quantum Bayesian error regions as a means to certify the quality of parameter point estimators have been reported in recent years. This task remains numerically formidable in practice for large dimensions and so far, no analytical expressions of the region size and credibility (probability of any given true parameter residing in the region) are known, which form the two principal region properties to be reported alongside a point estimator obtained from collected data. We first establish analytical formulas for the size and credibility that are valid for a uniform prior distribution over parameters, sufficiently large data samples and general constrained convex parameter estimation settings. These formulas provide a means to an efficient asymptotic error certification for parameters of arbitrary dimensions. Next, we demonstrate the accuracies of these analytical formulas as compared to numerically computed region quantities with simulated examples in qubit and qutrit quantum-state tomography where computations of the latter are feasible.
Standard computation of size and credibility of a Bayesian credible region for certifying any point estimator of an unknown parameter (such as a quantum state, channel, phase, etc.) requires selecting points that are in the region from a finite parameter-space sample, which is infeasible for a large dataset or dimension as the region would then be extremely small. We solve this problem by introducing the in-region sampling theory to compute both region qualities just by sampling appropriate functions over the region itself using any Monte Carlo sampling method. We take in-region sampling to the next level by understanding the credible-region capacity (an alternative description for the region content to size) as the average l_{p}-norm distance (p>0) between a random region point and the estimator, and present analytical formulas for p=2 to estimate both the capacity and credibility for any dimension and a sufficiently large dataset without Monte Carlo sampling, thereby providing a quick alternative to Bayesian certification. All results are discussed in the context of quantum-state tomography.
Bayesian error analysis paves the way to the construction of credible and plausible error regions for a point estimator obtained from a given dataset. We introduce the concept of region accuracy for error regions (a generalization of the point-estimator mean squared error) to quantify the average statistical accuracy of all region points with respect to the unknown true parameter. We show that the increase in region accuracy is closely related to the Bayesian region dual operations in Shang et al (2013 New J. Phys. 15 123026). Next with only the given dataset as viable evidence, we establish various adaptive methods to maximize the region accuracy relative to the true parameter subject to the type of reported Bayesian region for a given point estimator. We highlight the performance of these adaptive methods by comparing them with non-adaptive procedures in three quantum-parameter estimation examples. The results of and mechanisms behind the adaptive schemes can be understood as the region analog of adaptive approaches to achieving the quantum Cramér–Rao bound for point estimators.
Standard Bayesian credible-region theory for constructing an error region on the unique estimator of an unknown state in general quantum-state tomography to calculate its size and credibility relies on heavy Monte Carlo sampling of the state space followed by filtering to obtain the correct region sample. This conventional methodology typically gives negligible yield for very small error regions originating from large data sets. In this article, we discuss at length the in-region sampling theory for computing both size and credibility from region-average quantities that avoids this general problem altogether. Among the many possible numerical choices, we study the performance and properties of accelerated hit-and-run Monte Carlo algorithm for in-region sampling and provide its complexity estimates for quantum states. Finally with our in-region concept, by alternatively quantifying the region capacity with the region-average distance between two states in the region (measured for instance with either the Hilbert-Schmidt, trace-class, or Bures distance), we derive approximation formulas to analytically estimate both region capacity and credibility without Monte Carlo computation.
An estimator is a state that represents one's best guess of the actual state of the quantum system for the given data. Such estimators are points in the state space. To be statistically meaningful, they have to be endowed with error regions, the generalization of error bars beyond one dimension. As opposed to standard ad hoc constructions of error regions, we introduce the maximum-likelihood region—the region of largest likelihood among all regions of the same size—as the natural counterpart of the popular maximum-likelihood estimator. Here, the size of a region is its prior probability. A related concept is the smallest credible region—the smallest region with pre-chosen posterior probability. In both cases, the optimal error region has constant likelihood on its boundary. This surprisingly simple characterization permits concise reporting of the error regions, even in high-dimensional problems. For illustration, we identify optimal error regions for single-qubit and two-qubit states from computer-generated data that simulate incomplete tomography with few measured copies.
Precise measurement of physical quantities plays a crucial role in the development of science and technology.The main purpose of the dissertation is two-fold: to investigate the ultimate precision for estimation of physical quantities using Gaussian states and to propose an efficient method for certification of Bayesian error region in general quantum parameter estimation.In the first part, we begin with analyzing sensitivity for estimating a phase difference in an optical interferometer.Optical interferometry is widely used in science and industry for measuring small displacements.Recently, a large-scale optical interferometer socalled the Laser interferometer Gravitational-Wave Observatory (LIGO) has succeeded in detecting a gravitational wave, the signal of which is extremely small.On the other hand, it has been shown that a non-classical feature of quantum states can improve the sensitivity of estimation, such as in optical interferometer, including the LIGO.From a practical point of view, we inspect the practically achievable precision using non-classical Gaussian states in Mach-Zehnder interferometer with feasible measurements and realistic photon loss.We then investigate the precision of single-mode phase estimation using Gaussian measurement, which can be realized by using homodyne detection, and i show that non-Gaussian measurement is necessary to utilize the power of Gaussian input probes optimally.Finally, we find the optimal measurement for general Gaussian quantum metrology and identify three distinct optimal measurements corresponding to different circumstances.In the second part, we study the Bayesian error region, which is a crucial concept for a general estimation process.When estimating a physical quantity, one has to supply the error interval (single-parameter) or error region (multi-parameter) as well as the estimate.However, it has been shown that as the dimension of quantum systems of interest grows, it becomes intractable to calculate the size and credibility of Bayesian error regions.As an alternative, we derive an analytical expression for the properties, the size and credibility, of Bayesian error regions, in an asymptotic regime.We then propose an efficient numerical method to calculate them for high-dimensional quantum systems even in a non-asymptotic regime.
Quantum state tomography is the task of inferring the state of a quantum system by appropriate measurements. Since the frequency distributions of the outcomes of any finite number of measurements will generally deviate from their asymptotic limits, the estimates computed by standard methods do not in general coincide with the true state and, therefore, have no operational significance unless their accuracy is defined in terms of error bounds. Here we show that quantum state tomography, together with an appropriate data analysis procedure, yields reliable and tight error bounds, specified in terms of confidence regions-a concept originating from classical statistics. Confidence regions are subsets of the state space in which the true state lies with high probability, independently of any prior assumption on the distribution of the possible states. Our method for computing confidence regions can be applied to arbitrary measurements including fully coherent ones; it is practical and particularly well suited for tomography on systems consisting of a small number of qubits, which are currently in the focus of interest in experimental quantum information science.