Tutorial de GaAs (Parte A): Ejecutar una simulación 3D de deriva–difusión con un defecto vertical

1. Introducción

En este tutorial ejecutará una simulación de deriva–difusión 3D de un diodo pn de GaAs que contiene un defecto vertical intencionadamente evidente (una derivación de arriba abajo a través de la pila). Este es un problema deliberadamente “injusto” para modelos 1D: el defecto fuerza flujo lateral de corriente, amontonamiento de corriente y recombinación localizada que no pueden existir en una simulación puramente a través del espesor. Por qué esto importa: en dispositivos reales, el rendimiento suele estar limitado no por la pila ideal de capas sino por dónde fluye realmente la corriente. Defectos locales, pinholes, fugas en los bordes, contactos finitos y materiales no uniformes pueden dominar la JV medida. Un modelo 3D de deriva–difusión le permite vincular rasgos de la curva JV directamente con regiones específicas del dispositivo, respondiendo preguntas que 1D no puede.

Cómo usar esto en su propio trabajo: trate el defecto aquí como una plantilla para cualquier no idealidad localizada—por ejemplo un filamento de derivación, un pinhole a través de una capa bloqueante, un límite de grano, una región dañada por sonda o un electrodo de área finita. El flujo de trabajo es el mismo: construya una geometría 3D mínima, mantenga la malla lateral tan pequeña como sea posible, ejecute el solucionador y luego use curvas JV resueltas por contacto e instantáneas (jn.csv, \(\phi\), recombinación) para conectar el comportamiento eléctrico con la ubicación física.

El objetivo es doble: (1) darle un punto de entrada práctico al modelado eléctrico genuinamente 3D en OghmaNano, y (2) enseñar la decisión de modelado: cuándo el 3D es esencial (defectos, contactos finitos, variación lateral) y cuándo colapsa a 1D. En la Parte A ejecutará el dispositivo defectuoso, en la Parte B eliminará el defecto y mostrará que 3D/2D/1D coinciden una vez que se restaura la simetría, y en la Parte C habilitará la iluminación y verá dónde aparece la sensibilidad numérica a baja corriente.

2. Primeros pasos

Desde la ventana principal de OghmaNano, haga clic en New simulation para abrir la biblioteca de dispositivos (véase ??). Haga doble clic en GaAs demos. Esto abre una carpeta de ejemplos de arseniuro de galio (véase ??). Finalmente, haga doble clic en GaAs - 3D defect y guarde la simulación en una carpeta en la que tenga permisos de escritura.

Después de abrir el ejemplo, la ventana principal muestra una vista 3D del dispositivo (véase ??). La estructura es un diodo pn de GaAs que contiene dos bloques visibles que representan una región de defecto vertical intencionada (por ejemplo, un filamento de impureza o una trayectoria de pérdida tipo cortocircuito). El defecto va de arriba abajo a través del dispositivo, haciéndolo intrínsecamente tridimensional y capaz de producir amontonamiento de corriente y recombinación espacialmente localizada que estarían ausentes en un modelo puramente 1D.

Para inspeccionar los parámetros del material y del defecto, vaya a la pestaña Device structure y haga clic en Electrical parameters para abrir el editor de parámetros (véase ??). Encontrará conjuntos de parámetros separados para:

• GaAs p+ y GaAs n+ (las capas del diodo), y
• defect0 y defect1 (las regiones del defecto).

El ejemplo utiliza movilidades similares a las de GaAs y densidades de portadores físicamente razonables. La recombinación incluye un canal bimolecular simple libre-a-libre (a menudo denominado en contextos de dispositivos orgánicos como n³–p³) con un coeficiente representativo de 1 × 10⁻¹⁵ m³s⁻¹.

También observará que están habilitadas trampas de Shockley–Read–Hall (SRH) en equilibrio. En muchas simulaciones de dispositivos orgánicos se requieren trampas SRH dinámicas porque la ocupación de trampas almacena carga significativa y realimenta la electrostática. Para GaAs, sin embargo, las densidades de trampas suelen ser lo suficientemente bajas como para que las trampas puedan tratarse puramente como un mecanismo de pérdida por recombinación. Usar SRH en equilibrio captura esta física sin introducir variables de estado adicionales para la carga almacenada en trampas, manteniendo el modelo estable y eficiente.

3. Explorar la simulación: dopado y malla

La cinta Electrical proporciona acceso a los editores eléctricos principales utilizados para definir y diagnosticar la simulación (véase ??). En esta sección nos centraremos en dos herramientas: Doping/Ions, que define la unión pn, y Electrical mesh, que controla la resolución espacial de la solución de deriva–difusión.

Haga clic en Doping/Ions para abrir el editor de dopado (véase ??). Este ejemplo es un diodo pn deliberadamente simple, con regiones de tipo p y tipo n claramente separadas. Mantener aquí simple el perfil de dopado facilita atribuir cualquier comportamiento inusual posterior a la geometría (el defecto) y no a la electrostática.

