Cuanto de GaInNAs con bajo contenido de nitrógeno: Tutorial de 10 bandas sobre el modelo k·p para la estructura de zincblenda.
1. Introducción
En este tutorial, simularemos una heterounión cuántica de GaInNAs con una baja concentración de nitrógeno utilizando el modelo de diez bandas de OghmaNano para la estructura de tipo zincblenda. Abriremos un ejemplo existente, examinaremos la estructura de capas, ejecutaremos el cálculo del estado electrónico y analizaremos los espectros ópticos predichos. El ejemplo contiene una heterounión cuántica de GaInNAs de 7 nm dentro de una estructura más grande de GaAs/AlGaAs, en lugar de una heterounión de tres capas aislada.
¿Por qué es útil un modelo de diez bandas? La teoría k·p estándar de ocho bandas describe una banda de conducción y las bandas de valencia de agujeros pesados, agujeros ligeros y bandas de valencia separadas, incluyendo el espín. La introducción de una pequeña concentración de nitrógeno en una aleación III-V puede crear un nivel resonante relacionado con el nitrógeno que interactúa de manera especialmente intensa con la banda de conducción del material base. Una descripción de diez bandas añade dos componentes de espín a este estado relacionado con el nitrógeno. La interacción puede reducir significativamente la energía de transición fundamental y modificar la curvatura de la banda de conducción. Esto hace que los nitruros diluidos sean interesantes para dispositivos ópticos de longitudes de onda más largas.
2. Abrir el ejemplo del ecualizador de diez bandas.
Seleccione Archivo → Nueva simulación. En la ventana de nueva simulación, haga doble clic en los ejemplos de GaAs (??). La ventana siguiente contiene varios ejemplos de pozos cuánticos. Seleccione AlGaAs, aleación diluida de nitruro, 10 bandas k.p (??). Esto cargará la simulación preconfigurada de diez bandas, que incluye las definiciones de capa y los archivos de parámetros de material correspondientes.
3. Estructura del dispositivo y parámetros de los materiales.
La ventana principal muestra la estructura multicapa (??). A diferencia de las estructuras simples de GaAs/AlGaAs utilizadas en los tutoriales anteriores, este ejemplo incluye capas adicionales de revestimiento, capas de transición y capas de barrera alrededor de la capa de nitruro diluido. Para examinar la geometría, haga clic en «Editor de capas» en la barra de herramientas de la izquierda.
El editor de capas (??) muestra ocho entradas, que van desde p-GaAs hasta AlGaAs y la región activa, y luego hasta n-GaAs. En particular, la capa central de GaInNAs (pozo cuántico) tiene un grosor de 7,0 nm y su tipo de capa está configurado como pozo cuántico. El ejemplo visible también incluye capas de AlGaAs de 0,16 µm a cada lado de la región activa y una capa de barrera/SCL de GaAs de 0,02 µm. No es necesario reconstruir la estructura manualmente cuando se utiliza la plantilla proporcionada.
3.1. Revisar las definiciones de los materiales.
Haga clic en Parámetros eléctricos, seleccione la pestaña GaInNAs (QW) y abra el editor de código (??). La captura de pantalla muestra cómo el archivo de material contiene cantidades físicas, como la densidad de estados efectiva de la banda de valencia y la movilidad eléctrica. Por ejemplo, la movilidad de electrones a bajo campo que se observa es 0.05 m² V⁻¹ s⁻¹, equivalente a 500 cm² V⁻¹ s⁻¹. Estos son datos de entrada para el modelo del material, no valores inferidos de la dispersión del pozo cuántico que se muestra en el gráfico.
En un cálculo con bajo contenido de nitruro, es especialmente importante que la composición, las energías del borde de la banda y los parámetros de acoplamiento del nitrógeno sean mutuamente consistentes. Para esta primera ejecución, deje los valores predeterminados sin modificar; más adelante, al final de la página, analizaremos un experimento útil sobre la composición.
4. ¿Qué cambios se producen cuando introducimos nitrógeno?
Antes de ejecutar el modelo, es importante comprender por qué el hamiltoniano tiene diez bandas en lugar de ocho. Imaginemos que el límite de la banda de conducción del semiconductor huésped se desplaza hacia arriba en la escala de energía a medida que aumenta el vector de onda del electrón en el plano. Un estado de resonancia relacionado con el nitrógeno tiene su propia energía característica. Si los dos estados no interactuaran, sus energías podrían acercarse y cruzarse. Sin embargo, en un nitruro diluido, el nitrógeno se acopla a los estados de conducción del huésped, lo que modifica los niveles de energía resultantes.
Un ejemplo ilustrativo de este principio de la física es la matriz de dos estados que describe el fenómeno de evitación de cruce de bandas (BAC):
\[ H_{\mathrm{BAC}}(k)= \begin{pmatrix} E_C(k) & V_N\\ V_N & E_N \end{pmatrix}, \qquad V_N=\beta_N\sqrt{x_N}. \]
Aquí, \(E_C(k)\) representa la dispersión de la banda de conducción del material huésped, \(E_N\) es la energía del nivel relacionado con el nitrógeno, \(x_N\) es la fracción de nitrógeno y \(\beta_N\) controla la intensidad del acoplamiento. Al diagonalizar esta pequeña matriz, se obtienen dos ramas:
\[ E_{\pm}(k)=\frac{E_C(k)+E_N}{2} \pm\frac{1}{2}\sqrt{\big[E_C(k)-E_N\big]^2+4V_N^2}. \]
El término de la raíz cuadrada impide que las ramas acopladas se crucen cuando \(V_N\neq0\). En particular, la rama inferior, que se comporta como una banda de conducción, puede desplazarse hacia abajo en relación con la banda de conducción original, lo que reduce la banda prohibida fundamental. Dado que la interacción depende de la separación entre \(E_C(k)\) y \(E_N\), la dispersión también puede volverse fuertemente no parabólica.
Este ejemplo simplificado, que considera un sistema de dos por dos, explica la idea principal; no se trata de la solución completa del hamiltoniano del pozo cuántico que se presenta en este tutorial. El modelo de diez bandas k·p conserva la base más amplia de conducción/valencia de tipo zincblenda, incluyendo el espín, y añade los componentes relacionados con el nitrógeno. A continuación, se incluye el confinamiento cuántico a lo largo de la dirección de crecimiento. Para obtener el tratamiento teórico completo, consulte el capítulo sobre el modelo de diez bandas de zincblenda.
5. Ejecutar la simulación.
Vuelva a la pestaña «Archivo» y pulse «Ejecutar simulación». Seleccione la pestaña «Terminal» para seguir el cálculo (??). El programa ensambla un hamiltoniano disperso, obtiene los autoestados de conducción y valencia confinados y repite el cálculo en una secuencia de vectores de onda en el plano \(k_{\parallel}\).
En la captura de pantalla, el registro identifica el solucionador KP10, muestra los niveles de energía calculados e informa sobre la convergencia del bloque después de varias iteraciones. Los pares de energías muy próximas pueden incluir degeneraciones relacionadas con el espín; no deben interpretarse automáticamente como ramas distintas de anti-cruce del nitrógeno. Espere a que finalice la ejecución antes de examinar los archivos de salida.
6. Comprensión de la estructura de bandas calculada.
Los gráficos resultantes muestran la energía en función del vector de onda en el plano, y no la energía en función de la profundidad dentro del dispositivo. El potencial del pozo cuántico confina las cargas a lo largo de la dirección de crecimiento de la capa, lo que produce estados discretos; el movimiento paralelo a las capas le confiere a cada estado confinado una dispersión \(E_n(k_{\parallel})\). En un modelo parabólico simple, dicha dispersión es cuadrática. Un hamiltoniano de diez bandas acopladas permite que se desvíe de una parábola.
6.1 Estados de conducción y reducción de la separación.
El diagrama de conducción (??) muestra dos ramas calculadas que parten de diferentes niveles de energía. El nivel de conducción más bajo comienza aproximadamente en 0,89 eV en la referencia de energía utilizada para la simulación. No se debe confundir este valor con la banda prohibida del GaInNAs: se trata de una energía propia de esta estructura confinada en particular. Observe que las dos curvas representadas no tienen la misma curvatura. Este es uno de los aspectos que se deben analizar al comparar un cálculo de ocho bandas con el modelo de nitruro diluido.
Las ecuaciones simplificadas de la concentración de área activa (BAC) que se presentan anteriormente ayudan a explicar por qué la dispersión de la conductividad puede variar en un material que contiene nitrógeno. Sin embargo, por sí solas, no demuestran la proporción de nitrógeno en ninguna de las dos curvas coloreadas. Para ello, es necesario analizar los correspondientes autovectores o modificar \(x_N\) y observar cómo evolucionan los estados.
6.2 Estados de valencia
El diagrama de valencia (??) contiene varias ramas con formas claramente no parabólicas. Las características de los huecos pesados y los huecos ligeros pueden combinarse cuando el vector de onda en el plano es finito; las energías de confinamiento, la interacción espín-órbita y el acoplamiento entre los estados base influyen en las curvas. Un cambio ascendente o una modificación de la curvatura no indican necesariamente un fallo numérico. Por el contrario, una curva trazada de forma suave no es, por sí sola, prueba de que los parámetros del material sean correctos.
Es preferible considerar los colores individuales como ramas de autoestados representadas gráficamente, en lugar de asignarles de forma permanente las etiquetas HH o LH en todo el rango del vector de onda sin verificar sus características.
6.3 Analizando ambas bandas en conjunto.
Al activar la dispersión combinada (??), se colocan los estados electrónicos y de valencia en la misma escala de energía. La transición interbanda fundamental está relacionada con su separación, y no con la energía absoluta de ninguno de los estados representados:
\[ \hbar\omega_{nm}(k_{\parallel}) =E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}). \]
Esta expresión describe una transición óptica directa, aproximadamente vertical, en el vector de onda del plano: el momento del fotón es pequeño en relación con el momento del cristal electrónico. Sin embargo, el espectro observado no suele ser una línea nítida. Contribuyen muchas transiciones permitidas, con intensidades y poblaciones que varían de un estado a otro.
7. De los estados electrónicos a los espectros ópticos.
Con los estados calculados disponibles, OghmaNano puede construir espectros ópticos. La cuestión física importante no es simplemente si un electrón y una laguna tienen la diferencia de energía adecuada, sino también cuán fuertemente se acoplan ambos estados al campo electromagnético y si sus estados iniciales y finales están ocupados de manera apropiada.
7.1 ¿Por qué algunas transiciones tienen un impacto mayor que otras?
La regla de oro de Fermi proporciona un punto de partida útil para comprender las tasas de transición. Para una transición de un estado electrónico inicial \(i\) a un estado final \(f\), inducida por una perturbación óptica débil, su forma esquemática es:
\[ W_{i\rightarrow f}=\frac{2\pi}{\hbar} \left|\langle f|H_{\mathrm{int}}|i\rangle\right|^2 \delta(E_f-E_i-\hbar\omega). \]
El elemento de la matriz \(\langle f|H_{\mathrm{int}}|i\rangle\) describe la interacción de los estados con la luz. Depende, entre otras cosas, de su carácter orbital, las funciones de envolvente y la polarización de la luz. La función delta impone la conservación de la energía. En un espectro numérico, la transición idealizada y abrupta se representa mediante una función adecuada de ensanchamiento de la línea. Tras realizar la suma sobre las subbandas permitidas y los vectores de onda en el plano, el resultado es un espectro y no una única energía de transición.
En este caso, las ocupaciones de los portadores de carga permiten distinguir entre la absorción y la emisión. Si \(f_c\) representa la ocupación de un estado de conducción y \(f_v\) la ocupación de un estado de electrón de valencia, los factores característicos son:
\[ \begin{aligned} \text{absorption:}\quad &(1-f_c)f_v,\\ \text{stimulated emission:}\quad &f_c(1-f_v),\\ \text{net stimulated contribution:}\quad &f_c-f_v. \end{aligned} \]
La última expresión representa la diferencia entre los procesos ópticos ascendentes y descendentes. Para que haya una contribución positiva a la amplificación, es necesario que exista una inversión de población suficiente. Estos factores explican la interpretación de los siguientes gráficos; para obtener más detalles, consulte la sección teoría de la transición óptica.
7.2 Emisión y absorción estimuladas
Los dos gráficos (?? y ??) ilustran los procesos ópticos que compiten entre sí para una misma estructura de bandas de un pozo cuántico. Las múltiples curvas representan la secuencia de condiciones utilizadas en el ejemplo; sus magnitudes relativas varían porque la ocupación de los portadores determina qué transiciones son posibles.
Los espectros se extienden hasta el infrarrojo cercano, aproximadamente desde 1,1 hasta 1,5 µm. A modo de referencia, la energía de los fotones y la longitud de onda en el vacío están relacionadas mediante la siguiente ecuación:
\[ E_{\mathrm{photon}}=\frac{hc}{\lambda} \simeq \frac{1240\ \mathrm{eV\,nm}}{\lambda\,[\mathrm{nm}]}. \]
Por lo tanto, 1300 nm corresponde aproximadamente a 0,954 eV. Esto permite realizar una comprobación útil para verificar la coherencia entre la separación de energía en el diagrama de bandas combinado y la escala de los espectros ópticos, teniendo en cuenta el confinamiento, el rango de transiciones que contribuyen y la referencia de energía elegida.
7.3 Ganancia neta: se produce cuando la emisión supera la absorción.
La magnitud de interés para un medio amplificador es el equilibrio entre la emisión estimulada y la absorción. Utilizando la convención de signos empleada en el gráfico, se puede representar esquemáticamente de la siguiente manera:
\[ g_{\mathrm{net}}(\lambda) =g_{\mathrm{stim}}(\lambda)-\alpha(\lambda). \]
La gráfica de la ganancia neta calculada (??) incluye tanto curvas negativas como positivas. Cuando la absorción es el factor dominante, \(g_{\mathrm{net}}\) es negativo. A medida que cambian las poblaciones de electrones y huecos, la emisión estimulada puede llegar a ser lo suficientemente intensa como para que en algunas partes del espectro se desarrolle una ganancia positiva en el material. La región de transición representada en la gráfica presenta una estructura notable, en lugar de ser un único pico perfectamente suave, debido a que contribuyen múltiples transiciones en estados confinados.
Importante: obtener resultados positivos en cuanto a la ganancia de material no es, por sí solo, un indicador de que un láser completo funcionará por encima del umbral. También deben tenerse en cuenta las pérdidas en la cavidad, el confinamiento óptico, la reflectividad de los espejos y la inyección de portadores.
8. Localización de los resultados de la simulación.
bands y gain.Una vez que finalice el cálculo, abra la pestaña de resultados (??). La captura de pantalla muestra los archivos de resultados principales, que incluyen: Ec.csv (ramas de conducción), Ev.csv (ramas de valencia), Ec_Ev.csv (bandas combinadas), gain.csv, absorption.csv, net_gain.csv y spontaneous.csv. Los resultados adicionales se organizan en los directorios bands y gain. Haga doble clic en un archivo CSV para generar un gráfico o utilice los datos numéricos guardados para realizar análisis adicionales.
Aunque este tutorial ilustra la emisión estimulada, la absorción y la ganancia neta, el ejemplo también genera una salida de emisión espontánea. La emisión espontánea es una recombinación radiativa que no requiere la presencia de un fotón incidente, y su factor de ocupación característico es \(f_c(1-f_v)\). Sus coeficientes espectrales absolutos difieren de los de la emisión estimulada.
9. Pruébalo tú mismo: investiga la interacción del nitrógeno.
Experimento sugerido: Guarde una copia independiente de este ejemplo antes de modificarlo. A continuación, analice la concentración de nitrógeno en la definición del material GaInNAs, modificando un parámetro a la vez y manteniendo fija la geometría del pozo de 7 nm.
- Registre la energía propia más baja del electrón en el punto \(k_{\parallel}=0\).
- Modifique la proporción de nitrógeno dentro de un rango físicamente adecuado para la parametrización seleccionada.
- Vuelva a calcular la dispersión de la conducción y compare la forma de la rama inferior.
- Compare la posición del borde de absorción y las características principales de la ganancia.
La ilustración del modelo BAC predice que aumentar el acoplamiento del nitrógeno puede incrementar la separación de las ramas acopladas y modificar la energía de conducción más baja. Sin embargo, en un cálculo completo de la aleación, la composición también puede afectar a otros parámetros del material. Por lo tanto, interprete cualquier cambio resultante utilizando las funciones reales del material que estén habilitadas en su simulación, en lugar de asumir que la fórmula simple de dos estados es el modelo completo de la aleación.
✅ ¿Qué aspectos hay que tener en cuenta?
Investigue si el borde de la banda de conducción se desplaza y si su curvatura cambia al variar la concentración de nitrógeno. Compare los espectros ópticos en condiciones de ocupación idénticas, de modo que un cambio en las energías de transición no se confunda con un cambio no relacionado en la población de portadores. Si puede analizar la composición de los autovectores, será especialmente útil observar cómo cambian las contribuciones del material huésped y del nitrógeno a los estados.
Para un segundo experimento, mantenga los parámetros del material constantes y cambie el ancho de la capa cuántica de 7 nm a 6 nm. En general, aumentar el confinamiento cuántico tiende a elevar la energía del electrón confinado de nivel más bajo y modifica las energías de transición interbanda permitidas. La magnitud de este desplazamiento se determina mediante el cálculo de las diez bandas acopladas, y no simplemente asignando una única masa efectiva constante a cada estado.