arxiv: v1 [quant-ph] 29 Sep 2018

Transcription

1 Quantum interferometry via a coherent state mixed with a photon-added squeezed vacuum state Shuai Wang 1,, Xuexiang Xu 3, Yejun Xu 4, Lijian Zhang 1 School of Mathematics and Physics, Jiangsu University of Technology, Changzhou 13001, P.R. China Corresponding author: wshslxy@cczu.edu.cn College of Engineering and Applied Sciences, anjing University, anjing 10093, P.R. China 3 Department of Physics, Jiangxi ormal University, anchang 3300,P.R. China and 4 School of Mechanical and Electronic Engineering, Chizhou University, Chizhou 47000, P.R. China arxiv:181131v1 quant-ph] 9 Sep 018 We theoretically investigate the phase sensitivity with parity detection on a Mach-Zehnder interferometer with a coherent state combined with a photon-added squeezed vacuum state. When the phase shift approaches zero, the squeezed vacuum state is indeed the optimal state within a constraint on the average number of photons. However, when the phase shift to be estimated slightly deviates from zero, the optimal state is neither the squeezed vacuum state nor the photonsubtracted squeezed vacuum state, but the photon-added squeezed vacuum state when they carry many photons. Finally, we show that the quantum Cramér-Rao bound can be reached by parity detection. PACS number(s): 4.50.Dv,03.65.Ta I. ITRODUCTIO Phase estimation and optimal interferometry play a significant role for many precision measurement applications. For a Mach-Zehnder interferometer (MZI), when only coherent light is injected into one input port of the first beam splitter, the other input port is by default the vacuum of light, the sensitivity of the phase estimation is limited by the standard quantum noise limit (SL), i.e., φ 1/, where the average photon number in the input beam 1, ]. When using a nonclassical input state, for example, a coherent light and a squeezed vacuum light are injected into the two input ports of a MZI, Caves 3] find that the sensitivity of phase estimation below the SL. Since then, in order to go beyond the SL, many highly nonclassical states are employed to reduce the phase uncertainty 3 15], and to approach φ 1/, the so-called Heisenberg limit () 16, 17]. In general, the sensitivity of phase estimation within an interferometer crucially depends on input states as well as detection schemes. Coherent states and squeezed vacuum states as well as Fock states are useful for metrology under the present experimental technology. For example, via the analysis of the quantum Fisher information, for a MZI with a Fock state in one input port and an arbitrary state with the same total average photon number in the other input, same phase uncertainties will be achieved and can approach to the 18]. Following the theoretical work in Ref. 18], we analytically prove that the quantum Cramer-Rao bound can be reached via the parity detection in the limit ϕ 0 19]. By Bayesian analysis of the photon number statistics of the output state in a MZI, Pezzé and Smerzi show that the sensitivity of phase estimation can be achieved when the coherent light and squeezed vacuum light are mixed in roughly equal intensities 0]. As a simple alternative to the detection scheme of Ref.0], parity detection for a MZI with coherent and squeezed vacuum light has bee also investigated 1]. A parity measurement simply measures the even or odd number of photons in the output mode ]. Recently, Lang and Caves have considered the question: given that one input of an interferometer is entered by coherent light, what is the best state to inject the other input port for achieving high-sensitivity phase-shift measurements within a constraint on the average photon number that the state can carry? The answer, they find, is the squeezed vacuum state () 3]. As just pointed out in the conclusion in Ref. 3], it needs to further investigate whether the is the optimal state when the squeezing light carries many photons, even carries as many or more photons than the coherent input. On the other hand, when the phase shift to be estimated slightly deviates from zero, whether the with many photons is also the optimal state? For fixed initial squeezing parameter, Birrittella and Gerry claim that the corresponding sensitivity of the phase estimation can be increased via parity detection when the mixing of coherent states and photon-subtracted (PS) as input states of the MZI 4]. However, within a constraint on the total average photon number, whether the PS can indeed improve the phase sensitivity? In this work, we will answer the above questions by investigating the interferometry performed by mixing a coherent state with non-gaussian squeezed states, such as a photon-added (PA) and the PS. At present, the best experimentally realized non- Gaussian squeezed state in quantum optics is the photonsubtracted squeezed states 5, 6], and photon subtraction can be implemented by a beam splitter with high transmissivity 5]. For the photon addition operation,