A continuación, haga clic en Electrical mesh para abrir el editor de malla (véase ??). Verá que x, y y z están todas habilitadas, lo que significa que la simulación es genuinamente tridimensional. En este sistema de coordenadas:

• y es la dirección principal de transporte (a través del espesor del diodo),
• x y z son direcciones laterales utilizadas para resolver el defecto y el flujo lateral de corriente.

La malla en la dirección y está dividida deliberadamente entre las capas tipo n y tipo p para que cada región mantenga una resolución apropiada—esto se vuelve importante si más adelante cambia espesores de capa o introduce asimetría. En cambio, las mallas laterales (x y z) utilizan relativamente pocos puntos (por ejemplo, ~10 en x y ~5 en z). Esto es intencionado: en 3D, el coste computacional crece rápidamente. Si aumenta la densidad de malla en cada dimensión en un factor \(k\), el número total de incógnitas escala aproximadamente como \(k^3\). En términos prácticos, las simulaciones 3D se vuelven cúbicamente más lentas a medida que refina la malla, por lo que debería mantener la resolución lateral tan baja como sea posible mientras siga resolviendo la física que le interesa.

4. Ejecutar la simulación 3D y ver los resultados

Para ejecutar la simulación, haga clic en Run simulation (el triángulo azul). Debido a que este es un problema de deriva–difusión tridimensional genuino, tardará más que una ejecución 1D. En un portátil razonablemente moderno, el ejemplo debería completarse en aproximadamente 30 segundos. Si el tiempo de ejecución es significativamente mayor, compruebe primero que la simulación está guardada en una unidad local; las carpetas montadas en red o sincronizadas con la nube pueden ralentizar la ejecución sustancialmente debido al intenso I/O de archivos. Una vez finalizada la ejecución, abra la pestaña Output (véase ??). Este directorio contiene los archivos de salida estándar generados por una solución de deriva–difusión, incluidas instantáneas de variables internas y datos corriente–voltaje. En este punto es importante distinguir entre salidas agregadas y salidas resueltas por contacto.

Debido a que este dispositivo se simula en 3D con contactos de área finita, generalmente debería ignorar el archivo global jv.csv y utilizar en su lugar las curvas JV resueltas por contacto. Los archivos jv_contact0.csv y jv_contact1.csv informan de la corriente asociada a los contactos superior e inferior respectivamente. En simulaciones multidimensionales, estas curvas son la forma físicamente más significativa de interpretar el comportamiento eléctrico y de diagnosticar pérdidas específicas de contacto. Abra jv_contact0.csv para ver la curva JV del contacto superior (véase ??). Debido a que el dispositivo contiene un defecto vertical, no se espera que esta curva se parezca a una característica ideal de diodo pn 1D. En la siguiente sección (Parte B), deshabilitará el defecto y verá cómo cambia la respuesta JV—y cómo la misma estructura vuelve a colapsar a 2D y 1D una vez que la asimetría intrínsecamente tridimensional se elimina.

5. Explorar variables internas con el visor de instantáneas

Para inspeccionar lo que el solucionador está haciendo internamente, abra el directorio snapshots/ en la carpeta de salida. Esto da acceso a variables del simulador resueltas espacialmente en función del paso de simulación (normalmente la polarización aplicada). Al hacer doble clic en el directorio se abre el visor de instantáneas mostrado en ?? y ??.

En la ventana de instantáneas, haga clic en el botón azul + para añadir una gráfica y seleccione jn.csv en el menú desplegable File to plot. El archivo jn.csv corresponde a la densidad de corriente vertical de electrones. Usando el deslizador principal, puede explorar cómo evoluciona la distribución de corriente con el voltaje aplicado, mientras que los deslizadores y y z le permiten tomar cortes a través del volumen 3D.

Este ejemplo es un diodo en oscuridad, por lo que todo el flujo de corriente está impulsado eléctricamente y no es fotocorriente. La observación clave es que el defecto produce una trayectoria de conducción altamente localizada, que puede visualizarse directamente en 3D y vincularse sin ambigüedad al comportamiento JV no ideal. El visor de instantáneas puede utilizarse para explorar muchas otras magnitudes internas, incluidas Ec.csv y Ev.csv (bordes de banda), niveles cuasi-Fermi, tasas de recombinación y densidades de portadores. Explorar estos campos junto con la curva JV es una de las formas más eficaces de construir intuición sobre cómo se comportan los solucionadores de deriva–difusión en dispositivos multidimensionales.

Tenga en cuenta que, aunque el potencial electrostático \(\phi\) evoluciona con la polarización aplicada, la distribución de dopado (Nad.csv) está fijada por el editor de dopado y no cambia con el voltaje. Aun así, visualizar el dopado en 3D es útil—particularmente para tutoriales posteriores que impliquen variaciones laterales de dopado o defectos espacialmente localizados.