# Simple proof of the robustness of Gaussian entanglement in bosonic noisy channels

###### Abstract

The extremality of Gaussian states is exploited to show that Gaussian states are the most robust, among all possible bipartite continuous-variable states at fixed energy, against disentanglement due to noisy evolutions in Markovian Gaussian channels involving dissipation and thermal hopping. This proves a conjecture raised recently in [M. Allegra, P. Giorda, and M. G. A. Paris, Phys. Rev. Lett. 105, 100503 (2010)], providing a rigorous validation of the conclusions of that work. The problem of identifying continuous variable states with maximum resilience to entanglement damping in more general bosonic open system dynamical evolutions, possibly including correlated noise and non-Markovian effects, remains open.

###### pacs:

03.67.Mn, 03.65.YzContinuous-variable (CV) systems such as light modes and ultracold atomic ensembles COVAQIAL () provide advantageous resources to achieve unconditional implementations of quantum information processing brareview (), ranging from teleportation protocols Furusawa98 () to quantum key distribution cvcrypto () and one-way quantum computation 1waycv (). Gaussian states and Gaussian operations, that represent respectively the most natural and easily controllable light states as well as the set of manipulations efficiently realizable by linear optics, have traditionally occupied a privileged role in all such implementations. Furthermore, by virtue of their mathematical simplicity compared to general states living in infinite-dimensional Hilbert spaces, bosonic Gaussian states have been and are the preferred testgrounds for a broad variety of investigations into the structure, nature, and dynamics of CV entanglement and quantum correlations ourreview ().

Some prominent limitations of a Gaussian-only toolbox have been however recently exposed in several contexts nogos (), stimulating vigorous theoretical and experimental research into the realm of non-Gaussian state engineering and characterization belliniecc (), to assess and harness the potentially enhanced performance of de-Gaussified CV resources for quantum teleportation degausstelep (), entanglement distillation distillscheme (), parameter estimation souza (), universal quantum computation cvcomput (); 1waycv (), non-locality tests loopholefree (), etc. Still, one of the most powerful features of Gaussian quantum states is their extremality extra () in the space of all CV states, that allows to formulate valuable bounds on suitable entanglement measures and entropic degrees for a general non-Gaussian state , based on the corresponding (easier to compute) properties of the Gaussian state with the same first and second statistical moments as . This has important consequences for the security of CV quantum key distribution qkdextra ().

The transmission of one-mode and multi-mode, possibly entangled beams between distant locations, a basic necessity for the realization of a distributed quantum communication network qinternet (), is unavoidably affected by various types of noise. While phase diffusion (dephasing) dampens the coherences in the Fock basis transforming Gaussian states into non-Gaussian ones, the exposition to dissipative losses and thermal hopping yields an instance of a Gaussian channel, that preserves the Gaussian character of the states sent through it.

Recently, Allegra, Giorda and Paris have attempted a comparison (by means of several separability criteria) of the decoherence effects due to the latter types of channels, on special families of Gaussian as well as non-Gaussian two-mode entangled states, all initially having the same mean energy allegra (). Among the main findings of their work, they formulate a conjecture on the “robustness of Gaussian entanglement”, according to which initial Gaussian pure states should be the ones whose entanglement is the last to vanish, compared to other classes of initial non-Gaussian pure states with equal energy, when they are all evolving through a noisy Markovian channel with loss and thermal hopping.

In this Brief Report, we prove the conjecture of Ref. allegra () by resorting to the extremality of Gaussian states extra (), in particular to the fact that Gaussian states maximize entanglement among all CV states with fixed energy vanenk (). Our simple result establishes in full generality that no advantage can be gained by encoding information into entangled non-Gaussian states when they have to be transmitted through Gaussian noisy channels (under dissipation and/or amplification), and that Gaussian states are, hence, the most robust CV states indeed against photon losses and thermal hopping. An outlook on the problem of identification of the most robust CV entangled states for more general bosonic noisy evolutions, including dephasing channels, is provided in detail.

