## Abstract

It has been demonstrated that using two-mode squeezed vacuum state for phase estimation can break the Heisenberg limit. Our results reveal that the two-mode squeezed vacuum state is also applied to the optical rotation angle measurement. In our scheme, the resolution and sensitivity of the optical rotation angle signal are the same as the case of phase estimation. For the parameter estimation, phase or rotation angle, we discuss the influences of several imperfect factors on the resolution and sensitivity. First, the effect that the upper limit of photon-number resolving has on the maximum amount of available quantum Fisher information has been analyzed. Then, we have also studied the impacts of both the transmission efficiency in the transmission process and the detection efficiency on the detection results. Finally, conditions where all of the above imperfect elements are taken into account at the same time have also been explored. Additionally, other imperfect factors such as squeezing efficiency and dark counts are briefly discussed.

© 2017 Optical Society of America

## 1. Introduction

Quantum metrology [1–3] is the process of employing a quantum system or making use of quantum mechanical properties to estimate several parameters. As an emerging branch of quantum technology, quantum metrology is drawing plenty of attention in the community and showing its potential. The parameter estimation in interferometer has always been a significant problem in quantum metrology. In modern quantum technology, phase coding and polarization coding have been playing extraordinary roles in lots of fields such as quantum communication [4]. Therefore, phase and optical angle are two types of parameters with the most realistic value.

In 1981, Caves pointed out that the classical interferometry could not break through the shot noise limit due to the fact that one of the input ports is a vacuum state, and the vacuum fluctuation limits the measurement accuracy [5]. In order to improve the accuracy of estimation, several non-classical states that can achieve super-resolution and super-sensitivity are theoretically put forward and experimentally prepared [6–9], such as the N00N state and the squeezed state. And it has been studied that the two-mode squeezed vacuum (TMSV) state for phase estimation can break the Heisenberg limit [10].

In this paper, we demonstrate that the appropriate basis of polarized photons can be utilized to achieve the high-sensitivity estimation of the optical rotation angle under the premise of using the TMSV state. The ultimate estimation result of rotation angle is consistent with that of the phase estimation. This suggests that our scheme is applicable for both phase and optical rotating angle estimation to achieve the same accuracy. On the basis of this scene, we further explore the influence of the finite photon-number resolving on the maximum available quantum Fisher information of the system. We additionally derive an optimal relationship between the average photon number in the incident state and the finite photon-number resolving of the detector. The effects of both transmission efficiency and detection efficiency on estimation accuracy of system are also respectively analyzed. Finally, we simulate and discuss the estimation sensitivity of the system premeditating the above imperfect elements.

## 2. The principle of system

TMSV state is a dual-mode entangled state with forceful correlation between the two modes. Therefore, it has a major application in the researches of both basic issues in quantum mechanics and the development of the quantum information science and technology. The primary device used in our scheme is a Mach-Zehnder interferometer, which differs from the general phase measurement scheme in that we insert two quarter wave plates at two paths after the first beam splitter (see Fig. 1). The linearly polarized photons passing through the quarter wave plate will be converted to circularly polarized photons, and the right-handed circular polarization may be taken as an example in the following derivation. The two paths in interferometer are respectively placed with phase shifter and optical rotation medium, such as Faraday optical crystal [11] or Dove prism [12]. The input state can be written as a infinite superposition in two-mode Fock state basis $|{\psi}_{in}\u3009={\displaystyle \sum _{n=0}^{\infty}\sqrt{(1-t){t}^{n}}{|n,n\u3009}_{H/V}}$, where |*n*, *n*〉 ≡ |*n*〉* _{a}*|

*n*〉

*,*

_{b}*t*=

*N*/(

*N*+ 2),

*N*= 2sinh

^{2}

*s*the average photon number and

*s*is squeezing factor for the input TMSV state.

*H*/

