Tutorial sobre pozos cuánticos parabólicos de GaAs/AlGaAs
1. Introducción
En este tutorial, calcularemos las propiedades electrónicas y ópticas de un pozo cuántico simple de GaAs/AlGaAs utilizando el modelo de banda parabólica en OghmaNano. Se coloca una fina capa de GaAs entre dos barreras de AlGaAs con una banda prohibida más amplia. El cambio en la energía del borde de la banda confina los electrones y las huecos en la dirección de crecimiento, mientras que estos siguen siendo libres de moverse en paralelo a las capas. Esto produce energías de confinamiento discretas que se desarrollan en subbandas electrónicas bidimensionales.
Comenzaremos con una simulación ya preparada, analizaremos la estructura de las capas, ejecutaremos el modelo y observaremos las subbandas de conducción y valencia calculadas. A continuación, examinaremos cómo esos estados contribuyen a la absorción, la emisión estimulada, la ganancia óptica neta y la emisión espontánea. Si no está familiarizado con la física, página sobre la teoría de la banda parabólica ofrece una introducción más detallada.
2. Abrir el ejemplo del pozo cuántico.
Abra OghmaNano y seleccione «Nueva simulación» en la pestaña «Archivo». En el menú de selección de dispositivos, haga doble clic en «Ejemplos de GaAs» (??). Esto abrirá la selección de ejemplos de semiconductores y pozos cuánticos. Haga doble clic en «GaAs/AlGaAs QW parabólico» (??). El ejemplo preconfigurado se abrirá en la ventana principal de la simulación.
3. Comprensión de la estructura del dispositivo.
La ventana principal muestra la estructura de tres capas (??). La capa de GaAs es la capa activa; las dos capas de AlGaAs que la rodean forman las barreras. Haga clic en «Editor de capas» en el lado izquierdo de la ventana principal para consultar las dimensiones. El ejemplo proporcionado utiliza la siguiente estructura:
| Capa | Grosor | Propósito |
|---|---|---|
| AlGaAs | 40 nanómetros | Primera barrera |
| GaAs (heteroestructura de pozo cuántico) | 10 nanómetros | Pozo cuántico |
| AlGaAs | 40 nanómetros | Segunda barrera. |
Observe que a la capa central se le asigna explícitamente el tipo de capa de pozo cuántico, en lugar de la capa activa. Las entradas y la selección de materiales ópticos se pueden ver en ??. Las barreras de mayor banda prohibida crean los potenciales de confinamiento; los valores exactos de la diferencia de banda dependen de los parámetros de los materiales utilizados en el ejemplo.
¿Por qué es importante modificar el ancho del pozo?
El confinamiento a lo largo de la dirección del crecimiento permite únicamente la existencia de funciones de onda específicas para los electrones y las huecas. En la idealización más simple de un pozo cuántico unidimensional infinito de anchura \(L\), las energías de confinamiento de los electrones, medidas desde el fondo del pozo, son:
\[ E_n=\frac{n^2\pi^2\hbar^2}{2m^{*}L^2},\qquad n=1,2,3,\ldots \]
Aquí, \(m^{*}\) representa la masa efectiva de las portadoras y \(\hbar\) es la constante de Planck reducida. En un ejemplo real, las barreras de AlGaAs son finitas, por lo que OghmaNano calcula los estados confinados utilizando el perfil real del material, en lugar de asumir esta expresión de pozo infinito. No obstante, la ecuación proporciona una predicción útil: reducir el ancho del pozo de GaAs generalmente aumenta las energías de confinamiento de los electrones. Los huecos tienen su propia masa efectiva y las correspondientes energías de confinamiento.
Puede examinar o modificar las funciones de los materiales subyacentes a través del botón de parámetros eléctricos. En la captura de pantalla, la definición del material AlGaAs incluye un cálculo de la banda prohibida que depende de la composición y la temperatura (??). No es necesario modificar el código del material para completar este tutorial introductorio.
4. Ejecutar la simulación.
Vuelva a la pestaña «Archivo» y pulse «Ejecutar simulación» (el botón triangular azul). La pestaña «Terminal» mostrará los mensajes del solucionador mientras se ejecuta. En el ejemplo proporcionado, la terminal muestra las energías de los autoestados y repite los mensajes block converged after ... iterations (??). Espere a que termine la ejecución antes de abrir los gráficos.
El programa de resolución calcula los estados cuantificados en las diferentes capas y, a continuación, analiza cómo varían en función del vector de onda en el plano, \(k_{\parallel}\). Cuando se configura el cálculo óptico, también genera espectros para las diferentes condiciones de portadora en el ejemplo.
5. Subbandas electrónicas: interpretación del diagrama de bandas.
Abra el gráfico de dispersión de bandas calculado (??). Se pueden observar dos grupos de curvas claramente diferenciados: las subbandas de conducción, correspondientes a energías positivas, y las subbandas de valencia, que se encuentran cerca y por debajo de la referencia de la banda de valencia. La separación vertical refleja la banda prohibida del semiconductor, junto con las energías confinadas de los electrones y las huecas. Las diferentes energías iniciales dentro de cada grupo se deben a los estados cuantificados del pozo.
Aunque el movimiento a través del pozo está limitado, los portadores aún pueden desplazarse a lo largo de las capas. Por eso, una única energía confinada se representa como una parábola cuando se grafica en función del vector de onda paralelo:
\[ E_{e,n}(k_{\parallel})=E_{e,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_e^{*}}. \]
Aquí \(E_{e,n}(0)\) se representa la energía en el punto más bajo de la subbanda de electrones \(n\). Para las subbandas correspondientes de electrones de valencia, la forma parabólica simple se curva hacia abajo:
\[ E_{v,m}(k_{\parallel})=E_{v,m}(0) -\frac{\hbar^2 k_{\parallel}^2}{2m_h^{*}}. \]
Esta es la aproximación fundamental del tutorial: los estados tienen energías cuantizadas en la dirección de confinamiento, pero conservan una dispersión de masa efectiva en el plano. Los modelos multibanda k·p más avanzados permiten que las bandas se mezclen y dejen de ser parabólicas; la página teoría de la estructura de zincblenda con 8 bandas explica por qué esto puede ser importante.
6. De las subbandas a los espectros ópticos.
Los electrones pueden experimentar transiciones ópticas entre las subbandas de valencia y conducción. En condiciones normales de incidencia óptica, el momento del fotón es pequeño en comparación con el vector de onda electrónico, por lo que un modelo simplificado considera las transiciones que conservan el vector de onda en el plano. En el centro de la dispersión, la energía aproximada del fotón para una transición entre el nivel de electrón confinado \(n\) y el nivel de hueco \(m\) es:
\[ \hbar\omega\simeq E_g+E_{e,n}^{\mathrm{conf}}+E_{h,m}^{\mathrm{conf}}. \]
Las dos energías de confinamiento se miden a partir de los bordes de las bandas de energía correspondientes. Las transiciones posibles, la ocupación de sus estados iniciales y finales, y cualquier ensanchamiento espectral determinan la forma de la respuesta óptica calculada. Puede explorar estas ideas con más detalle en página sobre la teoría de las transiciones ópticas.
6.1 Emisión estimulada y ganancia neta.
El gráfico de emisión estimulada (??) muestra la contribución calculada a la emisión en función de la longitud de onda para la serie de condiciones de portadora utilizadas en el ejemplo. Los distintos pasos y picos espectrales corresponden a las energías en las que diferentes transiciones de subbanda contribuyen. La salida neta de ganancia separada (??) incluye la contribución de la absorción estimulada, que compite con la emisión, y por lo tanto puede ser positiva o negativa.
De forma esquemática, si consideramos un estado de conducción con una probabilidad de ocupación de \(f_c\) y un estado de electrón de valencia con una probabilidad de ocupación de \(f_v\), la contribución neta inducida contiene el factor de ocupación.
\[ \underbrace{f_c(1-f_v)}_{\text{stimulated emission}} -\underbrace{f_v(1-f_c)}_{\text{absorption}} =f_c-f_v. \]
Para lograr una contribución neta positiva a la estimulación, es necesario que se produzca una inversión suficiente de los estados que participan en la transición. El espectro total que se obtiene también depende de la intensidad de la transición óptica y del número de estados disponibles.
6.2 Absorción y emisión espontánea.
El gráfico de absorción (??) registra la contribución de la absorción de las transiciones permitidas. La absorción requiere un estado de valencia inicial ocupado y un estado de conducción final disponible. El gráfico de emisión espontánea (??) representa la recombinación radiativa de estados de conducción ocupados en estados de valencia vacíos, sin necesidad de un fotón incidente para estimular la transición. En consecuencia, el factor de ocupación de la emisión espontánea contiene \(f_c(1-f_v)\).
Observe que las características espectrales más marcadas se producen en regiones de longitud de onda similares en todos estos gráficos. Esto es de esperar: la misma estructura subbanda de electrones y huecos subyacente proporciona las posibles energías de transición óptica, mientras que la ocupación de los portadores diferencia la absorción de los distintos tipos de emisión. El gráfico de emisión espontánea utiliza una escala vertical logarítmica, por lo que tenga cuidado al comparar las alturas de las curvas.
7. Localización de los archivos calculados.
Una vez finalizada la simulación, abra la pestaña de resultados en la ventana principal (??). En el ejemplo, se generan archivos CSV separados para los distintos cálculos espectrales. Los nombres de archivo visibles incluyen entradas que siguen los patrones gain_00.csv con los sufijos de absorción, ganancia neta, emisión espontánea y emisión estimulada, seguidos de archivos con nombres similares para las condiciones calculadas posteriores. Haga doble clic en un resultado para examinarlo en la ventana de gráficos. Dado que los resultados se guardan como archivos CSV, también se pueden analizar o volver a graficar fuera de OghmaNano.
Los archivos numerados le permiten seguir cómo evoluciona el espectro calculado a lo largo de la serie de condiciones de la portadora, en lugar de limitarse a analizar un único cálculo óptico.
8. Pruébalo tú mismo: modifica el ancho del pozo.
Experimento sugerido: vuelva al editor de capas y cambie el grosor de la capa de GaAs (pozo cuántico) de 10 nm a 8 nm, dejando las capas de barrera de AlGaAs sin modificar. Vuelva a ejecutar la simulación y compare las energías de subbanda y los espectros ópticos resultantes con el ejemplo original.
- ¿Se desplaza el nivel de confinamiento electrónico más bajo hacia arriba en relación con el borde de la banda de conducción del GaAs?
- ¿Cómo cambian las energías de transición de los electrones y los huecos?
- ¿Se produce un desplazamiento del espectro óptico hacia longitudes de onda más cortas o más largas?
✅ Tendencia física prevista
En general, un pozo más estrecho aumenta las energías de confinamiento. Por lo tanto, la transición interbanda de menor energía tiende a desplazarse hacia una energía fotónica más alta, lo que corresponde a una longitud de onda más corta, aunque el desplazamiento exacto depende de las alturas finitas de la barrera y de las masas efectivas del electrón y el hueco. Como se indica en \(\lambda=hc/E_{\mathrm{photon}}\), una energía de transición más alta implica una longitud de onda de emisión o absorción más pequeña.
Una vez que se domine este ejemplo, la comparación lógica siguiente es con la versión de 8 bandas del modelo k·p de estructura de zincblenda aplicada a un pozo cuántico de GaAs/AlGaAs. Este modelo conserva el confinamiento cuántico, pero va más allá de las bandas parabólicas independientes para describir las interacciones y la mezcla entre los estados de conducción y valencia.