Nearly Optimal Measurement Scheme in a Noisy Mach-Zehnder Interferometer with Coherent and Squeezed Vacuum

The use of an interferometer to perform an ultra-precise parameter estimation under noisy conditions is a challenging task. Here we discuss nearly optimal measurement schemes for a well known,sensitive input state, squeezed vacuum and coherent light. We find that a single mode intensity measurement, while the simplest and able to beat the shot-noise limit, is outperformed by other measurement schemes in the low-power regime. However, at high powers, intensity measurement is only outperformed by a small factor. Specifically, we confirm, that an optimal measurement choice under lossless conditions is the parity measurement. In addition, we also discuss the performance of several other common measurement schemes when considering photon loss, detector efficiency, phase drift, and thermal photon noise. We conclude that, with noise considerations, homodyne remains near optimal in both the low and high power regimes. Surprisingly, some of the remaining investigated measurement schemes, including the previous optimal parity measurement, do not remain even near optimal when noise is introduced.


INTRODUCTION
Typical parameter estimation with the use of interferometric schemes aims to estimate some unknown parameter which is probed with the input quantum states of light. In principle, the sensitivity of these measurements depends on the chosen input states of light, the interferometric scheme, the noise encountered and the detection scheme performed at the output. For a real-world example, perhaps the most sensitive of these types of interferometers are the large scale interferometers used as gravitational wave sensors [1][2][3][4][5][6][7]. In general, if classical states of light are used, then the most sensitive measurement is limited to a classical bound, the shot-noise limit (SNL) [8][9][10]. Despite the remarkable precision possible with classical states, improvements are still possible. Here we discuss nearly optimal measurements achievable when one considers input states of coherent and squeezed vacuum [11,12], under many common noisy conditions and in realistic power regimes which are applicable to general interferometry.
It is of practical interest to consider the difficulty with implementing any particular measurement scheme as every additional optical element introduces further loss into the interferometer. It has been previously shown that the parity measurement is one example of an optimal measurement for coherent and squeezed vacuum input states under lossless conditions [13]. It was also previously shown that a more involved detection scheme is optimal under photon loss [14]. However, here we will discuss various common detection schemes, which are easier to implement in practice and perform nearly optimal. Discussion of a lossy MZI for Fock state inputs is also discussed in previous works [15,16].
While there are many technical challenges in using squeezed states of light, we show here that some of the measurement techniques commonly used in a classical setup are no longer near optimal. In addition, some measurements exhibit problems with effects such as phase drift and thermal photon noise. With the goal of choosing a simple, yet well-performing measurement, we investigate homodyne [17], parity measurements [18][19][20][21][22] and compare them to a standard intensity measurements. These measurements form a set that are either simple to implement, or are known to be optimal in the lossless case. Specifically, we confirm that, under lossless conditions, the parity measurement achieves the smallest phase variance. However, under noisy conditions, surprisingly the parity measurement suffers greatly, while the homodyne measurement continues to give a nearly optimal phase measurement. The parity measurement under losses was briefly discussed in the context of entangled coherent states by Joo et al. [23]. For the lossless case, we divide our results into two regimes, the low power regime (|α| 2 < 500, e.g. small scale sensors), in which different detection schemes can lead to significantly different phase variances, and the high power regime (|α| 2 > 10 5 , e.g. large scale, devoted interferometry), where all detection schemes are nearly optimal. While our scheme may hint at applications for setups like LIGO, a much more focused analysis, outside the scope of our investigation, would be required before drawing conclusions about LIGO's performance.

