Electric-field effects on defect migration energetics in GaN

Resumen Técnico: Efectos del campo eléctrico en la energética de la migración de defectos en GaN

Planteamiento del problema
La fiabilidad de los dispositivos basados en GaN y AlGaN, particularmente en entornos de alta potencia y radiación, se ve comprometida por la migración y reconfiguración de defectos puntuales. Si bien la formación de defectos bajo irradiación ha sido ampliamente estudiada, el comportamiento de los defectos existentes bajo la influencia combinada de campos eléctricos operativos y temperatura sigue siendo menos comprendido. En los transistores de alta movilidad electrónica (HEMT) y dispositivos de potencia, los campos operativos a menudo se acercan al límite de ruptura (~3 MV/cm), lo que puede impulsar el transporte de defectos que conduce a inestabilidades en el voltaje de umbral y deriva en el rendimiento. Los enfoques computacionales existentes enfrentan un compromiso: la Teoría del Funcional de la Densidad (DFT) ofrece una alta precisión pero está restringida a superceldas pequeñas, mientras que la Dinámica Molecular (MD) clásica con potenciales interatómicos estándar carece de la capacidad para modelar la transferencia de carga dinámica y el acoplamiento electrostático requeridos para simular la migración de defectos bajo campos eléctricos externos.

Metodología
Para abordar estas limitaciones, los autores desarrollaron un campo de fuerza reactivo ReaxFF específico para GaN. La metodología consistió en tres etapas principales:

1. Desarrollo del Campo de Fuerza: Se generó un conjunto de datos de entrenamiento exhaustivo utilizando cálculos DFT (VASP) dentro del marco del funcional PBE. Este conjunto de datos incluyó datos de la ecuación de estado (EOS) para varias fases de Ga y GaN (wurtzita, sal de roca, zincoblenda, etc.), perfiles moleculares de N2 y energías totales de superceldas que contienen defectos puntuales neutros (vacantes, antisitios, pares de Frenkel, intersticiales y defectos de Schottky).
2. Validación: Los parámetros de ReaxFF se optimizaron frente a este conjunto de datos DFT. El campo de fuerza resultante fue validado comparando los parámetros de red, energías de cohesión, módulo de compresibilidad y energías de formación/migración de defectos frente a resultados previos de DFT y experimentales.
3. Simulaciones de Campo Eléctrico: Utilizando el modelo ReaxFF validado dentro de LAMMPS, los autores realizaron cálculos de Banda Elástica de Nudged (CI-NEB) para determinar las barreras de migración para defectos puntuales intrínsecos ($V_{Ga}$, $V_N$, $Ga_{i1,oct}$ y $N_{i2,split}$) bajo campos eléctricos externos que variaron de 0 a 0.04 V/Å (0–40 MV/cm). Las simulaciones tuvieron en cuenta la equilibración de carga dinámica (QEq) y las distorsiones locales de la red.

Contribuciones Clave y Resultados

• Validación del Campo de Fuerza: El modelo ReaxFF desarrollado reproduce con éxito las propiedades estructurales de DFT (por ejemplo, constantes de red dentro de un 0.3% y módulo de compresibilidad dentro de un 1.5%) y la energética de los defectos. Aunque existen algunas desviaciones en las energías absolas de formación de defectos (hasta un 26% para ciertas configuraciones), el modelo captura las tendencias relativas y las barreras de migración con alta fidelidad respecto a la literatura.
• Respuesta de Campo No Lineal: Contrario a la simple aproximación lineal de dipolo-campo ($\Delta E_B = -\mu_{eff} \cdot E$), el estudio revela que los campos eléctricos modifican las barreras de migración de forma anisotrópica y no lineal. La respuesta depende fuertmente del tipo de defecto específico y su mecanismo de migración.
  • Vacantes: La migración de las vacantes de nitrógeno ($V_N$) exhibe una respuesta casi lineal y simétrica al campo, consistente con un modelo de dipolo rígido debido a la mínima deformación de la red y variación de carga. Las vacantes de galio ($V_{Ga}$) muestran una asimetría leve; la barrera hacia adelante aumenta mientras que la barrera inversa disminuye, pero con pendientes desiguales, atribuido al mayor tamaño atómico del Ga que induce deformación local y redistribución de carga.
  • Intersticiales: Los defectos intersticiales ($Ga_{i1,oct}$ y $N_{i2,split}$) muestran una marcada no linealidad y asimetría. Para $Ga_{i1,oct}$, el campo es lo suficientemente fuerte como para desplazar la configuración del punto de silla limitante de la velocidad, causando un cruce en qué pico de energía domina la barrera. Estos defectos involucran el movimiento cooperativo de múltiples átomos, lo que conduce a una compleja redistribución de carga parcial y respuestas de deformación que se desvían significamente de las predicciones lineales.
• Mecanismo de Desviación: La desviación del comportamiento lineal surge de la redistribución de carga parcial inducida por el campo y la distorsión de la red local a lo largo de la coordenada de reacción. El momento dipolar efectivo no es fijo, sino que evoluciona continuamente a medida que el defecto migra, un fenómeno capturado nativamente por el campo de fuerza reactivo pero omitido por los modelos de dipolo fijo.

Significancia y Reivindicaciones
El artículo sostiene que los campos eléctricos no solo sesgan el movimiento de los defectos de forma lineal, sino que pueden modificar de manera compleja el paisaje de migración de defectos a través del acoplamiento carga-red.

• Impacto Cinético: El estudio destaca que, si bien las reducciones de barrera inducidas por el campo para vacantes de alta barrera (>2 eV) pueden no alterar significamente las tasas de difusión en escalas de tiempo de dispositivos, los efectos son cinéticamente significativos para los defectos intersticiales con barreras cercanas a 1 eV. Una reducción de ~0.05 eV puede aumentar las tasas de migración en un orden de magnitud a temperatura ambiente, potencialmente acelerando la redistribución de defectos en regiones de alto campo de los dispositivos en operación.
• Fiabilidad de los Dispositivos: Estos hallazgos sugieren que en las regiones de alto campo (por ejemplo, bordes de compuerta en HEMTs), los defectos intersticiales pueden volverse sustancialmente más móviles, contribuyendo a la generación de trampas y a la deriva del voltaje de umbral. La asimetría direccional de la respuesta implica un desplazamiento neto de los defectos hacia regiones específicas del dispositivo, lo que potencialmente impulsa problemas de fiabilidad a largo plazo.
• Alcance Metodológico: Los autores señalan que sus resultados corresponden a defectos neutros debido a la naturaleza clásica del esquema de equilibración de carga QEq. Reconocen que los defectos altamente cargados en condiciones reales de operación (por ejemplo, $V_{Ga}^{3-}$) pueden exhibir respuestas de campo aún más fuertes, lo que sugiere que los hallazgos actuales representan un límite inferior de los efectos de campo. No obstante, este trabajo establece que el acoplamiento campo-defecto en GaN es anisotrópico y específico de cada defecto, lo que requiere simulaciones reactivas para una predicción precisa bajo condiciones operativas realistas.