Phase estimation of definite photon number states by using quantum circuits

We propose a method to map the conventional optical interferometry setup into quantum circuits. The unknown phase shift inside a Mach–Zehnder interferometer in the presence of photon loss is estimated by simulating the quantum circuits. For this aim, we use the Bayesian approach in which the likelihood functions are needed, and they are obtained by simulating the appropriate quantum circuits. The precision of four different definite photon-number states of light, which all possess six photons, is compared. The measurement scheme that we have considered is counting the number of photons detected after the final beam splitter of the interferometer, and photon loss is modeled by using fictitious beam splitters in the arms of the interferometer. Our results indicate that three of the four definite photon-number states considered can have better precision than the standard interferometry limit whenever the photon loss rate is in a specific range. In addition, the Fisher information for the four definite photon-number states in the setup is also estimated to check the optimality of the chosen measurement scheme.


I. INTRODUCTION
Precisely estimating an unknown parameter is important both in science and engineering.Depending on the parameter that one wants to estimate, one can use different approaches.In some cases, they may be able to measure the quantity directly, however, sometimes quantities are not directly accessible either in principle or due to experimental impediments, and they need to be estimated indirectly.
Different estimation approaches have different precisions in estimating the unknown quantity.Generally speaking, these approaches use the resources of atoms and/or photons to estimate the desired quantity, and by increasing the resources, they can reveal the unknown quantity with better precision.Nevertheless, increasing the resources as much as one wants is not possible in some cases (for example biological samples [1]), so it is essential to design an apparatus that gives the best possible estimation of the quantity by using a specified number of resources.
The optical phase is the most significant parameter in the context of quantum metrology, especially in the optical interferometry [2][3][4].Optical interferometers use photons as resources and are one of the most precise tools for measuring quantities such as length, index of refraction, velocity, etc [1,5], by mapping them into a phase shift.By using classical states of light, the best precision that one can estimate the phase shift scales like 1/ √ N [6], in which N is the number of photons used.This limit is called the standard quantum limit (SQL).However, this is not the ultimate precision limit; we can enhance the precision by a factor of 1/ √ N and reach the Heisenberg limit (1/N ) by employing quantum mechanical strategies like using entangled photonic states.
In a typical phase estimation setup, the input state together with the strategy that we use to extract the information and estimate the phase shift, determine the precision of the estimation [7][8][9][10][11][12][13][14][15][16][17][18].Choosing the best input state depends on the decoherence existing in the interferometer [6,19] and also on the number of resources (for example number of photons) that we can use.Apart from the above factors, our prior knowledge (before making any measurement) about the unknown quantity can also affect the optimal estimation strategy [20].
Consider a two-mode interferometer like Mach-Zehnder for instance.If prior to measurement we know that the unknown phase is around a specific value with slight deviation and the number of photons we can use is N , and there is no photon loss and decoherence in the interferometer, N OON states give the best possible precision and achieve the Heisenberg limit [21,22].However, when there is photon loss in the arms of the interferometer, N OON states become extremely fragile [21,[23][24][25][26][27] and there exist other states which outperform N OON states in different photon loss conditions.
A class of path-entangled photon Fock states, called mm states, were proposed in [28] to outperform N OON states in the presence of photon loss.Huver et al. have introduced a theoretically optimal detection operator and computed the precision of their phase estimation scheme by calculating the expectation value of the detection operator and using the error propagation formula.However, experimental realization of their optimal measurement scheme has not been reported yet.Performance of N OON and mm states by using parity detection, has also been investigated in the presence of photon loss [29] and phase fluctuation [30].Another class of definite photon number states that act better than the N OON states in the presence of loss are Holland-Burnett (HB) states [31].Unlike N OON states, HB states can be prepared easily by using a beam splitter and photon pairs from spontaneous parametric down-conversion (SPDC) [32].HB states are almost optimally robust not only to photon loss but also to inefficiencies in state preparation and photon detection and can beat SQL in realistic phase estimation scenarios [32][33][34][35][36][37].
In this work, we focus on four different definite photon number input states of light, which all possess six photons and belong to three different classes of non-classical states: N OON states, mm states, and HB states.We will compare their phase estimation precision in different photon loss conditions in a Mach-Zehnder interferometer by mapping the corresponding quantum optical setups into quantum circuits and simulating them to obtain the information needed.This is a new approach towards investigating phase estimation precision of definite photon number input states in the presence of loss.
We will choose a 50-50 beam splitter with two photon number-resolving detectors (PNRD) as our measurement scheme and use Bayesian estimation to estimate the unknown phase shift.Likelihood functions that we need in the Bayesian estimation for different input states will also be approximated by simulating appropriate quantum circuits rather than calculating them theoretically, which can be arduous in the presence of photon loss.
In detail, we proceed as follows: Section II explains our model and shows how to map the corresponding quantum optical setups into quantum circuits.In section III our results are presented and analyzed.Finally, a summary and conclusion are given in section IV.

II. DESCRIBING OUR MODEL AND DESIGNING CORRESPONDING QUANTUM CIRCUITS
The phase estimation setup considered here is similar to the ordinary Mach-Zehnder interferometer in which the first beam splitter is replaced by a black box that can generate any desired input photonic state.To model the photon loss, two extra fictitious beam splitters are used in both arms of the interferometer [28,29,[38][39][40].Different rates of photon loss can simply be created in the interferometer by adjusting the transmission coefficients of the beam splitters.To simulate the optical interferometry setup, an appropriate quantum circuit is designed.At first, we map the input photonic state into a corresponding state in qubits.There are two modes for input photons: The lower and the upper paths are the first and the second modes, respectively.So an arbitrary pure input state with a definite photon number (N ) would be like where 1) is in mode representation and |m l |N − m u means that m out of N photons are in the lower arm and the remaining N −m photons are in the upper arm.The input state is a superposition of states like |m l |N − m u , so if we can map these states into qubits of our circuit appropriately, we can represent the input state by the superposition of these corresponding qubit states.
By writing |m l |N − m u in the particle representation (and applying bosonic symmetry) we arrive at [41]: (2) Here, |l j (|u j ) means that the jth photon is in the lower (upper) arm and is an arbitrary permutation of 'l's and 'u's, and summation is over nontrivial permutations.Eq. ( 2) clarifies how to show the input state in qubits.We choose N qubits and label them from 1 to N by assuming that each qubit corresponds to a photon in Eq. ( 2).Then we create the bosonic symmetric state of these qubits.In each term of the total superposition state (Eq.( 2)), if the photon is in the lower arm, we put the analogous qubit in |0 , and if the photon is in the upper arm, we put the qubit in |1 .The interferometer's two arms are distinguishable, which means that considering the two orthogonal states (computational bases) in this mapping is legitimate.By this map, we can only devote two distinguishable modes for photons at each stage of the circuit, which can only show the case where there is no loss.When we introduce photon loss in both arms, we need to have two other distinguishable modes to each photon so that they can be in four distinguishable modes: lower arm, upper arm, lost from the lower arm, and lost from the upper arm.Therefore, we change the mapping, and two qubits for each photon are allocated.Our map in the presence of the photon loss is as follows: As a result, if the input photonic state has N photons, we allocate 2N qubits to them in our quantum circuit.
Figure .2 shows the design of our quantum circuit.We assume that all qubits are initially at |0 state.U input is a unitary evolution that transforms |0 1 |0 2 ... |0 2N to our desired input state.Considering that our input state is like Eq. (1) in photon representation, the analogous state of the qubits would be where arbitrarily permutes the corresponding qubit pairs' states of different photons and as in Eq. ( 2), means summation over nontrivial permutations.As mentioned earlier, U input transforms initial state |0 1 ... |0 2N to the state of Eq. (3).To implement the U input , we have to decompose it into products of gates that are available (we call them basis gates here) to us.Appendix (A) illustrates the list of the gates that we have used in our quantum circuits, and appendix (B) demonstrates how to do this decomposition.After preparing the input state, the phase shift, fictitious beam splitters, and the final beam splitter of the interferometer are implemented in order.As it is clear from Fig. 2, each of the above three parts (phase shift, fictitious beam splitters, and the final beam splitter) consist of N similar subroutines.Each particular subroutine acts only on the two qubits corresponding to a specific photon.Subroutines of the three parts transform the state of a qubit pair corresponding to a photon as below [42] In Eq. ( 4), φ is the phase shift in the interferometer, and t 0 and t 1 are the transmissivity of the lower and upper fictitious beam splitter, respectively.L 0 = 1 − t 0 and L 1 = 1 − t 1 determine the photon loss rate in the lower and upper arms respectively.We interpret the state after the final beam splitter in a way that |0 |0 and |1 |1 mean photon is in the path going to D 0 and D 1 respectively.Appendix (C) shows the sequence of the basis gates to implement these transformations.
At the end of the circuit, we measure the qubits in the computational basis and end up with a string of '0's and '1's.Finally, we convert this string to a result that shows us the number of photons detected by D 0 and D 1 .
Now our quantum circuit simulates the analogous quantum optical interferometer, and we can compare different input states' precision in estimating phase shift in different photon loss rates.First, we implement the appropriate U input to create our desired input state, then we put φ = φ * (which is the value of the quantity that the optical interferometer is going to estimate).Then we adjust t 0 and t 1 to create the photon loss rates that we want.Finally, we simulate the circuit for N r times, which means that the input state is sent through the interferometer for N r times to find a probability distribution P (φ) for phase.The variance of this probability distribution shows the error in estimating the phase shift; the lower the variance, the better the precision.
To estimate the properties of an object in a real experiment, it is placed in one of the arms of the interferometer and then the phase shift (induced by the object) is estimated by sending the input states (photons) through the interferometer for a specified number of times.The phase shift's actual value (φ * ) is unknown to the experimenter, but the number of photons that get detected at different detectors each time is known, and the phase shift can be estimated by this information.In this work, it is assumed that before sending the input states through the interferometer, the experimenter initially knows that the unknown phase shift is limited to a finite range, We show the measurement result of each run of the experiment by D(n 0 , n 1 ), which means that n 0 and n 1 photons are detected at D 0 and D 1 respectively; note that dφ P (D(n0,n1)|φ ) for different possible values of n 0 and n 1 , the probability distribution could be estimated by the following algorithm [43,44] 1. P (φ) = P pri (φ) 2. run the experiment and multiply P (φ) by P n0,n1 (φ) 3. P (φ) = normalized(P (φ) * P n0,n1 (φ)) 4. repeat steps 2 and 3 for N r − 1 times The output P (φ) of this algorithm is the probability distribution of the unknown phase shift (φ * ) after sending the input photonic state for N r times, and is a function of the N r measurement results ( D).We refer to it by P f inal (φ| D).We can estimate the φ * by the following formula: The standard deviation is given by ∆ φest is a function of the specific measurement results, D, that has been obtained by doing the experiment for N r times [45].Repeating the same experiment may yield different measurement results D , which could result in a slightly different ∆ φest .Therefore, in order to remove the dependency on specific D, we will calculate the average standard deviation Now we focus on how to find P (D(n 0 , n 1 )|φ) for each input state in specified condition of photon loss.We can find it by calculating it analytically or using the quantum circuit that we have designed to approximate it.Calculating P (D(n 0 , n 1 )|φ) analytically can sometimes be difficult, especially when there is photon loss, but our method avoids this difficulty and approximate the probability distribution up to the precision that we want.To do that, we prepare the input state, fix the values of t 0 and t 1 , which creates the desired photon loss condition, and choose n equally spaced points from φ i to φ f .Now by running the circuit for W times (W 1) at each of these n points, we can approximate P (D(n 0 , n 1 )|φ) by the following formula: where R(n 0 , n 1 ) is the number of times that n 0 and n 1 photons are detected at lower and upper detectors respectively.By increasing the number of points (n) and the number of times that we simulate the circuit at each point (W ), we can approach the exact P (D(n 0 , n 1 )|φ).
As mentioned in section I, we consider definite photon number input states of three different classes, which all possess the same number of photons.The states of N OON and mm classes for N photons are written as [29]: Here m and m are positive integers and without loss of generality we assume that m > m .m + m = N is the total number of photons that |m :: m state possesses.
The last input state that is investigated here is called the Holland-Burnett (HB) state, which has an even number of N photons and is given by the following formula [31,32] (before experiencing the phase shift): In our simulations we take N = 6 and the four input states considered are: |6 :: 0 , |5 :: 1 , |4 :: 2 and |HB(6) .

III. RESULTS
In this work, two different cases for photon loss are considered.Subsection III A shows the results for the case of symmetric photon loss where the probability of photon loss in the lower and upper arms of the interferometer are equal (L 0 = L 1 = L, t 0 = t 1 = t), and subsection III B shows the results for the asymmetric photon loss case where there is no photon loss in the lower arm of the interferometer (L 0 = 0).In our simulation of estimating the unknown phase shift, it is assumed that φ ∈ [φ i , φ f ], where φ f − φ i = π 6 ; this is done to avoid the ambiguity in estimating phase.The actual value of the unknown phase is , and it is chosen in a way that minimizes δφ for each input state.For |6 :: 0 , |5 :: 1 and |4 :: 2 , φ * is π/12, π/8 and π/4, respectively [30].For |HB (6) , approximate value of the φ * is found by simulating the quantum circuits with |HB(6) as input state for different values of phase shifts, and comparing their estimated average standard deviations.

A. Symmetric photon dissipation
Table.I shows the unnormalized likelihood functions (P n0,n1 (φ)) of the input states for t = t 0 = t 1 = 0.5.In TABLE I. Unnormalized likelihood functions for different input states in symmetric photon loss when t0 = t1 = 0.5.The vertical axis in each row shows Pn 0 ,n 1 (φ)/ max(Pn 0 ,n 1 (φ)) for the measurement result D(n0, n1), which is shown in the left side of each row.Each column in the table corresponds to one of the input states, which is shown at the top of the column.By moving down the table, the total number of detected photons (n0 + n1) generally decreases.
|6 this symmetric case, P n0,n1 (φ) for other values of t is similar; note that when t = 1, always the number of photons detected is 6, therefore P n0,n1 (φ) = 0 for n 0 + n 1 < 6.These probability distributions are calculated by simulating the corresponding quantum circuit at different phases in [0, 2π] for a large number of times and fitting an appropriate periodic function on them.They are plotted in a way that their maximum value is 1 (except P 3,3 (φ) for |5 :: 1 which is zero because this combination of photons never gets detected).They are not normalized, but this is not an issue because we normalize the final probability distribution at each step when we are estimating the unknown phase.Table .I is ordered in a way that the total number of detected photons (n 0 +n 1 ) generally decreases as we go below the table.As we can see, whenever the number of detected photons is less than 6 (meaning that at least one photon is lost during the process), the probability distribution for |6 :: 0 is a constant function, which means that we get no information about the phase shift in this case.Therefore, by increasing the photon loss probability, the precision of phase estimation with the N OON states will decrease more rapidly than the other states.In asymmetric photon loss, our simulation results show that likelihood functions change by altering the value of t 1 , and therefore we had to estimate them for different values of t 1 in [0, 1].Because of this, we have not shown them here.Fig. 4 Shows the average standard deviation of the input states in estimating the phase shift for N r = 400 in the asymmetric photon loss case.In this case when t 1 > 0.79, |6 :: 0 gives the best precision while for other values of t 1 , |HB(6) has the best precision.Unlike the symmetric loss case, in asymmetric case, |4 :: 2 never becomes the state with the lowest (∆φ) avg in any loss rate.
Comparing fisher information of our measurement scheme (F ρ φ , Mc ) with QFI (FQ(ρ φ )) for |6 :: 0 in symmetric (left plot) and asymmetric (middle plot) loss, and for |HB(6) in asymmetric loss (right plot).It can be seen that for the N OON state (|6 :: 0 ) our measurement scheme ( Mc) is optimal.However, for |HB(6) , F ρ φ , Mc at most points is a little lower than the QFI, which means that Mc is not the optimal measurement scheme at these points, but it is very close to the optimum value.

C. Estimating fisher information
Generally speaking, a phase estimation process can be divided into three steps: In the first step, one prepares a state (ρ φ ), which contains information about the unknown phase that they want to estimate.In the second step, a measurement scheme ( M ) is set up to extract the information of ρ φ , and in the last step, an appropriate estimator is chosen to estimate the unknown phase, based on the measurement results.The likelihood functions for the setup depends on ρ φ and M , so the fisher information of the phase estimation process depends on them can be calculated by the following formula [46,47]: where the summation is over all possible outcomes of the measurement ( D).In frequentist approach, the best possible precision that one can estimate (by using a locally unbiased estimator) the unknown phase is [5,45] ∆φ min = 1 and the estimator giving this precision is an efficient estimator.By changing the measurement scheme ( M ), one can change the fisher information of the setup, and as a result, the best possible precision in phase estimation.
Maximizing the fisher information over all generic quantum measurement schemes yields the Quantum Fisher Information (QFI) of ρ φ [5]: The measurement scheme that maximizes Fisher information is called the optimal measurement scheme for ρ φ .QFI depends only on ρ φ and can ba calculated by [5]: where L φ is the symmetric logarithmic derivative operator and (∆L φ ) 2 is its variance over ρ φ .
Regarding the interferometry setup that we are investigating (Fig. 1), ρ φ is the quantum state of light (after tracing out environment modes) before the final beam splitter, and the final beam splitter together with the PNRDs is the measurement scheme, which we show it by Mc .We estimate the fisher information for all the input states in both photon loss cases and discuss the measurement scheme's optimism.Instead of using Eq. ( 13), fisher information is estimated by using the Bayesian approach and simulating the quantum circuits a large number of times.
In the Bayesian estimation, when the number of measurements N r increases, the final probability distribution P f inal (φ| D) approaches a Gaussian distribution [45]: , (N r >> 1).
(17) This probability distribution is centered around the actual value of the unknown phase shift φ * , with the variance of 1 NrF (φ * ) in which F (φ * ) is the Fisher information of the setup.Using this result and simulating each quantum circuit a large number of times (N r >> 1), we have estimated the input states' Fisher information in symmetric and asymmetric photon loss conditions for the photon-counting measurement scheme.QFI has been calculated For N OON states in both symmetric and asymmetric photon loss [6,19,48] ) and asymmetric (the two lower figures) loss.By increasing photon loss rate, ∆φmin of parity detection surges drastically compared to the other two schemes, which means that parity detection is not a good measurement scheme for these mm states in the presence of loss.However, ∆φmin for our measurement scheme ( Mc), is close to the optimal scheme ( MA), which means that Mc is a near-optimal measurement scheme for the two mm states in the presence of loss.
asymmetric case [33].Fig. 5 shows comparison of these QFIs with the fisher information of our measurement scheme (F ρ φ , Mc ).According to this figure, for |6 :: 0 states in both loss cases, the photon-counting measurement scheme is the optimal measurement because F ρ φ , Mc is equal to QFI.For |HB(6) state in asymmetric loss, fisher information of our scheme is very close to QFI, and even at some points, it is equal to QFI; however, at most points, it is lower than the QFI which means that at these points our measurement scheme is not optimal.
To the best of our knowledge, QFI has not been calculated for mm states inside a lossy interferometer.However, the precision of phase estimation for this case has been calculated by two other measurement schemes.Bardhan et al. [30] by measuring the parity operator and using the error propagation formula have calculated the precision of phase estimation for parity detection.We show this measurement scheme by Mp here.FIG. 7. Comparing ∆φmin of our scheme with standard interferometry limit.Circle points are the values that are estimated by simulating the quantum circuits and the dashed lines are appropriate fits.The region between Heisenberg and standard classical limit is colored with light-gray.Whenever a line (corresponding to any input state) is in this region, it means that the corresponding input state has better precision than any classical interferometry setup, using the same number of resources.
According to these figures, |6 :: 0 , |5 :: 1 and |HB(6) can beat the classical limit in some loss rates in both symmetric (the left figure) and asymmetric case (the right figure); however, |4 :: 2 cannot beat the classical limit with our measurement scheme in any loss rates of symmetric and asymmetric losses.
They have shown that their operator is theoretically optimal for mm states, but they have not discussed the practical implementation of their optimal measurement scheme.We show this measurement scheme by MA .Figure. 6 shows comparison of ∆φ min of our measurement scheme with these two schemes, for |5 :: 1 and |4 :: 2 in symmetric and asymmetric losses.It reveals that, 1) when there is no photon loss, all three schemes have the same precision and ∆φ min = 1 m−m , which is the best possible precision for mm states in the absence of loss [28].2) By increasing the photon loss rate (in both symmetric and asymmetric cases), parity detection's precision decreases with a much greater slope than the other two schemes.It means that parity detection is not a good measurement scheme to extract the phase information in the presence of loss.3) The difference between the photon counting measurement scheme and the optimal scheme is tiny, especially for |5 :: 1 .Therefore the photon counting measurement scheme is near-optimal for these mm states and can be implemented in the lab by using a 50-50 beam splitter and two PNRD detectors.
At the final part, in Fig. 7, we have compared ∆φ min = 1 F ρ φ , Mc of our phase estimation setups with the Heisenberg limit, ∆φ HL = 1 N , and standard interferometry limit [6]: which is the best possible precision that a classical setup can reach.The region between the standard classical limit and the Heisenberg limit is colored with light-gray.Whenever a line is in this colored region, there is an improvement over classical interferometry.According to these results, |4 :: 2 cannot beat the classical limit in any loss rate.The other three input states can act better than the standard interferometry limit whenever the loss rate is lower than a specific value.

