Electron-acoustic phonon field induced tunnel scattering

Theory of electron-acoustic single phonon scattering has been reconsidered. It is assumed that the non-degenerate semiconductor has a spherical parabolic band structure. In the basis of the reconsideration there is a phenomenon of the tilting of semiconductor bands by the perturbing potential of an electric field. In this case, electron eigenfunctions are not plane waves or Bloch functions. In low-field regime, the expressions for electron intraband transition probability and scattering time are obtained under elastic collision approximation. Dependencies of scattering time on electron energy and uniform electric field are analyzed. The results of corresponding numerical computations for n-Si at 300 K are presented. It is established that there is no fracture on the curve of electron scattering time dependence on the electron energy.


Introduction
Current carrier (electron) mobility µ is an important parameter characterizing many transport phenomena in semiconductors under electric field F . Electron mobility is determined as µ = 〈eτ k /m〉 [1,2], where τ k is the electron quasi-momentum relaxation time, e is the electron charge magnitude, m is the electron effective mass, 〈· · · 〉 is the symbol of averaging over conduction zone quantum states. Relaxation time τ k is determined by the electron scattering by various dynamic and static imperfections of a crystall lattice such as lattice vibrations (optic and acoustic phonons), ionized and neutral impurity atoms, vacancies, etc. For theoretical consideration of scattering probability and relaxation time τ k , a flat-band semiconductor model is used, as a rule [1][2][3][4]. In this case, in low electric field region (F < F c , where F c is a characteristic field), the electron relaxation time τ k and, therefore, the mobility µ are field-independent quantities [1,4,5]. Particularly, the relaxation time related to electron-acoustic single phonon elastic scattering in non-degenerate n-type semiconductor with a spherical parabolic conduction band is given by [1][2][3][4] 1/τ k,ac = D 2 ac (2m) 3/2 k B T 2πħ 4 ρ r ν 2 0 ε k . F c ≈ 2.8 · 10 3 V/cm [5]. The dependence µ(F ) is explained by the phenomenon of electron gas heating-up under the effect of a high electric filed [1,4]. However, in recent work [6], a new mechanism of electron lattice scattering, referred to as electron-phonon FIT (field induced tunnel) scattering, is observed. In the basis of the FIT scattering there lies a phenomenon of tilting of semiconductor bands by the perturbing potential of an electric field. The effect of the electron-phonon FIT scattering is explained in terms of penetration of an electron wave function into a semiconductor band gap in the presence of an electric field. Contrary to a flat-band semiconductor, in a tilted-band semiconductor, a conduction electron transition in the band gap region is allowed. In [6], reconsidering the electron-phonon interaction theory, the case of electron intraband FIT scattering by non-polar optical phonon is analyzed. In the present work, electron-acoustic phonon FIT intraband scattering is considered. It is assumed that the non-degenerate n-type semiconductor has a parabolic conduction band.

Electron-acoustic phonon FIT transition probability
To theoretically characterize the carrier scattering, it is necessary to consider the scattering probability and evaluate the relaxation time. The task of the transition probability calculation is solved based on the perturbation theory (see, e.g. [1][2][3][4]). According to this theory, the probability per unit time of quantum system transition from λ state to λ state, to the first order in the perturbation, is determined as [1,2] Here, a λ ,λ (t ) is the perturbation matrix element: ψ λ (R, t ) = ψ λ (R) exp(−iE λ t /ħ), λ is the set of quantum numbers characterizing different states of a nonperturbed system, ψ λ (R) and E λ are the wave function and energy eigenvalues of stationary state of a non-perturbed quantum system, respectively,Ŵ (R, t ) is the perturbation operator, R is the set of the quantum system coordinates, V is the volume, ' * ' is the complex conjugate symbol.
Consideration of the electron scattering by phonons is based on (2) and (3) as well. In brief, the description of the probability calculation is as follows. In the present case, the quantum system consists of a conduction electron in the crystall periodic field and lattice normal vibrations. Then, λ should be replaced by an electron quasi-wave vector k and by the phonon occupation numbers of all possible states. The electron-phonon interaction Hamiltonian H e-ph is taken as perturbationŴ (R, t ); the r radius vector and the normal coordinates of lattice vibration are taken as a quantum system coordinate R. The wave function of non-perturbed state of an electron-phonon system is expressed as a product of one-electron wave function and harmonic oscillator wave functions.
In flat-band semiconductors, an electron state is described by Bloch functions. Therefore, for an electron scattered from an initial state |k〉 to a final state |k 〉, the transition probability per unit time W (k, k ) is evaluated based on the Bloch functions [1][2][3][4]. To simplify the calculations, a plane wave ψ k (r) = e ikr / L x L y L z (where L x , L y and L z are the sizes of a semiconductor, V = L x L y L z ) is used sometimes as an electron wave function [1][2][3][4] (nearly free electron approximation). Calculations of the probability W (k, k ) based on the Bloch function or plane wave are well known and reported in detail in numerous publications (see, for example, [1][2][3][4]). However, in the present work, a semiconductor, whose bands are tilted by the perturbing potential of a uniform electric field, is of interest. In this case, electron eigenfunctions are not plane waves or Bloch functions [1,3,7]. Here, based on such an assumption, the probability W (k, k ) is recalculated following the above-mentioned general approach. It is assumed that the perturbation Hamiltonian H e−ph , which is a harmonic function of time (harmonic perturbation), is determined using the deformation potential theory [1][2][3][4] of electron-phonon interaction. The results of our calculations show that: • the probability of electron transition with phonon absorption can be presented as follows: • the probability of electron transition with phonon emission can be presented as follows: Here, ψ k (r, t ) = ψ k (r) exp(−iε k t /ħ), ψ k (r) is the conduction electron wave function of stationary state k, ω q is the phonon angular frequency, n q is the occupation number of the equilibrium phonons, which is given by the Bose-Einstein distribution: The summations in (4) and (5) should be carried out in the range of the first Brillouin zone (BZ).
After integration over t , (4) and (5) are expressed as follows: Here and bellow, upper and lower symbols of double signs refer to electron transitions with the q phonon absorption and emission, respectively.
In the presence of a uniform electric field F (which is parallel to the z-axis), the electron wave-function is determined by the stationary Schrödinger equation which in the effective mass approximation is expressed as follows: Here and below, the index k of ε k is omitted to simplify the expressions, so we can substitute ε k → ε and ε k → ε in further expressions.