We consider general pure two-mode CV states with fixed initial mean energy , where , , are the ladder operators of mode satisfying the canonical commutation relations . Among all the possible input states, the two-mode squeezed state , where and , is the only Gaussian instance (up to local unitary operations). We let each possible initial state undergo a dissipative evolution through the channel described by the following master equation allegra (),

(1) |

that encompasses losses and thermal hopping in the presence of nonclassically fluctuating local environments. The dot denotes time-derivative and the Lindblad superoperator is defined as . Here is a loss coefficient and is the mean photon number in the stationary (thermal) state of each mode. The master equation (1) admits the operator solution allegra (); gmd94 (): , where denotes the evolution map corresponding to the noisy channel; are two additional fictitious modes in a thermal state ; and amounts to a beam-splitter transformation with .

The action of the noisy evolution in Eq. (1) can be thus simply modeled as letting each mode interfere with an ancillary thermal mode via a beam-splitter with transmissivity , and eventually tracing over the ancillae, as depicted schematically in Fig. 1. In this notation, () is the mean photon number of the ancillary mode (). In the Heisenberg picture, the beam-splitter transformation acts on the mode operators as follows nielsen ()

(2) |

where are the annihilation operators of the ancillary modes. We can calculate the mean energy of the generic state at any time during the evolution:

(3) | |||||

where all the additional cross-terms vanish because the ancillary thermal states have zero first moments, entailing . What is important here, is that the energy of the evolving states at any time only depends on the initial energy and on the channel parameters, i.e., if a set of states enter the channel all with the same energy, as it is the case in allegra () and in the present discussion, their energies will remain equal to each other (yet globally affected by the open system dynamics) throughout the whole evolution, regardless of the specific form of the states and of their Gaussian or non-Gaussian nature. Let us also remark once more that the channel in Eq. (1) is a Gaussian channel, i.e. the state evolved from the two-mode squeezed state will remain Gaussian during the whole evolution. The final important ingredient we need is the extremality of Gaussian states. If one adopts a bona fide measure of entanglement pv (), which is continuous and strongly superadditive, such as the distillable entanglement, or the squashed entanglement squashed (), then the Gaussian state with energy is the most entangled among all possible CV states with the same energy extra (); vanenk (), as a simple consequence of the maximum entropy property of Gaussian states holevo ().

Collecting all the above observations, we can conclude that at any time , in other words the Gaussian two-mode squeezed state starts in advantage, in terms of entanglement, over any other CV state with the same energy, and maintains its advantage throughout the whole noisy evolution. Given that (as soon as ) the channel in Eq. (1) takes a finite “separation time” to erase all the entanglement in any initial quantum state allegra (), the entanglement in the Gaussian state must be the last one to vanish, i.e., for all CV states with the same mean energy as . This proves the robustness of Gaussian entanglement conjecture of Ref. allegra (), for continuous and strongly superadditive bipartite entanglement monotones, in the general case of and without the need to restrict to be in the form of a photon-number-entangled-state (i.e., a two-mode state whose one-mode marginal density matrices are diagonal in the Fock basis). Actually, being essentially a consequence of extremality, the present result holds for any set of, not necessarily pure, iso-energetic input states of an arbitrary number of modes that include a Gaussian one , and for general Markovian Gaussian channels.

All the speculations of Ref. allegra () that were based on this conjecture can be now regarded as valid and truthful: Gaussian entanglement is extremal in terms of robustness against decoherence due to noise and dissipation. It is important to notice, however, that if entanglement is measured e.g. by the (logarithmic) negativity pv (), which does not satisfy the hypotheses of the theorem in extra (), then the extremality of Gaussian states can be violated by tiny amounts vanenk (); extra (), hence narrow regimes might be found where a particular non-Gaussian state could overcome in the entanglement hierarchy based on negativities. This is nonetheless to be regarded more as a pathological feature of the employed entanglement measure, rather than as a physically meaningful violation or counterexample to the robustness conjecture. Let us remark once more that our general proof in fact holds for any continuous and strongly superadditive bipartite entanglement monotone (such as the squashed entanglement squashed ()). It becomes quite clear, henceforth, how the degree of non-Gaussianity – measured as distance from a Gaussian reference state with the same covariance matrix nonG () – of different instances of ’s is set to play a key role in assessing their robustness to noise and their entanglement evolution, as evidenced by several examples in Ref. allegra ().

