Surmounting intrinsic quantum-measurement uncertainties in Gaussian-state tomography with quadrature squeezing

We reveal that quadrature squeezing can result in significantly better quantum-estimation performance with quantum heterodyne detection (of H. P. Yuen and J. H. Shapiro) as compared to quantum homodyne detection for Gaussian states, which touches an important aspect in the foundational understanding of these two schemes. Taking single-mode Gaussian states as examples, we show analytically that the competition between the errors incurred during tomogram processing in homodyne detection and the Arthurs-Kelly uncertainties arising from simultaneous incompatible quadrature measurements in heterodyne detection can often lead to the latter giving more accurate estimates. This observation is also partly a manifestation of a fundamental relationship between the respective data uncertainties for the two schemes. In this sense, quadrature squeezing can be used to overcome intrinsic quantum-measurement uncertainties in heterodyne detection.

Continuous-variable (CV) quantum tomography in phase space, or the diagnostics and reconstruction of observables of infinite Hilbert-space dimension in the continuous phase-space representation, is an indispensable technique for characterizing quantum light sources. These sources play a crucial role in practical optical quantum cryptography protocols 1-5 , many of which involve security analysis on Gaussian quantum states 6-10 -namely the coherent and squeezed states. Not only do Gaussian sources have properties that are intimately related to optomechanical phenomena 11,12 and extensively exploited in quantum metrology 13,14 , but they are also considered in the study of protocols such as state teleportation 15,16 , dense coding 17,18 and cloning 19,20 . The list of references is certainly non-exhaustive, and the reader is referred to two review articles on the subject of CV quantum optical communications 21,22 .
It therefore goes without saying that quantum tomography techniques for quantum states of light are of major interest in recent years. One of the most popular techniques, quantum homodyne detection [23][24][25] , is used in quantum optics to measure intensities of light signals from the outputs of a beam splitter that coherently merges the light mode from the source and that from a local oscillator, or reference coherent state. The result is an approximate measurement of the eigenstates of a rotated photonic quadrature whose phase angle depends on the phase of the local oscillator. The homodyne data obtained per binned angle constitute a distribution of points along a particular phase-space cut defined by this angle. This distribution is precisely the marginal distribution of the Wigner function of the quantum state of the source over the complementary quadrature. There is another well-known optical technique, heterodyne detection [26][27][28][29][30][31][32][33] , which has been extensively used to simultaneously probe a pair of optical beams of different frequencies in order to measure their relative phase. Compared to quantum homodyne detection, there are apparently relatively fewer published works on using its quantum variant to perform quantum tomography. This involves simultaneously measuring signal intensities of beams that are split from a single source signal by a beam splitter, thus realizing the approximate measurement of two rotated quadratures (position and momentum say) that are complementary to each other (see Fig. 1 for schematics). Because of the nature of such a quantum measurement, we shall understand the heterodyne scheme discussed here as a quantum double-homodyne tomography scheme, and henceforth, with common understanding, drop the adjective "quantum" when referring to these schemes. The heterodyne data obtained constitute a distribution of points according to the Husimi Q function of the state.
Manipulation of quantum tomograms have been carried out to address certain scalar quantities of interest (purity 34 for instance). Throughout the discussion, we shall focus on the reconstruction of covariance matrices that fully characterize the Wigner function of single-mode Gaussian states using the two CV schemes.
There exist elements that intrinsically affect the tomographic accuracy of quantum-state (or observable) reconstruction with typical CV tomography data. Two kinds of data disturbances that are ubiquitous in every experiment are statistical errors owing to the reconstruction with only a finite number of data points collected, as well as instrumental errors arising from errors in detector efficiency calibrations, measurement settings, etc. There are other tomographic elements that are particular to the two schemes of interest. In homodyne tomography, the phase space is sampled in terms of cuts resulted from marginalizations of the Wigner function over the complementary quadratures relative to the measured ones. The unknown quantum state, or any other field observable for that matter, is reconstructed by post-processing the homodyne data -also known as tomograms for historical reasons 35 -to reverse the marginalization in order to recover the full state. The heterodyne data, on the other hand, already consists of sample points distributed from the Q function in phase space. Nevertheless, it is well-known Figure 1. Schematics for homodyne and heterodyne detections. In our context, the heterodyne scheme shall be understood as a double-homodyne scheme (shown in (b)) that measures two complementary sets of quadrature eigenstates (those of position X and momentum P) simultaneously. that heterodyne tomography gives quantum measurement errors corresponding to the Arthurs-Kelly uncertainty relation, which originates from the simultaneous measurement of complementary quadrature observables 36,37 .
Two main goals are achieved in this discussion. Firstly, we reveal that for data collected along any pre-chosen phase-space direction, the data uncertainties for heterodyne detection typically become smaller relative to those for homodyne detection for highly-squeezed quantum states. This, as it turns out, is a result of an intimate relationship between the conditional variance of heterodyne data and the marginal variance of homodyne data, which defines the underlying statistical behavior for such data. In the context of covariance estimation, where covariances are the statistical quantities that completely characterize all Gaussian states, the difference in tomographic performance for covariance-matrix reconstruction between the two detection schemes is a sophisticated consequence of this fundamental relationship. This brings us to the second goal -the investigation of the optimal tomographic accuracy of unbiased covariance-matrix reconstruction for these two CV schemes. To appreciate the underlying physical framework, it will be shown that if the Arthurs-Kelly-type errors were absent, heterodyne detection on single-mode Gaussian states indeed always outperforms homodyne detection as far as the optimal limit of unbiased covariance-matrix reconstruction is concerned. In real scenarios where these errors are always present, this article shall demonstrate that heterodyne detection can still beat homodyne detection in optimal tomographic accuracy in many situations of practical relevance.

