# Entanglement Dynamics of Two Independent Cavity-Embedded Quantum Dots

## Abstract

We investigate the dynamical behavior of entanglement in a system made by two solid-state emitters, as two quantum dots, embedded in two separated micro-cavities. In these solid-state systems, in addition to the coupling with the cavity mode, the emitter is coupled to a continuum of leaky modes providing additional losses and it is also subject to a phonon-induced pure dephasing mechanism. We model this physical configuration as a multipartite system composed by two independent parts each containing a qubit embedded in a single-mode cavity, exposed to cavity losses, spontaneous emission and pure dephasing. We study the time evolution of entanglement of this multipartite open system finally applying this theoretical framework to the case of currently available solid-state quantum dots in micro-cavities.

###### pacs:

03.67.Mn, 03.65.Yz, 03.67.-a## 1 Introduction

Cavity quantum electrodynamics (CQED) deals with the interaction among photons confined in a reflective micro-cavity and atoms or other particles. When a two-level system (for example, a two-level atom) is strongly coupled to a cavity mode [1], it is possible to realize important quantum information processing tasks, such as controlled coherent coupling and entanglement of distinguishable quantum systems [2, 3, 4]. In this respect solid-state devices, and in particular semiconductor quantum dots (QDs), utilized as ”artificial atoms”, are one of the most promising architectures for the possibility of miniaturization, electrical injection, control and scalability. Recently, thanks to impressive progress in the technology of solid-state microcavities, substantial advances have been made towards these goals. The strong coupling regime has been reached for the excitonic transition of quantum dots (QDs) [5, 6, 7, 8, 9], and nanocrystals [10]coupled to optical semiconductor cavities, as well as for superconducting qubits coupled to microwave cavities [11]. In all of these systems, the cavity-mode quality factor (Q) can be very large while solid-state emitters are intrinsically coupled to the matrix they are embedded in. In fact, decoherence and phase relaxation unavoidably broaden any transition between the discrete states of such artificial atoms. High experimental performances are required to realize quantum processors and it is thus important to establish how long a sufficient degree of entanglement can be maintained in spite of losses, decoherence and noise. Implications are the possibility to store entangled states in solid-state memories and entanglement preservation during local operations in quantum algorithms [12]. Understanding how the entanglement is transferred from, e.g., a pair of independent initially entangled qubits to reservoirs has motivated several contributions in recent years [13, 14]. The aspect that has mostly drawn attention is the possibility of a complete disappearance of the entanglement between the qubits at finite times [15]. The occurrence of this phenomenon, termed entanglement sudden death (ESD) [16] and of entanglement revivals [17] have been shown in quantum optics experiments [18, 19]. Entanglement transfer from atoms to cavity modes leading to entanglement revivals has also been studied [20]. The effects on the dynamics of entanglement of cavity losses, spontaneous emission and pure dephasing for two-qubit systems have been extensively investigated [21, 22]. The general problem of the dynamics of entanglement in the simultaneous presence of more than one noise has been also studied finding that, for composite systems, the additivity of decay rates of single noises is not maintained [23]. The main aim of this paper is to investigate the dynamical behavior of entanglement for a realistic implementation for quantum computation of two qubits in separated cavities where several sources of noise are present.

Here we consider the entanglement dynamics in a system composed by two initially correlated solid-state emitters each strongly coupled to a lossy cavity interacting with its reservoir. We include the unavoidable losses due to spontaneous emission into external electromagnetic modes distinct from the lossy single-mode of the cavity. These, always present, additional losses arise from the coupling of the emitter to a continuum of leaky modes [24]. We also include pure dephasing noise which plays a relevant role in solid state quantum emitters. For example, a QD interacts with the phonons of the matrix it is embedded in, giving rise to sidebands in addition to the so-called zero-phonon line (ZPL). At sufficiently low temperature yet, the emission in the ZPL remains predominant, allowing to model these systems as effective two-level systems subject to additional pure dephasing [24]. We also study the entanglement transfer from the two-emitter system to the cavity modes. The input-output relations for optical cavities [25] show that the entanglement between cavity modes can in principle be measured experimentally by collecting photons escaping the cavities. Such two-cavity entanglement dynamics could thus be exploited to monitor the entanglement dynamics of cavity-embedded solid state emitters.