2 Agarwal and Tara 6] first theoretically studied the nonclassical properties of a photon-added coherent state. In 004, the photon addition operation was successfully demonstrated experimentally via a non-degenerate parametric amplifier with small coupling strength 8]. We also compared the phase sensitivity with another quantum limit, the quantum Cramér-Rao bound 9], which sets the ultimate limit for a set of probabilities that originated from measurements on a quantum system. The organization of this paper is as follows. In Sec. II, we make a brief review about both the PA and the PS. And then we describes the propagation of a twomode light, initially in the product state of PA and coherent light, through the MZI. Section III focuses on the parity detection scheme and provides the phase sensitivity. We will show that for giving a constraint on the average photon numbers of PA/PS/, almost the same phase sensitivity will be achieved by parity detection with such kinds of input states. Especially, when the phase shift slightly deviates from zero, the optimal state is neither the nor the PS, but the PA when they carry many photons. In Sec. IV, we prove that the quantum Cramer-Rao bound can be reached via the parity detection in the limit ϕ 0. II. RESULTED OUTPUT STATE OF THE IPUT FIELDS THROUGH THE ITERFEROMETER Previously, Birrittella and Gerry 4] studied the prospect of parity-based interferometry with mixing a PS and a coherent state. In contrast to that scheme, we mainly investigate the scheme when a PA and a coherent state ( z with the amplitude parameter z = z e θ ) are considered as the input state of a MZI. For our purposes, we first provide a brief review of both the PA and the PS which is a kind of non-gaussian squeezed vacuum states. And then, we derive the resulted output state when a PA and a coherent state are injected into a balanced MZI A. Photon-added and photon-subtracted squeezed vacuum states The is a Gaussian state, which is defined as 30] ( r = S (r) 0 = sech 1/ r exp 1 ) b tanh r 0, (1) where S (r) = exp r ( b b ) / ] with the squeezing parameter r. By repeatedly operating the photon addition operator b on a, one obtain the normalized PA which is defined as r, k = 1 k b k S (r) 0, () where k is the corresponding normalization factor 31] k = k t k τ k exp sinh r ( t + τ ) ] + tτ cosh r 4 t,τ=0 = k! cosh k rp k (cosh r). (3) The average photon number of the PA is given by n a = b b PA = k+1 k 1. (4) Different from that in Ref.4], we adopt the expression of the PS as following r, l = 1 Cl b l S (r) 0, (5) where b is the photon subtraction operator, and the normalization factor is 3] C l = l t l τ l exp sinh r ( t + τ ) ] + tτ sinh r t,τ=0 4 = l! ( i sinh r) l P l (i sinh r). (6) The average photon number of the PS is given by n s = b b PS = C k+1 C k. (7) Particularly, when k = l, the single-photon PA and the single-photon PS are the same non-gaussian squeezed state 33]. For the quantum metrology, the average photon number of input states is an important factor. Actually, both photon addition and subtraction can increase the average photon number of the state as shown in Fig. 1. Clearly, it is always better to perform addition rather than subtraction in order to increase the average photon number for given the initial squeezing. And hence, we expect that the PA can offer an improved phase resolution over both the PS-coherent input state and -coherent input state for a given initial squeezing parameter of the. The fact that the photon subtraction can increase the average photon number of the states with super- Poissonian statistics was explained in detail very recently by Stephen et al 34] and references involved in that. Indeed, the distribution for the single-mode, as is well known, is super-poissonian 35]. B. Resulted output state of the PA through the interferometer The balanced MZI considered here is mainly composed of two 50:50 beam splitters and two phase shifters. Generaly, the first beam splitter BS1 is described by the transformation U BS1 = exp iπ ( a b + ab ) /4 ].And the operator representation of the second beam splitter BS is taken as U BS = exp iπ ( a b + ab ) /4 ]. The operator U (ϕ) = exp iϕ ( a a b b ) / ] represents the two phase