Results
Notations. The covariance matrix that characterizes the Wigner function of a single-mode Gaussian quantum state can be represented by the positive 2 × 2 real matrix involving expectation values of functions of the respective position-and momentum-quadrature standard deviations Δ X and Δ P. The matrix elements of G w specify the size, shape, and orientation of the ellipsoidal level contours of the Wigner function.
Data uncertainties along a phase-space direction. In this section, we shall perform a generic analysis on uncertainty regions for the two schemes. To begin, we recall two different notions of probability distribution functions from statistics that can be derived from a joint probability distribution function of a set of parameters taking random values. The first kind of distribution functions, the conditional distribution function, is a probability distribution function, which values are defined by a subset of these parameters after fixing the values for the rest of the parameters. The second kind, the marginal distribution function, is a probability distribution function of a subset of these parameters by taking an average over the rest of the parameters. Data sampling according to these two kinds of distributions respectively produce data uncertainties that are quite different. The data uncertainty that goes with the marginal distribution are never less than that which goes with the conditional distribution. Figure 2 illustrates this point for a Gaussian joint distribution function and the proof is straightforward. Let us consider a Gaussian conditional probability distribution of two random variables characterized by a two-dimensional covariance matrix G, and denote the marginal variance by 2 σ θ , and the conditional variance by 2 Σ θ . We shall fix the reference direction, defined by the angle θ and shown graphically as dashed straight lines in Fig. 2 The simple inequality implies that the uncertainty for data acquired from the marginal distribution is always greater than those acquired from the full conditional distribution. The equality holds only in the following two cases: (i) The integration, or shadow projection, is done along principal axes of the ellipse; (ii) The ellipse is actually a circle (zero eccentricity) so that all directions are principal. The preceding background introduction sets the stage for clarifying an important operational difference between homodyne and heterodyne detection. In a homodyne measurement, the data obtained for a fixed angle θ is sampled according to the marginal distribution of the Wigner function. On the other hand, data gathered from a heterodyne measurement for the same angle are sampled directly from the Q function, which is the conditional distribution for the standard pair of complementary quadrature variables. Owing to the additional vacuum contribution present in the heterodyne measurement, the conditional covariance of the Q function is larger than that of the Wigner function. As a consequence, the heterodyne-data uncertainty region is not always smaller than the homodyne-data uncertainty region. Nevertheless, for Gaussian states with highly elongated Wigner functions, the uncertainty region of the heterodyne data is indeed smaller than that of the homodyne data. In other words, if one picks any phase-space direction to collect both homodyne and heterodyne data originating from a source described by a highly elongated Wigner function, then on average, the heterodyne data have a smaller spread than the homodyne data. Details of the uncertainty region analysis are given in the Methods section under Uncertainty regions.
Hence apparently, heterodyne detection is the better scheme for Gaussian states with highly elongated Wigner functions, and that profile asymmetry seems to be the crucial factor for data collected along a particular phase-space direction. However, as we shall soon witness, the two different statistical nature of the CV measurement schemes, as discussed previously, entangle with other data-processing factors in the schemes in a much more complicated way in quantum-state estimation.
To quantify the optimal tomographic performance between the two CV schemes, one needs to directly investigate the uncertainties of relevant tomographic quantities of interest, and measures of optimal tomographic accuracies generally depend on aspects of these uncertainty regions in a highly convoluted manner. Regardless, behind this underlying complexity lies a similarly elegant conclusion: Quadrature squeezing improves the tomographic performance of heterodyne data over homodyne data, thereby surmounting the intrinsic Arthurs-Kelly measurement uncertainties. That is, for highly-squeezed single-mode Gaussian states, heterodyne detection almost always gives a better tomographic performance than homodyne detection.

