# Phonon-assisted decoherence and tunneling in quantum dot molecules

###### Abstract

\abstcolWe study the influence of the phonon environment on the electron dynamics in a doped quantum dot molecule. A non-perturbative quantum kinetic theory based on correlation expansion is used in order to describe both diagonal and off-diagonal electron-phonon couplings representing real and virtual processes with relevant acoustic phonons. We show that the relaxation is dominated by phonon-assisted electron tunneling between constituent quantum dots and occurs on a picosecond time scale. The dependence of the time evolution of the quantum dot occupation probabilities on the energy mismatch between the quantum dots is studied in detail.

e-mail anna.grodecka-grad@nbi.dk, Phone: +45-353-25426, Fax: +45-353-25400

XXXX

## 1 Introduction

Quantum dots (QDs) are one of the promising candidates for building a feasible quantum computer. In particular, it has been proposed to use two coupled quantum dots called quantum dot molecules (QDMs) for various schemes of quantum computation, where the two ground states of a single confined electron as well as singlet and triplet states of doubly doped structures can be used as the logical qubit states. The electrical [1] and optical control [2, 3] of spin- and charge-based qubits have already been demonstrated. However, in such solid state systems, the phonon-assisted relaxation [5, 6] can strongly affect the coherent control. The energy difference between the two electron ground states, which is typically of order of a few meV, can lead to pure dephasing processes [7, 8, 9] and to electron tunneling [10, 11, 12].

In this paper, we present the full quantum kinetic description of the phonon-mediated relaxation in doped quantum dot molecules including non-Markovian effects. We show that the coupling to the phonon reservoir in quantum dot molecules can lead to a fast electron tunneling on a picosecond timescale, which strongly affects the coherent electron evolution. We employ the non-Markovian correlation expansion technique [13, 14] including up to three-particle correlations. Due to space constrains, we present here equations of motions with up to two-particle correlations. We include diagonal and off-diagonal couplings to the phonon reservoir representing virtual and real phonon-assisted processes, respectively. The dependence of the time evolution of the quantum dot occupation probabilities on the energy mismatch between the constituent quantum dots is studied in detail.

## 2 Model system

We consider a system consisting of a single quantum dot molecule doped with one electron. The relaxation between the two energetically lowest states of the electron in the left and in the right quantum dot is considered (see Fig. 1 with schematic plot of the energy levels in the quantum dot molecule). The free Hamiltonian of the electron reads

(1) |

where is the direct tunneling coupling element between the quantum dots and is the energy difference between the ground states in both quantum dots.

The Hamiltonian describing the free evolution of phonons is

(2) |

with and being phonon creation and annihilation operators, respectively. denotes phonon frequencies with phonon wave vector , where different phonon branches are labeled by (longitudinal and two transverse and ).

The interaction between the electron and phonon reservoir reads

(3) |

The relevant coupling elements for the acoustic phonons coupled via both piezoelectric coupling and deformation potential read:

(4) |

(5) |

with Gaussian electron wave functions

(6) |

We have numerically simulated the self-assembled GaAs quantum dots with the following parameters: nm, nm, and nm, where is the electron wave function size of the -th quantum dot, and the distance between the dots is . Here, kg/m is the crystal density, is the normalization volume of the phonon modes, is the speed of sound (longitudinal m/s or transverse m/s, depending on the phonon branch), C/m is the piezoelectric constant, eV is the deformation potential constant for the electrons and is the static dielectric constant. The function in the piezoelectric coupling element reads:

(7) | |||||

(8) | |||||

(9) |

We study the time evolution of the electron confined in a quantum dot molecule in the presence of the coupling to the phonon reservoir within the density matrix theory. We employ the second order correlation expansion method [13, 14], where it is assumed that correlations involving an increasing number of particles are of decreasing importance. This non-perturbative technique covers the memory effects in the non-Markovian regime.

With the help of Heisenberg equation we derive the equations of motion for the quantities of interest. The equation for the electron occupation reads:

(10) | |||||

It couples to the coherence as well as to the phonon-assisted coherences and and to phonon amplitudes . The factorization scheme has been used above with

(11) |

where the quantities have been decomposed into all possible lower-order factorizations. The occupation is directly affected by phonon reservoir in the case of the off-diagonal coupling and indirectly by the coherence , which evolves in the following way:

Here, a direct impact from phonons is present not only via real transitions but also for the pure dephasing described by the diagonal coupling elements and . If these coupling elements are equal, the two pure dephasing channels cancel each other out . In this case, the information cannot leak out, since the phonons is not delivering any which way information and no decoherence channel exists.

The evolution of phonons reveals the diagonal and off-diagonal couplings to the density and coherence, respectively,

The next step is to derive the equations of motion for all the phonon-assisted correlation quantities.

The dynamics of the phonon-assisted density, is described by

All these phonon-assisted two-particle correlations couple to the phonon-phonon correlations, phonon density and phonon coherence :

