k·p perturbation theory
In solid-state physics, the k·p perturbation theory is an approximated semi-empirical approach for calculating the band structure (particularly effective mass) and optical properties of crystalline solids. It is pronounced "k dot p", and is also called the "k·p method". This theory has been applied specifically in the framework of the Luttinger–Kohn model (after Joaquin Mazdak Luttinger and Walter Kohn), and of the Kane model (after Evan O. Kane).
Background and derivation
Bloch's theorem and wavevectors
According to quantum mechanics (in the single-electron approximation), the quasi-free electrons in any solid are characterized by wavefunctions which are eigenstates of the following stationary Schrödinger equation:
where k is a vector (called the wavevector), n is a discrete index (called the band index), and un,k is a function with the same periodicity as the crystal lattice.
For any given n, the associated states are called a band. In each band, there will be a relation between the wavevector k and the energy of the state En,k, called the band dispersion. Calculating this dispersion is one of the primary applications of k·p perturbation theory.
The periodic function un,k satisfies the following Schrödinger-type equation:
where the Hamiltonian is
Note that k is a vector consisting of three real numbers with dimensions of inverse length, while p is a vector of operators; to be explicit,
In any case, we write this Hamiltonian as the sum of two terms:
This expression is the basis for perturbation theory. The "unperturbed Hamiltonian" is H0, which in fact equals the exact Hamiltonian at k=0 (i.e., at the Gamma point). The "perturbation" is the term . The analysis that results is called "k·p perturbation theory", due to the term proportional to k·p. The result of this analysis is an expression for En,k and un,k in terms of the energies and wavefunctions at k=0.
Note that the "perturbation" term gets progressively smaller as k approaches zero. Therefore, k·p perturbation theory is most accurate for small values of k. However, if enough terms are included in the perturbative expansion, then the theory can in fact be reasonably accurate for any value of k in the entire Brillouin zone.
Expression for a nondegenerate band
For a nondegenerate band (i.e., a band which has a different energy at k=0 from any other band), with an extremum at k=0, and with no spin–orbit coupling, the result of k·p perturbation theory is (to lowest nontrivial order):
Since k is a vector of real numbers (rather than a vector of more complicated linear operators), the matrix element in these expressions can be rewritten as:
Therefore, you can calculate the energy at any k using only a few unknown parameters, namely En,0 and . The latter are called "optical matrix elements", closely related to transition dipole moments. These parameters are typically inferred from experimental data.
In practice, the sum over n often includes only the nearest one or two bands, since these tend to be the most important (due to the denominator). However, for improved accuracy, especially at larger k, more bands must be included, as well as more terms in the perturbative expansion than the ones written above.
Using the expression above for the energy dispersion relation, a simplified expression for the effective mass in the conduction band of a semiconductor can be found. To approximate the dispersion relation in the case of the conduction band, take the energy En0 as the minimum conduction band energy Ec0 and include in the summation only terms with energies near the valence band maximum, where the energy difference in the denominator is smallest. (These terms are the largest contributions to the summation.) This denominator is then approximated as the band gap Eg, leading to an energy expression:
The effective mass in direction ℓ is then:
Ignoring the details of the matrix elements, the key consequences are that the effective mass varies with the smallest bandgap and goes to zero as the gap goes to zero. A useful approximation for the matrix elements in direct gap semiconductors is:
which applies within about 15% or better to most group-IV, III-V and II-VI semiconductors.
In contrast to this simple approximation, in the case of valence band energy the spin–orbit interaction must be introduced (see below) and many more bands must be individually considered. The calculation is provided in Yu and Cardona. In the valence band the mobile carriers are holes. One finds there are two types of hole, named heavy and light, with anisotropic masses.
k·p model with spin–orbit interaction
where is a vector consisting of the three Pauli matrices. This Hamiltonian can be subjected to the same sort of perturbation-theory analysis as above.
Calculation in degenerate case
For degenerate or nearly degenerate bands, in particular the valence bands in certain materials such as gallium arsenide, the equations can be analyzed by the methods of degenerate perturbation theory. Models of this type include the "Luttinger–Kohn model" (a.k.a. "Kohn–Luttinger model"), and the "Kane model".
Electronic band structure
Notes and references
- P. Yu, M. Cardona (2005). Fundamentals of Semiconductors: Physics and Materials Properties (3rd ed.). Springer. Section 2.6, pp. 68 ff'. ISBN 3-540-25470-6.
- C. Kittel (1987). Quantum Theory of Solids (Second Revised Printing ed.). New York: Wiley. pp. 186–190. ISBN 0-471-62412-8.
- W. P. Harrison (1989) . Electronic Structure and the Properties of Solids (Reprint ed.). Dover Publications. pp. 158 ff. ISBN 0-486-66021-4.
- A direct gap semiconductor is one where the valence band maximum and conduction band minimum occur at the same position in k-space, usually the so-called Γ-point where k = 0.
- See Table 2.22 in Yu & Cardona, op. cit.
- See Yu & Cardona, op. cit. pp. 75–82
- Evan O. Kane (1957). "Band Structure of Indium Antimonide". Journal of Physics and Chemistry of Solids. 1: 249. Bibcode:1957JPCS....1..249K. doi:10.1016/0022-3697(57)90013-6.
- J. M. Luttinger, W. Kohn (1955). "Motion of Electrons and Holes in Perturbed Periodic Fields". Physical Review. 97: 869. Bibcode:1955PhRv...97..869L. doi:10.1103/PhysRev.97.869.