Simulación de un diodo de unión PN de InP: deriva-difusión, rectificación y recombinación de portadores.

1. Introducción: una unión rectificadora de InP.

El fosfuro de indio (InP) es un semiconductor III-V con una banda prohibida directa que se utiliza como base para dispositivos electrónicos de alta velocidad y tecnologías fotónicas integradas. Un diodo de unión PN es una de sus estructuras de dispositivo más simples y más ilustrativas: una región de tipo p se encuentra con una región de tipo n, lo que crea un campo eléctrico interno y una barrera de potencial que controla el movimiento de electrones y huecos. El diodo conduce fuertemente cuando se aplica una polarización directa y suprime la inyección de portadores cuando se aplica una polarización inversa, lo que produce rectificación eléctrica.

El InP tiene una banda prohibida a temperatura ambiente de aproximadamente 1,34 eV y una movilidad electrónica relativamente alta. Estas propiedades hacen que el transporte en la unión de InP sea relevante para la electrónica de semiconductores compuestos y para comprender las interfaces de semiconductores que se utilizan en las plataformas de dispositivos basadas en InP. El ejemplo de este tutorial demuestra el funcionamiento eléctrico de una unión homojunción compacta p-InP/n-InP con una capa de contacto en el lado n⁺.

En lugar de ajustar la curva J–V a una ecuación ideal de diodo, OghmaNano resuelve las ecuaciones acopladas de Poisson y de continuidad de deriva-difusión. Los datos resultantes revelan tanto la corriente en los terminales como la distribución interna de las energías de las bandas de conducción y de valencia, las densidades de electrones y huecos, los niveles cuasi-Fermi y la recombinación. Esto nos permite explicar por qué el diodo rectifica, en lugar de simplemente observar que lo hace.

La ventana principal de la simulación, que se muestra en ??, visualiza el dispositivo de InP con capas verticales. Aunque la geometría se muestra en tres dimensiones, el problema eléctrico que se resuelve aquí es unidimensional y se analiza a lo largo de la dirección de crecimiento.

2. Crear una nueva simulación.

En la ventana principal de OghmaNano, seleccione «Nueva simulación». Se abrirá el menú de selección de simulaciones, que mostrará las categorías de ejemplos disponibles (??). Abra las demostraciones de III–V/InP y seleccione el ejemplo de diodo PN de InP (??). La ubicación exacta del ejemplo puede variar según la versión de OghmaNano.

Una vez que se abra el ejemplo, seleccione la pestaña «Estructura del dispositivo» y abra el editor de capas. Verifique los nombres de las capas, sus espesores y los modelos de materiales eléctricos asociados antes de ejecutar la simulación. Las siguientes secciones explican la contribución de cada región a la respuesta general del diodo.

3. Análisis de la estructura del dispositivo, los contactos y el dopaje.

3.1 Estructura de la capa semiconductora

El dispositivo eléctrico contiene tres capas de InP con un espesor total de semiconductor de 700 nm. Las capas de p-InP y n-InP forman la unión homojunción rectificadora, mientras que la región de n⁺-InP proporciona una región del lado del cátodo con una mayor concentración de dopantes. La asignación de las capas se puede ver en el editor (??).

La región de aire de 200 nm que se muestra en el editor de la estructura forma parte de la visualización y de la geometría óptica circundante, y no de la pila eléctrica del semiconductor. Por lo tanto, su grosor no se incluye en el grosor total del dispositivo eléctrico, que es de 700 nm.

3.2 Dopaje y electrostática integrada.

Abra la sección de parámetros eléctricos → Dopaje/iones para ver las concentraciones de donantes y aceptores (??). El cambio neto en el dopaje invierte su signo en la interfaz entre el material p-InP y el material n-InP, aproximadamente a 250 nm de la parte superior de la pila de semiconductores activos. El segundo cambio, entre el material n-InP y el material n⁺-InP, modifica la densidad de electrones en dirección al cátodo, pero no crea otra unión PN.

En el estado de equilibrio, los electrones tienden a difundirse desde la región de tipo n hacia la región de tipo p, mientras que las «huecos» se difunden en la dirección opuesta. Los átomos dopantes ionizados que quedan sin cubrir generan un campo que se opone a una mayor difusión. La región de agotamiento resultante y la barrera de potencial interna constituyen la base de la respuesta rectificadora.

Aquí, \(\phi\) representa el potencial electrostático, \(\varepsilon\) la permitividad dieléctrica y \(\rho\) la densidad de carga local. Los átomos dopantes ionizados y los portadores móviles determinan el campo eléctrico local, mientras que el voltaje aplicado modifica la barrera efectiva para la inyección de portadores.