*V*indicates that the mode of the photons is linearly polarized at this time. After the first beam splitter, the TMSV state is naturally decoupled into two single-mode squeezed vacuum states. Note that the mode of photon is circularly polarized at this time. Then two modes go through phase shift and optical rotation, respectively. The phase shift and the optical rotation can be expressed by $\stackrel{\u02c6}{U}(\phi )=\mathrm{exp}(i{\stackrel{\u02c6}{n}}_{A}\phi )$ and $\stackrel{\u02c6}{U}(\theta )=\mathrm{exp}(i{\stackrel{\u02c6}{n}}_{B}\theta )$ under the circularly polarized basis. It is equivalent for circularly polarized photons to withstand optical rotation and phase shift owing to the fact that circularly polarized basis is the eigenstate of the rotational operator. The state after the phase shift and the optical rotation becomes

*φ*and

*θ*represent the phase and the optical rotation angle individually.

For the sake of ensuring the high-sensitivity, it is useful to introduce parity detection as detection strategy. Parity detection was originally proposed by Bollinger *et al.* [13] and later adopted for interferometry metrology by Gerry *et al.* [14, 15]. The parity operator at output port A is $\prod ^{\u02c6}={(-1)}^{\stackrel{\u02c6}{a}\u2020\stackrel{\u02c6}{a}}$. So the expectation value of parity detection is