METHOD
The interferometer considered here is a Mach-Zehnder interferometer (MZI) [24] as shown in Fig. 1 and is mathematically equivalent to a Michelson interferometer. Here, an input of a coherent state (|α ) and squeezed vacuum (|χ ) is used. With this input state, it is known that the phase sensitivity can be below the SNL, typically defined as ∆ 2 φ SNL = 1/N , where N is the mean number of photons entering the MZI [11].
FIG. 1: A general Mach-Zehnder interferometer with coherent |α and squeezed vacuum |χ states as input. Beam splitters (BS) mix the two spatial modes, while mirrors (M) impart a phase shift, which can be safely neglected since it it common to both modes. A phase shift φ represents the phase difference between the two arms of the MZI, which can be due to a path length difference. Our goal is to estimate the unknown parameter φ, which corresponds to the interaction of the quantum state with some process of interest.
For its close connection to the parity measurement, we shall describe our states in terms of Wigner functions. One can construct any Gaussian states Wigner function directly from the first and second moments by way of, where the covariance matrix, σ = σ ij = X i X j +X j X i − 2 X i X j , mean vector, d j = X j and X i , X j are orthogonal phase space variables.
We use this general form for our chosen input states of a coherent state (|α ) and squeezed vacuum (|χ = re iδ ) to define our states by, Here α, θ are the coherent amplitude and phase, respectively while r, δ denote the squeezing parameter and phase. As the input state we consider is a product state, it can be written in terms of the product [25], For simplicity, both states have equal initial phases, as this gives rise to the optimal phase sensitivity (discussed later) and are taken to be θ = δ = 0. This simply defines the coherent state to be displaced in the x 1 direction and the squeezed state to be squeezed along x 2 [26]. The average photon number in the coherent state is N coh = |α| 2 and in the squeezed vacuum state N sqz = sinh 2 r, which sets the SNL to be ∆ 2 φ SNL = 1/N tot = 1/(|α| 2 + sinh 2 r).
The propagation of this Wigner function is accomplished by the transformation of the phase space variables through the MZI, dictated by its optical elements. These transformations are described by where both beam splitters are fixed to be 50-50 and φ represents the unknown phase difference between the two arms of our MZI. We have chosen to use a symmetric phase model in order to simplify calculations as well as agree with previous results [27,28]. Our goal then will be minimizing our uncertainty in the estimation of the unknown parameter φ. Using these transforms, the total transform for the phase space variables is given by, where {x jf , p jf } represent the phase space variables, for each mode, after propagation through the MZI. We can also consider photon loss in the model by way of two mechanisms, photon loss to the environment inside the interferometer and photon loss at the detectors, due to inefficient detectors. Both of these can be modeled by placing a fictitious beam splitter in the interferometer with vacuum and a interferometer arm as input and tracing over one of the output modes, to mimic loss of photons to the environment [29]. This linear photon loss mechanism can be modeled with the use of a relatively simple transform, since these states are all of Gaussian form. Specifically this amounts to a transform of the covariance matrix according to σ L = (1 − L)I · σ + LI, 0 ≤ L ≤ 1 is the combined photon loss and I is the 4×4 identity matrix. Similarly the mean vector is transformed [14,30].

Quantum Cramér-Rao Bound
We consider an optimal measurement scheme with the meaning of saturating the quantum Cramér-Rao bound (QCRB) [31,32], which gives the best phase sensitivity possible for a chosen interferometer setup and input states. This optimality is independent of measurement scheme and it remains a separate task to show which measurement scheme achieves this optimal bound [12]. In what we call the classical version of this setup, a coherent state and vacuum state are used as input. With these two input states, the best sensitivity one can achieve is bounded by the SNL, which is achievable with many different detection schemes. Many interferometer models mainly focus on analytical analysis of Fisher information [33,34] when there is loss and phase drift. While this analysis is useful in that it demonstrates a 'best case scenario', it is unknown whether the optimal detection scheme is hard to realize in an actual experimental setup. Thus, in our analysis, we are more focused on Fisher information and how it compares with specific detection schemes, under noisy conditions. The benefit of using squeezed vacuum in place of vacuum is then that the phase measurement can now reach below the SNL. In order to compare various choices of measurement schemes, we not only need to calculate the various measurement outcomes, but also need to show the best sensitivity attainable with these input states. The best phase measurement one can do is given by the quantum Cramér-Rao bound [31] and is related to the quantum Fisher information (QFI, F ) [32], simply by ∆ 2 φ QCRB = F −1 . For the input states of a coherent and squeezed vacuum, one can use the Schwinger representation to calculate the QFI, since these are pure states [27,35]. Another option, and the method we use here, instead utilizes the Gaussian form of the states and can be calculated directly in terms of covariance and mean [36][37][38]. This method applies to pure and mixed states, as long as it maintains Gaussian form. Using this formalism, the QCRB for a coherent state and squeezed vacuum into an MZI can be found to be [12,13], While this gives us a bound on the best sensitivity obtainable with these given input states, it does not directly consider loss or even tell us which detection scheme attains this bound.