23702-3
where ε n = ħ 2 k 2 n /2m, the index n identifies the electron energy eigenvalues, k 2 ⊥ = k 2 x + k 2 y , k ⊥ is the electron wave-vector perpendicular to the electric field. Energy eigenvalues ε n (or k n ) are determined from the boundary conditions (see below) for the wave function χ n (z).
For a semiconductor of length L z in z-direction (i.e., −L z /2 z L z /2) from (9) and relation [1,3] 1 L z it follows Inserting (15) into (7) yields The Kronecker δ in this equation expresses the laws of conservation of perpendicular to the electric field quasi-momentum x, y components of the scattered particles. Here, only the normal N-processes of scattering are considered.
After summation with respect to q x and q y with the help of the Kronecker δ, (16) becomes as follows: Taking into account that ω q is an even function on q (and particularly on q z ), from (17) one obtains Using the following formal transformations (18) may be written as follows: For a sufficiently long time t , when the relation lim t →∞ sin(at )/a = πδ(a) can be used [1,2], (20) is expressed as follows: In this expression, the summation with respect to q 1,z is non-zero only when q z = q 1,z . Then, electron transition probability can be written as follows: Here, the Dirac δ-function indicates the energy conservation law. Note, if in (22) the plane wave (χ n (z) = e ik z z / L z ) or Bloch function is taken as an electron wave function, then the well-known classical expression of probability W a,e (k, k ) [1][2][3][4] (Fermi Golden rule) is derived. In the common case, the calculation of the integral and sum in (22) is complicated. On the other hand, the expression (1) for τ k,ac is derived (within the framework of the flat-band semiconductor model) by using the approximation of elasticity of electron-acoustic phonon scattering (elastic collision approximation, ħω q ε k ) with the assumption that Here, ω q = ν 0 q dispersion law of the long wavelength longitudinal acoustic phonon is used.
To reveal the difference between the results of flat and titled band approaches it is reasonable that these assumptions should be used here as well. Then, after simple summing with respect to q x and q y , (22) takes the following simpler form:
The result is as follows: Delta-integration over z 1 yields: Electron energy eigenvalues ε n are determined from boundary conditions to (10). Note, there is some difference between the peculiarities of the movement of an electron under electric field in vacuum and an electron in a semiconductor. Conduction band gap ∆E c of a semiconductor is a finite quantity. Contrary to the vacuum, the movement of a conduction electron in a semiconductor with perfect lattice has an oscillation character, as it is shown in figure 1. Electron 'el.1' oscillates between the bottom of the conduction band and semiconductor edge; electron 'el.2' periodically reflects from the conduction band top and bottom edges. The behavior of an 'el.2' is well-known as a Bloch oscillation. Boundary conditions for Schrödinger equation (10) and, therefore, electron energy eigenvalues in case of Bloch oscillations have been reported, for example, in [9]. Here, for definiteness, neglecting the Bloch oscillations, the case of 'el.1' is considered only, i.e., it is assumed that the electric field is low and the magnitude of ∆E c is very large. In other words, the model of a triangular quantum well with finite sizes is considered.
Here, we are interested in a large length semiconductor, particularly, in z-direction. Then, the allowed values of k n (or ε n ) are computed from the boundary condition χ n (z = −L z /2) = 0 as (see, [7] and figure 1) where a n are the zeros of the Airy function which are located in the negative part of the real axis [8].
They are well approximated as a n −(3π Solving (30) for k n and inserting the result into (13), one obtains energy eigenvalues [7]: The values ε in (31) are obtained as a function of the electron state quantum numbers. The quantities k x , k y (or k ⊥ ) and n (or k n ) are a set of quantum numbers which determine the conduction electron state in the presence of an electric field.
Normalization constant C n in (11) is determined as follows: It is known [7,8] that ds Ai 2 (s) = s Ai 2 (s) − Ai 2 (s), where Ai is the derivative of the Airy function. Consequently, normalization constant C n can be presented as follows: For a semiconductor with large length L z one has C −2 n = l Ai 2 (a n ).
Here, we used the fact that Ai(a n ) = 0; functions Ai(s) and Ai (s) exponentially vanish for positive large argument s [8]. On the other hand, the value Ai (a n ) is well approximated as [8] Ai (a n ) (−1) n−1 1 π (−a n ) 1/4 .

23702-7
Thus, using (11), (29) and taking into account that χ n (z) is the real function [see, (11), (12)] the expression (28) can be presented as follows: du Ai u − (k n l ) 2 2 Ai u − (k n l ) 2 2 where u = z/l is the dimensionless variable of integration. (38) describes electron transition probability in the presence of an electric field. Transition probability depends on the electric field and it is a symmetric function: W a,e (k, k ) = W a,e (k , k).

Electron-acoustic phonon FIT scattering time
The scattering time is defined by [3]: (38) shows that in the present case, the probabilities of phonon absorption and emission by electron are the same: W a (k , k) = W e (k, k ). Then, (39) can be written as follows: In the common case, the scattering time (the inverse of the scattering rate) differs from the relaxation time although sometimes both of them are equivalent (for example, the above-mentioned case of (1)) [3].
To calculate τ −1 k,sc in (40), we replace the sums over k x and k y by integrals over k x and k y , respectively. The transition from the sum to integral with the help of the relation dk x dk y = π dk 2 ⊥ can be presented as follows: Here, coefficient 2 in the numerator is related to the electron spin.
From (30) it follows that k 2 n+1 − k 2 n ∼ 1/l 2 . The distance between k n and k n+1 depends on n and it is small for large l . Therefore, at low-field regime, one can change the summation over n by an integral over k n : The derivative dn /dk n can be evaluated based on (30). The solution of (30) for n is given the following expression for n : Therefore, dn dk n = 2l 2 k n π L z 2l + (k n l ) 2 .
Based on (29), the derivative dn /dk n can be presented as follows: dn dk n = 2l 2 k n π −a n . 23702-8 Inserting (38) into (40), simultaneously taking into account (37) and transition (46), for the scattering time one obtains: Here, s = (l k n ) 2 is the dimensionless variable of integration, Brillouin zone is replaced by the infinite range: 0 k 2 ⊥ < ∞. The limits of integration over s are determined by the δ-function as: s min = −a 1 − L z /2l , s max = k 2 ⊥ − k 2 n l 2 , where a 1 = −(9π/8) 2/3 is the first zero of the Airy function [8]. δ-integration over s in (47) yields Here, s ⊥ = k 2 ⊥ l 2 is the dimensionless variable of integration, ∆ ≡ k 2 ⊥ l 2 + k 2 n l 2 = 2mεl 2 /ħ 2 is the dimensionless energy of an electron.
In (49), evaluations of the first and second the integrals over s ⊥ can be carried out based on (33) and following the integral of the product of two Airy functions [8], respectively, The result is as follows: where and τ −1 k,sc on electron dimensionless energy ε/k B T are plotted for n-Si at T = 300 K with the following parameters [5]: m = 0.32m 0 , ρ r = 2329 kg/m 3 , ν 0 = 8.43·10 3 m/s, D ac = 9 eV. The dependence τ −1 k,sc (ε/k B T ) (curve 'a') is calculated at F = 800 V/cm based on (61). The dependence τ −1 k,ac (ε/k B T ) (curve 'b') is calculated based on (1). As shown in figure 2  Taking into acount the problem solution reported in [6] it can be stated that the results of the present study can have a principial effect on the mobility fluctuation theory, especially. It should be noted that in the above-presented electron-acoustic phonon FIT scattering study, the electron quasi-momentum relaxation time and its relation to the scattering time is not included. It requires a separate consideration.