Covariance estimation -main result.
To formally compare the optimal tomographic accuracy of any unbiased estimation of the covariance matrix, we adopt the well-known tomographic measure  = − F Sp{ } 1 -the matrix trace of the inverse of the Fisher information matrix F that is scaled with the number of copies N -that sets the Cramér-Rao bound for all unbiased matrix estimators 38,39 , the best possible tomographic performance any unbiased estimator can achieve. By denoting the bounds for For a joint distribution function over a two-variable (two-dimensional) space, the marginal variance for data sampled along a particular reference direction (indicated by a double-headed arrow) can be understood as a shadow cast from the joint uncertainty region in the orthogonal direction (indicated by the corresponding perpendicular pair of dashed lines). On the other hand, the conditional variance is directly obtained by slicing the joint uncertainty region about its center along the same reference direction. The marginal uncertainty region is the region defined by pairs of points bounding the marginal variances in all directions, which is similar to how the conditional uncertainty region is related to the conditional variances.
Scientific RepoRts | 5:12289 | DOi: 10.1038/srep12289 homodyne and heterodyne detection respectively by HOM  and HET  , and the detector efficiency that is common to both schemes by η ≤ 1, the results are summarized by the following closed-form expressions: Since the measure  is invariant under a rotation of G w , as it should, we may parametrize this matrix using only two relevant parameters: μ ≥ 1 that controls the size (related to the temperature) and λ ≥ 1 that controls the shape or the ratio of the major to the minor axis (related to the extent of quadrature squeezing). After establishing Eq. (4), all subsequent comparisons shall revolve around these equations.
Since we are concerned with only profile estimation, the center of the Gaussian Wigner function, that is its mean value, is neglected.
Covariance estimation -revealing the physical consequences. To gain insights in the fundamental difference between homodyne and heterodyne detection, let us first consider the hypothetical situation where there are no quantum-mechanical consequences in simultaneously measuring two incompatible observables -the absence of all Arthurs-Kelly-type uncertainties. This entails the equalities 0 and G HOM = G HET for the detector-efficiency terms and covariance matrices. The ratio is then a monotonically increasing function of λ, μ and η. It turns out that this function has a maximum value of one, which is attained in the limit λ, μ → ∞. The ratio γ is smallest when λ = μ = 1, taking the minimal value of 3/10 for all minimum-uncertainty states (μ = 1) with circular Wigner-function profiles (λ = 1), i.e. the coherent states. Extreme elongation of the profiles as a result of huge photonic quadrature squeezing (λ ≫ 1) renders both CV schemes equivalent in tomographic performance since the significant regions of sampling approach phase-space lines of infinite length, and details of the two schemes in this hypothetical setting are irrelevant within such infinitesimally thin regions. Figure 3 illustrates all the observations made.
The physical implication of these findings can be succinctly written as the following message: For all single-mode Gaussian states, regardless of the detector efficiency, a direct sampling of the phase-space Q function would in principle give more accurate covariance-matrix estimators as compared to their counterparts obtained through tomogram processing of data that follow the marginalized distributions of the Wigner function. For this class of quantum states, this is consistent with an intuition that a direct inversion of a set of statistically consistent data points would yield more accurate Gaussian Wigner-function estimation than an indirect two-step post-processing procedure on marginally consistent data points that incorporates tomogram combination and data inversion. Covariance estimation -analysis for realistic scenarios with perfect detections. In the physical reality, heterodyne detection is always accompanied by the Arthurs-Kelly measurement uncertainty that comes from simultaneously measuring two complementary quadrature observables. Physically, the action of the beam splitter (apart from the other two used to carry out the joint homodyne measurements) that splits the incoming source signal into two in a heterodyne measurement set-up introduces vacuum fluctuation in the other input photon mode. The result is an additional vacuum noise that contributes to the overall Arthurs-Kelly-type measurement error induced by such a measurement, whose lower limit is known to be larger than the usual Heisenberg's uncertainty lower limit 36,37 . The formalism of quasi-probability distributions automatically accounts for the vacuum fluctuation by noting that the covariance matrix G w for the Gaussian Wigner function of a particular Gaussian state (associated with homodyne detection) is always less than that for the Q function, G Q = G w + 1/(2η), by a multiple of the identity that takes the usual one-half value when η = 1 for perfect detections. This additional beam splitter thus plays the fundamental role in introducing the Arthurs-Kelly uncertainty upon a subsequent measurement. It is the combination of this additive vacuum term and the details in handling different types of data that dictates the executive difference between homodyne and heterodyne detection in state estimation through the respective covariance matrices G HOM and G HET in Eq. (4).
To analyze the consequence of this quantum-mechanical uncertainty, we use the correct expressions for G HOM and G HET and first consider the ideal situation where the detections are perfect (η = 1). Figure 4 shows the surface plot generated for Eq. (4). From the plot, we note the maximal influence on γ as a manifestion of the Arthurs-Kelly uncertainty for minimum-uncertainty states (μ = 1), where 1 The tomographic accuracy associated with heterodyne detection takes the worst-case magnitude 6 5 HET H OM HOM    = / > for coherent states and remains greater than unity for all squeezed states (λ > 1). For Gaussian states of higher temperatures (μ > 1) that are sufficiently squeezed, γ would eventually be smaller than unity, since in the range λ ≫ 1, it can be shown that the gradient of HOM  in λ is always steeper than that of HET  . The ratio γ approaches unity as λ goes to infinity for all μ, in agreement with the previous discussion above.
We conclude from the short excursion above that in this perfect-detection scenario, as long as there exist slight perturbation on a Gaussian source, just as in any realistic setting, such that the resulting quantum state can no longer be of minimum uncertainty, one can always benefit from heterodyne detection with sufficiently large squeezing.
Covariance estimation -analysis for realistic scenarios with imperfect detections. In practice, detections are never perfect due to losses, which implies that the detector efficiency η is always less than unity. The surface plots for γ, as shown in Fig. 5, correspondingly possess rather different shapes and curvatures for different values of η. The response to η for homodyne and heterodyne tomography schemes turn out to be quite different. For example, in the limit of very small efficiency, η ≪ 1, the additional η-dependent factor HET δ η ( ) for the heterodyne scheme is almost twice as big as the factor HOM δ η ( ) for the homodyne scheme. One should expect homodyne tomography to perform better in this limit. Notice that for realistic detection efficiencies in the range 0.2 < η < 0.8, heterodyne tomography always outperforms homodyne tomography except for a small range of μ and λ parameters. Figure 6 illustrates the performance of the two schemes in terms of uncertainty-ellipse reconstructions of a Gaussian source. Details of its generation are given in Methods.
Numerical calculations show that in this η-range, the critical values of μ and λ beyond which γ < 1 are respectively ≈ 1.736 and ≈ 3.771. These values correspond to a quadrature squeezing of about There exists a small λ-μ region within which optimal unbiased covariance-matrix estimators obtained with homodyne detection are more accurate than those obtained with heterodyne. In typical experimental conditions where Gaussian states prepared can neither give rise to minimum uncertainties nor be truely coherent states, there exist a plethora of settings for which HET H OM   < .
− 3.369 dB and an anti-quadrature squeezing of about 8.1601 dB from the shot-noise (vacuum) level, which is well below the state-of-the-art squeezing/anti-squeezing levels, which are respectively ≈ − 12.7 dB and ≈ 19.9 dB 13 . Thus, physicists designing narrowband squeezed Gaussian states with error-suppression of such ratios can enjoy the benefit of tomographically accurate covariance-matrix estimators with heterodyne detection. However, decreasing η further to very small values would push the surface up to a level of 6/5 (the ultimate limit of γ for single-mode Gaussian states) eventually, since in the small-η limit, the ratio of the asymptotic value of 6 HET 2  η ≈ / to that of 5 HOM 2  η ≈ / becomes a constant.