Our analysis generalizes the one in Ref. vanenk (), where the two-mode squeezed Gaussian state had already been identified as “the most robust state of light” against evolutions in purely dissipative channels [i.e. with in Eq. (1)]. These results do not straightforwardly extend to more general scenarios in which other types of noisy channels are taken into account, that can violate Gaussianity and/or Markovianity. In these cases, the Gaussian nature of transmitted states is not preserved, and the state evolved from immediately loses its privileged role for , entering the competition among the ’s as a peer. This is the case, for instance, if purely local phase diffusive decoherence channels are considered, which are non-Gaussian channels modeled by a master equation of the form

(4) |

For such noise models, the coherences in an initial two-mode squeezed state are degraded undergoing a non-Gaussian evolution, eventually driving the two modes into an uncorrelated product of thermal states, each corresponding to the reduced one-mode density matrix of the input two-mode state. In this scenario, there is no apparent hierarchy regulating the dynamical comparison between Gaussian and non-Gaussian initial states , when energy or alternative resources (i.e., the squeezing degree) are fixed at time . More complicate is the situation when a two-mode dephasing channel, with correlated noises affecting each mode, is taken into account; in that case, not even the marginal states are preserved.

To formulate an outlook, we feel that even more challenging, yet surely important, would be to investigate the case of more realistic channels where dissipation, thermal hopping, as well as phase diffusion and possibly memory effects and correlated noise may all simultaneously take place. Some of these mechanisms are currently investigated, for instance, in the context of unveiling quantum coherence effects and noise-assisted energy transport phenomena in biological systems caruso (). The problem of entanglement evolution and robustness turns thus into a highly non-trivial arena, and since the powerful extremality results extra () are of little or no use at all in these more general cases, a numerical approach along the lines of Ref. allegra () should be the preferable avenue to pursue. This will be the subject of further work. The main result of this Brief Report should then be recovered in the limit of negligible dephasing, uncorrelated noises and memoryless channels, but we feel that no reliable prediction can be made at this stage about which type of state, in the uncountable arena of CV states, might gain the crown of the most robust state of light in the intermediate regime where all the diverse noise effects are comparable.

Fruitful discussions with Matteo Paris and Alex Monras are acknowledged.

## References

