Quantum multiparameter estimation with multi-mode photon catalysis entangled squeezed state

Abstract

We propose a method to generate the multi-mode entangled catalysis squeezed vacuum states (MECSVS) by embedding the cross-Kerr nonlinear medium into the Mach−Zehnder interferometer. This method realizes the exchange of quantum states between different modes based on Fredkin gate. In addition, we study the MECSVS as the probe state of multi-arm optical interferometer to realize multi-phase simultaneous estimation. The results show that the quantum Cramer−Rao bound (QCRB) of phase estimation can be improved by increasing the number of catalytic photons or decreasing the transmissivity of the optical beam splitter using for photon catalysis. In addition, we also show that even if there is photon loss, the QCRB of our photon catalysis scheme is lower than that of the ideal entangled squeezed vacuum states (ESVS), which shows that by performing the photon catalytic operation is more robust against photon loss than that without the catalytic operation. The results here can find applications in quantum metrology for multiparatmeter estimation.

Graphical abstract

Keywords

quantum metrology / squeezing vacuum / quantum catalysis / multi-parameter estimation

1 Introduction

Quantum metrology is one of the most important research fields in quantum information science which can provide significant quantum advantages over its classical counterpart. A fundamental task in quantum metrology is to improve the estimation precision of parameters to be measured through quantum resources allowing by the basic principles of quantum mechanics. Generally speaking, the typical quantum metrology includes three steps: the preparation of probe states, the evolution of probe states, and the readout of the evolved states. In quantum metrology, the quantum Cramer−Rao bound (QCRB) is usually used to quantify the estimation precision offered by quantum metrology, which gives the lower limit of the estimation precision that can be achieved using any possible detection methods [1-5]. For this purpose, prior works are focused on the estimation of a single parameter with superior quantum resources [6-20]. For instance, by adopting nonclassical or entangled states, such as single (two)-mode squeezed vacuum state (S(T)SVS) [11-13] and NOON state [14], the estimation precision can overcome the so-called standard quantum limit (SQL) scaling as $1/\sqrt{N}$ with $N$ being the mean photon number of the probe state, and in certain cases, the precision can even approach to the renowned Heisenberg limit (HL) with a scaling $1/N$. If we can perform certain non-Gaussian operations such as photon addition, subtraction and catalysis on the squeezing states, we may further enhance the precision of the metrology [15-20].

In the realistic scenarios, such as biological system measurement [21-23], the optical imaging and the sensor network [24-28], multiparameter quantum metrology is indispensable and thus has received a lot of increasing interest in recent years [29-43], as the number of parameters affecting a physical process is usually more than one. For instance, Humphreys et al. [29] treated the phase imaging problem regarded as a multiparameter estimation process, and showed the advantages of the multiparameter simultaneous estimation using the multi-mode NOON state, when comparing with the independent estimation scheme. The advantages remain even if there are photon losses, as studied by Yue et al. [30]. Apart from the photon losses, Ho et al. [31] estimated three components of an external magnetic field using the entangled Greenberger−Horne−Zeilinger state including the dephasing noise and showed that its sensitivity can beat the SQL. Besides, Yao et al. [32] investigated the multiple phase estimation problem for a natural parametrization of arbitrary pure states under the white noise. In order to further develop the quantum enhanced multiparameter simultaneous estimation, Hong et al. [33] proposed a method to generate the multi-mode NOON state, and they experimentally demonstrated that the QCRB can be saturated using the multi-mode NOON state. In addition to the NOON state, entangled coherent state is also widely used in the field of quantum metrology. To compare the performance of the entangled coherent state [44, 45] can be better with that of the multi-mode NOON state, Liu et al. [34] proposed a theoretical scheme of quantum enhanced multiparameter metrology with generalized entangled coherent state and showed that the entangled coherent state can indeed give better precision than that of the multi-mode NOON state. After that, Zhang et al. [43] investigated the quantum multiparameter estimation with generalized balanced multimode NOON-like states, including the entangled squeezed vacuum state (ESVS), the entangled squeezed coherent state, the entangled coherent state, and the NOON state. Comparing with other multimode NOON-like states, they found that the ESVS has the lowest QCRB if the mean photon number is the same. However, how to prepare the multimode ESVS has not yet been studied before. In addition, the maximum average photon number of the squeezing vacuum state achievable currently is still limited (about 15 dB for single mode squeezed vacuum [46] and 10.7 dB for two-mode squeezed vacuum [47]). A natural question is that can we improve the multiparameter estimation further under the same average photon number by performing certain non-Gaussian operation on the ESVS?

In this paper, we propose a scheme to generate the multimode ESVS or the entangled superposition of the photon-catalysis squeezed vacuum state which may find applications in the quantum information and quantum metrology. Due to the fact that the multi-photon catalysis operation [48, 49] can improve the fidelity in quantum teleportation [50, 51], extend the transmission distance in continuous variable quantum key distribution [52, 53], and enhance the sensitivity of phase estimation for a single-phase estimation [19] and undo the noise effect of the channel [54], we also propose a scheme to improve the multiparameter estimation precision by using the MECSVS. Our results clearly show that the multi-photon catalytic operation can further improve the precision of phase estimation compared with the result with ordinary ESVS as the probe state under the same average photon number. Moreover, the usage of multi-photon catalysis in the multiparameter quantum metrology is also more robust against the photon losses which can find applications in the practical scenarios.

The paper is organized as follows. In Section 2, we propose a scheme to generate the MECSVS. In Sections 3 and 4, we evaluate the QCRB of multiparameter estimation with the symmetric MECSVS under ideal and photon-loss cases, respectively. In Section 5, we study the QCRB when the anti-symmetric MECSVS is used. Finally, we summarize the results.

2 The generation of the MECSVS

In this section, we propose a scheme to generate the MECSVS which is the probe quantum state used for the quantum multiparameter estimation. The schematic diagram is shown in Fig.1 in which three steps are included. In the first step, we perform the $n$-photon catalysis (orange box) on the input SSVS $|\xi ⟩$. In the $n$-photon catalysis, an $n$-photon Fock state $|n⟩$ in an ancillary mode $\stackrel{\sim }{a}$ is injected into one port of a beam splitter with a transmissivity $T,$ and then is detected at the corresponding output of mode $\stackrel{\sim }{a}$. In the meantime, the input squeezed state $|\xi ⟩$ is injected into the other port $a$ of the beam splitter. If $n$ photons are detected at the output port $\stackrel{\sim }{a}$, the so called $n$-photon catalysis is performed and the output state of port $a$ is given by $|{\xi }^{\prime }⟩$ which is the single-mode $n$-photon catalyzed squeezed state. This multi-photon catalyzed process can be regarded as an effective operator [50], i.e.,

Fig.1 Generation of the $d$-mode MECSVS using quantum-optical Fredkin gate. The orange box is the multi-photon catalysis process, in which an $n$-photon Fock state is inputted into auxiliary mode of beam splitter, conditional measuring $n$ photons at one output. The dashed box is the entanglement generator, which creates two-mode entangled catalysis squeezed vacuum state by using the single-mode multi-photon catalysis squeezed vacuum state and vacuum state. Applying the entanglement generator $x$ times with each input of the following generator aligning with one of the output modes of the previous generator creates a $d$-mode MECSVS.

Full size|PPT slide

(1)$\begin{array}{rl}\hat{O}& \equiv ⟨n|\hat{B}\left(T\right)|n⟩\\ & =\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\{\frac{\exp \left[\mu (\tau )a^{†}a\right]}{1-\tau }\}{|}_{\tau =0},\end{array}$

where $\hat{B}\left(T\right)={\mathrm{e}}^{({\tilde{a}}^{†}a-\tilde{a}{a}^{†})\arccos \sqrt{T}}$ and the parameter

(2)$\begin{array}{r}\mu (\tau )=\ln \frac{\sqrt{T}-\tau /\sqrt{T}}{1-\tau }.\end{array}$

Considering that the input state of $a$-mode is a Gaussian SSVS $|\xi ⟩=\sqrt{\mathrm{s}\mathrm{e}\mathrm{c}\mathrm{h}r}\sum _{l=0}^{\mathrm{\infty }}{\left(-\frac{1}{2}\tanh r\right)}^l\sqrt{(2l)!}/l!|2l⟩$ as the input state [55], the corresponding single-mode multiphoton catalysis squeezed vacuum state is thus given by

(3)$\begin{array}{rl}|{\xi }^{\mathrm{\prime }}⟩& =\frac{1}{\sqrt{\mathbb{N}}}\hat{O}{|\xi ⟩}_a\\ & ={\frac{1}{\sqrt{\mathbb{N}}}\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}[\beta (\tau )\sum _{l=0}^{\mathrm{\infty }}{c}_l(\tau ){|2l⟩}_a]|}_{\tau =0},\end{array}$