Specific Measurements under lossless conditions
Now that we have a bound on the best possible sensitivity, we now seek to show how various choices of measurement compare to this bound. We consider some standard measurement choices including single-mode intensity, intensity difference, homodyne, and parity. While each of these measurements would require a significant reconfiguration of any interferometer, it is worthwhile to show how each choice impacts the resulting phase sensitivity measurement. We utilize the bosonic creation and annihilation operators (â † ,â), which obey the commutation relation, [â,â † ] = 1. We also utilize the quadrature operators (x,p) which are related to the creation and annihilation operators by the transformâ j = 1 √ 2 (x j + ip j ). These quadrature operators obey a similar commutator, [x,p] = i.
In terms of our output Wigner function, where "sym" indicates that this integral calculates the symmetric ordered expectation value of the operatorÔ. Each measurement operator, Ô , gives rise to a phase uncertainty by way of Starting with the simplest measurement, an intensity measurement is given by,Ô = â †â = x 2 +p 2 /2, which is implemented by simply collecting the outgoing light, directly onto a detector. For homodyne detection,Ô =x (we find the optimal homodyne measurement is taken along the x quadrature). For a balanced homodyne detection scheme, one would impinge one of the outgoing light outputs onto a 50-50 beam splitter, along with a coherent state of the same frequency as the input coherent state (usually this is derived from the same source) and perform intensity difference between the two outputs of this beam splitter. While there exist other implementations of homodyne than we describe here, we choose a standard balanced homodyne scheme, for simplicity. A standard intensity difference is defined asÔ =â †â −b †b . This particular measurement choice is also explored in Ref. [39]. Parity detection is defined to beÔ = (−1) â †â = πW (0, 0) ≡ Π . Parity detection has been implemented experimentally, though focusing on its ability for super-resolution [40]. While all chosen measurements can surpass the SNL, in the lossless case, to various degrees, in order of improving phase sensitivity, single-mode intensity performs the worst, followed by intensity difference, homodyne, and finally parity. The analytical forms of each detection scheme, at their respective minima, are listed below and we confirm that, under lossless conditions, the parity measurement matches the QCRB [13], homodyne attains, and intensity difference attains, while a single mode intensity measurement attains a minimum of, We can notice that for high coherent state powers (|α| 2 1), each detection scheme's leading term in its respective phase variance is given by ∆ 2 φ all ≈ (|α| 2 e 2r ) −1 , which is nearly optimal since the sinh 2 (r) term in the QCRB is negligible compared to large α. From these forms then, we can say that in the low-photon-number regime (|α| 2 < 500), the difference in these detection schemes can be significant, but in the high photon number regime (|α| 2 > 10 5 ), there is little difference between the various detection schemes.