The above derivation used the property of two-mode Fock state $\u3008n,n|\prod ^{\u02c6}|n,n\u3009\phantom{\rule{1em}{0ex}}={(-1)}^{n}{P}_{n}[\mathrm{cos}2(\phi +\theta )]$ [16] and the summation formula of Legendre polynomial $\sum _{n=0}^{\infty}{P}_{n}(x){t}^{n}}=1/\sqrt{1-2xt+{t}^{2}$ [17]. The details of the beam splitter and other optical components can be found in the work of Yurke *et al*. [18].

To clearly observe the behaviors of the signal resolution, in Fig. 2 we plot the curved surface of resolution as functions of *φ* and *θ*. Here we fix the average photon number *N* = 5 in the input state. As shown in Fig. 3, the full width at half maximum (FWHM) of signal is narrower than the curve of the Malus law, that is, our scheme is super-resolution for both phase and rotation angle. According to the definition of visibility [19], $V=\left({\u3008\prod ^{\u02c6}\u3009}_{\mathrm{max}}-{\u3008\prod ^{\u02c6}\u3009}_{\mathrm{min}}\right)/\left({\u3008\prod ^{\u02c6}\u3009}_{\mathrm{max}}+{\u3008\prod ^{\u02c6}\u3009}_{\mathrm{min}}\right)$, we calculate the visibility *N*/(*N* + 2) by ${\u3008\prod ^{\u02c6}\u3009}_{\mathrm{min}}=1/(N+1)$, ${\u3008\prod ^{\u02c6}\u3009}_{\mathrm{max}}=1$ and Eq. (2). Therefore, the signal visibility can be improved by raising average photon number.

For phase estimation of the TMSV state, it has been verified that parity detection is the optimal strategy by using quantum Fisher information [10]. As for our scheme, the standard deviation of the output signal can be estimated by calculating the error propagation.

The optimal value of above equation $1/\sqrt{N(N+2)}$ shows that the sensitivity of the parameter estimation is preferable to the Heisenberg limit. When the average photon number is large, the above equation is similar to the Heisenberg limit. Additionally, it can be seen that the two parameters have identical consequent in the form of evaluation. Therefore, all conclusions of phase estimation can also be transplanted to our scheme. For the sake of convenience, we uniformly use the parameter *δ* to indicate phase *φ* and rotation angle *θ* in the latter part of discussion.

## 3. Imperfect elements

The imperfect factor that we first discuss is derived from the finite photon-number resolving of photon counting detector. In the practical measurement process, the photon counting detector must have an upper limit with photon-number resolving, which will inevitably miss part of the measurement information. A copious amount of information in the case of loss will result in the poor measurement or even distortion. For balancing the amount of available information, we need to analyze the quantum Fisher information of the system. Since the input state is an infinite superposition of two-mode Fock state, the amount of the entire quantum Fisher information is equal to the sum of the different two-mode Fock states [20]. The transformation of the Mach-Zehnder interferometer for the estimated parameter can be simplified as

Where ${\stackrel{\u02c6}{J}}_{y}$ is the angular momentum operator of the Schwinger representation, therefore, the quantum Fisher information for pure state can be obtained by the formula $F=4\left(\u3008{\stackrel{\u02c6}{J}}_{y}^{2}\u3009-{\u3008{\stackrel{\u02c6}{J}}_{y}\u3009}^{2}\right)$ [21–23]. Quantum Fisher information of two-mode Fock state is *F _{n}* = 2

*n*(

*n*+ 1) through the calculation. Further, we can gain the quantum Fisher information of TMSV state

*F*=

_{TMSV}*N*(

*N*+ 2). This process depends on the orthogonality of the Fock state vector [24]. In the detection, the upper limit of the photon number will be truncated for this infinity sum, and we simulate the normalized quantum Fisher information of the system with distinct average photon number and upper limit of photon-number resolving. We can find that it can obtain more than 99% of the ideal Fisher information when the maximum photon-number resolving of detector is 6 times greater than average photon number in the input state, the detection signal can restore the real detection information at this time. This conclusion also has guiding and referential significance for numerical simulation.

Transmission efficiency is also a major imperfect factor for the detection results, the transmission efficiency in actual detection is affected by many elements, like the transmittance of optical devices, the absorption of photons in the environment and so on. The expectation value of parity signal can be offered by calculating the Wigner function of the output state. The Wigner function for the two-mode squeezed state is given by

Where the parameter *κ* is complex phase factor of squeezing parameter *ξ*, and *ξ* is generally a complex number to define the two-mode squeezed operator $\stackrel{\u02c6}{S}(\xi )=exp\left({\xi}^{*}\stackrel{\u02c6}{a}\stackrel{\u02c6}{b}-\xi {\stackrel{\u02c6}{a}}^{\u2020}{\stackrel{\u02c6}{b}}^{\u2020}\right)$. The modulus of *ξ* is squeezing factor *s* (*ξ* = *se ^{iκ}*). Usually

*κ*= 0 is chosen for simplicity.

*α*and

*β*are parameters characterizing two modes, respectively.

It is approximately equal for the optical path of two modes in the practical application of the interferometry metrology. So we presume that the two paths share the same transmission coefficient *T*, by making the substitutes

*W*(

_{in}*α*,

*β*) we obtain the Wigner function of the output state as ${W}_{out}(\alpha ,\beta )={W}_{in}(\tilde{\alpha},\tilde{\beta})$. Moreover, the expectation value of the parity operator can be supplied

It can be found that the visibility of the signal expectation is still *N/*(*N* + 2), that is, the transmission efficiency does not affect the visibility of the system. In accordance with the error propagation, the sensitivity can be written as

It is liable to check that Eq. (3) is equivalent to Eq. (8) when *T* = 1, and we conduct a simulation analysis with fixing *N* = 5 and *T* = 0.8, 0.9, 1, the results are demonstrated in Fig. 4(a). In the case of transmission efficiency less than 1, large transmission efficiency can still break through the shot noise limit and achieve super-sensitivity. We also investigate the relationship between the average photon number and the minimum transmission efficiency under the premise of achieving super-sensitivity. As can be seen from Fig. 4(b), smaller minimum transmission efficiency is required to achieve super-sensitivity when the average photon number is larger. In other words, it is more robust to resist losses for large average photon number. In particular, when *N* = 680, 10% of transmission efficiency is required to achieve the shot noise limit, from another point of view, at this time the system can tolerate 90% of the loss. This phenomenon is not bothersome to understand from the physical level. Because of the amount of information required for the shot noise limit is proportional to the *N*, and the Heisenberg limit information is proportional to *N*^{2}. Transmission efficiency *T* is a linear coefficient, the greater the *N*, the smaller the *T* that satisfies *N*^{2}*T* greater than *N*, super-sensitivity can be achieved under this circumstance. Of course, here the relationship (*N*^{2}*T* greater than the *N*) is a qualitative analysis based on the magnitude, not a quantitative calculation.

Detection efficiency of the detector is also an extremely critical defect ingredient. For the parity detection, the implementation of the detector is more demanding. Superconducting detectors can only achieve the detection of a few photons, which are obviously not conducive for practical applications. Therefore, most of the experiments to achieve the parity detection are based on the APD array [25–28], APD is a typical binary detector, which only response to the the presence of photons, even if more than one photon reach the APD can only count once a time. So the core idea of the APD array detector is to allocate an average photon number less than one in each APD by fiber distribution as possible.

This scheme takes advantage of readily available assemblies such as APD and time-multiplexed device [29]. In this kind of proposal, the importing photons are splitted into numerous spatially or temporally separate bins to make the presence of more than one photon per bin impossible. Subsequently, each bin is detected with single APD and the resolution results of photon number are obtained by inversely reckoning the number of outcomes from all of the bins.

This type of photon-number-resolving detector is characterized by positive operate valued measure (POVM) formalism [30] in which all photon-counting operations equivalent to POVM elements $\stackrel{\u02c6}{\mu}$.

Where *λ _{n}*

_{,}

*= [*

_{m}*C*·

_{W}*L*(

*η*)]

_{n}_{,}

*,*

_{m}*n*,

*m*and

*W*denote the detection pattern, photon number and the number of bins involved in time-multiplexed process, respectively. The other details about working principle of the detector can refer to the work of Achilles

*et al*. [28], where

*C*and

_{W}*L*(

*η*) represent weight and detection efficiency of detector individually. According to the above analysis, we simulate the influence of the detection efficiency on the detection results in the case of transmission efficiency is 1. Figures 5(a) and 5(b) are the results of resolution and sensitivity, respectively.

It can be viewed in Fig. 5(a) that the visibility of the detection signal increases with the increase of the detection efficiency. But when the detection efficiency is 1, the maximum detection expectation value cannot reach 1 due to the loss of optical fiber in the process of photon assignation. And this loss is inherent in the detection model. This also affects the sensitivity curve to an enormous extent. Moreover, through the calculation we can note that the effect of detection efficiency on the signal FWHM is extremely asthenic. This reveals that the super-resolution characteristic is favorable maintained.

As for the sensitivity, the results in Fig. 5(b) display that the optimal sensitivity will deteriorate as the detection efficiency decreases. The above result is built on the average photon number *N* = 5. It will increase the visibility and sensitivity if increasing the average photon number, however, the change trend with the detection efficiency will not change.

The squeezing efficiency of preparation stage is also an important imperfect factor. The cause of imperfect squeezing efficiency has diverse possibilities. For simplicity, an analytical scenario is to regard the input state as a mixed state containing TMVS state and minor coherent state (minor pump light that is not converted and filtered for some reasons). A well-known conclusion is that the estimation sensitivity of the coherent state is limited by the shot noise limit. Therefore, the mixture of coherent state pulls down the overall estimation sensitivity. This also means that the sensitivity decreases as the squeezing efficiency decreases. As for the resolution, the peak position of coherent state and that of TMSV state are respectively 0 and ±*π/*2 in the case of parity detection. Hence, the impact for visibility is not significant except for a bulge appears between the two signal peaks.

The numerical simulations that take all imperfect elements into account are discussed. We still fix the *N* =5 and the resolution limit of detector is 6*N*, which is to ensure the adequate amount of information. The simulation results are individually illustrated in Figs. 6(a) and 6(b). Figure 6(a) suggests that the effect of the two kinds of efficiencies on the visibility of the signal resolution. It is not exertive to find that the resolution is not affected by the transmission efficiency, which confirms our previous conclusion. In addition, the decrease in detection efficiency brings about resolution visibility decline. As to the sensitivity, Fig. 6(b) shows that the optimum sensitivity deteriorates with either the transmission efficiency or the detection efficiency. In contrast, this diversification is more conspicuous as the detection efficiency changes. The above results reveal that the entire detection upshot is more affected by the detection efficiency of the detector.

In addition, response time delay and dark counts of detector are momentous imperfect factors. Suppose that the rate of dark counts is *r*, the probability of *m* dark counts follows the Poisson distribution $D(m)={e}^{-r}{\displaystyle \sum _{m=0}^{\infty}{r}^{m}/m!}$. The parity is estimated just from the parity of the photon number at output port. Ideally, the sum of the photon number detected by two detectors is always even due to the property of TMSV state and the probability of even counts is ${P}_{e}=\left(1+\u3008\prod ^{\u02c6}\u3009\right)/2$. Therefore, the count is invalidated if the total number of detected photon is odd. The probability of even counts under such a protocol can be rewritten as

Where *D _{A}*

_{/}

*(*

_{B}*e*/

*o*) indicates that the probability of even or odd dark counts generated by detector