3 For the balanced MZI, applying the following transformation relations e iϕj a e iϕj = a cos ϕ + b sin ϕ, 3 e iϕj b e iϕj = b cos ϕ a sin ϕ (11) and the relation e iϕj 0 a 0 b = 0 a 0 b, in principle, one can obtain the explicit form of the output state the MZI. For the convenience of the later calculation, we rewrite the PA ψ r,k b in the basis of the coherent state as follows The average photon number PS l=1 PA l= l=3 k=1 l=6 k= k=3 k=6 k =l= FIG. 1. Plots of the average photon number of both the PA and the PS, as well as the, respectively. (a). The average photon number for a given (fixed squeezing parameter r) as a function of the number k or l; (b) The average photon number as a function of the squeezing parameter r for some values of k and l. shifters, the angle ϕ being the phase shift between the two arms to be estimated. The unitary transformation associated with such balanced MZI can be written as 14] r U (ϕ) = e i π J1 e iϕj3 e i π J1 = e iϕj, (8) where these operators consisted of two sets of Bose operator J 1 = 1 ( a b + ab ), J = 1 ( a b ab ), i J 3 = 1 ( a a b b ), (9) ψ r,k b = sech1/ r k d k d α α dg k π e +gα tanh r α α g=0, (1) ] where α = exp α / + αb 0 is a coherent state. When the product state ψ in = z a ψ r,k b is injected into the MZI, the resulted output state can be written as ψ out = sech1/ r d k d α z k dg k π e α +gα tanh r α e (α sin ϕ +z cos ϕ )a +(α cos ϕ z sin ϕ )b 0, 0 g=0 (13), which is the state of light at the output of the MZI. In the following section, we shall mainly present parity measurement scheme with calculations of the expected value of the parity operator and the corresponding phase sensitivity. III. PHASE ESTIMATIO WITH PARITY DETECTIO There are several detection methods for extracting phase information from the output states of the MZI, e.g., intensity detection, homodyne detection, and parity detection ]. Parity detection simply measures the even or odd number of photons in the output mode. In addition, as shown in Ref. 1], the parity detection saturates the quantum Cramér-Rao bound and in turn provides the phase sensitivity when the -coherent state are mixed in equal proportions. In experiments, the parity detection using a photon-number resolving detector with coherent states 37] has also been demonstrated. Here, we use the parity detection too. are the angular momentum operators in the well-known Schwinger representation 36]. They satisfy the commutation relation J i, J j ] = iɛ ijk J k (i, j, k = 1, (, 3), and commute with the Casimir operator J 0 = 1 a a + b b ), i.e., J 0, J i ] = 0. Propagation of the input fields (pure states) through these elements, the resulted output state can be written as out MZI = e iϕj ψ in. (10) A. The parity detection with the PA-coherent state For the detailed discussion of the parity detection in quantum optimal metrology, one can review that in Ref.]. Actually, the parity detection is to obtain the expectation value of the parity operator in the output state of the MZI. In order to calculate conveniently, using the operator identity

