# Entanglement and discord in two-mode Gaussian open quantum systems

###### Abstract

In the framework of the theory of open systems, we give a description of quantum entanglement and quantum discord for two non-interacting modes embedded in a thermal environment. We describe the evolution of entanglement in terms of the covariance matrix for Gaussian input states. For all values of the temperature of the thermal reservoir, an initial separable squeezed thermal state remains separable for all times. In the case of an entangled initial squeezed thermal state, entanglement suppression (entanglement sudden death) takes place, for non-zero temperatures of the environment. The Gaussian quantum discord, which is a measure of all quantum correlations in the state, including entanglement, decays asymptotically in time.

###### pacs:

03.65.Yz, 03.67.Bg, 03.67.Mn## 1 Introduction

In the framework of the theory of open systems based on completely positive quantum dynamical semigroups, we study the dynamics of quantum entanglement and quantum discord of a subsystem consisting of two uncoupled modes interacting with a common thermal bath. The initial state of the subsystem is taken of Gaussian form and the evolution under the dynamical semigroup assures the preservation in time of the Gaussian form of the state. In Sec. 2 we give the solution of the evolution equation for the covariance matrix corresponding to a generic two-mode Gaussian state of the two uncoupled modes interacting with the environment. We investigate in Sec. 3 the dynamics of entanglement, by analyzing the time evolution of the logarithmic negativity, which characterizes the degree of entanglement of the state. We study also the time evolution of the Gaussian quantum discord, which is a measure of all quantum correlations in the bipartite state including entanglement. Summary is given in Sec. 4.

## 2 Equations of motion for two modes interacting with an environment

We study the dynamics of a subsystem composed of two non-interacting modes in weak interaction with a thermal environment. In the axiomatic formalism based on completely positive quantum dynamical semigroups, the Markovian irreversible time evolution of an open system is described by the Lindblad-Kossakowski master equation. We are interested in the set of Gaussian states, therefore we introduce such quantum dynamical semigroups that preserve this set during time evolution of the system.

The two-mode Gaussian state is entirely specified by its covariance matrix, which is a real, symmetric and positive matrix:

(1) |

where , and are Hermitian matrices. and denote the symmetric covariance matrices for the individual reduced one-mode states, while the matrix contains the cross-correlations between modes. The equations of motion for the quantum correlations of the coordinates and momenta of the two modes are the following ( denotes the transposed matrix) [1]:

(2) |

where

(3) | |||

(4) |

is the mass, and are the frequencies of the non-resonant modes, and diffusion coefficients … and dissipation constant are real quantities.

## 3 Dynamics of continuous variable entanglement and discord

### 3.1 Time evolution of entanglement and logarithmic negativity

In order to quantify the degree of entanglement of the two-mode states it is suitable to use the logarithmic negativity. For a Gaussian density operator, the logarithmic negativity is given by where is the smallest of the two symplectic eigenvalues of the partial transpose of the 2-mode covariance matrix . In our model, the logarithmic negativity is calculated as [2, 3, 4, 5]

(7) |

where

(8) |

We assume that the initial Gaussian state is the two-mode squeezed thermal state, with the covariance matrix of the form [6]

(9) |

with the matrix elements given by

(10) | |||

(11) | |||

(12) |

where are the average numbers of thermal photons associated with the two modes and denotes the squeezing parameter. In the particular case and , (9) becomes the covariance matrix of the two-mode squeezed vacuum state. A two-mode squeezed thermal state is entangled when the squeezing parameter satisfies the inequality [6], where

(13) |

We suppose that the asymptotic state of the considered open system is a Gibbs state corresponding to two independent modes in thermal equilibrium at temperature Then the quantum diffusion coefficients have the following form (we put ) [7]:

(14) | |||

The calculations with separable initial squeezed thermal states show that these states remain separable for all values of the temperature and for all times. The evolution of entangled initial squeezed thermal states with the covariance matrix given by Eq. (9) is illustrated in Figure 1, where we represent the dependence of the logarithmic negativity on time and temperature for the case of an initial symmetric Gaussian state, when . We observe that for a non-zero temperature at certain finite moment of time, which depends on becomes zero and therefore the state becomes separable. This is the so-called phenomenon of entanglement sudden death. It is in contrast to the quantum decoherence, during which the loss of quantum coherence is usually gradual [8, 9]. For remains strictly positive for finite times and tends asymptotically to 0 for Therefore, only for zero temperature of the thermal bath the initial entangled state remains entangled for finite times and this state tends asymptotically to a separable one for infinitely large time. Dissipation favorizes the phenomenon of entanglement sudden death – with increasing the dissipation parameter the entanglement suppression happens earlier. The same qualitative behaviour of the time evolution of entanglement was obtained previously [10, 11] in the particular case and corresponding to an initial two-mode squeezed vacuum state.