*A*or detector

*B*. The situation, odd dark counts of one detector and even dark counts of the other detector, is invalid owing to violating the protocol. Similarly, we obtain ${P}_{o}^{\prime}=(1+{e}^{-4r})/4-\u3008\prod ^{\u02c6}\u3009{e}^{-2r}/2$. The expectation value of the output signal existing dark counts is

An intuitive phenomenon is that output signal remains the visibility of *N/N* + 2. The system sensitivity can be calculated and we plot Fig. 7 for observing the effect of dark counts on the system sensitivity. Under the current technology, the range of *r* is generally between 10^{−8} to 10^{−3}[31]. The response time delay forces the sampling detection gates to increase the width and introduces added dark counts. The role of the response time delay can be seen as a negative impact, which is to raise the dark counts. Take *N*=5 as an example with *r* = 10^{−2} (simulating the case of taking into account the response time delay and dark counts) and *r* = 10^{−3} (simulating the situation of only considering dark counts). It can be clearly observed that the effects of both response time delay and dark counts on sensitivity are slight and the sensitivity is still lower than the Heisenberg limit.

In general, response time delay and dark counts do not interfere the visibility of output signal, while the property of super-resolution is not changed. A relaxing phenomenon is that the influences of two imperfect factors on sensitivity are also feeble under the circumstance of actual detector parameter.