with $\beta (\tau )=1/\left(1-\tau \right)\sqrt{\cosh r},$ $c_l(\tau )={\left[-\frac{1}{2}{\mathrm{e}}^{2\mu (\tau )}\tanh r\right]}^l$ $\cdot \sqrt{\left(2l\right)!}/l!$ and $\mathbb{N}$ is the success probability given by

(4)$\begin{array}{r}\mathbb{N}={\frac{T^n}{{\left(n!\right)}^2}\frac{{\mathrm{\partial }}^{2n}}{\mathrm{\partial }{\tau }^n\mathrm{\partial }{\tau }^{\ast n}}\mathrm{\aleph }\left(\tau ,{\tau }^{\ast }\right)[1-\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right){]}^{-\frac{1}{2}}|}_{\tau ={\tau }^{\ast }=0},\end{array}$

with

(5)$\begin{array}{rl}& \mathrm{\aleph }\left(\tau ,{\tau }^{\ast }\right)=\frac{1}{\left(1-\tau \right)\left(1-{\tau }^{\ast }\right)\cosh r},\\ & \mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)={\mathrm{e}}^{2[\mu ({\tau }^{\ast })+\mu (\tau )]}{\tanh }^{2}r,\end{array}$

where ${\tau }^{\ast }$ is the complex conjugate of $\tau$. In particular, from Eq. (3), when $T=1$, the single-mode multi-photon catalysis squeezed vacuum state is reduced to the SSVS, as expected. Although the photon catalysis operation is probabilistic [56], it can be heralded and we choose to perform the quantum metrology only when the MECSVS is successfully generated.

After generating the single-mode multi-photon catalysis squeezed vacuum state as shown in Eq. (3), we next produce a two-mode multi-photon catalysis squeezed vacuum state using step 2 as shown in Fig.1. For this purpose, we propose to employ that quantum-optical Fredkin gate which consists of two Mach−Zehnder interferometers mediated by a cross-Kerr medium. To be more specific, the single-mode multi-photon catalysis squeezed vacuum states and a vacuum state $|0⟩$ are respectively injected into a cross-Kerr medium based Mach−Zehnder interferometer (MZI-2) from modes $a$ and $b$. Simultaneously, a vacuum state ${|0⟩}_u$ and a single photon state ${|1⟩}_v$ are used as inputs of the MZI-1. We assume that the four beam splitters ${\hat{B}}_j$ $(j=1,2,3,4)$ are chosen as 50:50, i.e., ${\hat{B}}_1={\mathrm{e}}^{\mathrm{i}\pi \left(u^{†}v+u{v}^{†}\right)/4}$ and ${\hat{B}}_3=$ ${\mathrm{e}}^{\mathrm{i}\pi \left(a^{†}b+a{b}^{†}\right)/4}$ with ${\hat{B}}_4={\hat{B}}_3^{†}$ and ${\hat{B}}_2={\hat{B}}_1^{†}$. After passing through the beams splitters $B_1$ and $B_3$, the photons in mode $a$ and mode $u$ pass through the corss-Kerr medium at the same time and the effective operator is given by [57, 58]

(6)$\begin{array}{r}{\hat{U}}_k=\exp \left(\mathrm{i}\chi t{a}^{†}a{u}^{†}u\right),\end{array}$

where $\chi$ is the nonlinear Kerr coupling coefficient and $t$ is the interaction time. For our purpose here, we choose $\chi t=\pi$. Unfortunately, the cross-Kerr nonlinearity in the natural medium is usually very small and $\chi t$ is usually much less than $\pi$. However, enhancement of the cross-Kerr nonlinearity is not impossible. In the past few decades, a number of methods have been proposed to achieve giant Kerr nonlinearity. For example, giant Kerr-nonlinearity such that $\chi >2\pi \times {10}^{10}\mathrm{H}\mathrm{z}$ can be achieved in an atomic ensemble by the electromagnetic induced transparency (EIT) and the slow light effect [59-63]. Besides, in Ref. [64], the author showed that by constructing a one-dimensional nonlinear photonic crystal from alternating layers of Kerr medium and linear dielectric medium, the phase of the wave function of the incident photons can be rotated by $\pi$ phase. In addition, by measurement-induced quantum operations, Costanzo et al. [65] demonstrated an experimental implementation of a strong Kerr nonlinearity where a $\pi$ phase shift is realized. Therefore, $\pi$ phase shift by the Kerr-type interaction is possible.

It is clearly seen that if the $u$-mode is vacuum, no phase shift for the $b$ mode. However, if the $u$-mode has 1 photon, there is a phase shift for the $b$ mode. Due to the phase shift, the $b$-mode may enter the upper path or lower path after passing through $B_4$. This is the basic principle for generating the entangled state in this setup. The combination of the operation ${\hat{U}}_F={\hat{B}}_4{\hat{U}}_k{\hat{B}}_3$ is actually the quantum-optical Fredkin gate which is given by

(7)$\begin{array}{r}{\hat{U}}_F={\mathrm{e}}^{\mathrm{i}\frac{\pi }{2}{u}^{†}u\left(a^{†}a+b^{†}b\right)}{\mathrm{e}}^{\frac{\pi }{2}{u}^{†}u\left(a{b}^{†}-a^{†}b\right)}.\end{array}$

This quantum gate can effectively entangle the two photon modes.

To be more specific, when $|{\xi }^{\mathrm{\prime }}⟩$ and $|0⟩$ are injected into the MZI-2 (${|0⟩}_u$ and ${|1⟩}_v$ are injected into the MZI-1), the unnormalized output state can be expressed as

