# Multipartite continuous-variable entanglement distillation using local squeezing and only one photon-subtraction operation

## Abstract

In this paper, we study entanglement distillation of multipartite continuous-variable Gaussian entangled states. Following Opatrný et al.’s photon subtraction (PS) scheme, the probability of successful distillation decreases exponentially with the number of parties . However, here, we shall propose an entanglement distillation scheme whose success probability scales as a constant with . Our protocol employs several local squeezers, but it requires only a single PS operation. Using the logarithmic negativity as a measure of entanglement, we find that both the success probability and the distilled entanglement can be improved at the same time. Moreover, an -mode transfer theorem (transferring states from phase space to Hilbert space) is presented.

###### pacs:

03.67.Mn, 03.67.Hk, 42.50.DvEntanglement, particularly multipartite entanglement state, is one of the most fundamental and puzzling aspect in quantum mechanics. However, entanglement is such a fragile resource that it may be easily degraded during its interaction with the environmental noise. To this point, entanglement distillation (always in a probabilistic way) has been proposed to increase the entanglement in the noise-disturbed entangled stateHage08 (); Dong (). Restricted by the famous No-Go theorem in continuous variable (CV) entanglement distillationnogo1 (); nogo2 (); nogo3 (), lots of efforts have been devoted to the non-Gaussian operations. As an example, photon subtraction(PS) operationOpatrny (), proposed by Opatrný et al in 2000, is principally simple and can be readily implemented with beamsplitter and photon detectors. Very recently, about 10 years after Opatrný et al’s pioneering work, an experiment which faithfully implements the PS-based two mode entanglement distillation has been reported NatPhotonic (). One of the challenge in this experiment is the extremely low successful probability, which is mainly due to the rather-high-transmittance beamsplitter used in PS operation—For one thing, the beamsplitter must own a relatively high transmittance to guarantee an entanglement-enhanced distillationdistill2copy (). For the other, high transmittance means low reflectiveness in beamsplitter and hence, low probability in PS operation. Assuming the beamsplitter’s transmittance is , as shown in Ref.NatPhotonic (), detectors ’s efficiency , the probability of each successful local PS operation is upper-bounded by . The successful probability of the whole entanglement distillation will be even lower, decreasing exponentially with the number of local PS operations: . This is really a serious problem if -partite continuous variable entanglement distillation is involved.

In this paper, we consider the distillation of -partite CV gaussian entanglement state with only one-time photon subtraction. Generally for partite (especially, the symmetric) Gaussian state, one-time PS will modify the permutation symmetry and the distilled entanglement, measured with logarithmic negativity (log-neg), will be even worse. However, if assisted with local squeezing, we show that the entanglement can be improvedcv2density (). Moreover, one-time PS can give a substantial increment in distilling probability, which keeps constants for arbitrary partite number .

Our entanglement distilling scheme is briefly shown in Fig.1, where Fig.1(a) describes the typical -time PS scheme which is a straightforward generalization of Opatrný et al’s scheme to partite case. Fig.1(b) is the our new distilling protocol with local squeezing ( described by symplectic operation ) and only one-time PS. Throughout this paper, for convenience, we assume that the second input-mode of beam-splitter are vacuum modes, representing a simple and concise expression of PS operation in Phase space(See AppendixA).

Preliminaries. Our result can be conveniently derived in Phase space. Let’s now introduce the basic facts and notation of CV modes state. First, it is convenient to express each mode, say mode, with the field quadrature operators , with being the mode annihilation and creation operators. By defining the vector of quadrature operators , the commutation relation can be written as , with The density matrix of arbitrary -mode system which resides in infinite-dimensional Hilbert space can be conveniently represented by the characteristic function in -dimensional real vector space, i.e., phase space . Gaussian state is a special kind of quantum state whose characteristic function is Gaussian in phase space, , inwhich is quadrature average and denotes the covariance matrix . In the following, we also use the wigner function, the Fourier transform of to express the state evolution. A standard normalized wigner function is defined by with and .

Multi-partite entanglement states. Multipartite CV Gaussian entanglement plays a prominent role in future quantum network and quantum communication protocoljzhang (); adesso (). In this paper, we are mainly interested in the distilling of a family of genuinely partite symmetric Gaussian entangled stateloockInequa (); loockTHEO (); aokiEXP (). The covariance matrix is given by

(1) |

inwhich is the unnormalized dimensional real vector , and is the identity matrix in dimensions. are diagonal matrix

(2) |

The genuinely partite state in Eq.(1) can be experimentally prepared with a particular sequence of phase-free beam splitters and squeezed input statesloockInequa (). To be simple, we will mainly focused on the unbiased states, namely .

Distillation with local squeezing and one-time PS. Let’s now derive the state evolution of the our one-time PS distilling protocol. As shown in Fig.1(b), local squeezing symplectic transformations are applied before PS. This corresponds in phase space to a transformation of covariance matrix . Without loss of generality, we assume that only PS operation is performed in the mode. Implemented with a beam-splitter (transmittance ), the PS operation (Fig.1(b) inset), couples the mode Gaussain state with the vacuum mode . The mode states now follows , with being the symplectic matrix Finally, a successful distillation is heralded if the detector register non-vacuum results. According to AppendixA, one can find the distilled state is a linear combination of two gaussian state

(3) |

with and being a normalized -paritie gaussian state with covariance matrix . The and are defined by partitioning of matrix , where are matrice respectively. The success probability of distilling follows . With the transfer theorem form phase space to Hilbert space (See Appendix B), one can easily compare the entanglement before and after distillation.