## 4. Conclusion

In conclusion, we manifest a quantum metrology scheme of rotating angle measurement based on the TMSV state. It can be used to estimate the rotation angle and obtain the same resolution and sensitivity compared to the phase estimation. On the basis of this theory, we also discuss the relationship between the upper limit *N _{up}* of photon-number resolving and the average photon number

*N*in the incident state from the perspective of quantum information theory. The results indicate that the detector can obtain 99% of information about the estimate parameter when the condition

*N*≥ 6

_{up}*N*is satisfied. Subsequently, the effects of transmission efficiency on sensitivity and resolution are also discussed. The outcomes demonstrate that the transmission efficiency will not affect the visibility of the system, but the sensitivity will decrease as the transmission efficiency decreases. In addition, increasing the average photon number will increase the ability of the system to achieve super-sensitivity in the presence of transmission loss. The effect of the detection efficiency is also considered, and decrease in detection efficiency will result in reduction in resolution and sensitivity. Additionally, all the imperfections are analyzed synthetically. Numerical results establish that only the detection efficiency has an influence on the system resolution, while the transmission efficiency and detection efficiency will contribute to the deterioration of the system sensitivity. The influence of the detection efficiency on the sensitivity is more serious. Finally, the impacts of multitudinous imperfect elements on the detection results have been studied, and the effect of response time delay and that of dark counts are inconspicuous in accordance with theoretical analysis. Overall, the function of TMSV state has also been expanded, which is of realistic connotation to the application of TMSV state.

