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

1. Problématique et Contexte

Les réflecteurs de Bragg distribués (DBR) en nitrure de gallium (GaN) sont des composants clés pour les dispositifs optoélectroniques (LED, lasers VCSEL, sources de photons uniques). Ils sont fabriqués par gravure électrochimique sélective de structures épitaxiales alternant des couches de GaN fortement dopées (n-type) et des couches non intentionnellement dopées (NID).

Le défi principal réside dans l'accès de l'agent de gravure aux couches dopées enfouies sous la surface (qui est généralement une couche NID non conductrice). Deux mécanismes ont été identifiés :

• Le modèle « Brochette » (Kebab model) : Chaque dislocation de vissage (TD) agit comme un pipeline vertical continu, permettant à l'agent de gravure de traverser toutes les couches NID pour atteindre chaque couche dopée, créant des pores alignés verticalement.
• Le modèle « Cascade » : Des données expérimentales récentes (tomographie FIB) montrent que le chemin de gravure est plus complexe. L'agent de gravure peut descendre verticalement le long d'une dislocation, puis se propager latéralement à travers une couche dopée pour atteindre une autre dislocation inactive, qui devient alors active. Ce mécanisme crée des structures poreuses en « cascade » où certaines dislocations s'activent et se désactivent à différents niveaux de profondeur.

La physique exacte régissant ces transitions (pourquoi certaines dislocations s'activent ou se désactivent) n'est pas entièrement comprise, et les modèles existants peinent à prédire la morphologie des pores et les courants de gravure (chronoampérométrie) en fonction de la tension appliquée.

2. Méthodologie

Les auteurs ont développé une simulation stochastique en Python (utilisant la bibliothèque Numba pour l'efficacité) pour modéliser ce processus de gravure.

• Modélisation : La structure DBR est représentée sous forme d'un tableau de pixels (5 nm x 5 nm). Les pixels peuvent être de l'agent de gravure, du GaN NID, du GaN dopé ou des dislocations (TD).
• Règles de probabilité :
  • Les pixels NID ont une probabilité de gravure nulle.
  • Les pixels dopés et les dislocations ont des probabilités de gravure finies et distinctes : $P_{doped}$ (probabilité de gravir la couche dopée) et $P_{disloc}$ (probabilité de gravir la région de la dislocation).
  • Une contrainte empêche la fusion de pores adjacents pour respecter les observations microscopiques de domaines de porosité non intersectés.
• Algorithme : À chaque itération, l'algorithme recherche les voisins des pixels déjà gravés. Si un pixel adjacent a une probabilité de gravure > 0, il est gravé selon une tirage aléatoire basé sur cette probabilité.
• Corrélation expérimentale :
  • Chaque pixel gravé génère une unité de charge simulée, permettant de tracer des courbes de courant en fonction du temps (simulant la chronoampérométrie).
  • La simulation est comparée à des échantillons réels gravés à différentes tensions (6 V à 12 V) et à des structures avec des épaisseurs de couches variables.
  • Des métriques quantitatives (« K5 » pour les dislocations traversant toutes les couches, « C » pour les cascades) sont utilisées pour analyser la morphologie simulée et la comparer aux données FIB expérimentales.

3. Contributions Clés

• Développement d'un modèle stochastique : Création d'un outil de simulation capable de reproduire à la fois la dynamique temporelle (courant) et la morphologie spatiale (pores) de la gravure médiée par les dislocations.
• Validation du modèle « Cascade » : La simulation démontre que le modèle « cascade » émerge naturellement lorsque les probabilités de gravure latérale et verticale sont ajustées, confirmant les observations expérimentales précédentes.
• Lien entre tension appliquée et probabilités : L'article établit une corrélation qualitative où l'augmentation de la tension de gravure expérimentale correspond à une augmentation de la probabilité de gravure des couches dopées ($P_{doped}$) dans la simulation.
• Généralisation du modèle : Démonstration que le même modèle peut être appliqué à des structures non uniformes (épaisseurs de couches variables) et à des géométries alternatives (gravure depuis des tranches latérales).

4. Résultats Principaux

A. Comparaison Chronoampérométrique

• Oscillations de courant : La simulation reproduit les oscillations observées expérimentalement dans les courbes de courant, qui correspondent au passage successif des couches dopées.
• Effet de la tension :
  • À haute tension (ou $P_{doped}$ élevé), les oscillations sont nettes et régulières, indiquant une gravure plus « couche par couche ».
  • À basse tension (ou $P_{doped}$ faible), les oscillations s'atténuent et disparaissent. La simulation montre que cela est dû à une gravure moins uniforme où les fronts de gravure de différentes dislocations ne sont pas synchronisés (certaines avancent plus vite, créant un mélange de gravure de plusieurs couches simultanément).
• Décalage initial : La simulation ne reproduit pas le pic de courant initial observé expérimentalement (dû à la gravure des bords/arrières de l'échantillon), mais ce phénomène est éliminé expérimentalement par un revêtement protecteur, validant ainsi le modèle de surface.

B. Morphologie des Pores et Statistiques

• Transition Kebab/Cascade : La simulation prédit correctement que l'augmentation de $P_{doped}$ (tension élevée) favorise la formation de structures « Kebab » (K5) et réduit la proportion de structures « Cascade » (C). Inversement, à basse probabilité, les cascades dominent.
• Limites du modèle :
  • La simulation ne prédit pas la formation de pores coniques dans les couches NID au-dessus des couches dopées (observés expérimentalement à haute tension), car elle suppose une probabilité nulle de gravure du NID loin des dislocations.
  • Elle ne reproduit pas parfaitement la réduction de l'efficacité de gravure des couches profondes observée expérimentalement à basse tension.

C. Structures à Épaisseurs Variables

• La simulation a été testée sur des échantillons où l'épaisseur des couches dopées varie linéairement (échantillons B et C).
• Elle reproduit avec succès les tendances des courbes de courant (amplitude des pics), montrant que le modèle est applicable au-delà des DBR standards et peut aider à interpréter des structures complexes.

5. Signification et Implications

• Compréhension des mécanismes : Bien que les détails physiques (tunneling Zener, états de pièges) restent à élucider, la simulation suggère que la tension appliquée affecte différemment la gravure des couches dopées et des dislocations. L'augmentation de la tension favorise la gravure des couches dopées, rendant le processus plus « couche par couche » et réduisant la complexité des cascades.
• Outil de prédiction : Ce modèle stochastique offre un moyen rapide et peu coûteux de prédire le comportement de la gravure électrochimique pour divers designs de structures et conditions de processus, sans avoir besoin de simulations physiques complexes à l'échelle atomique.
• Applicabilité étendue : Le modèle suggère que même dans les approches de gravure par tranches latérales (lithographie + RIE), les dislocations peuvent jouer un rôle actif dans la morphologie des pores, affectant potentiellement des propriétés comme la biréfringence du matériau.
• Optimisation des dispositifs : En permettant de comprendre comment contrôler le rapport entre structures « Kebab » et « Cascade », ce travail aide à optimiser la fabrication de DBR pour des applications spécifiques (extraction de lumière, confinement optique).

En résumé, cet article présente un outil de simulation robuste qui valide le modèle de « cascade » pour la gravure du GaN poreux et fournit un cadre pour interpréter les données expérimentales complexes en fonction des paramètres de processus.