4 4 exp ( λa a ) =: exp ( e λ 1 ) a a ] :, we rewrite the parity operator as follows Π b = ( 1) b b d = e iπb b =: e b b γ : = γ γ, π (14) where γ is a coherent state, and : : denotes the normally ordered form of Bose operators. Here, we consider performing parity detection on just one of the output modes, for instance, the b mode. The parity operator on an output mode b is described by Π B = ( 1) b b, then the expectation value of the parity operator is d γ Π b (ϕ) = out ψ π γ γ ψ out. (15) ow, we consider the corresponding the expectation value of the parity operator in the output state when the PA-coherent state is injected into the MZI. Then substituting Eqs. (13) and (14) into Eq. (15), and applying the integral formula d z π eζ z +ξz+ηz +fz +gz = 1 ζ 4fg e ζξη+ξ ζ 4fg g+η f (16) whose ( convergent ) condition is Re(ξ ± f ± g) < 0 and Re ζ 4fg < 0, after directly calculation, we obtain ξ±f±g Π b (ϕ) PA = Π b (ϕ) 0 k k h k g k exp cosh r cos ϕ 1 + sin ϕ sinh r gh + 4z cosh r sin ϕ + z sinh r sin ϕ 4 ( 1 + sin ϕ sinh r ) h + 4z cosh r sin ϕ + z sinh r sin ϕ 4 ( 1 + sin ϕ sinh r ) g sinh r cos ϕ 4 ( 1 + sinh r sin ϕ ) ( ] h + g ) h=g=0,(17) where Π b 0 is the corresponding expectation value of the parity operator for the input state with -coherent state 1], (cos ϕ 1 sinh r sin ϕ ) z sinh r sin ϕre(z ) Π b (ϕ) 0 = e (1+sinh r sin ϕ). 1 + sinh r sin ϕ (18) Because the goal of the interferometry is to estimate very small phase changes in quantum metrology, it may be interesting to expand Eq. (17) in the Taylor series around ϕ = 0. To write out the explicit Taylor expansion of Eq. (17) for general k in the limit ϕ 0 is a difficult task. However, for small k, the explicit Taylor expansions can be accessible. When k = 0, the Taylor expansions of Π b (ϕ) 0 is Π b (ϕ 0) 0 = 1 n z cos θ n a + n a + 1ϕ n z n a + n z + n a ϕ + O ( ϕ 4)], (19) where n z = z (z = z e iθ ) the average photon number of the coherent state, n a = ( k+1 / k 1) k=0 = sinh r the average photon number of the. While k = 1,, we can derive the analytical Taylor expansions of Π b (ϕ) PA as well, i.e., and Π b (ϕ 0) PA k=1 = ( 1) 1 n z cos θ n a + n a ϕ n z n a + n z + n a ϕ + O ( ϕ 4)], (0) Π b (ϕ 0) PA k= = 1 n z cos θ n a + n a 8 ( n a + 1 ) 3 ( n a 1 ) ϕ n z n a + n z + n a ϕ + O ( ϕ 4)], (1) where n a = k+1 / k 1 the average photon number of the PA. In order to obtain the good phase uncertainty, in the following we set θ = 0 (the phase of the coherent state) 1]. These analytical Taylor expansions of Π b (ϕ) PA mentioned above is useful for the following discussions about the parity detection. Different from that calculation method in Ref.4], for the normalized PS, using our method, we can also obtain the corresponding expectation value of the parity operator for the input state with PS-coherent states Π b (ϕ) PS = Π b (ϕ) 0 C l l h l g l exp sinh r cos ϕ 1 + sin ϕ sinh r gh z sinh r sin ϕ + z sinh r sin ϕ ( 1 + sin ϕ sinh r ) g z sinh r sin ϕ + z sinh r sin ϕ ( 1 + sin ϕ sinh r ) h sinh r cosh r ( 1 + sin ϕ sinh r ) ( ] h + g ) h=g=0 (), Although Eq. () is different from Eq. (16) in Ref.4], one can prove that two results are completely identical by numerical method. In addition, when k = l = 1, one can obtain Π b (ϕ) PS l=1 = Π b (ϕ) PA k=1, this is because

5 5 one photon-subtracted is identical to one photonadded 33]. When l =, the Taylor expansion of Π b (ϕ) PS in the limit ϕ 0 is Π b (ϕ) PS l= = 1 n z cos θ n s + n s 8 ( n s 1 ) 3 ( n s + 1 ) ϕ n z n s + n z + n s ϕ + O ( ϕ 4)] (3) where n s = C l+1 /C l the average photon number of the PS. By Eqs. (17) and (), we can investigate the expectation value of the parity operator as function of the phase shift ϕ. Given r = 0.3 and z =, in Fig. (a) we draft this expectation value against ϕ for k = 0, 1,, 3 and l = 0, 1,, 3. One can see that, for the central peak or trough of the Π b (ϕ), the narrowness of the maxima or minima at ϕ = 0 increases as k and l increase. In addition, we see that the central peak or trough of the Π b PA of the PA-coherent input states is narrower than that of the PS-coherent input states when k = l > 1.This result indicates that the photon addition can enhance supper-resolution better than photon subtraction for given initial coherent state amplitude and squeezing parameters. However, if given the same average photon number ( n a = n s ) of both the PA and the PS, the distributions of these central peaks or troughs of the Π b (ϕ) near ϕ = 0 are almost identical as shown in Fig. (b), which is also consistent with Eqs. (19-1) and (3). B. Phase sensitivity via the error propagation method The phase uncertainty ϕ is the main aspect of quantum optimal interferometry. The smallest phase uncertainty ϕ is the characteristic of the most sensitive measure. The smaller the value of Π b is, the higher the phase sensitivity is. From the error propagation method, the phase uncertainty ϕ of an interferometer can be determined as 1 Π b (ϕ) ϕ = (4) Π b (ϕ) / ϕ where we have used Π b = Π b (ϕ) Π b (ϕ) and the fact that Π b (ϕ) = 1. The phase sensitivity with parity detection for the interferometry with PAcoherent state is found to be best at ϕ = 0. Although it is difficult to write out the general explicit form of Eq. (4) when the PA-coherent state with the general value of k is considered as an interferometer state, based on Eqs. (19-1), the explicitly forms of ϕ for small k in the limit of ϕ 0 can be easily obtained. For example, P P PA k=1 k=0 k=3 r=0.3, z =4 k= l= l=0 l=3 PS PA l= n a =4, z =4 n s =4, z =4 f FIG.. The expectation value of the parity operator versus the phase shift ϕ for different PA-coherent and PScoherent input states of the MZI, respectively. according to Eq. (19), the explicit form of Eq. (4) in the case of ϕ 0 can be directly obtained 1], ( ϕ) PA k=0 = f 1 n z n a + n a n z n a + n z + n a, (5) where they have set θ = 0 (the phase of the coherent state) in order to obtain the good phase uncertainty. In addition, based on our Eqs. (0) and (1), we can also obtain the explicit form of Eq. (4) for k = 1,, respectively, ( ϕ) PA k=1 = and = 1 n z n a + n a + n z n a + n z + n a, (6) ( ϕ) PA k= 1 (7). n z n a + n a 8( 1+1 n a+1) 3( + n 1+1 n a 1) z n a + n z + n a