The paper is structured as follows. In Sec. 2 we present and solve a model for the physical configuration described above. In Sec. 3 we explore the dynamical behavior of entanglement for quite general values of the physical parameters of the system. In Sec. 4 we specialize to the case of currently available solid-state quantum dots in micro-cavities giving the characteristic lifetimes of entanglement. In Sec. 5 we summarize our results.

## 2 Model

Our system is composed by two noninteracting subsystems (), each consisting in a qubit (two-level emitter) coupled to a single-mode cavity in turn interacting with an external reservoir (see Fig. 1). The Hamiltonian of the total system is thus given by the sum of the Hamiltonians of the two noninteracting subsystems

(1) |

In each subsystem, we distinguish the bipartite system made by the qubit and the cavity from the reservoir . The Hamiltonian of each part reads like (we omit index )

(2) |

where describes, in the rotating wave approximation (RWA), the qubit-cavity system, is the coupling constant between qubit and cavity, denotes the usual diagonal Pauli matrix, are the two-level raising and lowering operators, and the annihilation and creation operators for the cavity mode. We observe that this Hamiltonian model is valid for values of smaller than [26]. describes the external environment responsible of the different noise sources which affect the qubit-cavity system and the interaction of the latter with the environment. The realistic conditions present in a system composed by quantum dots embedded in microcavities are modeled as three noise sources: cavity losses, qubit spontaneous emission and pure dephasing mechanisms.

In the usual master equation approach, considering the Markov approximation and an infinite number of bath oscillators, we can describe the dynamics of the qubit-cavity system by

(3) |

The Markovian processes are described by the Liouvillian term, that consist of 3 parts:

(4) |

describes the cavity losses of photon in the reservoir modes and is expressed by the following form [12]:

(5) |

being the photon escape rate from the cavity to free space.
In addition, the qubit is subjected to decay *via* spontaneous emission and losses of coherence.
describes spontaneous emission in the leaky modes: all the available light modes except the cavity one [12],

(6) |

and are based on Hamiltonians in RWA and thus are not valid for arbitrarily large values of the decay rates and [26]. Finally,

(7) |

describes pure dephasing processes [12]. Although the three noises are treated in the Markovian limit, the reduced dynamics of the qubit-cavity system can present non-Markovian behavior depending on the strength of the coupling constant with respect to the various decay rates. We shall comment quantitatively on this point after Eq. (2.1.1).

The optical cavity is an open quantum system, cavity photons can escape it and propagate into free space on along an optical fibre until they eventually reach another distant quantum system or can be detected. The quasimode approach is able to describe in a direct way such a behavior. A relationship between the external fields and the intracavity field may be obtained [25, 27] in the limit of continuous spectrum. If the coupling constant between cavity and external bosonic modes is independent of frequency over a band of frequencies about the characteristic frequency , :

(8) |

where the operators and are related to the reservoir operators and . Throughout this paper we will consider the case of no input photons, hence, once known the quantum state for the cavity mode, it will be possible to calculate expectation values and correlation functions for output photons that can be measured experimentally or be used as input to transmit the entanglement to distant quantum systems.

### 2.1 Procedure

Being the two subsystems noninteracting, they evolve independently so that we can analyze the dynamics of only one subsystem and in turn obtaining the evolution of the total four-partite open system [17]. We will consider initial states with zero cavity photons and at most one excitation in each qubit. Eq. (3) allows us to compute the joint evolution of the total four-partite open system starting from an arbitrary initial state where only one excitation is present in each atom. From the knowledge of the evolved density matrix, it will be possible to investigate the various reduced dynamics of the total system, for example that of the two qubits or of the two cavities. In the following we start showing how to compute for each part the time dependent density matrix elements, which in general may be different for the two subsystems.

#### Dynamics of subsystems

Here we consider the dynamics of a single subsystem . We choose the standard product basis , where () and () are the lower and upper state of the qubit (cavity). The dynamics of qubit under the effect of the master equation of Eq. (3) is described by the reduced density matrix

(9) |

