A stochastic simulation of the dislocation-mediated etching of porous GaN distributed Bragg reflectors

Título: Simulación estocástica del grabado mediado por dislocaciones de reflectores Bragg distribuidos (DBR) de GaN poroso

1. Planteamiento del Problema

El grabado electroquímico selectivo de conductividad de estructuras de nitruro de galio (GaN) permite fabricar reflectores Bragg distribuidos (DBR) porosos, esenciales para mejorar la extracción de luz en LEDs y láseres. El proceso tradicional requiere capas alternas de GaN dopado (n+) y no intencionalmente dopado (NID). Sin embargo, el desafío radica en cómo el agente químico (electrolito) accede a las capas dopadas enterradas bajo una capa superficial NID que no conduce corriente.

Históricamente, se han utilizado dos modelos para explicar este mecanismo:

• Modelo "Kebab": Asume que cada dislocación de torsión (TD) actúa como un canal vertical continuo que perfora todas las capas, permitiendo que el grabado ocurra en cada capa dopada alineada verticalmente.
• Modelo "Cascada" (Nuevo): Evidencia experimental reciente (mediante tomografía FIB y STEM) muestra que el camino de grabado no siempre es vertical. El electrolito puede viajar verticalmente por una TD, propagarse lateralmente a través de una capa dopada y activar una diferente TD para continuar hacia abajo. Esto crea estructuras complejas de "cascada" donde las TDs pueden activarse y desactivarse en diferentes profundidades.

La falta de un modelo físico completo que explique la transición entre estos comportamientos y su dependencia del voltaje aplicado motiva el desarrollo de una simulación estocástica.

2. Metodología

Los autores desarrollaron una simulación estocástica en Python (optimizada con la librería Numba) para modelar el proceso de porosificación.

• Modelo Geométrico: La estructura DBR se representa como una matriz de píxeles (5 nm x 5 nm). Incluye capas alternas de GaN dopado (azul) y NID (verde), con dislocaciones de torsión (TDs) modeladas como líneas verticales que atraviesan la estructura.
• Mecanismo de Grabado:
  • Las capas NID tienen una probabilidad de grabado ($P_{NID}$) de cero.
  • Las capas dopadas y las TDs tienen probabilidades de grabado finitas ($P_{doped}$ y $P_{disloc}$, respectivamente).
  • El algoritmo busca iterativamente píxeles adyacentes al electrolito y decide si se graban basándose en estas probabilidades.
  • Se implementó una restricción para evitar la fusión de poros adyacentes, imitando dominios de porosidad no interconectados observados experimentalmente.
• Correlación Experimental:
  • Se fabricaron muestras DBR (Muestra A) con capas uniformes y muestras con espesores variables (Muestras B y C) mediante epitaxia por fase de vapor de organometálicos (MOVPE).
  • Se realizaron experimentos de cronoamperometría a diferentes voltajes (6 V, 7 V, 8 V, 12 V) para medir la corriente en función del tiempo.
  • Se compararon los resultados de la simulación (corrientes simuladas y morfologías) con los datos experimentales y con tomografías FIB previas.

3. Contribuciones Clave

• Desarrollo de un Modelo Estocástico: Creación de una herramienta computacional capaz de reproducir tanto la morfología de los poros como las curvas de corriente experimental sin requerir una comprensión física completa y determinista de cada paso atómico.
• Validación del Modelo "Cascada": La simulación demuestra que, al variar las probabilidades de grabado, se pueden generar tanto estructuras tipo "kebab" como tipo "cascada", validando la hipótesis de que las TDs pueden activarse y desactivarse dinámicamente.
• Correlación Voltaje-Probabilidad: Establecen una relación cualitativa donde el aumento del voltaje aplicado en el experimento corresponde a un aumento en la probabilidad de grabado de las capas dopadas ($P_{doped}$) en la simulación.
• Aplicabilidad General: Demuestran que el modelo no solo funciona para DBRs estándar, sino que puede adaptarse a estructuras con espesores variables y geometrías de grabado desde paredes laterales (trenches), prediciendo cambios en la alineación de los poros.

4. Resultados Principales

• Reproducción de Curvas de Corriente: La simulación reproduce exitosamente las oscilaciones en las curvas de cronoamperometría.
  • A altos voltajes (alta $P_{doped}$), se observan picos de corriente claros y oscilaciones sostenidas, correspondientes a un grabado más uniforme capa por capa (tendencia al modelo "kebab").
  • A bajos voltajes (baja $P_{doped}$), las oscilaciones se atenúan y la corriente decae, reflejando un grabado menos ordenado y más dependiente de la propagación lateral (tendencia al modelo "cascada").
• Morfología de Poros:
  • La simulación predice correctamente que a voltajes bajos predominan las estructuras de cascada (TDs que se apagan y encienden).
  • A voltajes altos, aumenta la proporción de TDs que actúan como canales continuos ("K5" o "kebab" completo), coincidiendo con datos de tomografía FIB.
  • Se observa que la simulación no captura completamente la formación de poros cónicos en las capas NID cerca de las TDs a altos voltajes, ya que el modelo asume $P_{NID} = 0$, sugiriendo que mecanismos físicos adicionales (como dopantes no intencionales o estados en la banda prohibida) facilitan el grabado en NID.
• Efecto de Espesores Variables: En muestras con capas dopadas de espesor variable (Muestras B y C), la simulación predice correctamente las tendencias en la amplitud de los picos de corriente, aunque subestima la disminución de la eficiencia de grabado en las capas más profundas observada experimentalmente.

5. Significado e Impacto

• Herramienta de Diseño: La simulación ofrece una metodología rápida para predecir el comportamiento de grabado en diferentes estructuras de GaN, ayudando a optimizar parámetros de fabricación para dispositivos optoelectrónicos (LEDs, láseres VCSEL).
• Comprensión del Mecanismo: Aunque la física exacta del grabado mediado por TDs sigue siendo un misterio, el modelo sugiere que el voltaje afecta diferencialmente la probabilidad de grabado en las capas dopadas frente a las TDs en regiones NID. Esto apunta a que los mecanismos de conducción (tunelamiento Zener, dopantes no intencionales) juegan roles distintos en cada región.
• Relevancia en Geometrías Alternativas: El estudio sugiere que incluso en procesos que utilizan "trenches" (ranuras) litográficas para acceder a capas dopadas, el grabado mediado por TDs puede ocurrir a cierta distancia de la ranura, afectando la alineación de los poros y las propiedades ópticas (birrefringencia) del material resultante.

En conclusión, el trabajo proporciona un puente crucial entre la observación microestructural compleja y los datos macroscópicos de corriente, ofreciendo una simulación robusta que puede guiar la fabricación de futuros dispositivos basados en GaN poroso.