Efficient Analysis of Carrier Transport and TM-TE Emission in AlGaN UVC LEDs via Multi-band Localization Landscape Theory

Résumé technique : Analyse efficace du transport des porteurs et de l'émission TM-TE dans les LED UVC à base d'AlGaN via la théorie du paysage de localisation multi-bandes

Énoncé du problème
Les diodes électroluminescentes (LED) à base de nitrure d'aluminium et de gallium (AlGaN) émettant dans l'ultraviolet profond (UVC), opérant dans la plage 220–250 nm, sont confrontées à des défis d'efficacité importants. La haute teneur en aluminium (>60 %) requise pour ces longueurs d'onde entraîne un faible confinement des trous et une dominance de la bande de valence $|Z\rangle$ comme état fondamental. Cette structure de bandes résulte en une forte polarisation transverse magnétique (TM), ce qui entrave l'extraction de la lumière, et exacerbe les problèmes de transport de porteurs en raison de la faible mobilité des trous et de l'énergie d'ionisation élevée des accepteurs. De plus, les fluctuations aléatoires d'alliage ternaire dans l'AlGaN créent des variations spatiales tridimensionnelles (3D) des bords de bandes, de la déformation et des champs de polarisation, formant un paysage de potentiel accidenté qui ne peut être modélisé avec précision par des approximations de puits quantiques moyennés en une dimension.

Bien que les méthodes de simulation conventionnelles combinant des modèles $k \cdot p$ avec des solveurs de Schrödinger, de Poisson et de diffusion-dérive puissent traiter ces complexités, elles sont trop coûteuses en termes de calcul pour les simulations de dispositifs 3D à grande échelle, particulièrement pour la résolution du problème de valeur propre de la structure de bande de valence complète à 6 bandes sur l'ensemble de la pile du dispositif.

Méthodologie
Les auteurs proposent un cadre de simulation hybride qui intègre la théorie du paysage de localisation (LL) multi-bandes avec le formalisme de Wigner–Weyl (WW) et un solveur 3D de diffusion-dérive-continuité de courant (3D-DDCC).

1. Modèle de paysage de localisation (LL) : Au lieu de résoudre le problème de valeur propre $k \cdot p$ 3D complet, l'approche LL le remplace par une équation aux dérivées partielles elliptique linéaire ($\hat{H}u = 1$). Cela génère un potentiel de confinement effectif ($W = 1/u$) qui capture le confinement quantique, l'effet tunnel et la localisation induite par le désordre.
2. Intégration multi-bandes et de déformation : Pour traiter la physique spécifique des LED UVC, le modèle résout les bandes $|X \pm iY\rangle$ (polarisées TE) et $|Z\rangle$ (polarisées TM) séparément. Crucialement, il incorpore des potentiels de déformation induits par la contrainte dérivés d'un solveur de contrainte-déformation par éléments finis 3D. Cela permet au modèle de rendre compte de la manière dont la contrainte de compression déplace les positions énergétiques relatives des bandes de valence, modifiant l'état fondamental de $|Z\rangle$ vers $|X \pm iY\rangle$ sous certaines conditions.
3. Formalisme de Wigner–Weyl : Ce formalisme de l'espace des phases semi-classique est utilisé pour calculer les taux d'émission locaux et les rapports de polarisation optique. Il utilise les bandes interdites effectives spatialement variables dérivées des potentiels LL pour calculer les spectres d'électroluminescence (EL) sans nécessier de solutions globales d'états propres.
4. Validation : Les résultats LL-WW sont comparés à un solveur $k \cdot p$ à 6 bandes 3D complet (résolu uniquement pour la région active après convergence en raison des limites de calcul) et comparés à des données expérimentales publiées.

Contributions clés et résultats

• Efficacité de calcul : La contribution primaire est une réduction drastique du temps de simulation. Alors que la méthode $k \cdot p$ 3D nécessitait environ 2 240 heures pour résoudre 1 200 états propres de trous dans un seul puits quantique, la méthode LL-WW a complété la simulation complète du dispositif (incluant Poisson et la diffusion-dérive) en environ 0,5 heure sur un matériel comparable. Cela représente une réduction de plus d'un ordre de grandeur.
• Précision de la prédiction de la polarisation : Le modèle LL reproduit avec succès le rapport de polarisation dépendant de la longueur d'onde (croisement TE/TM) observé dans les expériences et prédit par la méthode $k \cdot p$. Il capture avec précision le comportement dépendant du tampon : les tampons AlN induisent une contrainte de compression plus forte, décalant l'ordre des bandes de valence pour favoriser l'émission TE aux longueurs d'onde plus courtes, tandis que les tampons AlGaN maintiennent des rapports de polarisation TM plus élevés sur une plage plus large.
• Analyse du transport de porteurs induit par la déformation : L'étude révèle que l'inclusion des potentiels de déformation modifie considérablement la physique du transport des trous.
  • Puits quantiques moins profonds : La contrainte de compression abaisse les bords des bandes de valence, rendant effectivement les puits quantiques de trous moins profonds.
  • Amélioration de l'injection des trous : Ce confinement moins profond, combiné à la masse effective plus légère de la bande $|Z\rangle$ dans la direction de croissance sous contrainte, facilite le transport vertical des trous à travers l'empilement des puits quantiques multiples (MQW).
  • Suppression du débordement d'électrons : Une meilleure injection de trous conduit à un meilleur équilibre électrons-trous, réduisant la fuite d'électrons (débordement) au-delà de la région active. La fraction de trous résidant dans la région de barrière est passée d'environ 48 % (sans les termes de déformation) à environ 75 % (avec les termes de déformation), indiquant une distribution de porteurs plus uniforme.
• Effets du désordre : Le modèle simule efficacement l'impact des fluctuations aléatoires d'alliage sur la localisation des porteurs et le croisement de polarisation TE-TM, montrant que les variations locales de composition peuvent déplacer le point de croisement de la polarisation par rapport aux prédictions du cristal virtuel.

Signification
L'article affirme que l'approche basée sur le paysage de localisation (LL) multi-bandes offre une solution rapide et fiable pour comprendre et optimiser les performances des LED UVC. En conservant les mécanismes physiques clés — tels que le mélange de bandes, la localisation induite par le désordre et le réordonnancement des bandes dépendant de la déformation — tout en évitant le coût de calcul élevé des solveurs de valeurs propres $k \cdot p$ complets, cette méthode permet une exploration rapide de l'espace de conception. Cela est particulièrement significatif pour l'optimisation des émetteurs UV profonds où la déformation, le désordre d'alliage et les interactions multi-bandes sont des déterminants critiques de l'efficacité du dispositif et des caractéristiques de polarisation. Le cadre permet aux chercheurs d'étudier l'interaction entre les caractéristiques électriques (structure de bandes, confinement des porteurs) et les propriétés optiques (spectres d'émission, polarisation) de manière informatiquement réalisable, facilitant l'optimisation des couches tampons et des structures de puits quantiques pour des LED UVC à haute efficacité.