IV. SUMMARY AND CONCLUSION
We have simulated a lossy Mach-Zehnder interferometer by using quantum circuits, for four definite photonnumber states (|6 :: 0 , |5 :: 1 , |4 :: 2 and |HB( 6) ) as its input states.The measurement scheme that we have considered here is counting the number of photons detected at the interferometer's two different output ports.Because of the current quantum hardware's noise, the quantum circuits have been simulated by a classical computer.We have compared the phase estimation precision of the four input states when there is symmetric and asymmetric photon loss in the interferometer.According to our results, |HB(6) state has the highest precision in most of the loss rates in both symmetric and asymmetric cases; in the symmetric loss whenever 0.68 > L > 0.19, and in the asymmetric loss whenever L 1 > 0.21, |HB(6) has the highest phase estimation precision.This means that |HB( 6) is better than the others for real-world quantum optical phase estimations.We have also checked the optimality of the photon-counting measurement scheme for these input states and found that it is optimum for |6 :: 0 in both loss cases and is near-optimum for the other three states when there is loss.Our results also show the range of photon loss rates, in which one can beat the standard interferometry limit by using |6 :: 0 , |5 :: 1 and |HB(6) as input states in the Mach-Zehnder interferometer and photon-counting as measurement scheme.
By introducing appropriate subroutines into our quantum circuits, we can also implement phase fluctuations, input state preparation inefficiencies, and detector inefficiency.Furthermore, we can also implement any other measurement scheme by adding the appropriate subroutine in the last part of our quantum circuits.These features help us compare different definite photon-number input states' precision in different decoherence conditions with different measurement schemes.We keep doing this until we fix the final projection coefficient (c n ), and afterward, our input state is ready.Figure 8 shows how these three operations are implemented in our circuit for qubit{1} and qubit{2}, which correspond to photon{1}.In implementing phase shift, we have used the fact that before the fictitious beam splitters, qubits corresponding to a photon can only be in state |0 |0 or |1 |1 .These gates are implemented to other qubits (related to other photons) similarly.