6 6 For the item 8( 1+1 n a+1) 3( 1+1 n a 1) in Eq. (7), one can easily obtain 8( 1+1 n a+1) 3( 1+1 n a 1) k= = 4 cosh4 r (3 cosh r 1), and prove that the inequality 8 3 8( 1+1 n a+1) 3( 6 for any values of 1+1 n a 1) the squeezing parameter r is satisfied. On the other hand, when the PS-coherent state with l = is considered as the interferometer state, according to Eq. (3), we obtain the corresponding phase uncertainty in the case of ϕ 0 as these corresponding phase uncertainties are indeed very small in the limit of ϕ 0. While the phase shift ϕ deviates from zero, the optimal state is neither the nor the PS, but the PA as shown in Fig. 3(b). On the other hand, for a given initial squeezing parameter r or given the same ratio of n a / n z (and n s / n z ) of the two input ports of the MZI, the PA has also the better performance in quantum precision measurement than both the PS and the when the phase shift ϕ slightly deviates from zero as shown in Fig. 4. ( ϕ) PS l= 1 = (8) n z n s + n s 8( 1+1 n s 1) 3( + n 1+1 n s+1) z n s + n z + n s n a = n s = 4 n z = 4 =0.15 where we have set θ = 0 as well. For the item 8( 1+1 n s 1) 3( 1+1 n s+1) l= = 4 sinh4 r in Eq. (8), one can (3 sinh r+1) easily prove that the inequality 0 4 sinh 4 r(3 sinh r + 1) 8 3 is satisfied. Obviously, Eqs. (5-8) indicate that, within constraints on the total average photon number ( = n a + n z or = n s + n z ) and the same ratios of n a / n z or n s / n z, the best phase uncertainty is obtained by the when one input port of a MZI is injected by a coherent state. When k = l = 0, both PA and PS reduce to a. Therefore, our results obtained via parity detection indeed support Lang and Caves s work that the is the optimal state for a MZI in the limit of ϕ 0 when one input is a coherent state 3]. On the other hand, given a larger value of the same average photon numbers of the PA, the PS and the, we can also see from Eqs. (5-8) that the difference among these corresponding phase uncertainties is very small. With the increases of the values of the k and l, we can numerically prove such conclusion is still true. As pointed in Refs. 0, 3], we can also prove that, in the case of the n a = n z (or n s = n z ), the best phase sensitivity will be obtained and reach the. However, when the phase shift ϕ to be estimated slightly deviates from zero, the with many photons may be not the optimal state for a MZI. Firstly, we investigate how the phase uncertainty varies with ϕ. Based on Eq. (4), we plot the phase uncertainty as a function of the phase shift ϕ in Fig. 3. Fig. 3(a) shows that, given somewhat small average photon numbers n a = n s = n z = 4, the MZI with one input port injected by a coherent state, the state is the optimal state for such MZI. Yet, interferometry with coherent states contained a few average number of photons is meaningless since such a measurement cannot be precise in principle. Moreover as pointed in the conclusion part in Ref. 3], it may interesting to explore how the phase sensitivity changes when the carry many photons. Thus, with the increases of the value of n a (or n s ) and n z, we can see from Fig. 3(b) that the differences among k=3 PA k= k=1 l=3 l= l=1 PS =15 =060 n a = n s = 16 j n z = FIG. 3. The phase uncertainty ϕ as a function of the phase shift ϕ when the PA-coherent state, the PS-coherent state and the -coherent state as interferometer states, respectively. (a) n a = n s = n z = 4; (b) n a = n s = n z = 16. The horizontal dashed line denotes the limit. Secondly, for given a constraint on the total average photon number of the PA, the PS and the, j