The time dependent coefficients and can be obtained analytically, however the presence of pure dephasing gives rise to very cumbersome and lengthy equations. In this section we present analytical results only for the case . Numerical results in presence of pure dephasing shall be included in next section. Analogously, for the cavity modes, the dynamics of the reduced density matrix can be expressed as

(10) |

where a zero-photon initial state ( and ) has been considered. The dynamics of the reduced density matrices is obtained using the master equation expressed by Eq. (3), that is appropriate when the reservoir is at zero temperature, the coupling between the cavity and the external modes of the reservoir has a flat spectrum in the range of involved frequencies and the qubit is resonant with the cavity [28, 29]. We will limit our investigation to this physical condition. In the absence of pure dephasing, from Eq. (3) one obtains,

(11) |

with and and where we have introduced the characteristic frequency . We point out that the function is not directly connected to the decay of the cavity excited state but is linked to the exchange of the initial qubit excitation between the qubit itself, the cavity and the reservoir. In this sense, the fact that goes to zero when reflects its dependence on the initial conditions and on the internal qubit-cavity dynamics. Non-Markovian features in the dynamics of the subsystems occur for values of large enough to make real, but such as not to compromise the validity of the RWA in our model.

## 3 Entanglement evolution

After obtaining in the previous section all the relevant dynamical coefficients, in this section we give the explicit expressions of concurrence for some couples of parties of the four-partite system. In particular we consider the two-qubit and two-cavity entanglement separately for two different initial configurations. The qubits are initially in one and two-excitation Bell-like states, while cavities are in their vacuum state.

We shall restrict our analysis of entanglement dynamics to the two-excitation entangled initial states

(12) |

where in each ket the first entry denotes a ( or ) state of subsystem , while the second entry a state of subsystem . Generalization of the results for other (eventually mixed) initial states is straightforward.

From the evolved state one finds the reduced density matrices of the bipartite system of interest tracing over the degrees of freedom of the noninvolved parties. We represent the density matrices in the standard computational basis . In this way, the two-qubit state at time is, e.g., given by

(13) |

In the following we will consider the case of two identical subsystems, hence . An analogous result can be obtained for the two-cavity system density matrix just replacing and with and respectively.

The concurrence corresponding to the density matrix of Eq. 13 in the case of identical subsystems is found to be [17]

(14) |

We present numerical calculations for a specific two-excitation entangled initial state with . Figure 2 displays the concurrence dynamics of the qubits and of the cavity modes in absence (2a) and presence (2b) of pure dephasing. Figure 2 also displays , providing information on the detectable output photon flux . In both cases , the other decay rates being fixed as and in Fig. 2a while , in Fig. 2b. For these values, non-Markovian features appear in each qubit-cavity subsystem dynamics, leading to similar effects in the dynamics of the plotted quantities. The figure shows that pure dephasing affects heavily the entanglement dynamics increasing the entanglement decay and enabling its sudden death. Although phase noise acts directly only on the emitter, owing to the strong coupling between the emitter and the cavity, it affects dramatically also the two-cavity concurrence.

Figure 3a shows the two-emitter concurrence as function of time and of the amount of phase noise . Figure 3b displays the two-cavity concurrence . We used , , . The detrimental effect of phase noise on both the dynamics of the entanglement of emitters and cavities is evident.

## 4 Application to quantum dots under realistic conditions

Here we specialize to the case where the considered system is implemented by using currently available quantum dots as quantum emitters, embedded in separated micro-cavities [9]. Solid state microcavities with three-dimensional photon confinement, high Q, and small volume mode can be realized by fabricating a photonic crystal slab structure with a nanocavity composed of one or more missing air holes. The slab incorporates a central layer of low density self-assembled InAs quantum dots [5, 9]. Another geometry of particular interest is that of micropillar cavities [6, 30]. In these systems the fundamental cavity mode can be coupled to and from the outside with a very high coupling efficiency. Moreover, they offer interesting perspectives for the implementation of quantum information protocols using charged quantum dots [31]. Figure 4a displays the concurrence dynamics of the quantum emitters and of the cavity modes obtained for two independent cavity-embedded quantum dots (). We consider typical system parameters for the state-of-art microstructures [9]: eV, eV,eV, eV. The figure shows that the two-dot entanglement after a rebirth survives up to about ps. Quantum operations based on all-optical implementations can be performed by means of ultrafast pulses. At optical frequencies, pulses of 20-100 fs time-width are currently available. Very recently structures displaying higher Q values but with a quite low coupling have been realized [30]. Figure 4b displays the concurrence dynamics of the quantum emitters and of the cavity modes obtained by using and parameters describing this novel micropillar structure [30]: eV, eV,eV, eV. In this case both the two-dot and two-cavity modes entanglement increase their lifetime. Figure 4 also displays the behavior of which is proportional to the detectable output photon flux .