## Funding

The Scientific Research Fund of Heilongjiang Provincial Education Department (12541839).

## References and links

**1. **V. Giovannetti, S. Lloyd, and L. Maccone, “Quantum metrology,” Phys. Rev. Lett. **96**(1), 010401 (2006). [CrossRef] [PubMed]

**2. **V. Giovannetti, S. Lloyd, and L. Maccone, “Advances in quantum metrology,” Nat. Photonics **5**(4), 222–229 (2011). [CrossRef]

**3. **C. F. Roos, M. Chwalla, K. Kim, M. Riebe, and R. Blatt, “‘Designer atoms’ for quantum metrology,” Nature **443**(7109), 316–319 (2006). [CrossRef] [PubMed]

**4. **J. W. Pan, Z. B. Chen, C. Y. Lu, H. Weinfurter, A. Zeilinger, and M. Zukowski, “Multiphoton entanglement and interferometry,” Rev. Mod. Phys. **84**(2), 777 (2012). [CrossRef]

**5. **C. M. Caves, “Quantum-mechanical noise in an interferometer,” Phys. Rev. D **23**(8), 1693 (1981). [CrossRef]

**6. **J. Joo, W. J. Munro, and T. P. Spiller, “Quantum metrology with entangled coherent states,” Phys. Rev. Lett. **107**(8) 083601 (2011). [CrossRef] [PubMed]

**7. **M. D. Lang and C. M. Caves, “Optimal quantum-enhanced interferometry using a laser power source,” Phys. Rev. Lett. **111**(17), 173601 (2013). [CrossRef] [PubMed]

**8. **C. E. Vollmer, C. Baune, A. Samblowski, T. Eberle, V. Handchen, J. Fiurasek, and R. Schnabel, “Quantum up-conversion of squeezed vacuum states from 1550 to 532 nm,” Phys. Rev. Lett. **112**(7), 073602 (2014). [CrossRef] [PubMed]

**9. **M. Lebugle, M. Grafe, R. Heilmann, A. Perez-Leija, S. Nolte, and A. Szameit, “Experimental observation of N00N state Bloch oscillations,” Nat. Commun. **6**(1), 8273 (2015). [CrossRef] [PubMed]

**10. **P. M. Anisimov, G. M. Raterman, A. Chiruvelli, W. N. Plick, S. D. Huver, H. Lee, and J. P. Dowling, “Quantum metrology with two-mode squeezed vacuum: parity detection beats the Heisenberg limit,” Phys. Rev. Lett. **104**(10), 103602 (2010). [CrossRef] [PubMed]

**11. **W. Zhang, Q. Qi, J. Zhou, and L. Chen, “Mimicking Faraday rotation to sort the orbital angular momentum of light,” Phys. Rev. Lett. **112**(15), 153601 (2014). [CrossRef] [PubMed]

**12. **J. Leach, M. J. Padgett, S. M. Barnett, S. Franke-Arnold, and J. Courtial, “Measuring the orbital angular momentum of a single photon,” Phys. Rev. Lett. **88**(25), 257901 (2002). [CrossRef] [PubMed]

**13. **J. J. Bollinger, W. M. Itano, D. J. Wineland, and D. J. Heinzen, “Optimal frequency measurements with maximally correlated states,” Phys. Rev. A **54**(6), R4649 (1996). [CrossRef] [PubMed]

**14. **C. C. Gerry, “Heisenberg-limit interferometry with four-wave mixers operating in a nonlinear regime,” Phys. Rev. A **61**(4), 043811 (2000). [CrossRef]

