We propose a new Auger-like mechanism for energy relaxation in quantum dots (QD) driven by resonant scattering of delocalized wetting layer (WL) carriers. It is demonstrated that resonant scattering leads to a considerable increase in the relaxation rate that can explain experimentally obtained relaxation rates. Analytical results for the relaxation rate are obtained for rectangular dots revealing a weak logarithmic dependence on the dot depth and level density. Comparing results for a rectangular and a parabolic QD model we conclude that the relaxation rate is not very sensitive to a chosen model.
1. Introduction
Capturing and subsequent relaxation of hot carriers (quasi-electrons and holes) in nanoscale QD are of central importance for practical applications of QD systems, for example, in lasers [1] single photon emitters [2] and photodetectors with unique properties [3, 4]. The discussion of the corresponding rate is controversial. It is long known that a singular spectrum in QD inhibits incoherent emission of phonons suppressing the most effective relaxation mechanism [5], so the term phonon bottleneck was coined. Indeed, experiments [6, 7] demonstrate a rapid enhancement of the relaxation rate when interlevel spacing coincides with phonon energy. However, in many experiments the bottleneck did not reveal itself and the relaxation rate was comparable or even faster than provided by phonons in bulk systems [8–15] possibly indicating the role of the continuous spectrum [16–19]. Further studies lead to a conclusion that the processes of capturing and relaxation can be separated and that the relaxation is likely to be defined by different mechanisms in specific situations for both capture and relaxation [13, 18, 20, 21].
Several relaxation scenarios were proposed to explain these experimental results: phonon-assisted tunnelling to impurity states [22], decay of the LO phonons [23], and several Auger-like mechanisms which included particle collisions [24], QD coupling with the WL plasma excitations [25], phonon-assisted intradot Auger interactions [26], and Auger interactions between excitons in the dot and in a localized state of the wetting layer [27]. Experimentally observed dependence of the relaxation rate on the number of created excitons prompted the conclusion that Auger-like processes are likely to play a role when exciton density is n>109 cm−2, while other mechanisms are dominant at smaller density [13, 20, 28]. Calculations of the rate due to WL carrier collisions [24, 29, 30] and plasma excitations [25] demonstrated that these mechanisms become important at n>1011 cm−2 leaving the phonon-assisted intradot Auger interaction mechanism [26] as the likely candidate to explain the experimentally observed relaxation at lower densities [13].
However, as we show below, the Coulomb interaction between delocalized WL and trapped carriers can provide an efficient relaxation mechanism yielding a high relaxation rate even for as low as n=109 cm−2 density when the resonant scattering of the WL carriers is considered. This mechanism is general and could be applied to describe relaxation to deep states in a general gaseous system. Here we apply it to the QD systems. In the discussion we calculate the rate for two models of the confining potential: a rectangular model for which analytical results are obtained and a parabolic potential for which we do numerical simulations. Qualitative agreement between the two models is demonstrated, showing the results are not very sensitive to details of the chosen model. By substituting standard parameters for the InGaAs dot systems we obtain relaxation rate which agrees with experiments.
2. Relaxation Rate
Scattering of delocalized WL carriers induces transitions of the trapped carrier between the QD states leading to relaxation in the system the rate of which is obtained from the scattering cross-section as(1)dεαdt=n∑β∫v2εβ-εασαβλnTλdλ,where σαβ(λ) is the scattering cross-section at energy λ of the WL carrier, α,β are initial and final two particle scattered states, where the energy of the trapped carrier is εα,β, respectively, v2 is initial velocity of the WL carrier, nT is temperature distribution of the wetting layer carriers, for which we take Boltzmann distribution, and the summation is done over all possible scattering states β. The cross-section σαβ(λ) includes electron-electron, hole-hole, and electron-hole scattering as well as carriers exchange. Hereafter all calculations are done in dimensionless Bohr units for InAs material.
Conventional calculations of σαβ(λ) within the Fermi rule approximation [24, 29, 30] overlook the fact that the energy dependence of the cross-section exhibits sharp resonances when E=λ+εα is close to energy of two particle virtual states where two carriers are temporarily confined in a QD during a time defined by the inverse width of the resonance. A schematic illustration of the resonant scattering as opposed to the direct scattering described by the Fermi rule is illustrated in Figure 1.
Schematic representation of direct (a) and resonant (b) scattering of a carrier with k2 momentum on a dot with a trapped carrier on a level εα. After collision the first carrier leaves carrying momentum K2 and the second one remains trapped on a new level εβ. Resonant scattering involves formation of a virtual level j.
High density of resonances for the scattering of WL carriers with energy in the temperature interval is the main reason why resonant scattering may dominate the relaxation. For weak Coulomb interaction this density is estimated as the density of two particle states obtained under the condition that transitions to the delocalized WL states is neglected. For example, if we assume that a QD has an equidistant (degenerate) single particle spectrum with the gap δe=30 meV for electrons and δh=15 meV for holes and U=200 meV depth, the number of respective single particle states is Ne=20 and Nh=80 while the total number of excitonic states is approximately NeNh=1600. The states in the energy interval -2U<E<-U are stable, whereas the states -U<E<0 are virtual if transitions to the WL states are taken into account. The number of such states gives the number of resonances in the cross-section. The degeneracy of the virtual states is lifted by the intradot Coulomb interaction and by deviations from the symmetric shape of the dot. While these factors have only a minor effect on the low lying stable states, the energy levels of the higher excited states especially of the virtual states become randomized and insensitive to details of the single particle spectrum. As a result, positions of the resonances can be described by an averaged density of virtual states obtained by assuming a quasi-continuous approximation for the single particle levels which is valid despite the point like single particle QD spectrum. For E=-U/2, that is, when a trapped carrier has energy ε≈-U/2, the mean distance between the resonances is estimated as 2U/NeNh=0.25 meV ≈3 K and, therefore, at T≫3 K many resonances contribute in (1) enhancing the relaxation.
In the vicinity of a resonance j, the cross-section is given by the 2D Breit-Wigner expression [31](2)σαβjλ=4m2v2γjαγjβλ+εa-Ej2+1/4γj2,where Ej is the energy of the resonance, γjα is the partial width of the resonance, γj=∑αγjα is its total width, and λ and εα are the energy of the scattered and the trapped carriers, respectively. Partial widths are found as(3)γjα=2π∑kVjα2δEj-εα-λ,where Vjα is the transition matrix element between the virtual state j and a scattering state α and Ej=εj+λj is the zeroth-order approximation for the resonance energy, where εj and λj are single particle QD levels for temporarily trapped carriers. Equation (1) is valid if neighboring resonances do not overlap; that is, γj<Ej-Ej+1. If the temperature is not too low γj≪T the summation over λ in (1) for the cross-section (2) is carried out yielding for the rate(4)dεadt=8πnm2∑jβγjαγjβγjεα-εβnTEj-εα,where the condition εα<Ej holds. In the quasi-continuous approximation (1) can be further simplified if the density of resonances does not change within the δE=T interval. The number of contributing resonances in the sum in (1) can be approximated as T/Δ¯, where Δ is the average gap between neighboring resonances. This new T factor cancels the normalization T-1 in the Boltznmann distribution and the latter can be approximated by the delta function nT=δ(λ).
3. Coulomb Interaction Matrix Elements
For a rectangular shaped l×l QD single particle eigenfunctions are plane waves defined by a 2D wave vector k. Assigning indexes 1 and 2 to the trapped and free carrier, respectively, the zeroth-order two particle states are defined as pairs of single particle states. Coulomb transition matrix elements and the resonance widths are calculated with the following assumptions: (1) the confinement length in a direction perpendicular to the wetting layer is small and is disregarded in the calculations, which are strictly two-dimensional, (2) the eigenstates are taken as exp(ikx), which is equivalent to periodic boundary conditions, and (3) the Coulomb interaction is weak and acts only inside the dot, which allows us to expand the corresponding integrations over the entire space (with a proper normalization for the eigenfunctions). The latter assumption yields the already mentioned quasi-continuous approximation for the matrix elements which will be used in the sums over the discrete states in (1). With these assumptions (3) yields(5)γjαq=2π3l41q2δεj+λj-εα-λ,where q=K1-k1=k2-K2. It is convenient to introduce the quantity(6)γ¯jα=l22π2∫γjαδεj-εα-Δdq=2πl21εα-εj2v1jv2jsinΘv1j2+v2j2-2v1jv2jcosΘ,where Θ is angle between k1j and k2j and v1(2)j=k1(2)j/m1(2) is carrier velocity. Equation (6) defines the averaged width that corresponds to a particular energy change, Δ, of the trapped carrier. In deriving (6) the momentum transfer was assumed small; |k1j-k1α|≪k1(2)j. As what follows from (6) transitions to neighboring states dominate the relaxation. Using (6) we write the rate (1), averaged over the degenerate initial and final states α and β, as(7)dεαdt=16π2nm1m2l2∑j,εβγ¯jαγ¯jβγjεβ-εαδεj+λj-εα.The summation over εβ and j is done in the quasi-continuous approximation where sums are reduced to integrals. The limits of integration are obtained from the conditions εj+λj≥εα(β) and εα(β)<0. The final expression for the rate reads (8)dεαdt=-8n∫εα+δ0dεε-εaε+Uεα+UlnU+εε-εα-1×lnv1+v2v1-v2,v1=2U+εm1,v2=2U+εa-εm2,where δ is cutoff equal to the QD single particle energy gap between εα and the next lower level. The logarithmic singularity at v1=v2 is an artifact of the approximation used in calculating the transition matrix elements; however, it disappears after integration. If a lighter WL electron scatters on a heavier trapped hole, the logarithmic singularity is far from ε=εα most contributing to the integral and reduces to a small factor v1/v2<1 in the resulting rate.
The integration in (8) can be easily done numerically; nevertheless, it is instructive to evaluate its largest part on the lower integration limit analytically which yields(9)dεαdt=-8nlnv1+v2v1-v2π26+lnUδ-εαU+εα×12ln-εαδUU+εα+1,where v1=2(U+εα)/m1 and v2=2U/m2. In the limit of U/δ≫1 and when the trapped carrier is still on the higher QD level, that is, -εα≪U, (9) yields(10)dεαdt∝-lnUδ2,where U/δ is roughly equal to the total number of single particle QD states. One can see that the largest contribution to the relaxation rate depends logarithmically on the QD depth and the density of its single particle levels and is not sensitive to other details of the QD model. This gives an indication that the obtained result has an applicability range beyond the rectangular model QD.
To compare the resonant scattering contribution to that of the direct scattering (Figure 1(a)) we calculate the cross-section using the Fermi golden rule approximation. Substituting the result into (1) and using the same assumptions as before we obtain(11)dεadt=-nlogv1+v2v1-v2logU+εaδ∝lnUδ.As (10), this expression also depends on the QD depth and level density logarithmically; however, here the logarithm is not squared. Consequently for deep dots with a large number of levels, U/δ≫1, resonant scattering always dominates the relaxation.
Now we take a parabolic QD model where the confining potential is defined by V(r)=-U+m1,2ω1,22r2/2 at r<a=2U/m1,2ω1,22 and V=0 otherwise. As above we assume that the interaction acts mainly inside a dot, which allows using harmonic oscillator eigenfunctions for both QD and WL states. The two particle states are defined by four quantum numbers j=(n,m;N,M) and α=(r,s;R,S), where n, N, r, R are the main quantum numbers and m, M, s, S are the angular momenta of entering single particle states. The number R which describes the scattered particle is given by the energy conservation condition R=(U+λ1,2)/ω1,2 and can be, therefore, noninteger, leading to a divergency of the corresponding wave function at r→∞. However, in the calculation of the matrix elements this divergency is suppressed by the wave function with the main number N that corresponds to a lower energy. The transition matrix element can be expressed as an integral(12)Vnm;NMrs;RS=-iM-Sim-sδm+s-M-S×∫0∞dqϕn+mr-n+s-mq24m1ω1ϕnr-nq24m1ω1×ϕN+MR-N+S-Mq24m2ω2ϕNR-Nq24m2ω2,where ϕnm(x)=exp(x/2)Lnm(x)xm/2 is the radial part of the oscillator wave function (Lnm(x) denotes the Laguerre functions). Equation (12) is substituted into (3) and (1) that are evaluated numerically.
In the above analysis we assumed the scattering process at which the carriers are of a different type. For the carriers of the same type the cross-section is obtained by carrier permutation in the scattering amplitude f(1,2) as σ=f1,22+f2,12±2|f(1,2)f∗(2,1)|. It follows from the calculations f1,22 and f2,12 provide an equal contribution to the rate, while the interference term 2|f(1,2)f∗(2,1)| leads to a much smaller contribution because it involves transitions with larger energy transfer. Consequently, the rate acquires an additional factor 2. Another factor 2 comes from spin degeneracy of WL carriers. It must be noted that when carriers are of the same type the Coulomb blockade introduces an additional exponential factor exp(-Ec/T) into the rate. The value of Ec for the InAs dots can be estimated as few meV and, therefore, the Coulomb blockade is noticeable at low temperatures.
4. Results and Discussion
Results of calculations for the relaxation rate are summarized in Figure 2. Direct comparison between considered models is not possible since they have a different structure of single particle levels. However, we utilize the fact that the level structure dependence in (9) and (10) comes only via the energy gap between neighboring states and that the rate is mainly defined by transitions to the nearest levels. It has to be noted that in the parabolic potential zeroth-order two particle states eigenvalues are degenerate, while in real dots this degeneracy is lifted as discussed above. To bypass this difficulty we formally assume a slightly different energy gap for the second particle δ1=1.1δ2, which was enough to break the artificial degeneracy of the resonances imposed by the degeneracy of single particle QD states. We note that the quasi-continuous approximation is not needed for the parabolic model, where the summation in (1) is done numerically. The temperature dependence of the rate is weak as long as Δ¯≪T≪ω, for the calculations we chose T=40 K.
Relaxation rate versus energy of a trapped carrier, calculated for parabolic (symbols) and rectangular (continuous lines) dots for electron-electron (e-e) and hole-hole (h-h) scattering.
The results in Figure 2 demonstrate that the relaxation due to resonant scattering is considerably faster for deep QD levels. One can also see that resonant scattering driven relaxation for both considered models is qualitatively similar, while quantitatively the models differ only by a factor 1.5. The rate increases in the interval εα<δ, then saturates, then abruptly drops, and becomes zero at εα>-δ, where virtual states nor longer exist. The relaxation rate of heavier holes, which have smaller energy gaps, is faster than that of electrons, as expected from (11) and (10).
An effective relaxation time, defined as the time when a carrier with initial energy -δ reaches the state εα=-U/2, is found by solving (1) with the right-hand side given by the data plotted in Figure 2. Thus we obtain τ≈5 ps for the hole-hole scattering and τ≈15 ps for the electron-electron scattering when the WL carrier density is n=109 cm−2.
We have demonstrated that resonant scattering considerably enhances relaxation of a single carrier already captured in a QD below a critical level -δ. This initial capturing may be considered a bottleneck of this relaxation process. Few capturing mechanisms can bring the carrier onto those levels, for example, phonon emission by a WL carrier [32]. One may note, however, that the bottleneck appears entirely due to the assumption that the single particle states energy gap remains constant for higher QD levels. This assumption is not applicable for many realistic QD models, where the gap decreases for the upper states, diminishing the bottleneck interval and thus enhancing the capture rate without assistance by phonons. Moreover, the presence of localized states near the QD makes the capture onto the shallow QD levels more efficient [27].
When the energy of a trapped carrier is close to the ground state the relaxation rate rapidly decreases. The resonance widths become smaller leading to a larger life time of two particle virtual states. However, in this case scattering of WL carriers on a QD with two carriers inside, which is also of a resonant character, ensures a fast relaxation of the two particle trapped states. Similar reasoning can be also applied to relaxation of a general multiexciton complex [12]. Estimation of the corresponding rate may be complicated since the Breit-Wigner formula is likely inapplicable because of the high density of resonances. A detailed discussion of these questions is beyond the aims of this work.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
BimbergD.GrundmannM.LedentsovN. N.1999New York, NY, USAWileyMichlerP.KirazA.BecherC.SchoenfeldW. V.PetroffP. M.ZhangL.HuE.ImamogluA.A quantum dot single-photon turnstile device200029055002282228510.1126/science.290.5500.2282WangW. W.GuoF. M.LiY. Q.Modeling and simulation of a resonant-cavity-enhanced InGaAs/GaAs quantum dot photodetector20152015684751010.1155/2015/847510BhattacharyaP.MiZ.Quantum-dot optoelectronic devices2007959172317402-s2.0-4554909316010.1109/JPROC.2007.900897BockelmannU.BastardG.Phonon scattering and energy relaxation in two-, one-, and zero-dimensional electron gases19904214894789512-s2.0-000093205810.1103/PhysRevB.42.8947FafardS.LeonardD.MerzJ. L.PetroffP. M.Selective excitation of the photoluminescence and the energy levels of ultrasmall InGaAs/GaAs quantum dots199465111388139010.1063/1.112060FarfadS.LeonR.LeonardD.MerzJ. L.PetroffP. M.Phonons and radiative recombination in self-assembled quantum dots199552HeitzR.VeitM.LedentsovN. N.HoffmannA.BimbergD.UstinovV. M.Kop’evP. S.AlferovZ. I.Energy relaxation by multiphonon processes in InAs/GaAs quantum dots19975616104351044510.1103/PhysRevB.56.10435SosnowskiT. S.NorrisT. B.JiangH.SinghJ.KamathK.BhattacharyaP.Rapid carrier relaxation in In0.4Ga0.6As/GaAs quantum dots characterized by differential transmission spectroscopy199857R942310.1103/PhysRevB.57.R9423RaymondS.HinzerK.FafardS.MerzJ. L.Experimental determination of Auger capture coefficients in self-assembled quantum dots20006124R16331R1633410.1103/PhysRevB.61.R16331CruzJ.de SalesF.da SilvaS.SolerM.MoraisP.da SilvaM.QuivyA.LeiteJ.CW photoluminescence determination of the capture cross-section of self-assembled InAs quantum dots20031710710810.1016/S1386-9477(02)00721-XNguyenH. S.SallenG.AbbarchiM.FerreiraR.VoisinC.RoussignolP.CassaboisG.DiederichsC.Photoneutralization and slow capture of carriers in quantum dots probed by resonant excitation spectroscopy2013871111530510.1103/PhysRevB.87.1153052-s2.0-84874882139KlimovV. I.McBranchD. W.LeatherdaleC. A.BawendiM. G.Electron and hole relaxation pathways in semiconductor quantum dots199960, article 1374010.1103/PhysRevB.60.13740MarcinkeviciusS.LeonR.Carrier capture and emission in self-assembled InGaAs/GaAs quantum dots1999797982MaZ.PierzK.HübnerJ.RühleW. W.Nanosecond scale carrier dynamics of self-assembled InAs/AlAs quantum dots studied by time-resolved photoluminescence20052832032082-s2.0-2124449863410.1016/j.physe.2005.02.009TodaY.MoriwakiO.NishiokaM.ArakawaY.Efficient carrier relaxation mechanism in InGaAs/GaAs self-assembled quantum dots based on the existence of continuum states19998220411441172-s2.0-034499202910.1103/PhysRevLett.82.4114FerreiraR.BastardG.Unbound states in quantum heterostructures2006121201362-s2.0-3384686704410.1007/s11671-006-9000-1LingH. S.LeeC. P.Lo JournalM. C.Energy-dependent carrier relaxation in self-assembled InAs quantum dots200810312431110.1063/1.2947599RihaniJ.SedrineN. B.SalletV.OueslatiM.ChtourouR.Observation of many-body Coulomb interaction effects on the photoluminescence spectra of InAs/GaAs quantum dots20082541031253129MarcinkeviciusS.GaardnerA.LeonR.Rapid carrier relaxation by phonon emission in In0.6Ga0.4As/GaAs quantum dots200164115307ShamirzaevT. S.AbramkinD. S.NenashevA. V.ZhuravlevK. S.TrojánekF.DzuňrákB.MalýP.Carrier dynamics in InAs/AlAs quantum dots: lack in carrier transfer from wetting layer to quantum dots2010211515570310.1088/0957-4484/21/15/1557032-s2.0-77949892305SercelP. C.Multiphonon-assisted tunneling through deep levels: a rapid energy-relaxation mechanism in nonideal quantum-dot heterostructures19955120145321454110.1103/PhysRevB.51.145322-s2.0-35949007644LiX.-Q.ArakawaY.Phonon bottleneck in quantum dots: role of lifetime of the confined optical phonons199959506910.1103/PhysRevB.59.5069UskovA. V.AdlerF.SchweizerH.PilkuhnM. H.Auger carrier relaxation in self-assembled quantum dots by collisions with two-dimensional carriers19978112789578992-s2.0-000117815110.1063/1.365363BockelmannU.EgelerT.Electron relaxation in quantum dots by means of Auger processes1992462315574155772-s2.0-000100285210.1103/PhysRevB.46.15574EfrosA. L.KharchenkoV. A.RosenM.Breaking the phonon bottleneck in nanometer quantum dots: role of Auger-like processes199593428128410.1016/0038-1098(94)00760-82-s2.0-0029197612SyperekM.BaranowskiM.SękG.MisiewiczJ.LöfflerA.HöflingS.ReitzensteinS.KampM.ForchelA.Impact of wetting-layer density of states on the carrier relaxation process in low indium content self-assembled (In,Ga)As/GaAs quantum dots201387125305SunK. W.ChenJ. W.LeeB. C.LeeC. P.KechiantzA. M.Carrier capture and relaxation in InAs quantum dots2005169153015352-s2.0-2344443993310.1088/0957-4484/16/9/021NielsenT.GartnerP.JahnkeF.Many-body theory of carrier capture and relaxation in semiconductor quantum-dot lasers200469235314NilssonH. H.ZhangJ.-Z.GalbraithI.Homogeneous broadening in quantum dots due to Auger scattering with wetting layer carriers200572205331LandauL. D.LifshitzE. M.1972Oxford PressMR0120782FereiraR.BastardG.Phonon-assisted capture and intradot Auger relaxation in quantum dots199974281810.1063/1.124024