UDC 538.958
V G. Tyuterev, S. V Obukhov
ON THE RADIATIONLESS DECAY OF A DIRECT EXCITON IN AN INDIRECT GAP SEMICONDUCTOR
The theory of a radiationless dissociation of exciton in a semiconductor with the indirect forbidden gap is developed. It is shown that in GaAs. GaP and Ge under a pressure the main mechanism of exciton decay is the short-wave phonon-assisted scattering inside the conduction band.
Key words: semiconductors, exciton, electron-phonon interaction, density functional method.
1. Introduction
In several optical experiments the discrete levels of a direct exciton in the diamond and zincblende structure semiconductors are observed at the background of the uniform spectrum of single-electronic states [1-3]. In these cases they turn out to be resonance levels with finite lifetimes and hence the radiationless dissociation of a direct exciton becames possible. This is just a case of direct exciton in the indirect band gap semiconductor GaP [1]. In the direct band gap GaAs the lifetime of direct exciton sharply decreases under the external uniform pressure under the threshold where GaAs became the indirect one [2] . In indirect Ge the lifetime of direct exciton strongly depends on the external uniform pressure [3]. The decay mechanism of exciton’s state in all these cases is invoked by emission or absorption of short-wave phonons with large wave-vectors close to the crystal Brillouine zone edges.
A theory of exciton-phonon interaction is developed in the book [4] but it concerns to the long-wave phonons only. We present here a generalization of this theory [4] for the case of the exciton interaction with short-wave phonons.
2. General theory of the radiationless exciton decay
A multi-particle wave function of a free non-interacting electron-hole pair according to [4] can be written as an antisymmetrized product of single-particle states [4]
0™'(k h. k e ,{x}) =
AVvM (X1 ) ^vklP (X2 ) Vvk2a (X3 ) Vvk2b (X4 ) ...
...Vvkho (X>n) Vcke(r' (X>m)-.
•"V vk N a (X N )v vk nP(X N ).
Bloch states. Band index 1 has values (c,v) for conduction and valence bands correspondingly; s is a spin variable; a spin projection o has two possible values a
and p. The operator A means antisymmetrization by arguments X1,X2,.-,XN. The removal of this valence electron is equivalent to the creation of a hole with the spin (- o) and a wave vector (-k h). Certainly due to the double occupation of energy levels in the multiparticle ground state the single electronic state Vvk,-afe ) is still occupied by another valence electron in the valence band with the same vector kh but with the opposite spin projection (- o).
Following to [4] let us introduce the basis of elec-tron-hole excitations in the “excitonic” representation:
®"'(-,R l , R}):
I
e-kRL O'
'K - K, k, {y).
Here K is the wave vector of the electron-hole couple as a whole. The center of gravity of the couple is placed into the lattice vector RL. By definition k h = k - K and k e = k so K = k e - k h. In what follows the indices c,v are omitted. In addition the electron-phonon interaction can not change a spin so a = a' and the spin indices also can be omitted.
The bounded i-th state of an exciton which moves with a quasi-momentum K can be presented as a decomposition onto the basis of the states of virtual elec-tron-hole pairs.
k h . R L .
R})5.
The weights UK ( R L ) of virtual electron-hole pair 0(ke - kh.RL.{£,}) inside the bounded state can be
h? L ’
found from the equation [4]
Here one of the valence electron from a singleelectron state Yvk CT(with a wave-vector k h and spin o is removed from a multi-particle crystal ground state and instead another one in a conduc- = E (K )UK (R ).
tion state V ckea'(^m) is added with a wave-vector ke aind e spin o'. Here £ = (r,s);
V^k yin ) = ^k (r )Xa{s) are the single-electronic
Z(®(K.Rl) H„ ®(K,RlO)UK(Rl ) =
Here H0 is the multi-particle Hamiltonian of electron subsystem in the equilibrium lattice. Et (K) is
L
the energy of an exciton with a wave vector K . For the Wannier exciton the envelope approximation is justified and it is shown in [4] that
TH
ul (R,) = expo---------KR,)F„m (R,).
m„ + m.
The envelope function Fnlm (RL ) is equal to the hydrogenic one tynlm (p) at lattice vectors
Fnlm (RL ) = %lm (p)L
P=Ri
Here.p = (mere + mhrh)/( + mh); i;.rh are the coordinates of a conduction electron and a hole in the real space. m
the energy of a bound state one has
e , mh are their effective masses. For
Enlm (K) Eband +
*4
+Eg -
+
h2 k 2
2n2h2 2 (me + mh )
Here ^ 1 = mel + m, 1
*2
h ’
= e2 / a 0; s0 is a dielectric constant of a crystal. Eband is the sum of valence band energies. (-^e*4/in2%2) is the discrete energy levels of bounded exciton state. The last term is the kinetic energy. The normalization condition for
the envelope is £ |Fs (Rl )| = 1. Next it is conve-
Rl
nient to rewrite the multi-particle wave functions in the form [4]:
f \
m„
K
Fnlm (Q) = 4; Z Fnlm (Ri )
VN Rl
-iQRL
i,q
i,q
V (K,{^}) = ZFnlm k--
Rl
xO(k - K,k,{^}).
Here Fnlm (Q ) irier transformation of the
envelope
The normalization condition for the envelope is reduced to £ \Fnm (Q)|2 = 1
QeBZ
3. Exciton-phonon interaction in indirect gap band structure
For the Wannier exciton due to weak link between the electron and the hole one can treat their interaction with phonons independently. Hence the multi-particle operator of the electron-phonon interaction is:
^=S8h (X )=£«*, (r.).
Here &i (r ) is a single-particle operator of interaction of phonon with the wave vector q and the i-th electron. We are looking for the dissociation that is for
a decay of exciton to the uncoupled electron and hole pair. The initial state here is a ground s-state of exciton
¥s (K. with its energy Es (K). The final state is
uncoupled state of electron and hole with wave vectors ke and kh correspondingly. This last state is described
by a pair wave function 0( \, k ' e, .The energy
of Ms state is Efree = Eband + ec (ke) - sv (k;) ± .
The energy conservation rule for the phonon-assisted dissociation should be written as:
Es (K) = Eband + SC (ke ) -Sv (K ) ± .
Here is the phonon energy. The matrix element of the phonon-assisted transition from ¥, (K. {!;} to ®(k (,k'„,{);} is
M (k'h. k',. K ) = J>T, (K .{X})x
x^siq, (X, )$(k'h. k 'e .{X}) d NX.
i.q ' '
At absolute temperatures close to zero we can neglect the kinetic energy of exciton and take K = 0. The calculation reduces to the calculation of the matrix elements between the virtual electron-hole (e - h)
states Ojkh.ke.{&,} with “vertical” excitation
kh = ke = k which enters into the bounded states with a weight Fnm (k) and the final state of free pair
®(k (, k 'e ,{&,}. Keeping in mind that J.I J. ,dr1dr2...drN . using the nor-
sls2---sN
malization condition and the momentum conservation law in the scattering of Bloch electron at the phonon
{V*nk (r)^h±q (r)Vn'k'(r)dr =
= 5k'.k±q j V*nk (r)^h±q (r) Vn'k±q (r)dr
one gets after some algebra
M (k h. k e. K ) = 4. (k h. k e. K) + Acc (k h. k e. K).
Here
Av (kh. ke. K) = Sk kSl>k±q X
xF,(k. (r)8h,qV„,(r)dr
is the matrix element of a transition with the emission of phonon with the wave vector q between the virtual (e-h) pair with “vertical” excitation of electron having the quasi-momentum ke = kh = k and the uncoupled free electron with the wave vector k' = k and the
hole with kh = k ± q. For a Wannier exciton the weight Fs (k) is a delta-shaped function and so the lengths of ke and k are restricted by values k| << 2n / a . It is easy to find by inspection (Fig. 1) that the configurations of a valence band in the indirect semiconductors do not give an opportunity to find a wave vector k'h in order to satisfy both energy and momentum conservation rules. So (h-h) scattering inside a valence band is impossible in this configuration. A second term
Acc ( k h, k e, K ) = 5khk 5
Here r means that the states belong to the center of a Brillouin zone and it is used that energy in T-point of a conduction band is s cr = s,,
Jv r + EG .
k^k kek ±q
XFs (kh )jv*vke (r)5h± q Vvkh (r) dr
is a matrix element of transition assisted by phonon with a wave vector q between a virtual “vertically” ke = kh = k excited (e-h) couple and a free hole with kh = k and conduction electron with a wave vector k e = k ± q . In the indirect semiconductor this process appears to be possible (Fig. 1) and it involves a shortwave phonon. Consequently the total probability includes only the scattering inside the conduction band and it is assisted by shortwave phonons. Phonon energies %raq. exciton’s levels and its kinetic energy are small in comparison with a forbidden gap. Neglecting them one gets a energy conservation condition in the form ec (ke) - sv (k'h) = EG. The total probability is written as:
W = Zl Fs (k )|2 |j Vlk±q (r )5h±q (r ) Vck (r ) dr X
kq
X5(sc (k ± q )-sv (k)-Eg ).
Let us use the integral mean value theorem. then
W = (ZI F (k )|2 )x XZ| j vl-±q ( r )§h±q ( r ) Vck (r) dr\2 X
q
XS(s c ( k ± q )-s v (k )-Eg ).
Due to the delta shape of Fs (k) the mean value k
i |2
is close to zero. Keeping in mind Zl F (k ) = 1
k
and using the Fermi golden rule one gets the probability per second of the exciton dissociation
W = Z | j Vc±q (r ) 5h±q (r ) Vcr (r ) dr
q
X5(sc (±q )-scr).
Fig. 1. Scheme of the direct exciton decay in the indirect gap semiconductor
4. Results and conclusion
The last expression means that the problem of the calculation of the decay probability of a direct exciton in the indirect gap material reduces to the calculation of the “horizontal” short-wave phonon-assisted transition from a single-electronic Bloch state in the centre of a Brillouin zone into aside minima (Fig.1). This theory was applied in [5, 6] to calculate the exciton decay in GaP and GaAs. Electron and phonon spectra were calculated in the density functional perturbation theory [7]. According this our calculation the line-width of exciton in GaP is 9.2 meV in a reasonable agreement with experimental values 16 meV [2] and 12 meV [8]. The calculated life time of direct exciton in GaAs t = 1.5±0.2 ps [5, 6] also is in a good agreement with experimental 1.4 ps [9], 2.0±0.5 ps [10]. The calculated pressure dependence [11] of exciton lifetimes and the width of exciton spectral lines in Ge under the external pressure also agree satisfactorily with an existing experiment [3]. We believe, that the method developed in [5, 6, 11] for the calculating of electron-phonon scattering probabilities can be useful in the investigation of relaxation of highly excited electrons [12] and the problem of high-temperature superconductivity [13].
References
1. Goni A. R., Cantarero A., Syassen K., Cardona M. // Phys. Rev. B. V. 41, 10111 (1990).
2. Zollner S., Kircher J., Cardona M., Gopalan S. // Solid-State Electron. V. 32, 1585 (1989).
3. Li G. H., Goni A. R., Syassen K., Cardona M. // Phys.Rev. B. V. 49, 8017 (1994).
4. Knox R. S. Theory of excitons. N.Y., L., Academic Press, 1963.
5. Sjakste J., Vast N., Tyuterev V. // Phys. Rev. Lett. V. 99, 236405 (2007).
6. Sjakste J., Tyuterev V., Vast N. // Phys. Rev. B. V. 74, 235216 (2006).
7. Baroni S., Gironcoli S. de, Corso A. D., Giannozzi P. // Rev. Mod. Phys. V. 73, 515 (2001).
8. Sell D. D., Lawaetz P. // Phys. Rev. Lett. 26, 311 (1971).
9. Young J., Gong T., Faucher P., Kelly P. // Phys.Rev. B. V. 50, 2208 (1994).
10. Shah J., Deveaud B., Damen T C. et al. // Phys. Rev. Lett. V. 59, 2222 (1987).
11. Tyuterev V., Obukhov S. V., Vast N., Sjakste J. // Phys. Rev. B. V. 84, 035201 (2011).
12. Tyuterev V. G. // Tomsk State Pedagogical University Bulletin. 2006. Issue 6. P. 7-11.
13. Manca P., Mula G., Kuz'menko A. B. et al. // Tomsk State Pedagogical University Bulletin. 2002. Issue 2(30). P. 3-10.
Tyuterev V. G.
Tomsk State Pedagogical University.
Ul. Kievskaya. 60. Tomsk. Russia. 634061.
E-mail: [email protected]
Obukhov S. V.
Tomsk State Pedagogical University.
Ul. Kievskaya. 60. Tomsk. Russia. 634061.
E-mail: [email protected]
Received 14.03.2011.
В. Г. Тютерев, С. В. Обухов БЕЗИЗЛУЧАТЕЛЬНЫЙ РАСПАД ПРЯМОГО эКСИТОНА В НЕПРЯМОЗОННОМ ПОЛУПРОВОДНИКЕ
Построена теория безизлучательной диссоциации экситона в полупроводнике с непрямой запрещенной зоной. Показано, что в GaAs, GaP, а также в Ge под давлением основным механизмом распада экситона является рассеяние электрона на коротковолновых фононах в зоне проводимости.
Ключевые слова: полупроводники, экситон, электрон-фононное взаимодействие, метод функционала плотности.
Тютерев В. Г., докотор физико-математических наук, професор, зав. кафедрой.
Томский государственный педагогический университет.
Ул. Киевская, 60, Томск, Россия, 634061.
E_mail: [email protected]
Обухов С. В., научный сотрудник отдела исследований и разработок.
Томский государственный педагогический университет.
Ул. Киевская, 60, Томск, Россия, 634061.
E-mail: [email protected]