3.3 Contactos eléctricos

El modelo utiliza una barrera óhmica con predominio de huecos en el ánodo y una barrera óhmica con predominio de electrones en el cátodo. Estas condiciones de contorno son adecuadas para estudiar la respuesta de la unión sin introducir una barrera de inyección adicional en los contactos. En los dispositivos fabricados, la metalización y el procesamiento de los contactos determinan qué tan cerca se acercan los contactos al comportamiento óhmico ideal.

4. Modelos de materiales y parámetros de transporte.

4.1 Modelos de materiales de Lua: origen de los parámetros de InP.

Una de las características importantes de este ejemplo de diodo de InP es que los parámetros de las propiedades eléctricas del material se definen mediante modelos de materiales Lua. Esto establece una conexión clara entre las propiedades físicas de un semiconductor y la solución de difusión-deriva. Las regiones p-InP, n-InP y n⁺-InP hacen referencia a la misma definición de material, Inorganic/InP/InP.lua. Su diferente comportamiento eléctrico se debe principalmente al perfil de dopaje, y no a la asignación de tres conjuntos distintos de constantes de material.

Para analizar el modelo, abra la sección Estructura del dispositivo → Parámetros eléctricos (el editor de materiales eléctricos) y seleccione una de las capas de InP. La pestaña Editor de código muestra el script de material Lua asociado a esa región; la pestaña Parámetros muestra las propiedades del material calculadas. Estas dos vistas deben utilizarse conjuntamente. El código define cómo se obtienen las propiedades, mientras que los valores calculados revelan qué es lo que realmente recibe el solucionador eléctrico.

Las imágenes de ?? y ?? muestran cómo se accede al editor de código. La configuración de la capa de la simulación es la que determina el modelo de material asignado: en este ejemplo, cada región activa está asociada al archivo Lua de InP. Los nombres de los archivos de imagen se heredan de la secuencia de captura de pantalla y no determinan qué semiconductor se está simulando.

En la simulación proporcionada, los interruptores db_overrides se configuran en false para los parámetros clave de transporte y de banda. Por lo tanto, OghmaNano está configurado para obtener los valores correspondientes de la base de datos de materiales/modelo Lua en lugar de utilizar los valores predeterminados que se guardan en otro lugar del archivo JSON. Esto es importante al examinar la configuración: los valores efectivos que se muestran en la pestaña de Parámetros son más informativos que los valores predeterminados almacenados en un campo de anulación inactivo.

4.2 Estructura de bandas, densidad de estados y movilidad.

Los primeros parámetros que se deben verificar en la pestaña de parámetros evaluados son la banda prohibida \(E_g\) y la afinidad electrónica \(\chi\). La banda prohibida establece la separación entre las bandas de conducción y de valencia, mientras que la afinidad electrónica indica la posición del borde de la banda de conducción en relación con una referencia de energía de vacío. Para el InP, la banda prohibida a temperatura ambiente es de aproximadamente \(1.34\,\mathrm{eV}\). Una comprobación inicial útil para verificar la coherencia es determinar si la separación entre \(E_c\) y \(E_v\) en el diagrama de bandas resultante coincide con la banda prohibida evaluada.

Las densidades de estados efectivas de las bandas de conducción y de valencia, \(N_c\) y \(N_v\), respectivamente, influyen en la relación entre la concentración de portadores y la posición del nivel cuasi-Fermi. Por lo tanto, contribuyen a las poblaciones de equilibrio de electrones y huecos, así como a la respuesta a la inyección bajo polarización directa. La permitividad relativa \(\varepsilon_r\) interviene en la ecuación de Poisson e influye en el potencial electrostático y en el ancho de la zona de agotamiento. Estas cantidades deben analizarse como un conjunto coherente, y no como parámetros numéricos independientes que se ajustan por separado.

Las movilidades del electrón y del hueco, \(\mu_n\) y \(\mu_p\), determinan la respuesta a la deriva de los portadores de carga ante un campo eléctrico; además, la difusión depende de los gradientes de densidad de portadores. En el modelo unidimensional, las densidades de corriente se expresan de la siguiente manera:

Aquí, \(n\) y \(p\) representan las densidades locales de electrones y huecos, \(E\) es el campo eléctrico y \(D_n\) y \(D_p\) son los coeficientes de difusión. Los dos componentes de corriente pueden variar considerablemente a lo largo del diodo, incluso cuando la corriente total en estado estacionario es aproximadamente independiente de la posición. Por este motivo, las secciones siguientes analizan \(J_n\) y \(J_p\) junto con los perfiles de densidad de portadores.

