Spontaneous Emission Spectrum in Double Quantum Dot Devices

See allHide authors and affiliations

Science  30 Oct 1998:
Vol. 282, Issue 5390, pp. 932-935
DOI: 10.1126/science.282.5390.932


A double quantum dot device is a tunable two-level system for electronic energy states. A dc electron current was used to directly measure the rates for elastic and inelastic transitions between the two levels. For inelastic transitions, energy is exchanged with bosonic degrees of freedom in the environment. The inelastic transition rates are well described by the Einstein coefficients, relating absorption with stimulated and spontaneous emission. The most effectively coupled bosons in the specific environment of the semiconductor device used here were acoustic phonons. The experiments demonstrate the importance of vacuum fluctuations in the environment for quantum dot devices and potential design constraints for their use for preparing long-lived quantum states.

Electronic quantum devices allow the quantum mechanical properties of electrons confined to small regions in a solid to be explored. Existing devices include semiconductor resonant tunneling diodes (1) (based on quantum mechanical confinement), superconducting Josephson junction circuits (2) (based on macroscopic phase coherence), metallic single-electron transistors (3) (based on quantization of charge), and molecular electronic devices (4). The principle of operation in circuits of these devices is based on controlling energy states, for instance, by means of an external (gate) voltage. Thermal energy is always a source for unwanted transitions and errors. Even at zero temperature, however, vacuum fluctuations in the environment can give rise to transitions between states of nonequal energy by spontaneous emission of an energy quantum. Such inelastic transitions cause errors in many proposed schemes for quantum circuits. Here we studied inelastic transitions in a fully controllable, two-level quantum system realized in a double quantum dot device. We can relate the transition rates involving emission to absorption rates by the Einstein coefficients over the full energy and temperature range we studied. At the lowest temperature (23 mK), we directly measured the energy-dependent rate for spontaneous emission and determined that, in our specific semiconductor device, this energy is emitted into the environment formed by acoustic phonons.

Our double quantum dot (Fig. 1A) is fabricated in the two-dimensional electron gas (2DEG) of an AlGaAs-GaAs semiconductor heterostructure (5). The source and drain are large 2DEG regions that serve as leads for contacting current and voltage wires. The two dots, L and R, are separated from each other and from the leads by potential barriers induced by negative voltages applied to the three metallic gates. Tunneling between the different regions is sufficiently strong to detect current but weak enough that the number of electrons in each dot is a well-defined integer. The energy states in such fully confined regions are discrete, 0D states, resembling discrete atomic states (6, 7). The discrete energies include contributions from single-electron charging energies (arising from Coulomb interactions) and from quantum-mechanical confinement. The lowest energy state for one additional electron in the L dot is labeled in Fig. 1, B to D, asE L, and similarlyE R for the R dot. Figure 1C illustrates the resonance condition, E L =E R, in which case an electron can tunnel elastically from an occupied state in the source viaE L and E R to an empty state in the drain. Such tunneling sequences of single electrons are regulated by the Coulomb charging energies (3, 7). When the two states are not aligned, E LE R, only inelastic transitions are allowed for which some energy needs to be exchanged with the environment. A measured off-resonance current, therefore, directly provides information about the coupling between electrons on the dots to degrees of freedom in the environment. The inelastic rates can be analyzed with well-developed methods in quantum optics (8, 9).

Figure 1

(A) Double quantum dot device defined in the 2DEG of a AlGaAs-GaAs heterostructure by focused ion beam implantation. The narrow channel connects the large 2D source and drain leads. Negative voltages (V GL,V GC, and V GR) applied to the metal gates (G L,G C, and G R; widths are 40 nm) induce three tunable tunnel barriers in the wire. The two quantum dots, L and R, respectively, contain ∼15 and ∼25 electrons; charging energies are ∼4 and ∼1 meV; and the measured average spacing between single-particle states is ∼0.5 and ∼0.25 meV. (B, through D) Energy diagrams (vertical axis) along the spatial axis through the dots (horizontal axis) for the tunnel situations: absorption, elastic, and emission. Thick vertical lines denote tunnel barriers. The continuous electron states in the leads are filled up to the Fermi energies μS and μD. The external voltage V SDbetween leads opens a transport window of sizeeV SD = μS − μD. The energy, ɛ ≡ E LE R, is defined as the difference between the topmost filled discrete state of the L dot, E L, and the lowest discrete state for adding an extra electron to the R dot,E R. (The interdot capacitance prevents that E L and E R are simultaneously occupied.) An elastic current can flow when ɛ = 0, otherwise a nonzero current requires absorption (ɛ < 0) or emission of energy (ɛ > 0). T c is the tunnel coupling and Γi is the inelastic rate between the two dots. ΓL and ΓR are the tunnel rates across the left and the right barriers. (E) Typical measurement of the current (solid) versus ɛ at 23 mK. The measured current is decomposed in an elastic (dashed) and an inelastic (dotted-dashed) part.