Figure. 1
shows the corresponding optical interferometry setup in which the detectors (D 0 and D 1 ) are photon-number-resolving detectors (PNRD).

FIG. 1 .
FIG. 1.The optical interferometry setup.Black box prepares the desired input state.φ indicates the phase shift between two arms and photon loss is modeled by adding a fictitious beam splitter (red box) in each arm.The blue box represents a 50-50 beam splitter and D0 and D1 are photon-numberresolving detectors.

FIG. 2 .
FIG.2.Our quantum circuit to simulate setup of Fig.1.2N qubits are chosen to represent an input state corresponding to N photons.Uinput (black box) implements the unitary evolution to prepare the desired input state.Orange, red and blue boxes implement unitary transformations of phase shift, photon loss and the final 50-50 beam splitter respectively.Finally qubits are measured in computational bases.

FIG. 8 .
FIG. 8. Implementation of phase shift, fictitious beam splitter and the final beam splitter on the first photon.
state and in a similar manner, by implementing following unitary evolution on the subspace spanned by |ψ 2 and |ψ 3 we can fix the coefficient of projection onto |ψ 2 : Appendix C: Implementation of phase shift, Fictitious Beam Splitter and the final Beam Splitter by Basis Gates we know the conditional likelihood functions P n0,n1 (φ) ≡ FIG.3.Plot of (∆φ)avg as a function of t for different input states for Nr = 400.Circle points are the values that are calculated by simulating the quantum circuits and the dashed lines are appropriate functions fitted to the data.
|6 :: 0 is the best input state and gives the least error in estimating the phase shift.However, when 0.32 < t < 0.81, |HB(6) has the highest precision and when t < 0.32, |4 :: 2 becomes the best state to estimate the phase shift.When t = 0, the photon loss rate is one, and this means that no photons get detected, and we gain no information about the phase shift, regardless of the input state.So, in this case, the information about the phase shift is the initial probability distribution for the phase shift, and at this point, (∆φ) avg of all input states approach to a fixed value.FIG.4.Plot of (∆φ)avg as a function of t1 for different input states for Nr = 400.Circle points are the values that are calculated by simulating the quantum circuits and the dashed lines are appropriate functions fitted to the data.
is less than 2. This means that at high loss rates, we expect |4 :: 2 and |HB(6) states to estimate phase better than |5 :: 1 .Figure.3 shows the average standard deviations of the input states in estimating the phase shift, for N r = 400 in the symmetric photon loss case.Obviously, for t > 0.81, , and for HB states in