Feb 27, 2004 - Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, USA. Received 3 September 2003; ...
PHYSICAL REVIEW B 69, 075214 共2004兲
High-field transport of electrons and radiative effects using coupled force-balance and Fokker-Planck equations beyond the relaxation-time approximation Danhong Huang, T. Apostolova, P. M. Alsing, and D. A. Cardimona Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, USA 共Received 3 September 2003; revised manuscript received 17 October 2003; published 27 February 2004兲 The dynamics of a many-electron system under both dc and infrared fields is separated into a center-of-mass and a relative motion. The first-order force-balance equation is employed for the slow center-of-mass motion of electrons, and the Fokker-Planck equation is used for the ultrafast relative scattering motion of degenerate electrons. This approach allows us to include the anisotropic energy-relaxation process which has been neglected in the energy-balance equation in the past. It also leads us to include the anisotropic coupling to the incident infrared field with different polarizations. Based on this model, the transport of electrons is explored under strong dc and infrared fields by going beyond the relaxation-time approximation. The anisotropic dependence of the electron distribution function on the parallel and perpendicular kinetic energies of electrons is displayed with respect to the dc field direction, and the effect of anisotropic coupling to an incident infrared field with polarizations parallel and perpendicular to the applied dc electric field is shown. The heating of electrons is more accurately described beyond the energy-balance equation with the inclusion of an anisotropic coupling to the infrared field. The drift velocity of electrons is found to increase with the amplitude of the infrared field due to a suppressed momentum-relaxation process 共or frictional force兲 under parallel polarization but decreases with the amplitude due to an enhanced momentum-relaxation process under perpendicular polarization. DOI: 10.1103/PhysRevB.69.075214
PACS number共s兲: 72.20.Ht, 61.82.Fk, 61.80.Az, 72.10.Di
I. INTRODUCTION
High-field transport of electrons is usually referred to when the electric field applied to the system is no longer small and the current-voltage relation deviates from linearresponse theory. There are many theories that have been proposed to describe this phenomenon since 1930. However, the semiclassical Boltzmann transport equation with a drift term still seems to be the only one amenable to practical calculations. In practice, however, solution of the Boltzmann transport equation also becomes a laborious task when the system goes beyond the linear regime. Fro¨hlich and Paranjape1 first utilized a displaced Maxwellian distribution function to describe the electron transport in insulators and semiconductors. Later, Arai2 applied a similar model to describe electron transport in metals at low temperatures but by means of a Fermi-Dirac distribution function, in which a finite electron temperature was predicted under an applied field even when the lattice temperature became zero. Shortly after that, Lei and Ting3 proposed the coupled first-order force-balance and energy-balance equations to describe electrons in semiconductors and metals at both low and high temperatures by assuming an isotropic quasi-thermal-equilibrium 共FermiDirac兲 distribution for hot electrons with a temperature different from the lattice temperature. The effect of acceleration of electrons under a dc field can be included in the Boltzmann transport equation using a field-induced drift term after an acceleration ansatz has been employed for uncorrelated scatters.4 This approach was recently generalized to study the drifting motion of electrons under a low-frequency ac bias.5 There appears to be much less justification for treating electron transport with an ac field oscillating at optical frequencies by directly adding a drift term to the Boltzmann transport equation even though the resulting conductivity is inversely proportional to the 0163-1829/0/69共7兲/075214共12兲/$22.50
square of the optical frequency, since the time-average velocity of an electron under an optical field is in fact zero between two successive scattering events. For weak external fields, linear-response theory6 or a linear-field expansion of the Boltzmann transport equation4 could be applied. In this weak-field limit, a complete approach for treating electronphonon drag effects of any degree under a magnetic field was proposed by simultaneously solving Boltzmann equations for both electrons and phonons.7 However, these theories fail to produce correct results for high-field electron transport. It has been known for quite a long time that the drift of electrons under a dc field can be regarded as a field-driven center-of-mass 共collective兲 motion of many electrons.8 The scattering of electrons by the lattice or by impurities can then be considered as a motion relative to this center of mass. A uniform external optical field is also known to couple only to the center-of-mass motion. The ions on the lattice and ionized impurity atoms in the material system remain stationary with respect to the moving center of mass. Therefore, the center-of-mass motion couples to the relative scattering motion. This causes a frictional force to act on each drifting electron due to lattice or impurities.3 As a result, the acceleration of electrons under the dc field will be opposed by this friction. The classically measurable frictional force is a quantum-statistical average of the friction acting on each electron, involving a many-electron relative scattering motion, and therefore depends on the occupation of electrons in different states. On the other hand, the distribution of electrons in various states is determined by the relative elastic and inelastic scattering motions of electrons, including Coulomb scattering, phonon scattering, and impurity scattering. The existence of a dc or an optical field couples the center-of-mass motion of electrons with their relative scattering motion. In this situation, the dc drift modifies the energy conservation in electron scattering with phonons by adding a
69 075214-1
PHYSICAL REVIEW B 69, 075214 共2004兲
HUANG et al.
Doppler shift along the dc field direction, causing the scattering process to be anisotropic.3 The time-dependent drift, on the other hand, anisotropically modifies the strength of the electron scattering as well as the energy conservation at the same time.9,10 For a spatially uniform optical field incident on semiconductors, direct photon absorption is forbidden when the photon energy is smaller than the band gap. However, the incident photons can still be absorbed by electrons with the help of their scattering with phonons, which provides momentum to the electrons. The simplest way to treat a steady-state distribution of electrons is the introduction of an electron temperature for a Fermi-Dirac function of hot electrons. The value of the electron temperature, different from the lattice temperature, can be obtained by employing a first-order energy-balance equation for the relative motion of electrons.3 However, the assumption of an isotropic Fermi-Dirac function for hot electrons cannot be justified even in steady state when the external field exists for all time. In the limits of strong optical field and zero drift velocity, the exact solution of the Boltzmann scattering equation including the electronelectron interaction shows the existence of multiple side peaks in the electron distribution as a function of electron kinetic energy.10 Thus, the electron distribution is nonequilibrium and the energy-balance formalism fails.11 After the optical field has been turned off, the electron distribution quickly relaxes into a quasiequilibrium Fermi-Dirac function through electron-electron interaction, bringing the energybalance equation back into validity. The energy-balance equation cannot be justified when either the drift velocity is large 共the Doppler shift is comparable to the phonon energy兲 or the optical field is strong. In addition, the energy-balance equation excludes the effect of electron-electron scattering on the electron temperature. A more rigorous method involves the use of the Boltzmann scattering equation without a drift term for anisotropic scattering of electrons in relative motion. In the presence of an intense optical field only, this approach has been successfully applied to study laser damage in dielectric and semiconductor materials.9,10,12 However, the incoherent Boltzmann scattering theory cannot be applied to study the dynamics of electrons if a dipole coupling of electrons to the external optical field exists when either the photon energy becomes greater than the band gap or a spatially nonuniform optical field is incident on the system after diffraction by a surface grating. In this case, a full coherent density-matrix theory13 with off-diagonal matrix elements is required. One advantage of the current approach is that it includes a generalization of the definition of electron temperature for non-equilibrium electrons. In addition, the extension of the current approach to one that incorporates coherent density-matrix equations13 for the relative motion of electrons is straight forward when pump-laser-induced interband excitation exists. In this paper, we use a first-order force-balance equation to describe electron drift due to a slow center-of-mass motion under a dc electric field. At the same time, we adopt the Boltzmann scattering equation without a drift term for the ultrafast process of anisotropic scattering of electrons with phonons and impurities in a relative motion by assuming that
no dipole coupling to the external optical field exists in the system. In order to highlight the physics and simplify the calculation of high-field electron transport in the presence of both dc and optical fields, we expand the Boltzmann scattering equation up to second order for degenerate electrons to obtain the Fokker-Planck equation including Joule heating and antidiffusion terms.12 In this way, the anisotropic distribution of electrons and the anisotropic scattering of electrons under an incident optical field with different polarizations can be fully incorporated into our model. In addition, the heating of electrons can be more precisely described beyond the isotropic Fermi-Dirac function of hot electrons used in a first-order energy-balance equation. The organization of the paper is as follows. In Sec. II, we introduce our model and derive the coupled force-balance and Boltzmann scattering equations by separating the centerof-mass and relative motions of electrons. After the expansion of the Boltzmann scattering equation we obtain the Fokker-Planck equation to simplify the calculation of highfield electron transport beyond the relaxation-time approximation. Analytical expressions for the expansion coefficients are presented in Sec. III. Numerical results are displayed in Sec. IV for the dependence of anisotropic distributions of electrons on the parallel and perpendicular kinetic energies with respect to the dc field direction, as well as for the dependence of drift velocities and heating of electrons on polarizations of the incident optical field parallel and perpendicular to an applied dc electric field. The paper is briefly concluded in Sec. V.
II. MODEL
Let us consider a bulk semiconductor, e.g., GaAs, which is doped with electron concentration 3D to form a threedimensional 共3D兲 electron gas. In this section, we first separate the dynamics of a many-electron system into a centerof-mass motion plus a relative motion under both dc and infrared fields. The equation for the center-of-mass motion of electrons is built after a quantum-statistical average is taken to obtain a classical frictional force. The relative motion of electrons is studied by using the Boltzmann scattering equation including anisotropic scattering of electrons with phonons and impurities beyond the relaxation-time approximation. After that, the Fokker-Planck equation for degenerate electrons including Joule heating and antidiffusion terms is derived by expanding the Boltzmann scattering equation to second order, which significantly simplifies the calculation of high-field electron transport.
A. Center-of-mass and relative motions
In the presence of a spatially uniform infrared field with a ជ (t) and a dc field Eជ dc , the Hamiltonian vector potential A of an interacting three-dimensional electron gas can be written as
075214-2
PHYSICAL REVIEW B 69, 075214 共2004兲
HIGH-FIELD TRANSPORT OF ELECTRONS AND . . .
1
H⫽
2m * ⫺e ⫺
兺i
ˆ ជ 兲2⫹ 共 pជ i ⫺eA
兺 i⬍ j
e2 4 ⑀ 0 ⑀ r兩 rជ i ⫺rជ j 兩
ជ a兲 U imp共 rជ i ⫺R 兺i rជ i •Eជ dc⫹ 兺 i,a
ជrជ U ion共 rជ i ⫺Rជ ᐉ 兲 , uជ ᐉ •ⵜ 兺 i,ᐉ
共1兲
i
where i⫽1,2, . . . ,N e is the index of N e electrons, a ⫽1,2, . . . ,N a , is the index for N a impurity atoms, ᐉ ⫽1,2, . . . ,N ᐉ is the index for N ᐉ lattice ions, rជ i is the posiជ a and Rជ ᐉ are the position tion vector for the ith electron, R vectors of impurity atoms and lattice ions, uជ ᐉ represents the ion displacement from the thermal equilibrium position, m * is the effective mass of electrons, and ⑀ r is the relative dielectric constant of the host material. The single-electron ជrជ , and both the impurity momentum operator is pជˆ i ⫽⫺iបⵜ
† tion兲 operator of electrons and bˆ qជ (bˆ qជ ) to denote the creation 共annihilation兲 operator of phonons. U imp(q) is the Fourier transform of the impurity potential and C q is the electron-phonon coupling constant 共given in Sec. II B兲. The existence of the infrared field causes the ions in the lattice 共or ionized impurity atoms兲 to oscillate with time relative to the electrons. This driven oscillation transfers kinetic energy from the lattice to the electrons through the electron-phonon coupling. The quantum mechanics for the electron system gives rise to the following operator commutations14 with N eⰇ1:
ˆ ជ c兴 ⫽ 关 Rជ c, ˆ pជ ⬘ i 兴 ⫽0, 关 rជ i ⬘ , P 关 R ␣c , Pˆ c 兴 ⫽iប ␦ ␣ ,
冋
关 r i⬘␣ ,pˆ ⬘j  兴 ⫽iប ␦ ␣ ␦ i j ⫺
i
ជ a) and the ion potential U ion(rជ i ⫺Rជ ᐉ ) are potential U (rជ i ⫺R included. We first define the center-of-mass momentum and position vectors by imp
ˆ ជ c⫽ P
Ne
兺
i⫽1
pជˆ i ,
Rជ c⫽
1 Ne
Ne
兺 rជ i ,
i⫽1
ជ c. rជ i ⬘ ⫽rជ i ⫺R
共2兲
d ˆc 1 ˆc ˆ rel兴 ជ ⫽ 关P ជ ,Hcm⫹H P dt iប
ជ dc⫺i ⫽N eeE
共3兲
⫺i
By using the center-of-mass and relative momentum and position vectors defined above, we can separate the total Hamiltonian, including the Hamiltonians of electrons and phonons, into one center-of-mass Hamiltonian Hcm and anˆ rel , given by other relative Hamiltonian H Hcm⫽
ˆ rel⫽ H
ˆ ជ c兲 2 共P 2N em *
⫹
N ee 2 A 2 2m *
⫺
e m*
ជ •ˆ ជ c⫺N eeEជ dc•Rជ c, A P
共4兲
† † k aˆ kជ aˆ kជ ⫹ 兺 ប q bˆ qជ bˆ qជ 兺 kជ , qជ ,
⫹ ⫹ ⫹
1 2
兺 兺qជ
kជ ,kជ ⬘ , , ⬘
c
ជ ជ ជ U imp共 q 兲 qជ e iq •(R ⫺R ) ˆ qជ , 兺 qជ ,a c
N em *
c
uជ d⬅
where the infrared field is treated classically, V is the volume of the system, ប q is the phonon energy with wave number q for mode 共totally three modes兲, k ⫽ប 2 k 2 /2m * is the kinetic energy of electrons with wave number k, and the index ⫽⫾1 is for up and down spin states for electrons. † Here, we use aˆ kជ (aˆ kជ ) to represent the creation 共annihila-
共7兲
共8兲
† where ˆ qជ ⫽ 兺 kជ , aˆ kជ ⫹qជ aˆ kជ is the density-wave operator for ជ (t)⫽(Eជ op /⍀ f)cos(⍀ft), ជ op(t) and E electrons, A ជ ⫽E opsin(⍀ft). Here, ⍀ f is the infrared-field angular frequency and E op is its amplitude. Applying a quantumstatistical average 具具 ••• 典典 to both Eqs. 共7兲 and 共8兲, we define the following quantities:
† † aˆ ជ ជ aˆ ជ ជ aˆ kជ aˆ kជ ⑀ 0 ⑀ rq 2 V k ⫹q k ⬘ ⫺q ⬘ ⬘ ⬘
共5兲
a
ˆ ជc P e ˆ d ជc 1 ជc ˆ rel兴 ⫽ ជ, ⫽ 关 R ,Hcm⫹H ⫺ A uជ ⫽ R dt iប N em * m *
ជ ជ † † C q 共 bˆ qជ ⫹bˆ ⫺qជ 兲 e iq •R aˆ kជ ⫹qជ aˆ kជ 兺 兺 kជ , qជ ,
兺 兺 kជ , qជ ,a
ជ ជ † C q qជ e iq •R 共 bˆ qជ ⫹bˆ ⫺qជ 兲 ˆ qជ 兺 ជq ,
e2
ជ ជc ជ † U imp共 q 兲 e iq •(R ⫺R a) aˆ kជ ⫹qជ aˆ kជ ,
共6兲
where ␣ ,  ⫽x,y,z. By using these operator commutations, we get two Heisenberg equations for the center-of-mass motion of electrons:
and those for the relative motion by 1 ˆc ˆ ˆ ជ , pជ ⬘ i ⫽pជ i ⫺ P Ne
册
1 ⬇iប ␦ ␣ ␦ i j , Ne
d uជ ⬅ dt 0
冓冓 冔冔 d c Rជ dt
冓冓 冔冔
⫽uជ 0 ⫺
d ˆc ជ P dt e
,
ជ op cos共 ⍀ ft 兲 , E m *⍀ f
共9兲
共10兲
ˆ ជ c/N em * 典典 is the part of the drift velocity rewhere uជ 0 ⫽ 具具 P lated to the center-of-mass canonical momentum. The center-of-mass motion of electrons can be regarded as a slow motion in comparison with the fast relative scattering motion of electrons. Combining Eqs. 共7兲–共10兲 and using the Fermi 共Bose兲 statistics for electron 共phonon兲 operators † † aˆ kជ ,aˆ kជ (bˆ qជ ,bˆ qជ ), we get the force-balance equation3
075214-3
PHYSICAL REVIEW B 69, 075214 共2004兲
HUANG et al.
d d ជ op共 t 兲 , uជ d⫽N em * uជ 0 ⫹N eeE dt dt
共11兲
d ជ dc⫹Fជ i关 uជ d兴 ⫹Fជ p关 uជ d兴 , uជ ⫽N eeE dt 0
共12兲
N em * N em *
where Fជ i (Fជ p ) are the frictional forces due to impurity 共phonon兲 scattering. For steady state or cases in which the time period of an optical field is much smaller than the momentum-relaxation time of electrons, the last term in Eq. 共11兲 has no contribution. By considering the slow drifting motion of the center-of-mass in comparison with the ultrafast energy-relaxation process of electrons in the relative motion and keeping only the leading-order interaction of electrons with phonons and impurities, we can introduce a quantumstatistical average11 to the frictional forces acting on all electrons which occur on the right-hand side of Eq. 共12兲 due to impurity and phonon scattering. We know that the electron temperature should be determined by the ultrafast relative 共internal兲 motion of electrons. In contrast, the drifting motion of electrons is very slow compared to the relative scattering process. As a result, when the frictional forces acting on all electrons are quantum-statistically averaged, one can use the Fermi-Dirac function as a first-order approximation with an effective temperature adiabatically evolving from T * (t⫽ ⫺⬁)⫽T to T * (t⭓0)⫽T e . The anisotropic scattering in the relative motion of electrons produces nonzero frictional forces given by
ជ i关 uជ d兴 ⫽⫺2 n iV F ⬁
⫻
兺
M ⫽⫺⬁
兺kជ 兺qជ qជ 兩 U imp共 q 兲 兩 2共 n kជ ⫹qជ ⫺n kជ 兲 J 兩2M 兩 共 兩 Mqជ 兩 兲 ␦ 共 k⫹q ⫺ k ⫹បqជ •uជ 0
⫹M ប⍀ f兲 ,
ជ p关 uជ d兴 ⫽⫺4 F
兺kជ 兺 qជ ,
qជ 兩 C q 兩 2 共 n kជ ⫹qជ ⫺n kជ 兲
共13兲 ⬁
兺
M ⫽⫺⬁
J 兩2M 兩 共 兩 Mqជ 兩 兲
⫻关 N 0 共 q 兲 ⫺N 0⬘ 共 q ⫹qជ •uជ 0 ⫹M ⍀ f兲兴 ⫻ ␦ 共 k⫹q ⫺ k ⫹ប q ⫹បqជ •uជ 0 ⫹M ប⍀ f兲 .
共14兲
Here, M ⫽0,⫾1,⫾2, . . . , J M (x) is the M th-order Bessel function, n i is the impurity concentration, and N 0 ( q ) ⫽ 关 exp(បq /kBT)⫺1 兴 ⫺1 is the Bose-Einstein function for equilibrium phonons with lattice temperature T. The electron-phonon coupling strength 兩 C q 兩 2 and electronimpurity coupling strength 兩 U imp(q) 兩 2 have essentially been modified by the multiplicative factor J 兩2M 兩 ( 兩 Mqជ 兩 ). Using Eq. 共13兲, we get the same resistivity obtained from the Boltzmann transport equation to leading order for both nondegenerate and degenerate electron gases when impurity scattering is considered and E op⫽0. At high temperatures (T ⭓300 K), it is known that the scattering from impurities becomes much less significant than that from phonons, and we will ignore it hereafter. N 0⬘ ( q ⫹qជ •uជ 0 ⫹M ⍀ f) in Eq. 共14兲 is obtained from N 0 ( q ⫹qជ •uជ 0 ⫹M ⍀ f) by replacing the lattice temperature T with the nonequilibrium electron
temperature T e , which will be explicitly given3 in Sec. III. The force-balance equation describes the momentum dissipation process for drifting electrons. It is easy to show from Eq. ជ p关 uជ d兴 ⫽0 if uជ 0 ⫽0ជ . Here, the center-of共14兲 that Fជ i关 uជ d兴 ⫽F mass drift causes a nonzero frictional force by producing a Doppler shift qជ •uជ 0 relative to the phonon frequency q and introducing anisotropic phonon scattering. On the other hand, the infrared field mainly modifies the electron-phonon coupling strength and drives T e away from T by changing the average kinetic energy of the electrons. n kជ in Eq. 共14兲 is the occupation probability of electrons, and Mqជ ⫽
ជ op eqជ •E m * ⍀ f2
.
共15兲
B. Boltzmann and Fokker-Planck equations
As we noted earlier, both uniform infrared and dc fields couple only to the center-of-mass motion of electrons but not directly to their relative motions. The drift of electrons, however, couples the center-of-mass and relative motions. At high temperatures, phonon scattering in the system will dominate impurity scattering. The first-order Coulomb effects can be included within the Hartree-Fock approximation. The Hartree energy shifts only the band edges and is irrelevant to our current work. The Fock energy between two electrons can be neglected compared to the average kinetic 1/3 a B* Ⰷ1 with effective Bohr energy of each electron if 3D 2 2 * radius a B ⫽4 ⑀ 0 ⑀ rប /m * e . For the concentration 3D cho1/3 sen in this paper, 3D a B* Ⰷ1 is satisfied. The second-order Coulomb effects include screening and pair scattering. Screening has been considered in this paper within the static Thomas-Fermi approximation. A discussion about using a more sophisticated random-phase approximation can be found in a recent paper.15 Pair scattering does not directly contribute to the drift velocity. For time scales much longer than the pair scattering time, this effect can be included phenomenologically by introducing a homogeneous lifetime 1/⌫ ee to electrons, where ⌫ ee is the pair scattering rate on the order of 10 ps⫺1 . This will broaden the loss ( ␦ ) functions in both the force-balance and Boltzmann scattering equations. Obviously, this will not qualitatively change the features predicted in this paper. The force-balance equation3 for the center-of-mass motion takes into account the momentumdissipation effect. The energy-relaxation effect is included in the relative scattering motion. If we could assume an isotropic Fermi-Dirac function for hot electrons,3 the electron temperature would be found from a first-order energy-balance equation for the relative motion of electrons. However, there is no justification for this assumption if the external field exists for all time. On the other hand, the Boltzmann scattering equation has been very successfully applied to explain laser damage of dielectric and semiconductor materials10,12 under an intense laser field. We consider the situation in which the infrared photon energy is much smaller than the band gap of bulk GaAs and the amplitude of the infrared field is moderate. In this case, we can neglect the multiphoton interband excitation process. For the relative motion of
075214-4
PHYSICAL REVIEW B 69, 075214 共2004兲
HIGH-FIELD TRANSPORT OF ELECTRONS AND . . .
electrons under both dc and infrared fields, we use the Boltzmann scattering equation without the drift term to describe the relative scattering motion with phonons at high temperatures, given by d (in) (out) n ជ ⫽W kជ 共 1⫺n kជ 兲 ⫺W kជ n kជ , dt k
共16兲
where the scattering-in and scattering-out rates of electrons by phonons are12 (in) 2 W kជ ⫽ ប
⬁
兩 C q 兩 兺 J 兩2M 兩 共 兩 Mqជ 兩 兲关 n kជ ⫺qជ N qជ ␦ 共 k 兺 M ⫽⫺⬁ qជ , 2
⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫺M ប⍀ f兲
d d f 共 kជ 兲 ⫹ 关 V T共 kជ 兲 ⫹V F共 kជ 兲兴 f 共 kជ 兲 ⫺ 关 D T共 kជ 兲 ⫹D F共 kជ 兲兴 dt d kជ
⫹n kជ ⫹qជ 共 N qជ ⫹1 兲 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫹M ប⍀ f兲兴 , (out)
W kជ
⫽
2 ប
兺 qជ ,
共17兲 ⬁
兩 C q 兩 2
兺
M ⫽⫺⬁
J 兩2M 兩 共 兩 Mqជ 兩 兲关共 1⫺n kជ ⫹qជ 兲
⫻N qជ ␦ 共 k⫹q ⫺ k ⫺ប q ⫹បqជ •uជ 0 ⫺M ប⍀ f兲 ⫹ 共 1 ⫺n kជ ⫺qជ 兲共 N qជ ⫹1 兲 ␦ 共 k⫺q ⫺ k ⫹ប q ⫺បqជ •uជ 0 ⫹M ប⍀ f兲兴 .
共18兲
The terms containing N qជ ⫹1 represent the contributions from phonon emission, while the terms containing N qជ are the contributions from phonon absorption. In the absence of both dc and infrared fields, the steady-state solution of Eq. 共16兲 gives rise to a Fermi-Dirac function with the lattice temperature. The pair scattering effect can be phenomenologically included by changing the ␦ functions to Lorenzian functions with a homogeneous broadening ប⌫ ee , where the scattering rate ⌫ ee is on the order of 10 ps⫺1 . However, this does not qualitatively change the features predicted in this paper. The modification of electron-phonon scattering with M ⫽0 results from the drifting electrons interacting with the infrared field through phonon scattering. This effect is also
A F共 kជ 兲 ⫽
historically known as free-carrier absorption of photons. The free-carrier absorption of photons occurs only close to the Fermi surface if the Fermi energy of electrons is much larger than 共i.e., degenerate electrons兲 the electron temperature and phonon energies. In this case, by assuming small energies, e.g., ប q , ប⍀ f , បqជ •uជ 0 relative to the electron kinetic energy k ⬇ F with F⫽ប 2 (3 2 3D) 2/3/2m * at T⫽0 K for weak dc and moderate infrared fields, we can expand the Boltzmann scattering equation in Eq. 共16兲 to second order with respect to these energies. As a result, the Fokker-Planck equation for drifting electrons under dc and infrared fields formally becomes12
2ប
⫻
d2 2
d kជ
f 共 kជ 兲 ⫽ 关 A T共 kជ 兲 ⫹A F共 kជ 兲兴 f 共 kជ 兲 .
共19兲
Here, f ( kជ )⫽ ( k )n kជ is the continuous distribution function ( k ) of electrons with density of states ⫽ 冑(2m * ) 3 k /2 2 ប 3 . We note that this equation is different from the standard Fokker-Planck equation16 and cannot be written as a conservation law for particle number and current within the energy space even in the absence of source terms on the right-hand side of Eq. 共19兲. But we will still call Eq. 共19兲 the Fokker-Planck equation for simplicity. This equation is linear with respect to f ( kជ ) and, thus, greatly simplifies the calculation of the scattering dynamics of electrons. On the right-hand side of Eq. 共19兲, we include two source terms. The coefficient for the thermal spontaneous phonon emission is A T共 kជ 兲 ⫽
2 ប
兩 C q, 兩 2 关 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲 兺 qជ ,
⫺ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲兴 ,
共20兲
and the coefficient for the field-induced spontaneous phonon emission is
2 兩 C q, 兩 2 M qជ 关 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫺ប⍀ f兲 兺 ជq ,
⫺ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲 ⫺ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫺ប⍀ f兲兴 ⫺
ប
2 兩 C q, 兩 2 M qជ 关 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲 ⫺ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲兴 . 兺 ជq ,
共21兲
In addition, we define in Eq. 共19兲 the thermal energy transfer rate due to phonon scattering by V T共 kជ 兲 ⫽
2 ប
兩 C q, 兩 2 共 ប q ⫺បqជ •uជ 0 兲关 N qជ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲 ⫺ 共 N qជ ⫹1 兲 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲兴 , 兺 qជ ,
and the field-induced energy transfer rate by 075214-5
共22兲
PHYSICAL REVIEW B 69, 075214 共2004兲
HUANG et al.
V F共 kជ 兲 ⫽
2ប
2 兩 C q, 兩 2 M qជ 兵 N qជ 共 k ⫺ k⫺q 兲关 ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫺ប⍀ f兲兴 兺 ជq ,
⫹ 共 N qជ ⫹1 兲共 k ⫺ k⫹q 兲关 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫺ប⍀ f兲兴 其 ⫺
ប
2 兩 C q, 兩 2 M qជ 共 ប q ⫺បqជ •uជ 0 兲关 N qជ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲 ⫺ 共 N qជ ⫹1 兲 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲兴 . 兺 ជq ,
共23兲 Finally, we define in Eq. 共19兲 the thermal energy diffusion coefficient due to phonon scattering by D T共 kជ 兲 ⫽
ប
兩 C q, 兩 2 共 ប q ⫺បqជ •uជ 0 兲 2 关 N qជ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲 ⫹ 共 N qជ ⫹1 兲 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲兴 , 兺 qជ ,
共24兲
and the field-induced energy diffusion coefficient by D F共 kជ 兲 ⫽
4ប
2 兩 C q, 兩 2 M qជ 兵 N qជ 共 k ⫺ k⫺q 兲 2 关 ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫺ប⍀ f兲兴 兺 ជq ,
⫹ 共 N qជ ⫹1 兲共 k ⫺ k⫹q 兲 2 关 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫺ប⍀ f兲兴 其 ⫺
2ប
2 兩Cq,兩2M qជ 共 ប q ⫺បqជ •uជ 0 兲 2 关 N qជ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲 ⫹ 共 N qជ ⫹1 兲 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲兴 . 兺 ជq,
共25兲 The electron-phonon scattering depends not only on the electron distribution but also on the phonon distribution. The nonequilibrium phonon distribution can be computed from another scattering equation due to their coupling to electrons:17 d N qជ ⫺N 0 共 q 兲 em abs , N ជ ⫽⌰ qជ 共 N qជ ⫹1 兲 ⫺⌰ qជ N qជ ⫺ dt q q
共26兲
where qជ is the relaxation time of phonons from the boundary scattering. In Eq. 共26兲, the phonon-emission rate due to electron-phonon scattering to leading order is found to be10 em
⌰ qជ ⫽
4 兩 C q, 兩 2 ប
f 共ជ兲
兺kជ 共 kk 兲
再
冉
冊
1 2 1 2 关 ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲兴 1⫺ M qជ ⫹ M qជ 关 ␦ 共 k 2 4
⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫺ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫺ប⍀ f兲
冎
⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲兴 ⫹
4 兩 C q, 兩 2 ប
兺kជ
冋
共 k⫹q ⫺ k 兲 d f 共 kជ 兲 共 k兲 d kជ
册再
1 2 M 关 ␦ 共 k ⫺ k⫹q ⫹ប q 4 qជ
冉
1 2 ⫺បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫺ប⍀ f兲兴 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲 1⫺ M qជ 2
冊冎
共27兲
,
and the phonon-absorption rate is abs
⌰ qជ ⫽
4 兩 C q, 兩 2 ប
f 共ជ兲
兺kជ 共 kk 兲
再
冉
1 2 关 ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 兲兴 1⫺ M qជ 2
冊
1 2 ⫹ M qជ 关 ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫺ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q 4
冎
⫹បqជ •uជ 0⫺ប⍀f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប q ⫺បqជ •uជ 0 ⫹ប⍀ f兲兴 ⫹
4 兩 C q, 兩 2 ប
兺kជ
冋
共 k⫺q ⫺ k 兲 d f 共 kជ 兲 共 k兲 d kជ
册再
冉
1 2 M 关 ␦ 共 k ⫺ k⫺q 4 qជ
1 2 ⫺ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 ⫺ប⍀ f兲兴 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប q ⫹បqជ •uជ 0 兲 1⫺ M qជ 2 075214-6
冊冎
.
共28兲
PHYSICAL REVIEW B 69, 075214 共2004兲
HIGH-FIELD TRANSPORT OF ELECTRONS AND . . .
Equations 共11兲, 共12兲, 共16兲, and 共26兲 can be combined together to study the electron-phonon drag effect.7 The scattering interaction between drifting electrons and phonons constitutes a frictional force on the center-of-mass motion of electrons. The force-balance equation determines the center-of-mass drift velocity from Eq. 共11兲, where upon ignoring the frictional force due to impurity scattering in Eq. 共12兲 the remainder of the frictional force from phonons in Eq. 共14兲 can now be written as Fជ p [uជ d]⫽⫺4
1
兺kជ 共 k 兲
冋
d d kជ
f 共 kជ 兲
册兺
qជ 兩 C q, 兩 2 ( k⫹q ⫺ k )
qជ ,
再冉
冊
1 2 1⫺ M qជ [N 0 ( q )⫺N ⬘0 ( q ⫹qជ •uជ 0 )] 2
1 2 ⫻ ␦ ( k⫹q ⫺ k ⫹ប q ⫹បqជ •uជ 0 )⫹ M qជ 兵 [N 0 共 q 兲 ⫺N 0⬘ 共 q ⫹qជ •uជ 0 ⫹⍀ f兲 ] ␦ 共 k⫹q ⫺ k ⫹ប q ⫹បqជ •uជ 0 ⫹ប⍀ f兲 4
冎
⫹ 关 N 0 共 q 兲 ⫺N 0⬘ 共 q ⫹qជ •uជ 0 ⫺⍀ f兲兴 ␦ 共 k⫹q ⫺ k ⫹ប q ⫹បqជ •uជ 0 ⫺ប⍀ f兲 其 .
From Fro¨hlich electron-phonon coupling, we find the coupling matrix to be12 兩 C q兩 2⫽
冉 冊冉 ប LO 2V
冊
1 e2 1 ⫺ , ⑀ ⬁ ⑀ s ⑀ 0 共 q 2 ⫹Q s2 兲
共30兲
where LO is the frequency of dominant longitudinal-optical 共LO兲 phonon modes, ⑀ ⬁ and ⑀ s are the high-frequency and static dielectric constants of the host material, and Q s2 ⫽(e 2 / ⑀ 0 ⑀ r)(m * / 2 ប 2 )(3 2 3D) 1/3 represents the ThomasFermi screening effect. For acoustic phonon scattering, on the other hand, we use the deformation-potential approximation.17 This leads to 兩 C qᐉ 兩 2 ⫽
ប qᐉ 2 c 2ᐉ V
兩 C qt 兩 ⫽ 2
冋
D 2⫹
ប qt
9 32q
13
2 c 2t V 64q
共 eh 14兲 2 2
共 eh 14兲 2
2
册冉
q2 q 2 ⫹Q s2
冉 冊 q2
q 2 ⫹Q s2
冊
phonons under a weak electron-phonon coupling in the relative scattering motion of electrons 共neglecting the electronphonon drag effect兲. By using these approximations, analytical expressions for the expansion coefficients in the FokkerPlanck equation can be obtained. For the case with both dc and infrared fields, the thermal quantities A T( kជ ), V T( kជ ), and D T( kជ ) can be obtained from the functions defined in Eqs. 共A1兲–共A3兲 共in the appendix兲
A T共 k⬜ , k 储 兲 ⫽
4␣
共 ប LO兲 3/2
⫺ 关G⫹ 0 共 k⬜ ,k 储 兲 ⫺G 0 共 k⬜ ,k 储 兲兴 ,
共33兲
2
, 共31兲
V T共 k⬜ , k 储 兲 ⫽⫺
2
,
共29兲
共32兲
4␣
冑ប LO
⫺ 关G⫹ 1 共 k⬜ ,k 储 兲 ⫺G 1 共 k⬜ ,k 储 兲兴 ,
共34兲
where ⫽ᐉ,t corresponds to one longitudinal and two transverse acoustic-phonon modes, c ᐉ and c t are the sound velocities for these modes, is the ion mass density, D is the deformation-potential coefficient, and h 14 is the piezoelectric constant. Applying the Debye model to low-energy acoustic phonons, we get q ⫽c q with ⫽ᐉ,t. III. RADIATIVE EFFECT AND ELECTRON TRANSPORT UNDER DC AND INFRARED FIELDS
It is known that the scattering between LO phonons and electrons becomes dominant at high temperatures in polar semiconductors, such as GaAs. In this section, we discuss the radiative effect and electron transport under both an infrared field and a dc electric field by including the opticalphonon scattering at high temperatures in GaAs. We consider the indirect absorption of the infrared field by electrons in the relative motion and the resulting change in the nonequilibrium electron temperature for an anisotropic distribution of electrons. We assume an outside heat bath connected to the GaAs system to keep a constant lattice temperature. Under this assumption we further assume N qជ ⬇N 0 ( q ) for
FIG. 1. Illustration for orientations of electron (kជ ) and phonon ជ (q ) wave vectors. Eជ dc and uជ 0 lie in the 共z兲 direction. The electron wave vector kជ can be decomposed as a parallel component (kជ 储 ) 共in the z direction兲 and a perpendicular component kជ⬜ 共in the x direction兲. The phonon wave vector qជ can also be decomposed in the same way. The angle between qជ ⬜ and the x direction within the x-y plane is denoted as , and the angle between kជ and the z direction within the x-z plane is represented by .
075214-7
PHYSICAL REVIEW B 69, 075214 共2004兲
HUANG et al.
D T共 k⬜ , k 储 兲 ⫽
2 ␣ 冑ប LO ⫹ 关 G 2 共 k⬜ ,k 储 兲 ⫹G ⫺ 2 共 k⬜ ,k 储 兲兴 , 共35兲
2 where ␣ ⫽(e 2 /16 ⑀ 0 ) 冑2m * LO (1/⑀ ⬁ ⫺1/⑀ s), k⬜ ⫽ 冑2m * k⬜ /ប, and k 储 ⫽ 冑2m * k 储 /ប. The anisotropic dependence on k⬜ and k 储 in Eqs. 共33兲–共35兲 is a consequence of the applied dc electric field, as shown in Fig. 1. The field-induced quantities A F( kជ ), V F( kជ ), and D F( kជ ), 2 , can be obtained from the functions deproportional to E op fined in Eqs. 共A4兲–共A6兲:
A F共 k⬜ , k 储 兲 ⫽
␣ 共 ប LO兲 3/2
(⫹1) (⫺1) 关R⫹ 共 k⬜ ,k 储 兲 ⫹R ⫹ 共 k⬜ ,k 储 兲
(⫹1) (⫺1) ⫺R ⫺ 共 k⬜ ,k 储 兲 ⫺R ⫺ 共 k⬜ ,k 储 兲兴
⫺
2␣
共 ប LO兲
3/2
(0) (0) 关R⫹ 共 k⬜ ,k 储 兲 ⫺R ⫺ 共 k⬜ ,k 储 兲兴 ,
共36兲 V FT共 k⬜ , k 储 兲 ⫽⫺
␣ 冑ប LO
(⫹1) (⫺1) 关S⫹ 共 k⬜ ,k 储 兲 ⫹S ⫹ 共 k⬜ ,k 储 兲
(⫹1) (⫺1) ⫺S ⫺ 共 k⬜ ,k 储 兲 ⫺S ⫺ 共 k⬜ ,k 储 兲兴 ,
V J共 k⬜ , k 储 兲 ⫽
2␣
冑ប LO
共37兲
Fp N em * ⫻
⫽
冋
3D
f 共 k⬜ , k 储 兲
⫹
k⬜
⫺⬁
k⬜ dk⬜
0
f 共 k⬜ , k 储 兲 k储
⫹⬁
册
1 共 k⬜ ⫹ k 储 兲
关 H 共 k⬜ ,k 储 兲 ⫹W 1 共 k⬜ ,k 储 兲
共41兲
The anisotropic scattering contributes to the frictional force in Eq. 共41兲 in two respects: 共i兲 through the functions H, W 1 , defined in Eqs. 共A7兲–共A9兲 and 共ii兲 through the and W ⫾1 2 distribution function f ( k⬜ , k 储 ) obtained from the FokkerPlanck equation for degenerate electrons. Having calculated these expansion coefficients in Eqs. 共33兲–共40兲 and the frictional force in Eq. 共41兲, the Fokker-Planck equation in Eq. 共19兲 can be solved simultaneously with the force-balance equation in Eq. 共11兲 along with the boundary conditions f (0, 0)⫽ f (0, k 储 )/ k⬜ ⫽ f ( k⬜ ,0)/ k 储 ⫽0. When electrons stay in a nonequilibrium state, we cannot define an electron temperature through the Fermi-Dirac distribution.3 However, we can still define a general nonequilibrium electron temperature through the average energy of each electron in the relative motion. By calculating (T e) for fixed T e and 3D, T e can be found from the following equation:2
m * 3D
冕
⫹⬁
0
再冕
共38兲
⫽ e
␣ 冑ប LO (⫹1) (⫺1) D FT共 k⬜ , k 储 兲 ⫽ 关 T ⫹ 共 k⬜ ,k 储 兲 ⫹T ⫹ 共 k⬜ ,k 储 兲 2
⫹
共39兲
and D F( kជ )⫽D FT( kជ ) V F( kជ )⫽V FT( kជ )⫹V J( kជ ), ⫹D A( kជ ). V FT( kជ ) and D FT( kជ ) with M ⫽⫾1 represent the infrared-field-induced corrections to relevant thermal quantities V T( kជ ) and D T( kជ ), respectively. V J( kជ )⬎0 for M ⫽0 comes from the effect of Joule heating due to power absorption from the infrared field by electrons, and D A( kជ )⬍0 for M ⫽0 is a result of the antidiffusion of electrons in energy space due to a correction to spontaneous phonon emission by Joule heating.12 The anisotropic dependence on k⬜ and k 储 in Eqs. 共36兲–共40兲 reflects different radiative effects for two polarizations 共parallel and perpendicular to uជ 0 ) of the incident infrared field in addition to the existence of a dc field 共parallel to uជ 0 ), as shown in Fig. 1. With the help of functions defined in Eqs. 共A7兲–共A9兲, the frictional force from phonon scattering can be expressed as
2m *
dk 储
⫹W (⫹1) 共 k⬜ ,k 储 兲 ⫹W (⫺1) 共 k⬜ ,k 储 兲兴 . 2 2
(0) (0) 关S⫹ 共 k⬜ ,k 储 兲 ⫺S ⫺ 共 k⬜ ,k 储 兲兴 ,
␣ 冑ប LO (0) (0) D A共 k⬜ , k 储 兲 ⫽⫺ 关 T ⫹ 共 k⬜ ,k 储 兲 ⫹T ⫺ 共 k⬜ ,k 储 兲兴 , 共40兲
⫹⬁
1
3
2ប2
(⫹1) (⫺1) ⫹T ⫺ 共 k⬜ ,k 储 兲 ⫹T ⫺ 共 k⬜ ,k 储 兲兴 ,
冑 冕 冕
4␣
k 4 dk 1⫹exp兵 关 k ⫺ 共 T e兲兴 /k BT e其
p
0
1 3D
1 ជ dc⫹Eជ op共 t 兲兴 •uជ d共 t 兲 dt⫺ m * 兩 uជ d共 t 兲 兩 2 关E 2
冕
⫹⬁
0
f 共 kជ 兲 k d kជ ,
冎 共42兲
where (T e) is the chemical potential for free electrons at the temperature T e . The combination of the first and second terms on the right-hand side of Eq. 共42兲 represents the net increase of internal energy of electrons from frictional forces, and the last term comes from the change in the average kinetic energy of electrons in the relative motion mainly due to indirect photon absorption. This T e depends on both external and system parameters, such as 3D , ប LO , m * , E dc , E op , and ⍀ f , but not on the electron kinetic energy, and it is different from the lattice temperature T. For very high lattice temperatures at which the quantum degeneracy can be neglected, the left side of Eq. 共42兲 can simply be replaced by 23 k BT e In Eq. 共42兲, the momentumrelaxation time p is given by 1 1 ⫽ p 3D
冕
⫹⬁
0
冋
1
T共 kជ 兲
⫹
1
F共 kជ 兲
册
f 共 kជ 兲 d kជ ,
共43兲
where the state-dependent nonlocal quantities can be expressed through
075214-8
PHYSICAL REVIEW B 69, 075214 共2004兲
HIGH-FIELD TRANSPORT OF ELECTRONS AND . . .
1
T共 kជ 兲 1
R共 kជ 兲
⫽
2ប
⫽
2 ប
兺qជ 兩 C q兩 2 关 Rk,k⫹q ␦ 共 k ⫺ k⫹q ⫹ប LO⫺បq 储 u 0 兲 ⫹Rk,k⫺q ␦ 共 k ⫺ k⫺q ⫺ប LO⫹បq 储 u 0 兲兴 ,
共44兲
兺qជ 兩 C q兩 2 M q2ជ 兵 Rk,k⫹q 关 ␦ 共 k ⫺ k⫹q ⫹ប LO⫺បq 储 u 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫹q ⫹ប LO⫺បq 储 u 0 ⫺ប⍀ f兲兴
⫹Rk,k⫺q 关 ␦ 共 k ⫺ k⫺q ⫺ប LO⫹បq 储 u 0 ⫹ប⍀ f兲 ⫹ ␦ 共 k ⫺ k⫺q ⫺ប LO⫹បq 储 u 0 ⫺ប⍀ f兲兴 其 ⫺
ប
兺qជ 兩 C q兩 2 M q2ជ 关 Rk,k⫹q ␦ 共 k ⫺ k⫹q ⫹ប LO⫺បq 储 u 0 兲 ⫹Rk,k⫺q ␦ 共 k ⫺ k⫺q ⫺ប LO⫹បq 储 u 0 兲兴 .
The magnitude of 1/ p reflects the strength of the frictional force. In the above two equations, we have defined the angular-distribution factor for electron momentum dissipation18 Rk,k⫾q ⫽1⫺
kជ • 共 kជ ⫾qជ 兲 . 兩 kជ • 共 kជ ⫾qជ 兲 兩
共46兲
IV. NUMERICAL RESULTS AND DISCUSSION
In our numerical calculations, we have chosen GaAs as ជ dc is assumed in the z direction, as shown the host material. E ជ in Fig. 1, and E op is chosen to be either parallel or perpenជ dc . For GaAs we have taken parameters as foldicular to E lows: m * ⫽0.067m 0 with free-electron mass m 0 , 3D⫽1 ⫻1018 cm⫺3 , ប LO⫽36 meV, ⑀ r⫽12, ⑀ s⫽11, ⑀ ⬁ ⫽13, ប⍀ f⫽25 meV, and T⫽300 K. Other parameters, such as E dc and E op , will be given directly in the figure captions. In the numerical calculations below, we will show results only for steady-state cases. In Fig. 2 the scaled distribution function ¯f ( k⬜ , k 储 ) 共in units of 3D / F) is exhibited as a function of k⬜ 共dashed curve with k 储 ⫽0) and of k 储 共solid curve with k⬜ ⫽0) for ជ dc with E dc⫽0.58 kV/cm and E op⫽20 kV/cm. A subEជ op储 E stantial difference between our model and the energy-balance approach11 can be seen clearly here; energy-balance equation give the dashed curve. It is evident from this figure that the distribution function of electrons exhibits an anisotropic dependence on k 储 and k⬜ , which is absent in the energybalance approach. This anisotropy is attributed to the exisជ dc and Eជ op applied to the GaAs system. tence of E We display in Fig. 3 the scaled distribution function ¯f ( k , k ) as a function of parallel kinetic energy k of 储 储 ⬜ ជ op储 Eជ dc 共solid curve兲 and Eជ op⬜Eជ dc electrons with k ⫽0 for E ⬜
共dashed curve兲 with E dc⫽0.58 kV/cm and E op⫽20 kV/cm. From this figure we find that the distribution function of electrons evidently depends on the polarization of the incident infrared field. This is a consequence of the different modifications to electron-phonon scattering in the presence of an incident infrared field with parallel and perpendicular
共45兲
ជ dc . This anisotropic coupling to the polarization relative to E infrared field will contribute differently to u 0 共compare the solid curves in Figs. 4 and 5兲. In Fig. 4 we show the calculated drift velocity u 0 共solid curve and left axis兲 and the nonequilibrium electron temperature T e 共dashed curve and right axis兲 as functions of the ជ op储 Eជ dc with amplitude of the incident infrared field E op for E E dc⫽0.58 kV/cm. From the figure we find that u 0 increases sublinearly with E op in the range of E op shown, similar to results reported by Lei at high E op , T⫽10 K, and terahertz frequencies.11 However, the higher temperature T⫽300 K used in our calculation brings u 0 down below 107 cm/s. In addition, the increase of T e is greatly limited at T⫽300 K and the range for the ratio T e /T is significantly reduced. In addition, we present the calculated drift velocity u 0 共solid curve and left axis兲 and nonequilibrium electron temperature T e 共dashed curve and right axis兲 as functions of E op
FIG. 2. Scaled distribution function ¯f ( k⬜ , k 储 ) in units of 3D / F as a function of parallel kinetic energy k 储 / F 共solid curve with k⬜ ⫽0) and of perpendicular kinetic energy k⬜ / F 共dashed ជ dc . Here, k ⫽ប 2 k 2储 /2m * , k curve with k ⫽0) for Eជ op储 E 储
⫽ប 2 k⬜2 /2m * , F⫽ប 2 (3 2 3D) 2/3/2m * ⫽0.58 kV/cm, and Eជ op⫽20 kV/cm.
075214-9
⬜
储
at
T⫽0 K,
Eជ dc
PHYSICAL REVIEW B 69, 075214 共2004兲
HUANG et al.
FIG. 3. Scaled distribution function ¯f ( k⬜ , k 储 ) in units of 3D / F as a function of parallel kinetic energy k 储 / F of electrons with k ⫽0 for parallel polarization Eជ op储 Eជ dc 共solid curve兲 and per⬜
ជ op⬜Eជ dc 共dashed curve兲. Here, Eជ dc pendicular polarization E ⫽0.58 kV/cm and Eជ op⫽20 kV/cm.
ជ op⬜Eជ dc with E dc⫽0.58 kV/cm. By comparing in Fig. 5 for E this figure with Fig. 4 we find that u 0 decreases with E op instead of increasing with E op in this case. Although u 0 decreases with E op , T e still increases with E op due to the increased average kinetic energy of electrons by power absorption from the applied infrared field. The effect of the polarization of the incident infrared field can be further visualized by looking at the momentumជ op⬜Eជ dc relaxation rate 1/ p in Fig. 6 as a function of E op for E ជ ជ 共solid curve兲 and E op储 E dc 共dashed curve兲 with E dc ⫽0.58 kV/cm. From this figure we find that the momentum-
FIG. 4. Calculated drift velocity u 0 in units of 106 cm/s 共solid curve and left axis兲 and nonequilibrium electron temperature T e in units of K 共dashed curve and right axis兲 as functions of the amplitude of the incident infrared field E op . Here, E dc is set to be 0.58 ជ op is assumed to be parallel to Eជ dc . kV/cm and E
FIG. 5. Calculated drift velocity u 0 in units of 106 cm/s 共solid curve and left axis兲 and nonequilibrium electron temperature T e in units of K 共dashed curve and right axis兲 as functions of the amplitude of the incident infrared field E op . Here, E dc is set to be 0.58 ជ op is assumed to be perpendicular to Eជ dc . kV/cm and E
relaxation process is suppressed by the incident infrared field for parallel polarization but enhanced for perpendicular polarization. This result provides an explanation for Figs. 4 and 5, where u 0 for parallel polarization in Fig. 4 increases with E op due to the suppression of momentum relaxation or frictional force, while u 0 for perpendicular polarization in Fig. 5 decreases with E op due to the enhancement of momentum relaxation.
V. CONCLUSIONS
In conclusion, we have developed a model, in which electrons are interacting with an infrared field in a bulk semicon-
FIG. 6. Calculated momentum-relaxation rate 1/ p in units of THz as a function of the amplitude of the incident infrared field E op ជ dc 共solid curve兲 and parallel for perpendicular polarization Eជ op⬜E ជ ជ polarization E op储 E dc 共dashed curve兲. Here, E dc is set to be 0.58 kV/cm.
075214-10
PHYSICAL REVIEW B 69, 075214 共2004兲
HIGH-FIELD TRANSPORT OF ELECTRONS AND . . .
ductor while they are undergoing transport due to an applied strong dc electric field, by coupling the force-balance equation for the slow center-of-mass motion of electrons and the Fokker-Planck equation for the ultrafast relative scattering motion of degenerate electrons. This allows us to fully include in our model the anisotropic energy-relaxation process, which was neglected in the past when the first-order energybalance equation was used. At the same time, it also allows us to include the anisotropic coupling of electrons to the incident infrared field which gives rise to different frictional forces for polarizations of the infrared field perpendicular or parallel to the applied dc field. Based on this model, we can study high-field transport of electrons under strong dc and infrared fields beyond the relaxation-time approximation. From our numerical study, we have discovered the anisotropic dependence of the electron distribution function on parallel and perpendicular kinetic energies of electrons. We have also quantified the effect of anisotropic coupling to the incident infrared field with polarizations parallel and perpendicular to the applied dc electric field. We have found that the drift velocity of electrons increases with the amplitude of the infrared field for the parallel polarization due to a suppression of momentum relaxation but decreases with the infrared-field amplitude for the perpendicular polarization due to an enhancement of momentum relaxation. The heating of electrons has been precisely described beyond the firstorder energy-balance equation in addition to the inclusion of anisotropic coupling to an incident infrared field having different polarizations. The first requirement for the Fokker-Planck equation to be a good approximation to the Boltzmann scattering equation in the relative scattering motion of electrons is that the Fermi energy F must be greater than k BT for equilibrium electrons. The other requirements include FⰇប q , ប⍀ f , and បq 储 u 0 . When a dipole coupling of electrons to an external optical field exists, i.e., using a grating to produce a nonuniform field or with photon energy greater than the band gap of the host materials, the incoherent Boltzmann scattering equation should be replaced by full coherent density-matrix equations with off-diagonal matrix elements.
⫾ 共 kជ ,qជ ,M 兲 ⫽
⫾共kជ,qជ ,M 兲⫽␦
冋
The authors are grateful to Professor S. W. Koch and Professor J. V. Moloney for their helpful discussions and comments on the relative scattering motion of electrons. APPENDIX
We can decompose the electron and phonon wave vectors ជ dc and perpendicular to it, as shown along the direction of E in Fig. 1. Here, for optical-phonon scattering of electrons at high temperatures in polar semiconductors we first define three functions that determine A T( kជ ), V T( kជ ), and D T( kជ ) ជ dc⫽0 but Eជ op⫽0 关see Eqs. 共33兲–共35兲兴: in the case of E G⫾ 0 共 k⬜ ,k 储 兲 ⫽
冑2m * ប LO ប ⫻
G⫾ 1 共 k⬜ ,k 储 兲 ⫽
⫹⬁
⫺⬁
dq 储
册
⫹⬁
0
dq⬜ ⫾ 共 kជ ,qជ ,0 兲
ប
冕
⫹⬁
⫺⬁
q⬜ 2 q储 ⫹q⬜2 ⫹Qs2
共A1兲
冉
dq 储
1⫺
冕
⫹⬁
0
q 储u 0 LO
dq⬜ ⫾ 共 kជ ,qជ ,0 兲
冊冋
册
1 1 N 0 共 LO兲 ⫹ ⫾ , 2 2 共A2兲
G⫾ 2 共 k⬜ ,k 储 兲 ⫽
冑2m * ប LO ប ⫻
冕
⫺⬁
q⬜ q2储 ⫹q⬜2 ⫹Qs2
冕
⫹⬁
冉
dq 储
1⫺
⫹⬁
0
q 储u 0 LO
dq⬜ ⫾ 共 kជ ,qជ ,0 兲
冊冋 2
册
1 1 N 0 共 LO兲 ⫹ ⫾ . 2 2 共A3兲
In the above three equations, we have introduced the following function for k⬜ q⬜ ⫽0:
冑4k⬜2 q⬜2 ⫺ 关 q⬜2 ⫹q 2储 ⫾2k 储 q 储 ⫿2m * 共 LO⫺q 储 u 0 ⫹M ⍀ f兲 /ប 兴 2
2m* 共LO⫺q 储 u 0 ⫹M ⍀ f兲 . ប
冕
q⬜ , 2 q 储 ⫹q⬜2 ⫹Q s2
冑2m * ប LO
⫻
冕
共 2k⬜ q⬜ ⫺ 兩 q⬜2 ⫹q 2储 ⫾2k 储 q 储 ⫿2m * 共 LO⫺q 储 u 0 ⫹M ⍀ f兲 /ប 兩 兲
where (x) is a unit step function. On the other hand, for k⬜ q⬜ ⫽0 we define
q⬜2 ⫹q2储 ⫾2k储q储⫿
ACKNOWLEDGMENTS
,
M ⫽0,⫾1 关see Eqs. 共36兲–共40兲兴:
(M ) R⫾ 共 k⬜ ,k 储 兲 ⫽
Moreover, we define another three functions that will determine A F( kជ ), V F( kជ ), and D F( kជ ) in the case of Eជ op⫽0 with 075214-11
␥ 冑2m * ប LO ប ⫻
冕
⫹⬁
⫺⬁
dq 储
q⬜ 共 ␣⬜ q⬜2 ⫹ ␣ 储 q 2储 兲 q 2储 ⫹q⬜2 ⫹Q s2
,
冕
⫹⬁
0
dq⬜ ⫾ 共 kជ ,qជ ,M 兲 共A4兲
PHYSICAL REVIEW B 69, 075214 共2004兲
HUANG et al.
(M ) S⫾ 共 k⬜ ,k 储 兲 ⫽
␥ 冑2m * ប LO ប ⫻
冕
⫹⬁
⫺⬁
dq 储
q⬜ 共 ␣⬜ q⬜2 ⫹ ␣ 储 q 2储 兲
冋
q 2储 ⫹q⬜2 ⫹Q s2
冉 册
冕
⫹⬁
0
q 储u 0 M ⍀ f 1⫺ ⫹ LO LO
冕
⫹⬁
⫺⬁
dq 储
q⬜ 共 ␣⬜ q⬜2 ⫹ ␣ 储 q 2储 兲 ⫻ q 2储 ⫹q⬜2 ⫹Q s2
冋
冉 册
⫻ 关 N 0 共 LO兲 ⫺N 0⬘ 共 LO⫹q 储 u 0 兲兴 1⫹
冕
⫻
0
1⫺
q 储u 0 M ⍀ f ⫹ LO LO
冊
⫹⬁
⫺⬁
dq 储
冕
⫹⬁
0
共 kជ ,qជ ,M 兲 ⫽
共A6兲
⫺⬁
冕
⫹⬁
0
dq⬜ 共 kជ ,qជ ,0 兲
q⬜ q 储 共 ␣⬜ q⬜2 ⫹ ␣ 储 q 2储 兲 q 2储 ⫹q⬜2 ⫹Q s2
q⬜ q 储 2 q 储 ⫹q⬜2 ⫹Q s2
␥ 4
冕
⫻
⫹⬁
⫺⬁
dq 储
冕
⫹⬁
0
冉
关 N 0 共 LO兲
冊
q 储u 0 , LO
共A8兲
dq⬜ 共 kជ ,qជ ,M 兲
q⬜ q 储 共 ␣⬜ q⬜2 ⫹ ␣ 储 q 2储 兲 q 2储 ⫹q⬜2 ⫹Q s2
冉
⫹q 储 u 0 ⫹M ⍀ f兲兴 1⫹
关 N 0 共 LO兲 ⫺N 0⬘ 共 LO
冊
q 储 u 0 ⫹M ⍀ f . LO
共A9兲
In the above three equations, we have introduced the following function for k⬜ q⬜ ⫽0:
共2k⬜q⬜⫺兩q⬜2 ⫹q2储 ⫹2k储q储⫹2m*共q储u0⫹LO⫹M ⍀ f兲 /ប 兩 兲
冑4k⬜2 q⬜2 ⫺ 关 q⬜2 ⫹q 2储 ⫹2k 储 q 储 ⫹2m * 共 q 储 u 0 ⫹ LO⫹M ⍀ f兲 /ប 兴 2
On the other hand, for k⬜ q⬜ ⫽0 we define
冋
共kជ,qជ ,M 兲⫽␦ q⬜2 ⫹q2储 ⫹2k储q储⫹
1
dq 储
2
) W (M 2 共 k⬜ ,k 储 兲 ⫽
dq⬜ 共 kជ ,qជ ,0 兲
⫹⬁
⫺N 0⬘ 共 LO⫹q 储 u 0 兲兴 1⫹
where ␥ ⫽(eE op /m * ⍀ f2 ) 2 . We take ␣⬜ ⫽0 and ␣ 储 ⫽1 for ជ op储 Eជ dc , but we set ␣ 储 ⫽0 and ␣⬜ ⫽1/2 for Eជ op⬜Eជ dc 共unpoE larized infrared field兲. Finally, for the frictional force in Eq. 共41兲 we define the following functions with M ⫽0,⫾1:
冕
冕
dq⬜ ⫾ 共 kជ ,qជ ,M 兲
1 1 ⫻ N 0 共 LO兲 ⫹ ⫾ , 2 2
H 共 k⬜ ,k 储 兲 ⫽
␥ 2
W 1 共 k⬜ ,k 储 兲 ⫽⫺ 共A5兲
⫹⬁
冊
q 储u 0 , LO 共A7兲
冊
1 1 ⫻ N 0 共 LO兲 ⫹ ⫾ , 2 2
␥ 冑2m * ប LO (M ) T⫾ 共 k⬜ ,k 储 兲 ⫽ ប
冉
dq⬜ ⫾ 共 kជ ,qជ ,M 兲
H. Fro¨hlich and B.V. Paranjape, Proc. Phys. Soc. London, Sect. B 69, 21 共1956兲. 2 M.R. Arai, Appl. Phys. Lett. 42, 906 共1983兲. 3 X.L. Lei and C.S. Ting, Phys. Rev. B 32, 1112 共1985兲. 4 J. M. Ziman, Principles of the Theory of Solids, 1st ed. 共Cambridge University Press, Cambridge, England, 1964兲, pp. 179– 186. 5 X.-G. Zhao, G.A. Georgakis, and Q. Niu, Phys. Rev. B 56, 3976 共1997兲. 6 R. Kubo, J. Phys. Soc. Jpn. 12, 570 共1957兲. 7 F.G. Bass and Yu.G. Gurevich, Sov. Phys. Usp. 14, 113 共1971兲; Yu.G. Gurevich and O.L. Mashkevich, Phys. Rep. 181, 327 共1989兲. 8 C.S. Ting, S.C. Ying, and J.J. Quinn, Phys. Rev. B 14, 4439 共1976兲. 9 T. Apostolova, D.H. Huang, and D.A. Cardimona, Phys. Rev. B 67, 205323 共2003兲.
.
册
2m* 共q储u0⫹LO⫹M ⍀ f兲 . ប
10
A. Kaiser, B. Rethfeld, H. Vicanek, and G. Simon, Phys. Rev. B 61, 11 437 共2000兲. 11 X.L. Lei, J. Appl. Phys. 84, 1396 共1998兲. 12 T. Apostolova, D.H. Huang, P.M. Alsing, J. McIver, and D.A. Cardimona, Phys. Rev. B 66, 075208 共2002兲. 13 M. Lindberg and S.W. Koch, Phys. Rev. B 38, 3342 共1988兲; J.V. Moloney, R.A. Indik, J. Hader, and S.W. Koch, J. Opt. Soc. Am. B 11, 2023 共1999兲. 14 L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 3rd ed. 共Pergamon Press, Oxford, 1976兲, pp. 1– 45. 15 X.L. Lei and C.S. Ting, J. Phys. C 18, 77 共1985兲. 16 L.H. Holway and D.W. Fradin, J. Appl. Phys. 46, 279 共1975兲. 17 S.K. Lyo and D.H. Huang, Phys. Rev. B 66, 155307 共2002兲. 18 M. Sparks, D.L. Mills, R. Warren, T. Holstein, A.A. Maradudin, L.J. Sham, E. Loh, Jr., and D.F. King, Phys. Rev. B 24, 3519 共1981兲.
075214-12