text:
# Entanglement induced by spontaneous emission in spatially extended two-atom systems

## I Introduction

The subject of generation of entangled states has attracted much interest since it became clear that entanglement lies at the heart of many new applications which come under the general heading of quantum information and quantum computation. Several methods of creating entanglement have been proposed involving trapped and cooled ions or neutral atoms . Of particular interest is generation of entangled states in two-atom systems, since it is generally believed that entanglement of only two microscopic quantum systems (two qubits) is essential to implement quantum protocols such as quantum computation. It has been shown that entangled states in a two-atom system can be created by a continuous driving of the atoms with a coherent or chaotic thermal field , or by spontaneous emission from two distant atoms initially prepared in a coherent superposition state or in an entangled state . In particular, Cabrillo et al.  demonstrated that two three-level atoms initially prepared by a weak pulse in an entangled state can evolve under spontaneous emission into an entangled state of the ground states of the atoms.

The effect of spontaneous emission on entanglement creation has also been discussed by others . These studies, however, have been limited to the small sample (Dicke) model . The disadvantage of the Dicke model is that it does not include the dipole-dipole interaction among the atoms and does not correspond to realistic experimental situations of atoms located (trapped) at different positions. In fact, the model corresponds to a very specific geometrical configuration of the atoms confined to a volume much smaller compared with the atomic resonant wavelength (the small-sample model). The present atom trapping and cooling techniques can trap two atoms at distances of order of a resonant wavelength , which makes questionable the applicability of the Dicke model to physical systems.

In this paper we study what may be termed ”spontaneously induced transient entanglement” in a system of two interacting atoms. In related works, Kim et al.  and Zhou et al.  have shown that a transient entanglement can be created in the Dicke model driven by a thermal (chaotic) field. Sørensen et al.  have proposed a method to produce a transient entanglement in Bose-Einstein condensate excited by a single pulse. Turchette et al.  have recently realised experimentally a transient entanglement in two trapped ions. Unlike previous works , we consider entanglement creation by spontaneous emission from initially uncorrelated atoms and without the presence of external coherent or incoherent fields. We are particularly interested in three aspects of entanglement creation by spontaneous emission: (1) The dependence of the entanglement creation on specific arrangements of initial uncorrelated states, (2) the role of the antisymmetric state in the entanglement creation, and (3) sharing and transfer of entanglement between two entangled states. We do not make the small-sample approximation, so that our results are valid for arbitrary interatomic separations. The antisymmetric state and its dynamics are not neglected, even when dealing with small interatomic separations. We show that spontaneous emission from two spatially separated atoms can lead to a transient entanglement of initially unentangled atoms. This result contrasts the Dicke model where spontaneous emission cannot produce entanglement from initially unentangled atoms . Moreover, we show that the entanglement creation relies crucially on the population distribution between the entangled symmetric and antisymmetric states and attains maximal values when the population of the symmetric state becomes zero. This is a rather surprising prediction, since the symmetric state is an example of maximally entangled state and one might conclude that its participation in the atomic dynamics would enhance entanglement.

## II Collective two-atom systems

We consider a system of two non-overlapping two-level atoms with ground states $`|g_i`$ and excited states $`|e_i(i=1,2)`$ connected by dipole transition moments $`\stackrel{}{\mu }_i`$. The atoms are located at fixed positions $`\stackrel{}{r}_1`$ and $`\stackrel{}{r}_2`$ and coupled to all modes of the electromagnetic field, which we assume are in the vacuum state. We consider spontaneous emission from identical as well as non-identical atoms prepared in two different initial states. In the case of nonidentical atoms, we assume that atoms have equal dipole moments $`\stackrel{}{\mu }_1=\stackrel{}{\mu }_2=\stackrel{}{\mu }`$, but different transition frequencies $`\omega _1`$ and $`\omega _2`$, such that $`\omega _2\omega _1\omega _0=(\omega _1+\omega _2)/2`$, so that the rotating-wave approximation can be applied to calculate the dynamics of the system.

The time evolution of the system of atoms coupled through the vacuum field is given by the following master equation