4.3 Cómo inspeccionar y comparar la configuración de los materiales

En el editor eléctrico, seleccione p-InP y examine su código Lua y los parámetros calculados. A continuación, cambie a n-InP y, por último, a n⁺-InP. Compare el modelo de material asignado, la banda prohibida, la afinidad electrónica, las movilidades, las densidades de estados efectivas y los coeficientes de recombinación. Para esta unión homojunción, se utiliza el mismo modelo de InP subyacente en las tres regiones; el principal cambio intencional entre las regiones es su concentración de aceptores o donantes.

Es importante distinguir dos conceptos: el modelo de material describe las propiedades del InP, mientras que el perfil de dopaje establece qué portadores de carga predominan localmente y dónde se encuentra la unión. Esta distinción resulta especialmente útil al pasar de una unión PN homopolar a una heterounión basada en InP, donde las diferencias entre los modelos de material pueden generar, por sí mismas, desplazamientos de banda. En este diodo, la flexión de las bandas de conducción y valencia que se observa en las imágenes refleja principalmente la electrostática a través de la transición p/n.

Para obtener un estudio reproducible, mantenga el modelo de material de Lua constante al modificar las concentraciones de dopaje. Por el contrario, si está investigando la sensibilidad del resultado a la movilidad o al coeficiente de recombinación, modifique ese parámetro del material de forma intencionada y compare la curva J–V resultante, las densidades de portadores y los perfiles de corriente interna. El editor de parámetros y las vistas del código Lua proporcionan el vínculo entre los supuestos físicos de entrada y la salida simulada.

4.4 Parámetros de recombinación

Las regiones de InP incluyen la recombinación de Shockley-Read-Hall (SRH) y un canal de recombinación de portadores libres, tal como se describe en el modelo de material utilizado. La recombinación de SRH describe la captura de electrones y huecos facilitada por defectos; la recombinación de portadores libres es especialmente relevante en el InP de banda prohibida directa, ya que las transiciones electrón-hueco pueden producirse de forma radiativa. En el presente ejercicio, ambos procesos son útiles como mecanismos de pérdida de portadores que modifican la corriente eléctrica.

El coeficiente \(B\), los tiempos de vida de las trampas \(\tau_n\) y \(\tau_p\), y las densidades auxiliares \(n_1\) y \(p_1\) determinan la magnitud de cada proceso. El aspecto clave que se analizará en este tutorial es cómo cambia la distribución espacial de la recombinación a medida que el sesgo impulsa un mayor número de electrones y huecos hacia la misma región.

4.5. Malla y solución numérica.

El cálculo eléctrico utiliza una malla unidimensional de 150 puntos a lo largo de 700 nm, lo que corresponde a una separación promedio de aproximadamente 4,7 nm. Las ecuaciones no lineales acopladas se resuelven mediante la iteración de Newton y el solucionador lineal disperso UMFPACK. Esta malla constituye un punto de partida útil para analizar la unión y las regiones cuasi-neutras.

La salida incluye información sobre la convergencia y perfiles de los dispositivos que dependen del sesgo. Una buena comprobación inicial consiste en analizar la suavidad de la curva J–V y, a continuación, confirmar que las corrientes de electrones y huecos cumplen con la ley general de conservación de la corriente en estado estacionario. El editor de malla se muestra en ??.

5. Ejecutar la simulación e interpretar la respuesta de la curva J–V.

Vuelva a la ventana principal y seleccione «Ejecutar simulación». OghmaNano aplica diferentes voltajes y registra la solución obtenida en cada punto. El archivo de ejemplo guardado utiliza un rango de voltaje de 0 a 3 V con incrementos de 25 mV. La ventana de resultados registra el progreso del solucionador durante la simulación, mientras que la pestaña «Salida» permite acceder a los datos resultantes.

Una vez finalizada la simulación, los resultados se muestran en la pestaña de resultados (??). Abra jv.csv para generar un gráfico de la densidad de corriente en función del voltaje aplicado. Las capturas de pantalla incluyen tanto una representación lineal de la curva J–V (??) como una representación logarítmica de la corriente (??).

5.1. Comprensión de las características.