A typical current spectrum versus ɛ ≡ E LE R at our lowest lattice temperatureT = 23 mK (10) is shown in Fig. 1E. The gate voltages V GR andV GL are swept simultaneously such that the respective energies are like those illustrated in Fig. 1, B to D; that is, ɛ = 0 occurs in the middle between the Fermi energies of source and drain, μS and μD, and |ɛ| =eV SD (its maximum) corresponds to having the states E L and E Raligned to one of the Fermi energies. To analyze the large asymmetry, we decomposed the total current I tot(ɛ) =I el(ɛ) + I inel(ɛ > 0) into a symmetric part I el(ɛ) =I el(−ɛ) (dashed curve) and the remaining asymmetric part I inel(ɛ > 0) (dotted-dashed curve). At T = 0, Iel(ɛ) is due to elastic tunneling and has a Lorentzian line shapeI el(ɛ) =I el,maxω2/(ω2 + ɛ2) (11). The full width at half maximum, 2ω, can be tuned by the central gate voltageV GC roughly from 4 to 20 μeV. From measurements of I el(ɛ) at positive and negative V SD, it is possible to extract values for the tunnel couplings ΓL, ΓR, and T c (11,12).

The remaining current, I inel(ɛ > 0), which is nonzero only for ɛ > 0, is due to inelastic tunneling. In Fig. 1E,I inel is nonzero over an energy range of ∼100 μeV, despite the thermal energy kT (23 mK) = 2 μeV being much smaller (k is Boltzmann's constant). (The irregular fine structure is discussed below.) In general, we found thatI inel vanishes when one of the levels,E L or E R, crosses one of the two Fermi energies. In the specific case of Fig. 1E,E L and E R cross the Fermi energies simultaneously, implying that I inel is cut off at ɛ = eV SD. Below this cutoff, the value of I inel was not influenced by the value of V SD (13). ForT = 0, we can write the condition for a nonzero inelastic current as μS > E L >E R > μD = μSeV SD. The amount of inelastic current depends on the transition rates as I inel(ɛ) =eL −1 + Γi −1(ɛ) + ΓR −1]−1. When the inelastic rate Γi(ɛ) from E L toE R is much smaller than the rates through the outer barriers, then I inel(ɛ) =eΓi(ɛ).

The effect of a nonzero temperature on the current is shown in Fig. 2A. A higher temperature Tenhances I tot on both the emission (ɛ > 0) and the absorption (ɛ < 0) side. The absorption spectrum shows an exponential temperature dependence,e ɛ/kT (dashed lines), for absolute energies larger than the elastic current measured at 23 mK, that is |ɛ| > w.

Figure 2

(A) Measured current versus ɛ forT = 23 to 300 mK. The current is measured foreV SD = 140 μeV while sweepingV GR and V GLsimultaneously in opposite directions such that we change the energy difference ɛ. Gate voltage is translated to energy ɛ by a calibration better than 10% with photon-assisted tunneling measurements (16). Dashed lines indicate exponential dependence, e ɛ/kT, for |ɛ| ≫ kT. Arrows point at steplike structure on the emission side (ɛ > 0) and a shoulder on the absorption side (ɛ < 0). From fits (11) to the elastic current part at 23 mK, we obtainedhΓRhT c ∼ 1 μeV and hΓL ∼ 5 μeV for this data set. (B) Reconstructed current, I tot(ɛ) = I el(23 mK) +I inel(ɛ,T) for differentT. The spontaneous emission spectra derived from the measured data at 23 mK and Eq. 1 are used to reconstruct the full temperature and energy dependence. (C) The absorption rateW a (open symbols) and emission rateW e (closed symbols) normalized by the spontaneous emission rate A versus kT/|ɛ|. Circles, squares, upper and lower triangles, and diamonds are taken at |ɛ| = 18, 24, 40, 60, and 80 μeV, respectively [see also symbols in (A)]. The solid line indicates the Bose-Einstein distribution, 〈n〉, whereas the dashed line shows 〈n〉 + 1.

To analyze the temperature dependence, we assumed boson statistics for the degrees of freedom in the environment. The average occupation number 〈n〉 of environmental modes at energy ɛ is given by the Bose-Einstein distribution function: 〈n〉 = 1/(e ɛ/kT − 1). The rates for absorption, W a, and emission, W e, can be expressed very generally by W a =B aρ and W e= A + B eρ, where the Einstein coefficients stand for spontaneous emission (A) and stimulated emission (B e) and absorption(B a), and ρ is the energy density (8). From the Einstein relations, B a= B e = A〈n〉/ρ (8), we obtainEmbedded Image(1) Embedded ImageTo test whether the inelastic current follows emission and absorption statistics, we calculated the full current spectrum from Eq. 1. First, we obtained the spontaneous emission rate fromA(ɛ) = I inel(ɛ > 0,T = 23 mK)/e. The trace at 23 mK is effectively at zero temperature for ɛ ≫ 2 μeV because then 〈n〉 ≪ 1. The emission current at higher temperatures follows from I inel(ɛ > 0,T) = e(〈n〉 + 1)A(ɛ), whereas the absorption current follows fromI inel(ɛ < 0, T) =enA(−ɛ). The reconstructed current spectrum is shown in Fig. 2B. The central part of the curves (|ɛ| < 10 μeV) is kept blank because Eq. 1 does not include the T dependence of I el. The calculated current reproduces the measured current well up to 200 mK. Even the small steplike feature seen at ɛ ∼ 30 μeV is reflected by a shoulderlike feature at ɛ ∼−30 μeV in the measured and in the calculated absorption spectra (indicated by arrows). ForT > 200 mK, the measured current significantly exceeds the calculated current, which is probably due to thermally excited electrons. (Equation 1 only describes the T dependence of the environment. The thermal excitations in the electron leads are not included.) Further confirmation of the applicability of the Einstein relations to our quantum dot system follows from the prediction [I inel(ɛ > 0) −I inel(ɛ < 0)]/eA(|ɛ|) = [W eW a]/A = 1, which is valid independent of T. A plot of the normalized rates,W a/A andW e/A, versuskT/|ɛ| for various ɛ and T up to 200 mK is shown in Fig. 2C. The measured data closely follow the prediction [W eW e]/A = 1; that is, the normalized rates, W a/A andW e/A, differ by 1 over the temperature range T < 200 mK without fitting any parameter.

The inelastic rate for a two-level system coupled to a bosonic environment at T = 0 is expected to have aT c 2 dependence (14,15). Still, without identifying the bosonic environment, we can test this dependence on the elastic tunnel couplingT c between the two dots. The inelastic current increases with T c, as shown inFig. 3A. For the largest coupling, we obtain a saturation where the elastic current peak can no longer be distinguished. By fitting the elastic current part to a Lorentzian line shape (11) we can obtain rough estimates forT c as long as the current is less than the saturation value. We found that with these fitted values, the inelastic current scales as T c α with an exponent α = 2.5 ∼ 3, perhaps somewhat greater than expected. The effect of the increased coupling on the symmetric part of the current at low T is shown in Fig. 3B. For small tunnel coupling, we always obtain Lorentzian line shapes. For increasing couplings, the data still fit to a Lorentzian tail on the absorption side. However, we generally found significant deviations for small ɛ, implying that for large coupling the elastic and inelastic rates can become of the same order. This effect may form a significant limitation for the coherence time in coupled quantum devices (16).

Figure 3

Current spectrum for different coupling energies at 23 mK. (A) The magnetic field is 1.6 T for (i) and 2.4 T for the other curves (10). The curves have an offset, and curve (i) is multiplied by 5. Rough estimates for the coupling energies are (i) hT c (∼0.1 μeV;h is Planck's constant) ≪ hΓR(∼10 μeV), (ii) hT chΓR (∼1 μeV), (iii)hT c > hΓR (∼0.1 μeV), and (iv) hT chΓR (∼0.01 μeV), and ΓL ≳ ΓR for all curves. The two dotted curves are the derivatives−dI/dɛ in arbitrary units for curves (i) and (ii) to enhance the bump structure. (B) Logarithmic-linear plots for (i) and (ii). Dashed lines are Lorentzian fits. For (ii) we chose parameters that fit the tail for negative ɛ. (C) Logarithmic-logarithmic plots of the emission spectrum for two different samples. The solid lines are based on the focused ion beam sample in Fig. 1A [upper trace is the same as (ii) in (A); lower trace is for coupling energies between (i) and (ii) in (A)]. The dotted line is based on a surface gate sample with a distance between left and right barriers of 600 nm (12). The dashed lines indicate a 1/ɛ and 1/ɛ2 dependence expected for piezoelectric interaction with 3D and 2D phonons, respectively.

The importance of fluctuations in the environment on electron tunneling through quantum devices has been recognized for a long time. Environmental studies on Coulomb blockade devices with two or more tunnel junctions have only discussed effects due to absorption (3, 17). For emission it is required that electrons are first pumped to a higher energy state, as has recently been done in a superconducting Cooper pair transistor under microwave irradiation (18). In the case of a double dot, pumping occurs whenE L > E R and an electron tunnels in from the left reservoir toE L. A double dot thus offers a unique two-level system that is pumped by a dc voltage without the induction of heating currents. It is therefore possible to reach an out-of-equilibrium situation so close to T = 0 that vacuum fluctuations become the main source for generating electron transport.

To identify whether photons, plasmons, or phonons form the bosonic environment (19), we measured spontaneous emission spectra while placing the double dot in different electromagnetic environments. In the regime from 10 to 100 μeV, the typical wavelengths are 1 to 10 cm for photons and 0.3 to 30 cm for 2DEG plasmons. We tested the coupling to the photonic environment by placing the sample in microwave cavities of different size (20). To check the coupling to plasmons, we measured different types of devices with largely different dimensions of the 2DEG leads, gate pads, and bonding wires. Both types of variation had no effect on the emission spectra; even the fine structure was reproduced.

The third option of acoustic phonons is the most likely possibility (21). Phonon emission rates have been calculated for single dots (22). For a double dot system, we can obtain the general energy dependence (15). For a deformation potential, we expect a rate dependence of ɛD−2 (ɛ for 3D phonons and constant for 2D phonons) and for piezoelectric interaction of ɛD−4 (1/ɛ for 3D phonons and 1/ɛ2 for 2D surface acoustic waves) (23). InFig. 3C traces measured on two different types of devices are compared. Here, the emission current is plotted versus ɛ on a log-log scale. Ignoring the bumps, we found an energy dependence between 1/ɛ and 1/ɛ2. This result implies that the dominant emission mechanism is the piezoelectric interaction with 2D or 3D acoustic phonons. Note that a 1/ɛ or 1/ɛ2 dependence should be avoided in coherent devices, because the inelastic rate becomes large near resonance (ɛ ∼ 0) (16).

The bumps observed in both type of devices suggest the existence of resonances, for instance, due to a finite size in the phonon environment. The bumps are particularly clear in the derivative of the current to energy (dotted curves in Fig. 3A). The large bump in Fig. 3A at ɛ = 30 μeV corresponds to a frequency of f = ɛ/h = 7.3 GHz. For 3D phonons, this yields a wavelength λ3D =s 3D/f = 640 nm (s 3D = 4800 m/s is the 3D sound velocity), whereas for 2D surface acoustic waves λ2D =s 2D/f = 380 nm (s 2D = 2800 m/s). Both wavelengths, more or less, fit with the dimensions of the two quantum dot devices. We have not yet been able to control these resonances by studying devices with a variety of gate dimensions. However, we predict that it will be possible to gain control over the phonon environment by making 3D phonon cavities in hanging bridges (24) or by creating a 2D phonon band gap with a periodic gate geometry (25).

  • * To whom correspondence should be addressed. E-mail: leo{at}


View Abstract

Stay Connected to Science

Navigate This Article