(8)$\begin{array}{rl}& {\hat{B}}_2{\hat{B}}_4{\hat{U}}_k{\hat{B}}_3{\hat{B}}_1{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v\\ =& \frac{1}{2}[\left({|0⟩}_a{|{\mathrm{e}}^{\mathrm{i}\pi l}{\xi }^{\mathrm{\prime }}⟩}_b+{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b\right){|1⟩}_u{|0⟩}_v\\ & +\mathrm{i}\left(-{|0⟩}_a{|{\mathrm{e}}^{\mathrm{i}\pi l}{\xi }^{\mathrm{\prime }}⟩}_b+{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b\right){|0⟩}_u{|1⟩}_v].\end{array}$

Here, the state ${|{\mathrm{e}}^{\mathrm{i}\pi l}{\xi }^{\mathrm{\prime }}⟩}_b$ describes the state similar to Eq. (3) but the coefficient $c_l(\tau )$ is multiplied by a phase factor ${\mathrm{e}}^{\mathrm{i}\pi l}$. To remove this additional phase shift, we place a quarter wave plate in the output route of $b$ mode. When $2l$ photons pass through this quarter wave plate, an additional $\pi l$ phase shift is accumulated which exactly cancels out the previous phase factor. Two single photon detectors $D_1$ and $D_2$ are placed in the output routes of $u$ and $v$ modes. According to Eq. (8), the output state depends on the detection results of $u$ and $v$ modes. Finally, the output state is given by

(9)$\begin{array}{r}{|\psi ⟩}_{ab}=\frac{1}{\sqrt{{\tilde{N}}_2}}(\pm {|0⟩}_a{|{\xi }^{\mathrm{\prime }}⟩}_b+{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b),\end{array}$

where ${\tilde{N}}_2$ is a normalization factor. The symmetric (antisymmetric) state is obtained when a single photon is detected in $D_1$ ($D_2$) and no photon is detected in $D_2$ ($D_1$). It is clearly seen that the two-mode ECSVS is generated. Both the symmetric and the antisymmetric ECSVSs can be used for improving the precision of the metrology. In the following, we mainly consider the symmetric case and disscuss the antisymmetric case in Section 5.

After preparing the two-mode ECSVS state, we use it as inputs of two MZIs and repeat the procedures as those in step 2. By repeating these procedures for a number of time, we can in principle generate the MECSVS. If all the detection results of the auxiliar qubits are one photon in the mode $u$ and zero photon in the mode $v$, the output state is then given by

(10)$\begin{array}{r}|\mathrm{\Psi }⟩=\frac{1}{\sqrt{\tilde{\mathbb{N}}}}\sum _{j=0}^d{|0⟩}_0{|0⟩}_1{|0⟩}_2{|0⟩}_3...{|{\xi }^{\mathrm{\prime }}⟩}_j...{|0⟩}_d,\end{array}$

where $\tilde{\mathbb{N}}=\left(d+1\right)\left(1+d\left|⟨0\right|{{\xi }^{\mathrm{\prime }}⟩|}^2\right)$ is the normalized coefficient, which can be calculated as

(11)$\begin{array}{rl}\tilde{\mathbb{N}}=& \frac{T^n}{{\left(n!\right)}^2}\frac{{\mathrm{\partial }}^{2n}}{\mathrm{\partial }{\tau }^n\mathrm{\partial }{\tau }^{\ast n}}\left(d+1\right)\mathrm{\aleph }\left(\tau ,{\tau }^{\ast }\right)\\ & {\times \left\{{\left[1-\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)\right]}^{-1/2}+d\right\}|}_{\tau ={\tau }^{\ast }=0}.\end{array}$

We note that $\tilde{\mathbb{N}}\ge d+1$ because $\left|⟨0\right|{{\xi }^{\mathrm{\prime }}⟩|}^2\ge 0$. Therefore, the average photon number of the symmetric superposition of the MECSVS $\overline{n}=(d+1)⟨{\xi }^{\mathrm{\prime }}|a^{†}a|{\xi }^{\mathrm{\prime }}⟩/\tilde{\mathbb{N}}$ is usually smaller than the average photon number of the single mode catalysis squeezed vacuum state $⟨{\xi }^{\mathrm{\prime }}|a^{†}a|{\xi }^{\mathrm{\prime }}⟩$. When ${\xi }^{\mathrm{\prime }}$ is the squeezed vacuum, we can see that the average photon number of the symmetric superposition of ESVS is usually smaller than that of the single-mode SVS, i.e., ${\sinh }^{2}r$. The quantum state shown in Eq. (10) is a symmetric MECSVS which is a highly entangled state. In particular, when $T=1$, the MECSVS is reduced to the multi-mode ESVS, as expected. In the following, we shall use the MECSVS as the inputs of a multi-arm interferometer in order to effectively improve the precision of multiparameter estimations of multiple optical phases at the same time with and without photon losses.

3 The QCRB of multiparameter estimation without photon losses

In this section, we investigate the QCRB of multiple independent phases estimation simultaneously by a multi-arm interferometer with inputting the MECSVS shown in Eq. (10) under the ideal case. The schematic diagram of $d$ parameters estimation is shown in Fig.2 where the mode $0$ is the reference beam and the modes from $1$ to $d$ are the parametric beams. The $d$ estimated parameters are assumed to be mutually independent linear phases and their transformations through the multi-arm interferometer can be represented by a unitary operator [29, 43]

Fig.2 Schematic of multiparameter estimation with d phase shifts under an ideal case.

Full size|PPT slide

(12)$\begin{array}{r}{\hat{U}}_{\theta }=\exp \left(\mathrm{i}\sum _{j=1}^d{\theta }_j\cdot {\hat{N}}_j\right),\end{array}$

where ${\theta }_j$ and ${\hat{N}}_j=a_j^{†}{a}_j$ represent the phase shift and the photon number operator for the $j$th mode, respectively. Since all the phases are independent to each other, the photon number operators ${\hat{N}}_j$ for different $j$ are commutable. When the input MECSVS $|\mathrm{\Psi }⟩$ goes through the interferometer, the output state can be expressed as $|{\mathrm{\Psi }}_{out}⟩={\hat{U}}_{\theta }|\mathrm{\Psi }⟩$. According to the definition of QCRB, the estimation precision of $\hat{\theta }$ is inversely proportional to the quantum Fisher inforamtion (QFI) of the output state $|{\mathrm{\Psi }}_{out}⟩$, i.e.,

(13)$\begin{array}{r}{\left|\mathrm{\Delta }\hat{\theta }\right|}^2\ge {\left|\mathrm{\Delta }\hat{\theta }\right|}_{QCRB}^2=\text{Tr}[{\mathit{F}}_{\theta }^{-1}],\end{array}$

where $|\mathrm{\Delta }\hat{\theta }{|}^2$ is the total variance as a figure of merit for the multiparameter estimation [34], Tr$[\cdot ]$ represents the trace operation and ${\mathit{F}}_{\theta }^{-1}$ is the inverse of the $d\times d$ QFI matrix. Given an arbitrary pure state, it is possible to saturate the QCRB if $\text{Im}⟨\psi |L_j{L}_k|\psi ⟩=0$ is satisfied for all $j,k$ and $\theta ,$ in which $L_{j(k)}$ is the symmetric logarithmic derivative given by $L_{j(k)}$ = $2(|{\mathrm{\partial }}_{j(k)}{\mathrm{\Psi }}_{out}⟩⟨{\mathrm{\Psi }}_{out}|+|{\mathrm{\Psi }}_{out}⟩$ $⟨{\mathrm{\partial }}_{j(k)}{\mathrm{\Psi }}_{out}|)$ with $|{\mathrm{\partial }}_{j(k)}{\mathrm{\Psi }}_{out}⟩$ = ${\mathrm{\partial }}_{j(k)}|{\mathrm{\Psi }}_{out}⟩/{\mathrm{\partial }}_{j(k)}{\theta }_{j(k)}$ [66, 67]. For our scheme, based on Eqs. (10) and (12), the elements of the QFI matrix are given by

(14)$\begin{array}{r}{f}_{jk}=4(⟨\mathrm{\Psi }|{\hat{N}}_j{\hat{N}}_k|\mathrm{\Psi }⟩-⟨\mathrm{\Psi }|{\hat{N}}_j|\mathrm{\Psi }⟩⟨\mathrm{\Psi }|{\hat{N}}_k|\mathrm{\Psi }⟩).\end{array}$

As a result, the QFI matrix reads as [68]

(15)$\begin{array}{r}{\mathit{F}}_{\theta }=4\frac{1}{\tilde{\mathbb{N}}}⟨{\xi }^{\mathrm{\prime }}|{\hat{N}}^2|{\xi }^{\mathrm{\prime }}⟩I-4\frac{1}{{\tilde{\mathbb{N}}}^2}⟨{\xi }^{\mathrm{\prime }}|\hat{N}|{\xi }^{\mathrm{\prime }}{⟩}^2\overline{I},\end{array}$

where I represents the identity matrix and $\overline{I}$ denotes the matrix with the elements ${\overline{I}}_{jk}=1$ for all $j$ and $k$. From Eqs. (13) and (15), we can finally obtain the expression of the QCRB for our scheme [34], i.e.,

(16)$\begin{array}{r}|\mathrm{\Delta }\hat{\theta }{|}_{QCRB}^2=\frac{d}{4⟨{\xi }^{\mathrm{\prime }}|{\hat{N}}^2|{\xi }^{\mathrm{\prime }}⟩}\left(\tilde{\mathbb{N}}+\frac{1}{\mathbb{k}-d/\tilde{\mathbb{N}}}\right),\end{array}$

where $\mathbb{k}=⟨{\xi }^{\mathrm{\prime }}|{\hat{N}}^2|{\xi }^{\mathrm{\prime }}⟩/⟨{\xi }^{\mathrm{\prime }}|\hat{N}|{\xi }^{\mathrm{\prime }}{⟩}^2$ with

(17)$\begin{array}{rl}⟨{\xi }^{\mathrm{\prime }}|{\hat{N}}^2|{\xi }^{\mathrm{\prime }}⟩=& \frac{T^n}{{\left(n!\right)}^2}\frac{{\mathrm{\partial }}^{2n}}{\mathrm{\partial }{\tau }^n\mathrm{\partial }{\tau }^{\ast n}}\mathrm{\aleph }\left(\tau ,{\tau }^{\ast }\right)\\ & \times \frac{\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)\left[\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)+2\right]}{{\left[1-\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)\right]}^{5/2}}{|}_{\tau ={\tau }^{\ast }=0},\\ ⟨{\xi }^{\mathrm{\prime }}|\hat{N}|{\xi }^{\mathrm{\prime }}⟩=& \frac{T^n}{{\left(n!\right)}^2}\frac{{\mathrm{\partial }}^{2n}}{\mathrm{\partial }{\tau }^n\mathrm{\partial }{\tau }^{\ast n}}\mathrm{\aleph }\left(\tau ,{\tau }^{\ast }\right)\\ & \times \frac{\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)}{{\left[1-\mathrm{\Im }\left(\tau ,{\tau }^{\ast }\right)\right]}^{3/2}}{|}_{\tau ={\tau }^{\ast }=0}.\end{array}$

Note that $|\mathrm{\Delta }\hat{\theta }{|}_{QCRB}^2$ is positive definite according to Eq. (16). Especially, when $T=1$, one can obtain $⟨{\xi }^{\mathrm{\prime }}|{\hat{N}}^2|{\xi }^{\mathrm{\prime }}⟩={\sinh }^{4}r+2{\sinh }^{2}r{\cosh }^{2}r$ and $⟨{\xi }^{\mathrm{\prime }}|\hat{N}|{\xi }^{\mathrm{\prime }}⟩={\sinh }^{2}r$, which is consistent with the result of the $d+1$ modes ESVS case, as expected [43]. When $T\ne 1$, the QCRB can be different from that of the ESVS case.