Cuando se aplica un voltaje bajo, la barrera de la unión integrada restringe considerablemente la inyección de portadores. Al aumentar la polarización directa, se reduce esa barrera, por lo que los electrones y las huecos entran más fácilmente en el lado opuesto y la corriente total aumenta. En el régimen intermedio, la corriente suele ser aproximadamente exponencial con respecto al voltaje; con corrientes más elevadas, el transporte y las pérdidas resistivas pueden modificar la pendiente.

Esta expresión del diodo ideal proporciona una referencia útil para analizar un intervalo aproximadamente exponencial, pero no es impuesta por el modelo de deriva-difusión. El factor de idealidad \(n\) y el prefijo \(J_0\) deben determinarse únicamente en un intervalo de polarización adecuado; pueden variar en función de la recombinación, la resistencia en serie y el comportamiento del contacto.

La gráfica logarítmica ayuda a identificar el inicio de una corriente significativa y los cambios en la pendiente. La gráfica lineal facilita la interpretación de la magnitud y la forma de la curva de corriente elevada. Para obtener una imagen más clara del funcionamiento del diodo, examine ambas gráficas junto con los perfiles de la banda interna y la densidad de portadores.

6. Funcionamiento interno: cómo funciona el diodo.

El directorio «snapshots» contiene perfiles espaciales guardados en los puntos de polarización seleccionados durante el barrido J–V. Las capturas de pantalla proporcionadas muestran cuatro posiciones a lo largo de ese barrido. Estas imágenes muestran una progresión desde una inyección baja hacia una polarización directa más intensa. Las cuatro imágenes que componen cada grupo corresponden a las mismas posiciones seleccionadas en el control deslizante, lo que permite comparar directamente las diferentes magnitudes físicas.

6.1 Bandas de conducción y de valencia, y niveles cuasi-Fermi.

Abra el visor de imágenes y seleccione Ec.csv, Ev.csv, Fn.csv y Fp.csv. Estos representan el borde de la banda de conducción, el borde de la banda de valencia, el nivel cuasi-Fermi de los electrones y el nivel cuasi-Fermi de las huecas. Su separación espacial y sus pendientes revelan la barrera de potencial, la respuesta a la polarización aplicada y la desviación del equilibrio.

A baja inyección, los bordes de la banda se curvan a través de la unión como resultado del campo electrostático inherente. A medida que aumenta la polarización directa, la barrera efectiva se reduce y los niveles cuasi-Fermi de electrones y huecos se separan a medida que los portadores se desplazan a través del dispositivo. El punto clave es que el InP sigue siendo el mismo semiconductor en ambos lados de la unión: los perfiles cambiantes reflejan su electrostática y la inyección de portadores, que están definidas por el dopaje.

La comparación entre ?? y ?? ilustra cómo el sesgo terminal modifica la estructura de bandas internas. En particular, los niveles cuasi-Fermi proporcionan un indicador directo del estado de no equilibrio que se establece durante la conducción en estado estacionario.

6.2 Densidades de electrones y huecos.

A continuación, se muestran los gráficos Q_nfree.csv y Q_pfree.csv, que representan la distribución espacial de los electrones libres y las «huecos». A baja polarización, los electrones predominan en el lado de tipo n y los «huecos» en el lado de tipo p. Sus densidades disminuyen considerablemente alrededor de la unión, lo que refleja la región de agotamiento.

A medida que aumenta la polarización directa, se inyectan portadores minoritarios en los lados opuestos. Esto aumenta la superposición de las poblaciones de electrones y huecos y modifica la forma de ambos perfiles de densidad. En niveles de inyección más altos, las poblaciones de portadores libres a lo largo de la unión difieren considerablemente de sus valores correspondientes a una polarización baja.

Compare la variación de la densidad con la posición de la transición de dopaje y, a continuación, con las ubicaciones donde se produce una fuerte flexión de la banda. Esto relaciona los portadores suministrados por el contacto con la barrera electrostática generada por la unión PN.

6.3 Densidades de corriente de electrones y huecos.

Represente gráficamente Jn.csv y Jp.csv. Estos perfiles muestran la contribución de los electrones y las «huecos» a la corriente eléctrica que atraviesa el dispositivo. Los componentes individuales pueden variar según la posición, ya que la recombinación transfiere portadores entre las poblaciones de electrones y «huecos», pero la corriente total en estado estacionario debería permanecer aproximadamente constante cuando no se acumule ninguna carga neta.

