Simulación de un diodo PIN de GaN: modelos de deriva-difusión, curvas J-V y transporte en materiales de banda prohibida amplia.
1. Introducción
Un diodo PIN de nitruro de galio (GaN) es un dispositivo semiconductor de banda prohibida amplia que consta de una región de tipo p, una región con una concentración de impurezas relativamente baja (nominalmente intrínseca) y una región de tipo n. A diferencia de una unión PN convencional, un diodo PIN introduce deliberadamente una región de deriva considerable entre los terminales dopados. Bajo polarización inversa, esta región puede soportar un campo eléctrico distribuido; bajo polarización directa, los electrones y los huecos inyectados permiten que la estructura conduzca la corriente. Estas propiedades hacen que la arquitectura PIN sea relevante para los principios subyacentes de la rectificación de potencia, el bloqueo de la tensión inversa y la ingeniería de dispositivos de banda prohibida amplia.
El GaN tiene una banda prohibida a temperatura ambiente de aproximadamente 3,4 eV y puede soportar campos eléctricos más altos en comparación con el silicio. Sin embargo, los diodos de potencia de GaN prácticos requieren un control cuidadoso del grosor de la capa de deriva y de la concentración de impurezas de fondo, la resistencia de contacto y los defectos cristalinos, así como de la terminación de los bordes. En este ejemplo, la región intrínseca de GaN de 800 nm está diseñada para que sea compacta: esto facilita el estudio de la electrostática central y la inyección de portadores, pero por sí sola no representa un rectificador de potencia de alta tensión con las características especificadas.
Este tutorial utiliza las ecuaciones de Poisson y de deriva-difusión acopladas de OghmaNano para resolver el problema de un diodo PIN de GaN unidimensional. En lugar de tratar la corriente como una única magnitud de circuito, el modelo determina el perfil de la banda interna, los niveles cuasi-Fermi, las concentraciones de electrones y huecos y las tasas de recombinación a medida que se modifica el sesgo. El objetivo es relacionar una curva J–V rectificadora intuitiva con los procesos físicos que tienen lugar en el interior del semiconductor.
2. Crear una nueva simulación.
Abra OghmaNano y seleccione «Nueva simulación». Elija las demostraciones de semiconductores III-V y abra el ejemplo del diodo PIN de GaN. Las categorías de simulación disponibles y la ventana de selección de demostraciones se muestran en las figuras ?? y ??. Si la demostración ya está instalada, al abrirla debería mostrarse la estructura de capas en la ventana principal del dispositivo.
La vista principal de la simulación (??) muestra una representación tridimensional del dispositivo. Aunque la visualización es en 3D, el ejemplo numérico que se presenta aquí solo resuelve el transporte eléctrico a lo largo de la dirección de crecimiento de la capa. Abra el editor de capas para ver el grosor real y la asignación de materiales de cada región.
3. Estructura en capas y dopaje.
El prototipo finalizado utiliza la siguiente estructura semiconductora activa eléctricamente. La capa intrínseca de 800 nm es la configuración utilizada para la simulación que se muestra aquí; las versiones anteriores del ejemplo pueden contener una capa intrínseca de 400 nm y deben actualizarse antes de intentar reproducir estos resultados.
| Capa | Grosor | Concentración nominal de impurezas. | Propósito |
|---|---|---|---|
| p-GaN | 200 nm | Concentración de átomos aceptores: 5 × 1017 cm−3 | inyección en el lado p y región del ánodo |
| GaN intrínseco | 800 nm | No se utilizó ningún tipo de sustancia prohibida de forma intencionada. | Región de deriva y región de agotamiento con una baja concentración de impurezas. |
| nitruro de galio | 200 nm | Donantes: 5 × 1017 cm−3 | región de unión de tipo n |
| GaN dopado con n⁺ | 100 nanómetros | Donantes: 1 × 1018 cm−3 | Región de contacto posterior. |
La capa activa de semiconductores tiene un grosor de 1,3 µm. En el editor de capas, puede aparecer una región de aire separada de 200 nm para la visualización del dispositivo y la geometría óptica; no forma parte de la pila PIN activa desde el punto de vista eléctrico. Las cuatro regiones de semiconductores utilizan GaN, por lo que el propósito de las interfaces en este caso es modificar la concentración de impurezas, y no la composición de la aleación de semiconductores.
Abra el editor de contactos (véase la figura que aparece a continuación) para examinar el ánodo y el cátodo, y luego abra el editor de dopaje (véase la figura que aparece a continuación). La región p dopada con aceptores y las regiones n/n⁺ dopadas con donantes deben aparecer en lados opuestos de una amplia región central no dopada. El centro de la estructura no es una heterounión de materiales separada: su principal característica distintiva es un bajo nivel de dopaje intencional.
En la aproximación idealizada de agotamiento, una región con una baja carga neta puede soportar un campo eléctrico que se extiende a lo largo de gran parte de su espesor. El voltaje soportado por una región está relacionado con la integral espacial de su campo eléctrico:
\[V_{\mathrm{region}}=-\int_{y_0}^{y_1} E(y)\,\mathrm{d}y.\]Por lo tanto, una región más amplia con una ligera concentración de impurezas puede distribuir un potencial inverso aplicado a lo largo de una distancia mayor. Esta es una de las razones principales por las que se estudian las estructuras PIN para el bloqueo de voltaje. Esto no significa que el voltaje de ruptura se pueda obtener simplemente multiplicando el grosor por un campo crítico supuesto: los dispositivos reales tienen campos no uniformes, carga de fondo, contactos y concentración de campo en los bordes.
4. Modelos eléctricos, contactos y configuración numérica.
4.1 Parámetros del material GaN
Abra el editor de parámetros eléctricos para examinar la configuración de transporte y recombinación utilizada para cada región de GaN. La simulación está configurada para cargar las propiedades del GaN desde la base de datos de materiales OghmaNano, con las opciones de anulación de la base de datos desactivadas. Esto es preferible a asignar valores arbitrarios de movilidad efectiva, banda prohibida o afinidad electrónica en el archivo JSON del dispositivo. Verifique los datos de materiales cargados si necesita una comparación cuantitativamente reproducible con los resultados experimentales.
En el caso del GaN, las consideraciones principales incluyen el ancho de la banda prohibida, la movilidad de los electrones y las huecos, las estadísticas de los portadores libres, la recombinación radiativa y la recombinación asistida por trampas. La concentración nominal de aceptores que se muestra en el editor de dopaje no debe identificarse automáticamente con la densidad de huecos libres: los aceptores de magnesio en el GaN de tipo p no están completamente ionizados a temperatura ambiente. Se debe considerar en qué medida el modelo eléctrico gestiona la ionización de los dopantes antes de asignar una resistencia de contacto física o de obtener una densidad de portadores medida.
4.2 Contactos y condiciones de contorno
La demostración utiliza contactos óhmicos ideales para el ánodo del lado p y el cátodo del lado n. Esto simplifica la interpretación al evitar un segundo elemento rectificador en la interfaz metal/semiconductor. En los dispositivos de GaN fabricados, especialmente en los contactos de p-GaN, lograr un comportamiento óhmico de baja resistencia puede ser un desafío e implica la pila de metales, la preparación de la superficie y el dopaje.
4.3 Red eléctrica.
El editor de la malla eléctrica (que se muestra a continuación) indica la discretización a lo largo del espesor del semiconductor. Al modificar la región intrínseca de 400 nm a 800 nm, asegúrese de que la malla cubra la pila activa actualizada de 1,3 µm. Una elección inicial de aproximadamente 250 a 260 puntos proporciona un espaciamiento del orden de 5 nm; el software puede volver a generar la malla automáticamente, según la configuración. Para un estudio más cuantitativo, repita la simulación con una malla más fina y compruebe que la curva J–V y los perfiles internos no cambien de forma significativa.
4.4 Transporte y recombinación
OghmaNano combina la ecuación de Poisson electrostática con las ecuaciones de continuidad para electrones y huecos. En una dimensión, las ecuaciones en estado estacionario pueden expresarse esquemáticamente de la siguiente manera:
\[\frac{\mathrm{d}}{\mathrm{d}y}\!\left(\varepsilon\frac{\mathrm{d}\phi}{\mathrm{d}y}\right)=-\rho,\qquad \frac{\mathrm{d}J_n}{\mathrm{d}y}=q(R-G),\qquad \frac{\mathrm{d}J_p}{\mathrm{d}y}=-q(R-G),\]donde \(\phi\) es el potencial electrostático, \(\rho\) la densidad local de carga espacial, \(J_n\) y \(J_p\) las densidades de corriente de electrones y huecos, \(G\) la tasa de generación y \(R\) la tasa neta de recombinación. Con estas convenciones de signos, la densidad de corriente total \(J_n+J_p\) es constante en un dispositivo unidimensional que ha alcanzado un estado estacionario y convergente.
El modelo de recombinación puede incluir contribuciones de recombinación directa entre bandas (de electrones libres a electrones libres) y contribuciones del modelo de Shockley-Read-Hall (SRH). Su importancia relativa depende de la calidad del material, el nivel de inyección y los parámetros del modelo. La banda prohibida directa del GaN hace que la recombinación radiativa sea físicamente posible, pero su inclusión en un modelo no implica que un dispositivo específico representado en el modelo sea un LED eficiente.
5. Ejecución de la simulación J–V en condiciones de baja iluminación.
Establezca el rango de barrido de voltaje de -2,0 V a +5,0 V, con incrementos de 0,025 V. El voltaje negativo polariza inversamente la estructura PIN, y el voltaje positivo la polariza directamente, permitiendo la conducción. La configuración de polarización deseada es:
"Vstart": -2.0,
"Vstop": 5.0,
"Vstep": 0.025Desde la ventana principal, seleccione «Ejecutar simulación». Mientras se realiza el cálculo, utilice la pestaña «Terminal» para supervisar el progreso y la convergencia. Una vez finalizado, abra la pestaña «Resultados» y seleccione «jv.csv» para generar un gráfico de la densidad de corriente en función del voltaje aplicado. A continuación, se muestra el visor de resultados, la terminal del solucionador y la curva J–V resultante.
Cuando se aplica una polarización inversa, la corriente simulada es pequeña en comparación con su valor cuando se aplica una polarización directa. A medida que aumenta la polarización positiva, la inyección de portadores aumenta rápidamente y el diodo se activa. Con voltajes directos más altos, la curva se vuelve más lineal en un eje de densidad de corriente lineal, lo que indica que la corriente se ve cada vez más afectada por el transporte resistivo y otras no idealidades en la estructura simulada. Una pendiente que se aplana en un eje de corriente logarítmico no debe interpretarse automáticamente como una saturación brusca de la corriente: un aumento de corriente aproximadamente lineal también se vuelve progresivamente más plano en una gráfica logarítmica.
Para una comparación sencilla con la teoría de los diodos, la característica de polarización directa baja a moderada puede relacionarse con:
\[J\approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right],\]donde \(J_0\) es una corriente de saturación efectiva y \(n\) un factor de idealidad efectivo. Un diodo PIN puede desviarse considerablemente de un régimen ideal de Shockley único, ya que la recombinación en la región de agotamiento, la inyección en la capa intrínseca, la inyección a niveles altos y la resistencia en serie pueden influir en la corriente. Por lo tanto, la expresión simplificada debe considerarse más bien como una aproximación diagnóstica que como la ecuación fundamental que resuelve OghmaNano.
La capacidad de bloqueo en polarización inversa no es equivalente a la ruptura por avalancha observada. Una prueba que reduce la tensión hasta -2 V solo evalúa el comportamiento con una pequeña polarización inversa. La tensión a la que se produce la ruptura de un dispositivo de potencia de GaN fabricado depende de la física en campos elevados, los defectos del material, las vías de fuga, la geometría y el tratamiento de los bordes, y no se puede deducir únicamente de esta demostración.
6. Análisis de los resultados de la simulación interna.
Tras la ejecución completa del barrido J–V, el visor de imágenes de OghmaNano permite representar gráficamente los resultados con resolución espacial en los puntos de polarización seleccionados. Las capturas de pantalla que se muestran a continuación ilustran cuatro posiciones en el barrido almacenado, desde la primera hasta la última instantánea seleccionada. Sus nombres (-0 a -3) identifican las posiciones de las capturas de pantalla, no los voltajes. Consulte el valor real del voltaje seleccionado en el visor de imágenes al realizar comparaciones cuantitativas.
6.1 Bordes de la banda y niveles cuasi-Fermi.
Abra el visor de imágenes y represente gráficamente Ec.csv, Ev.csv, Fn.csv y Fp.csv en el mismo gráfico. Estos corresponden a los bordes de las bandas de conducción y de valencia, así como a los niveles cuasi-Fermi de electrones y huecos. Las figuras ?? a ?? muestran su evolución a lo largo del barrido de polarización.
En condiciones de polarización inversa, se espera que gran parte de la variación del potencial se produzca en la región intrínseca ligeramente dopada, lo que da como resultado una pendiente de banda relativamente amplia. La polarización directa reduce la barrera de inyección efectiva y crea una separación de los niveles cuasi-Fermi de electrones y huecos que depende de la polarización. Los gradientes espaciales de estos niveles son importantes porque determinan el flujo de corriente; la curvatura de la banda y la separación de los niveles cuasi-Fermi deben considerarse en conjunto, en lugar de interpretarse de forma independiente.
Como comprobación física básica, la separación entre \(E_c\) y \(E_v\) en las capas de GaN debe ser coherente con la banda prohibida de GaN que figura en la base de datos de materiales, en lugar de ser un valor predeterminado residual de otro semiconductor. Dado que se trata de un dispositivo PIN de un solo material, no se esperan cambios bruscos en la banda inducidos por la composición en las transiciones p/i/n.
6.2 Densidades de electrones y huecos.
Trace las curvas de Q_nfree.csv y Q_pfree.csv (o los resultados de la densidad de portadores de electrones y huecos que muestra la versión instalada). A bajas densidades de inyección, la mayoría de los portadores se concentran en sus respectivas regiones dopadas. La región central no dopada contiene densidades de portadores de equilibrio mucho más bajas. En condiciones de polarización directa intensa, los portadores inyectados ocupan una mayor proporción de la región intrínseca y la recombinación puede aumentar considerablemente.
El cambio brusco en la densidad nominal de impurezas no implica un perfil idéntico y brusco de portadores libres. La densidad de portadores se determina de forma autoconsistente a partir de las energías de la banda, los niveles de Fermi, la densidad de estados y la neutralidad de la carga. Esto es especialmente relevante en el lado del p-GaN si el modelo del material utiliza un tratamiento simplificado de la ionización incompleta de los átomos aceptores.
6.3 Recombinación en la región intrínseca.
Observe los gráficos de recombinación espacial en los mismos puntos de polarización. Las figuras de recombinación de este ejemplo muestran cómo evoluciona la región de superposición y recombinación significativas de electrones y huecos a medida que aumenta la inyección directa. Interprete la altura y el ancho de cada perfil trazado en conjunto: un pico local grande no necesariamente representa la mayor contribución a la tasa de recombinación total del dispositivo.
Para un proceso de transición directa entre bandas, la tasa neta de emisión radiativa suele seguir la siguiente forma:
\[R_{\mathrm{rad}}=B\left(np-n_i^2\right),\]donde \(B\) es el coeficiente de recombinación del material. La recombinación SRH depende de los parámetros de los defectos y puede presentar una distribución espacial diferente. Para comparar los mecanismos, represéntelos por separado si ambos están disponibles y habilitados en el dispositivo cargado. No dé por sentado que un gráfico de recombinación total sin etiquetar identifica un único mecanismo microscópico.
6.4. Tasas de recombinación/continuidad de electrones y huecos.
Los gráficos adicionales que muestran la velocidad de los electrones y los huecos son útiles para analizar dónde cambian más rápidamente las poblaciones de portadores de carga, ya sea por recombinación o por inyección. Las cuatro capturas de pantalla que se muestran a continuación indican cómo cambian estos perfiles calculados en función del sesgo. Si la velocidad mostrada tiene signo o incluye tanto contribuciones de generación como de recombinación, utilice la leyenda del gráfico y el encabezado del archivo de datos para determinar la interpretación física correcta antes de comparar las magnitudes.
6.5. La situación actual de conservación como modelo de referencia.
Para cualquier punto de polarización en estado estacionario, represente gráficamente la densidad de corriente de electrones \(J_n\) y la densidad de corriente de huecos \(J_p\) si se dispone de esos datos. Cada componente puede variar según la posición a medida que los portadores se recombinan, pero el total \(J_n+J_p\) debe mantenerse constante en el espacio en un dispositivo 1D convergente, sin términos adicionales de acumulación de carga. Esto resulta especialmente útil para verificar si una pequeña ondulación aparente en una señal de baja corriente podría confundirse con un efecto físico.
7. Ejercicios recomendados.
8. Supuestos y limitaciones físicas.
Este ejemplo está diseñado como una simulación educativa a gran escala del transporte de electrones y la difusión en materiales de nitruro de galio (GaN). Ofrece una primera visión, con base física, de la electrostática de los dispositivos PIN y de la inyección de portadores, pero no se trata de un modelo comercial calibrado. En particular, los resultados deben interpretarse teniendo en cuenta las siguientes consideraciones.
- Contactos idealizados: las barreras de contacto metal/GaN y la resistencia de contacto práctica del contacto p-GaN no se validan explícitamente.
- Dopaje nominal: la activación del aceptor de tipo p y la carga residual de fondo en el GaN, que nominalmente es intrínseco, pueden alterar la distribución física de la carga.
- Comportamiento en campos magnéticos intensos: la variación de la curva J–V no determina el fenómeno de ruptura por avalancha, el voltaje de bloqueo limitado por la corriente de fuga ni la fiabilidad en condiciones de polarización inversa elevada.
- Transporto unidimensional: el modelo no tiene en cuenta la concentración lateral de la corriente, las fugas en los bordes, las placas de campo y la geometría de terminación del dispositivo.
- Defectos y recombinación: para obtener predicciones cuantitativas realistas sobre las fugas y la recombinación, es necesario calibrar los parámetros de los defectos volumétricos y de la interfaz.
- Efectos térmicos y dinámicos: el autocalentamiento, las transiciones durante el cambio de estado y los efectos relacionados con la acumulación de carga de los portadores minoritarios deben evaluarse por separado cuando sea pertinente.
Estas limitaciones no disminuyen el valor del ejemplo como punto de partida. La disposición de las capas y la curva J–V ilustran lo que ocurre cuando se introduce una región de deriva con baja concentración de impurezas en un diodo semiconductor de banda prohibida amplia. También proporcionan una base de referencia reproducible para probar modelos más avanzados de transporte y recombinación.
9. Próximos pasos.
Tras analizar el diodo PIN de GaN, compare su distribución de campo con la de una unión PN de GaN convencional y, a continuación, investigue cómo el grosor de la capa intrínseca y el dopaje residual influyen en la rectificación. Un trabajo más avanzado podría incorporar estados de trampa calibrados, contactos no ideales y mecanismos de transporte en campos elevados para examinar los límites del modelo simple de deriva y difusión.