A flavor of the k·p perturbation theory used for calculating the structure of multiple, degenerate electronic bands in bulk and quantum well semiconductors. The method is a generalization of the single band k·p theory.
In this model the influence of all other bands is taken into account by using Löwdin's perturbation method.[1]
Background
All bands can be subdivided into two classes:
- Class A: six valence bands (heavy hole, light hole, split off band and their spin counterparts) and two conduction bands.
- Class B: all other bands.
The method concentrates on the bands in Class A, and takes into account Class B bands perturbatively.
We can write the perturbed solution
as a linear combination of the unperturbed eigenstates
:
![{\displaystyle \phi =\sum _{n}^{A,B}a_{n}\phi _{n}^{(0)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a9ecfb1344f261dcb1cfd82268e5f1f10747ab96)
Assuming the unperturbed eigenstates are orthonormalized, the eigenequation are:
,
where
.
From this expression we can write:
,
where the first sum on the right-hand side is over the states in class A only, while the second sum is over the states on class B. Since we are interested in the coefficients
for m in class A, we may eliminate those in class B by an iteration procedure to obtain:
,
![{\displaystyle U_{mn}^{A}=H_{mn}+\sum _{\alpha \neq m}^{B}{\frac {H_{m\alpha }H_{\alpha n}}{E-H_{\alpha \alpha }}}+\sum _{\alpha ,\beta \neq m,n;\alpha \neq \beta }{\frac {H_{m\alpha }H_{\alpha \beta }H_{\beta n}}{(E-H_{\alpha \alpha })(E-H_{\beta \beta })}}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b53728672be3201e74252f679be36960543efde)
Equivalently, for
(
):
![{\displaystyle a_{n}=\sum _{n}^{A}(U_{mn}^{A}-E\delta _{mn})a_{n}=0,m\in A}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55da82e280b4704b0f1b8fbf08bb1a9dda9fe59e)
and
.
When the coefficients
belonging to Class A are determined so are
.
Schrödinger equation and basis functions
The Hamiltonian including the spin-orbit interaction can be written as:
,
where
is the Pauli spin matrix vector. Substituting into the Schrödinger equation in Bloch approximation we obtain
,
where
![{\displaystyle \mathbf {\Pi } =\mathbf {p} +{\frac {\hbar }{4m_{0}^{2}c^{2}}}{\bar {\sigma }}\times \nabla V}](https://wikimedia.org/api/rest_v1/media/math/render/svg/64c966e661f24a201effa5d6c33228903a762fd6)
and the perturbation Hamiltonian can be defined as
![{\displaystyle H'={\frac {\hbar }{m_{0}}}\mathbf {k} \cdot \mathbf {\Pi } .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49a766edf4f797692403fdd2c609562002714b53)
The unperturbed Hamiltonian refers to the band-edge spin-orbit system (for k=0). At the band edge, conduction band Bloch waves exhibit s-like symmetry, while valence band states are p-like (3-fold degenerate without spin). Let us denote these states as
, and
,
and
respectively. These Bloch functions can be pictured as periodic repetition of atomic orbitals, repeated at intervals corresponding to the lattice spacing. The Bloch function can be expanded in the following manner:
,
where j' is in Class A and
is in Class B. The basis functions can be chosen to be
![{\displaystyle u_{10}(\mathbf {r} )=u_{el}(\mathbf {r} )=\left|S{\frac {1}{2}},{\frac {1}{2}}\right\rangle =\left|S\uparrow \right\rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec06dbc6a1dcbc56e4103145b36d86d2938dcd92)
![{\displaystyle u_{20}(\mathbf {r} )=u_{SO}(\mathbf {r} )=\left|{\frac {1}{2}},{\frac {1}{2}}\right\rangle ={\frac {1}{\sqrt {3}}}|(X+iY)\downarrow \rangle +{\frac {1}{\sqrt {3}}}|Z\uparrow \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/44d3e479b040d22e6118e96e0a04fd44b2056c4c)
![{\displaystyle u_{30}(\mathbf {r} )=u_{lh}(\mathbf {r} )=\left|{\frac {3}{2}},{\frac {1}{2}}\right\rangle =-{\frac {1}{\sqrt {6}}}|(X+iY)\downarrow \rangle +{\sqrt {\frac {2}{3}}}|Z\uparrow \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/c8b75e7d7af530f20c6ffb1bfd0d7c5bd3ec9a42)
![{\displaystyle u_{40}(\mathbf {r} )=u_{hh}(\mathbf {r} )=\left|{\frac {3}{2}},{\frac {3}{2}}\right\rangle =-{\frac {1}{\sqrt {2}}}|(X+iY)\uparrow \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/6bf3390c6aef9c628a534929e9c87d91a76ac9e9)
![{\displaystyle u_{50}(\mathbf {r} )={\bar {u}}_{el}(\mathbf {r} )=\left|S{\frac {1}{2}},-{\frac {1}{2}}\right\rangle =-|S\downarrow \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f2b130262637b63da14dfe34dec45a2b565fbb2)
![{\displaystyle u_{60}(\mathbf {r} )={\bar {u}}_{SO}(\mathbf {r} )=\left|{\frac {1}{2}},-{\frac {1}{2}}\right\rangle ={\frac {1}{\sqrt {3}}}|(X-iY)\uparrow \rangle -{\frac {1}{\sqrt {3}}}|Z\downarrow \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/3465b331c9caee83d17333eef0986f63c76f687e)
![{\displaystyle u_{70}(\mathbf {r} )={\bar {u}}_{lh}(\mathbf {r} )=\left|{\frac {3}{2}},-{\frac {1}{2}}\right\rangle ={\frac {1}{\sqrt {6}}}|(X-iY)\uparrow \rangle +{\sqrt {\frac {2}{3}}}|Z\downarrow \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd04b6ded6e9b532fcb1a79859e5d185e1604245)
.
Using Löwdin's method, only the following eigenvalue problem needs to be solved
![{\displaystyle \sum _{j'}^{A}(U_{jj'}^{A}-E\delta _{jj'})a_{j'}(\mathbf {k} )=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/61cb1cde970b03a1d63d849ec416ef426fdec706)
where
,
![{\displaystyle H_{j\gamma }^{'}=\left\langle u_{j0}\right|{\frac {\hbar }{m_{0}}}\mathbf {k} \cdot \left(\mathbf {p} +{\frac {\hbar }{4m_{0}c^{2}}}{\bar {\sigma }}\times \nabla V\right)\left|u_{\gamma 0}\right\rangle \approx \sum _{\alpha }{\frac {\hbar k_{\alpha }}{m_{0}}}p_{j\gamma }^{\alpha }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a8213d5fd7f6f023ab1b1bbacb292ea0caf6d4d)
The second term of
can be neglected compared to the similar term with p instead of k. Similarly to the single band case, we can write for
![{\displaystyle D_{jj'}\equiv U_{jj'}^{A}=E_{j}(0)\delta _{jj'}+\sum _{\alpha \beta }D_{jj'}^{\alpha \beta }k_{\alpha }k_{\beta },}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ab32e22f9986e3a5fbf86ebc7666d580398fc10)
![{\displaystyle D_{jj'}^{\alpha \beta }={\frac {\hbar ^{2}}{2m_{0}}}\left[\delta _{jj'}\delta _{\alpha \beta }+\sum _{\gamma }^{B}{\frac {p_{j\gamma }^{\alpha }p_{\gamma j'}^{\beta }+p_{j\gamma }^{\beta }p_{\gamma j'}^{\alpha }}{m_{0}(E_{0}-E_{\gamma })}}\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad2027289efc750ad8091e5570a7c3d76899068d)
We now define the following parameters
![{\displaystyle A_{0}={\frac {\hbar ^{2}}{2m_{0}}}+{\frac {\hbar ^{2}}{m_{0}^{2}}}\sum _{\gamma }^{B}{\frac {p_{x\gamma }^{x}p_{\gamma x}^{x}}{E_{0}-E_{\gamma }}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0146310c1bc4e9b88ff67b1af1233212a7c00106)
![{\displaystyle B_{0}={\frac {\hbar ^{2}}{2m_{0}}}+{\frac {\hbar ^{2}}{m_{0}^{2}}}\sum _{\gamma }^{B}{\frac {p_{x\gamma }^{y}p_{\gamma x}^{y}}{E_{0}-E_{\gamma }}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e92b69b8b9b4e992f122c2509968adf245c4ea5)
![{\displaystyle C_{0}={\frac {\hbar ^{2}}{m_{0}^{2}}}\sum _{\gamma }^{B}{\frac {p_{x\gamma }^{x}p_{\gamma y}^{y}+p_{x\gamma }^{y}p_{\gamma y}^{x}}{E_{0}-E_{\gamma }}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2890670679d99a6c081d82deb419ab4fd0b373f1)
and the band structure parameters (or the Luttinger parameters) can be defined to be
![{\displaystyle \gamma _{1}=-{\frac {1}{3}}{\frac {2m_{0}}{\hbar ^{2}}}(A_{0}+2B_{0}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/401f8dbc9d190a3e29d68412eeb89268cbdde17a)
![{\displaystyle \gamma _{2}=-{\frac {1}{6}}{\frac {2m_{0}}{\hbar ^{2}}}(A_{0}-B_{0}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e2372c6536d34d0d814df128b5eed38aab163488)
![{\displaystyle \gamma _{3}=-{\frac {1}{6}}{\frac {2m_{0}}{\hbar ^{2}}}C_{0},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0ea1b30a93d22742c31825f1de52856618b059b9)
These parameters are very closely related to the effective masses of the holes in various valence bands.
and
describe the coupling of the
,
and
states to the other states. The third parameter
relates to the anisotropy of the energy band structure around the
point when
.
Explicit Hamiltonian matrix
The Luttinger-Kohn Hamiltonian
can be written explicitly as a 8X8 matrix (taking into account 8 bands - 2 conduction, 2 heavy-holes, 2 light-holes and 2 split-off)
![{\displaystyle \mathbf {H} =\left({\begin{array}{cccccccc}E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\P_{z}^{\dagger }&P+\Delta &{\sqrt {2}}Q^{\dagger }&-S^{\dagger }/{\sqrt {2}}&-{\sqrt {2}}P_{+}^{\dagger }&0&-{\sqrt {3/2}}S&-{\sqrt {2}}R\\E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\E_{el}&P_{z}&{\sqrt {2}}P_{z}&-{\sqrt {3}}P_{+}&0&{\sqrt {2}}P_{-}&P_{-}&0\\\end{array}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e67941bfa1b5d464564a18adac792dc3813aa42a)
Summary
![[icon]](//upload.wikimedia.org/wikipedia/commons/thumb/1/1c/Wiki_letter_w_cropped.svg/20px-Wiki_letter_w_cropped.svg.png) | This section is empty. You can help by adding to it. (July 2010) |
References
2. Luttinger, J. M. Kohn, W., "Motion of Electrons and Holes in Perturbed Periodic Fields", Phys. Rev. 97,4. pp. 869-883, (1955). https://journals.aps.org/pr/abstract/10.1103/PhysRev.97.869