We notice that in the systems considered in this section, dots and cavities frequencies are eV so that the validity limits of RWA are well satisfied.

## 5 Conclusions

We have studied the dynamics of entanglement in a system made by two solid-state emitters, as two quantum dots, embedded in two separated micro-cavities. In addition to the coupling with cavity mode, the emitter is subject to spontaneous emission, due to the coupling with a continuum of leaky modes, and to phonon-induced pure dephasing mechanisms. We have modeled this physical system as a multipartite system composed by two independent parts each containing a qubit exposed to cavity losses, spontaneous emission and pure dephasing. The numerical results presented here for arbitrary values of the physical parameters put forward the impact of pure dephasing on the entanglement dynamics of the quantum emitters and of the two-cavity modes. Experimental information about the latter can be gathered by detection of the collected cavity output field. We have finally applied this theoretical framework to the specific case of currently available solid-state quantum dots in micro-cavities.

## References

### References

- Scully M O and Zubairy M S 1997 Quantum Optics (Cambridge University Press)
- Imamoglou A et al. 1999 Phys. Rev. Lett. 83 4204
- Monroe C 2002 Nature 416 238
- Michler P et al. 2000 Science 290 2282
- Reithmaier J P et al. 2004 Nature 432 197
- Yoshie T et al. 2004 Nature 432 200
- Peter E et al. 2005 Phys. Rev. Lett. 95 067401
- Press D, Gotzinger S, Reitzenstein S, Hofmann C, A Loffler M Kamp A F and Yamamoto Y 2007 Phys. Rev. Lett. 98 117402
- Hennessy K et al. 2008 Nature 445 896
- Thomas N L et al. 2006 Nano Lett. 6 557
- Wallraff A et al. 2004 Nature 431 162
- Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge University Press)
- López C E, Romero G, Lastra F, Solano E and Retamal J C 2008 Phys. Rev. Lett. 101 080503
- Bai Y K, Ye M Y and Wang Z D 2009 Phys. Rev. A 80 044301
- Yu T and Eberly J H 2009 Science 323 598
- Yu T and Eberly J H 2004 Phys. Rev. Lett. 93 140404
- Bellomo B, Lo Franco R and Compagno G 2007 Phys. Rev. Lett. 99 160502
- Almeida M P et al. 2007 Science 316 579–582
- Xu J S et al. 2010 Phys. Rev. Lett. 104 100502
- López C E, Romero G and Retamal J C 2010 Phys. Rev. A 81 062114
- Huang J H and Zhu S Y 2007 Phys. Rev. Lett. 76 062322
- Bellomo B, Lo Franco R and Compagno G 2008 Phys. Rev. A 77 032342
- Yu T and Eberly J H 2006 Phys. Rev. Lett. 97 140403
- Aufféves A, Gérard J M and Poizat J P 2007 Phys. Rev. A. 79 053838
- Walls D F and Milburn G J 2008 Quantum Optics (Berlin: Springer-Verlag)
- Werlang T, Dodonov A V, Duzzioni E I and Villas-Bôas C J 2008 Phys. Rev. A 78 053805
- Gardiner C W and Collett M J 1985 Phys. Rev. A 31 3761
- Scala M, Militello B, Messina A, Piilo J and Maniscalco S 2007 Phys. Rev. A 75 013811
- Scala M and al 2007 J. Phys. A 40 14527
- Loo V et al. 2010 arXiv: 1011.1155v1 [cond–mat.mes–hall]
- Hu C J et al. 2009 Phys. Rev. B 80 205326