Discussion
Quantum homodyne and heterodyne (double-homodyne) techniques are well-known quantum diagnostic schemes that are accessible in quantum-optics laboratories. The respective tomographic capabilities of these two continuous-variable schemes involve an intricate combination of both the statistical and the data-processing characteristics in each of the schemes. We revealed these inherent tomographic capabilities of heterodyne detection by comparing it with the more popular homodyne detection and analyzing the optimal tomographic accuracy of unbiased covariance-matrix estimation of single-mode Gaussian states with the help of the scaled Fisher information matrix.  The uncertainty ellipses plotted are reconstructions from particular sets of data that are respectively from measuring N = 50 (red), 100 (green) and 150 (blue) copies or sampling events in phase space, which are compared to the true uncertainty ellipse (black, dashed). The true Gaussian squeezed state considered here is characterized by the parameters μ = 2 and λ = 10, with the detection efficiency η = 0.5. The thickness of each closed curve is proportional to the distance between the reconstructed ellipse and the true ellipse. These plots are representatives of the average performance of the two schemes, a signature of the existence of a wide range of settings for which heterodyne detection is tomographically better than homodyne detection even for moderate values of N.
Despite the existence of Arthurs-Kelly uncertainties as a consequence of simultaneous measurements of incompatible observables, heterodyne detection exhibits significantly better tomographic performance when probing sufficiently squeezed Gaussian states that are not of minimum-uncertainty, which is the case in all experiments, especially when the detection efficiency is within the practical range in quantum-optics experiments. The amount of quadrature squeezing considered here can be observed with present-day technology.
The seemingly counterintuitive fact that quadrature squeezing can improve heterodyne tomography over homodyne tomography by effectively mask the detrimental effects of intrinsic quantum measurement uncertainties is an important physical consequence that can be easily overlooked if one bases his or her understanding of optimal tomographic performance solely on the manifestation of these uncertainties. The intricate amalgam of statistical and data-processing characteristics in a highly sophisticated fashion is what determines the ultimate tomographic limit for these schemes. This interesting and important revelation justifies the relevance of future experimental work with quantum heterodyne detections to reveal possible novel tomographic enhancements. While we acknowledge the current technical complexities in carrying out heterodyne detection, on which we are in no position to comment, we hope that the conclusions drawn in this article would kick-start some interesting prototypical experimental work.

