Abstract
Two-mode squeezed states are scalable and robust entanglement resources for continuous-variable and hybrid quantum information protocols that are realized at a distance. We consider the effect of a linear cross talk in the multimode distribution of two-mode squeezed states propagating through parallel similar channels. First, to reduce degradation of the distributed Gaussian entanglement, we show that the initial two-mode squeezing entering the channel should be optimized already in the presence of a small cross talk. Second, we suggest simultaneous optimization of relative phase between the modes and their linear coupling on a receiver side prior to the use of entanglement, which can fully compensate the cross talk once the channel transmittance is the same for all the modes. For the realistic channels with similar transmittance values for either of the modes, the cross talk can be still largely compensated. This method relying on the mode interference overcomes an alternative method of entanglement localization in one pair of modes using measurement on another pair and feed-forward control. Our theoretical results pave the way to more efficient use of multimode continuous-variable photonic entanglement in scalable quantum networks with cross talk.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Photonic quantum entanglement [1] is not only a puzzling physical phenomenon, but as well a resource for quantum information processing and quantum communication tasks, in particular for quantum key distribution (QKD) [2,3], quantum metrology [4] or quantum computing [5]. Continuous-variable (CV) entanglement, using generally multiphoton states of light, has experimentally enabled quantum communication and information processing with large information capacity [6–14]. It can be quantified using logarithmic negativity (LN) [15], being a measure of negativity of a state’s partial transpose [16,17]. Advantageously, CV entanglement can be deterministically generated reaching two-mode squeezing bellow $-10 dB$ (corresponds to logarithmic negativity up to $LN=4.3$) [18]. In CV quantum communication, one of the possible ways to simultaneously transfer many entangled states during the same period of time is to use frequency or spatial multiplexing of quantum states residing in different modes [19–22]. Such an approach was in particular used to improve quantum communication with multiplexed QKD [23] or with multiplexed quantum teleportation using frequency pulse modes [24]. The techniques for preparation and detection of multimode quantum states of light were drastically improved in the past years, which enabled generation of highly multimode frequency comb states [25,26] with the focus on quantum networks [27,28] or cluster states, which are multiplexed in the time domain [29–32].
In a massive mode multiplexing, it is challenging to avoid cross talk between the modes during preparation and distribution [33–35], which can reduce or even destroy logarithmic negativity of entangled pairs and undermine applicability of shared entanglement. It was, in particular, shown that cross talk effects in the multimode homodyne detection can undermine security of mode-non-discriminating CV QKD [36–38]. However, the cross talk already in the state preparation before the distribution may substantially influence also other entanglement-based protocols, especially, if they are implemented over attenuating (lossy) channels [26,27].
In the current paper, we theoretically study the role of linear cross talk between the signal modes on the state preparation side in distribution of multimode CV entangled states of light through an imperfect (lossy and noisy) channel (differently from the recent study of cross talk from the co-propagating classical signals in multimode CV QKD [39]). We consider Gaussian bipartite two-mode squeezed vacuum states (TMSV) [40], also known as twin beams, and show that even a small cross talk substantially degrades Gaussian entanglement contained in the mode pairs and makes the states more sensitive to excess noise in quantum channels used for entanglement distribution. Importantly, amount of initial two-mode squeezing should be optimally adjusted to maximize its transmission in the presence of cross talk. To overcome this basic passive method, we suggest active control of phase between the modes and controllable linear coupling of the output modes prior to their use on a remote side in order to compensate the cross talk, inspired by the method used to remove correlations in quantum memory channels [41] and in the squeezed state generation [42]. We show that such an active method, if used optimally, can completely eliminate the negative effect of cross talk once the noiseless channel transmittance is the same for all the modes. For a realistic case of asymmetrical quantum channels with slightly different transmittance for different modes, the cross talk can be still largely compensated for. To demonstrate the advantage we compare our method of active cross talk elimination to a deterministic entanglement localization scheme [43,44] using homodyne detection on one pair of modes and feed-forward control on another pair. We show that the active scheme provides better results in a wide range of realistic parameters, particularly, at relatively small cross talk and similar channel transmittance values for different modes, while remarkably preserving multimode structure of the signal. We also address stability of our method and show that it is robust against deviations of the linear coupling used to compensate the cross talk.
2. Effect of cross talk on multimode continuous-variable entanglement distribution
We consider multimode entangled idler and multimode entangled signal beams constituting TMSV emitted by a source [40]. Such states are entangled so that quadratures $\hat {x}=\hat {a}^{\dagger}+\hat {a}$ and $\hat {p}=i(\hat {a}^{\dagger}-\hat {a})$, defined for a given mode, are strongly respectively correlated and anti-correlated within a pair of modes belonging to the signal and the idler beams. To verify the entanglement, a sender (Alice) is performing homodyne quadrature measurement on each of the modes of a respective beam (e.g., signal) and a receiver (Bob) is measuring quadratures of each of many modes of another twin beam (e.g, idler) using another homodyne detector, after the modes experience cross talk and travel through generally attenuating and noisy quantum channels. In the relevant experimental examples [26,27,45] the crosstalk appears already in the source prior to the lossy and noisy channel.
We assume the most common linear cross talk between two neighbouring modes, which can be represented by a beamsplitter coupling $t_{c}$ between the modes, as it is schematically shown in Fig. 1. Changing $t_c$ from 0 to 1 then means a transition from very strong inter-mode coupling (hence strong cross talk) to the absence of cross talk. Here we analyse weak nearest-neighbour coupling, although in general case multimode coupling is more complex, the model allows us to obtain necessary conditions for improvement to further numerically apply these methods to a more complex multiple mode cross talk. We parametrize the quantum channel with transmittance $T_i$ for a given $i$-th mode and with the amount of phase-insensitive excess noise $\varepsilon$ added to the quadrature variance, which is assumed (and typically is) the same for all the modes. Excess noise in general depends on properties of the channel , but can be also concerned with unaccounted detection noise or imperfect state generation. For the fiber channels the observed excess noise is typically below 1 % shot noise unit (SNU), being the level of vacuum quadrature fluctuations, at the channel output [46,47]. We first analyse the role of cross talk in a basic case of two entangled mode pairs (so that the overall number of modes in signal and idler beams is four), then extend the results to three pairs.
We describe the Gaussian multimode TMSV states [40] by covariance matrices with elements $\gamma _{ij}=\langle r_i r_j \rangle$, where $r_i=\{x_i,p_j\}$ is a set of quadratures of a given $i$-th mode, taking into account their zero mean values, $i \in [1,2N]$, here $2N$ is the overall number of modes (hence, $N$ is the number of mode pairs, or, equivalently, the number of modes in each signal and idler beams). Then the system of four pairwise perfectly entangled modes $A_{\{1,2\}},B_{\{1,2\}}$, shown in Fig. 1, each having quadrature variance $V \geq 1$, is described by the 8x8 covariance matrix of the form
Here $\mathbb {I} = { diag[1,1]}$ is the unity matrix, $\mathbb {Z}= {diag[1,-1]}$ is a Pauli matrix. Quadrature variance $V$ is related to the squeezing parameter $r$ of TMSV, which defines the amount of two-mode squeezing applied to the two-mode vacuum in both quadratures $\textrm {Var}[(x_{A_i}+x_{B_i})/\sqrt {2}]=e^{-2r}$ and $\textrm {Var}[(p_{A_i}-p_{B_i})/\sqrt {2}]=e^{-2r}$, and $V=\cosh (2r)$, so that the larger $r$, the stronger is quadrature correlations between the modes, given by $\langle (x_{A_i}+x_{B_i}) (p_{A_i}-p_{B_i}) \rangle =\sinh {r}$. After a linear cross talk $t_c$ and lossy and noisy quantum channel, the covariance matrix changes to
Here $r_c \equiv 1-t_c$. We characterize the bipartite Gaussian entanglement using LN [15] of a state $\rho$ is defined as
where $\rho ^\Gamma$ is a partial transpose of $\rho$, $||\rho ||_1$ is the trace norm of the operator $\rho$, that is equal to the sum of the absolute values of the negative eigenvalues of $\rho ^\Gamma$, quantifying the degree to which a partially transposed state fails to be positive [16,17]. For a Gaussian states with covariance matrix $\gamma$ the Eq. (3) becomes the sum of all symplectic eigenvalues of the partially transposed that are less than one: $LN(\gamma ) = \sum _k \max \lbrace 0, - \log _2 \nu _k\rbrace$. For the TMSV states with the symplectic eigenvalues $\{\nu _+, \nu _-\}$ the larger eigenvalue is $\nu _+\geqslant 1$ [48]. This way the LN of the $i$-th pair of modes simplifies to where $\nu _-$ is the smallest symplectic eigenvalue of the covariance matrix of the partially transposed state of either the first or the second pair. We limit our study to the Gaussian TMSV states as the typical case of bipartite quadrature-entangled states of light. As the task is to deterministically transmit all pairs through the channel, we evaluate logarithmic negativity for each pair of modes separately. We also define the initial logarithmic negativity of a TMSV state (1), before the cross talk and before sharing the state over a noisy and lossy channel, see Fig. 1, as which is the same for the both pairs of modes. It represents the maximum of Gaussian entanglement in the pairs of modes, which can be shared in the perfect case of no cross talk and ideal (lossless and noiseless) channels. In the limit of small state variance ($V\to 1$), it behaves as $LN_0\sim \frac {1}{\log 2} \sqrt {2(V-1)}$. In the limit of large $V \to \infty$, $LN_0 \sim \log _2 (2V)$.In the absence of cross talk ($t_c = 1$) the logarithmic negativity monotonously grows with increasing TMSV quadrature variance $V$, but its limit depends on the channel parameters as
For noiseless channels ($\varepsilon =0$) the entaglement never vanishes for any $T_i>0$. For lossless channels ($T_i=1$) entanglement is lost at $\varepsilon > 2$. In the presence of cross talk between the signal modes and after a lossy and noisy channel, as shown in Fig. 1, the logarithmic negativity becomes
For channels with high loss ($T_i\ll 1$), expanding (7) around $T_i=0$ gives
For a very small transmittance, $T_i\to 0$, LN is well described by the first term in (8), $LN \sim \frac {T_i [2-\varepsilon -(1-t_c) (V+1)] }{\log (2)}$. In this limit, the shared entanglement has linear dependence on the transmittance $T_i$. Provided that the cross talk is absent ($t_c=1$), the entanglement is independent on the initial variance $V$ and is destroyed by the excess noise $\varepsilon =2$. The presence of cross talk introduces a new term $-T_i(1-t_c) V$ that reduces the shared entanglement, which then monotonically decreases with the increase of $V$.
Beyond the limit of small initial entanglement ($V \rightarrow 1$), sensitivity of Gaussian entanglement of TMSV to cross talk can be seen in Fig. 2, left, where logarithmic negativity is plotted versus the initial entanglement at different cross talk coupling $t_c$ values and compared to the case when the cross talk is absent ($t_c=1$). It is numerically evident from (7) and from the plots in Fig. 2, left, that the amount of entanglement shared over the noiseless channel in the presence of cross talk becomes sensitive to the initial entanglement even for weak channel attenuation $T_i=0.9$ and can be broken once the initial entanglement is too strong. This is concerned with the fact that due to cross talk the signal from an adjacent mode is coupled to the signal in another mode and thus effectively acts an uncorrelated noise appearing in one of the beams if all modes are treated independently. Amount of such a noise is larger for the higher state variance, i.e., for a larger initial entanglement. This destructive extra noise effect from the crosstalk is also dominant for channels with excess noise. Moreover, the presence of cross talk degrades robustness of the shared entanglement to the channel noise, as can be seen from Fig. 2, right. The shared entanglement is then destroyed by the level of noise, which is independent of the channel transmittance $T_i$ (hence being valid for any $T_i>0$), and is given by
Equivalently, the bound on the shared entanglement (variance $V$ for which entanglement turns to $0$) is the same for any $T_i>0$, and it reads
Importantly, the maximum tolerable initial variance given by $V_{max}$ (or, equivalently, maximum tolerable initial entanglement) also does not depend on the channel transmittance and, for a purely attenuating channel, depends only on the cross talk coupling $t_c$. The necessary condition $V<V_{max}$ therefore simplifies optimization if the cross talk is present only in the source. Note that if the channel is present also before the crosstalk, optimization becomes more complex and involves the channel transmittance values, however, the limitation $V<V_{max}$ remains.
The above given results were obtained in the assumption that the cross talk occurs between two modes. In the case, when each of the multiple TMSV modes interacts with two neighbouring modes, i.e., the cross talk dominantly occurs between three modes, the Eqs. (7) and (10) hold up to the substitution $t_c \to t_c^2$, which leads to stronger sensitivity of the shared entanglement to the linear cross talk between the signal modes. In this case, $V_{max}$ is even more limited and therefore less entanglement can be transmitted.
The shared entanglement, in the presence of cross talk, can be maximized by optimizing the initial entangled resource. The optimal state variance $V$, which maximizes the shared entanglement, reads
By optimizing the initial TMSV variance as in (11) we both reduce the negative influence of the cross talk, and rise tolerance to the excess noise. The maximal tolerable level of excess noise then becomes $\varepsilon _{max}=2 t_c$ SNU.
For a small cross talk $t_c \to 1$ the optimal variance can be approximated as
If we neglect the excess noise, the expression (12) simplifies to
Putting $T_i=1$ in (13) allows to easily demonstrate the existence of the upper bound on the optimal state variance $V_{opt}(T=1)=\frac {1}{\sqrt {1-t_c}}$. It illustrates the fact that for any cross talk, even without noise and attenuation, the maximal possible shared entanglement is rather limited. For a strongly attenuating channel ($T_i \to 0$) together with excess noise $\varepsilon$, the optimal variance $V_{opt}$ is approximately
For a channel with strong attenuation ($T_i \to 0$), the optimized entanglement can be then approximated as
In the limit of a very small cross talk ($t_c >0.98$) for relatively strong excess noise ($\varepsilon >0.2$ SNU) it is possible to find better approximation expanding optimized entanglement $LN(V_{opt})$ into series for a very low cross talk ($t_c \sim 1$) in the the limit of strongly attenuating channel ($T_i \to 0$), then the optimized entanglement can be better analytically approximated as
Figures 2 and 3 illustrate that although optimization of the state variance during the state preparation helps to somewhat mitigate the negative effects of cross talk (allowing to reach maximal logarithmic negativity for given cross talk and excess noise that corresponds to the maximums in Fig. 2), nevertheless such optimization fails to fully compensate for the cross cross talk especially in the presence of a stronger excess noise. Comparison of the left and the right panel in the Fig. 3 and the Eqs. (15)–(17) show that in the strongly attenuating channel while for a weak cross talk logarithmic negativity decreases as $LN\sim -\sqrt {1-t_c}$ (Eq. (15) and dashed lines), very soon it gets suppressed more by higher orders of $t_c$ (Eq. (15)–(17) and dotted lines). Thus, even though initial entanglement can be optimized to achieve the maximum shared entanglement in the presence of cross talk, the shared entanglement will be still largely reduced by the cross talk, especially if combined with the channel loss and/or noise.
We therefore consider the local active manipulations on either sender or receiver side in order to improve and possibly fully restore the shared entanglement by compensating the cross talk. If the mode coupling happens only in the source during the state preparation, it can be compensated by manipulating either the modes $A_1$ and $A_2$, or $B_1$ and $B_2$ on the sender side or the modes $B_1$ and $B_2$ on the remote side after the states are shared through an imperfect channel. If an additional mode coupling also happens in the channel, manipulations on the receiver side are needed in addition to the manipulations at the sender side in order to simultaneously compensate both cross talk in the state preparation and in the channel. Generally both, the cross talk and the losses, can occur on both sender and receiver side and also in the source itself during the state preparation and distribution. In our study we assume that the shared part of the states (hence, modes $B_1$ and $B_2$) is the subject to cross talk. We then consider the manipulations on the receiver side, i.e., after the imperfect channel, in order to compensate the cross talk. The case of an ideal cross talk compensation on the sender side would be equivalent to the absence of cross talk, as can be seen in the next section.
3. Cross talk compensation by optical interference
We consider the feasible state operations prior to using the entangled states in order to compensate the cross talk, i.e., to reduce or completely eliminate its negative effect on the entanglement distribution. Since the cross talk is typically caused by an energetically passive photon exchange between the signal modes, we consider the beam-splitter type interaction based on a variable coupling $t_r$ between the signal modes. The compensating interaction can in principle be performed on either sender or receiver side. If no signal loss occurs on the sender side, it is sufficient to implement the beam-splitter with $t_r=t_c$ coupling the modes $A_1$ and $A_2$ at the sender. To take into account the effects caused by the lossy channels we further concentrate on the case where the compensation is implemented on the modes $B_1$ and $B_2$ by the remote party of the entanglement distribution scheme. In this case the beam splitter needs to be preceded by an optimal phase shift (in the case of linear cross talk between two modes it is given by $\pi$) on one of the modes (e.g., mode 2), as shown in Fig. 1. Then for a pair $A_1,B_1$ the covariance matrix reads
Choosing $t_r=\frac {T_2 t_c}{T_2 t_c+ T_1 r_c}$ the states change to
We have therefore eliminated the contribution of the second mode $B_2$ to the first one $B_1$, but not vice versa. If, on the other hand, we chose $t_r=\frac {T_2 t_c}{T_2 t_c+T_1 r_c}$, the cross talk will be eliminated from the second mode, but not from the first one. It is easy to see, that if $T_1=T_2 \equiv T$, (22) and (23) turn into $|\alpha ''\rangle _{B_1} = \left | \sqrt {T}\alpha \right \rangle _{B_1}$ and $|\beta ''\rangle _{B_2}= \left |\sqrt {T }\beta \right \rangle _{B_2}$ as it should be for two independent modes attenuated by a channel. Therefore for the perfectly balanced channels ($T_1=T_2$) the cross talk is fully eliminated for a pair of arbitrary pure two-mode state.
When the channel transmittance values $T_1,T_2$ for modes $B_1,B_2$ are different, the cross talk cannot be fully compensated in both modes and entanglement cannot be fully restored in both transmitting channels. While transmittance values for different modes are typically similar in fiber channels [49], they can vary, e.g. for frequency modes in atmospheric channels [50]. We therefore consider the efficiency of the cross talk compensation in the unbalanced channels. In order to maximally reconstruct entanglement, one has to optimize the coupling $t_r$, as shown in Fig. 4. It is evident from the plots, that entanglement can still be largely compensated for weaker crosstalk ($t_c\to 1$) between asymmetrical lossy channels and that the method of optical interference is stable to the setting $t_r$. There is a difference between the optimal $t_r$ for the first pair and for the second one (it corresponds to the maximums of the curves on the left and right plots in Fig. 4), although, as it noticeable from the plots, their maximums are very close. In the second pair $A_2B_2$ we even can overcome the entanglement before the cross talk because we get stronger signal in the pair $A_2B_2$ that is taken from the pair $A_1B_1$ due to the cross talk.
In a general case, $t_r$ should be optimized numerically to result in the maximally restored entanglement. However, the optimal setting can be derived analytically in the relevant limits of small and large initial entanglement (or, equivalently, small and high large initial state variance $V$). For a purely lossy channel ($\varepsilon =0$), the optimal transmittance, that maximizes the entanglement in the first pair of modes $A_1B_1$, reads
andAnd vice versa ($T_1\leftrightarrow T_2$) for another pair. It is important to note that the optimal $t_r$ for any value of the state variance $V$ always lies between the bounds given by (24–25) and, for realistically close $T_1$ and $T_2$, this interval is quite narrow. For $t_r$ defined by (25) the logarithmic negativity is an increasing function of $V$ and for large initial state variance it approaches the limit:
Note, that the above given Eqs. (24)–(26) are applicable for maximization of the enanglement for the pair $A_1,B_1$. If, on the other hand, the goal is to maximize the logarithmic negativity in the second pair $A_2,B_2$ , the same Eqs. (24)–(26) can be applied with replacement $T_1 \leftrightarrow T_2$. This also implies, that in the case of unbalanced channels $T_1 \neq T_2$, the entanglement cannot be optimally restored in both modes. While applying the optical interference method to compensate the cross talk in the channels with different transmittance for each mode ($T_1 \neq T_2$), one has to chose one mode in which to maximize the entanglement at the expense of the rest of the modes. The Eqs. (24)–(26) refer to the first pair of modes $A_1B_1$. However, since in practical situations strong unbalancing between transmittance values for different modes is unlikely, i.e. $T_1$ is typically close to $T_2$, optimal settings for $t_r$ are also close and choosing $t_r$ in between the optimal settings for either of the pairs will give nearly optimal results in terms of entanglement for both the pairs.
Advantage of the decoupling by interference is that cross talk from all entangled pairs can be nearly perfectly removed. Moreover, technical imperfections, as the modes match imperfectly at the decoupling BS, can be incorporated to the full optimization. Despite such a principal advantage, one may be interested in maximizing entanglement in only one pair of modes. In this case entanglement localization [43,44] can be applied by performing only an optimized measurement on one pair and feed-forward control on another pair of modes based on the measurement outcomes, which we consider in the next Section and compare it to the above suggested method of reversed coupling.
4. Entanglement localization by measurement and feed-forward control
In order to compensate the negative effect of cross talk in only one entangled pair in the scheme given in Fig. 1, we now consider a possibility to apply optimized Gaussian measurement of one pair of modes and feed-forward control of the other pair on the receiving sides and thus to increase the shared entanglement at the cost of losing one of the mode pairs. This method relies on the highly efficient and low-noise homodyne measurement and the high fidelity feed-forward control and coherent displacement [9,51]. Similarly as in the previous Section, we theoretically consider a simple two-mode case with the feed-forward control applied by Alice and Bob on the first pair of modes $A_1,B_1$ after the measurement on the second pair of modes $A_2,B_2$, as shown in Fig. 5. Instead of a direct interference of multimode signals on a coupler with variable $t_r$, we interfere the signal modes $A_2, B_2$ with the local oscillator beams in a generalized Gaussian measurement [52] on both Alice and Bob sides of the scheme, as shown in Fig. 5. Photocurrents from this detection control the modulation units in the other pair of modes $A_1, B_1$. We therefore assume that both the sender and the receiver perform a general Gaussian measurement [52] on the second pair of modes ($A_2,B_2$), i.e., a generalized heterodyne (also known as double-homodyne) measurement, using unbalanced beam-splitters with transmittance values $t_A$ and $t_B$ in Alice’s and Bob’s sides respectively. The beam-splitters divide each of the measured modes $A_2,B_2$ (into auxiliary detection modes denoted as $C_A,C_B,D_A,D_B$) and the outputs are then measured in conjugate quadratures using two homodyne detectors (with no loss of generality we assume $x$-quadrature measured on $C_A,C_B$ and $p$-quadrature measured on $D_A,D_B$). It is advantageous to optimize the measurement shown in Fig. 5, we are therefore searching for optimal $t_A$ and $t_B$ that maximize the logarithmic negativity of the conditional state in the first pair of modes $A_1,B_1$. We evaluate conditional matrix of the state in modes $A_1,B_1$ after a homodyne measurement of a quadrature $r$ (being either $x$ or $p$) in mode $K$ (being one of $C_A,C_B,D_A,D_B$). The covariance matrix of the modes $A_1, B_1$ and $K$ before the measurement is $\gamma _{A_1B_1 K}= \left ( \begin {array}{ccc} {\gamma }_{A_1} & c_{A_1 B_1} & c_{A_1 K}\\ c_{A_1 B_1} & {\gamma }_{B_1} & c_{B_1 K}\\ c_{A_1 K} & c_{B_1 K} & {\gamma }_K \end {array} \right )$. The measurement resulting in an outcome $r_K$ transforms the covariance matrix as [53]
First, we observe that the best possible result for any given channel transmittance values and any cross talk level is achieved by using homodyne measurement on the receiver side. It corresponds to $t_B=1$ or $t_B=0$, i.e. measuring either only $C_B$ (homodyne measurement of $x$-quadrature) or only $D_B$ (homodyne measurement of $p$-quadrature) mode in Fig. 5. We further assume with no loss of generality that $t_B=1$, which is equivalent to Bob measuring $x$-quadrature of the mode $A_2$ with a homodyne detector. On the other hand, an optimal $t_A$ on the receiver side generally depends on the TMSV state variance $V$, channel parameters (transmittance values and excess noise) and the cross talk coupling $t_c$. In the general case $t_A$ can be found only numerically, but in the limit of large $V$ and noiseless channels, we can find the $t_A$ that maximizes the logarithmic negativity in $\gamma _{A_1B_1|x_{C_A},x_{B_2},p_{D_A}}$ analytically. Then for $V\to \infty$ and $\varepsilon =0$ in the limit of the weak cross talk $t_c\lesssim 1$ the optimal $t_A$ on the sender side is independent of the TMSV variance and can be approximated as
In the case, when additionally losses are low ($T_1 \to 1, T_2 \to 1$), the optimal setting at the sender side is $t_A=0$. This means that for the highly transmitting channels or low cross talk the best strategy for the feed-forward method of cross talk compensation is to use homodyne measurement on Alice’s side and measure complementary quadrature to the one that the receiver (Bob) measures. This phenomena can be explained using continuous-variable quantum erasing, where Alice by measurement prepares a squeezed state in front of the cross talk which, using homodyne measurement of antisqueezed variable and feed-forward control, can be used to eliminate beam-splitter type of cross talk [54,55]. The covariance matrix for a noiseless channel then turns to
In a simple case of a perfect channel ($T_i=1$ and $\varepsilon =0$), the logarithmic negativity of the state in modes $A_1,B_1$ (conditioned on homodyne measurement outcomes on both sides) can be expressed analytically as
It is always larger than the logarithmic negativity in the same pair of modes after the crosstalk and in the same perfect channels, but without conditioning, which reads $LN = - \log _2 \left [V-\sqrt {t_c \left (V^2-1\right )}\right ]$. The latter equation reaches maximum when $V_{opt}=\frac {1}{\sqrt {1-t_c}}$, while (33) is a constantly growing function of $V$. It illustrates the fact that the conditional measurement overcomes the limit on the initial entanglement given by (10) and generally any need to optimize the initial state. Therefore it allows to use as high initial entanglement as it is experimentally achievable. When the state is conditioned on the optimal homodyne measurement in a highly transmitting channel, the logarithmic negativity of the state in modes $A_1,B_1$ is not bounded and in the limit of arbitrarily strong TMSV variance $V \to \infty$ , assuming a noiseless channel ($\varepsilon =0$) tends to
In the case, when the channels are strongly attenuating ($T_{1,2}\to 0$), the optimal $t_A$ given by (31) turns to $1/2$, i.e., the balanced heterodyne detection on the sender side becomes the optimal measurement. In this case, when we also chose the optimal $t_B=1$, the conditional covariance matrix of the state in modes $A_1,B_1$ becomes (assuming no excess noise, $\varepsilon =0$)
The logarithmic negativity is then as well not bounded and the optimization of the initial entanglement (11) is not needed. In the limit of very strong initial entanglement $V \rightarrow \infty$ (assuming no excess noise, $\varepsilon =0$) it tends to
In the general case of intermediate channel attenuation and level of cross talk an optimized heterodyne measurement on the sender side should be chosen by evaluating the respective logarithmic negativities, but entanglement in the remaining pair of modes after the measurement and feed-forward control is always no less, than before these operations. Therefore, optimized measurement and feed-forward can always improve (and fully restore in the case of almost perfectly balanced channels) entanglement in a single pair of modes at the cost of loosing another pair. We illustrate the efficiency and stability of the method with respect to the sender measurement settings in Fig. 6 and show that it can largely restore the entanglement in the remaining pair of modes and is stable with respect to the measurement setting $t_A$. As it is evident from Fig. 6, any kind of Gaussian measurement on the pair $A_2,B_2$, even not an optimized one, largely restores entanglement in the pair $A_1,B_1$. The best result is achieved with the homodyne measurement of the $x$-quadrature on the sender side ($t_A=1$) in the low loss channels and with the heterodyne measurement ($t_A=1/2$) in the strong loss channels, which agrees with approximate analytical results given above in (31). Remarkably, the optimized measurement allows to recover entanglement that was completely destroyed by the cross talk.
We also compare the method of optimized measurement and feed-forward to the method of optical interference, suggested in the previous section, as shown in Fig. 7 for low and high channel losses.
It is evident from the plots, that the method of optimized optical interference is always more efficient in restoring the Gaussian entanglement and, besides, it preserves the multimode structure. Our further analysis shows, that this superiority holds always, unless the channels are strongly unbalanced (with either of the channels being strongly attenuating while another one being almost lossless, $T_i \to 1$ and $T_j \to 0$), which is, however, very unlikely in practical situations. While being less efficient, the measurement strategy can be sometimes nearly optimal and become very close to the optimized measurement strategy, however, the measurement inevitably destroys one pair modes. Both of the methods increase the amount of entanglement, that can be transferred, and remove the necessity to optimize the initial entanglement. The lines corresponding to interference method (purple) and to the conditional measurement (green and yellow) continue to the right approaching the asymptotes given by (26), (34), (36), and the lines corresponding to the ideal case without cross talk (blue) approach the asymptote (6). We also compare the methods in Fig. 8 for the same initial entanglement in the large region of transmittance values and in the presence of unbalancing. In Fig. 8 we plot the Gaussian entanglement for the both pairs of modes in case of the optimized optical interference and show that method allows to restore the entanglement in both the modes contrary to the method of conditional measurement. For a weaker unbalancing the entanglement of both pairs of modes after the optimized interference would overlap. Only one (blue) line in Fig. 8 corresponds to the conditional measurement method, it demonstrates that although this method is quite effective for restoring entanglement in one pair, it also destroys the second pair of modes. We have not presented the case with no cross talk on the plot because for the given set of parameters it will be almost fully overlapping with the results for the optimized optical interference method. It illustrates the fact that even for relatively strong cross talk, the optimized interference approach allows to recover the entanglement almost completely. Both active methods give substantial gain compared to initial state before recovery. For any cross talk $t_c$, state variance $V$, and channel transmittance $T_1$, $T_2$ both the proposed methods outperform the simple optimization of the initial entanglement shared (or equivalently the variance $V$) as proposed in the Sec. 2.
Conditional measurement with feed-forward control gives guaranteed entanglement gain no matter how strong is the cross talk and does not rely on our knowledge of cross talk coupling $t_c$, while discarding one of the modes. However, this benefit is heavily paid for by a low probability of success of this asymptotic method. On the other hand, the optical interference scheme allows to preserve all the modes. It gives better results than the entanglement concentration with the conditional measurement scheme for most realistic values of $t_c$, especially in channels with higher loss. In the case of the state having only several modes, both of the methods perform comparably. For multimode states with higher number of modes the advantage of the optimal interference method will be more pronounced as it preserves all the modes, while the conditional measurement method is based on discarding most of them.
5. Conclusion
Entanglement distribution is a key ingredient of modern quantum technology. We considered the effect of cross talk on entanglement distribution using Gaussian multimode twin-beam states and shown that initial entanglement inserted to a multimode link has to be optimized to reach its maximal transmission already if cross talk is weakly present. We then suggest the method of cross talk compensation based on optimized optical interference using adjustable phase control and linear coupling of modes, and show that the method can completely eliminate the negative effect of cross talk once the channel transmittance is balanced, i.e., is the same for all the modes. The method is still very efficient for small unbalancing between the channel transmittance values for different modes and is stable with respect to the linear coupling setting, but requires knowledge of the cross talk strength. This methodology can be extended to many multiplexed modes, however, in such a case the numerical analysis and optimization are needed for specific real channels, as the model of the cross talk will be more complicated and vary in practical situations. We therefore leave the analysis of the real cross talk in the source or channel for future experimental tests. As an alternative, we suggest the method of cross talk compensation for one of the pairs of modes by optimized Gaussian measurement on another pair and feed-forward control, which does not rely on knowledge of cross talk strength and can restore entanglement in the remaining pair, while reducing the multimode structure, and is also limited in the unbalanced channels. The method of optical interference remains more efficient unless the channels are very strongly unbalanced, which is however not likely in the practical situations, and preserves the amount of entangled modes transmitted through the channel. Our methods can be prospective for realization of multimode quantum communication in the presence of cross talk in the sources and channels. Note that the described cross talk compensation methods distribute entanglement between the modes, but do not increase it beyond the initial entanglement in the source before the cross talk, hence not contradicting the impossibility of entanglement distillation by local Gaussian unitary operations [53,57]. Our results demonstrate basic principles in a simple case of only nearest-mode cross talk. The further investigation for larger number of modes in the presence of cross talk will require multi-parameter numerical optimizations and possibly application of modern methods of deep machine learning applicable to quantum optics [58] to find efficient strategy of cross talk compensation. Such numerical optimizations in application to complex multimode cross-talk effects will be the subject of future studies.
Funding
Univerzita Palackého v Olomouci (IGA-PrF-2021-006); Ministerstvo Školství, Mládeže a Tělovýchovy (CZ 02.1.01/0.0/0.0/16_026/0008460); Horizon 2020 Framework Programme (731473 (project 8C20002 ShoQC), 820466 CIVIQ); Grantová Agentura České Republiky (19-23739S).
Disclosures
The authors declare that there are no conflicts of interest related to this article.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. W. T. Tittel and G. Weihs, “Photonic entanglement for fundamental tests and quantum communication,” QIC 1(2), 3–56 (2001). [CrossRef]
2. N. Gisin, G. Ribordy, W. Tittel, and H. Zbinden, “Quantum cryptography,” Rev. Mod. Phys. 74(1), 145–195 (2002). [CrossRef]
3. S. Pirandola, U. L. Andersen, L. Banchi, M. Berta, D. Bunandar, R. Colbeck, D. Englund, T. Gehring, C. Lupo, C. Ottaviani, J. L. Pereira, M. Razavi, J. S. Shaari, M. Tomamichel, V. C. Usenko, G. Vallone, P. Villoresi, and P. Wallden, “Advances in quantum cryptography,” Adv. Opt. Photonics 12(4), 1012–1236 (2020). [CrossRef]
4. P. Kómár, E. M. Kessler, M. Bishof, L. Jiang, A. S. Sørensen, J. Ye, and M. D. Lukin, “A quantum network of clocks,” Nat. Phys. 10(8), 582–587 (2014). [CrossRef]
5. H. J. Kimble, “The quantum internet,” Nature 453(7198), 1023–1030 (2008). [CrossRef]
6. R. Ukai, S. Yokoyama, J.-i. Yoshikawa, P. van Loock, and A. Furusawa, “Demonstration of a controlled-phase gate for continuous-variable one-way quantum computation,” Phys. Rev. Lett. 107(25), 250501 (2011). [CrossRef]
7. X. Su, H. Shuhong, X. Deng, L. Ma, M. Wang, X. Jia, C. Xie, and K. Peng, “Gate sequence for continuous variable one-way quantum computation,” Nat. Commun. 4(1), 2828 (2013). [CrossRef]
8. R. Filip, “Continuous-variable quantum nondemolishing interaction at a distance,” Phys. Rev. A 69(5), 052313 (2004). [CrossRef]
9. S. Yokoyama, R. Ukai, J.-i. Yoshikawa, P. Marek, R. Filip, and A. Furusawa, “Nonlocal quantum gate on quantum continuous variables with minimum resources,” Phys. Rev. A 90(1), 012311 (2014). [CrossRef]
10. O. Morin, K. Huang, J. Liu, H. Le Jeannic, C. Fabre, and J. Laurat, “Remote creation of hybrid entanglement between particle-like and wave-like optical qubits,” Nat. Photonics 8(7), 570–574 (2014). [CrossRef]
11. H. Jeong, A. Zavatta, M. Kang, S.-W. Lee, L. S. Costanzo, S. Grandi, T. C. Ralph, and M. Bellini, “Generation of hybrid entanglement of light,” Nat. Photonics 8(7), 564–569 (2014). [CrossRef]
12. Z. Yan, L. Wu, X. Jia, Y. Liu, R. Deng, S. Li, H. Wang, C. Xie, and K. Peng, “Establishing and storing of deterministic quantum entanglement among three distant atomic ensembles,” Nat. Commun. 8(1), 718 (2017). [CrossRef]
13. M. Huo, J. Qin, J. Cheng, Z. Yan, Z. Qin, X. Su, X. Jia, C. Xie, and K. Peng, “Deterministic quantum teleportation through fiber channels,” Sci. Adv. 4(10), eaas9401 (2018). [CrossRef]
14. K. Huang, H. Le Jeannic, O. Morin, T. Darras, G. Guccione, A. Cavailles, and J. Laurat, “Engineering optical hybrid entanglement between discrete- and continuous-variable states,” New J. Phys. 21(8), 083033 (2019). [CrossRef]
15. G. Vidal and R. F. Werner, “Computable measure of entanglement,” Phys. Rev. A 65(3), 032314 (2002). [CrossRef]
16. A. Peres, “Separability criterion for density matrices,” Phys. Rev. Lett. 77(8), 1413–1415 (1996). [CrossRef]
17. M. Horodecki, P. Horodecki, and R. Horodecki, “Separability of mixed states: necessary and sufficient conditions,” Phys. Lett. A 223(1-2), 1–8 (1996). [CrossRef]
18. T. Eberle, V. Händchen, and R. Schnabel, “Stable control of 10 db two-mode squeezed vacuum states of light,” Opt. Express 21(9), 11546–11553 (2013). [CrossRef]
19. B. Hage, A. Samblowski, and R. Schnabel, “Towards Einstein-Podolsky-Rosen quantum channel multiplexing,” Phys. Rev. A 81(6), 062301 (2010). [CrossRef]
20. M. Heurs, J. G. Webb, A. E. Dunlop, C. C. Harb, T. C. Ralph, and E. H. Huntington, “Multiplexed communication over a high-speed quantum channel,” Phys. Rev. A 81(3), 032325 (2010). [CrossRef]
21. T. Kouadou, L. L. Volpe, S. De, C. Fabre, V. Parigi, and N. Treps, “Single-pass generation of spatial and spectral multimode squeezed states of light,” in Frontiers in Optics + Laser Science APS/DLS, (Optical Society of America, 2019), p. FTh3B.7.
22. S. Shi, L. Tian, Y. Wang, Y. Zheng, C. Xie, and K. Peng, “Demonstration of channel multiplexing quantum communication exploiting entangled sideband modes,” Phys. Rev. Lett. 125(7), 070502 (2020). [CrossRef]
23. Z. Qu and I. B. Djordjevic, “Four-dimensionally multiplexed eight-state continuous-variable quantum key distribution over turbulent channels,” IEEE Photonics J. 9(6), 1–8 (2017). [CrossRef]
24. A. Christ, C. Lupo, and C. Silberhorn, “Exponentially enhanced quantum communication rate by multiplexing continuous-variable teleportation,” New J. Phys. 14(8), 083007 (2012). [CrossRef]
25. G. de Valcarcel, G. Patera, N. Treps, and C. Fabre, “Multimode squeezing of frequency combs,” Phys. Rev. A 74(6), 061801 (2006). [CrossRef]
26. O. Pinel, P. Jian, R. M. De Araujo, J. Feng, B. Chalopin, C. Fabre, and N. Treps, “Generation and characterization of multimode quantum frequency combs,” Phys. Rev. Lett. 108(8), 083601 (2012). [CrossRef]
27. J. Roslund, R. M. De Araujo, S. Jiang, C. Fabre, and N. Treps, “Wavelength-multiplexed quantum networks with ultrafast frequency combs,” Nat. Photonics 8(2), 109–112 (2014). [CrossRef]
28. Y. Cai, Y. Xiang, Y. Liu, Q. He, and N. Treps, “Versatile multipartite Einstein-Podolsky-Rosen steering via a quantum frequency comb,” Phys. Rev. Res. 2(3), 032046(R) (2020). [CrossRef]
29. S. Yokoyama, R. Ukai, S. C. Armstrong, C. Sornphiphatphong, T. Kaji, S. Suzuki, J.-i. Yoshikawa, H. Yonezawa, N. C. Menicucci, and A. Furusawa, “Ultra-large-scale continuous-variable cluster states multiplexed in the time domain,” Nat. Photonics 7(12), 982–986 (2013). [CrossRef]
30. J.-i. Yoshikawa, S. Yokoyama, T. Kaji, C. Sornphiphatphong, Y. Shiozawa, K. Makino, and A. Furusawa, “Invited article: Generation of one-million-mode continuous-variable cluster state by unlimited time-domain multiplexing,” APL Photonics 1(6), 060801 (2016). [CrossRef]
31. F. Arzani, C. Fabre, and N. Treps, “Versatile engineering of multimode squeezed states by optimizing the pump spectral profile in spontaneous parametric down-conversion,” Phys. Rev. A 97(3), 033808 (2018). [CrossRef]
32. X. Zhu, C.-H. Chang, C. González-Arciniegas, A. Pe’er, J. Higgins, and O. Pfister, “Hypercubic cluster states in the phase modulated quantum optical frequency comb,” arXiv preprint arXiv:1912.11215.
33. T. Kudo and T. Ishigure, “Analysis of interchannel crosstalk in multimode parallel optical waveguides using the beam propagation method,” Opt. Express 22(8), 9675–9686 (2014). [CrossRef]
34. L. Szostkiewicz, M. Napierala, A. Ziolowicz, A. Pytel, T. Tenderenda, and T. Nasilowski, “Cross talk analysis in multicore optical fibers by supermode theory,” Opt. Lett. 41(16), 3759–3762 (2016). [CrossRef]
35. B. Nada and T. Berceli, “Crosstalk reduction in fiber links using double polarization,” Opt. Quantum Electron. 52(4), 200 (2020). [CrossRef]
36. V. C. Usenko, L. Ruppert, and R. Filip, “Entanglement-based continuous-variable quantum key distribution with multimode states and detectors,” Phys. Rev. A 90(6), 062326 (2014). [CrossRef]
37. V. C. Usenko, L. Ruppert, and R. Filip, “Quantum communication with macroscopically bright nonclassical states,” Opt. Express 23(24), 31534–31543 (2015). [CrossRef]
38. O. Kovalenko, K. Y. Spasibko, M. V. Chekhova, V. C. Usenko, and R. Filip, “Feasibility of quantum key distribution with macroscopically bright coherent light,” Opt. Express 27(25), 36154 (2019). [CrossRef]
39. T. A. Eriksson, B. J. Puttnam, G. Rademacher, R. S. Luís, M. Fujiwara, M. Takeoka, Y. Awaji, M. Sasaki, and N. Wada, “Crosstalk impact on continuous variable quantum key distribution in multicore fiber transmission,” IEEE Photonics Technol. Lett. 31(6), 467–470 (2019). [CrossRef]
40. C. Weedbrook, S. Pirandola, R. García-Patrón, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, “Gaussian quantum information,” Rev. Mod. Phys. 84(2), 621–669 (2012). [CrossRef]
41. C. Lupo, L. Memarzadeh, and S. Mancini, “Removing correlations in signals transmitted over a quantum memory channel,” Phys. Rev. A 85(1), 012320 (2012). [CrossRef]
42. R. Filip, “Squeezing restoration by a noisy probe from a classically correlated environment,” Phys. Rev. A 81(3), 032330 (2010). [CrossRef]
43. F. Verstraete, M. Popp, and J. I. Cirac, “Entanglement versus correlations in spin systems,” Phys. Rev. Lett. 92(2), 027901 (2004). [CrossRef]
44. F. Sciarrino, E. Nagali, F. De Martini, M. Gavenda, and R. Filip, “Entanglement localization after coupling to an incoherent noisy system,” Phys. Rev. A 79(6), 060304 (2009). [CrossRef]
45. M. Chen, N. Menicucci, and O. Pfister, “Experimental realization of multipartite entanglement of 60 modes of a quantum optical frequency comb,” Phys. Rev. Lett. 112(12), 120505 (2014). [CrossRef]
46. P. Jouguet, S. Kunz-Jacques, E. Diamanti, and A. Leverrier, “Analysis of imperfections in practical continuous-variable quantum key distribution,” Phys. Rev. A 86(3), 032309 (2012). [CrossRef]
47. D. Huang, P. Huang, D. Lin, and G. Zeng, “Long-distance continuous-variable quantum key distribution by controlling excess noise,” Sci. Rep. 6(1), 19201 (2016). [CrossRef]
48. G. Adesso, A. Serafini, and F. Illuminati, “Extremal entanglement and mixedness in continuous variable systems,” Phys. Rev. A 70(2), 022318 (2004). [CrossRef]
49. V. Alwayn, Optical Network Design and Implementation (Cisco Press, 2014).
50. J. L. Streete, “Infrared measurements of atmospheric transmission at sea level,” Appl. Opt. 7(8), 1545–1549 (1968). [CrossRef]
51. Y. Shiozawa, J.-i. Yoshikawa, S. Yokoyama, T. Kaji, K. Makino, T. Serikawa, R. Nakamura, S. Suzuki, S. Yamazaki, W. Asavanant, S. Takeda, P. Loock, and A. Furusawa, “Quantum nondemolition gate operations and measurements in real time on fluctuating signals,” Phys. Rev. A 98(5), 052311 (2018). [CrossRef]
52. K. Miyata, H. Ogawa, P. Marek, R. Filip, H. Yonezawa, J.-i. Yoshikawa, and A. Furusawa, “Implementation of a quantum cubic gate by an adaptive non-gaussian measurement,” Phys. Rev. A 93(2), 022301 (2016). [CrossRef]
53. J. Eisert, S. Scheel, and M. B. Plenio, “Distilling gaussian states with gaussian operations is impossible,” Phys. Rev. Lett. 89(13), 137903 (2002). [CrossRef]
54. R. Filip, “Continuous-variable quantum erasing,” Phys. Rev. A 67(4), 042111 (2003). [CrossRef]
55. U. Andersen, O. Glöckl, S. Lorenz, G. Leuchs, and R. Filip, “Experimental demonstration of continuous variable quantum erasing,” Phys. Rev. Lett. 93(10), 100403 (2004). [CrossRef]
56. M. Sabuncu, R. Filip, G. Leuchs, and U. Andersen, “Environmental assisted quantum information correction for continuous variables,” Phys. Rev. A 81(1), 012325 (2010). [CrossRef]
57. J. Fiurášek, “Gaussian transformations and distillation of entangled gaussian states,” Phys. Rev. Lett. 89(13), 137904 (2002). [CrossRef]
58. V. Dunjko and H. J. Briegel, “Machine learning and artificial intelligence in the quantum domain: a review of recent progress,” Rep. Prog. Phys. 81(7), 074001 (2018). [CrossRef]