- (1) N. Cerf, G. Leuchs, and E. S. Polzik eds., Quantum Information with Continuous Variables of Atoms and Light (Imperial College Press, London, 2007).
- (2) S. L. Braunstein and P. van Loock, Rev. Mod. Phys. 77, 513 (2005).
- (3) A. Furusawa, J. L. Sørensen, S. L. Braunstein, C. A. Fuchs, H. J. Kimble, and E. S. Polzik, Science 282, 706 (1998).
- (4) F. Grosshans, G. Van Assche, J. Wenger, R. Brouri, N. J. Cerf, and P. Grangier, Nature 421, 238 (2003).
- (5) R. Ukai, N. Iwata, Y. Shimokawa, S. C. Armstrong, A. Politi, J. Yoshikawa, P. van Loock, and A. Furusawa, arXiv:1001.4860.
- (6) G. Adesso and F. Illuminati, J. Phys. A 40 7821, (2007).
- (7) J. Eisert, S. Scheel, and M. B. Plenio, Phys. Rev. Lett. 89, 137903 (2002); J. Fiurás̆ek, ibid. 89, 137904 (2002); J. Niset, J. Fiurasek, and N. J. Cerf, ibid. 102, 120501 (2009).
- (8) A. I. Lvovsky, H. Hansen, T. Aichele, O. Benson, J. Mlynek, and S. Schiller, Phys. Rev. Lett. 87, 050402 (2001); A. Zavatta, S. Viciani, and M. Bellini, Science 306, 660 (2004); J. Wenger, R. Tualle-Brouri, and P. Grangier, Phys. Rev. Lett. 92, 153601 (2004); J. S. Neergaard-Nielsen, B. Melholt Nielsen, C. Hettich, K. Mølmer, and E. S. Polzik, Phys. Rev. Lett. 97, 083604 (2006); K. Wakui, H. Takahashi, A. Furusawa, and M. Sasaki, Opt. Express 15, 3568 (2007); A. Ourjoumtsev, H. Jeong, R. Tualle-Broui, and P. Grangier, Nature 448, 784 (2007).
- (9) T. Opatrny, G. Kurizki, and D.-G. Welsch, Phys. Rev. A 61, 032302 (2000); P. T. Cochrane, T. C. Ralph, and G. J. Milburn, ibid. 65, 062306 (2002); S. Olivares, M. G. A. Paris, and R. Bonifacio, ibid. 67, 032314 (2003); A. Kitagawa, M. Takeoka, M. Sasaki, and A. Chefles, Phys. Rev. A 73, 042310 (2006); F. Dell’Anno, S. De Siena, L. Albano, and F. Illuminati, Phys. Rev. A 76, 022301 (2007).
- (10) J. Eisert, D. E. Browne, S. Scheel, and M. B. Plenio, Ann. Phys. 311, 431 (2004); A. Ourjoumtsev, A. Dantan, R. Tualle-Brouri, and P. Grangier, Phys. Rev. Lett. 98, 030502 (2007). H. Takahashi, S. Neergaard-Nielsen, M. Takeuchi, M. Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, Nature Photon. 4, 178 (2010).
- (11) G. Adesso, F. Dell’Anno, S. De Siena, F. Illuminati, and L. A. M. Souza, Phys. Rev. A 79, 040305(R) (2009).
- (12) S. Lloyd and S. L. Braunstein, Phys. Rev. Lett. 82, 1784 (1999); N. C. Menicucci, P. van Loock, M. Gu, C. Weedbrook, T. C. Ralph, and M. A. Nielsen, Phys. Rev. Lett. 97, 110501 (2006).
- (13) H. Nha and H. J. Carmichael, Phys. Rev. Lett. 93, 020401 (2004); R. Garcia-Patron, J. Fiurasek, N. J. Cerf, J. Wenger, R. Tualle-Brouri, and P. Grangier, ibid. 93, 130409 (2004).
- (14) M. M. Wolf, G. Giedke, and J. I. Cirac, Phys. Rev. Lett. 96, 080502 (2006).
- (15) M. Navascues, F. Grosshans, and A. Acin, Phys. Rev. Lett. 97, 190502 (2006); R. Garcia-Patron and N. J. Cerf, ibid. 97, 190503 (2006).
- (16) H. J. Kimble, Nature 453, 1023 (2008).
- (17) M. Allegra, P. Giorda, and M. G. A. Paris, Phys. Rev. Lett. 105, 100503 (2010).
- (18) S. J. van Enk and O. Hirota, Phys. Rev. A 71, 062322 (2005).
- (19) G.M. D’Ariano, Phys. Lett. A 187, 231 (1994).
- (20) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cambridge, 2000).
- (21) M. B. Plenio and S. Virmani, Quant. Inf. Comp. 7, 1 (2007).
- (22) M. Christandl and A. Winter, J. Math. Phys. 45, 829 (2004).
- (23) A. S. Holevo, M. Sohma, and O. Hirota, Phys. Rev. A 59, 1820 (1999).
- (24) M. G. Genoni, M. G. A. Paris, and K. Banaszek, Phys. Rev. A 78, 060303(R) (2008).
- (25) F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and M. B. Plenio, J. Chemi. Phys. 131, 105106 (2009).