Electric-field effects on defect migration energetics in GaN

Résumé technique : Effets du champ électrique sur l'énergétique de la migration des défauts dans le GaN

Énoncé du problème
La fiabilité des dispositifs à base de nitrure de gallium (GaN) et d'AlGaN, particulièrement dans les environaux de haute puissance et de rayonnement, est compromise par la migration et la reconfiguration des défauts ponctuels. Si la formation de défauts sous irradiation est bien étudiée, le comportement des défauts existants sous l'influence combinée de champs électriques opérationnels et de la température reste moins bien compris. Dans les transistors à haute mobilité d'électrons (HEMT) et les dispositifs de puissance, les champs opérationnels approchent souvent la limite de claquage (~3 MV/cm), ce qui peut induire un transport de défauts menant à des instabilités de la tension de seuil et à une dérive des performances. Les approches de calcul existantes sont confrontées à un compromis : la théorie de la fonctionnelle de la densité (DFT) offre une grande précision mais est limitée à de petites supercellules, tandis que la dynamique moléculaire (MD) classique avec des potentiels interatomiques standards manque de la capacité de modéliser le transfert de charge dynamique et le couplage électrostatique requis pour simuler la migration des défauts sous champs électriques externes.

Méthodologie
Pour remédier à ces limitations, les auteurs ont développé un champ de force réactif ReaxFF spécifique pour le GaN. La méthodologie comprenait trois étapes principales :

1. Développement du champ de force : Un ensemble de données d'entraînement complet a été généré à l'aide de calculs DFT (VASP) dans le cadre de la fonctionnelle PBE. Cet ensemble comprenait des données d'équation d'état (EOS) pour diverses phases de Ga et de GaN (wurtzite, sel gemme, zinc-blende, etc.), des profils moléculaires de N2, et les énergies totales de supercellules contenant des défauts ponctuels neutres (vacances, antisites, paires de Frenkel, interstitiels et défauts de Schottky).
2. Validation : Les paramètres ReaxFF ont été optimisés par rapport à cet ensemble de données DFT. Le champ de force résultant a été validé en comparant les paramètres de réseau, les énergies de cohésion, le module de compressibilité et les énergies de formation/migration des défauts avec des résultats DFT et expérimentaux antérieurs.
3. Simulations de champ électrique : En utilisant le modèle ReaxFF validé dans LAMMPS, les auteurs ont effectué des calculs de chemin de transition (CI-NEB) pour déterminer les barrières de migration pour les défauts ponctuels intrinsèques ($V_{Ga}$, $V_N$, $Ga_{i1,oct}$ et $N_{i2,split}$) sous des champs électriques externes allant de 0 à 0,04 V/Å (0–40 MV/cm). Les simulations ont pris en compte l'équilibrage de charge dynamique (QEq) et les distorsions locales du réseau.

Contributions clés et résultats

• Validation du champ de force : Le modèle ReaxFF développé reproduit avec succès les propriétés structurelles de la DFT (par exemple, les constantes de réseau à 0,3 % près et le module de compressibilité à 1,5 % près) ainsi que l'énergétique des défauts. Bien que certaines déviations existent dans les énergies absolues de formation des défauts (jusqu'à 26 % pour certaines configurations), le modèle capture fidèlement les tendances relatives et les barrières de migration par rapport à la littérature.
• Réponse non linéaire au champ : Contrairement à la simple approximation linéaire dipôle-champ ($\Delta E_B = -\mu_{eff} \cdot E$), l'étude révèle que les champs électriques modifient les barrières de migration de manière anisotrope et non linéaire. La réponse dépend fortement du type de défaut spécifique et de son mécanisme de migration.
  • Vacances : La migration des vacances d'azote ($V_N$) présente une réponse presque linéaire et symétrique au champ, cohérente avec un modèle de dipôle rigide en raison d'une déformation de réseau et d'une variation de charge minimales. Les vacances de gallium ($V_{Ga}$) montrent une légère asymétrie ; la barrière directe augmente tandis que la barrière inverse diminue, mais avec des pentes inégales, ce la est attribué à la taille atomique plus grande du Ga induisant une contrainte locale et une redistribution de charge.
  • Interstitiels : Les défauts interstitiels ($Ga_{i1,oct}$ et $N_{i2,split}$) affichent une non-linéarité et une asymétrie prononcées. Pour $Ga_{i1,oct}$, le champ est assez fort pour déplacer la configuration du point de selle limitant la vitesse, provoquant un croisement de l'énergie de pic dominante de la barrière. Ces défauts impliquent un mouvement coopératif de plusieurs atomes, conduisant à une redistribution complexe des charges partielles et à des réponses de contrainte qui s'écartent considérablement des prédictions linéaires.
• Mécanisme de déviation : Le départ du comportement linéaire provient de la redistribution de charge partielle induite par le champ et de la distorsion du réseau local le long de la coordonnée de réaction. Le moment dipolaire effectif n'est pas fixe mais évolue continuellement à mesure que le défaut migre, un phénomène capturé nativement par le champ de force réactif mais manqué par les modèles à dipôle fixe.

Signification et affirmations
L'article affirme que les champs électriques ne se contentent pas de biaiser le mouvement des défauts de manière linéaire, mais peuvent modifier de manière complexe le paysage de migration des défauts via le couplage charge-réseau.

• Impact cinétique : L'étude souligne que si les réductions de barrière induites par le champ pour les vacances à haute barrière (> 2 eV) pourraient ne pas modifier significativement les taux de diffusion aux échelles de temps des dispositifs, les effets sont cinétiquement importants pour les défauts interstitiels dont les barrières sont proches de 1 eV. Une réduction d'environ 0,05 eV peut augmenter les taux de migration d'un ordre de grandeur à température ambiante, potentiellement en accélérant la redistribution des défauts dans les régions de haut champ des dispositifs en fonctionnement.
• Fiabilité des dispositifs : Ces conclusions suggèrent que dans les régions de haut champ (par exemple, les bords de grille dans les HEMTs), les défauts interstitiels peuvent devenir nettement plus mobiles, contribuant à la génération de pièges et à la dérive de la tension de seuil. L'asymétrie directionnelle de la réponse implique un dérive nette des défauts vers des régions spécifiques du dispositif, pouvant entraîner des problèmes de fiabilité à long terme.
• Portée méthodologique : Les auteurs notent que leurs résultats concernent des défauts neutres en raison de la nature classique du schéma d'équilibrage de charge QEq. Ils reconnaissent que les défauts fortement chargés dans les conditions réelles de fonctionnement (par exemple, $V_{Ga}^{3-}$) pourraient présenter des réponses au champ encore plus fortes, suggérant que les résultats actuels représentent une limite inférieure des effets de champ. Néanmoins, ce travail établit que le couplage champ-défaut dans le GaN est anisotrope et spécifique au défaut, nécessitant des simulations réactives pour des prédictions précises sous des conditions de fonctionnement réalistes.