Pozo cuántico de GaN/AlGaN: Tutorial de 8 bandas sobre la estructura de Wurtzite k·p.
1. Introducción
En este tutorial, calcularemos los estados electrónicos y la respuesta óptica de un pozo cuántico de GaN situado entre dos barreras de AlGaN, utilizando el modelo de ocho bandas de wurtzita k·p de OghmaNano. El ejemplo incluye un pozo de GaN de 10 nm y dos capas de AlGaN de 20 nm. Primero, ejecutaremos el ejemplo proporcionado y, a continuación, analizaremos la dispersión de las subbandas de conducción y valencia, así como los espectros ópticos resultantes.
La diferencia entre la estructura cristalina de wurtzita y la de zincblenda es importante en este caso. Ambas son estructuras cristalinas semiconductoras, pero la wurtzita presenta una orientación hexagonal preferente a lo largo del eje c. La diferente simetría cristalina modifica la división de la banda de valencia, el acoplamiento entre los estados electrónicos y la intensidad de la luz polarizada en diferentes direcciones. En consecuencia, no se puede simplemente reutilizar un cálculo de ocho bandas basado en la estructura de zincblenda con los parámetros de los materiales de GaN.
2. Abrir el ejemplo de wurtzita con ocho bandas.
En la pestaña «Archivo», haga clic en «Nueva simulación». En la ventana que aparece, haga doble clic en «Ejemplos de GaAs» (??). A pesar del nombre de la categoría, también contiene ejemplos de heterouniones cuánticas de nitruro y silicio. En la siguiente ventana, haga doble clic en «GaN Wurtzita 8 bandas k.p» (??). Esto abrirá la pila de materiales preparada y seleccionará la demostración adecuada.
3. Inspección del pozo cuántico de GaN/AlGaN.
La ventana principal muestra la estructura de tres capas (??). En el cálculo real del pozo cuántico, el confinamiento se produce a lo largo del grosor de esta estructura: las funciones de envolvente electrónicas varían entre las capas, mientras que los diagramas de dispersión describen el movimiento paralelo a ellas.
Haga clic en «Editor de capas» en la barra de herramientas de la izquierda. Como se muestra en ??, el ejemplo utiliza la siguiente geometría:
| Capa | Grosor | Tipo |
|---|---|---|
| AlGaN | 20 nanómetros | Activo / barrera |
| GaN (heteroestructura de pozo cuántico) | 10 nanómetros | Pozo cuántico |
| AlGaN | 20 nanómetros | Activo / barrera |
Para su primera prueba, no modifique la configuración de los materiales proporcionados. La captura de pantalla muestra el grosor de las capas, pero no indica la fracción molar de aluminio, la configuración de la deformación ni los parámetros de polarización; estos parámetros deben verificarse en la configuración de los materiales del ejemplo antes de utilizarlo para representar una muestra experimental específica.
4. ¿Por qué utilizar un hamiltoniano de wurtzita con ocho bandas?
Un modelo parabólico sencillo asigna una masa efectiva a cada portador de carga y asume que su energía varía aproximadamente como el cuadrado de su vector de onda en el plano. Esto es útil cerca de un extremo de banda, pero no puede describir cómo los diferentes estados de valencia se mezclan entre sí a medida que aumenta el vector de onda. En los materiales de tipo wurtzita, el entorno cristalino hexagonal también distingue el movimiento y la polarización de la luz a lo largo del eje c de los que se dan en el plano perpendicular.
El modelo de ocho bandas conserva una órbita similar a la de conducción y tres órbitas similares a las de valencia, cada una con dos componentes de espín. Conceptualmente, un estado propio tiene la forma:
\[ \Psi_{n,\mathbf{k}_{\parallel}}(\mathbf{r})= e^{i\mathbf{k}_{\parallel}\cdot\mathbf{r}_{\parallel}} \sum_{\alpha=1}^{8}F_{n\alpha}(z)u_{\alpha}(\mathbf{r}), \]
donde las funciones \(F_{n\alpha}(z)\) describen el confinamiento a lo largo del pozo, y \(u_\alpha\) son los ocho estados base de Bloch. Por lo tanto, un estado confinado dado puede tener contribuciones de más de una orbital de la banda de valencia; por esta razón, describir cada curva coloreada utilizando una única etiqueta fija para agujeros pesados o ligeros puede ser engañoso, especialmente cuando nos alejamos del borde de la banda.
El ordenamiento de la banda de valencia de la wurtzita se ve afectado por las interacciones del campo cristalino y del acoplamiento espín-órbita. En una notación concisa, el cálculo resuelve:
\[ H_{8\times8}\!\left(k_x,k_y,-i\frac{\partial}{\partial z}\right) \mathbf{F}_n(z)=E_n(\mathbf{k}_{\parallel})\,\mathbf{F}_n(z). \]
Obsérvese que el vector de onda en el plano se mantiene explícitamente, mientras que el vector de onda en la dirección de crecimiento se sustituye por un operador diferencial. Al repetir este problema de autovalores para diferentes valores de \(k_{\parallel}\), se obtienen las curvas de dispersión de las subbandas que se muestran a continuación. Una explicación más detallada de la base, los términos del campo cristalino y los coeficientes del material se ofrece en el capítulo teoría de Wurtzite k·p.
5. Ejecute la simulación.
Vuelva a la pestaña «Archivo» y pulse «Ejecutar simulación». Abra la pestaña «Terminal» para seguir el progreso. En ??, la terminal identifica el solucionador WZ KP8, informa sobre el ensamblaje de la matriz dispersa, muestra información sobre la convergencia y avanza a través de los vectores de onda muestreados. Resuelve los estados electrónicos y de valencia por separado en cada punto.
Dado que se trata de un problema de función envolvente con múltiples componentes, es más complejo que resolver una ecuación parabólica escalar independiente para cada banda. Espere a que el solucionador de autovalores y los cálculos ópticos posteriores hayan finalizado antes de examinar los archivos resultantes.
5.1. Inspección de las intensidades de las transiciones ópticas.
Más adelante, durante el proceso, el terminal muestra la intensidad de los elementos de la matriz óptica (??). Para una estructura de wurtzita convencional en el plano c, TE se refiere al campo eléctrico óptico perpendicular al eje c, mientras que TM se refiere a su componente paralela a ese eje. Estas polarizaciones no tienen por qué acoplarse con la misma intensidad a una transición electrónica determinada.
Una forma útil de analizar la intensidad óptica es mediante el elemento de la matriz de transición elevado al cuadrado:
\[ M_{cv}^{(\mathbf{e})}(\mathbf{k}_{\parallel})= \left\langle\Psi_{c,\mathbf{k}_{\parallel}}\middle| \mathbf{e}\!\cdot\!\mathbf{p}\middle| \Psi_{v,\mathbf{k}_{\parallel}}\right\rangle, \qquad S_{cv}^{(\mathbf{e})}\propto\left|M_{cv}^{(\mathbf{e})}\right|^2. \]
Aquí, \(\mathbf{e}\) especifica la polarización de la luz y \(\mathbf{p}\) es el operador de momento. Un estado puede estar presente en la estructura de bandas electrónicas, pero contribuir muy poco a una polarización óptica específica si su elemento de matriz es pequeño. Por lo tanto, los valores relativos de TE/TM en el terminal proporcionan información física adicional más allá de las posiciones de los niveles de energía.
6. Comprensión de las bandas electrónicas calculadas.
El eje horizontal de cada diagrama de bandas representa el vector de onda en el plano \(k_{\parallel}\), y no la posición dentro del pozo cuántico. El eje vertical representa la energía de una subbanda confinada. En el punto \(k_{\parallel}=0\), cada curva comienza en una de las energías de confinamiento discretas. Al desplazarse hacia la derecha, se observa cómo se dispersa el mismo estado electrónico cuando la carga se mueve en paralelo a las capas.
6.1 Estados de la banda de conducción
Active la salida de la banda de conducción e inspeccione la dispersión en ??. Se pueden observar varios niveles confinados, que comienzan en diferentes energías cercanas al borde de la banda de GaN. Para valores pequeños del vector de onda, una aproximación útil para una subbanda de conducción es:
\[ E_{c,n}(k_{\parallel})\simeq E_{c,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_{e,\parallel}^{*}}. \]
El término \(E_{c,n}(0)\) contiene la energía de confinamiento, mientras que \(m_{e,\parallel}^{*}\) describe la respuesta del portador al movimiento en el plano. Esta relación ayuda a explicar las curvas que muestran un aumento gradual, aunque el cálculo completo con ocho bandas mantiene el acoplamiento interbanda en lugar de imponer una forma parabólica fija.
6.2 Estados de la banda de valencia.
El diagrama de la banda de valencia (??) es más complejo. Su rama superior se eleva inicialmente alejándose de \(k_{\parallel}=0\) antes de curvarse hacia abajo; varias ramas inferiores presentan curvaturas notablemente diferentes. Estas formas son posibles en un modelo de banda de valencia de wurtzita acoplado, en el que la división del campo cristalino, el acoplamiento espín-órbita, el confinamiento y la mezcla de bandas actúan en conjunto.
No se debe interpretar esta representación gráfica como prueba de que cada curva coloreada tenga una identidad fija en las bandas A, B o C en todo el rango del vector de onda. Para realizar una asignación definitiva, examine las contribuciones de los componentes del vector propio subyacente. De manera similar, los extremos inusuales que no se encuentran en el centro merecen ser comparados con el material y los parámetros de deformación elegidos, si se utiliza el modelo de forma cuantitativa.
7. Transiciones ópticas, emisión estimulada y absorción.
Las bandas electrónicas determinan las energías de las posibles transiciones ópticas, pero la intensidad de estas transiciones también depende de las funciones de onda de los electrones y las huecas, la polarización y la ocupación de los portadores. Para un vector de onda común en el plano, la energía del fotón para una transición vertical es:
\[ \hbar\omega=E_{c,n}(\mathbf{k}_{\parallel}) -E_{v,m}(\mathbf{k}_{\parallel}). \]
La regla de oro de Fermi explica por qué los elementos de la matriz óptica que se muestran en la pantalla son importantes. En forma esquemática, la tasa de una transición es:
\[ W_{v\rightarrow c}^{(\mathbf{e})}= \frac{2\pi}{\hbar}\left|H_{cv}^{\prime(\mathbf{e})}\right|^2 \delta\!\left(E_c-E_v-\hbar\omega\right). \]
El elemento de la matriz de perturbación \(H'_{cv}\) depende de la polarización de la luz y de los estados que participan en la interacción. La función delta impone la conservación de la energía; al construir un espectro práctico, se suman las contribuciones de las transiciones permitidas sobre los estados muestreados y se amplían. Los factores de ocupación de las portadoras determinan entonces si la absorción o la emisión estimulada predominan en una longitud de onda determinada. Un factor neto esquemático útil es \(f_c-f_v\), donde \(f_c\) y \(f_v\) son las ocupaciones de electrones de los estados de conducción y de valencia, respectivamente.
Dado que la capa de GaN presenta una amplia banda prohibida, sus espectros calculados se manifiestan en el rango ultravioleta cercano, principalmente alrededor de 330 a 370 nm en el ejemplo proporcionado. La longitud de onda y la energía del fotón están relacionadas mediante la ecuación \(E_{\gamma}=hc/\lambda\).
7.1 Emisión estimulada y ganancia neta.
El gráfico de emisión estimulada (??) contiene una serie de espectros para diferentes condiciones de la señal portadora. El cálculo de la ganancia neta (??) tiene en cuenta la competencia entre la emisión y la absorción. Una ganancia positiva del material indica una amplificación dentro de este modelo electrónico, pero por sí sola no demuestra que una cavidad láser completa haya alcanzado el umbral: también deben considerarse el factor de confinamiento y las pérdidas del resonador.
7.2 Absorción
Los resultados de la absorción se muestran en ??. Sus bordes y características espectrales están determinados por las separaciones entre los estados de valencia ocupados y los estados de conducción disponibles, ponderados por sus elementos de matriz óptica. En un pozo de tipo wurtzita, los espectros TE y TM correspondientes pueden diferir porque el contenido orbital de un estado determina la intensidad con la que se acopla a cada polarización.
Las características que se observan en estos espectros son el resultado del modelo configurado y del ensanchamiento. No deben identificarse automáticamente con picos excitónicos específicos: un cálculo de transición interbanda de un solo electrón, basado en el modelo k·p, no constituye, por sí solo, un tratamiento completo de los excitones ligados por la interacción de Coulomb entre electrones y huecos.
8. Emisión espontánea
A diferencia de la emisión estimulada, la emisión espontánea se produce sin que un fotón incidente inicie la transición. Depende de la disponibilidad de un electrón en el estado de conducción y de un estado vacío en la banda de valencia (un hueco). Por lo tanto, su factor de ocupación electrónica es proporcional a
\[ R_{\mathrm{sp}}(\omega)\ \propto\ \sum_{c,v,\mathbf{k}_{\parallel}} \left|M_{cv}(\mathbf{k}_{\parallel})\right|^2 f_c(1-f_v)\,L\!\left(E_c-E_v-\hbar\omega\right), \]
donde \(L\) indica el ensanchamiento espectral elegido y la proporcionalidad omite los prefactores fotónicos. El resultado se representa gráficamente en ??. La escala logarítmica permite comparar espectros que abarcan varios órdenes de magnitud. Observe las características ópticas cerca de la parte del espectro correspondiente a longitudes de onda más cortas y el cambio en la magnitud general en las diferentes condiciones representadas.
9. Pruébalo tú mismo.
💡 Dos cálculos adicionales que pueden ser útiles:
- Modifique el ancho de la capa de GaN. Intente reducir la capa central de 10 nm a 7 nm, dejando las barreras sin cambios. Vuelva a calcular los autoestados y compare los niveles de conducción más bajos y las energías de transición óptica. Un confinamiento más fuerte generalmente aumentará la separación entre los niveles confinados, aunque la respuesta detallada de la banda de valencia es multibanda.
- Compare las intensidades de TE y TM. Utilice las intensidades ópticas registradas por el terminal para identificar las transiciones que favorecen fuertemente una polarización en particular. Analice cómo esto se relaciona con el carácter del estado de valencia, en lugar de asumir que todas las transiciones contribuyen por igual al espectro.
✅ ¿Qué aspectos debe tener en cuenta?
En general, un pozo más estrecho tiende a aumentar las energías de confinamiento de los electrones y a modificar los niveles de los huecos, lo que a su vez altera las energías disponibles para las transiciones interbanda. Las relaciones TE/TM pueden variar considerablemente de un par de estados a otro, ya que examinan diferentes componentes orbitales de las funciones de onda de la estructura wurtzita. Para un estudio cuantitativo de los materiales, también es importante analizar los efectos de la composición de la aleación, la tensión y los campos de polarización.