To elaborate the advantages of the MECSVS as inputs of the multi-arm interferometer, we plot the $|\mathrm{\Delta }\hat{\theta }{|}_{QCRB}^2$ as a function of the mean photon number for several different catalytic photon numbers $n\in \{1,2,3\}$, as shown in Fig.3 where $d=5$ and $T=0.7,0.8,0.9$. As a comparison, the black solid line corresponds to the result using the multi-mode ESVS case. It is clearly seen that the QCRBs in the case using MECSVS as inputs are obviously lower than those using the normal ESVS for all three catalytic photon numbers (i.e., $n=1,2$ or $3$) especially when the average photon number is small. This indicates that by catalyzing the SSVS before inputting into the interferometer we can significantly improve the phase detection precision. The QCRB is the lowest when $T=0.9$ comparing with $T=0.7$ and $T=0.8$ for all three catalytic photon numbers. We also note that the mean photon number of the MECSVS increases with the squeezing parameter $r$ but is saturated when $r$ is large enough (see Fig.4). The maximum reachable mean photon number increases when T or the catalysis photon number increases. In Fig.3 we mainly show the improvement of precision in the range of low photon number. This may find applications in the situation where the disturbance to the system needs to be minimized when the phases are measured. Despite that we mainly show the improvement in the low photon number range, the improvement when the photon number is relatively large can also be achieved if the transmissivity $T$ or the number of photon catalysis increase (see Fig.4) or the antisymmetric state is used for the metrology (see Fig.5 in Section 5).

Fig.3 The QCRB as a function of mean photon number for the ESVS (black) and the MECSVSs with $T=0.9$ (blue), $T=0.8$ (red), and $T=0.7$ (green). (a) The case of single-photon catalytic, (b) the case of two-photon catalytic, (c) the case of three-photon catalytic.

Full size|PPT slide

Fig.4 The mean photon number as a function of mean photon number for the MECSVSs with $T=0.9$ (blue), $T=0.8$ (red), and $T=0.7$ (green). (a) the case of single-photon catalytic, (b) the case of two-photon catalytic, (c) the case of three-photon catalytic. The corresponding mean photon number of the SSMS is also plotted (see black line).

Full size|PPT slide

Fig.5 The QCRB as a function of the mean photon number (a) for given $T=0.9$ and different catalysis photon number $n=1,2,3$ in ideal case and (b) for given $T=0.9$ and $\eta =0.9$ with photon losses.

Full size|PPT slide

4 The QCRB of the multiparameter estimation with photon losses

Under the realistic environment, photon losses and phase diffusions are usually unadvoidable, which can affect the ultimate precision limit of the optical interferometry. Here we focus on studying the QCRB of multiparameter estimation with photon losses, as shown in Fig.6. To describe the influences of noise environment on the phase probing system $S$, additional degrees of freedom should be introduced, i.e., the environment denoted as E. Here we assume that the environment is in the vacuum state which is a reasonable assumption for the optical regime. In such circumstance, the initial state of the system and the environment are separable and can be written as ${|\mathrm{\Psi }⟩}_S{|\mathbf{0}⟩}_E$. The evolution of the whole combined system is unitary and can be denoted as ${\hat{U}}^{S+E}\left(\theta \right)$. Thus, the final state before the measurement can be expressed as [69]

Fig.6 Schematic of multiparameter estimation with d phase shifts under the photon losses.

Full size|PPT slide

(18)$\begin{array}{r}{|\mathrm{\Psi }\left(\theta \right)⟩}_{S+E}\text{=}{\hat{U}}^{S+E}\left(\theta \right){|\mathrm{\Psi }⟩}_S{|\mathbf{0}⟩}_E\text{=}\sum _{\mathit{k}}{\hat{\mathrm{\Pi }}}_{\mathit{k}}\left(\theta \right){|\mathrm{\Psi }⟩}_S{|\mathit{k}⟩}_E,\end{array}$

where ${\hat{U}}^{S+E}\left(\theta \right)$ = $U_0^{S_0+E_0}{U}_1^{S_1+E_1}\left({\theta }_1\right)$....$U_d^{S_d+E_d}\left({\theta }_d\right)$ and the environment vacuum state ${|\mathbf{0}⟩}_E$ = ${|\mathbf{0}⟩}_{E_0}{|\mathbf{0}⟩}_{E_1}$....${|\mathbf{0}⟩}_{E_d}$. In the second identity, ${|\mathit{k}⟩}_E$ = ${|k_0⟩}_{E_0}{|k_1⟩}_{E_1}$...${|k_d⟩}_{E_d}$ are orthonormal states of the environment, and ${\hat{\mathrm{\Pi }}}_{\mathit{k}}\left(\theta \right)$ = ${\hat{\mathrm{\Pi }}}_{k_0}\left({\theta }_0\right)\otimes {\hat{\mathrm{\Pi }}}_{k_1}\left({\theta }_1\right)$....$\otimes {\hat{\mathrm{\Pi }}}_{k_d}\left({\theta }_d\right)$ is the tensor product of all Kraus operator ${\hat{\mathrm{\Pi }}}_{k_j}\left({\theta }_j\right)$ [70, 71], each of which is defined as

(19)$\begin{array}{r}{\hat{\mathrm{\Pi }}}_{k_l}\left({\theta }_l\right){\text{=}}_{E_l}⟨k|{\hat{U}}_l^{S_l+E_l}\left({\theta }_l\right){|0⟩}_{E_l}.\end{array}$

From Eq. (18), we can then calculate the QFI matrix of the whole system including the noisy environment and the QFI matrix elements is given by

(20)$\begin{array}{r}ϝ(\theta ,{\hat{\mathrm{\Pi }}}_{\mathit{k}}\left(\theta \right){)}_{lj}=4(⟨{\hat{\pounds }}^{\left(lj\right)}⟩-⟨{\hat{\mathrm{\Gamma }}}^{\left(l\right)}⟩⟨{\hat{\mathrm{\Gamma }}}^{\left(j\right)}⟩),\end{array}$

where $⟨\cdot ⟩$ standing for $⟨\mathrm{\Psi }|\cdot |\mathrm{\Psi }⟩$ and

(21)$\begin{array}{r}{\hat{\mathrm{\Gamma }}}^{\left(l\right)}=\sum _{k_l}\mathrm{i}\frac{\mathrm{d}{\hat{\mathrm{\Pi }}}_{k_l}^{†}\left({\theta }_l\right)}{\mathrm{d}{\theta }_l}{\hat{\mathrm{\Pi }}}_{k_l}\left({\theta }_l\right),\end{array}$