In addition, the two-particle correlations couple to three particle correlations, e.g. . and these couple to up to four particle correlations, etc. Thus, in order to get a closed set of equations, one needs to truncate the hierarchy by neglecting higher order correlations. In the present paper, up to three particle correlations were included.

## 3 Phonon-mediated tunneling

Initially, one electron is injected into the quantum dot with higher energy, . The time of the injection is assumed to be much shorter than the response of the phonon reservoir, thus at the initial time we can set all correlations to zero.

The time evolution of the occupation probability of the quantum dot with higher energy in the presence of the coupling to the phonon reservoir is shown in Fig. 2 for different values of the energy difference between the two electron ground states. The case of the ideal evolution, where the electron-phonon coupling is absent, is shown in the inset of Fig. 2. The electron tunnels coherently between the two quantum dots with a period of these oscillations that is determined by the energy difference between the ground states and the tunneling coupling and scales like . The maximum occupation also depends on these parameters. If the energy difference is twice the value of the tunneling coupling, the electron tunnels between the quantum dots with maximum occupation of of the quantum dot with lower energy. For the larger energy differences, the tunneling coupling is too weak to drive the electron into the adjacent quantum dot and the electron stays in the initial quantum dot.

The dynamics strongly changes if the coupling between the electron and the phonon reservoir is present (see Fig. 2), which results in the electron tunneling to the neighboring quantum dot on a short picosecond time scale. This phonon-assisted tunneling is only properly modeled if the off-diagonal electron-phonon coupling term in the Hamiltonian is included. Both coupling via deformation potential and piezoelectric effect are important and contribute to the relaxation. In the case of the energy difference of meV, the thermalization process takes about ps and finishes with a delocalized state of the electron with a probability of for being in the left quantum dot and in the second dot. For larger energy mismatches of and meV, the relaxation takes about and picoseconds, respectively, and in the final state, the electron is localized in the quantum dot with lower energy. In these cases, the probability of phonon emission or/and absorption processes is the highest since the energies of the relevant acoustic phonons lie in this parameter regime. If we increase the energy difference between the dots, the probability of phonon-mediated relaxation will decrease and the tunneling is slower and less efficient. The presented results were calculated for the low temperature of K. The phonon-assisted tunneling will get stronger and faster at higher temperatures and will also affect the final state, since the ratio of probabilities of phonon emission and absorption changes with temperature.

## 4 Conclusion

A full description of the electron dynamics in the presence of the electron-phonon coupling in a quantum dot molecule doped with a single electron has been presented. The two couplings to the relevant acoustic phonons via deformation potential and piezoelectric coupling with diagonal and off-diagonal interactions have been taken into account. It has been shown that the phonon-mediated relaxation is a fast process on a picosecond timescale strongly modifying the coherent evolution of the electron. It is dominated by off-diagonal electron-phonon coupling. We analyzed the dependence of the relaxation on the energy difference between the two quantum dots and indicated the vales of parameters when the phonon-mediated tunneling is most efficient.

The authors acknowledge support from the Emmy Noether Program of the DFG (Grant No. FO 637/ 1-1) and the DFG Research Training Group GRK 1464.

## References

- [1] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, Science 309, 2180 (2005).
- [2] M. Bayer, P. Hawrylak, K. Hinzer, S. Fafard, M. Korkusinski, Z. R. Wasilewski, O. Stern, and A. Forchel, Science 291, 451 (2001).
- [3] H. E. Türeci, J. M. Taylor, and A. Imamoglu, Phys. Rev. B 75, 235313 (2007).
- [4] D. S. Golubev and A. D. Zaikin, Phys. Rev. Lett. 81, 1074 (1998).
- [5] T. Brandes and B. Kramer, Phys. Rev. Lett. 83, 3021 (1999).
- [6] S. Debald, T. Brandes, and B. Kramer, Phys. Rev. B 66, 041301 (2002).
- [7] A. Grodecka, C. Weber, P. Machnikowski, and A. Knorr, Phys. Rev. B 76, 205305 (2007).
- [8] P. Machnikowski and L. Jacak, Phys. Rev. B 69, 193302 (2004).
- [9] A. Grodecka, P. Machnikowski, and J. Förstner, Phys. Rev. A 79, 042331 (2009).
- [10] A. Grodecka, P. Machnikowski, and J. Förstner, Phys. Rev. B 78, 085302 (2008).
- [11] A. Grodecka-Grad and J. Förstner, Phys. Rev. B 81, 115395 (2010).
- [12] K. Gawarecki, M. Pochwala, A. Grodecka-Grad, and P. Machnikowski, Phys. Rev. B 81, 245312 (2010).
- [13] F. Rossi and T. Kuhn, Rev. Mod. Phys. 74, 895 (2002).
- [14] J. Förstner, C. Weber, J. Danckwerts, and A. Knorr, Phys. Rev. Lett. 91, 127401 (2003).