$`{\displaystyle \frac{\widehat{\rho }}{t}}`$ $`=`$ $`i{\displaystyle \underset{i=1}{\overset{2}{}}}\omega _i[S_i^z,\widehat{\rho }]i{\displaystyle \underset{ij}{\overset{2}{}}}\mathrm{\Omega }_{ij}[S_i^+S_j^{},\widehat{\rho }]`$ (1)
$``$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{i,j=1}{\overset{2}{}}}\mathrm{\Gamma }_{ij}\left(\widehat{\rho }S_i^+S_j^{}+S_i^+S_j^{}\widehat{\rho }2S_j^{}\widehat{\rho }S_i^+\right),`$ (2)
where $`S_i^+(S_i^{})`$ are the dipole raising (lowering) operators and $`S^z`$ is the energy operator of the $`i`$th atom. In Eq. (2), $`\mathrm{\Gamma }_{ij}(i=j)`$ are the spontaneous emission rates of the atoms, equal to the Einstein $`A`$ coefficient for spontaneous emission, whereas $`\mathrm{\Gamma }_{ij}`$ and $`\mathrm{\Omega }_{ij}(ij)`$ describe the interatomic coupling , and are the collective damping and the dipole-dipole interaction potential defined, respectively, by

$`\mathrm{\Gamma }_{ij}=\mathrm{\Gamma }_{ji}`$ $`=`$ $`{\displaystyle \frac{3}{2}}\mathrm{\Gamma }\{[1(\overline{\mu }\overline{r}_{ij})^2]{\displaystyle \frac{\mathrm{sin}\left(k_0r_{ij}\right)}{k_0r_{ij}}}`$ (4)
$`+[13(\overline{\mu }\overline{r}_{ij})^2][{\displaystyle \frac{\mathrm{cos}\left(k_0r_{ij}\right)}{\left(k_0r_{ij}\right)^2}}{\displaystyle \frac{\mathrm{sin}\left(k_0r_{ij}\right)}{\left(k_0r_{ij}\right)^3}}]\},`$
and

$`\mathrm{\Omega }_{ij}`$ $`=`$ $`{\displaystyle \frac{3}{4}}\mathrm{\Gamma }\{[1(\overline{\mu }\overline{r}_{ij})^2]{\displaystyle \frac{\mathrm{cos}\left(k_0r_{ij}\right)}{k_0r_{ij}}}`$ (6)
$`+[13(\overline{\mu }\overline{r}_{ij})^2][{\displaystyle \frac{\mathrm{sin}\left(k_0r_{ij}\right)}{\left(k_0r_{ij}\right)^2}}+{\displaystyle \frac{\mathrm{cos}\left(k_0r_{ij}\right)}{\left(k_0r_{ij}\right)^3}}]\},`$
where $`k_0=\omega _0/c`$, $`r_{ij}=\left|\stackrel{}{r}_j\stackrel{}{r}_i\right|`$ is the distance between the atoms, $`\overline{\mu }`$ is unit vector along the atomic transition dipole moments, that we assume are parallel to each other, and $`\overline{r}_{ij}`$ is the unit vector along the interatomic axis.

The master equation (2) has been used for many years to study a wide variety of problems involving the interaction of collective atomic systems with the radiation field . Using the master equation (2), we can write down the equations of motion for the components of the density matrix of the two-atom system in the basis of the product states $`|e_1|e_2`$, $`|e_1|g_2`$, $`|g_1|e_2`$ and $`|g_1|g_2`$ of the individual atoms. However, the problem simplifies by working in the basis of the collective states of the system which contains symmetric and antisymmetric combinations of the product states. For identical atoms $`(\omega _1=\omega _2)`$ the collective states are

$`|e`$ $`=`$ $`|e_1|e_2,`$ (7)
$`|s`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left(|e_1|g_2+|g_1|e_2\right),`$ (8)
$`|a`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left(|e_1|g_2|g_1|e_2\right),`$ (9)
$`|g`$ $`=`$ $`|g_1|g_2.`$ (10)
In the collective state representation, the two-atom system behaves as a single four-level system with the ground state $`|g`$, the upper state $`|e`$, and two intermediate states: the symmetric $`|s`$ and antisymmetric $`|a`$ states. The most important property of the collective states is that the symmetric and antisymmetric states are maximally entangled states. The states are linear superpositions of the product states which cannot be separated into product states of the individual atoms.

For non-identical atoms, the collective states of the system contain non-maximally entangled states, which can be written as linear combinations of the maximally entangled states

$`|e`$ $`=`$ $`|e_1|e_2,`$ (11)
$`|s^{}`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left[\left(\alpha +\beta \right)|s+\left(\beta \alpha \right)|a\right],`$ (12)
$`|a^{}`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left[\left(\alpha \beta \right)|s+\left(\alpha +\beta \right)|a\right],`$ (13)
$`|g`$ $`=`$ $`|g_1|g_2,`$ (14)
where $`\alpha =d/\sqrt{d^2+\mathrm{\Omega }_{12}^2},\beta =\mathrm{\Omega }_{12}/\sqrt{d^2+\mathrm{\Omega }_{12}^2}`$, $`d=\mathrm{\Delta }+\sqrt{\mathrm{\Omega }_{12}^2+\mathrm{\Delta }^2}`$, and $`\mathrm{\Delta }=(\omega _2\omega _1)/2`$.

Thus, in both cases of identical or non-identical atoms, we can limit the considerations to the basis of the collective states (10). In this basis, the density matrix elements satisfy the following set of simple differential equations

$`\dot{\rho }_{ee}`$ $`=`$ $`2\mathrm{\Gamma }\rho _{ee},`$ (15)
$`\dot{\rho }_{ss}`$ $`=`$ $`\left(\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)\left(\rho _{ss}\rho _{ee}\right)+i\mathrm{\Delta }\left(\rho _{as}\rho _{sa}\right),`$ (16)
$`\dot{\rho }_{aa}`$ $`=`$ $`\left(\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)\left(\rho _{aa}\rho _{ee}\right)i\mathrm{\Delta }\left(\rho _{as}\rho _{sa}\right),`$ (17)
$`\dot{\rho }_{as}`$ $`=`$ $`\left(\mathrm{\Gamma }+2i\mathrm{\Omega }_{12}\right)\rho _{as}+i\mathrm{\Delta }\left(\rho _{ss}\rho _{aa}\right),`$ (18)
$`\dot{\rho }_{se}`$ $`=`$ $`\left[{\displaystyle \frac{1}{2}}\left(3\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)i\left(\omega _0\mathrm{\Omega }_{12}\right)\right]\rho _{se}+i\mathrm{\Delta }\rho _{ae},`$ (19)
$`\dot{\rho }_{ae}`$ $`=`$ $`\left[{\displaystyle \frac{1}{2}}\left(3\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)i\left(\omega _0+\mathrm{\Omega }_{12}\right)\right]\rho _{ae}+i\mathrm{\Delta }\rho _{se},`$ (20)
$`\dot{\rho }_{gs}`$ $`=`$ $`\left[{\displaystyle \frac{1}{2}}\left(\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)i\left(\omega _0+\mathrm{\Omega }_{12}\right)\right]\rho _{gs}+\left(\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)\rho _{se}i\mathrm{\Delta }\rho _{ga},`$ (21)
$`\dot{\rho }_{ga}`$ $`=`$ $`\left[{\displaystyle \frac{1}{2}}\left(\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)i\left(\omega _0\mathrm{\Omega }_{12}\right)\right]\rho _{ga}\left(\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)\rho _{ae}i\mathrm{\Delta }\rho _{gs},`$ (22)
$`\dot{\rho }_{eg}`$ $`=`$ $`\left(\mathrm{\Gamma }+2i\omega _0\right)\rho _{eg}.`$ (23)
Equations (23) show that all transitions rates to and from the symmetric state are equal to $`(\mathrm{\Gamma }+\mathrm{\Gamma }_{12})`$. On the other hand, all transitions rates to and from the antisymmetric state are equal to $`(\mathrm{\Gamma }\mathrm{\Gamma }_{12})`$. Thus, the symmetric state decays with an enhanced (superradiant) rate, while the antisymmetric state decays with a reduced (subradiant) state. Hence, the population of the antisymmetric state experiences a variation on a time scale of order $`(\mathrm{\Gamma }\mathrm{\Gamma }_{12})^1`$, which can lead to interesting effects not observed in the Dicke model. These effects result from the fact that the set of equations (23) has two different solutions depending on whether $`\mathrm{\Gamma }_{12}=\mathrm{\Gamma }`$ or $`\mathrm{\Gamma }_{12}\mathrm{\Gamma }`$. The case of $`\mathrm{\Gamma }_{12}=\mathrm{\Gamma }`$ corresponds to the small sample (Dicke) model, whereas the case of $`\mathrm{\Gamma }_{12}\mathrm{\Gamma }`$ corresponds to spatially extended atomic systems. The existence of two different solutions of Eq. (23) is connected with conservation of the total spin $`S^2`$, that $`S^2`$ is a constant of motion for the Dicke model and $`S^2`$ not being a constant of motion for a spatially extended system of atoms . We can explain it by expressing the square of the total spin of the two-atom system in terms of the density matrix elements of the collective system as

$`S^2\left(t\right)=22\rho _{aa}\left(t\right).`$ (24)
It is clear from Eq. (24) that $`S^2`$ is conserved only in the Dicke model, in which the antisymmetric state is ignored. For a spatially extended system the antisymmetric state participates fully in the dynamics and $`S^2`$ is not conserved. The Dicke model evolves between the triplet states $`|e`$, $`|s`$, and $`|g`$, while the spatially extended two-atom system evolves between the triplet and the antisymmetric states.

## III Transient entanglement

The entanglement creation by spontaneous emission is illustrated most clearly if one assumes that a system of two atoms decays spontaneously from initially unentangled (uncorrelated) states. Several different measures have been proposed to identify entanglement between two atoms, and we choose the Wootters entanglement measure , the concurrence $`C`$, defined as

$`C=\mathrm{max}(0,\sqrt{\lambda _1}\sqrt{\lambda _2}\sqrt{\lambda _3}\sqrt{\lambda _4}),`$ (25)
where $`\lambda _1,\mathrm{},\lambda _4`$ are the eigenvalues of the matrix $`\stackrel{~}{\rho }=\rho (\sigma _y\sigma _y)\rho ^{}(\sigma _y\sigma _y)`$ and $`\sigma _y`$ is the Pauli matrix. The range of the concurrence is from 0 to 1. For unentangled atoms $`C=0`$ whereas $`C=1`$ for the maximally entangled atoms. The concurrence involves eigenvalues of the complicated matrix $`\stackrel{~}{\rho }`$ and, in general, is difficult to calculate analytically. Therefore, for the understanding and explanation of the entanglement creation via spontaneous emission, we will use the Peres-Horodecki (negativity) measure for entanglement . The negativity criterion is given by the quantity

$`E=\mathrm{max}(0,2{\displaystyle \underset{i}{}}\mu _i),`$ (26)
where the sum is taken over the negative eigenvalues $`\mu _i`$ of the partial transposition of the density matrix $`\widehat{\rho }`$ of the system. The value $`E=1`$ corresponds to maximum entanglement between the atoms whilst $`E=0`$ describes completely separated atoms.

The two entanglement measures, the concurrence $`C`$ and negativity $`E`$, give the same results for criteria for entanglement, but they can give different results for a degree of entanglement .

### A Identical atoms

We begin with spontaneous emission from two identical atoms and consider entanglement creation for two different sets of initial atomic conditions at $`t=0`$. In the first, one of the atoms is in its excited state and the other is assumed to resides in its ground state. In the second, both atoms are assumed to reside in their excited states.

#### 1 Initial state of only one atom excited: $`|\mathrm{\Phi }_0=|e_1|g_2`$.

In the basis of the collective states of the system, the initial condition of only one atom excited corresponds to the initial condition with non-zero density matrix elements $`\rho _{ss}(0)=\rho _{aa}(0)=\rho _{sa}(0)=\rho _{as}(0)=1/2`$. In Fig. 1, we plot the time evolution of the populations $`\rho _{ss}(t)`$, $`\rho _{aa}(t)`$ and the concurrence $`C`$, which we have found solving the equations of motion (23) with the condition that initially at $`t=0`$ one of the atoms was prepared in its excited state and the other in the ground state. One can see from Fig. 1 that the concurrence is zero at $`t=0`$; in other words there is no entanglement in the system at $`t=0`$. The time evolution of the concurrence reflects the time evolution of the entanglement between the atoms. It is seen from Fig. 1 that the concurrence builds up as time develops indicating that spontaneous emission can create entanglement between the initially unentangled atoms. As time progresses, the concurrence rapidly increases and reaches a maximum at time close to the point where the symmetric state becomes depopulated. At later times, the concurrence slowly decreases and overlaps with the time evolution of the population of the antisymmetric state.

The physical understanding of the creation of the transient entanglement can be achieved by considering the time-dependent density matrix of the system, which for the initial condition of only one atom excited has the following form

$`\rho (t)`$ $`=`$ $`\left(\begin{array}{cccc}0& 0& 0& 0\\ 0& \rho _{ss}(t)& \rho _{sa}(t)& 0\\ 0& \rho _{as}(t)& \rho _{aa}(t)& 0\\ 0& 0& 0& \rho _{gg}(t)\end{array}\right),`$ (31)
where the non-zero time-dependent density matrix elements are given by

$`\rho _{ss}(t)`$ $`=`$ $`{\displaystyle \frac{1}{2}}\mathrm{exp}\left[\left(\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)t\right],`$ (32)
$`\rho _{aa}(t)`$ $`=`$ $`{\displaystyle \frac{1}{2}}\mathrm{exp}\left[\left(\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)t\right],`$ (33)
$`\rho _{gg}(t)`$ $`=`$ $`1\mathrm{exp}\left(\mathrm{\Gamma }t\right)\mathrm{cosh}\mathrm{\Gamma }_{12}t,`$ (34)
$`\rho _{as}(t)`$ $`=`$ $`\rho _{sa}^{}(t)={\displaystyle \frac{1}{2}}\mathrm{exp}\left[\left(\mathrm{\Gamma }+2i\mathrm{\Omega }_{12}\right)t\right].`$ (35)
One can see that the density matrix is not diagonal due to the presence of coherences $`\rho _{sa}(t)`$ and $`\rho _{as}(t)`$. The density matrix can be rediagonalized to give new ”diagonal” states

$`|\mathrm{\Psi }_1`$ $`=`$ $`\left\{\left[P_1(t)\rho _{ss}(t)\right]|a+\rho _{as}(t)|s\right\}/\left\{\left[P_1(t)\rho _{ss}(t)\right]^2+\left|\rho _{as}(t)\right|^2\right\}^{\frac{1}{2}},`$ (36)
$`|\mathrm{\Psi }_2`$ $`=`$ $`\left\{\rho _{as}(t)|a+\left[P_2(t)\rho _{aa}(t)\right]|s\right\}/\left\{\left[P_2(t)\rho _{aa}(t)\right]^2+\left|\rho _{as}(t)\right|^2\right\}^{\frac{1}{2}},`$ (37)
$`|\mathrm{\Psi }_3`$ $`=`$ $`|g,`$ (38)
$`|\mathrm{\Psi }_4`$ $`=`$ $`|e,`$ (39)
where the diagonal probabilities (populations of the new states) are

$`P_1(t)`$ $`=`$ $`\mathrm{exp}\left(\mathrm{\Gamma }t\right)\mathrm{cosh}\mathrm{\Gamma }_{12}t,`$ (40)
$`P_2(t)`$ $`=`$ $`0,`$ (41)
$`P_3(t)`$ $`=`$ $`1\mathrm{exp}\left(\mathrm{\Gamma }t\right)\mathrm{cosh}\mathrm{\Gamma }_{12}t,`$ (42)
$`P_4(t)`$ $`=`$ $`0.`$ (43)
Thus, the coherences $`\rho _{sa}(t)`$ and $`\rho _{as}(t)`$ cause the system to evolve effectively only between two states: the ground product state $`|g`$ and the entangled state $`|\mathrm{\Psi }_1`$, which is a linear combination of the states $`|s`$ and $`|a`$. In this case, the density matrix of the system is diagonal for all times $`t`$, and is given by

$`\rho (t)=P_1(t)|\mathrm{\Psi }_1\mathrm{\Psi }_1|+\rho _{gg}(t)|gg|.`$ (44)
It is easy to find from Eq. (35) that $`|\rho _{as}(t)|^2=\rho _{aa}(t)\rho _{ss}(t)`$, and then the state $`|\mathrm{\Psi }_1`$ can be written as

$`|\mathrm{\Psi }_1`$ $`=`$ $`\left(\sqrt{\rho _{ss}(t)}|s+\sqrt{\rho _{aa}(t)}|a\right)/\sqrt{\rho _{aa}(t)+\rho _{ss}(t)},`$ (45)
The state $`|\mathrm{\Psi }_1`$ reduces to a nonentangled state when the states $`|s`$ and $`|a`$ have equal populations. On the other hand, the state $`|\mathrm{\Psi }_1`$ reduces to a maximally entangled state when either $`\rho _{ss}(t)`$ or $`\rho _{aa}(t)`$ are equal to zero.

We are now in a position to understand quantitatively the method of creation of entanglement via spontaneous emission. The spontaneously induced entanglement results from unequal populations of the symmetric and antisymmetric states. This observation is supported by the Peres-Horodecki measure for entanglement. It is easy to show that the eigenvalues of the partial transposition of the density matrix (31) are

$`\mu _1`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left[\rho _{ss}(t)+\rho _{aa}(t)\rho _{as}(t)\rho _{sa}(t)\right],`$ (46)
$`\mu _2`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left[\rho _{ss}(t)+\rho _{aa}(t)+\rho _{as}(t)+\rho _{sa}(t)\right],`$ (47)
$`\mu _3`$ $`=`$ $`{\displaystyle \frac{1}{2}}\{\rho _{gg}(t)+\left[\rho _{gg}^2(t)+\left(\rho _{ss}(t)\rho _{aa}(t)\right)^2\left(\rho _{as}(t)\rho _{sa}(t)\right)^2\right]^{\frac{1}{2}}\},`$ (48)
$`\mu _3`$ $`=`$ $`{\displaystyle \frac{1}{2}}\{\rho _{gg}(t)\left[\rho _{gg}^2(t)+\left(\rho _{ss}(t)\rho _{aa}(t)\right)^2\left(\rho _{as}(t)\rho _{sa}(t)\right)^2\right]^{\frac{1}{2}}\}.`$ (49)
It is obvious from Eqs. (49) and (35) that $`\mu _1,\mu _2`$ and $`\mu _3`$ are always positive. The eigenvalue $`\mu _4`$ becomes negative if and only if

$`|\rho _{ss}(t)\rho _{aa}(t)|>|\rho _{as}(t)\rho _{sa}(t)|,`$ (50)
which, according to Eq. (35), is satisfied if $`\rho _{ss}(t)\rho _{aa}(t)`$.

Since the population of the symmetric state decays faster than the antisymmetric state (see Eq. (35)), at time when the state $`|s`$ becomes depopulated, the state $`|\mathrm{\Psi }_1`$ reduces to the maximally entangled antisymmetric state $`|a`$. The above analysis give clear evidence that a transient entanglement created by spontaneous emission can appear only in spatially extended two-atom systems where the antisymmetric state fully participates in the dynamics of the system.

#### 2 Initial state of both atoms excited: $`|\mathrm{\Phi }_0=|e_1|e_2`$.

The role of the antisymmetric state in the entanglement creation is more evident when we choose the initial state of both atoms in their excited states. Here, $`\rho _{ee}(0)=1`$ and the initial values of the remaining density matrix elements are zero. Note, that in this case there are no any initial coherences between the atoms, and also there are no any initial coherences between the collective states. Recent analysis of the entanglement creation in the Dicke model have shown that there is no entanglement creation in spontaneous emission when the atoms are initially prepared in their excited states . In Fig. 2, we show the concurrence $`C`$ and the populations $`\rho _{ss}(t)`$ and $`\rho _{aa}(t)`$ for the case of two spatially separated atoms prepared in their excited states at $`t=0`$. Since, at $`t=0`$, only the upper state $`|e`$ is occupied, one finds from Eq. (23) that the spontaneous emission populates the symmetric and antisymmetric states with different rates. Thus, at early times the symmetric state is more populated than the antisymmetric state. Again, there is no entanglement between the atoms at $`t=0`$ as the initial state $`|e`$ is an unentangled state. However, the figure clearly demonstrates that contrary to the Dicke model, a transient entanglement can be generated in spontaneous emission from two spatially separated atoms. Note that the entanglement is not generated until the population $`\rho _{ss}(t)`$ becomes zero.

This effect can be explained by the Peres-Horodecki measure for entanglement. Consider the time evolution of the density operator of the system, with the initial condition of $`\rho _{ee}(0)=1`$, has the diagonal form for all times

$`\rho (t)=\rho _{ee}(t)|ee|+\rho _{ss}(t)|ss|+\rho _{aa}(t)|aa|+\rho _{gg}(t)|gg|,`$ (51)
where $`\rho _{ii}(t)`$ $`(i=e,s,a,g)`$ are the time dependent populations of the collective atomic states. According to the Peres-Horodecki criterion for entanglement, the two-atom system represented by the density matrix (51) is entangled when

$`|\rho _{ss}(t)\rho _{aa}(t)|>2\sqrt{\rho _{ee}(t)\rho _{gg}(t)}.`$ (52)
Thus, in contrast to the case of only one atom excited, it is not enough to produce an unbalanced population distribution between the states $`|s`$ and $`|a`$ to obtain an entanglement in the system. The reason is the presence of a population in the unentangled state $`|e`$. In order to analyse the dependence of the entanglement on the populations of the collective states, we solve Eq. (23) with the initial condition $`\rho _{ee}(0)=1`$ and find that the time-dependent populations are given by

$`\rho _{ee}(t)`$ $`=`$ $`\mathrm{exp}(2\mathrm{\Gamma }t),`$ (53)
$`\rho _{ss}(t)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Gamma }+\mathrm{\Gamma }_{12}}{\mathrm{\Gamma }\mathrm{\Gamma }_{12}}}\left\{\mathrm{exp}\left[\left(\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)t\right]\mathrm{exp}\left(2\mathrm{\Gamma }t\right)\right\},`$ (54)
$`\rho _{aa}(t)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Gamma }\mathrm{\Gamma }_{12}}{\mathrm{\Gamma }+\mathrm{\Gamma }_{12}}}\left\{\mathrm{exp}\left[\left(\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)t\right]\mathrm{exp}\left(2\mathrm{\Gamma }t\right)\right\}.`$ (55)
One can see from Eq. (55) that in general the populations decay with different decay rates. However, for small interatomic separations, $`\mathrm{\Gamma }_{12}\mathrm{\Gamma }`$, and then the upper state $`|e`$ and the symmetric state $`|s`$ decay with the same rate $`(2\mathrm{\Gamma })`$, whereas the antisymmetric state decays with a significantly reduced rate $`\mathrm{\Gamma }\mathrm{\Gamma }_{12}`$. At early times the population is mostly in the state $`|e`$, and then the inequality (52) is not satisfied. The inequality (52) is not satisfied until $`\rho _{ee}(t)0`$. Since the population of the symmetric state decays with the same rate as $`\rho _{ee}(t)`$, at time where $`\rho _{ee}(t)0`$ the population $`\rho _{ss}(t)0`$. Thus, the lack of the entanglement for $`\rho _{ss}(t)0`$ can be attributed to a large population of the product state $`|e`$ which, in turn, decays on the same time scale as $`\rho _{ss}(t)`$.

At time $`t_s`$ when $`\rho _{ss}(t)=0`$, the upper state $`|e`$ is also depopulated, but there is still some population accumulated in the antisymmetric state, as a result of a slow decay of the population with the reduced rate $`(\mathrm{\Gamma }\mathrm{\Gamma }_{12})`$. Therefore, for times larger than $`t_s`$, the system behaves effectively as a two-level system, whose the density matrix can be written as

$`\rho (t)=\rho _{aa}(t)|aa|+\rho _{gg}(t)|gg|.`$ (56)
Following the Peres-Horodecki measure for entanglement, we see that the atoms are entangled until the antisymmetric state is populated, i.e. when $`\rho _{aa}(t)0`$. The entanglement persists for a longer time and slowly decays to zero on a time scale of order $`(\mathrm{\Gamma }\mathrm{\Gamma }_{12})^1`$, which for small interatomic separations is much longer than the single atom decay rate $`\mathrm{\Gamma }^1`$.

### B Non-identical atoms

Although most of the work on entanglement is concerned with identical atoms, there can be interesting information available on entanglement creation with non-identical atoms. In this section, we give illustrative examples of both the entanglement creation and entanglement transfer between two nonidentical atoms with different transition frequencies.

Figure 3 shows the concurrence $`C`$ and the populations $`\rho _{ss}(t)`$ and $`\rho _{aa}(t)`$ as a function of time for the initial condition of only one atom excited. We take $`\mathrm{\Delta }=\mathrm{\Gamma }`$ that the atom $`\mathrm{"}2\mathrm{"}`$ has a higher frequency (energy) than the atom $`\mathrm{"}1\mathrm{"}`$, i.e. $`\omega _2>\omega _1`$, and assume that initially the atom of the higher transition frequency was in its excited state while the other atom of lower transition frequency was in the ground state. As in the case of identical atoms, there is no initial entanglement between the atoms, and at early times the entanglement builds up rapidly to a maximum appearing at short time $`\mathrm{\Gamma }t<1`$. However, comparing with the entanglement for identical atoms, shown in Fig. 1, we see that the maximum of the entanglement obtained with nonidentical atoms is greater than that obtained with identical atoms. It is interesting to note that after passing through the maximum, the entanglement decays with two different time scales. This effect is more pronounced if we choose the atom of the lower frequency to be initially in its excited state and the other atom of higher frequency in the ground state. We illustrate this in Fig. 4, where we plot the concurrence $`C`$ and the populations $`\rho _{ss}(t)`$ and $`\rho _{aa}(t)`$ as a function of time for the same parameters as in Fig. 3. The two decay time scales are well resolved. At early times the entanglement decays on a time scale of order $`\mathrm{\Gamma }^1`$, and next the time decay ”jumps” into the time scale of order $`(\mathrm{\Gamma }\mathrm{\Gamma }_{12})^1`$, which coincides with the decay of the population of the antisymmetric state.

The enhanced entanglement and the two different decay time scales can be explained as a consequence of the coherent transfer of the population from the symmetric to the antisymmetric state, that is absent for identical atoms. In order to show this, consider the dynamics of the entangled states $`|s`$ and $`|a`$, described by Eq. (23), that in the case of initially only one atom excited simplify to the following equations of motion

$`\dot{\rho }_{ss}`$ $`=`$ $`\left(\mathrm{\Gamma }+\mathrm{\Gamma }_{12}\right)\rho _{ss}+i\mathrm{\Delta }\left(\rho _{as}\rho _{sa}\right),`$ (57)
$`\dot{\rho }_{aa}`$ $`=`$ $`\left(\mathrm{\Gamma }\mathrm{\Gamma }_{12}\right)\rho _{aa}i\mathrm{\Delta }\left(\rho _{as}\rho _{sa}\right),`$ (58)
$`\dot{\rho }_{as}`$ $`=`$ $`\left(\mathrm{\Gamma }+2i\mathrm{\Omega }_{12}\right)\rho _{as}+i\mathrm{\Delta }\left(\rho _{ss}\rho _{aa}\right).`$ (59)
Equations (51) are formally identical to the optical Bloch equations for a two-level atom driven by a coherent field . Here, the detuning $`\mathrm{\Delta }`$ plays the role of a Rabi frequency that coherently transfers the population between the symmetric and antisymmetric states. Note that the interaction between the states does not involve the ground state $`|g`$, and therefore is not accompanied by spontaneous emission. Thus, the coherent transfer of the population between the states is a decoherence free process. In this process, the population is efficiently transferred from the more populated symmetric state to the antisymmetric state before it decays to the ground state leading to the enhancement of the entanglement.

Finally, we point out one more difference between the entanglement creation in identical and nonidentical atoms. Figure 5 shows the concurrence $`C`$ and the populations $`\rho _{ss}(t)`$ and $`\rho _{aa}(t)`$ for the same parameters as in Fig. 3, but with the new initial condition $`\rho _{ss}(0)=1`$. In this case there is perfect $`(C=1)`$ initial entanglement between the atoms. Since, at $`t=0`$, only the state $`|s`$ is occupied, one finds that at early times the initial entanglement decays with the enhanced rate $`(\mathrm{\Gamma }+\mathrm{\Gamma }_{12})`$ corresponding to the decay rate of the population of the symmetric state. The entanglement decays in time until the atoms become disentangled. This happens at time $`\mathrm{\Gamma }t2`$, where the population is equally distributed over the entangled states, $`\rho _{ss}(t)=\rho _{aa}(t)`$. However, as time develops, the entanglement emerges again. This is remarkable as spontaneous emission is essentially an irreversible process, and one might expect that spontaneous emission merely degrades the initial entanglement. The revival effect is absent for identical atoms, and is due to the coherent transfer of the population from the symmetric to the antisymmetric state. The revival of the entanglement at time $`\mathrm{\Gamma }t2`$ is more easily understood by reference to the density matrix of the system. It is seen from Fig. 5 that for times $`\mathrm{\Gamma }t>2`$ the population of the symmetric state is negligible, and therefore the system behaves as a two-level system whose the density matrix is in the diagonal form $`\rho (t)=\rho _{aa}(t)|aa|+\rho _{gg}(t)|gg|`$. Since $`\rho _{aa}(t)0`$ for $`\mathrm{\Gamma }t>2`$, the atoms are entangled until the population $`\rho _{aa}(t)`$ decays eventually to the ground state.

In contrast, if the atoms were initially prepared in the antisymmetric state, $`\rho _{aa}(0)=1`$, the initial entanglement remains in the antisymmetric state for all times $`t`$ even if a part of the population is transferred to the symmetric state. This feature is shown in Fig. 6, where we plot the time evolution of the concurrence and the population $`\rho _{aa}(t)`$ for the same parameters as in Fig. 5, but with the initial condition $`\rho _{aa}(0)=1`$. It is clear from Fig. 6 that the entanglement decays with the reduced rate $`(\mathrm{\Gamma }\mathrm{\Gamma }_{12})`$ corresponding to the decay rate of the population of the antisymmetric state, and the coherent coupling does not transfer much of the entanglement to the symmetric state.

Our calculations clearly demonstrate that the transient entanglement induced from initially unentangled atoms depends crucially on the presence of the antisymmetric state, that is characteristic of spatially extended atomic systems.

## IV Summary

In this paper, we have analysed the role of the antisymmetric state in entanglement creation by spontaneous emission from two spatially separated atoms. The results show that the generation of the entanglement is strongly dependent upon the population distribution between the symmetric and antisymmetric states of the system. The entanglement is maximal when the population of the symmetric state becomes zero. Thus, our results show the participation of the symmetric state in the atomic dynamics has a destructive effect on the entanglement. We have also considered the entanglement creation in two non-identical atoms. We have found that the entanglement can be enhanced by the process of the coherent transfer of the population between the symmetric and antisymmetric states. In addition, we have shown that the entanglement can decay with two different time scales and exhibits a revival behavior due to the coherent transfer of the population from the rapidly decaying symmetric to the slowly decaying antisymmetric state before it decays to zero. Although the entanglement created by spontaneous emission appears only in the transient regime, it provides an useful information about entanglement sharing and entanglement transfer between two entangled states.