Lossy Inteferometer
We now consider the effects of loss and calculate the lossy QCRB. This is done following the same loss procedure described previously. The lossy QCRB of this mixed state becomes [14] ∆ 2 φ Loss QCRB = L(e 2r − 1) + 1 (1 − L){|α| 2 e 2r + sinh 2 (r)[L(e 2r − 1) Note that this QCRB with loss only considers linear photon loss caused by photon loss inside the interferometer and photon loss due to inefficient detectors. In reality, there may be more specific sources of noise one needs to consider, but our method's purpose is to show a preliminary case when simple loss models are considered. We note that a measurement scheme proposed by Ono and Hofmann is exactly optimal (thus it is able to achieve the bound given by Eq. 11) under loss [14], but we wish to explore how simpler measurement schemes perform when compared with this bound.
In the case of losses, the forms of each phase variance necessarily becomes much less appealing. For this reason, we only list the analytical form of the homodyne measurement under loss, as it is our prime candidate for a nearly optimal measurement. The phase variance of homodyne in a lossy interferometer is given by, where we have fixed the optimal phase to φ = π to obtain the phase variance minimum. We can note several interesting comparisons from this form, including the obvious ∆ 2 φx(L) ≥ ∆ 2 φx and ∆ 2 φx(0) = ∆ 2 φx. However, if we investigate Eq. 11 for high powers (large |α|), we find, This expansion illustrates the fact that homodyne is nearly optimal and approaches the QCRB in the large power limit. One can see in Fig. 2, which shows when loss is considered, parity detection suffers greatly, while other detection schemes are still able to achieve sub-SNL phase variances. In all but the intensity and parity measurement schemes, the optimal phase (the point at which each curve achieves its minimum) has a constant value and therefore should not prove overly difficult to stabilize. In the case of intensity and parity measurement however, this optimal phase depends on both the squeezing strength r and the amplitude of the coherent state |α|. Therefore, fluctuations in the source will actually affect the optimal phase setting and in general degrade the phase measurement in this measurement scheme. Note that, in practice, typical experiments use an offset to remain near these optimum values, but purposely remain slightly away from the minimum, due to noise considerations. At this point we can note, that current technological limits enforce |α| 2 sinh 2 (r), as generally it is relatively easier to increase laser power, than to increase squeezing power. Just as it was in the lossless case, under lossy conditions then, Fig. 2 shows that homodyne remains nearly optimal in the low power regime. In contrast, the previously optimal measurement, parity, is now not able to even reach sub-SNL.
We can also plot the phase variance as a function of average photon number, shown in Fig. 3 for large powers, which can be related to the light's optical frequency Note that a homodyne measurement nearly, but not exactly, reaches the QCRB as shown by the inset. and power by |α| 2 = P/( ω 0 ) [41]. The phase variances shown in Fig. 3 are at their respective minima in terms of optimal phase. In this form, it's clear that a parity measurement suffers greatly, under lossy conditions and at high powers. Parity may also be difficult to implement in certain inteferometric setups as it either involves number counting (which is not feasible at very large powers) or several homodyne measurements [42]. Alternatively, a single homodyne measurement is nearly optimal in this lossy case, still only requires measurement on a single mode, is simpler to implement than parity, and is not nearly as sensitive to loss. We note that while homodyne appears to meet the QCRB in Fig. 2, it actually doesn't exactly reach the QCRB (as indicated in the inset of Fig 2). While intensity difference is also close in phase variance to a homodyne measurement (when |α| 2 > 100) it requires utilization of both output modes for phase measurement, which may not be feasible in some setups. We note that while we have chosen typical parameters for |α| and r, the trend of homodyne achieving near optimal measurement generalizes to other parameter choices as well.

Phase Drift
Returning to Fig. 2, it is clear at which value of phase the various measurements attain their lowest value. It is this value of phase that one attempts to always take measurements at with the use of a control phase inside the interferometer. The width of each of curve then can be interpreted as the chosen measurement scheme's resistance to phase drift. The mechanism of phase drift comes about due to the limited ability to set control phases in the interferometer with infinite precision. In general, the control phase value will vary around the optimal phase We have assumed one can set the control phase to its optimal value, to obtain the best phase variance in each measurement choice. Squeezing strength in the squeezed state is set to r = 1. Note that Parity is now not able to achieve even sub-SNL, due to loss, while homodyne and intensity difference quickly approach the QCRB (appear on top of one another). SNL and QCRB are also plotted with the same lossy parameters.
setting. For this reason we aim to show this phase drift in a more rigorous way. We therefore will use the analytical forms of the various measurement phase variances, as a function of unknown phase φ. We simulate phase drift by computing a running average of the phase variance, with a pseudo-randomly chosen phase, near the optimal phase, for each measurement. This pseudo-random choice is made from a Gaussian distribution, whose mean is fixed at each measurements respective optimal phase choice and has a chosen variance of σ = 0.15. As predicted in the previous discussion, this gives a clearer picture of each measurement's behavior under phase drift. For simplicity, we focus on the lossless case for this treatment of phase drift. Shown in Fig. 4, we see the phase variance ratio to the QCRB for each measurement scheme, as a function of the number of measurements. As the number of measurements is increased, the phase variance asymptotes to the ideal measurement case, given by the phase variance at the optimal phase. This is an illustration of the law of large numbers, as the number of measurement increases, each scheme approaches its true average. However, as is clear from Fig. 4, each scheme approaches its average at significantly different rates. We can see that parity performs fairly poorly, as compared to the other measurement schemes. In the case of intensity, homodyne, and intensity difference measurements, it's clear that homodyne and intensity difference attain a small phase variance, while also being more tolerant of phase drift. This confirms a special case of Genoni et al. [43], who showed that homodyne measurement is resistant to phase diffusion in pure Gaussian states. In principle, all of the different measurement schemes will each attain their respective phase variance minimum, as the number of measurements increases to infinity, but it is instructive to see how quickly a finite number of measurements approaches the ideal phase variance minimum.

Thermal Photon Noise
In addition to photon loss, detector efficiency, and phase drift, we also model the inevitable interaction with thermal photon noise from the environment. This is accomplished much in the same way as a photon loss model, but here we consider a thermal photon state incident on a fictitious beam splitter, on both arms of the interferometer, and trace out one of its output modes. This allows a tunable amount of thermal photon noise (by changing the average photon number in the thermal state), into the interferometer. The effects of this unwanted thermal photon noise, to the various measurement's phase variance, is shown in Fig. 5. From this, we can see that even in the regime of a relatively low photon number of thermal photon noise, such noise significantly degrades the phase variance of each scheme, but drastically affects the par-ity scheme, making it significantly above the SNL. The SNL and QCRB under this noise model do not directly incorporate the additional thermal photons. Also in this regime, a standard single-mode intensity measurement now does not acheive sub-SNL phase variance, but homodyne and intensity difference continue to reach sub-SNL. We also note that the advantage of homodyne over intensity difference measurement is significantly decreased in the presence of thermal photon noise, but homodyne still maintains its superiority. The introduction of larger average thermal photon number continues to degrade all measurements so that they no longer beat the SNL, but this example showcases their behavior under this noise model. It should be noted that in the optical regime, the occupation of a thermal state, at room temperature is approximately n th ≈ 10 −20 and therefore, some interferometric schemes do not deal with significant contribution from this model of thermal photon noise, but experiments in the microwave frequencies can have n th ≈ 1, where this model is more applicable. We recommend that a homodyne measurement is the simplest, nearly optimal measurement choice for a setup as discussed here. Homodyne is a typical measurement choice in interferometer experiments, as well as being a single mode measurement, likely resistant to photon loss, detector efficiency, and phase drift. It shows its main benefits in the low power regime, but performs nearly optimal in both the low and high power regimes.

CONCLUSION
In this paper, we have seen the performance of many common interferometric measurement schemes. While all are able to achieve a sub-SNL phase variance measurement in the lossless case, for the choice of a coherent and squeezed vacuum input state, all are outperformed by a homodyne measurement when loss is in-troduced. While these measurements each come with their own challenges in implementing, we have shown that each measurement's performance can vary significantly under different noise models. We have also shown that in the high-photon regime, with loss, most measurement schemes approach the QCRB except for parity which suffers significantly. Our results may imply that simpler measurement schemes are overall appealing when using large powers. The behavior of each measurement scheme under phase drift and thermal photon noise is also discussed, and we find that homodyne and intensity difference measurement behave best within these models. This should be expected as both homodyne and intensity difference measurements operate in a similar way, subtracting intensities between two modes, removing common noise sources. Therefore, when considering ease of implementation as well as near optimal detection, we conclude that homodyne is nearly optimal under loss, phase drift, and thermal photon noise, for the specific choice of input states of coherent and squeezed vacuum and in both power regimes.
B.T.G. would like to acknowledge support from the National Physical Science Consortium & National Institute of Standards and Technology graduate fellowship program as well as helpful discussions with Dr. Emanuel Knill at NIST-Boulder. C.Y. would like to acknowledge support from an Economic Development Assistantship from the National Science Foundation and the Louisiana State University System Board of Regents. D.K.M. would like to acknowledge support from University Grants Commission, New Delhi, India for Raman Fellowship. T.R.C. would like to acknowledge support from the National Science Foundation grants PHY-1150531. This document has been assigned the LIGO document number LIGO-P1600084. J.P.D would like to acknowledge support from the Air Force Office of Scientific Research, the Army Research Office, the Boeing Corporation, the National Science Foundation, and the Northrop Grumman Corporation. We would all also like to thank Dr. Haixing Miao and the MQM LIGO group for helpful discussions. All authors contributed equally to this work.

Declarations Appendix 1: Measurement
We have shown that different choices of measurements lead to varied ability to perform parameter estimation. Here we show the details of each detections analytical calculation for the ideal, lossless case. For each detection choice,Ô, we need to calculate ∆ 2 φ = ∆ 2Ô /|∂ Ô /∂φ| 2 = ( Ô 2 − Ô 2 )/|∂ Ô /∂φ| 2 (13) that is, we need the variance of the chosen measurement and the derivative of its first moment. Therefore, in general we need the first and second moments for each chosen detection scheme. We also need the Wigner function at the output, this is obtained by following the transformation of phase space variables described in the text and results in a final output Wigner function of, W (x 1 , p 1 , x 2 , p 2 ) = 1 π 2 Exp{− 1 2 e −2r [p 2 2 + x 2 1 + e 4r (p 2 1 + x 2 2 ) + e 2r (p 2 1 + p 2 2 + x 2 1 + x 2 2 + 4|α| 2 + 4 √ 2|α|(x 2 cos (φ/2) − p 1 sin (φ/2))) + 2e 2r sinh r((p 2 2 − x 2 1 + e 2r cos φ(p 1 − x 2 )(p 1 + x 2 )) + 2 sin φ(p 2 x 1 + e 2r p 1 x 2 ))]}, (14) where subscripts label spatial modes, α is the amplitude in the coherent state and r the squeezing strength in the squeezed state and we have chosen θ coh = δ sqz = 0 for simplicity. For an intensity measurement, in terms of phase space quadrature operatorsx,p, for the second moment, we have, where "sym" denotes the symmetric operator form which is calculated from (â †â ) 2 sym = ∞ −∞ 1 4 (x 2 + p 2 ) 2 × W (X)dX and W (X) is the output Wigner function, given by Eq. 14. For the variance calculation then, For homodyne detection, since the optimal homodyne measurement is alongx, we simply have, For intensity difference we have, and finally for parity measurement, utilizing Π 2 = 1 and W (0, 0) is the value of the Wigner function at the origin in phase space.