Una muestra de la teoría de perturbación k · p utilizada para calcular la estructura de múltiples bandas electrónicas degeneradas en semiconductores de pozos cuánticos y de volumen . 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 separada y sus contrapartes de giro) y dos bandas de conducción.Clase B : todas las demás bandas.El método se concentra en las bandas de la Clase A y tiene en cuenta las bandas de la Clase B de forma perturbadora.
Podemos escribir la solución perturbada como una combinación lineal de los estados propios no perturbados : ϕ {\ Displaystyle \ phi _ {} ^ {}} ϕ I ( 0 ) {\ Displaystyle \ phi _ {i} ^ {(0)}}
ϕ = ∑ norte A , B a norte ϕ norte ( 0 ) {\displaystyle \phi =\sum _{n}^{A,B}a_{n}\phi _{n}^{(0)}} Suponiendo que los estados propios no perturbados están ortonormalizados, las ecuaciones propias son:
( E − H m m ) a m = ∑ n ≠ m A H m n a n + ∑ α ≠ m B H m α a α {\displaystyle (E-H_{mm})a_{m}=\sum _{n\neq m}^{A}H_{mn}a_{n}+\sum _{\alpha \neq m}^{B}H_{m\alpha }a_{\alpha }} ,donde
H m n = ∫ ϕ m ( 0 ) † H ϕ n ( 0 ) d 3 r = E n ( 0 ) δ m n + H m n ′ {\displaystyle H_{mn}=\int \phi _{m}^{(0)\dagger }H\phi _{n}^{(0)}d^{3}\mathbf {r} =E_{n}^{(0)}\delta _{mn}+H_{mn}^{'}} .A partir de esta expresión podemos escribir:
a m = ∑ n ≠ m A H m n E − H m m a n + ∑ α ≠ m B H m α E − H m m a α {\displaystyle a_{m}=\sum _{n\neq m}^{A}{\frac {H_{mn}}{E-H_{mm}}}a_{n}+\sum _{\alpha \neq m}^{B}{\frac {H_{m\alpha }}{E-H_{mm}}}a_{\alpha }} ,donde la primera suma en el lado derecho está sobre los estados de la clase A solamente, mientras que la segunda suma está sobre los estados de la clase B. Como estamos interesados en los coeficientes de m en la clase A, podemos eliminar los de la clase B mediante un procedimiento de iteración para obtener: a m {\displaystyle a_{m}}
a m = ∑ n A U m n A − δ m n H m n E − H m m a n {\displaystyle a_{m}=\sum _{n}^{A}{\frac {U_{mn}^{A}-\delta _{mn}H_{mn}}{E-H_{mm}}}a_{n}} , U m n A = H m n + ∑ α ≠ m B H m α H α n E − H α α + ∑ α , β ≠ m , n ; α ≠ β H m α H α β H β n ( E − H α α ) ( E − H β β ) + … {\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 } De manera equivalente, para ( ): a n {\displaystyle a_{n}} n ∈ A {\displaystyle n\in A}
a n = ∑ n A ( U m n A − E δ m n ) a n = 0 , m ∈ A {\displaystyle a_{n}=\sum _{n}^{A}(U_{mn}^{A}-E\delta _{mn})a_{n}=0,m\in A} y
a γ = ∑ n A U γ n A − H γ n δ γ n E − H γ γ a n = 0 , γ ∈ B {\displaystyle a_{\gamma }=\sum _{n}^{A}{\frac {U_{\gamma n}^{A}-H_{\gamma n}\delta _{\gamma n}}{E-H_{\gamma \gamma }}}a_{n}=0,\gamma \in B} .Cuando se determinan los coeficientes pertenecientes a la Clase A, también lo son . a n {\displaystyle a_{n}} a γ {\displaystyle a_{\gamma }}
Ecuación de Schrödinger y funciones básicas El hamiltoniano, incluida la interacción espín-órbita, se puede escribir como:
H = H 0 + ℏ 4 m 0 2 c 2 σ ¯ ⋅ ∇ V × p {\displaystyle H=H_{0}+{\frac {\hbar }{4m_{0}^{2}c^{2}}}{\bar {\sigma }}\cdot \nabla V\times \mathbf {p} } ,donde está el vector de matriz de espín de Pauli . Sustituyendo en la ecuación de Schrödinger obtenemos σ ¯ {\displaystyle {\bar {\sigma }}}
H u n k ( r ) = ( H 0 + ℏ m 0 k ⋅ Π + ℏ 2 k 2 4 m 0 2 c 2 ∇ V × p ⋅ σ ¯ ) u n k ( r ) = E n ( k ) u n k ( r ) {\displaystyle Hu_{n\mathbf {k} }(\mathbf {r} )=\left(H_{0}+{\frac {\hbar }{m_{0}}}\mathbf {k} \cdot \mathbf {\Pi } +{\frac {\hbar ^{2}k^{2}}{4m_{0}^{2}c^{2}}}\nabla V\times \mathbf {p} \cdot {\bar {\sigma }}\right)u_{n\mathbf {k} }(\mathbf {r} )=E_{n}(\mathbf {k} )u_{n\mathbf {k} }(\mathbf {r} )} ,donde
Π = p + ℏ 4 m 0 2 c 2 σ ¯ × ∇ V {\displaystyle \mathbf {\Pi } =\mathbf {p} +{\frac {\hbar }{4m_{0}^{2}c^{2}}}{\bar {\sigma }}\times \nabla V} y la perturbación hamiltoniana se puede definir como
H ′ = ℏ m 0 k ⋅ Π . {\displaystyle H'={\frac {\hbar }{m_{0}}}\mathbf {k} \cdot \mathbf {\Pi } .} El hamiltoniano imperturbable 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 una simetría similar a la s, mientras que los estados de la banda de valencia son similares a p (degenerados 3 veces sin espín). Denotemos estos estados como , y , y respectivamente. Estas funciones de Bloch se pueden representar como una repetición periódica de orbitales atómicos, repetidos a intervalos correspondientes al espaciado de la red. La función Bloch se puede ampliar de la siguiente manera: | S ⟩ {\displaystyle |S\rangle } | X ⟩ {\displaystyle |X\rangle } | Y ⟩ {\displaystyle |Y\rangle } | Z ⟩ {\displaystyle |Z\rangle }
u n k ( r ) = ∑ j ′ A a j ′ ( k ) u j ′ 0 ( r ) + ∑ γ B a γ ( k ) u γ 0 ( r ) {\displaystyle u_{n\mathbf {k} }(\mathbf {r} )=\sum _{j'}^{A}a_{j'}(\mathbf {k} )u_{j'0}(\mathbf {r} )+\sum _{\gamma }^{B}a_{\gamma }(\mathbf {k} )u_{\gamma 0}(\mathbf {r} )} ,donde j ' está en la Clase A y está en la Clase B. Las funciones base se pueden elegir para ser γ {\displaystyle \gamma }
u 10 ( r ) = u e l ( r ) = | S 1 2 , 1 2 ⟩ = | S ↑ ⟩ {\displaystyle u_{10}(\mathbf {r} )=u_{el}(\mathbf {r} )=\left|S{\frac {1}{2}},{\frac {1}{2}}\right\rangle =\left|S\uparrow \right\rangle } u 20 ( r ) = u S O ( r ) = | 1 2 , 1 2 ⟩ = 1 3 | ( X + i Y ) ↓ ⟩ + 1 3 | Z ↑ ⟩ {\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 } u 30 ( r ) = u l h ( r ) = | 3 2 , 1 2 ⟩ = − 1 6 | ( X + i Y ) ↓ ⟩ + 2 3 | Z ↑ ⟩ {\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 } u 40 ( r ) = u h h ( r ) = | 3 2 , 3 2 ⟩ = − 1 2 | ( X + i Y ) ↑ ⟩ {\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 } u 50 ( r ) = u ¯ e l ( r ) = | S 1 2 , − 1 2 ⟩ = − | S ↓ ⟩ {\displaystyle u_{50}(\mathbf {r} )={\bar {u}}_{el}(\mathbf {r} )=\left|S{\frac {1}{2}},-{\frac {1}{2}}\right\rangle =-|S\downarrow \rangle } u 60 ( r ) = u ¯ S O ( r ) = | 1 2 , − 1 2 ⟩ = 1 3 | ( X − i Y ) ↑ ⟩ − 1 3 | Z ↓ ⟩ {\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 } u 70 ( r ) = u ¯ l h ( r ) = | 3 2 , − 1 2 ⟩ = 1 6 | ( X − i Y ) ↑ ⟩ + 2 3 | Z ↓ ⟩ {\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 } u 80 ( r ) = u ¯ h h ( r ) = | 3 2 , − 3 2 ⟩ = − 1 2 | ( X − i Y ) ↓ ⟩ {\displaystyle u_{80}(\mathbf {r} )={\bar {u}}_{hh}(\mathbf {r} )=\left|{\frac {3}{2}},-{\frac {3}{2}}\right\rangle =-{\frac {1}{\sqrt {2}}}|(X-iY)\downarrow \rangle } .Utilizando el método de Löwdin, solo es necesario resolver el siguiente problema de valores propios
∑ j ′ A ( U j j ′ A − E δ j j ′ ) a j ′ ( k ) = 0 , {\displaystyle \sum _{j'}^{A}(U_{jj'}^{A}-E\delta _{jj'})a_{j'}(\mathbf {k} )=0,} donde
U j j ′ A = H j j ′ + ∑ γ ≠ j , j ′ B H j γ H γ j ′ E 0 − E γ = H j j ′ + ∑ γ ≠ j , j ′ B H j γ ′ H γ j ′ ′ E 0 − E γ {\displaystyle U_{jj'}^{A}=H_{jj'}+\sum _{\gamma \neq j,j'}^{B}{\frac {H_{j\gamma }H_{\gamma j'}}{E_{0}-E_{\gamma }}}=H_{jj'}+\sum _{\gamma \neq j,j'}^{B}{\frac {H_{j\gamma }^{'}H_{\gamma j'}^{'}}{E_{0}-E_{\gamma }}}} , H j γ ′ = ⟨ u j 0 | ℏ m 0 k ⋅ ( p + ℏ 4 m 0 c 2 σ ¯ × ∇ V ) | u γ 0 ⟩ ≈ ∑ α ℏ k α m 0 p j γ α . {\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 }.} 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 una sola banda, podemos escribir para Π {\displaystyle \Pi } U j j ′ A {\displaystyle U_{jj'}^{A}}
D j j ′ ≡ U j j ′ A = E j ( 0 ) δ j j ′ + ∑ α β D j j ′ α β k α k β , {\displaystyle D_{jj'}\equiv U_{jj'}^{A}=E_{j}(0)\delta _{jj'}+\sum _{\alpha \beta }D_{jj'}^{\alpha \beta }k_{\alpha }k_{\beta },} D j j ′ α β = ℏ 2 2 m 0 [ δ j j ′ δ α β + ∑ γ B p j γ α p γ j ′ β + p j γ β p γ j ′ α m 0 ( E 0 − E γ ) ] . {\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].} Ahora definimos los siguientes parámetros
A 0 = ℏ 2 2 m 0 + ℏ 2 m 0 2 ∑ γ B p x γ x p γ x x E 0 − E γ , {\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 }}},} B 0 = ℏ 2 2 m 0 + ℏ 2 m 0 2 ∑ γ B p x γ y p γ x y E 0 − E γ , {\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 }}},} C 0 = ℏ 2 m 0 2 ∑ γ B p x γ x p γ y y + p x γ y p γ y x E 0 − E γ , {\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 }}},} y los parámetros de la estructura de la banda (o los parámetros de Luttinger ) se pueden definir para ser
γ 1 = − 1 3 2 m 0 ℏ 2 ( A 0 + 2 B 0 ) , {\displaystyle \gamma _{1}=-{\frac {1}{3}}{\frac {2m_{0}}{\hbar ^{2}}}(A_{0}+2B_{0}),} γ 2 = − 1 6 2 m 0 ℏ 2 ( A 0 − B 0 ) , {\displaystyle \gamma _{2}=-{\frac {1}{6}}{\frac {2m_{0}}{\hbar ^{2}}}(A_{0}-B_{0}),} γ 3 = − 1 6 2 m 0 ℏ 2 C 0 , {\displaystyle \gamma _{3}=-{\frac {1}{6}}{\frac {2m_{0}}{\hbar ^{2}}}C_{0},} 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 . γ 1 {\displaystyle \gamma _{1}} γ 2 {\displaystyle \gamma _{2}} | X ⟩ {\displaystyle |X\rangle } | Y ⟩ {\displaystyle |Y\rangle } | Z ⟩ {\displaystyle |Z\rangle } γ 3 {\displaystyle \gamma _{3}} Γ {\displaystyle \Gamma } γ 2 ≠ γ 3 {\displaystyle \gamma _{2}\neq \gamma _{3}}
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 luces y 2 de escisión) D j j ′ {\displaystyle \mathbf {D_{jj'}} }
H = ( E e l P z 2 P z − 3 P + 0 2 P − P − 0 P z † P + Δ 2 Q † − S † / 2 − 2 P + † 0 − 3 / 2 S − 2 R E e l P z 2 P z − 3 P + 0 2 P − P − 0 E e l P z 2 P z − 3 P + 0 2 P − P − 0 E e l P z 2 P z − 3 P + 0 2 P − P − 0 E e l P z 2 P z − 3 P + 0 2 P − P − 0 E e l P z 2 P z − 3 P + 0 2 P − P − 0 E e l P z 2 P z − 3 P + 0 2 P − P − 0 ) {\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)} Resumen Esta sección está vacía. Puedes ayudar
agregando más .
( Julio de 2010 )
Referencias