Cuando la corriente es muy baja, la anulación entre los términos de deriva y difusión puede hacer que los perfiles de corriente individuales parezcan ruidosos o con una estructura muy marcada. Las capturas de pantalla obtenidas con una baja polarización muestran este comportamiento de forma más evidente que las capturas obtenidas con una corriente más alta. Es útil examinar los perfiles en una escala adecuada y comparar la corriente total con el valor terminal J–V antes de interpretar los detalles numéricos.

6.4 Perfil de recombinación

La recombinación determina en qué punto se eliminan los electrones y las «huecos» inyectados de sus respectivas poblaciones de portadores. Abra el visor de imágenes con R_nfree_to_pfree.csv y R_ss_srh.csv para comparar los canales de recombinación directa y los canales asistidos por defectos, respectivamente. Ambos procesos dependen de la disponibilidad de electrones y huecos, por lo que sus distribuciones espaciales cambian a medida que cambian las poblaciones de portadores inyectados.

Al comienzo de la conducción, la recombinación puede concentrarse cerca de la unión, donde se superponen las poblaciones de electrones y huecos. Con una inyección directa más intensa, la superposición se extiende a una parte mayor del semiconductor, y los perfiles de recombinación pueden volverse más amplios. Compare la posición y el ancho de los picos de recombinación con las densidades de portadores a la misma polarización.

Al comparar los dos canales, hay que tener en cuenta tanto el valor máximo como la integral espacial. Un pico local pronunciado no tiene por qué ser el factor dominante en la recombinación integrada si otro canal está activo en un volumen mucho mayor. Esto es especialmente importante al relacionar los perfiles de recombinación con la corriente eléctrica total.

6.5 Perfiles de la tasa neta de generación de electrones y huecos.

Las imágenes etiquetadas como R_n.csv y R_p.csv ofrecen información adicional sobre los valores que se introducen en las ecuaciones de continuidad de electrones y huecos. Estas imágenes pueden indicar dónde se concentra con mayor intensidad el equilibrio de los términos locales de generación, captura y recombinación. Interprete sus signos de acuerdo con la convención de representación gráfica utilizada por OghmaNano.

En conjunto, la estructura de bandas, las densidades de portadores, las corrientes y los perfiles de velocidad proporcionan una explicación completa de las características del diodo: la barrera interna controla la inyección, el voltaje aplicado modifica la electrostática interna y el transporte y la recombinación determinan la corriente resultante en estado estacionario.

7. Ejercicios e interpretación.

8. Interpretación del modelo y análisis complementario.

Este ejemplo utiliza un modelo de deriva y difusión de estado estacionario unidimensional de InP, una estructura semiconductora compacta p/n/n⁺ y terminales óhmicos idealizados. Es ideal para demostrar cómo surge la rectificación eléctrica de la barrera de la unión y cómo este comportamiento se relaciona con la distribución espacial de los portadores de carga. Estudios más detallados del dispositivo pueden incorporar las características de contacto medidas, parámetros de material más precisos, variación de la temperatura y detalles geométricos adicionales.

Las propiedades del material InP se obtienen del modelo Lua seleccionado. Para estudios cuantitativos, confirme la banda prohibida, la afinidad electrónica, la movilidad y los coeficientes de recombinación evaluados en el editor de materiales o en los datos de salida. Verifique la conservación de la corriente y los residuos del solucionador, y repita los puntos de polarización representativos con una malla más fina para establecer la estabilidad numérica.

En el ejemplo proporcionado, el rango estándar de barrido va de 0 a 3 V. Un experimento útil que se puede realizar a continuación consiste en incluir un rango modesto de polarización inversa, lo que permite comparar directamente la conducción directa y la corriente de fuga inversa. El rango de polarización inversa y las características de fuga que son relevantes desde el punto de vista físico dependen del dispositivo y de los mecanismos incluidos en la simulación.

Para obtener más información, varíe las concentraciones de dopaje, una a la vez. En general, aumentar la concentración nominal de dopaje reduce la región de agotamiento y modifica la electrostática en equilibrio; cambiar los coeficientes de recombinación influye en la corriente, que depende de la inyección. Estas son comparaciones controladas útiles porque los materiales y la geometría de la capa pueden permanecer constantes.

9. Próximos pasos.

El diodo PN de InP constituye un punto de partida útil para estructuras semiconductoras más avanzadas basadas en InP. Una extensión lógica sería introducir una región ligeramente dopada y comparar el comportamiento electrostático y de transporte de una unión PN y una unión PIN. Otra opción sería examinar heterostructures que involucren InGaAs/InP, donde las diferencias en la composición introducen desplazamientos de banda adicionales y un nuevo comportamiento en la inyección de portadores.