Methods
Uncertainty regions. To analyze the uncertainty regions more thoroughly, let us focus on the uncertainty regions of the two tomography schemes. These regions are centered at the origin of the phase space, and in terms of the angle θ, their respective boundaries are described by the functions σ θ and Σ θ . These are precisely the standard deviations for the two schemes of a fixed angle θ. For simplicity we assume ideal detector efficiency η = 1. Then, with Gauss-Weierstrass transform that turns Wigner functions into Q functions, and vice versa, the respective covariance matrices G w and G Q are related by the expression G Q = G w + 1/2.
For minimum uncertainty states, where The uncertainty region for homodyne tomography is more complicated, but for states of highly elongated regions, that is λ ≫ 1 or λ ≪ 1, this region looks like two circles in contact. Plots of these uncertainty regions are given below. Figure 7 shows that homodyne detection seems to give smaller errors for nearly circular states, while for states of highly elongated Wigner functions, uncertainty suppression is observed with heterodyne detection for almost all angles. Naturally, if we integrate/cut approximately along the direction of a principal axis, homodyne data is less noisy. However, the interval of angles for which σ θ < Σ θ shrinks as the Wigner function becomes more elongated. Figure 8 shows the uncertainty regions for the two schemes. The area of the heterodyne uncertainty ellipse is given by whereas that of the quartic homodyne uncertainty region is The latter is / λ 2 /( ) λ [1 2 ] times more than the former in the limit of λ ≫ 1 [λ ≪ 1]. The areas become equal when One solution is the reciprocal of the other as it should be, for the values λ = λ crit and λ = 1/λ crit essentially correspond to the same Gaussian state apart from a rotation by a π/2 angle. Generalization of the above to inefficient detections -η < 1 -is possible and closed form expressions for the uncertainty areas are available. Figure 9 shows the plot of the critical value λ crit for which the homodyne and heterodyne areas are equal as a function of the detector efficiency η.
Smaller detection efficiencies mean larger differences between the corresponding Wigner and Q functions and more state asymmetry (larger λ crit in other words) is needed to see the effect. A realistic choice of η = 0.8 gives λ = 0.149 which requires the smallest standard deviation to be squeezed by a factor of 2.6 below the shot noise, well within the experimentally feasible range for optical quadrature squeezing.
Needless to say, the minimum-uncertainty states discussed above are the worst possible states for this kind of comparison since the extra noise from the Q function stands out. With a lot of excess noise (mixed states), as we always have in actual experiments, where the anti-squeezed quadrature is never controlled as well as the squeezed one, this extra vacuum term matters less and the uncertainty suppression in heterodyne tomography over homodyne tomography will be easy to observe.
Covariance estimation -background. We are interested in estimating the true covariance matrix of the Wigner function that is parametrized as can be written as a linear combination of trace-orthonormal Hermitian basis matrices Γ k -Sp{Γ k Γ l } = δ kl .
Let G ŵ be an estimator of G w . We quantify the performance of this estimator using the scaled mean squared-error (MSE) with N that is defined to be the mean Hilbert-Schmidt distance between the true covariance matrix and the estimator, where F is the scaled Fisher information matrix that is represented by a 3 × 3 positive real matrix. The right-hand side of (14) is known as the Cramér-Rao bound. By denoting the row of parameters as    is the scaled Fisher matrix of a single quadrature, which is really a special case of Eq. (15), and the normalization by d is there to account for splitting the total ensemble among d quadratures. Finally, for a fair comparison in phase space, we take the limit of infinitely many quadrature angle settings and replace the sum by the integral. This also closely follows the experimental practice, where rather than measuring a fixed number of quadratures, the phase of the local oscillator is continuously being changed during the experiment, so that effectively a very large number of quadratures (up to 10 6 ) is realized. Thereafter, We calculate the scaled Fisher matrix and scaled MSE for both homodyne and heterodyne schemes. The parametrization of Eq. (12) is used. The easiest way to do this is to set g 3 = 0, and after which, employ the invariance properties of the scaled Fisher matrix. Once the final scaled MSE formulas are cast in basis-independent form, that is, in terms of matrix invariants, they become general.
Indeed, in heterodyne detection, any 2 × 2 covariance matrix can be diagonalized by applying a rotation R, according to which a new trace-orthonormal operator basis is generated: Γ k → Γ′ k = RΓ k R T , RR T = R T R = 1 and Sp{Γ′ k Γ′ l } = δ kl . Expressing the old basis in terms of the new basis elements, Γ k = ∑ l w kl Γ′ l , the scaled Fisher matrix transforms into F → W T FW, where orthonormality of the new basis implies that WW T = W T W = 1. Therefore, the scaled MSE, which is given by the matrix trace of F −1 , does not change. Alternatively, one can rotate the covariance matrix, rather than the coordinate system, and see that the right-hand side of Eq. (15) does not change upon the mapping G → R T GR.
In the case of homodyne detection, the invariance with respect to a change in the orientation of the uncertainty ellipse follows from integrating the scaled Fisher information over all quadratures (angles).
For homodyne measurements over the entire phase space, there is, unfortunately, no closed-form expression for the LIN/ML estimator, even for large N. However, one can still compute the scaled MSE in the large-N limit from the scaled Fisher matrix in Eq. (20). This requires calculating integrals over all angles between zero and π, which can be carried out elegantly using contour integration techniques.
Upon recognizing the form of C(θ) in Eq. (16)   For this type of tomography, the optimal LIN estimator is well-known. It is essentially given by the total sample covariance matrix of all the collected data two-tuples {(x j , p j )}, up to an additive η-dependent multiple of the identity. The scaled MSE of this LIN estimator (also of the ML estimator for large N) can easily be computed either by calculating the scaled Fisher matrix (g 3 = 0) using Eq. and then take its inverse, or by directly performing the average over all data of the Hilbert-Schmidt distance in Eq. (13) under the consideration of 2-D Gaussian statistics. The two approaches are equivalent since direct-sampling of the Q function results in 2-D Gaussian statistics. Either way, the closed-form expression, written in the basis-independent form, is the one given as the second equation of (4).
Covariance estimation -Comparisons of reconstructed uncertainty ellipses. Figure 6 compares reconstructed uncertainty ellipses for a given true squeezed Gaussian state using the two CV schemes. The parameters for this state are chosen so that γ < 1 in the large N limit. Here, the potential of heterodyne detection can be witnessed even for moderate values of N.
To generate the curves, we investigate unbiased maximum-likelihood (ML) covariance-matrix estimators G ML that are in principle asymptotically optimal with respect to the Cramér-Rao bound for the which we take as the measure for tomographic accuracy.
The ML estimator for the homodyne data is obtained by simply taking all data obtained from various angles and performing a maximization of the likelihood, which is a Gaussian distribution function characterized by elements of the covariance matrix, to obtain the covariance matrix that maximizes this likelihood. To obtain the ML estimator for the heterodyne data, the sampled phase-space points are gathered and the corresponding sample covariance matrix is calculated from these points. This sample covariance matrix is, by definition, the maximum-likelihood estimator for Gaussian states. From these covariance matrices, uncertainty ellipses can be directly obtained by computing the eigenvalues and eigenvectors of these matrices, where the two orthonormal eigenvectors of each matrix represent the basis vectors for the two principle axes, and the eigenvalues represent the lengths of these axes.