El modelo de Luttinger-Kohn es una versión de la teoría de perturbación k·p utilizada para calcular la estructura de múltiples bandas electrónicas degeneradas en semiconductores de pozo cuántico y masivo . El método es una generalización de la teoría k · p de banda única .
En este modelo, la influencia de todas las demás bandas se tiene en cuenta mediante el método de perturbación de Löwdin . [1]
Fondo
Todas las bandas se pueden subdividir en dos clases:
- Clase A : seis bandas de valencia (agujero pesado, agujero ligero, banda dividida y sus contrapartes de espín) y dos bandas de conducción.
- Clase B : todas las demás bandas.
El método se concentra en las bandas de Clase A y tiene en cuenta las bandas de Clase B de forma perturbativa.
Podemos escribir la solución perturbada, como una combinación lineal de los estados propios no perturbados :![{\displaystyle \phi _{}^{}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \phi _{i}^{(0)}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \phi =\sum _ {n}^{A,B}a_ {n}\phi _ {n}^{(0)}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Suponiendo que los estados propios no perturbados estén ortonormalizados, las ecuaciones propias son:
,
dónde
.
De esta expresión podemos escribir:
,
donde la primera suma del lado derecho cubre los estados de la clase A únicamente, mientras que la segunda suma cubre los estados de la clase B. Como estamos interesados en los coeficientes para m en la clase A, podemos eliminar los de la clase B mediante un procedimiento de iteración para obtener:![{\ Displaystyle a_ {m}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
,
![{\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 }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
De manera equivalente, para ( ):![{\ Displaystyle a_ {n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle n\en A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle a_{n}=\sum _{n}^{A}(U_{mn}^{A}-E\delta _{mn})a_{n}=0,m\in A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y
.
Cuando se determinan los coeficientes pertenecientes a la Clase A, también se determinan .![{\ Displaystyle a_ {n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle a _ {\ gamma}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ecuación de Schrödinger y funciones básicas
El hamiltoniano , incluida la interacción espín-órbita, se puede escribir como:
,
¿Dónde está el vector de matriz de espín de Pauli ? Sustituyendo en la ecuación de Schrödinger en la aproximación de Bloch obtenemos![{\displaystyle {\bar {\sigma }}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
,
dónde
![{\displaystyle \mathbf {\Pi } =\mathbf {p} +{\frac {\hbar }{4m_{0}^{2}c^{2}}}{\bar {\sigma }}\times \ nabla V}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y la perturbación hamiltoniana se puede definir como
![{\displaystyle H'={\frac {\hbar }{m_{0}}}\mathbf {k} \cdot \mathbf {\Pi } .}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El hamiltoniano imperturbado se refiere al sistema de órbita de giro de borde de banda (para k = 0). En el borde de la banda, las ondas de Bloch de la banda de conducción exhiben simetría tipo s, mientras que los estados de la banda de valencia son tipo p (3 veces degenerados sin espín). Denotemos estos estados como , y , y respectivamente. Estas funciones de Bloch pueden representarse como repeticiones periódicas de orbitales atómicos, repetidas a intervalos correspondientes al espaciamiento de la red. La función Bloch se puede ampliar de la siguiente manera:![{\displaystyle |S\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |X\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |Y\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |Z\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
,
donde j' está en la Clase A y está en la Clase B. Las funciones base se pueden elegir para que sean![{\displaystyle \gamma}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u_{10}(\mathbf {r} )=u_{el}(\mathbf {r} )=\left|S{\frac {1}{2}},{\frac {1}{2 }}\right\rangle =\left|S\uparrow \right\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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 \ rango}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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 }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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 }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u_{50}(\mathbf {r} )={\bar {u}}_{el}(\mathbf {r} )=\left|S{\frac {1}{2}},- {\frac {1}{2}}\right\rangle =-|S\downarrow \rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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 }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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 }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
.
Usando el método de Löwdin, solo es necesario resolver el siguiente problema de valores propios
![{\displaystyle \sum _{j'}^{A}(U_{jj'}^{A}-E\delta _{jj'})a_{j'}(\mathbf {k} )=0,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde
,
![{\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 }.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El segundo término de puede despreciarse en comparación con el término similar con p en lugar de k . De manera similar al caso de banda única, podemos escribir para![{\displaystyle \Pi }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle U_{jj'}^{A}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle D_{jj'}\equiv U_{jj'}^{A}=E_{j}(0)\delta _{jj'}+\sum _{\alpha \beta }D_{jj'}^ {\alpha \beta }k_{\alpha }k_{\beta },}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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].}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ahora definimos los siguientes parámetros.
![{\displaystyle A_{0}={\frac {\hbar ^{2}}{2m_{0}}}+{\frac {\hbar ^{2}}{m_{0}^{2}}}\ suma _{\gamma }^{B}{\frac {p_{x\gamma }^{x}p_{\gamma x}^{x}}{E_{0}-E_{\gamma }}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle B_{0}={\frac {\hbar ^{2}}{2m_{0}}}+{\frac {\hbar ^{2}}{m_{0}^{2}}}\ suma _{\gamma }^{B}{\frac {p_{x\gamma }^{y}p_{\gamma x}^{y}}{E_{0}-E_{\gamma }}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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 }}} ,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y los parámetros de estructura de banda (o los parámetros de Luttinger ) se pueden definir para ser
![{\displaystyle \gamma _{1}=-{\frac {1}{3}}{\frac {2m_{0}}{\hbar ^{2}}}(A_{0}+2B_{0}) ,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \gamma _{2}=-{\frac {1}{6}}{\frac {2m_{0}}{\hbar ^{2}}}(A_{0}-B_{0}) ,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \gamma _{3}=-{\frac {1}{6}}{\frac {2m_{0}}{\hbar ^{2}}}C_{0},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Estos parámetros están muy relacionados con las masas efectivas de los huecos en varias bandas de valencia. y describir el acoplamiento de los estados , y con los otros estados. El tercer parámetro se relaciona con la anisotropía de la estructura de la banda de energía alrededor del punto cuando .![{\displaystyle \gamma _{1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \gamma _{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |X\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |Y\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |Z\rangle }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \gamma _{3}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Gamma}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \gamma _{2}\neq \gamma _{3}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Matriz hamiltoniana explícita
El hamiltoniano de Luttinger-Kohn se puede escribir explícitamente como una matriz de 8X8 (teniendo en cuenta 8 bandas: 2 de conducción, 2 de agujeros pesados, 2 de agujeros ligeros y 2 de separación)![{\displaystyle \mathbf {D_{jj'}} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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^{\daga }/{\sqrt {2}}&-{\sqrt {2}}P_{+}^{\daga }&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)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Resumen
Referencias
- ^ SL Chuang (1995). Física de dispositivos optoelectrónicos (Primera ed.). Nueva York: Wiley. págs. 124-190. ISBN 978-0-471-10939-6. OCLC 31134252.
2. Luttinger, JM Kohn, W., "Movimiento de electrones y agujeros en campos periódicos perturbados", Phys. Apocalipsis 97,4. págs. 869-883, (1955). https://journals.aps.org/pr/abstract/10.1103/PhysRev.97.869