7 r = 0.9 n z = 100 ture of the expectation value of the parity operator, but there are other photon numbers where the uncertainty is still below the SL for both the PA and the PS. n a = n s = 4 SL PA k=1 k= k=3 PS l=1 l= l= j n a = n s = 5 SL n a = n s = 16 n z = j n a = n s = 49 SL c FIG. 4. The phase uncertainty ϕ as a function of the phase shift ϕ when the PA-coherent state, the PS-coherent state and the -coherent state as interferometer states. (a) r = 0.9, n z = 100; (b) n a = n s = 16, n z = 100. it is important to investigate how the phase sensitivities change with the total average photon number. In Fig. 5, for the PA (k = 0, 1,, 3, 6) and the PS (l = 0, 1,, 3, 6), we give some same values of both n a and n s and plot the phase sensitivities versus the total number of photons. Obviously, from Fig. 5 we can see that almost the same phase uncertainties can be obtained in the limit of ϕ 0. In addition, when the values of both n a and n s increase, these phase uncertainties more quickly approach to the. Of course, when n a = n z (or n s = n z ), the optimal phase uncertainties will be obtained. In Fig. 6, we repeat these graphs for ϕ = 5. In the latter case especially, we see that the optimal state is the PA. When n a and n s increase, the phase uncertainties blow up due to the periodic na FIG. 5. Plots of the phase uncertainty versus total average photon at ϕ = 10 4 for fixed the average photon number of the PA and the PS, as well as the, along with the corresponding curves for the SL and the limits. Only the parameter is being changed. Finally, we further investigate how the photon addition and subtraction affect the phase sensitivity for a given initial squeezing parameter of the. Comparing with that results in Ref. 4], we find that, for the same initial squeezing parameter r, it is also better to perform photon addition rather than photon subtraction

8 8 in respect to reducing the phase uncertainty. This is because that it is always better to perform photon addition rather than photon subtraction in order to increase the average photon number for given the initial squeezing as shown in Fig. 1. IV. QUATUM FISHER IFORMATIO OF THE MZI ITERFEROMETER In this section, we will prove that the parity detection is the optimal measurement for our considered interferometric scheme. ow, we use the quantum Fisher information to find the maximum level of the phase sensitivity by the Cramér Rao bound as given by 17] ϕ min = 1 FQ, (9) For pure states injected into a MZI, the quantum Fisher information F Q can be obtained by 38] F Q = 4 ψ (ϕ) ψ (ϕ) ψ (ϕ) ψ (ϕ) ], (30) where ψ (ϕ) = e iϕj3 e iπj1/ ψ in is the state just before the second beam splitter of the MZI, and ψ (ϕ) = ψ (ϕ) / ϕ. In terms of the input state, the quantum Fisher information becomes F Q = 4 in J in in J in ], (31) and thus the quantum Fisher information is, up to factor of 4, the variance of the operator J. In the present work, the PA-coherent state as the interferometer state. In order to obtain the quantum Fisher information, based on Eqs. () and (3), we first obtain the following expectation values, i.e., b PA = b = 0, and PA k+,k b PA as well as = k+ t k+ τ k e sinh r 4 (t +τ )+tτ cosh r t,τ=0, (3) k,k+ b PA = k+ t k τ k+ e sinh r 4 (t +τ )+tτ cosh r t,τ=0.(33) Then, according to Eqs. (9) and (31), we can directly obtain the Quantum Fisher information of the MZI as F QA = n z n a + n z + n a n z k+,k k, (34) where we have also set θ = 0 (the phase of the coherent state) for obtaining the good phase uncertainty and used the relation k+,k = k,k+. For k = 0, 1,, combining Eqs. (5-8) and Eq. (34), we can analytically prove that the quantum Cramér-Rao bound can be reached via the parity detection in the limit ϕ 0. For general k, we can numerically prove this is still true. Using the similar method, we derive the quantum Fisher information for the PS-coherent interferometer state where F QS = n z n s + n z + n s n z C l+,l C l, (35) C l+,l b PS = k+ t k+ τg k e sinh r 4 (t +τ )+tτ sinh r t,τ=0, C l,l+ b PS = k+ t k τ k+ e sinh r 4 (t +τ )+tτ sinh r t,τ=0, (36) and we have used the relation C l+,l = C l,l+. Similarly, we can check that the quantum Cramér-Rao bound can be reached via the parity detection in the limit ϕ 0 for the PS-coherent considered as the interferometer state. V. COCLUSIOS In summary, we have studied the quantum optimal interference by mixing a coherent state with a PA. Given a constraint on the total average number of photons and in the limit of ϕ 0, for a MZI with PAcoherent and PS-coherent as well as -coherent input states, almost the same phase uncertainties can be obtained. However, when the phase shift ϕ somewhat deviates from zero, the optimal state is neither the nor the PS, but the PA when these three states contain many photons. On the other hand, for fixed the initial squeezing r, it is better to perform photon addition rather than photon subtraction for improving the phase sensitivity of the MZI. This may be because that it is always better to perform addition rather than subtraction in order to increase the average photon number of the for given the initial squeezing. Finally, we show that the quantum Cramér-Rao bound can be reached via the parity detection in the limit ϕ 0. ACKOWLEDGMETS This work is supported by the ational atural Science Foundation of China (o , o , and o ) and sponsored by Qing Lan Project of the Higher Educations of Jiangsu Province of China.