(22)$\begin{array}{r}{\widehat{\mathit{\text{£}}}}^{\left(lj\right)}=\{\begin{array}{c}\sum _{k_l}\frac{\mathrm{d}{\hat{\mathrm{\Pi }}}_{k_l}^{†}\left({\theta }_l\right)}{\mathrm{d}{\theta }_l}\frac{\mathrm{d}{\hat{\mathrm{\Pi }}}_{k_l}\left({\theta }_l\right)}{\mathrm{d}{\theta }_l},l=j,\\ {\hat{\mathrm{\Gamma }}}^{\left(l\right)}{\hat{\mathrm{\Gamma }}}^{\left(j\right)},l\ne j.\end{array}\end{array}$

As an example of the noise environment, the photon loss can be simulated using an optical beam splitter with a transmissivity ${\eta }_l$ (${\eta }_l=1$ corresponds to lossless case, and ${\eta }_l=0$ corresponds to complete photon loss). From Eqs. (19), (21) and (22), we can obtain ${\hat{\mathrm{\Gamma }}}^{\left(l\right)}={\mu }_l{\hat{N}}_l$ and ${\widehat{\mathit{\text{£}}}}^{\left(ll\right)}={\mu }_l{\hat{N}}_l^2+v_l{\hat{N}}_l$ when $l=j$ (and ${\widehat{\mathit{\text{£}}}}^{\left(lj\right)}={\hat{\mathrm{\Gamma }}}^{\left(l\right)}{\hat{\mathrm{\Gamma }}}^{\left(j\right)}$ when $l\ne j$), with ${\mu }_l=1-\left(1+{\epsilon }_l\right)\left(1-{\eta }_l\right)$ and $v_l=(1+$ ${\epsilon }_l{)}^2{\eta }_l\left(1-{\eta }_l\right)$. As shown in Fig.6, if such photon losses exist in the each mode of the multi-arm interferometer, a series of Kraus operators in each mode is given by [30]

(23)$\begin{array}{r}{\hat{\mathrm{\Pi }}}_{k_l}\left({\theta }_l\right)=\sqrt{\frac{{\left(1-{\eta }_l\right)}^{k_l}}{k_l}}{\mathrm{e}}^{\mathrm{i}{\theta }_l\left({\hat{N}}_l-{\epsilon }_l{k}_l\right)}{\eta }_l^{\frac{{\hat{N}}_l}{2}}{a}_l^{k_l},\end{array}$

where ${\epsilon }_l$ is an arbitrary real number. ${\epsilon }_l$ need to be optimized to make the lower bound $ϝ(\theta ,{\hat{\mathrm{\Pi }}}_{\mathit{k}})$ as tight as possible. According to Eqs. (20) and (23), the lower bound for the optimal precision of multiparameter estimation is given by [72, 73]

(24)$\begin{array}{r}\text{cov}\left(\theta \right)\ge \frac{1}{ϝ(\theta ,{\hat{\mathrm{\Pi }}}_{\mathit{k}})}.\end{array}$

Based on Eq. (24), let us take the MECSVS as the input of multi-arm interferometer. In this case, all (off-) diagonal elements of $ϝ(\theta ,{\hat{\mathrm{\Pi }}}_{\mathit{k}})$ are the same, denoted as ${ϝ}_d$ (${ϝ}_o$). For simplicity, here we can make a reasonable assumption that ${\eta }_l=\eta$ and ${\epsilon }_l=\epsilon$ since all modes are symmetric for the probe state and then we can obtain

(25)$\begin{array}{l}{ϝ}_d=4({\mu }_l^2⟨\mathrm{\Delta }{\hat{N}}_l^2⟩+v_l⟨{\hat{N}}_l⟩),\end{array}$

(26)$\begin{array}{l}{ϝ}_o=4[{\mu }_l^2(⟨{\hat{N}}_l{\hat{N}}_j⟩-⟨{\hat{N}}_l⟩⟨{\hat{N}}_j⟩)],\end{array}$

where $⟨\mathrm{\Delta }{\hat{N}}_l^2⟩=⟨{\hat{N}}_l^2⟩-⟨{\hat{N}}_l{⟩}^2$ and $⟨{\hat{N}}_l{\hat{N}}_j⟩=0,$ $⟨{\hat{N}}_l⟩⟨{\hat{N}}_j⟩=$ $⟨{\hat{N}}_l{⟩}^2$ for the MECSVS. Then we can calculate the analytical expression of the QCRB with MECSVS as inputs of the multi-arm interferometer under the photon losses [30]

(27)$\begin{array}{r}{\left|\mathrm{\Delta }\hat{\theta }\right|}_{QCRB}^2=\underset{{\hat{\mathrm{\Pi }}}_{\mathit{k}}}{max}\text{Tr}[{ϝ}^{-1}(\theta ,{\hat{\mathrm{\Pi }}}_{\mathit{k}})],\end{array}$

where

(28)$\begin{array}{r}\text{Tr}[{ϝ}^{-1}(\theta ,{\hat{\mathrm{\Pi }}}_{\mathit{k}})]\text{=}\frac{d-1}{{ϝ}_d-{ϝ}_o}+\frac{1}{{ϝ}_d+\left(d-1\right){ϝ}_o}.\end{array}$

In Fig.7(a), we compare the QCRBs of the ideal (solid lines) and the photon-loss (dashed lines) cases with different input states (i.e., the ESVS and the single-photon catalyzed state) where the strength of the photon loss is chosen to be $\eta =0.9$. We can see that with photon loss, the QCRBs increase for both inputs. However, we can see that the QCRB of using the single-photon catalyzed state is still lower than that using the ESVS case which indicates that the single-photon catalyzed state still has better performance than the ESVS case under the photon losses. We also find that when the mean photon number is less than $0.32$, the QCRB for the MECSVS with $n=1$ with photon losses can be still lower than that for the multi-mode ESVS without photon losses. In addition, we also consider the QCRB with the input MECSVS for several different catalytic photon numbers $n\in \{1,2,3\},$ as shown in Fig.7(b) where $\eta =0.9$. It is found that, the QCRBs for all the three catalysis photon numbers are about the same with photon losses and all of them perform better than that with the ESVS input.

Fig.7 In the presence of photon losses, the phase uncertainty $|\mathrm{\Delta }\hat{\theta }{|}_{QCRB}^2$ as a function of the mean photon number. (a) Comparison of the ideal case and the loss case for the ESVS and the single-photon catalysis operation. (b) Comparison of the phase uncertainty $|\mathrm{\Delta }\hat{\theta }{|}_{QCRB}^2$ with the number of the photon catalysis $n=1,2,$ and $3$, and the ESVS (black dashed line) for 10% photon loss ($\eta$ $=0.9$).

Full size|PPT slide

We also investigate the QCRB for different strength of the photon loss which are shown in Fig.8 where we set $T=0.9$. From the three subfigures, we can see that with the decrease of $\eta \in \{1,0.98,0.96,0.94,0.92,0.90\},$ the QCRB for the MECSVS increases, meaning that the precision of multiparameter estimation decreases with the increase of photon-loss intensity. More interestingly, wih the increase of $n=1,2,3$, the usage of the MECSVS with respect to the QCRB canbe superior to that of the multi-mode ESVS at a certain range of the mean photon number. To be more specific, for $\eta =0.94$ (green line), when the mean photon number is less than $0.72,0.90$ and $0.95$, the QCRB for the MECSVS of $n=1,2,3$ with photon losses can be lower than that for the multi-mode ESVS without photon losses, and such a phenomenon is also true when $\eta =0.96$ and $0.98$. Thus, the results above clearly show that using the MECSVS can improve the phase estimation precision over the ESVS in both the noisy and noisy-free cases.

Fig.8 The QCRB as the function of the mean photon number in the presence of photon losses. The dot line represents the ideal ESVS, the dotdash line represents the single-photon catalytic ESVS, and the colored line represents the QCRB corresponding to different dissipation factors. (a) The case of single-photon catalysis, (b) the case of two-photon catalysis, (c) the case of three-photon catalysis.

Full size|PPT slide

5 The QCRB of the other multi-mode entangled states

In Section 2, we show that if the measurement results for $D_1$ and $D_2$ detectors are zero and one photon, respectively, antisymmetric instead of symmetric superposition of vacuum and ESVS is obtained. In this section, we show that this antisymmetric state can be also used for improving the phase estimation precision. As an example, we assume that all the $D_1$ and $D_2$ detectors detect zero and one photon, respectively. In this case, the output state for the 6 modes example (i.e., $d=5$) can be calculated as

(29)$\begin{array}{rl}|{\mathrm{\Psi }}^{\mathrm{\prime }}⟩=& \frac{1}{\sqrt{{\tilde{\mathbb{N}}}^{\mathrm{\prime }}}}({|\xi ⟩}_0{|0⟩}_1{|0⟩}_2{|0⟩}_3{|0⟩}_4{|0⟩}_5\\ & -{|0⟩}_0{|\xi ⟩}_1{|0⟩}_2{|0⟩}_3{|0⟩}_4{|0⟩}_5\\ & -{|0⟩}_0{|0⟩}_1{|\xi ⟩}_2{|0⟩}_3{|0⟩}_4{|0⟩}_5\\ & +{|0⟩}_0{|0⟩}_1{|0⟩}_2{|\xi ⟩}_3{|0⟩}_4{|0⟩}_5\\ & -{|0⟩}_0{|0⟩}_1{|0⟩}_2{|0⟩}_3{|\xi ⟩}_4{|0⟩}_5\\ & +{|0⟩}_0{|0⟩}_1{|0⟩}_2{|0⟩}_3{|0⟩}_4{|\xi ⟩}_5),\end{array}$

where ${\tilde{\mathbb{N}}}^{\mathrm{\prime }}=(d+1)(1-\left|⟨0\right|{{\xi }^{\mathrm{\prime }}⟩|}^2)$ is the normalization coefficient. Different from the symmetric MECSVS, the normalization factor of the antisymmetric MECSVS ${\tilde{\mathbb{N}}}^{\mathrm{\prime }}\le d+1$ and therefore, the average photon number of the antisymmetric MECSVS can be larger than that of the single-mode CSVS. When ${\xi }^{\mathrm{\prime }}$ is the squeezed vacuum, we can see that the average photon number of the antisymmetric superposition of ESVS can be larger than that of the single-mode SVS, i.e., ${\sinh }^{2}r$. Based on Eq. (29), we can calculate the QCRB of multiparameter estimations with and without photon losses.

In order to more see whether this entangled state [Eq. (29)] can be also used to improve the precision of multiphase estimation or not, we plot the QCRB as a function of the mean photon number for given $T=0.9$ and different catalysis photon number $n=1,2,3$ in the ideal case [see Fig.5(a)] and in the case with photon losses with $\eta =0.9$ [see Fig.5(b)]. For comparison, the QCRBs of ESVS as input state with and without photon losses are also shown in the figures. From Fig.5(a), we can see that, for the entangled state given in Eq. (29), single-photon catalysis and two-photon catalysis have obvious advantages in improving the precision of multiphase estimation. However, when the number of catalytic photons $n\ge 3$, the QCRB of the catalytic entangled states does not surpass that of the ESVS when the mean photon number is large (i.e, larger than 4.13 in this example). From Fig.5(b), we can see that the QCRB in the case of single-photon catalysis is still lower than that in the ESVS case with photon losses. Hence, using the MECSVS in the antisymmetric case can also improve the phase estimation precision over using the ESVS.

6 Success probability

In this section, we discuss the success probability of our scheme. As known to us, the photon catalysis operation is probabilistic. The success probability depends on the transmissivity of the beam splitter and the number of photon catalysis which is given by Eq. (4). The higher the transmissivity, the higher the success probability. Fortunately, in our scheme the photon-catalysis is only performed once in the step one shown in Fig.1 and we find that the improvement of precision is the best when the transmissivity is around 0.9. For this transmissivity, the success probability of the photon catalysis is very high. For example, the success probability is about 70%−90% for the single-photon catalysis, 60%−80% for the two-photon catalysis, and 50%−70% for the three-photon catalysis when the squeezing parameter is in the range between 0 to 1 (Fig.9).

Fig.9 The success probability of the photon catalysis process as a function of squeezed parameter $r$ for three different number of photon catalysis. The transmissivity $T=0.9$.

Full size|PPT slide

On the other hand, the quantum-optical Fredkin gate used for the step 2 shown in Fig.1 is also probabilistic and different measurement results yield different output states. For example, when the photon detector $D_1$ detects a photon, the output state is the symmetric state $|\psi {⟩}_{ab}=(|0{⟩}_a|{\xi }^{\prime }{⟩}_b+|{\xi }^{\prime }{⟩}_a|0{⟩}_b)/\sqrt{N}$, while it is the antisymmetric state $|\psi {⟩}_{ab}=(-|0{⟩}_a|{\xi }^{\prime }{⟩}_b+|{\xi }^{\prime }{⟩}_a|0{⟩}_b)/\sqrt{N^{\prime }}$ when a single photon is detected in $D_2$. Despite the different outcomes, both of them can be applied to enhance the sensitivity of the phase measurement. For example, the improvement of phase estimation precision using the symmetric state is shown in Sections 3 and 4 and the improvement of phase estimation precision using the antisymmetric state is discussed in Section 5.

7 Conclusion

In this paper, we first propose a scheme to prepare the ESVS and MECSVS. Then we show that using MECSVS as input of a multi-arm optical interferometer, we can improve the QCRB of the multi-phase measurement over that using the corresponding usual ESVS in both cases with and without photon losses. We also compared the performance of the single-, two-, and three-photon catalysis states. The results shown that in both cases with and without photon losses for the symmetric MECSVS, the three-photon catalysis entangled state have better performance than the other two, which indicates that multi-photon catalysis has more resilient to the environment noises in this case. Additionally, for the antisymmetric MECSVS, it turns out that the single-photon catalyzed entangled state gives the best QCRB, which implies that the MECSVS with antisymmetric case can be also used to improve the phase estimation precision. Our results can find applications in the quantum metrology for multiparameter estimation.

8 Appendix A: Derivation of the two-mode ECSVS [Eq. (9)]

Assume that the initial input state is the single-mode multi-photon catalysis squeezed vacuum state ${|{\xi }^{\mathrm{\prime }}⟩}_a$ [see Eq. (3)] in the $a$ mode and the vacuum state ${|0⟩}_b$ in the $b$ mode. The output state of the step 2 shown in Fig.1 is given by

(A1)$\begin{array}{}& \begin{array}{rl}|\psi ⟩& ={\hat{B}}_2{\hat{B}}_4{\hat{U}}_k{\hat{B}}_3{\hat{B}}_1{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v\\ & =\frac{1}{\sqrt{2}}{\hat{B}}_2{\hat{U}}_F{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b\left({|1⟩}_u{|0⟩}_v+\mathrm{i}{|0⟩}_u{|1⟩}_v\right)\\ & =\frac{1}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\sum _{l=0}^{\mathrm{\infty }}{c}_l\left(\tau \right)\\ & \cdot {\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_0}{\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_2}{|2l⟩}_a{|0⟩}_b\left({|1⟩}_u{|0⟩}_v+\mathrm{i}{|0⟩}_u{|1⟩}_v\right){|}_{\tau =0}.\end{array}\end{array}$

Since $|n⟩=a^{†n}/\sqrt{n!}|0⟩$, Eq. (A1) can be rewritten

(A2)$\begin{array}{}& \begin{array}{rl}|\psi ⟩=& \frac{1}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}{\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_0}\\ & \cdot {\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_2}{a}^{†2l}{|0⟩}_a{|0⟩}_b\left({|1⟩}_u{|0⟩}_v+\mathrm{i}{|0⟩}_u{|1⟩}_v\right){|}_{\tau =0}.\end{array}\end{array}$

By using the formulas ${\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_i}{|1⟩}_u=$ ${\mathrm{e}}^{\mathrm{i}\chi J_i}{|1⟩}_u,$ ${\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_i}{|0⟩}_u={|0⟩}_u$ with $i=0,2$, ${\mathrm{e}}^{\mathrm{i}\lambda J_2}{a}^{†}{\mathrm{e}}^{-\mathrm{i}\lambda J_2}=a^{†}\cos \frac{\lambda }{2}-$ $b^{†}\sin \frac{\lambda }{2}$ and ${\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_2}{|0⟩}_a{|0⟩}_b={|0⟩}_a{|0⟩}_b$, we can obtain

(A3)$\begin{array}{}& \begin{array}{rl}|\psi ⟩=& \frac{1}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}\\ & \cdot {\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_0}{\left(a^{†}\cos \frac{\chi }{2}-b^{†}\sin \frac{\chi }{2}\right)}^{2l}\\ & \cdot {|0⟩}_a{|0⟩}_b\left({|1⟩}_u{|0⟩}_v+\mathrm{i}{|0⟩}_u{|1⟩}_v\right){|}_{\tau =0}\\ & +\frac{\mathrm{i}}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}{a}^{†2l}\\ & \cdot {|0⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v{|}_{\tau =0}.\end{array}\end{array}$

Taking $\chi =\pi$ yields

(A4)$\begin{array}{}& \begin{array}{rl}|\psi ⟩=& \frac{1}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\times \sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}{\mathrm{e}}^{\mathrm{i}\chi u^{†}u{J}_0}\\ & \cdot {\left(-b^{†}\right)}^{2l}{|0⟩}_a{|0⟩}_b\left({|1⟩}_u{|0⟩}_v+\mathrm{i}{|0⟩}_u{|1⟩}_v\right){|}_{\tau =0}\\ & +\frac{\mathrm{i}}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}{a}^{†2l}\\ & \cdot {|0⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v{|}_{\tau =0}.\end{array}\end{array}$

Finally, substituting $J_0=\frac{1}{2}\left(a^{†}a+b^{†}b\right)$ into Eq. (A4) and using ${\mathrm{e}}^{\lambda b^{†}b}{b}^{†}{\mathrm{e}}^{-\lambda b^{†}b}=b^{†}{e}^{\lambda }$, we can get

(A5)$\begin{array}{}& \begin{array}{rl}|\psi ⟩=& \frac{1}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\\ & \cdot \sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}{\mathrm{e}}^{\mathrm{i}\pi l}{b}^{†2l}{|0⟩}_a{|0⟩}_b{|1⟩}_u{|0⟩}_v{|}_{\tau =0}\\ & +\frac{\mathrm{i}}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\\ & \cdot \sum _{l=0}^{\mathrm{\infty }}\frac{c_l\left(\tau \right)}{\sqrt{\left(2l\right)!}}{a}^{†2l}{|0⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v{|}_{\tau =0}\\ =& \frac{1}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\\ & \cdot \sum _{l=0}^{\mathrm{\infty }}{c}_l\left(\tau \right){\mathrm{e}}^{\mathrm{i}\pi l}{|0⟩}_a{|2l⟩}_b{|1⟩}_u{|0⟩}_v{|}_{\tau =0}\\ & +\frac{\mathrm{i}}{\sqrt{2\mathbb{N}}}{\hat{B}}_2\frac{T^{n/2}}{n!}\frac{{\mathrm{\partial }}^n}{\mathrm{\partial }{\tau }^n}\beta \left(\tau \right)\\ & \cdot \sum _{l=0}^{\mathrm{\infty }}{c}_l\left(\tau \right){|2l⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v{|}_{\tau =0}.\end{array}\end{array}$

In the path of mode b, we place a quarter waveplate to compensate the $\pi$$l$ phase shift and the output state is

(A6)$\begin{array}{}& \begin{array}{r}|\psi ⟩=\frac{1}{\sqrt{2}}{\hat{B}}_2\left({|0⟩}_a{|{\xi }^{\mathrm{\prime }}⟩}_b{|1⟩}_u{|0⟩}_v+\mathrm{i}{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b{|0⟩}_u{|1⟩}_v\right).\end{array}\end{array}$

Then, after the transformations of the second beam splitter B2 of MZI-1, Eq. (6) becomes

(A7)$\begin{array}{}& \begin{array}{rl}|\psi ⟩=& \frac{1}{2}[\left({|0⟩}_a{|{\xi }^{\mathrm{\prime }}⟩}_b+{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b\right){|1⟩}_u{|0⟩}_v\\ & +\mathrm{i}\left(-{|0⟩}_a{|{\xi }^{\mathrm{\prime }}⟩}_b+{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b\right){|0⟩}_u{|1⟩}_v].\end{array}\end{array}$

From Eq. (A7), we can see that if a single photon is detected in the $u$ mode or in the $v$ mode, we can obtain

(A8)$\begin{array}{}& \begin{array}{r}{|\psi ⟩}_{out}=\frac{1}{\sqrt{{\tilde{N}}_2}}\left(\pm {|0⟩}_a{|{\xi }^{\mathrm{\prime }}⟩}_b+{|{\xi }^{\mathrm{\prime }}⟩}_a{|0⟩}_b\right)\end{array}\end{array}$

which is the two-mode ECSVS [Eq. (9)].

References

[1] C.W. Helstrom, Quantum Detection and Estimation Theory, Academic Press, New York, 1976

[2] S. L. Braunstein, C. M. Caves. Statistical distance and the geometry of quantum states. Phys. Rev. Lett., 1994, 72(22): 3439

[3] J. Joo, W. J. Munro, T. P. Spiller. Quantum metrology with entangled coherent states. Phys. Rev. Lett., 2011, 107(8): 083601

[4] J. Liu, X. X. Jing, X. G. Wang. Phase-matching condition for enhancement of phase sensitivity in quantum metrology. Phys. Rev. A, 2013, 88(4): 042316

[5] L. J. Fiderer, J. M. E. Fraïsse, D. Braun. Maximal quantum Fisher information for mixed states. Phys. Rev. Lett., 2019, 123(25): 250502

[6] W. Zhong, L. Zhou, Y. B. Sheng. Double-port measurements for robust quantum optical metrology. Phys. Rev. A, 2021, 103(4): 042611

[7] J. D. Zhang, C. L. You, C. Li, S. Wang. Phase sensitivity approaching the quantum Cramer−Rao bound in a modified SU(1, 1) interferometer. Phys. Rev. A, 2021, 103(3): 032617

[8] M. Eaton, R. Nehra, A. Win, O. Pfister. Heisenberg limited quantum interferometry with multiphotonsubtracted twin beams. Phys. Rev. A, 2021, 103(1): 013726

[9] R. Okamoto, T. Tahara. Precision limit for simultaneous phase and transmittance estimation with phase shifting interferometry. Phys. Rev. A, 2021, 104(3): 033521

[10] S. Y. Lee, Y. S. Ihn, Z. Kim. Optimal entangled coherent states in lossy quantum-enhanced metrology. Phys. Rev. A, 2020, 101(1): 012332

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

[12] L. Pezzé, A. Smerzi. Mach−Zehnder interferometry at the Heisenberg limit with coherent and squeezed-vacuum light. Phys. Rev. Lett., 2008, 100(7): 073601

[13] C. Oh, S. Y. Lee, H. Nha, H. Jeong. Practical resources and measurements for lossy optical quantum metrology. Phys. Rev. A, 2017, 96(6): 062304

[14] S. Slussarenko, M. M. Weston, H. M. Chrzanowski, L. K. Shalm, V. B. Verma, S. W. Nam, G. J. Pryde. Unconditional violation of the shot-noise limit in photonic quantum metrology. Nat. Photonics, 2017, 11(11): 700

[15] L. L. Guo, Y. F. Yu, Z. M. Zhang. Improving the phase sensitivity of an SU(1, 1) interferometer with photon-added squeezed vacuum light. Opt. Express, 2018, 26(22): 29099

[16] R. Birrittella, C. C. Gerry. Quantum optical interferometry via the mixing of coherent and photon-subtracted squeezed vacuum states of light. J. Opt. Soc. Am. B, 2014, 31(3): 586

[17] S. Wang, X. X. Xu, Y. J. Xu, L. J. Zhang. Quantum interferometry via a coherent state mixed with a photonadded squeezed vacuum state. Opt. Commun., 2019, 444: 102

[18] D. Braun, P. Jian, O. Pinel, N. Treps. Precision measurements with photon-subtracted or photon-added Gaussian states. Phys. Rev. A, 2014, 90(1): 013821

[19] H. Zhang, W. Ye, C. P. Wei, Y. Xia, S. K. Chang, Z. Y. Liao, L. Y. Hu. Improved phase sensitivity in a quantum optical interferometer based on multi-photon catalytic two-mode squeezed vacuum states. Phys. Rev. A, 2021, 103(1): 013705

[20] Y. Ouyang, S. Wang, L. J. Zhang. Quantum optical interferometry via the photon added two-mode squeezed vacuum states. J. Opt. Soc. Am. B, 2016, 33(7): 1373

[21] M. A. Taylor, W. P. Bowen. Quantum metrology and its application in biology. Phys. Rep., 2016, 615: 1

[22] N. P. Mauranyapin, L. S. Madsen, M. A. Taylor, M. Waleed, W. P. Bowen. Evanescent single-molecule biosensing with quantum-limited precision. Nat. Photonics, 2017, 11(8): 477

[23] M. A. Taylor, J. Janousek, V. Daria, J. Knittel, B. Hage, H. A. Bachor, W. P. Bowen. Biological measurement beyond the quantum limit. Nat. Photonics, 2013, 7(3): 229

[24] M. Tsang, R. Nair, X. M. Lu. Quantum theory of superresolution for two incoherent optical point sources. Phys. Rev. X, 2016, 6(3): 031033

[25] L. J. Fiderer, T. Tufarelli, S. Piano, G. Adesso. General expressions for the quantum Fisher information matrix with applications to discrete quantum imaging. PRX Quantum, 2021, 2(2): 020308

[26] C. Lupo, Z. X. Huang, P. Kok. Quantum limits to incoherent imaging are achieved by linear interferometry. Phys. Rev. Lett., 2020, 124(8): 080503

[27] G. Brida, M. Genovese, I. Ruo Berchera. Experimental realization of sub-shot-noise quantum imaging. Nat. Photonics, 2010, 4(4): 227

[28] L. Pezzé. Entanglement-enhanced sensor networks. Nat. Photonics, 2021, 15(2): 74

[29] P. C. Humphreys, M. Barbieri, A. Datta, I. A. Walmsley. Quantum enhanced multiple phase estimation. Phys. Rev. Lett., 2013, 111(7): 070403

[30] J. D. Yue, Y. R. Zhang, H. Fan. Quantum-enhanced metrology for multiple phase estimation with noise. Sci. Rep., 2014, 4(1): 5933

[31] L. B. Ho, H. Hakoshima, Y. Matsuzaki, M. Matsuzaki, Y. Kondo. Multiparameter quantum estimation under dephasing noise. Phys. Rev. A, 2020, 102(2): 022602

[32] Y. Yao, L. Ge, X. Xiao, X. G. Wang, C. P. Sun. Multiple phase estimation for arbitrary pure states under white noise. Phys. Rev. A, 2014, 90(6): 062113

[33] S. L. Hong, J. Rehman, Y. S. Kim, Y. W. Cho, S. W. Lee, H. Jung, S. Moon, S. W. Han, H. T. Lim. Quantum enhanced multiple-phase estimation with multi-mode NOON states. Nat. Commun., 2021, 12(1): 5211

[34] J. Liu, X. M. Lu, Z. Sun, X. G. Wang. Quantum multiparameter metrology with generalized entangled coherent state. J. Phys. A, 2016, 49(11): 115302

[35] Z. B. Hou, R. J. Wang, J. F. Tang, H. D. Yuan, G. Y. Xiang, C. F. Li, G. C. Guo. Control-enhanced sequential scheme for general quantum parameter estimation at the Heisenberg limit. Phys. Rev. Lett., 2019, 123(4): 040501

[36] H. Kwon, Y. Lim, L. Jiang, H. Jeong, C. Oh. Quantum metrological power of continuous-variable quantum networks. Phys. Rev. Lett., 2022, 128(18): 180503

[37] S. S. Pang, A. N. Jordan. Optimal adaptive control for quantum metrology with time-dependent Hamiltonians. Nat. Commun., 2017, 8(1): 14695

[38] J. Yang, S. S. Pang, Y. Y. Zhou, A. N. Jordan. Optimal measurements for quantum multiparameter estimation with general states. Phys. Rev. A, 2019, 100(3): 032104

[39] R. Nichols, P. Liuzzo-Scorpo, P. A. Knott, G. Adesso. Multiparameter Gaussian quantum metrology. Phys. Rev. A, 2018, 98(1): 012114

[40] L. Pezzé, M. A. Ciampini, N. Spagnolo, P. C. Humphreys, A. Datta, I. A. Walmsley, M. Barbieri, F. Sciarrino, A. Smerzi. Optimal measurements for simultaneous quantum estimation of multiple phases. Phys. Rev. Lett., 2017, 119(13): 130504

[41] W. Ge, K. Jacobs, Z. Eldredge, A. V. Gorshkov, M. Foss-Feig. Distributed quantum metrology with linear networks and separable inputs. Phys. Rev. Lett., 2018, 121(4): 043604

[42] Z. B. Hou, J. F. Tang, H. Z. Chen, H. D. Yuan, G. Y. Xiang, C. F. Li, G. C. Guo. Zero-trade-off multiparameter quantum estimation via simultaneously saturating multiple Heisenberg uncertainty relations. Sci. Adv., 2021, 7(1): eabd2986

[43] L. Zhang, K. W. C. Chan. Quantum multiparameter estimation with generalized balanced multimode NOON like states. Phys. Rev. A, 2017, 95(3): 032321

[44] X. X. Jing, J. Liu, W. Zhong, X. G. Wang. Quantum Fisher information of entangled coherent states in a lossy Mach–Zehnder interferometer. Commum. Theor. Phys., 2014, 61(1): 115

[45] J. Joo, K. Park, H. Jeong, W. J. Munro, K. Nemoto, T. P. Spiller. Quantum metrology for nonlinear phase shifts with entangled coherent states. Phys. Rev. A, 2012, 86(4): 043828

[46] H. Vahlbruch, M. Mehmet, K. Danzmann, R. Schnabel. Detection of 15 dB squeezed states of light and their application for the absolute calibration of photonelectric quantum efficiency. Phys. Rev. Lett., 2016, 117(11): 110801

[47] Y. J. Wang, W. H. Zhang, R. X. Li, L. Tian, Y. H. Zheng. Generation of −10.7 dB unbiased entangled states of light. Appl. Phys. Lett., 2021, 118(13): 134001

[48] T. J. Bartley, G. Donati, J. B. Spring, X. M. Jin, M. Barbieri, A. Datta, B. J. Smith, I. A. Walmsley. Multi-photon state engineering by heralded interference between single photons and coherent states. Phys. Rev. A, 2012, 86(4): 043820

[49] A. I. Lvovsky, J. Mlynek. Quantum-optical catalysis: Generating nonclassical states of light by means of linear optics. Phys. Rev. Lett., 2002, 88(25): 250401

[50] L. Y. Hu, Z. Y. Liao, M. S. Zubairy. Continuousvariable entanglement via multi-photon catalysis. Phys. Rev. A, 2017, 95(1): 012310

[51] X. X. Xu. Enhancing quantum entanglement and quantum teleportation for two-mode squeezed vacuum state by local quantum-optical catalysis. Phys. Rev. A, 2015, 92(1): 012318

[52] Y. Guo, W. Ye, H. Zhong, Q. Liao. Continuousvariable quantum key distribution with non-Gaussian quantum catalysis. Phys. Rev. A, 2019, 99(3): 032327

[53] W. Ye, H. Zhong, Q. Liao, D. Huang, L. Y. Hu, Y. Guo. Improvement of self-referenced continuous-variable quantum key distribution with quantum photon catalysis. Opt. Express, 2019, 27(12): 17186

[54] A. E. Ulanov, I. A. Fedorov, A. A. Pushkina, Y. V. Kurochkin, T. C. Ralph, A. I. Lvovsky. Undoing the effect of loss on quantum entanglement. Nat. Photonics, 2015, 9: 764

[55] F. Jia, W. Ye, Q. Wang, L. Y. Hu, H. Y. Fan. Comparison of nonclassical properties resulting from non-Gaussian operations. Laser Phys. Lett., 2019, 16(1): 015201

[56] C. Kumar, Rishabh Arora. Realistic non-Gaussian-operation scheme in parity-detection-based Mach−Zehnder quantum interferometry. Phys. Rev. A, 2022, 105(5): 052437

[57] C. C. Gerry, R. A. Campos. Generation of maximally entangled photonic states with a quantum-optical Fredkin gate. Phys. Rev. A, 2001, 64(6): 063814

[58] N. Imoto, H. A. Haus, Y. Yamamoto. Quantum nondemolition measurement of the photon number via the optical Kerr effect. Phys. Rev. A, 1985, 32(4): 2287

[59] H. Schmidt, A. Imamoglu. Giant Kerr nonlinearities obtained by electromagnetically induced transparency. Opt. Lett., 1996, 21(23): 1936

[60] D. Vitali, M. Fortunato, P. Tombesi. Complete quantum teleportation with a Kerr nonlinearity. Phys. Rev. Lett., 2000, 85(2): 445

[61] A. B. Matsko, I. Novikova, G. R. Welch, M. S. Zubairy. Enhancement of Kerr nonlinearity by multiphoton coherence. Opt. Lett., 2003, 28(2): 96

[62] H. Kang, Y. F. Zhu. Observation of large Kerr nonlinearity at low light intensities. Phys. Rev. Lett., 2003, 91(9): 093601

[63] S. Rebić, J. Twamley, G. J. Milburn. Giant Kerr nonlinearities in circuit quantum electrodynamics. Phys. Rev. Lett., 2009, 103(15): 150503

[64] H. Azuma. Quantum computation with Kerr-nonlinear photonic crystals. J. Phys. D, 2008, 41(2): 025102

[65] L. S. Costanzo, A. S. Coelho, N. Biagi, J. Fiurášek, M. Bellini, A. Zavatta. Measurement-induced strong Kerr nonlinearity for weak quantum states of light. Phys. Rev. Lett., 2017, 119(1): 013601

[66] A. Fujiwara. Estimation of SU(2) operation and dense coding: An information geometric approach. Phys. Rev. A, 2001, 65(1): 012316

[67] K. Matsumoto. A new approach to the Cramer–Rao type bound of the pure-state model. J. Phys. Math. Gen., 2002, 35(13): 3111

[68] M. G. A. Paris. Quantum Estimation for quantum technology. Int. J. Quant. Inf., 2009, 7(supp01): 125

[69] B. M. Escher, R. L. de Matos Filho, L. Davidovich. General framework for estimating the ultimate precision limit in noisy quantum-enhanced metrology. Nat. Phys., 2011, 7(5): 406

[70] M. Sarovar, G. J. Milburn. Optimal estimation of one parameter quantum channels. J. Phys. Math. Gen., 2006, 39(26): 8487

[71] H. Zhang, W. Ye, C. P. Wei, C. J. Liu, Z. Y. Liao, L. Y. Hu. Improving phase estimation using number conserving operations. Phys. Rev. A, 2021, 103(5): 052602

[72] J. Liu, H. D. Yuan, X. M. Lu, X. G. Wang. Quantum Fisher information matrix and multiparameter estimation. J. Phys. A Math. Theor., 2020, 53(2): 023001

[73] T. Baumgratz, A. Datta. Quantum enhanced estimation of a multidimensional field. Phys. Rev. Lett., 2016, 116(3): 030801

Acknowledgements

This work was supported by the National Key R&D Program of China (Grant No. 2021YFA1400800), the Key-Area Research and Development Program of Guangdong Province (Grant No. 2018B030329001), the Natural Science Foundation of Guangdong (Grant No. 2021A1515010039), the National Natural Science Foundation of China (No. 11964013), and the Major Discipline Academic and Technical Leaders Training Program of Jiangxi Province (No. 20204BCJL22053).

RIGHTS & PERMISSIONS

2023 Higher Education Press