In Fig.2, we evaluate the entanglement after and before distillation with the log-neg logneg () as the figure of merit for entanglement. For simplicity, we consider three-mode entanglement distillation as the example. For a fixed initial squeezing , as shown in Fig.2(a), we plot the entanglement as a function of . Through our simulation, we assume the transmittance of PS beamsplitter is which is available in recent experimentsNatPhotonic (). The probability of success is briefly shown in Fig.2(b). It should be noted that the success probability which is about . This is mainly due to the rather low initial squeezing (), which results extremely low photon number in each transmission mode, which certainly decrease the probability of being photon subtracted. Our method can be applied for even stronger squeezing. In Fig.(c), we increase and find the optimal squeezing which may maximize the log-neg of output entanglement state. The numerical results support the linearity reliance of upon the increasing . Also, we plot the corresponding optimized log-neg and success probability in Fig.2(d) and Fig.3. The success probability (Fig.3) is about which is an pronounced improvement compared with the -time PS strategy().

Discussions and Prospectives. We presented here a photon-subtraction based entanglement distillation for arbitrary partite continuous variable entanglement state. As an example, now in this paper, only the three-partite symmetric Gaussian state is involved. This method is applicable for arbitrary partite CV state. Indeed, even for , this improvement in both log-neg and success probability also applies. As an auxiliary result, we also derive the transfer theorem for partite Gaussian state from Phase space to Hilbert space. We can envisage that this theorem could find more application in the entanglement evaluation tasks, such as entanglement swapping and entanglement distribution.

Acknowledgements. The authors acknowledge financial support from National Fundamental Research Program, also by National Natural Science Foundation of China (Grant No. 10674128 and 60121503) and the Innovation Funds and “Hundreds of Talents” program of Chinese Academy of Sciences and Doctor Foundation of Education Ministry of China (Grant No. 20060358043). SZ acknowledges support by the Max Planck Gesellschaft, Chinese Academy of Sciences Joint Doctoral Promotion Programme (MPG-CAS-DPP). PvL acknowledges support from the Emmy Noether Program of the DFG.

## Appendix A Wigner function description of Photon Subtraction

In ideal cases, a perfect photon subtraction is described with the annihilation operation: in Hilbert space. However, this is not an unitary operation and cannot be implemented deterministically. A convenient way is to use beamsplitter and photon detectorsOpatrny (). For ease, we consider a single partite state as input, the PS operation (see Fig.1) can be represented with a completely-positive map from to normalized output state , with

(4) |

where denotes the positive operators projecting to non-vacuum subspace and denotes the Beam-splitting operation between and modes.

In our calculation, indeed, it is convenient to use the wigner function to describe the PS process above (Eq.(4)). In fact, the operator is a difference of two operation whose wigner function are both GaussianPRLCVCHSH (), i.e.,

(5) |

In case that the input state is Gaussian, the wigner-function of distilled entanglement can be easily formulated with the linear combination of a series of Gaussian function, each of which can be conveniently expressed with the covariance marices.

## Appendix B Quantum State from phase-space to Hilbert Space

In this section, we give the detailed techniques we use in the processing of multi-partite Gaussian quantum state from phase-space to Hilbert Space. For a -partite CV state, the density matrix follows

with being the mode displacement operator.

The matrix entries of can be conveniently obtained by observing the equationcv2density ():

where is the -dimensional real vector.

By noticing the fact that

(6) | |||

one obtains that (after integrating the )

(7) |

inwhich

(8) | |||||

and is the Pauli matrices. Then, one can check that the state (Eq.(7)) is now automatically normalized, i.e., .

### References

- B. Hage, Aiko Samblowski, James Diguglielmo, Alexander Franzen, Jaromír Fiurášek and Roman Schnabel, Nat. Phys. 4, 915 (2008).
- R. Dong, Mikael Lassen, Joel Heersink, Christoph Marquardt, Radim Filip, Gerd Leuchs and Ulrik L Andersen, Nat. Phys. 4, 919 (2008).
- J. Eisert, S. Scheel, M. B. Plenio, Phys. Rev. Lett, 89, 137903 (2002).
- J. Fiurášek, Phys. Rev. Lett. 89, 137904 (2002).
- G. Giedke, J. I. Cirac, Phys. Rev. A 66, 032316 (2002).
- T. Opatrný, G. Kurizki and D. G. Welsch, Phys. Rev. A 61, 032302 (2000).
- H. Takahashi, J. S. Neergaard-Nielsen, M. Takeuchi, M. Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, Nature Photonic 4, 178 (2010).
- ShengLi Zhang, Peter van Loock, Phys. Rev. A 82, 062316 (2010).
- ShengLi Zhang, Peter van Loock, arXiv:1103.4500v1 .
- Jing Zhang, Gerardo Adesso, Changde Xie and Kunchi Peng, Phys. Rev. Lett. 103, 070501(2009).
- Gerardo Adesso, and Fabrizio Illuminati, Phys. Rev. A. 78, 042310(2008).
- Peter van Loock and Akira Furusawa, Phys. Rev. A 67, 052315(2003).
- P. van Loock, and Samuel L. Braunstein, Phys. Rev. Lett. 84, 3482 (2000).
- Takao Aoki, Nobuyuki Takei, Hidehiro Yonezawa, Kentaro Wakui, Takuji Hiraoka, Akira Furusawa and Peter van Loock, Phys. Rev. Lett. 91, 080404 (2003).
- G. Vidal and R.F. Werner, Phys. Rev. A 65, 032314 (2002).
- R. García-Patrón, J. Fiurášek, N. J. Cerf, J. Wenger, R. Tualle-Brouri, and Ph. Grangier, Phys. Rev. Lett. 93, 130409 (2004).