9 9 1] V. Giovannetti, S. Lloyd, and L. Maccone, Advances in quantum metrology, at. Phot. 5, 9 (011). ] L. Pezzè, A Smerzi, M.K. Oberthaler, R. Schmied, and P. Treutlein, Quantum metrology with nonclassical states of atomic ensembles, arxiv: v3 (018). 3] C.M. Caves, Quantum-mechanical noise in an interferometer, Phys. Rev. D 3, (1981). 4] Z.Y. Ou, Fundamental quantum limit in precision phase measurement, Phys. Rev. A 55, (1997). 5] J.J. Bollinger, W.M. Itano, D.J. Wineland, and D.J. Heinzen, Optimal frequency measurements with maximally correlated states, Phys. Rev. A 54, R4649 R465 (1996). 6] J. Dunningham and T. Kim, Using quantum interferometers to make measurements at the Heisenberg limit, J. Mod. Opt. 53, (006). 7] G. Gilbert, M. Hamrick, and Y.S. Weinstein, Use of maximally entangled -photon states for practical quantum interferometry, J. Opt. Soc. Am. B 5, (008). 8] B.L. Higgins, D.W. Berry, S.D. Bartlett, H.M. Wisman, and G.J. Pryde, Entanglement-free Heisenberg-limited phase estimation, at. Lett. 450, (007). 9] H. Lee, P. Kok, and J.P. Dowling, A quantum Rosetta stone for interferometry, J. Mod. Opt. 49, (010). 10] C.C. Gerry, A. Benmoussa, and R.A. Campos, onlinear interferometer as a resource for maximally entangled photonic states: Application to interferometry, Phys. Rev. A 66, (00). 11] M.W. Mitchell, J.S. Lundeen, and A.M. Steinberg, Super-resolving phase measurements with a multiphoton entangled state, ature 49, (004). 1] G.Y. Xiang, B.L. Higgins, D.W. Berry, H.M. Wiseman, and G.J. Pryde, Entanglement-enhanced measurement of a completely unknown optical phase, at. Phot. 5, (011). 13] R. Demkowicz-Dobrzanski, M. Jarzyna, and J. Kolodynski, Quantum Limits in Optical Interferometry, Prog. Opt. 60, (015). 14] Q.S. Tan, J.Q. Liao, X.G. Wang, and F. ori, Enhanced interferometry using squeezed thermal states and even or odd states, Phys. Rev. A 89, 0538 (014). 15] S.Y. Lee, C.W. Lee, J. Lee, and H ha, Quantum phase estimation using path-symmetric entangled states, Sci. Rep. 6, (016). 16] M.J. Holland and K. Burnett, Interferometric detection of optical phase shifts at the Heisenberg limit, Phys. Rev. Lett. 71, (1993). 17] Z.Y. Ou, Complementarity and Fundamental Limit in Precision Phase Measurement, Phys. Rev. Lett. 77, (1996). 18] L. Pezzè and A. Smerzi, Ultrasensitive Two-Mode Interferometry with Single-Mode umber Squeezing, Phys. Rev. Lett. 110, (013). 19] S. Wang, Y.T. Wang, L.J Zhai, and L.J. Zhang, Twomode quantum interferometry with a single modefock state and parity detection, J. Opt. Soc. Am. B, 35, (018). 0] L. Pezzé and A. Smerzi, Mach-Zehnder Interferometry at the Heisenberg Limit with Coherent and Squeezed- Vacuum Light, Phys. Rev. Lett. 100, (008). 1] K.P. Seshadreesan, P.M. Anisimov, H. Lee and J.P. Dowling, Parity detection achieves the Heisenberg limit in interferometry with coherent mixed with squeezed vacuum light, ew J. Phys. 13, (011). ] C.C. Gerry and J. Mimih, The parity operator in quantum optical metrology, Contemp. Phys. 51, (010). 3] M.D. Lang and C.M. Caves, Optimal Quantum- Enhanced Interferometry Using a Laser Power Source, Phys. Rev. Lett. 111, (013). 4] R. Birrittella and C.C. Gerry, Quantum optical interferometry via the mixing of coherent and photon-subtracted squeezed vacuum states of light, J. Opt. Soc. Am. B 31, (014). 5] J. Wenger, R. Tualle-Brouri, and P. Grangier, on- Gaussian Statistics from Individual Pulses of Squeezed Light, Phys. Rev. Lett. 9, (004). 6] A. Ourjoumtsev, A. Dantan, R. Tualle-Brouri, and P. Grangier, Increasing Entanglement between Gaussian States by Coherent Photon Subtraction, Phys. Rev. Lett. 98, (007). 7] G.S. Agarwal and K. Tara, onclassical properties of states generated by the excitations on a coherent state, Phys. Rev. A 43, (1991). 8] A. Zavatta, S. Viciani, and M. Bellini, Quantum-to- Classical Transition with Single-Photon-Added Coherent States of Light, Science 306, (004). 9] C.W. Helstrom, Quantum Detection and Estimation Theory (Academic Press, ew York, 1976). 30] M.O. Scully and M.S. Zubairy, Quantum Optics (Cambridge University Press, ew York,1997). 31] S. Wang, H.C. Yuan, and X.X. Xu, Conditional generation of Fock states and Schrödinger-cat states via adding multiple photons to a squeezed vacuum state, Eur. Phys. J. D 67, 10 (013). 3] X.G. Meng, Z. Wang, H.Y. Fan, and J.S. Wang, onclassicality and decoherence of photon-subtracted squeezed vacuum states, J. Opt. Soc. Am. B 9, (01). 33] M.S. Kim, Recent developments in photon-level operations on travelling light fields, J. Phys. B, 41, (008). 34] M.B. Stephen, F. Gergely, R.G. Claire, and C.S. Fiona, Statistics of photon-subtracted and photonadded states, Phys. Rev. A, 98, (018). 35] See, for example C. C. Gerry and P. L. Knight, Introductory Quantum Optics (Cambridge University Press, ew York, 005), Chap ] B. Yurke, S.L. McCall, and J.R. Klauder, SU() and SU(1,1) interferometers, Phys. Rev. A 33, (1986). 37] L. Cohen, D. Istrati, L. Dovrat, and H.S. Eisenber, Super-resolved phase measurements at the shot noise limit by parity measurement, Opt. Express, (014). 38] Y. Ben-Aryeh, Phase estimation by photon counting measurements in the output of a linear Mach-Zehnder interferometer, J. Opt. Soc. Am. B 9, (01).

10 10 j=5 n a = n s = 4 SL j=5 n a = n s = 5 SL j=5 n a = n s = 49 SL c FIG. 6. Plots of the phase uncertainty versus the total average photon at ϕ = 5 for fixed the average photon number of the PA and the PS, as well as the, along with the corresponding curves for the SL and the limits. The thick red line represents the PA(or PS)-coherent state with k = l = 1; the orange line represents the PAcoherent state with k = ; the orange dashed line denotes the PA-coherent state with k = 3; the orange dotteddashed line denotes k = 6. While, the blue line represents the PS-coherent state with l = ; the blue dashed line represents l = 3; the blue dotted-dashed line denotes l = 6.

11 k=l=0 n a = n s = 1.05 s Ö s Ö SL PA and PS k=l=1 n a = n s = k=l=0 n a = n s = 1.05 SL PS PA l= n s = 6.8 k= s n a = 7.68 Ö s Ö H k=l=0 n a = n s = 1.05 PS SL c PASV l=3 n s = 8.89 k=3 n s a = 11.5 Ö s Ö FIG. 7. For r = 0.9, plots of the phase uncertainty versus the total average photon number at ϕ = 10 4, along with the corresponding curves for the SL and limits. Only the parameter is being changed.