Abstract
Quantum metrology can achieve enhanced sensitivity for estimating unknown parameters beyond the standard quantum limit. Recently, multiple-phase estimation exploiting quantum resources has attracted intensive interest for its applications in quantum imaging and sensor networks. For multiple-phase estimation, the amount of enhanced sensitivity is dependent on quantum probe states, and multi-mode N00N states are known to be a key resource for this. However, its experimental demonstration has been missing so far since generating such states is highly challenging. Here, we report generation of multi-mode N00N states and experimental demonstration of quantum enhanced multiple-phase estimation using the multi-mode N00N states. In particular, we show that the quantum Cramer-Rao bound can be saturated using our two-photon four-mode N00N state and measurement scheme using a 4 × 4 multi-mode beam splitter. Our multiple-phase estimation strategy provides a faithful platform to investigate multiple parameter estimation scenarios.
Similar content being viewed by others
Introduction
Quantum metrology has attracted intensive interest in recent years, as it allows us to estimate an unknown parameter with enhanced sensitivity over classical approaches. Developments are now directed towards various applications such as microscopy1,2,3,4,5, biological imaging6,7,8,9,10 and sensor network11,12,13. In such practical applications, quantum metrology is naturally extended to multiple parameter scenario, since the number of parameters that affect a physical process is generally more than one. For estimating multiple parameters, simultaneous estimation is desirable as it can achieve higher precision than an individual estimation of each parameter using the same amount of resources14,15. Furthermore, if a physical system of interest has time dynamics, estimation of multiple parameters has to be done simultaneously. However, unlike single parameter estimation, it is non-trivial to optimize multiple parameter estimation to achieve the maximum sensitivity. Various strategies have been thus proposed to enhance the sensitivity by optimizing either the probe states or the measurement scheme14,15,16,17,18,19,20,21,22. For example, quantum strategies exploiting various quantum probe states such as Greenberger–Horne–Zeilinger states23, single-photon Fock states24, squeezed states12,25,26,27,28, Holland-Bernett states29,30,31, and N00N states14,32,33 have been extensively studied.
In particular, N00N states have been outstanding to investigate the fundamental quantum limit of quantum metrology given by the Heisenberg uncertainty principle with a fixed number of particles34. In a single-phase estimation scheme, N00N states can saturate the Heisenberg limit thanks to its largest number variance between the two modes34. Enhanced sensitivity beyond the standard quantum limit has been experimentally demonstrated with N00N states35,36. Recently, the concept of N00N states has been extended to generalized multi-mode N00N states to investigate the quantum enhancement in multiple-phase estimation14. Moreover, for multiple-phase estimation with limited resources, it has been known that multi-mode N00N states allow achieving the enhanced sensitivity outperforming the other quantum probe states or classical strategies for a multi-mode interferometer14,15,32. However, since generation of such multi-mode quantum probe states is challenging, experimental demonstrations of multiple parameter quantum metrology have been limited to utilizing quantum states other than multi-mode N00N states12,15,23,24,27,29.
In this work, we experimentally demonstrate quantum enhanced multiple-phase estimation using a multi-mode N00N state with photon number N = 2 and mode number m = 4. To this end, we propose a scheme for generating multi-mode N00N states. At first, we observe that the measured interference fringes exhibit phase super-resolution as a function of phase differences φi (i = 1, 2, 3) of the generated 4-mode 2002 state with a measurement scheme using a 4 × 4 multi-mode beam splitter, so-called a quarter37,38. Then, we demonstrate quantum enhanced phase sensitivity of the 4-mode 2002 state by analyzing the Cramer-Rao bound (CRB) and the quantum Cramer-Rao bound (QCRB)14,16,17,39. Moreover, we show that the measured sensitivity is better than a coherent state \(|\alpha \rangle\), which is a classical probe state, as well as another quantum probe state prepared by single-photon Fock states24. Our results can motivate an investigation into the quantum strategies using multi-mode and multi-particle entanglement to develop the quantum enhanced multiple parameter metrology.
Results
Multiple-phase estimation scenario
Let us begin by introducing the system model for multiple-phase estimation schemes with multi-mode N00N states defined as
where N is the number of photons, and they are distributed along m modes14,32,33. The multi-mode N00N state is a coherent superposition of all possibilities where N photons in one mode and none in any of the other m − 1 modes34. Then, we theoretically analyze the sensing scheme of d = m − 1 unknown multiple phases and the fixed photon number N = 2 with different probe states. A quantum probe state undergoes individual phase shifts φ = {φ1, φ2, ⋯ , φd} and the phase shifted state is combined using an m × m multi-mode beam splitter and then detected by photon number-resolving detectors (PNRDs) at each mode. Here, the goal is to minimize the total uncertainty of the phase estimation governed by the CRB and the QCRB. The lower bound of the sum of the variance of each phase estimation given by the CRB and the QCRB is14,16,17,39
where FC(φ) is the classical Fisher information matrix (CFIM), FQ(φ) is the quantum Fisher information matrix (QFIM), and μ is the number of measurements. Note that F−1 refers to the inverse of the Fisher information matrix F. The first and second inequalities are direct consequences of the CRB and the QCRB, respectively. Note that in our experiments, the CRB can always be saturated, asymptotically in μ, by using maximum likelihood estimator, and the mean of phases are expected to converge to the true values17,39.
We theoretically analyze the CRB and the QCRB of d multiple-phase estimation scheme with two different probe states. The first probe state we consider is a classical state prepared by injecting a coherent state \(\left|\alpha \right\rangle\) with an average photon number \(\overline{N}=2\) into one of the input ports of an m × m multi-mode beam splitter (Fig. 1a), and the other probe state is the two-photon m-mode N00N state (Fig. 1b). The measurement scheme is identical for two probe states. Each mode of the probe states is combined at another m × m multi-mode beam splitter after undergoing phase shifts, and then measured by PNRDs. We provide the CRB and the QCRB values of the total variance \(\sum {\left|{{\Delta }}\varphi \right|}^{2}\) for two probe states of coherent states and multi-mode N00N states with various number of phases d in Fig. 1c. Note that multi-mode N00N states always have lower phase uncertainty than coherent states, i.e., classical probe states.
Generation of the 4-mode 2002 state
The conceptual diagram of our proposed scheme for preparing the multi-mode N00N state with N = 2 and m = 4 is shown in Fig. 2a. The generation process of 4-mode 2002 state \(|{\psi }_{4}^{2}\rangle\) is the following:
Here, \({|}{{{\Phi }}}^{+}\rangle =\frac{1}{\sqrt{2}}(|{1}_{{\alpha }_{0}}^{H}{1}_{{\alpha }_{1}}^{H}\rangle +|{1}_{{\alpha }_{0}}^{V}{1}_{{\alpha }_{1}}^{V}\rangle )\) is the triplet Bell state where, for example, \(|{1}_{{\alpha }_{0}}^{H}\rangle\) denotes the horizontally polarized single-photon state in the mode α0. Our scheme can generate \(\left|{\psi }_{4}^{2}\right\rangle\) with a unity conversion probability from the pre-selected initial \(\left|{{{\Phi }}}^{+}\right\rangle\) state into \(\left|{\psi }_{4}^{2}\right\rangle\), thus two-photons always move together, see Methods for the detailed information. Then, prepared \(\left|{\psi }_{4}^{2}\right\rangle\) undergoes phase encoding φi (i = 1, 2, 3) and it is combined by a quarter as shown in Fig. 2b. The output states are measured by single-photon detectors located at each mode. Note that one output port of the quarter is split into two (\({b}_{0^{\prime} }\) and b0″) using a 50/50 fiber beam splitter to implement a pseudo-PNRD. Post-selected two-photon coincidence counts of both output port \({b}_{0^{\prime} }\) and b0″ correspond to the two-photon counting at b0 with success probability of 1/2. The schematic of experimental setup is shown in Fig. 2c.
As a first step to verify the generation of \(\left|{\psi }_{4}^{2}\right\rangle\), we investigate the coherence among all four modes of \(\left|{\psi }_{4}^{2}\right\rangle\) by observing the interference fringes between the reference mode (a0) and one of the other modes (ai). We measure the two-photon coincidence counts on \({b}_{0^{\prime} }\) and b0″, \({{{\mbox{C}}}}_{b0^{\prime} b0^{\prime\prime} }\), of two input modes while we block the other two input modes. The phase shift φi is realized by adjusting the optic axis angle θi of the half waveplate (HWP) located between two quarter waveplates (QWPs). Note that φi = 2θi. As shown in Fig. 3a–c, results of \({{{\mbox{C}}}}_{b0^{\prime} b0^{\prime\prime} }\) clearly reveal two times faster sinusoidal modulations compared to the single-photon input case due to the λ/2 photonic de Broglie wavelength of two-photon N00N states40,41,42.
Then, in order to directly confirm the generation of \(\left|{\psi }_{4}^{2}\right\rangle\), we measure two-photon probability distributions for the output states when we use all of the input modes a0, a1, a2, and a3 of \(\left|{\psi }_{4}^{2}\right\rangle\) while we vary one of the phase encoding φi. The experimental results of the interference fringes from all of the modes are shown in Fig. 3d–f, and they are compared with the theoretical calculations. We obtain the theoretical predictions of \({{{\mbox{C}}}}_{{b}_{0}{b}_{0}}\), \({{{\mbox{C}}}}_{{b}_{0}{b}_{1}}\), \({{{\mbox{C}}}}_{{b}_{0}{b}_{2}}\), and \({{{\mbox{C}}}}_{{b}_{0}{b}_{3}}\) based on the ideal \(\left|{\psi }_{4}^{2}\right\rangle\) and the experimentally reconstructed quarter transition matrix. Then we experimentally measured post-selected coincidence counts of Ci for \(i={b}_{0^{\prime} }{b}_{0^{\prime\prime} },{b}_{0^{\prime} }{b}_{1},{b}_{0^{\prime} }{b}_{2}\) and \({b}_{0^{\prime} }{b}_{3}\). See Supplementary Note 1 for the detailed information on a quarter transition matrix and theoretical calculations. The results of Fig. 3d–f shows that the experimentally obtained interference fringes without normalization are very well-matched to the theoretical calculations, thus we can confirm that the prepared input state is \(\left|{\psi }_{4}^{2}\right\rangle\).
Experimental multiple-phase estimation
Then, we investigate the sensitivity bound of multiple-phase estimation using our prepared \(\left|{\psi }_{4}^{2}\right\rangle\) as a probe state. At first, we theoretically calculate the CRB with an ideal \(\left|{\psi }_{4}^{2}\right\rangle\) and an ideal quarter when the CRB saturates the QCRB. We obtain the theoretical two-photon detection probability set of {Pl(φ)} (l = 0, 1, . . . , 9) with a set of projectors \(\{{\hat{{{\Pi }}}}_{l}\}=\left\{\left|2000\right\rangle \right.\left\langle 2000\right|\), \(\left|0200\right\rangle \left\langle 0200\right|\), \(\left|0020\right\rangle \left\langle 0020\right|\), \(\left|0002\right\rangle \left\langle 0002\right|\), \(\left|1100\right\rangle \left\langle 1100\right|\), \(\left|0011\right\rangle \left\langle 0011\right|\), \(\left|1010\right\rangle \left\langle 1010\right|\), \(\left|0101\right\rangle \left\langle 0101\right|\), \(\left|1001\right\rangle \left\langle 1001\right|\), \(\left|0110\right\rangle \left.\left\langle 0110\right|\right\}\) satisfying the normalizing condition ∑lPl(φ) = 1. Note that P0 = P1 = P2 = P3, P4 = P5, P6 = P7, and P8 = P9 for all φ (See Methods). Then CFIM is given by
where j and k can be 1, 2, and 3. The minimum value of the CRB is obtained to be \(\,{{\mbox{Tr}}}\,[{{{{{{{{\bf{F}}}}}}}}}_{C}^{-1}({{{{{{{\boldsymbol{\varphi }}}}}}}})]\) = 1.5 where φ1 ≃ π/2, φ2 = 0, and φ3 = π/2, and it saturates the QCRB = 1.5. In order to experimentally estimate the CRB value, we obtained interference fringes by scanning φ1 near the point where we expected both the CRB and the QCRB to be saturated (φ2 = −0.07π and φ3 = 0.52π) for the prepared \({\left|{\psi }_{4}^{2}\right\rangle }_{\exp }\) probe state (See Methods). Two-photon detection probabilities \({P}_{{b}_{0}{b}_{0}}^{\,{{\mbox{m}}}\,}\), \({P}_{{b}_{0}{b}_{1}}^{\,{{\mbox{m}}}\,}\), \({P}_{{b}_{0}{b}_{2}}^{\,{{\mbox{m}}}\,}\) and \({P}_{{b}_{0}{b}_{3}}^{\,{{\mbox{m}}}\,}\) are then obtained from the measured post-selected coincidence counts with μ ≃ 8, 144, Ci, (\(i={b}_{0^{\prime} }{b}_{0^{\prime\prime} },{b}_{0^{\prime} }{b}_{1},{b}_{0^{\prime} }{b}_{2}\), and \({b}_{0^{\prime} }{b}_{3}\)), which were appropriately normalized assuming the following relations \({P}_{0}^{\exp }={P}_{1}^{\exp }={P}_{2}^{\exp }={P}_{3}^{\exp }={P}_{{b}_{0}{b}_{0}}^{\,{{\mbox{m}}}\,}\), \({P}_{4}^{\exp }\)\(={P}_{5}^{\exp }={P}_{{b}_{0}{b}_{1}}^{\,{{\mbox{m}}}\,}\), \({P}_{6}^{\exp }={P}_{7}^{\exp }={P}_{{b}_{0}{b}_{2}}^{\,{{\mbox{m}}}\,}\), and \({P}_{8}^{\exp }={P}_{9}^{\exp }={P}_{{b}_{0}{b}_{3}}^{\,{{\mbox{m}}}\,}\). Note that the assumed relations are always satisfied for Pl(φ) with an ideal quarter and the experimentally reconstructed quarter transition matrix is close to an ideal quarter. Then obtained \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\) are functions of unknown phases φ and used to calculate the derivatives in Eq. (4). The detailed relation between \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\) and \({P}_{{b}_{i}{b}_{j}}^{\,{{\mbox{m}}}\,}\) are provided in Supplementary Note 3.
Experimentally obtained two-photon detection probabilities are shown in Fig. 4a where our experimental data are well-matched to our fitting function \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\), which are obtained from \({\left|{\psi }_{4}^{2}\right\rangle }_{\exp }\) and \({{{{{{{{\bf{U}}}}}}}}}_{q,\exp }\) (see “Methods”). The CFIM can be obtained by using \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\), and the diagonal terms of the CFIM at various φ1 are plotted in Fig. 4b (see “Methods” for the detailed information). Note that the maximum values of all diagonal terms of the CFIM are 3 at φ1 ≃ 0.5π with an ideal \(\left|{\psi }_{4}^{2}\right\rangle\) state and an ideal quarter. Then we numerically find the minimum CRB of \(\,{{\mbox{Tr}}}\,[{({{{{{{{{\bf{F}}}}}}}}}_{C}^{\exp })}^{-1}({{{{{{{\boldsymbol{\varphi }}}}}}}})]=1.85\pm 0.01\) when φ1 ≃ 0.47π, φ2 = −0.07π, and φ3 = 0.52π from \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\), which is clearly smaller than the CRB of the coherent state \(\,{{\mbox{Tr}}}\,[{({{{{{{{{\bf{F}}}}}}}}}_{C,{{{{{{{\rm{coh}}}}}}}}})}^{-1}({{{{{{{\boldsymbol{\varphi }}}}}}}})]=3\) with the same average photon number \(\overline{N}=2\) as shown in Fig. 4c (see Supplementary Note 3 for the detailed calculation on estimating the CRB and results with respect to φ2 and φ3.). Note that the ideal CRB value can achieve 1.5 to saturate the QCRB with an ideal \(\left|{\psi }_{4}^{2}\right\rangle\) state and an ideal quarter at φ1 ≃ 0.5π, φ2 = 0, and φ3 = 0.5π.
In our experiment, experimental errors are mainly from non-unity visibility of the observed interference, a phase fluctuation in each arm of an interferometer, a normalization assumption due to lack of superconducting nanowire single-photon detector (SNSPD) channels, and the fact that the quarter (\({{{{{{{{\bf{U}}}}}}}}}_{q,\exp }\)) used in our experiment is slightly different from an ideal quarter Uq, see Methods for comparison between \({{{{{{{{\bf{U}}}}}}}}}_{q,\exp }\) and Uq. In particular, in Eq. (4), with an ideal quarter, ∂Pl(φ)/∂φj = 0 should be at the same φj value for all Pl(φ), however, a non-ideal quarter transition matrix \({{{{{{{{\bf{U}}}}}}}}}_{q,\exp }\) makes \(\partial {P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})/\partial {\varphi }_{j}=0\) happens at slightly different φj values. One can notice it from Fig. 4a that \({P}_{{b}_{0}{b}_{0}}^{\,{{\mbox{m}}}\,}\), \({P}_{{b}_{0}{b}_{1}}^{\,{{\mbox{m}}}\,}\), and \({P}_{{b}_{0}{b}_{3}}^{\,{{\mbox{m}}}\,}\) (\({P}_{{b}_{0}{b}_{2}}^{\,{{\mbox{m}}}\,}\)) do not have their minimum (maximum) value at φ1 = 0.5π. This is the reason why the minimum value of the CRB is not obtained at φ1 ≃ 0.5π but φ1 ≃ 0.47π.
Furthermore, we emphasize that the enhanced sensitivity can be obtained by the prepared \(\left|{\psi }_{4}^{2}\right\rangle\) probe state compared to other quantum probe states using single-photon Fock states \(\left|{\psi }_{{{\mbox{Fock}}}}\right\rangle\) proposed in ref. 24, where \(\left|1100\right\rangle\) state is used as an input state instead of \(\left|\alpha \right\rangle\) in Fig. 1a (detailed calculations are provided in Supplementary Note 2). Moreover, we theoretically compare the sensitivity bounds between our multi-mode N00N state and an amplitude-unbalanced multi-mode N00N state \(\left|{\psi }_{u}\right\rangle\), which is proposed in refs. 14,15,39. The amplitude-unbalanced multi-mode N00N state has the form of \(\left|{\psi }_{u}\right\rangle =\alpha \left|N0\cdots 0\right\rangle +\beta (\left|0N0\cdots 0\right\rangle +\cdots +\left|0\cdots 0N\right\rangle )\) with α2 + dβ2 = 1 and \(\beta =1/\sqrt{d+\sqrt{d}}\). \(\left|{\psi }_{u}\right\rangle\) is known to have the minimum QCRB among the multi-mode N00N states14,15,39. In a measurement scheme using a quarter and PNRDs, the QCRB and the CRB values of \(\left|{\psi }_{u}\right\rangle\) with N = 2 and d = 3 are theoretically calculated to be 1.4 and 1.62, respectively. Here, we find that even though \(\left|{\psi }_{u}\right\rangle\) has the lower QCRB of 1.4 than 1.5 of our \(\left|{\psi }_{4}^{2}\right\rangle\), \(\left|{\psi }_{4}^{2}\right\rangle\) can provide a better sensitivity (smaller CRB) 1.5 than 1.62 of \(\left|{\psi }_{u}\right\rangle\) with a realistic measurement scheme using a quarter. Note that an optimal measurement saturating the QCRB may not be experimentally feasible14,17.
In Table 1, we summarize the ideal QCRB and CRB values for various probe states as well as the CRB values obtained from our experimental results. We emphasize that the experimentally obtained CRB value of 1.85 ± 0.01 provides a better sensitivity than the ideal CRB values of 3 for \(\left|\alpha \right\rangle\) and 2.44 for \(\left|{\psi }_{{{\mbox{Fock}}}}\right\rangle\).
Our experiments use the post-selection technique and does not consider the photonic losses due to experimental imperfection43 and lack of high-efficiency PNRDs36. By considering the post-selection probability and photon losses, the CRB of our \(\left|{\psi }_{4}^{2}\right\rangle\) state cannot outperform the classical strategy. However, a genuine quantum enhancement can be achieved if one uses the state-of-art technologies such as high-efficiency PNRDs and optimized low-loss optical components with our \(\left|{\psi }_{4}^{2}\right\rangle\) state. Note that the post-selection technique does not affect the proof-of-concept verification of quantum enhancement and post-selection is standard technique used in almost previous quantum metrology experiment15,23,24,29. The corresponding quantitative analysis as well as the detailed discussions on photonic losses in our experimental setup are provided in Supplementary Note 4.
Discussion
In conclusion, we proposed a scheme for generating a multi-mode N00N state and experimentally demonstrate that the prepared quantum probe state is the 4-mode 2002 state by observing various interference fringes shown in Fig. 3 using a quarter and photon number resolving detection using post-selective pseudo-PNRDs. Moreover, we exploit the prepared 4-mode 2002 state as a quantum probe state for simultaneously estimating three phases of a 4-mode interferometer with quantum enhanced sensitivity. Then, we confirm that the CRB obtained by our 4-mode 2002 state and measurement scheme can saturate the QCRB. Our results provide a practical platform to investigate intriguing issues in the field of quantum multiple parameter metrology. At first, we emphasize that our scheme can be extended to generation of higher-mode N00N states. For instance, one can exploit a multiple-path Sagnac interferometer44 to increase the number of modes from 4 to 4n, i.e., \(\left|{\psi }_{4n}^{2}\right\rangle\), where n is the number of Sagnac interferometers, and then one can estimate up to 4n − 1 phases simultaneously. Note that another scheme for generating multi-mode N00N states with N ≥ 2 has been theoretically proposed33. However, experimental demonstration seems to be challenging within current technology since it requires extremely strong nonlinearity or deterministic generation of multi-photon states. Another interesting future direction would be finding a realistic measurement scheme for minimizing the CRB. In general, the CRB obtained by measurement scheme using a balanced multi-mode beam splitter cannot saturate the QCRB, see Fig. 1c with d = 2, 4, 5. Since the optimal measurement saturating the QCRB involves complex multi-photon states, it may not be experimentally feasible14,17. Hence, finding an experimentally realistic measurement scheme minimizing the CRB is essential for practical applications. Our results have direct applications for the quantum enhanced phase object imaging requiring a low photon flux7,45,46. In addition, our results can pave the way for demonstrating distributed quantum enhanced multiple-phase estimation by increasing the number of phases in local interferometers consisting of distributed quantum sensors12,23,27.
Methods
4-mode 2002 state preparation
We used a CW single frequency laser operating at a center wavelength of 780 nm. The polarization of pump laser is set to \(\left|D\right\rangle =\left(\left|H\right\rangle +\left|V\right\rangle \right)/\sqrt{2}\) polarization. The 10 mm-thick type-II periodically poled KTiOPO4 (PPKTP) crystal with 46.15 μm poling period is located at the center of the Sagnac interferometer, which consists of a dual wavelength polarizing beam splitter (PBS), a dual wavelength HWP whose optic axis is oriented at 45∘, and two dual wavelength mirrors as shown in Fig. 2c47. Here, dual wavelength optical components are designed for working at both 780 and 1560 nm photons. The horizontal (vertical) polarization component of the pump laser is transmitted (reflected) at dual wavelength PBS. The vertically polarized pump laser is changed to the horizontal polarization after transmitting a dual wavelength HWP. Then, the both of clockwise and counter-clockwise propagating pump laser have horizontal polarization in front of the PPKTP crystal and they probabilistically create photon pairs having orthogonal polarizations of \(\left|H\right\rangle\) and \(\left|V\right\rangle\), respectively, via degenerate spontaneous parametric down conversion process. Then, after passing through the dual wavelength PBS, the counter propagating photon pair states are emerging to \((\left|HV\right\rangle +{e}^{i\phi }\left|VH\right\rangle )/\sqrt{2}\), where ϕ is relative phase between two states. Then, we can prepare \(\left|{{{\Phi }}}^{+}\right\rangle\) state using a set of WPs. The measured heralding efficiency of \(\left|{{{\Phi }}}^{+}\right\rangle\) state is 37% with our SNSPD whose detection efficiency is 80%.
Two photons prepared in \(\left|{{{\Phi }}}^{+}\right\rangle\) state simultaneously entered at both input ports of a lateral beam splitter (LBS), then we can observe the Hong-Ou-Mandel (HOM) interference48. The path length difference between two photons are matched to minimize the coincidence count rate detected at both output ports. See Supplementary Note 1 for the HOM interference results with \(\left|{{{\Phi }}}^{+}\right\rangle\) input state. After the LBS, we can prepare the two-photon N00N state with different polarization (See the second line of Eq. (3)). The horizontally (vertically) polarized photons are transmitted (reflected) at the PBS so that they split into 4 path modes depending on the polarization states. Note that only two-photon states can exist in each mode. By rotating the polarization state of the vertically polarized photons into the horizontal polarization, we can prepare \(\left|{\psi }_{4}^{2}\right\rangle\).
Theoretical analysis of \(\left|{\psi }_{4}^{2}\right\rangle\)
Our probe state is the 4-mode 2002 state of the form
where \({\hat{a}}_{{{{{{{{\rm{0}}}}}}}}}^{{{{\dagger}}} }\) is a creation operator, which creates a single-photon in the input mode a0 of a quarter. A quarter has four input modes (a0, a1, a2, a3) and four output modes (b0, b1, b2, b3), respectively, as shown in Fig. 2c. The unitary matrices for the phase encoding Uφ and the ideal quarter transition Uq are given by37,38,
and
respectively. After the initial state \(\left|{\psi }_{4}^{2}\right\rangle\) undergoes the phase encoding and the quarter transformation, then the output state \(\left|{\psi }_{{{\mbox{out}}}}\right\rangle ={{{{{{{{\bf{U}}}}}}}}}_{q}{{{{{{{{\bf{U}}}}}}}}}_{\varphi }\left|{\psi }_{4}^{2}\right\rangle\) becomes
with c0 = c1 = c2 = c3, c4 = c5, c6 = c7, and c8 = c9. Then, one can obtain the two-photon detection probability of Pl(φ) = ∣cl∣2 (l = 0, 1, ..., 9), and it satisfies the normalizing condition ∑lPl(φ) = 1. Pl(φ) is a function of φ, and used for theoretical calculations in Fig. 3. See Supplementary Note 3 for the detailed information on Pl(φ) = ∣cl∣2.
Analysis considering experimental errors
In order to include the errors caused by our experimental imperfection, we consider the generation of 4-mode 2002 state \({\left|{\psi }_{4}^{2}\right\rangle }_{\exp }\) from an imperfect Bell state of \(|{{{\Phi }}}_{\exp }^{+}\rangle\), which is given by49
where \({\hat{\alpha }}_{j}^{H{{{\dagger}}} }\) (\({\hat{\alpha }}_{j}^{V{{{\dagger}}} }\)) creates a photon at αj with a horizontal (vertical) polarization shown in Fig. 2c, and \({\hat{\alpha }}_{j}^{^{\prime} H{{{\dagger}}} }\) (\({\hat{\alpha }}_{j}^{^{\prime} V{{{\dagger}}} }\)) creates a distinguishable photon with \({\hat{\alpha }}_{j}^{H{{{\dagger}}} }\) (\({\hat{\alpha }}_{j}^{V{{{\dagger}}} }\)) at same position due to experimental imperfection. Note that photons created from \({\hat{\alpha }}_{j}^{H{{{\dagger}}} }\) (\({\hat{\alpha }}_{j}^{V{{{\dagger}}} }\)) and \({\hat{\alpha }}_{j}^{^{\prime} H{{{\dagger}}} }\) (\({\hat{\alpha }}_{j}^{^{\prime} V{{{\dagger}}} }\)) are distinguishable and they do not interfere each other, and ϵ is real parameter varying from 0 to 1. Then after LBS, PBS, and HWP at 45∘, the experimentally prepared \({\left|{\psi }_{4}^{2}\right\rangle }_{\exp }\) becomes
where \({\hat{a}}_{i}^{{{{\dagger}}} }\) creates a photon at input port of quarter ai and \({\hat{a}}_{i}^{^{\prime} {{{\dagger}}} }\) creates a distinguishable photon due to experimental imperfection. Note that photons created from \({\hat{a}}_{i}^{{{{\dagger}}} }\) and \({\hat{a}}_{i}^{^{\prime} {{{\dagger}}} }\) are distinguishable and they do not interfere each other. After \({\left|{\psi }_{4}^{2}\right\rangle }_{\exp }\) undergoes the phase encoding Uφ and the experimentally reconstructed quarter transition matrix \({{{{{{{{\bf{U}}}}}}}}}_{q,\exp }\), then the output state becomes \({\left|{\psi }_{{{\mbox{out}}}}\right\rangle }_{\exp }={{{{{{{{\bf{U}}}}}}}}}_{q,\exp }{{{{{{{{\bf{U}}}}}}}}}_{\varphi }{\left|{\psi }_{4}^{2}\right\rangle }_{\exp }\). Here, \({{{{{{{{\bf{U}}}}}}}}}_{q,\exp }\) is given by
Then experimental two-photon detection probabilities are obtained as \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\), it is used for fitting curves in Fig. 4a. See Supplementary Note 3 for the detailed calculations on our error analysis.
Quantum Fisher information matrix and CRB
The QFIM is given by
By calculating FQ, one can obtain the QCRB from \(\,{{\mbox{Tr}}}\,[{{{{{{{{\bf{F}}}}}}}}}_{Q}^{-1}]\)14,16,17,39. For a \(\left|{\psi }_{4}^{2}\right\rangle\) probe state, the QCRB is calculated to be \(\,{{\mbox{Tr}}}\,[{{{{{{{{\bf{F}}}}}}}}}_{Q}^{-1}]\) = 1.5, and the minimum value of the CRB is obtained at φ1 ≃ 0.5π, φ2 = 0, and φ3 = 0.5π, and it saturates the QCRB with \(\,{{\mbox{Tr}}}\,[{{{{{{{{\bf{F}}}}}}}}}_{Q}^{-1}]\) = \(\,{{\mbox{Tr}}}\,[{{{{{{{{\bf{F}}}}}}}}}_{C}^{-1}({{{{{{{\boldsymbol{\varphi }}}}}}}})]\) = 1.5. The ideal FC(φ) and FQ matrices are obtained to be
where φ1 ≃ 0.5π, φ2 = 0, and φ3 = 0.5π.
The experimentally obtained \({{{{{{{{\bf{F}}}}}}}}}_{C}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\) matrix using \({P}_{l}^{\exp }({{{{{{{\boldsymbol{\varphi }}}}}}}})\) is given as below:
where φ1 ≃ 0.47π, φ2 = −0.07π, and φ3 = 0.52π. The experimentally obtained CRB is evaluated to be 1.85 ± 0.01. The detailed calculations are provided in Supplementary Note 3.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Code availability
The code used to generate the figures within this paper and other findings of this study are available from the corresponding author upon reasonable request.
Change history
06 December 2021
A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-27483-2
References
Moreau, P.-A., Toninelli, E., Gregory, T. & Padgett, M. J. Imaging with quantum states of light. Nat. Rev. Phys. 1, 367–380 (2019).
Brida, G., Genovese, M. & Berchera, I. R. Experimental realization of sub-shot-noise quantum imaging. Nat. Photon. 4, 227–230 (2010).
Guerrieri, F. et al. Sub-Rayleigh imaging via N-photon detection. Phys. Rev. Lett. 105, 163602 (2010).
Cui, J.-M., Sun, F.-W., Chen, X. D., Gong, Z.-J. & Guo, G.-C. Quantum statistical imaging of particles without restriction of the diffraction limit. Phys. Rev. Lett. 110, 153901 (2013).
Edgar, M. P. et al. Imaging high-dimensional spatial entanglement with a camera. Nat. Commun. 3, 984 (2012).
Schwartz, O. et al. Superresolution microscopy with quantum emitters. Nano Lett. 13, 5832–5836 (2013).
Ono, T., Okamoto, R. & Takeuchi, S. An entanglement-enhanced microscope. Nat. Commun. 4, 2426 (2013).
Taylor, M. A. & Bowen, W. P. Quantum metrology and its application in biology. Phys. Rep. 615, 1–59 (2016).
Schirhagl, R., Chang, K., Loretz, M. & Degen, C. L. Nitrogen-vacancy centers in diamond: nanoscale sensors for physics and biology. Annu. Rev. Phys. Chem. 65, 83–105 (2014).
Mauranyapin, N. P., Madsen, L. S., Taylor, M. A., Waleed, M. & Bowen, W. P. Evanescent single-molecule biosensing with quantum-limited precision. Nat. Photon. 11, 477–481 (2017).
Komar, P. et al. A quantum network of clocks. Nat. Phys. 10, 582–587 (2014).
Xia, Y. et al. Demonstration of a reconfigurable entangled radio-frequency photonic sensor network. Phys. Rev. Lett. 124, 150502 (2020).
Ge, W., Jacobs, K., Eldredge, Z., Gorshkov, A. V. & Foss-Feig, M. Distributed quantum metrology with linear networks and separable inputs. Phys. Rev. Lett. 121, 043604 (2018).
Humphreys, P. C., Barbieri, M., Datta, A. & Walmsley, I. A. Quantum enhanced multiple-phase estimation. Phys. Rev. Lett. 111, 070403 (2013).
Polino, E., Valeri, M., Spagnolo, N. & Sciarrino, F. Photonic quantum metrology. AVS Quantum Sci. 2, 024703 (2020).
Gessner, M., Pezzè, L. & Smerzi, A. Sensitivity bounds for multiparameter quantum metrology. Phys. Rev. Lett. 121, 130503 (2018).
Pezzè, L. et al. Optimal measurements for simultaneous quantum estimation of multiple phases. Phys. Rev. Lett. 119, 130504 (2017).
Goldberg, A. Z. et al. Multiphase estimation without a reference mode. Phys. Rev. A 102, 022230 (2020).
Gebhart, V., Smerzi, A. & Pezzè, L. Heisenberg-limited Bayesian multiphase estimation algorithm. Preprint at https://arxiv.org/abs/2010.09075 (2020).
Jachura, M., Chrapkiewicz, R., Demkowicz-Dobzański, R., Wasilewski, W. & Banaszek, K. Mode engineering for realistic quantum-enhanced interferometry. Nat. Commun. 7, 11411 (2016).
Vidrighin, M. D. et al. Joint estimation of phase and phase diffusion for quantum metrology. Nat. Commun. 5, 3532 (2014).
Roccia, E. et al. Multiparameter approach to quantum phase estimation with limited visibility. Optica 5, 1171–1176 (2018).
Liu, L. Z. et al. Distributed quantum phase estimation with entangled photons. Nat. Photon. 15, 137–142 (2021).
Polino, E. et al. Experimental multiphase estimation on a chip. Optica 6, 288–295 (2019).
Schnabel, R. Squeezed states of light and their applications in laser interferometers. Phys. Rep. 684, 1–51 (2017).
Lawrie, B. J., Lett, P. D., Marino, A. M. & Pooser, R. C. Quantum sensing with squeezed light. ACS Photon. 6, 1307–1318 (2019).
Guo, X. et al. Distributed quantum sensing in a continuous-variable entangled network. Nat. Phys. 16, 281–284 (2020).
Gessner, M., Smerzi, A. & Pezzé, L. Multiparameter squeezing for optimal quantum enhancements in sensor networks. Nat. Commun. 11, 3817 (2020).
Su, Z. et al. Multiphoton interference in quantum fourier transform circuits and applications to quantum metrology. Phys. Rev. Lett. 119, 080502 (2017).
Ciampini, M. A. et al. Quantum-enhanced multiparameter estimation in multiarm interferometers. Sci. Rep. 6, 28881 (2016).
Xiang, G. Y., Higgins, B. L., Berry, D. W., Wiseman, H. M. & Pryde, G. J. Entanglement-enhanced measurement of a completely unknown optical phase. Nat. Photon. 5, 43–47 (2011).
Zhang, L. & Chan, K. W. C. Quantum multiparameter estimation with generalized balanced multimode NOON-like states. Phys. Rev. A 95, 032321 (2017).
Zhang, L. & Chan, K. W. C. Scalable generation of multi-mode NOON states for quantum multiple-phase estimation. Sci. Rep. 8, 11440 (2018).
Dowling, J. P. Quantum optical metrology-the lowdown on high-N00N states. Contemp. Phys. 49, 125–143 (2008).
Slussarenko, S. et al. Unconditional violation of the shot-noise limit in photonic quantum metrology. Nat. Photon. 11, 700–703 (2017).
You, C. et al. Scalable multiphoton quantum metrology with neither pre- nor post-selected measurements. Preprint at https://arxiv.org/abs/2011.02454 (2020).
Peruzzo, A., Laing, A., Politi, A., Rudolph, T. & O’Brien, J. L. Multimode quantum interference of photons in multiport integrated devices. Nat. Commun. 2, 224 (2011).
Pryde, G. J. & White, A. G. Creation of maximally entangled photon-number states using optical fiber multiports. Phys. Rev. A 68, 052315 (2003).
Liu, J., Yuan, H., Lu, X. M. & Wang, X. Quantum Fisher information matrix and multiparameter estimation. J. Phys. A: Math. Theor. 53, 023001 (2019).
Edamatsu, K., Shimizu, R. & Itoh, T. Measurement of the photonic de Broglie wavelength of entangled photon pairs generated by spontaneous parametric down-conversion. Phys. Rev. Lett. 89, 213601 (2002).
D’Angelo, M., Chekhova, M. V. & Shih, Y. Two-Photon ––graphy. Phys. Rev. Lett. 87, 013602 (2001).
Boto, A. N. et al. Quantum interferometric optical lithography: exploiting entanglement to beat the diffraction limit. Phys. Rev. Lett. 85, 2733–3726 (2000).
Datta, A., Zhang, L., Dorner, U., Smith, B. J. & Walmsley, I. A. Quantum metrology with imperfect states and detectors. Phys. Rev. A 83, 063836 (2011).
Hu, X. M. et al. Beating the channel capacity limit for superdense coding with entangled ququarts. Sci. Adv. 4, eaat9304 (2018).
Wolfgramm, F., Vitelli, C., Beduini, F. A., Godbout, N. & Mitchell, M. W. Entanglement-enhanced probing of a delicate material system. Nat. Photon. 7, 28–32 (2013).
Morris, P. A., Aspden, R. S., Bell, J. E., Boyd, R. W. & Padgett, M. J. Imaging with a small number of photons. Nat. Commun. 6, 1–6 (2015).
Kim, T., Fiorentino, M. & Wong, F. N. Phase-stable source of polarization-entangled photons using a polarization Sagnac interferometer. Phys. Rev. A 73, 012316 (2006).
Hong, C. K., Ou, Z. Y. & Mandel, L. Measurement of subpicosecond time intervals between two photons by interference. Phys. Rev. Lett. 59, 2044 (1987).
Lyons, A. et al. Attosecond-resolution Hong-Ou-Mandel interferometry. Sci. Adv. 4, eaap9416 (2018).
Acknowledgements
This work was supported by the National Research Foundation of Korea (NRF) (2019M3E4A1079777, 2019M3E4A1078660, 2020M3E4A1079939, 2021R1C1C1003625), the Institute for Information and Communications Technology Promotion (IITP) (2020-0-00947, 2020-0-00972), and the KIST research program (2E31021).
Author information
Authors and Affiliations
Contributions
H.-T.L. initiated and led the project. The experimental scheme was designed by S.H., Y.-S.K., Y.-W.C., H.-T.L., and performed by S.H. Theoretical calculations were carried out by S.H., J.R., S.-W. L., and H.-T.L. S.H. and H.-T.L. analyzed the experimental data. H.J., S.M., and S.-W.H. helped to analyze the data and discuss the results. All authors contributed to discussions and writing the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hong, S., ur Rehman, J., Kim, YS. et al. Quantum enhanced multiple-phase estimation with multi-mode N00N states. Nat Commun 12, 5211 (2021). https://doi.org/10.1038/s41467-021-25451-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-021-25451-4
This article is cited by
-
Distributed quantum sensing of multiple phases with fewer photons
Nature Communications (2024)
-
Entanglement-enhanced optomechanical sensing
Nature Photonics (2023)
-
Quantum multiparameter estimation with multi-mode photon catalysis entangled squeezed state
Frontiers of Physics (2023)
-
Simultaneous multiple angular displacement estimation precision enhanced by the intramode correlation
Quantum Information Processing (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.