Simulation d’une diode à jonction PN en GaN : modèles de dérive-diffusion, courbes J-V et recombinaison.
1. Introduction
Une diode à jonction PN en nitrure de gallium (GaN) est constituée de régions de type p et de type n du même semi-conducteur, qui sont réunies pour former une jonction de redressement. Les électrons et les trous se redistribuent à l’interface, ce qui crée une zone de déplétion et une barrière de potentiel intégrée. Sous polarisation directe, la barrière diminue et le courant augmente brusquement ; sous polarisation inverse, la jonction empêche l’injection de porteurs. La réponse asymétrique résultante en termes de courant et de tension constitue le principe de fonctionnement de base d’un redresseur à semi-conducteurs.
Le GaN est un semi-conducteur de type III-V à large bande interdite, avec une bande interdite d’environ 3,4 eV à température ambiante. Son champ électrique critique élevé et sa faible concentration de porteurs intrinsèques en font un matériau intéressant pour la conversion de puissance, la rectification électrique et l’électronique à haute température. La jonction homojonction compacte p-GaN/n-GaN présentée dans ce tutoriel offre un moyen clair d’étudier comment la barrière intégrée contrôle le flux de courant et comment les propriétés du matériau GaN influencent une jonction de rectification.
Dans ce tutoriel, OghmaNano est utilisé pour étudier la rectification électrique, le transport des porteurs par dérive et diffusion, la courbure des bandes d’énergie et la recombinaison. Ses équations de Poisson et de continuité des électrons/trous, résolues de manière couplée, permettent de calculer le courant qui traverse la diode, ainsi que les profils internes des porteurs et le champ électrostatique en fonction de la position et de la tension appliquée. Cela permet de relier la courbe caractéristique courant-tension (J–V) aux processus physiques qui se produisent au niveau de la jonction.
Les fines couches, disposées verticalement les unes sur les autres, sont affichées dans la fenêtre principale de l’appareil OghmaNano (??). La visualisation est en trois dimensions, mais la simulation électrique présentée ici est unidimensionnelle, le long de la direction de croissance des couches.
2. Créer une nouvelle simulation.
Dans la fenêtre principale d’OghmaNano, sélectionnez « Nouvelle simulation ». La première boîte de dialogue affiche les catégories de simulation disponibles (??). Ouvrez les exemples III–V/GaN et sélectionnez la démonstration de la diode PN en GaN (??). L’ordre exact des exemples peut varier selon les versions d’OghmaNano.
Après avoir chargé l’exemple, ouvrez l’éditeur de couches à partir du ruban de l’appareil. Cela permet d’afficher les régions semi-conductrices, leurs épaisseurs et les entrées sélectionnées dans la base de données des matériaux. L’appareil est une hétérojonction PN, les couches p-GaN et n-GaN formant la jonction électriquement active.
3. Analyse de la structure du dispositif, des contacts et du dopage.
3.1 Structure en couches
L’empilement électrique comprend trois régions en GaN. Les couches de GaN de type p et de GaN de type n forment la jonction de redressement ; la couche de GaN⁺ fournit une région plus fortement dopée, adjacente à la cathode. L’épaisseur totale du semi-conducteur est de 700 nm, sans compter l’air environnant utilisé pour la visualisation. L’éditeur de couches est présenté dans ??.
| Couche | Épaisseur | Dopage nominal | Fonction dans le domaine de l’électricité. |
|---|---|---|---|
| GaN dopé au magnésium (p-GaN) | 250 nm | \(N_A=5\times10^{17}\,\mathrm{cm^{-3}}\) | Côté de la jonction PN où la concentration de trous est la plus élevée. |
| n-GaN (dopé au silicium) | 350 nm | \(N_D=5\times10^{17}\,\mathrm{cm^{-3}}\) | Côté de la jonction PN où les électrons sont majoritaires. |
| GaN dopé de type n⁺ | 100 nanomètres | \(N_D=1\times10^{18}\,\mathrm{cm^{-3}}\) | Région du côté cathode avec une concentration de dopants plus élevée. |
3.2 Contacts
L’exemple utilise des conditions aux limites ohmiques idéales au niveau de l’anode de type p et de la cathode de type n. Ces conditions aux limites permettent de se concentrer sur la réponse électrique interne de la jonction. Dans les diodes GaN fabriquées, la résistance de contact et l’activation des accepteurs de magnésium influencent la chute de tension et l’injection de courant ; ces deux paramètres sont utiles pour des études plus approfondies.
3.3 Profil de dopage et zone de déplétion
Ouvrez l’éditeur « Doping / Ions » dans le ruban « Electrical » pour afficher les distributions d’accepteurs et de donneurs (??). Une seule transition abrupte entre le p-GaN dopé avec des accepteurs et le n-GaN dopé avec des donneurs définit la jonction PN active. La région n⁺ introduit une deuxième étape dans la concentration de donneurs près de la cathode.
À l’équilibre thermique, la jonction développe une zone de déplétion, car les porteurs de charge mobiles se diffusent et laissent derrière eux des dopants ionisés. Dans le cadre de l’approximation de la zone de déplétion, la largeur caractéristique dépend de la permittivité, du dopage de chaque côté et du potentiel maintenu par la jonction. Pour une jonction abrupte classique, la relation bien connue est la suivante :
\[ W \simeq \sqrt{\frac{2\varepsilon_s}{q}\left(\frac{1}{N_A}+\frac{1}{N_D}\right)(V_{\mathrm{bi}}-V)}. \]Ici, \(V\) représente la tension appliquée à la jonction, \(V_{\mathrm{bi}}\) représente le potentiel intégré et \(\varepsilon_s\) représente la permittivité du GaN. Cette équation constitue une référence utile pour interpréter le champ électrique auto-cohérent calculé par OghmaNano, en particulier lorsqu’on étudie l’influence du dopage et de la tension appliquée sur la jonction.
4. Paramètres électriques et modèles de matériaux.
4.1 Modèle de matériau GaN Lua
Les trois couches semi-conductrices actives électriquement font référence au modèle de matériau Inorganic/GaN/GaN.lua. La base de données des matériaux électriques fournit des données telles que la largeur de bande interdite, l’affinité électronique, la densité d’états, la permittivité, les mobilités et les paramètres de recombinaison. Dans l’exemple JSON, les paramètres db_overrides pertinents sont désactivés, de sorte que les valeurs utilisées doivent être celles renvoyées par le modèle Lua plutôt que les valeurs inactives visibles ailleurs dans le fichier enregistré. La largeur de bande interdite du GaN doit être d’environ 3,4 eV à température ambiante ; vérifiez cette valeur dans les profils de bande affichés ou dans les paramètres de matériau calculés.
Les captures d’écran des vues relatives au code de matériau et à l’éditeur de paramètres sont regroupées ci-dessous. Certains noms de fichiers des captures d’écran d’origine conservent les anciennes étiquettes GaAs/AlGaAs. Cependant, la simulation active décrite dans ce tutoriel utilise exclusivement le GaN. Dans le logiciel, vérifiez le chemin d’accès Lua du matériau GaN plutôt que de déduire le matériau à partir du nom du fichier image.
4.2 Transport par transporteur
OghmaNano résout l’équation de Poisson pour les semi-conducteurs, ainsi que les équations de continuité pour les électrons et les trous, en utilisant des expressions de courant basées sur les phénomènes de dérive et de diffusion. Les densités de courant des électrons et des trous dépendent du champ local, de la mobilité et du gradient de concentration des porteurs :
\[ J_n=q\mu_n n E+qD_n\frac{dn}{dy},\qquad J_p=q\mu_p p E-qD_p\frac{dp}{dy}. \]Ces équations décrivent deux mécanismes complémentaires : la dérive due au champ électrique interne et la diffusion due aux gradients de concentration des porteurs de charge. À l’état stationnaire, leur interaction détermine le courant net de la diode. Le modèle permet également de déterminer les niveaux quasi-Fermi des électrons et des trous, dont la variation offre une perspective différente sur les forces motrices des porteurs de charge.
4.3 Recombinaison
Les électrons et les trous injectés peuvent se recombiner à l’intérieur du GaN. L’exemple comprend les processus de recombinaison configurés dans le modèle de matériau, notamment les termes de Shockley–Read–Hall (SRH) médiatisés par des défauts, lorsque ces derniers sont activés par les données de matériau pertinentes. L’expression SRH idéale habituelle est la suivante :
\[ R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}. \]La durée de vie des porteurs dans les régions de recombinaison dépend de la densité des pièges, de la vitesse thermique des porteurs et de la section efficace de capture. Ces paramètres déterminent l’efficacité avec laquelle les porteurs injectés se recombinent et, par conséquent, influencent la forme et l’amplitude de la courbe caractéristique courant-tension. L’analyse du taux de recombinaison spatial permet de déterminer où se concentre la perte de porteurs au sein de la jonction.
4.4 Maillage numérique
La grille électrique s’étend sur la structure en GaN de 700 nm et comporte 150 points le long de la direction de croissance, ce qui donne un espacement moyen d’environ 4,7 nm. Il s’agit d’une résolution de départ acceptable pour une démonstration pédagogique. L’éditeur de grille (??) vous permet d’examiner la distribution des nœuds et de l’affiner lors de l’étude de gradients abrupts dans la zone d’appauvrissement.
5. Lancer la simulation J–V.
L’exemple est configuré pour une variation de tension allant de 0 V à +5 V, par incréments de 0,025 V. Lancez la simulation à l’aide de la commande « Exécuter la simulation ». Pendant la simulation, la fenêtre affiche l’avancement numérique et les informations sur la convergence (??). Une fois la simulation terminée, l’onglet « Résultats » contient les résultats générés, notamment jv.csv et les données relatives au dispositif.
Ouvrez jv.csv pour examiner la densité de courant en fonction de la tension appliquée. La vue linéaire (??) et la vue logarithmique (??) présentent la même réponse électrique, mais de manière complémentaire. Le graphique linéaire montre comment le courant augmente à forte polarisation directe ; le graphique logarithmique permet de mieux visualiser le début de l’injection et la réponse à faible courant.
5.1 Interprétation de la courbe de polarisation directe.
À faible polarisation directe, la barrière de potentiel limite fortement l’injection des porteurs minoritaires, de sorte que le courant net reste faible. Lorsque la tension appliquée augmente, la barrière de jonction diminue, les électrons et les trous traversent plus facilement la jonction PN et le courant augmente rapidement. Sur un graphique semi-logarithmique, une région quasi exponentielle est souvent décrite localement à l’aide de l’équation de la diode.
\[ J\simeq J_0\left[\exp\left(\frac{qV}{n_{\mathrm{id}}k_{\mathrm B}T}\right)-1\right], \]où \(J_0\) représente le courant de saturation effectif et \(n_{\mathrm{id}}\) le facteur d’idéalité. Cette équation est un outil d’interprétation utile, mais elle ne constitue pas un modèle indépendant de l’ensemble du dispositif simulé. En particulier, une variation de la pente aux courants plus élevés peut indiquer une résistance série, des effets d’injection ou des changements de recombinaison. La forte augmentation du courant observée dans les captures d’écran fournies est cohérente avec la rectification en polarisation directe, mais un facteur d’idéalité spécifique ne doit être mentionné qu’après avoir ajusté un intervalle de tension approprié.
5.2 Étude du comportement de la jonction en polarisation inverse.
La courbe de référence courant-tension (J–V) s’étend de 0 V à +5 V et illustre l’évolution de la conduction directe. Pour un exercice complémentaire, réglez une tension de départ négative et examinez les modifications de la courbure de la bande, de la largeur de la zone de déplétion et de la distribution du champ électrique en polarisation inverse. Cela permet d’établir une comparaison utile avec l’injection directe et illustre la manière dont la jonction PN réagit à la polarité de la tension appliquée.
6. Analyse des captures d’écran de la simulation : bandes, densités de porteurs, courants et recombinaison.
Le visualiseur d’instantanés affiche les données internes aux points sélectionnés le long de la courbe de balayage de tension. Les quatre images de chaque groupe ci-dessous correspondent aux positions sélectionnées de manière progressive dans la courbe de balayage enregistrée ; elles sont étiquetées « instantané 0 à 3 » plutôt que par les valeurs de tension exactes, car la position du curseur à elle seule ne détermine pas la valeur de polarisation. La tension appliquée sélectionnée est affichée dans le visualiseur et doit être lue directement pour une interprétation quantitative.
6.1 Bandes de conduction et de valence ; niveaux quasi-Fermi
Ouvrez les fichiers Ec.csv, Ev.csv, Fn.csv et Fp.csv à partir des données enregistrées. Les limites de la bande de conduction et de la bande de valence décrivent la barrière de jonction électrostatique. Les quasi-niveaux de Fermi des électrons et des trous indiquent comment les populations de porteurs s’écartent de l’équilibre. En l’absence de polarisation, le niveau de Fermi à l’équilibre doit être le même dans tout le dispositif connecté, tandis qu’une tension directe appliquée crée une séparation croissante entre les potentiels électrochimiques des porteurs.
Dans une hétérojonction GaN, les deux côtés utilisent le même semi-conducteur, de sorte que la variation de la bande de conduction à la jonction PN est principalement due à la courbure électrostatique de la bande. Comparez la barrière près de l’interface p-n à faible polarisation avec le profil obtenu après une augmentation de l’injection directe. La transition donneur n/n⁺ peut affecter la courbure locale sans modifier la bande interdite fondamentale du GaN.
6.2 Densités de porteurs d’électrons et de trous.
Les profils de densité des électrons libres et des trous libres enregistrés (Q_nfree.csv et Q_pfree.csv) donnent une représentation directe des porteurs majoritaires et minoritaires. Dans la première image, les électrons prédominent sur le côté de type n et les trous sur le côté de type p, avec des variations importantes près de la zone de déplétion. Sous polarisation directe, l’injection de porteurs minoritaires augmente la concentration de porteurs des deux côtés de la jonction, ce qui finit par modifier la forme spatiale des profils de densité.
L’échelle logarithmique de la densité est particulièrement utile pour comparer les porteurs minoritaires et majoritaires, qui peuvent différer de plusieurs ordres de grandeur. Il ne faut pas assimiler la concentration nominale des accepteurs de magnésium à la concentration des trous libres sans tenir compte de l’ionisation des accepteurs ; l’activation incomplète du magnésium est une caractéristique importante du GaN de type p obtenu expérimentalement.
6.3 Densités de courant des électrons et des trous.
Tracez les courbes Jn.csv et Jp.csv afin d’examiner la répartition du courant entre les électrons et les trous. À l’état stationnaire, les courants d’électrons et de trous peuvent échanger de l’énergie le long du dispositif lorsque les porteurs se recombinent. En tenant compte des conventions de signe, leur courant électrique total devrait être approximativement indépendant de la position en l’absence d’accumulation nette de charges.
La capture d’écran obtenue avec un faible courant présente des caractéristiques numériques bien définies. Une telle structure peut apparaître lorsque de faibles composantes de courant sont évaluées à partir de flux de dérive et de diffusion opposés, ou lorsqu’une transition abrupte n’est pas suffisamment prise en compte. Il ne faut pas l’interpréter automatiquement comme une pointe de courant physique. Il est nécessaire d’évaluer la convergence du solveur et de comparer le courant total traversant le dispositif avant de formuler une affirmation quantitative.
6.4 Mécanismes de recombinaison
Les captures d’écran de la recombinaison combinent les taux de recombinaison libre-libre et ceux liés au mécanisme SRH (R_nfree_to_pfree.csv et R_ss_srh.csv dans la séquence de captures d’écran enregistrées). À faible injection, la recombinaison médiatisée par les défauts peut être concentrée autour de la jonction, où les deux types de porteurs sont présents ; à une polarisation directe plus forte, la distribution spatiale de la recombinaison s’élargit à mesure que la population de porteurs minoritaires augmente. L’étendue et la taille relative de ces canaux dépendent des paramètres numériques du matériau et ne doivent pas être déduites uniquement de la hauteur d’un pic.
En comparant les profils de recombinaison obtenus pour différentes tensions de polarisation directe, on peut déterminer les zones où les électrons et les trous injectés interagissent le plus fortement. Ces informations spatiales complètent la courbe J–V et permettent de mettre en relation les variations de courant avec les populations de porteurs à l’intérieur de la jonction.
6.5 Profils de la vitesse des électrons et des trous.
Les valeurs de sortie des taux d’électrons et de trous (R_n.csv et R_p.csv) permettent d’affiner le diagnostic de l’équilibre local des porteurs de charge. Leurs signes et leurs composantes individuelles doivent être interprétés conformément aux conventions de sortie d’OghmaNano. Comparez les pics marqués observés à faible polarisation avec les distributions plus larges obtenues à une polarisation plus forte ; l’interprétation physique doit se concentrer sur les tendances cohérentes, et non sur l’interprétation de chaque variation numérique étroite comme un mécanisme distinct.
7. Exercices et interprétation.
8. Hypothèses utilisées pour la modélisation et analyses complémentaires.
La simulation utilise un modèle unidimensionnel et à l’état stationnaire de dérive-diffusion d’une structure semi-conductrice GaN compacte de 700 nm, avec des conditions aux limites ohmiques au niveau des contacts. Cette approche est particulièrement utile pour étudier l’influence de la concentration d’impuretés à la jonction, des paramètres des matériaux et de la recombinaison sur la structure de bande interne et la réponse courant-tension. Des études sur des configurations complètes de dispositifs de puissance peuvent s’appuyer sur ce modèle en intégrant la géométrie du dispositif, la résistance de contact et le comportement thermique.
La couche de p-GaN est caractérisée par une densité nominale d’accepteurs de magnésium. Étant donné que les accepteurs de magnésium dans le GaN sont relativement profonds, leur activation électrique et la densité de trous qui en résulte nécessitent une attention particulière. Les détails dépendent du modèle de matériau utilisé pour les statistiques des porteurs et l’ionisation des impuretés. De même, les défauts d’interface, les pièges dans le matériau, la mobilité des porteurs et les coefficients de recombinaison influencent la réponse électrique calculée.
La plage de balayage de 0 à 5 V permet d’observer en détail la conduction directe, depuis le début de l’injection de porteurs jusqu’à la plage de courant plus élevé. En étendant le balayage à la polarisation inverse, on obtient une autre approche utile pour analyser la zone de déplétion et le comportement du champ électrique. Les études plus avancées sur les dispositifs de puissance peuvent également inclure les effets de champ élevé et le blocage de la tension en fonction de la géométrie.
Pour obtenir des prédictions quantitatives de niveau recherche, vérifiez les propriétés évaluées du matériau GaN, examinez le résidu du solveur et la continuité du courant, affinez le maillage spatial et comparez les résultats avec les mesures expérimentales ou les données disponibles dans la littérature pour une structure présentant une correspondance étroite. Ces vérifications permettent d’améliorer le modèle initial en vue d’une comparaison quantitative avec les dispositifs fabriqués.
9. Prochaines étapes
Une fois l’analyse initiale de la courbe J–V terminée, l’étape suivante consiste à faire varier le type de dopage p et n et à comparer les variations de la courbure de la bande et le comportement de la zone de déplétion. Une étude complémentaire consiste à introduire une zone centrale faiblement dopée et à comparer les architectures GaN PN et PIN. Ces modifications permettent de démontrer comment la conception des couches affecte l’électrostatique, l’injection en polarisation directe et la capacité à supporter une tension inverse.