Simulation de la diode PIN en GaN : modèles de dérive-diffusion, courbes J-V et transport dans les matériaux à large bande interdite.
1. Introduction
Une diode PIN en nitrure de gallium (GaN) est un dispositif semi-conducteur à large bande interdite, doté d’une région de type p, d’une région faiblement dopée (théoriquement intrinsèque) et d’une région de type n. Contrairement à une jonction PN ordinaire, une diode PIN introduit délibérément une zone de dérive importante entre les bornes dopées. Sous polarisation inverse, cette zone peut supporter un champ électrique réparti ; sous polarisation directe, les électrons et les trous injectés permettent à la structure de conduire le courant. Ces propriétés rendent l’architecture PIN pertinente pour les principes fondamentaux de la rectification de puissance, du blocage de la tension inverse et de la conception de dispositifs à large bande interdite.
Le GaN présente, à température ambiante, une largeur de bande interdite d’environ 3,4 eV et peut supporter des champs électriques plus élevés que le silicium. Toutefois, les diodes de puissance GaN pratiques nécessitent un contrôle précis de l’épaisseur de la couche de dérive et du dopage de fond, de la résistance de contact, des défauts du cristal et de la terminaison des bords. Dans cet exemple, la région intrinsèque de GaN de 800 nm est délibérément compacte : elle facilite l’étude de l’électrostatique centrale et de l’injection de porteurs, mais ne constitue pas en soi un redresseur de puissance haute tension performant.
Ce tutoriel utilise les équations couplées de Poisson et de diffusion-dérive d’OghmaNano pour résoudre le problème d’une diode PIN en GaN unidimensionnelle. Au lieu de considérer le courant comme une simple grandeur de circuit, le modèle analyse le profil interne de la bande, les niveaux quasi-Fermi, les concentrations d’électrons et de trous, ainsi que les taux de recombinaison en fonction de la variation de la tension de polarisation. L’objectif est de relier une courbe J–V intuitive, représentant le comportement de redressement, aux processus physiques qui se déroulent à l’intérieur du semi-conducteur.
2. Créer une nouvelle simulation.
Ouvrez OghmaNano et sélectionnez « Nouvelle simulation ». Choisissez les exemples de démonstration sur les semi-conducteurs III–V, puis ouvrez l’exemple de la diode PIN en GaN. Les catégories de simulation disponibles et la fenêtre de sélection des démonstrations sont présentées dans les figures ?? et ??. Si la démonstration est déjà installée, son ouverture devrait afficher la structure en couches dans la fenêtre principale de l’appareil.
La vue principale de la simulation (??) affiche une représentation tridimensionnelle du dispositif. Bien que la visualisation soit en 3D, l’exemple numérique présenté ici ne prend en compte que le transport électrique dans la direction de la croissance de la couche. Ouvrez l’éditeur de couches pour afficher l’épaisseur réelle et la composition de chaque région.
3. Structure en couches et dopage
La démonstration achevée utilise la structure semi-conductrice électriquement active suivante. La couche intrinsèque de 800 nm est la configuration utilisée pour la simulation illustrée ici ; les versions antérieures de l’exemple peuvent contenir une couche intrinsèque de 400 nm et doivent être mises à jour avant de tenter de reproduire ces résultats.
| Couche | Épaisseur | Concentration nominale de dopant | Objectif |
|---|---|---|---|
| p-GaN | 200 nm | Concentration d’accepteurs : 5 × 1017 cm−3 | injection dans la zone p et région anodique |
| GaN intrinsèque | 800 nanomètres | Aucun dopage intentionnel. | Zone de dérive et zone d’appauvrissement faiblement dopées. |
| GaN de type n | 200 nm | Concentration des donneurs : 5 × 1017 cm−3 | région de jonction de type n |
| GaN dopé de type n⁺ | 100 nanomètres | Concentration des donneurs : 1 × 1018 cm−3 | Zone de contact arrière. |
L’empilement actif de semi-conducteurs a une épaisseur de 1,3 µm. Une zone d’air distincte de 200 nm peut apparaître dans l’éditeur de couches pour la visualisation du dispositif et la géométrie optique ; elle ne fait pas partie de l’empilement PIN électriquement actif. Les quatre régions de semi-conducteurs utilisent du GaN, de sorte que le but des interfaces ici est de modifier le dopage, et non la composition de l’alliage de semi-conducteurs.
Ouvrez l’éditeur de contacts (voir la figure ci-dessous) pour examiner l’anode et la cathode, puis ouvrez l’éditeur de dopage (voir la figure ci-dessous). La région p dopée avec des accepteurs et les régions n/n⁺ dopées avec des donneurs doivent apparaître de part et d’autre d’une large région centrale non dopée. Le centre de la structure ne constitue pas une hétérojonction matérielle distincte : sa principale caractéristique est un faible niveau de dopage intentionnel.
Dans le cadre de l’approximation idéale de l’épuisement, une région à faible charge nette peut supporter un champ électrique qui s’étend sur une grande partie de son épaisseur. La tension supportée par une région est liée à l’intégrale spatiale de son champ électrique :
\[V_{\mathrm{region}}=-\int_{y_0}^{y_1} E(y)\,\mathrm{d}y.\]Ainsi, une zone plus large et faiblement dopée peut répartir une tension inverse appliquée sur une plus grande distance. C’est l’une des raisons principales pour lesquelles les structures PIN sont étudiées pour le blocage de la tension. Cela ne signifie pas que la tension de claquage peut être obtenue simplement en multipliant l’épaisseur par un champ critique supposé : les dispositifs réels présentent des champs non uniformes, une charge de fond, des contacts et une concentration du champ aux bords.
4. Modèles électriques, contacts et configuration numérique.
4.1 Paramètres des matériaux GaN
Ouvrez l’éditeur des paramètres électriques pour examiner les paramètres de transport et de recombinaison utilisés pour chaque région de GaN. La simulation est configurée pour charger les propriétés du GaN à partir de la base de données de matériaux OghmaNano, les indicateurs de remplacement de la base de données étant désactivés. Il est préférable de ne pas attribuer de valeurs arbitraires aux mobilités, aux largeurs de bande ou aux affinités électroniques dans le fichier JSON du dispositif. Vérifiez les données de matériaux chargées si vous avez besoin d’une comparaison quantitativement reproductible avec les résultats expérimentaux.
Pour le GaN, les principaux éléments à prendre en compte sont la large bande interdite, la mobilité des électrons et des trous, les statistiques des porteurs libres, la recombinaison radiative et la recombinaison assistée par des pièges. La concentration nominale des accepteurs indiquée dans l’éditeur de dopage ne doit pas être automatiquement associée à la densité des trous libres : les accepteurs de magnésium dans le GaN de type p ne sont pas complètement ionisés à température ambiante. Il convient de déterminer dans quelle mesure le modèle électrique prend en compte l’ionisation des dopants avant d’attribuer une résistance de contact physique ou d’estimer une densité de porteurs mesurée.
4.2 Contacts et conditions aux limites
La démonstration utilise des contacts ohmiques idéaux pour l’anode de type p et la cathode de type n. Cela simplifie l’interprétation en évitant la présence d’un deuxième élément rectificateur à l’interface métal/semiconducteur. Dans les dispositifs GaN fabriqués, en particulier au niveau des contacts p-GaN, l’obtention d’un comportement ohmique à faible résistance peut être difficile et implique la structure métallique, la préparation de la surface et le dopage.
4.3 Grille électrique
L’éditeur de maillage électrique (illustré ci-dessous) affiche la discrétisation le long de l’épaisseur du semi-conducteur. Lorsque vous modifiez la région intrinsèque de 400 nm à 800 nm, assurez-vous que le maillage couvre la nouvelle pile active de 1,3 µm. Un choix initial d’environ 250 à 260 points permet d’obtenir un espacement de l’ordre de 5 nm ; le logiciel peut automatiquement recalculer le maillage en fonction de la configuration. Pour une étude plus quantitative, répétez la simulation avec un maillage plus fin et vérifiez que la courbe J–V et les profils internes restent effectivement inchangés.
4.4 Transport et recombinaison
OghmaNano couple l’équation de Poisson électrostatique aux équations de continuité des électrons et des trous. Dans une dimension, les équations à l’état stationnaire peuvent être représentées schématiquement comme suit :
\[\frac{\mathrm{d}}{\mathrm{d}y}\!\left(\varepsilon\frac{\mathrm{d}\phi}{\mathrm{d}y}\right)=-\rho,\qquad \frac{\mathrm{d}J_n}{\mathrm{d}y}=q(R-G),\qquad \frac{\mathrm{d}J_p}{\mathrm{d}y}=-q(R-G),\]où \(\phi\) représente le potentiel électrostatique, \(\rho\) la densité de charge spatiale locale, \(J_n\) et \(J_p\) les densités de courant des électrons et des trous, \(G\) le taux de génération et \(R\) le taux de recombinaison net. En utilisant ces conventions de signe, la densité de courant totale \(J_n+J_p\) reste constante dans un dispositif unidimensionnel, à l’état stationnaire et après convergence.
Le modèle de recombinaison peut inclure des contributions de recombinaison directe bande à bande (libre à libre) et de recombinaison de Shockley–Read–Hall (SRH). Leur importance relative dépend de la qualité du matériau, du niveau d’injection et des paramètres du modèle. La bande interdite directe du GaN rend la recombinaison radiative physiquement possible, mais sa présence dans un modèle n’implique pas qu’un dispositif particulier représenté sur un graphique soit une LED efficace.
5. Exécution de la simulation J–V en mode sombre.
Réglez la plage de variation de la tension de −2,0 V à +5,0 V, par incréments de 0,025 V. Une tension négative polarise inversement la structure PIN, tandis qu’une tension positive la polarise en mode direct. La configuration de polarisation souhaitée est la suivante :
"Vstart": -2.0,
"Vstop": 5.0,
"Vstep": 0.025Dans la fenêtre principale, sélectionnez « Exécuter la simulation ». Pendant que le calcul s’effectue, utilisez l’onglet « Terminal » pour surveiller l’avancement et la convergence. Une fois le calcul terminé, ouvrez l’onglet « Résultats » et sélectionnez « jv.csv » pour afficher la densité de courant en fonction de la tension appliquée. Le visualiseur de résultats, le terminal du solveur et la courbe J–V résultante sont présentés ci-dessous.
En polarisation inverse, le courant simulé est faible par rapport à sa valeur en polarisation directe. Lorsque la polarisation positive augmente, l’injection de porteurs augmente rapidement et la diode s’active. Pour des tensions directes plus élevées, la courbe devient plus linéaire sur un axe de densité de courant linéaire, ce qui indique que le courant est de plus en plus influencé par le transport résistif et d’autres non-idéalités dans la structure simulée. Une pente qui s’aplatit sur un axe de courant logarithmique ne doit pas être interprétée automatiquement comme une saturation du courant : une augmentation du courant qui est à peu près linéaire apparaît également progressivement plus plate sur un graphique logarithmique.
Pour une comparaison simple avec la théorie des diodes, la caractéristique de polarisation directe faible à modérée peut être mise en relation avec
\[J\approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right],\]où \(J_0\) représente le courant de saturation effectif et \(n\) le facteur d’idéalité effectif. Une diode PIN peut s’écarter considérablement d’un régime idéal de Shockley unique, car la recombinaison dans la zone de déplétion, l’injection dans la couche intrinsèque, l’injection à fort niveau et la résistance série peuvent influencer le courant. L’expression simple doit donc être considérée comme une approximation diagnostique plutôt que comme l’équation fondamentale résolue par OghmaNano.
La capacité de blocage en tension inverse n’est pas équivalente à la rupture en avalanche observée. Un balayage jusqu’à -2 V ne permet de tester que le comportement en faible polarisation inverse. La tension à laquelle un dispositif de puissance en GaN fabriqué se rompt dépend de la physique des champs élevés, des défauts du matériau, des chemins de fuite, de la géométrie et de la terminaison des bords, et ne peut être déduite de cette seule démonstration.
6. Analyse des résultats de la simulation interne.
Après l’exécution du balayage J–V, le visualiseur d’instantanés d’OghmaNano permet d’afficher les résultats en fonction de la position, pour des points de polarisation sélectionnés. Les captures d’écran ci-dessous illustrent quatre positions dans le balayage enregistré, en commençant par l’instantané le plus ancien et en terminant par le plus récent. Leurs noms (-0 à -3) indiquent les positions des captures d’écran, et non les tensions. Pour effectuer des comparaisons quantitatives, consultez la valeur de la tension sélectionnée dans le visualiseur d’instantanés.
6.1 Bords de bande et niveaux quasi-Fermi
Ouvrez le visualiseur d’instantanés et affichez les courbes de Ec.csv, Ev.csv, Fn.csv et Fp.csv ensemble. Elles représentent respectivement les bords des bandes de conduction et de valence, ainsi que les niveaux quasi-Fermi des électrons et des trous. Les figures ?? à ?? illustrent leur évolution en fonction de la variation de la tension appliquée.
En polarisation inverse, la majeure partie de la variation du potentiel devrait se produire dans la région intrinsèque faiblement dopée, ce qui entraîne une pente de bande relativement étendue. En polarisation directe, la barrière d’injection effective est réduite et une séparation des quasi-niveaux de Fermi des électrons et des trous, qui dépend de la polarisation, est créée. Les gradients spatiaux de ces niveaux sont importants car ils déterminent le flux de courant ; la courbure de la bande et la séparation des quasi-niveaux de Fermi doivent être considérées ensemble plutôt que d’être interprétées séparément.
À titre de vérification physique de base, la distance entre \(E_c\) et \(E_v\) dans les couches de GaN doit correspondre à la valeur de la bande interdite du GaN indiquée dans la base de données des matériaux, et non à une valeur par défaut résiduelle provenant d’un autre semi-conducteur. Étant donné qu’il s’agit d’un dispositif PIN constitué d’un matériau homogène, on ne s’attend pas à ce que des décalages de bande induits par une variation de composition se produisent brusquement aux transitions p/i/n.
6.2 Densités d’électrons et de trous.
Tracez les courbes Q_nfree.csv et Q_pfree.csv (ou les résultats relatifs à la densité des porteurs électroniques et des trous, tels qu’ils sont présentés par la version installée). À faible niveau d’injection, les densités des porteurs majoritaires sont concentrées dans leurs régions dopées respectives. La région centrale non dopée contient une concentration d’équilibre des porteurs beaucoup plus faible. En cas de forte polarisation directe, les porteurs injectés occupent une plus grande partie de la région intrinsèque, et la recombinaison peut augmenter considérablement.
Le changement brutal de la concentration nominale des dopants n’implique pas nécessairement un profil tout aussi brutal de la concentration des porteurs de charge libres. La concentration des porteurs de charge est déterminée de manière cohérente à partir des énergies de bande, des niveaux de Fermi, de la densité d’états et de la neutralité de charge. Ceci est particulièrement pertinent pour la couche de type p en GaN si le modèle de matériau utilise une approche simplifiée de l’ionisation incomplète des accepteurs.
6.3 Recombinaison dans la région intrinsèque
Examinez les graphiques de recombinaison spatiale aux mêmes points de polarisation. Les figures de recombinaison présentées dans cet exemple illustrent comment la zone de chevauchement et de recombinaison significatifs entre les électrons et les trous se développe à mesure que l’injection directe augmente. Interprétez la hauteur et la largeur de chaque profil représenté ensemble : un pic local important ne représente pas nécessairement la contribution la plus importante au taux de recombinaison global du dispositif.
Dans le cas d’une transition directe entre deux bandes d’énergie, le taux d’émission radiative net suit généralement la formule suivante :
\[R_{\mathrm{rad}}=B\left(np-n_i^2\right),\]où \(B\) est le coefficient de recombinaison du matériau. La recombinaison de type SRH dépend des paramètres des défauts et peut présenter une distribution spatiale différente. Pour comparer les mécanismes, représentez-les séparément s’ils sont tous les deux disponibles et activés dans l’appareil chargé. Ne présumez pas qu’une représentation graphique de la recombinaison totale sans légende indique un seul mécanisme microscopique.
6.4 Taux de recombinaison/de continuité des électrons et des trous.
Les graphiques supplémentaires illustrant les taux d’électrons et de trous sont utiles pour déterminer les zones où les populations de porteurs varient le plus rapidement en raison de la recombinaison et de l’injection. Les quatre captures d’écran ci-dessous montrent comment ces profils calculés évoluent en fonction de la polarisation. Si le taux affiché est signé ou contient à la fois des contributions de génération et de recombinaison, utilisez la légende du graphique et l’en-tête du fichier de données pour déterminer la bonne interprétation physique avant de comparer les valeurs.
6.5. La méthode de conservation actuelle comme modèle de vérification.
Pour tout point de polarisation en régime permanent, représentez graphiquement la densité de courant des électrons \(J_n\) et la densité de courant des trous \(J_p\), si ces valeurs sont disponibles. Chaque composante peut varier en fonction de la position, car les porteurs se recombinent, mais le total \(J_n+J_p\) doit rester constant dans l’espace dans un dispositif 1D convergent, sans termes supplémentaires d’accumulation de charges. Il s’agit d’une vérification particulièrement utile lorsque une petite ondulation apparente dans une courbe de faible courant pourrait être interprétée à tort comme un effet physique.
7. Exercices proposés.
8. Hypothèses et limites physiques
Cet exemple est conçu comme une simulation éducative à grande échelle du phénomène de dérive-diffusion dans le GaN. Il offre une première approche, pertinente sur le plan physique, de l’électrostatique des structures PIN et de l’injection de porteurs, mais il ne s’agit pas d’un modèle d’appareil commercial calibré. En particulier, les résultats doivent être interprétés en tenant compte des réserves suivantes.
- Contacts idéalisés : les barrières de contact métal/GaN et la résistance de contact pratique du p-GaN ne sont pas explicitement validées.
- Dopage nominal : l’activation des accepteurs de type p et la charge résiduelle de fond dans le GaN, qui est théoriquement intrinsèque, peuvent modifier la distribution physique des charges.
- Dépassement du champ électrique élevé : la variation de la tension et du courant ne permet pas de déterminer le seuil de claquage en avalanche, la tension de blocage maximale ou la fiabilité sous forte polarisation inverse.
- Transport unidimensionnel : le modèle ne prend pas en compte l’effet de concentration du courant latéral, les fuites en périphérie, les plaques de champ et la géométrie de terminaison du dispositif.
- Défauts et recombinaison : pour obtenir des prédictions quantitatives réalistes concernant les pertes et la recombinaison, il est nécessaire de calibrer les paramètres des défauts présents dans le matériau et à l’interface.
- Effets thermiques et dynamiques : l’auto-échauffement, les transitoires de commutation et les effets liés au stockage de charges par les porteurs minoritaires doivent être évalués séparément, le cas échéant.
Ces limites ne diminuent pas la valeur de cet exemple en tant que point de départ. Le schéma de superposition des couches et la courbe J–V illustrent ce qui se passe lorsqu’une région de dérive faiblement dopée est insérée dans une diode à semi-conducteur à large bande interdite. Ils fournissent également une base de référence reproductible pour tester des modèles de transport et de recombinaison plus avancés.
9. Prochaines étapes
Après avoir étudié la diode PIN en GaN, comparez sa distribution de champ avec celle d’une jonction PN en GaN classique, puis étudiez l’influence de l’épaisseur de la couche intrinsèque et du dopage résiduel sur la rectification. Des travaux plus approfondis pourraient introduire des états de piège calibrés, des contacts non idéaux et des mécanismes de transport à champ élevé afin d’examiner les limites du modèle simple de dérive-diffusion.