### 3.2 Gaussian quantum discord

Quantum discord was introduced [12] as a measure of all quantum correlations in a bipartite state, including – but not restricted to – entanglement. Originally the quantum discord was defined and evaluated mainly for finite dimensional systems. Very recently [13, 14] the notion of discord has been extended to the domain of continuous variable systems, in particular to the analysis of bipartite systems described by two-mode Gaussian states. Closed formulas have been derived for bipartite thermal squeezed states [13] and for all two-mode Gaussian states [14]. The Gaussian quantum discord of a general two-mode Gaussian state can be defined as the quantum discord where the conditional entropy is restricted to generalized Gaussian positive operator valued measurements (POVM) on the mode 2, and in terms of symplectic invariants it is given by (the symmetry between the two modes 1 and 2 is broken) [14]

(15) |

where

(16) |

(17) |

(18) |

and are the symplectic eigenvalues of the state, given by

(19) |

where Notice that Gaussian quantum discord only depends on , i.e., entangled () and separable states are treated on equal footing.

The evolution of the Gaussian quantum discord is illustrated in Figure 2, where we represent the dependence of on time and temperature for an entangled initial symmetric Gaussian state, taken of the form of a two-mode squeezed thermal state (9), for such values of the parameters which satisfy for all times the first condition in formula (17). For entangled initial states the Gaussian discord remains strictly positive in time and in the limit of infinite time it tends asymptotically to zero, corresponding to the thermal product (separable) state, with no correlations at all. The Gaussian discord has nonzero values for all finite times and this fact certifies the existence of nonclassical correlations in two-mode Gaussian states, either separable or entangled.

From Figure 2 we notice that, in concordance with the general properties of the Gaussian quantum discord [14], the states can be either separable or entangled for and all the states above the threshold are entangled. We also notice that the decay of quantum discord is stronger when the temperature is increasing.

## 4 Summary

In the framework of the theory of open systems we investigated the Markovian dynamics of quantum correlations for a subsystem composed of two non-interacting modes embedded in a thermal bath. We have described the time evolution of the logarithmic negativity in terms of the covariance matrix for squeezed thermal states, for the case when the asymptotic state is a Gibbs state corresponding to two independent quantum harmonic oscillators in thermal equilibrium. For all values of the temperature of the thermal reservoir, an initial separable Gaussian state remains separable for all times. In the case of an entangled initial squeezed thermal state, entanglement sudden death takes place for non-zero temperatures of the environment. Only for a zero temperature of the thermal bath the initial entangled state remains entangled for finite times, but in the limit of infinite time it evolves asymptotically to an equilibrium state which is always separable. The Gaussian quantum discord decreases asymptotically in time. This is contrast with the sudden death of entanglement. The time evolution of quantum discord is very similar to that of the entanglement before the sudden suppression of the entanglement. After the sudden death of the entanglement, the nonzero values of quantum discord quantify the nonclassical correlations for separable mixed states.

## References

## References

- [1] Sandulescu A, Scutaru H and Scheid W 1987 J. Phys. A: Math. Gen. 20 2121
- [2] Isar A 2008 Eur. J. Phys. Special Topics 160 225
- [3] Isar A 2008 Int. J. Quantum Inf. 6 689
- [4] Isar A 2009 Phys. Scr., Topical Issue 135 014033
- [5] Isar A 2009 Open Sys. Inf. Dynamics 16 205
- [6] Marian P, Marian TA and Scutaru H 2003 Phys. Rev. A 68 062309
- [7] Isar A, Sandulescu A, Scutaru H, Stefanescu E and Scheid W 1994 Int. J. Mod. Phys. E 3 635
- [8] Isar A and Scheid W 2007 Physica A 373 298
- [9] Isar A 2007 J. Russ. Laser Res. 28 439
- [10] Isar A 2011 Phys. Scr. T143 014012
- [11] Isar A 2011 Open Sys. Inf. Dynamics 18 1
- [12] Zurek WH 2000 Annalen der Physik (Leipzig) 9 853
- [13] Giorda P and Paris MGA 2010 Phys. Rev. Lett. 105 020503
- [14] Adesso G and Datta A 2010 Phys. Rev. Lett. 105 030501