**15. **C. C. Gerry and R. A. Campos, “Generation of maximally entangled photonic states with a quantum-optical Fredkin gate,” Phys. Rev. A **64**(6), 063814 (2001). [CrossRef]

**16. **R. A. Campos, C. C. Gerry, and A. Benmoussa, “Optical interferometry at the Heisenberg limit with twin Fock states and parity measurements,” Phys. Rev. A **68**(2), 023810 (2003). [CrossRef]

**17. **M. Abramowitz and I. Stegun, *Handbook of Mathematical functions* (Dover, 1965).

**18. **B. Yurke, S. L. McCall, and J. R. Klauder, “SU (2) and SU (1, 1) interferometers,” Phys. Rev. A **33**(6) 4033 (1986). [CrossRef]

**19. **J. P. Dowling, “Quantum optical metrology-the lowdown on high-N00N states,” Contemp. Phys. **49**(2), 125–143 (2008). [CrossRef]

**20. **P. Liu, P. Wang, W. Yang, G. R. Jin, and C. P. Sun, “Fisher information of a squeezed-state interferometer with a finite photon-number resolution,” Phys. Rev. A **95**(2), 023824 (2017). [CrossRef]

**21. **S. Knysh, V. N. Smelyanskiy, and G. A. Durkin, “Scaling laws for precision in quantum interferometry and the bifurcation landscape of the optimal state,” Phys. Rev. A **83**(2), 021804 (2011). [CrossRef]

**22. **J. Liu, X. Jing, and X. Wang, “Phase-matching condition for enhancement of phase sensitivity in quantum metrology,” Phys. Rev. A **88**(4), 042316 (2013). [CrossRef]

**23. **Q. S. Tan, J. Q. Liao, X. Wang, and F. Nori, “Enhanced interferometry using squeezed thermal states and even or odd states,” Phys. Rev. A **89**(5), 053822 (2014). [CrossRef]

**24. **M. Takeoka, K. P. Seshadreesan, C. L. You, S. Izumi, and J. P. Dowling, “Fundamental precision limit of a Mach-Zehnder interferometric sensor when one of the inputs is the vacuum,” arXiv:1705.09506, (2017).

**25. **I. Afek, O. Ambar, and Y. Silberberg, “High-NOON states by mixing quantum and classical light,” Science , **328**(5980), 879–881 (2010). [CrossRef] [PubMed]

**26. **Y. Israel, S. Rosen, and Y. Silberberg, “Supersensitive polarization microscopy using NOON states of light,” Phys. Rev. Lett. **112**(10), 103604 (2014). [CrossRef] [PubMed]

**27. **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,” Nat. Photonics **5**(1), 43–47 (2011). [CrossRef]

**28. **D. Achilles, C. Silberhorn, C. Sliwa, K. Banaszek, I. A. Walmsley, M. J. Fitch, B. C. Jacobs, T. B. Pittman, and J. D. Franson, “Photon-number-resolving detection using time-multiplexing,” J. Mod. Opt. **51**(9–10), 1499–1515 (2004). [CrossRef]

**29. **D. Achilles, C. Silberhorn, C. Sliwa, K. Banaszek, and I. A. Walmsley, “Fiber-assisted detection with photon number resolution,” Opt. Lett. **28**(23), 2387–2389 (2003). [CrossRef] [PubMed]

**30. **J. S. Lundeen, A. Feito, H. Coldenstrodt-Ronge, K. L. Pregnel, C. Silberhorn, T. C. Ralph, J. Eisert, M. B. Plenio, and I. A. Walmsley, “Tomography of quantum detectors,” Nat. Phys. **5**(1), 27–30 (2009). [CrossRef]

**31. **R. H. Hadfield, “Single-photon detectors for optical quantum information applications,” Nat. Photonics **3**, 696–705 (2009). [CrossRef]