NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2021, 12 (6), P. 680-689
Temperature dependence of recombination radiation in semiconductor nanostructures with quantum dots containing impurity complexes
V. D. Krevchik1, A. V. Razumov1, M. B. Semenov1, I. M. Moyko1, A. V. Shorokhov1,2
1Penza State University, Krasnaya, 40, Penza, 440026, Russia ^University of Jyvaskyla, PO Box 35, Jyvaskyla, FI-40014, Finland
physics@pnzgu.ru, alex.shorokhov@gmail.com
PACS 73.40.Gk, 03.65.Xp DOI 10.17586/2220-8054-2021-12-6-680-689
Temperature dependence of the spectral intensity of recombination radiation in a quasi-zero-dimensional structure, containing impurity complexes “A+ + e” (a hole localized on a neutral acceptor, interacting with an electron localized in the ground state of a quantum dot), has been investigated in an external electric field in the presence of tunneling decay of a quasistationary A+-state. Probability of dissipative tunneling of a hole has been calculated in the one-instanton approximation, and the influence of tunneling decay and of an external electric field on the A+ -state binding energy and on the spectra of recombination radiation, associated with the optical transition of an electron from the ground state of a quantum dot to the A+ -state of the impurity center, has been investigated in the adiabatic approximation. “Dips” in the temperature dependence of the SIRR have been revealed, which are associated with the presence of resonant tunneling at certain values of temperature and strength of the external electric field, for which the double-well oscillatory potential becomes symmetric.
Keywords: spectral intensity of recombination radiation, quasi-zero-dimensional structure, impurity complexes, quantum dots.
Received: 15 November 2021 Revised: 19 November 2021
1. Introduction
The relevance of studies of the spectral intensity of recombination radiation (SIRR) temperature dependence in quasi-zero-dimensional structures is determined by the fact that, first, electro-optical systems based on quantum dots (QD)s have significantly better parameters compared to similar devices based on quantum wells. Second, the temperature quenching of luminescence is usually associated with the presence of a nonradiative channel caused by a defect located in a QD [1,2] or at the barrier boundary. In this work, we would like to draw attention to the possible existence of one more channel of temperature quenching of luminescence associated with tunneling processes, in particular, with the process of dissipative tunneling of a hole localized at the A+-center into the matrix, surrounding the QD. As will be shown below, the temperature “spreading” of the wave function for the quasi-stationary A+-state under tunnel decay conditions is accompanied by an increase in the energy of the radiative transition of an electron and a corresponding decrease in the overlap integral of the wave functions of a hole and an electron, which leads to temperature quenching of the recombination radiation. The aim of this work is to study theoretically the temperature effect, associated with the electron - phonon interaction, on the binding energy of a hole in the “A+ + e”-impurity complex in a spherically symmetric QD, as well as on the spectral intensity of the recombination radiation of a quasizero-dimensional structure in the presence of a dissipative tunneling.
2. Binding energy of a quasi-stationary A+-state in a quantum dot in the presence of dissipative tunneling in an external electric field
Let us consider the problem of quasi-stationary states of a hole in an impurity complex A+ + e in a semiconductor spherically symmetric QD. The interaction of an electron in the QD ground state with a hole, localized at the A0-center, will be considered in the framework of the adiabatic approximation [3]. In this case, the electron potential Vn,l,m(r) acting on the hole can be considered averaged over the electron motion:
Vnlm (r)
№nlm (Ге) |2
4пєо£J |r - Ге|
0
dre
2
e
(1)
where e is the the electron charge; є is the dielectric permittivity of QD material; є0 is the electrical constant; m = 0, ±1, ±2,... is the magnetic quantum number; l = 0,1,2... is the orbital quantum number, Фп1т (re) is the wave
function of an electron in a QD determined by the expression (2):
^nlm (re) — Ylm (Өе? +e)
Jl + З (Knlre)
л/2пД0VJ 3 (KniRo)
(2)
where Knl is defined by an expression of the following form:
Kn
nl
lXn,i kT
R02 + Eh ln
4 sinh
----^=1 sinh
VbAs/kT J
2- У
vtaVKT )
(3)
Here vLA, vTA are the velocities of longitudinal and transverse phonons; - — n^G / JL, g is the overlap
a у pV
integral, Ө is the deformation potential, a is the lattice constant, p is the QD material density, V is the QD volume, Eh is the Bohr hole energy, JCn,l is the half-integer root of the Bessel function l + 1/2.
Let us consider the case corresponding to the p-state of an electron, in this case l — 1, m — 0, then expression (1) can be represented in the form
Vnio(rh)
2n
4пє0є
sin Өе dӨe / dpe
rh 2 2 Ro 2 2
f |^n10 (re, Өe, +e) | re2 dre f |^n10 (reA,+e)| re2 dre
+
\Jrh2 + re2 - 2rhre cos Өе J \Jrh2 + r 2
^o rh
re2 - 2rhre cos Өe
(4)
Or, taking into account the explicit expression for the wave function:
3e2
2n
Vnio(rh) — -
32n3£oeR‘oJl (KniRo)
2 o
sin Өe cos2 ӨedӨe dpe
J5 (Knire) r-edr-e 2
Vrh2 + re2 - 2rhre cos Өе
Ro
+
J5 (Knl re) re dre
Vrh2 + re2 - 2rh,re cos Өе
(5)
7T
2
e
After integration over the angular variables <pe, Өe and over the radial coordinate of the electron re, we obtain:
Vnio(rh) —
27n3eoeRo6Knirh4 J1 (Kni r*h)
[3R04 + 6KniR04rh2 + 9rh4 + 6KniR02rh4
- 4KniRo4rh4 - 9rh4 cos (2KniR0) + 12K^R0 2rh4 cos (2K^R0)
- 3Ro4 cos (2K^R0) + 4KniRo4rh4 (Ci (2K^R0) - Ci (2Kn+h))
-4KniR04rh4 ln ( R ) - (Qrh3 + 3R03 - KniR0 3rf) 2K^R0 rh sin(2KniR0)
(6)
Here, rh — rh/ah, R0 — Ro/ah,
2
e
Kn —
X2
n,l kT — + ^ „.,.o ln
R02 EhR02
. - \
4 sinh I ----I sinh
vlav kT /
2- У
vta
VkTA
Ci (x) is the cosine integral.
In the case under consideration, p-state with m — 0, the potential energy minimum is shifted relative to the QD center. Position of the minimum rh min is determined from the solution of the transcendental equation, which is obtained by equating the first derivative of expression (6) to zero:
dVn,i,o(rh)
drh
Expanding expression (6) in a Taylor series near rh min and limiting ourselves to the quadratic term in the radial coordinate of the hole rh, we obtain:
Vnio(rh)
e22-7n-3
n2min( ^nWrhLn (3R*2+зг^ - ^nWrhLn))+
Д*
o
_ e0£R06Kn51rh6minJ52/2 (K*5R0 )
+ 4K:5R04rh4mi^Ci (2КПіД0) - Ci (2K>h min) - In - 2КПдҺ minR0 (K;5R03rh2min sin (2К*5^ min) - ^fmin Sm (2K*5R0) - 3R3 ^ (2К*5^ min) ) +
гі1
+ (rh - rh min) 2R04rh min ( - 6 - 6K^min + 2K*5rh4min +6 (1 - К*^^ COS (2К*^ min) + + K*irh min (K*5rh2min - 12) Sin (2К*^ min) )-- (rh - rh min)2 2R04 ( - 15 - 9K*5rh2min + ^4^ +
+ (15 - 21K*5rh2min + K:5rh4min) COS (2K* prh min) + (30K*^ min - T^rfm^) sin (2K* ^ min) )
where rh min = rh min/ah.
Within the framework of the adiabatic approximation for the p-state of an electron, we obtain:
e2 ~ тһ< (rh - rT)
г5Гһ min)
* r
Vn50(rh) = T-,* в*
eR;
0
where values в*, and rT are defined as follows:
2-7n-3
вп R * 5 k * 5 r * 6 j2 (K * R * )
R0 Kn1 rh min J5 /2 (Kn1R0 )
rh2miJ2K*5R02rh2min (3R02 + 3rh2min - 2K*5l02rhLJ ) +
0 Kn5 'h min J5/2 (Kn5R0) |_ V 7
+ 3R04 + 9rh4min + 3rh4min COS (2K*^) (6K*5R02 - 3) - 3R04 COS (2K*prh min) +
+ 4K:lR*4rh4min Ci (2K*iR0) - Ci (2K*^ min) - ln
R
n
h min
(8)
(9)
- 2K*prhminR^K*5R03rh2min Sin (2K*і'Һmin) - 9rh3min Sin (2K*iR0) - 3R03 Sin (2K*і'Һmin) ) + + 2R*4rh2min ( - 6 - 6K*5rh2min + 2^^ + 6 (1 - K*5^) COS (2K*і'Һ min) +
+ K*5'Һ min (K*5rh2min - 12) Sin (2K*5^ min) ) X (15 + 9K*5rh2min - ^l^min-
- COS (2K*5'Һ min) (15 - 21K*5rh2min + K*5r^min) - (30K*5'Һ min - 7K*?rh3min) Sin (2K*prh min) ) ^
(10)
2
hwn
h2e22 5n 3
mh£0£R02K;5rh6min J/2 (Kn5І0)
(15 + 9K*5rh2min - K*5r
4 4
h min -
- COS (2K*prh min) (15 - 21K*5rh2min + K^in) - (30K*5'Һ min - 7K*3rh3min) Sin (2K*5'Һ min) )
and
't = — r
h min
( - 6 - 6K*5'h2min + 2K:5rh4min + 6 (1 - K*5^^) COS (2K*5^ min) + + K*5'h min (K*5'h2min - 12) Sin (2K*5'h min) ) X [ - 15 - 9K*5^п + K*5^п +
r
h min 2
1 5/2
(11)
+ (15 - 21K*5'h2min + K*5rh4miJ COS (2K*5'h min) + (30K*5'Һ min - 7K*3rh3min) Sin (2K*5'h min)
T - 5
Since the confining potential of a QD, generally speaking, should have a finite depth, then in our model of the hole potential of confinement (9), the amplitude of the potential U0 is an empirical parameter and satisfies the relation
Uo = -е2в„/4пєоєДо + ш*һш2п R0/2 = ш^2Д2/2,
whence w0 = ё2вП/2пё0ётҺд3 is the characteristic frequency of the hole confining potential of a QD within
the adiabatic approximation, and in this case U0/ (Һш0) > 1.
Usage of the adiabatic approximation makes it possible to take into account the effect of an external electric field on the bound state of the hole. Let the electric field strength vector E0 be directed along the x coordinate axis, then the energy levels of the oscillatory potential (9) will have the form:
En,0,0 (T) = _ e e IеI E0
En\п2пз (T )= £R0 РҺ 2mh^
+ bwn
ni + П2 + П3 + 2 ) +
+ kT ln
4sinh
2П
(vL^/kT J \vTAsZkT J _
z ) sinh ( -
and the corresponding one-particle wave functions are written as:
(x - X0)2 + y2 + z2 \
ФПШ2П3 (x,y,z) = CneXP --
2 a?
Hn
x - x0
an
HnJ— Hn. -
(13)
(14)
Here, Cn
2пі+п2+пЗщ і П2! П3! n3/2a,n
-1/2
; an = у Һ/ (m*hUn)\X0 = |e| E0/ (mh^n); Hn (x) are the Hermite polynomials; nl, n2, n3 are the quantum numbers corresponding to energy levels of a harmonic oscillator. We will assume that the process of decay of the quasi-stationary level of the A+-center is due to dissipative tunneling. It should be noted that in the one-instanton approximation, the decay probability Г0 (dissipative tunneling), under conditions of an external electric field, can be represented in the form Г0 = B exp (-S), where the expressions for S and B were obtained in our previous papers (in the Bohr units) [4,5]:
s = 1( b0 + x0
2 M0 + x0
'(1 - x2)
+1
3 -
b0 + x0 a0 + x0
1
b0 + x0
2 в’ \ a0 + :
+ 1) T02 -^
1
b0 + x0
2ү* \a0 + x0
+1
coth (h-Va) - cosh((e'-T°> M -f”hPVx^> + cosh (,e- - t0) 77)
Y v V sinh (вN /A) Vv^ 02 v V
(1 - xl)
coth (ТРЩ - COSh((e’ - ^ + cosh ((в' - T- ) )
(15)
B=
2 Ed
ущ (
b0 + x’
+1
where
Һ ypn \a0 + x’
+ A’ ^1 -1 A
£T X < A’
в’ cosh ( у ) - 1
+ D’
в2 cosh - 1
el
■ T0i
2 sinh ^ У )
в'
+ D’
в’ cosh ( ^ - Г02'
sinM 7?
cosh ( e— t0 l'
sinM 7?
1
+D
в’
1
cosM - T02
+
sinM
1
x l ,2
1
" 2
*
Y
+1+
Pc a ’ 2
T
' є*2 a ’ 2
+1+
4U’ 4є*ь2Щ ' У V 4U’ 4є*ь2Щ J U’
£*L2a ’2
^(eL2a ’2/ (4U’) + 1 + єС 4a’2/ ^є^Щ))2 - є*2a’2/U’,
arsinh
1
b0 + x '
sinh (в ’) / ( 1 +
b0 + x 0
+ в’
a'0 + x ’ J \ a'0 + x 0
єТ = kT/Eh, є* = fe^/Eh, єС = hVC/Ed, в’ = УЩ /a ’ At , b0 = b0/ah, a0
l/2
(16)
a0/ah, x0 x0/ah;
3
a
a
n
n
2
2
2
2
2
2
є ’2a ’2
’2
Eh and ah are the Bohr energy and the hole radius correspondingly;
A* = (2e*2a*2 - x*) / ((x* - x*) X), D* = (2e*2a*2 - x*) / ((x* - x*) x*),
e* = V2 xi/ (a* £t), = V2 x 2 / (a*£t), Toi = To V~x[/E2, T02 = "biVE/V'2.
Using the zero-range potential procedure (see, for example, [3]), in combination with the one-instanton approximation [5], we obtain an equation that determines the dependence of the hole energy in the complex A+ + e on temperature T, QD parameters, and dissipative tunneling.
The short-range potential of an impurity is simulated by a potential of zero radius with a power ү = 2nh2 / (am*h), which has the form:
Vs (x, y, z; xa, ya, Za) = Y6 (x - xa) 6 (y - ya ) 6 (z - Za)
' , , д
1 + (r - ra) dr
(17)
where a is determined by the binding energy Ei of the same A+ -center in a bulk semiconductor.
In the effective mass approximation, the wave function ФЛҺ (x, y, z; xa, ya, za) of an electron localized at a short-range potential satisfies the Schrodinger equation:
(EXh - H) ФЛҺ (x, y, z; xa, ya ,za) = Vs (x, y, z; xa, ya, za) ^ЛҺ (x, y, z; xa, ya, za) , (18)
where EQD = -h2X2 / (2mh) - the eigenvalues of the Hamilton operator Hs = H + Vs (x, y, z; xa ,ya, za); H = -h2/(2m*h) V2 + mh^2 (x2 + y2 + z2) /2 - |e| E°x.
To determine the binding energy of a hole in a complex A+ + e in the adiabatic approximation, it is necessary to construct a one-particle Green’s function G (x, y, z; xa, ya, za; EЛп) to the Schrodinger equation with a Hamiltonian, containing potential (17):
G (x,y,z; xa,ya,za; Eлh)
E
n1,n2,n3
фп *
Пі,П2 ,П3
(xa,ya,za) Ф
п
Пі,П2,Пз
(x, y, z)
-eQd + ihro + En;0n°2,n3 (T)
(19)
The Lippmann-Schwinger equation for a A+-state in a QD with a parabolic confinement potential can be written
as:
1-00 + TO -bOO
фh (x,y,z; xa,ya,za)= j J j dxidyidziG (x,y,z; xa,ya ,z^; Елһ) x
X Vs (x,y,z; Xa,Уa,Za)^Лh (xi, yi, zi, xa, ya, za) . (20)
-TO —TO —TO
Substituting (17) into (20), we obtain:
фh (x, y, z; xa ,ya, za) = YG (x, y, z; xa, ya, za'; Елһ) (T Ф h) (x, y, z; xa, ya, za) ,
where:
(T ^ЛҺ )(x,y,z; xa,ya,za) = ДҺп
V^Va
ө^өа
л . ,д
1+(r - ra ) дГ
^ЛҺ (x,y,z; xa,ya,za) .
(21)
(22)
Acting by the operator T on both sides of relation (21), we obtain an equation that determines dependence of the energy Елһ of the bound state of the A+ -center on the parameters of the QD, on the position of the impurity
Ra = (xa,ya,za), and on the temperature T:
a = (TG) (xa, ya, za'; xa, ya, za'; Елһ) , (23)
mh
here, a is determined by the energy Ei of the bound state of the same A+ -center in a massive semiconductor.
Then, for Green’s function (29), taking into account (23) and (24), we obtain in the Bohr units:
G ix,y,z,xa,ya,za'; Г)Лһ
1 (x* - x0)2 + у *2 + z *2 + (x*a - x0)2 + y*a + z>a
exp I X
\2 I * 2 I ^ * 2 '
n3/2a?nEh
E
Пі,П2,Пз
Hni (x* - x*) Hm (xa - x°) Hn2 (y*) Hn2 (ya) Hn3 (z*) Hn3 (za).
2ni !ni!
2n2 !n2!
2n3 !n3!
2
-Vh - вҺ - 4ф—А + і4Г° + в i ( ni + n2 + n3 ^ о ) ^ ln
3\ kT. 2) + Eh
4 sinh (--) sinh ( 2E_ )!}
. \vLAVkT ) yvTA'/kT J J J
i
where the next notations are introduced: nh = EQD/Eh; R— = R0/ah; 3 = Eh/hwn; a*n = an/ah;
3— = е2вһ/£ЕҺ)Еһаһ; Г— = bFo/4Eh.
Further, given that:
-1
*2 0 — 2
-nh - 3Һ -
x0 3
kT
4 ahh + ,4Гһ + Eh 1,1
3 1 J dt exp
o
4 sinh
----^=1 sinh
VLAs/kT J
2-
xhh 3—2
-t( -3(nh + 3һ + -һОһ—г
VTAsfkT
- *4Г0 ) +
3
+ ni + П2 + П3 + 2
3 kT3
+П1 + П2 + П3 + - + 1n
2 Eh
the expression for the Green’s function takes the following form:
4sinh Г---^1=] sinh f 2-_")
. \VLAs/kT J \VTAVkT J _l_
(25)
G (X,y,Z,Xa,Va,Za,nh) = 3
1 exp I (x— -x—)2 + (xa- x—)2+y2+ya+z2 + zhV,
exp I ^ I X
n2 аПвЕһ
2an
xj dt exp
0
-t -3 % + 3— +
x—23 2 , kT3
4ah2
- І4Г0 +
Eh
1n
4 sinh
- \ / 2- VI 3
= I Sinh I ---;-- I +--
VLAVkTj {vTAs/kT J \ 2
ni =0
e—^niHni{нт[“) ^ fe—t\n>Man Man
ni!
E V
n2 = 0
E(V
пз=0
П2!
-^пз Hn z) Hn Za
П3!
. (26)
Making summation over the quantum numbers n1, n2, n3, and separating the divergent part of the expression for the Green’s function (26), we obtain:
G (r,Ra,nh) = -
1
n2 an 3Eh
+
0
kT3
~Eh
dt exp -t[ -3 ( nh + 3h + 4a—2
( - ) (
4 sinh -------sinh
2 x—23 2 з
2 1 «* 1 0 - *4Г0 +-----+
1n
2-
r2 + Д2
X ^ (1 - exp(-t)) — 2 exp(---2^a
exp
\VLAs/kT J V vtaVeT
f2(rRa) e—t - (r2 + Ra)
)])
a) e—2t
(1 - e—2t)
tit
exp
(r - Rg)2 2alt
exp ( V2 (-nh - 3— - Pf + -4Г— + g 1n [4sinh (sinh (+ з) ^
-------—----------------------------------3---------------------------------------------------. (27)
2n2 an3Eh |r - Ra|
Substituting (27) into (23) we obtain the dispersion equation that determines the dependence of the binding energy of a hole Eqd in the complex A+ + e on the parameters of the QD, on temperature T, and on the electron quantum number n. Using the procedure of the zero-range potential method (see, for example, [3]) in combination with the one-instanton approximation [5], we obtain the dispersion equation that determines dependence of the hole binding energy in the complex A+ + e on the temperature T and on parameters of dissipative tunneling:
— 2 0-2
x023
kT
ni = V -nh- 3h-iiar +і4Г— + Eh1n
4 sinh
-
VLA
VkT
sinh
2-
VTA
VkT
3
+ 2 +
+ v 3t.fdt exp{-t3 (-nh - 3h - ^
*2я—2.л^ 3 kT
+ *4Г— + — + — X
1n
4 sinh
-
VLA
л/kT
sinh
( 2-
\VTAs/~kT J _
1 - (1 - e—2t) — 2
2t%/2t
23 Eh
exp f-^ x '
a
where R*a = Ra/an; Пі = V\EiI /E, Пһ = \/E<fD)/Eh.
It should be noted that the binding energy of a hole in the case under consideration is a complex quantity. Its
real part determines the average binding energy of the resonant state of the A+-center EhQD) = Re E(QD), and the doubled imaginary part determines the broadening of the corresponding energy level AEh = 2Im E^QD). Fig. 1 shows the result of a numerical analysis of the dispersion equation for the case of a centered A+-center (R* = 0) at different values of the QD radius R*. It was taken into account that the binding energy of the A+-state is measured from the level of the ground state of the adiabatic oscillatory well.
Fig. 1. Dependence of the binding energy of quasi-stationary states of a hole in the complex A+ + e on the temperature T and the magnitude of the external electric field E0, for the following values of the QD parameters and dissipative tunneling R* = 1 (R* = R0/ah - QD radius R0 in units of the effective Bohr radius of the hole ah); U* = 350 (U* = U0/Eh is the amplitude of the adiabatic potential U0 of an electron in units of the effective Bohr energy of a hole Eh), n = 3.
As can be seen from Fig. 1, that in the field dependence of the binding energy of the A+-state, there are “dips” at a certain temperature. This is due to the effect of “tuning” the starting energy level of the A+-state to the states caused by the hole-phonon interaction in the matrix surrounding the quantum dot, i.e. with the effect of resonant tunneling. The depth of the dip increases with increasing temperature, which is due to the dynamics of the temperature-dependent peak in the field dependence of the probability of dissipative tunneling [5]. A decrease in the binding energy of the A+-state with an increase in the external electric field strength is associated with the Stark shift in energy and with the polarization of the A+-center, and with an increase in temperature - with the broadening of energy levels and the corresponding “spreading” of the wave function of the A+-state.
3. Temperature and field dependences of the spectral intensity of recombination radiation in a quasi-zero-dimensional structure with impurity complexes
Let us consider the process of radiative transition of an excited electron to the level of A+-center. The Coulomb interaction of an electron with a hole is accompanied by a radiative transition of an electron to the energy level of A+-center under conditions of tunneling decay of a quasi-stationary A+-state. The energy spectrum of an electron in the size-quantized band can be represented as [6,7]:
E
n,l
Xn,i Eh R*2
+ kT ln
4 sinh
о
VLA
VkT
sinh
20
VTA
VkT
(29)
here XXn,l is the root of the Bessel function of half-integer order l + 1/2.
The wave function of an electron is given by expression (2).
The SIRR, taking into account the dispersion of QD sizes and the finite lifetime of the resonant A+-state, is determined by an expression of the next form [8]:
4(u2 V^e2 c3 V
Peh e0
m о
/Е\MI' P (u)
nlm
Го
у0 + (Enlm — EXh —
where m0 - the free electron mass; Peh is the matrix element of the momentum operator at the Bloch amplitudes of band carriers; w is the frequency of the emitted electromagnetic wave with polarization e0; V is the QD volume; P (u) - the Lifshitz-Slezov function [9]:
34eu2 exp [-1/ (1 — 2u/3)]
P (u)=<| 25/3 (u + 3)7/3 (3/2 — u)11/3’ u< 2’ o, u > 2.
(31)
The wave function of A+ -state, in the case of a central location of the A+ center in the QD, has the form
Фh (r) = —c j dt exp j —fit^ —Һ — вҺ — + І4Г5 +
kT
Eh в
ln
4 smh I --smh I -------+
vlav kT J \VTAV kT ) _
+2в)}x (1—exp(—1)) 3/2 exp —
(1 + e-2t) r2 \
x к , (32)
where C is the normalization factor, determined by an expression of the next form:
c
2Vn
rf1 + щп\.an r
[nh (t )]2 W
nh (T) vh (T)
2
2
+ 1 — Ф
vh (T) 1
2
2 1
(33)
Here, 'rj\h (T) is determined by the dispersion equation (28).
The matrix element of the radiative transition of an excited electron to the level of A+ -center in the dipole approximation has the form:
M = i\o\l 2па І0 (En,l,m — Eh) (Фһ (r) |(eA,r)| (p, Ө, p)) .
(34)
Taking into account (2), (28) and (32), the matrix element of the radiative recombination transition of an electron from the ground state of the size-quantized band to the level of A+-center in the QD can be represented as:
M
a-12 2 n 4 V3® 2 V2kRQ4Jз (Xn,i
r(1 + ^ )
[nh (T )]2 Г (
Щ1 (^ҚЩТІ +л — ФЛШ — 1, — 1
і п 2п
■ 2
r*2dr* cos Ө sin ӨdӨdpx
0 0 0
x dt exp
0
—м — %в — eeh—
х02в , ,лот» , кТв
4
+ і4вГ0 +
Eh
ln
4 sinh
tin \ , ( 2Hn
vLAVkT) sin4VTA^kT
3
+ 2
rJl + 3 \Xn,l
x (1 — exp(—1)) 2 ex^— (1 + e_2t) x T-^j Yl,m (Ө,Р) 1+2 (R°
where RQ = R0/ah.
Calculation of (35) leads to integrals giving selection rules for quantum numbers m and l:
2n s
(35)
exp (imp) dp =
2n, if m = 0,
o, if m = o,
(36)
j Pl (cos Ө) cos Ө sin ӨdӨ = < 3
2, if l = 1, 3’ ’
0, if l = 1.
(37)
Thus, the radiative transition of an electron to the level of the A+-center is possible only from states with the values of quantum numbers l = 0 and m = 0.
1
2
The remaining integral over the radial coordinate r* has the form:
dr*r*3/2Jl+ 3 .
l+3 R
exp -
(1 + e-2t)
-------------r
2(1 - e-2t)
2
iXn.i I 1(1 - e-2^ Xn.i
exp —
R
*
2 V 1 + e-2t
R
*
1 - e-2‘ V3/2
1 + e-2t
. (38)
Taking into account (36), (37), and (38), for the square of the matrix element (35), we have:
|M |2
2-2п-5/2вһХп,і
a2 R*9 ahR 0
Jз (Xn.i) J5 (Xfn.i
г (1 + f-
[nh ]2 г
dt exp
-t ( -Vh( - fifth —+ i4er* +
-1
~ж
I И nh+0 - *( nh -1' -1
i
ln
4 sinh
" sinh
VLA
VkT/
2"
VTA
VkT
3
+ 2
-2t \ 3/2
(1 - exp(-t))-3/2exp (-2 (2+И (XRf)) (тг^)
2
. (39)
Taking into account (31) and (39) for the spectral intensity of recombination radiation (SIRR) in QD (30), we can write:
2
2
Ф(Х,Т) = Фо x
аһвһХ n.iX 2
R* 12 R0
J 3 (Xn.i)
2 V . J
г (1 + П
3/2
x j duP (u) x
(nh )2 г( ^v1)
f К f+0 - *( nh - 2> -1
dt exp
х*2в , , кТв
-t -%e - eeh —+ i4er0 +
Eh
ln
4 sinh
(1-exp(-t))-3/2exp(-2(1+^)fe)) (f^)
---"^=1 sinh
vLA\fkT J
1- e-2‘ \3/2
2"
VTA'fkT
2
x
3
+ 2
Г*2 + У xn,i + kT in
to + ( r2 + Eh ln
4 sinh
iVkT
sinh
( 2Я \vTAVkT
- nlh - X
2
(40)
2
2
i
г
where X = frw/Eh; Ф0 = /єе2 |Pehe0| /4п2 h3c3m0.
Figure 2 shows the SIRR dependence on the magnitude of the external electric field E0. It can be seen that the decrease in the SIRR value with increasing of E0 is accompanied by “dips”, that appear at certain values of the external electric field strength and temperature. In [5], it is shown that variation of the strength of the external electric field can lead to transformation of the shape of the double-well oscillatory potential, which simulates the system “QD - surrounding matrix”, while the transition to the symmetric shape of the double-well oscillatory potential is accompanied by the appearance of a peak in the field dependence of tunneling probability. Thus, the nature of the dip appears to be related to the effect of resonant tunneling, when the double-well oscillator potential becomes symmetric. An increase in the SIRR value with the temperature increasing (Fig. 2) is associated with an increase in the overlap integral of the wave functions of the initial and final states due to temperature smearing of energy levels. It should be noted that the presence of dissipative tunneling makes the optics of quantum dots very sensitive to the parameters of the surrounding matrix, which determine, respectively, the frequency of the phonon mode єL and the constant of interaction with the contact medium (with the heat-bath) єс. With an increase in the value of eL, the wave function of the A+ -state “spreads” due to the hole-phonon interaction, which is accompanied by a decrease in the SIRR value. An increase in the parameter єс leads to an increase in the “viscosity” of the surrounding matrix, i.e. to a decrease
Fig. 2. The SIRR dependence on the inverse temperature 1000/T and the external electric field strength E0 at RQ = 1; UQ = 350; n = 3.
in the probability of dissipative tunneling. As a result, the binding energy of the A+-state increases, and the overlap integral of the wave functions of the initial and final states decreases, which leads to a decrease in the SIRR value.
4. Conclusion
The possibility of the existence of a luminescence temperature quenching channel in semiconductor quasi-zerodimensional structures with “A+ + e” impurity complexes associated with dissipative tunneling of a hole into the matrix, surrounding the QD, has been theoretically substantiated. It is shown that this channel “triggers” efficiently at the temperature and strength of the external electric field, for which the double-well oscillatory potential, simulating the “quantum dot - surrounding matrix” system, becomes symmetric (resonant tunneling effect).
Acknowledgements
The present study was supported by the Ministry of Education and Science of the Russian Federation (Project No. 0748-2020-0012).
References
[1] Shamirzaev T.S., Klochikhin A.A, et. al. Photoluminescence of germanium quantum dots grown in silicon on a SiO2 monolayer. Physics of the Solid State, 2005, 47 (1), P. 82-85.
[2] Reznitsky A.N., Seksenbaev K.S, Permogorov S.A. Temperature dependence of the photoluminescence intensity of self-organized CdTe quantum dots in a ZnTe matrix under different excitation conditions. Physics of the Solid State, 2012, 54 (1), P. 123-133.
[3] Krevchik V.D., Levashov A.V. Energy spectrum of the complex A+ + e in a quantum dot in the adiabatic approximation. Physics of the Solid State, 2006, 48 (3), P. 589-592.
[4] Kusmartsev F.V., Krevchik V.D, et al. Phonon assisted resonant tunnelling and its phonons control. JETP Letters, 2016, 104, P. 392-397.
[5] Leggett A.J., Krevchik V.D, et al. Controllable dissipative tunneling. Tunnel transport in low-dimensional systems. Fizmatlit, Moscow, 20112012.496 p.
[6] Vainshtein I.A., Zatsepin A.F, Kortov V.S. Applicability of the empirical Varshni relation for the temperature dependence of the width of the band gap. Physics of the Solid State, 1999, 41 (6), P. 905-908.
[7] Ridley B. Quantum processes in semiconductors. Mir, Moscow, 1986.
[8] Zhou Hai-Yang, Gu Shi-Wei, Shi Yao-Ming. Electronic and shallow impurity states in semiconductor heterostructures under an applied electric field Commun. Theor. Phys., 2005, 44, P. 375-380.
[9] Lifshits I.M., Slezov V.V. On the kinetics of diffusion decomposition of supersaturated solid solutions. JETP: Journal of Experimental and Theoretical Physics (ZhETF), 1958, 35 (2), P. 479-492.