Simulation de la diode à hétérojonction AlGaAs/GaAs : alignement des bandes d’énergie, courbes J–V et transport des porteurs.
1. Introduction : pourquoi simuler une hétérojonction AlGaAs/GaAs ?
Une diode hétérojonction AlGaAs/GaAs est une diode semi-conductrice formée à l’endroit où deux matériaux différents de la famille III-V se rencontrent : l’arséniure de gallium et d’aluminium (AlxGa1−xAs) et l’arséniure de gallium (GaAs). Contrairement à une diode PN conventionnelle en GaAs, cet appareil combine une modification du dopage avec une modification de la structure de bande du semi-conducteur. La jonction présente donc à la fois un potentiel électrostatique intrinsèque et des décalages entre les bandes de conduction et de valence.
Dans ce tutoriel, nous étudions une jonction p-Al0,30Ga0,70As / n-GaAs à l’aide du solveur unidimensionnel de Poisson et de diffusion-dérive d’OghmaNano. La région AlGaAs, qui présente une bande interdite plus large, est en contact avec la région GaAs, qui présente une bande interdite plus étroite ; leur interface modifie l’injection d’électrons et de trous et constitue une introduction simple à l’alignement des bandes dans une hétérojonction. L’ensemble de la structure semi-conductrice ne fait que 700 nm d’épaisseur, ce qui facilite la visualisation des profils de transport internes.
Les hétérostructures AlGaAs/GaAs sont des éléments essentiels dans l’électronique et l’optoélectronique des matériaux du groupe III-V, notamment dans les dispositifs à hétérojonction, les couches de fenêtre photovoltaïques et les structures électroluminescentes. Cet exemple est spécifiquement une démonstration de diode/redresseur : il n’est pas présenté comme une photodiode, une diode laser ou un composant commercial à haute vitesse optimisés. Nous utiliserons le modèle pour établir un lien entre une caractéristique mesurable de la densité de courant d’obscurité en fonction de la tension (J–V) et le diagramme de bande microscopique ainsi que la physique de la recombinaison.
La vue principale est affichée dans ??. Le reste du tutoriel suit la même progression que les autres exemples de diodes OghmaNano : charger un dispositif, examiner ses entrées, exécuter une simulation et interpréter les résultats.
2. Créer une nouvelle simulation.
Dans la fenêtre principale d’OghmaNano, sélectionnez « Nouvelle simulation ». Le navigateur de démonstrations regroupe les exemples de semi-conducteurs connexes. Ouvrez les démonstrations GaAs/III–V et sélectionnez l’exemple de diode à hétérojonction AlGaAs/GaAs, ou chargez directement le projet d’exemple fourni si son libellé de menu diffère selon les versions. Les deux étapes du navigateur sont présentées dans ?? et ??.
Une fois l’exemple ouvert, examinez la pile verticale dans la vue de l’appareil. Bien que la visualisation soit en trois dimensions, ce tutoriel traite du transport électrique unidimensionnel à travers l’épaisseur de la couche. Le modèle ne prend pas en compte la concentration latérale du courant ni les géométries réelles des contacts métalliques.
3. Analyse de la structure en couches, de la composition de l’alliage et du dopage.
3.1 Empilement des couches
Ouvrez l’éditeur de calques à partir du ruban de l’appareil. La diode active est composée de trois couches semi-conductrices, avec une zone d’air au-dessus pour la représentation de l’appareil. Le problème électrique concerne l’ensemble de la pile semi-conductrice de 700 nm plutôt que la zone d’air de 200 nm.
| Région (de haut en bas) | Épaisseur | Utilisation de substances dopantes chez les mannequins. | Objectif |
|---|---|---|---|
| Air (non conducteur d’électricité) | 200 nm | — | Environnement situé au-dessus de l’appareil. |
| p-Al0,30Ga0,70As | 200 nm | NA = 5 × 1017 cm−3 | Face de type p à large bande interdite d’une hétérojonction. |
| n-GaAs | 400 nm | ND = 5 × 1017 cm−3 | Côté de type n à faible largeur de bande interdite. |
| n+-GaAs | 100 nanomètres | ND = 1 × 1018 cm−3 | Région du côté cathode avec une concentration de dopants plus élevée. |
La jonction p-AlGaAs/n-GaAs constitue à la fois la jonction métal-semiconducteur PN et l’interface entre les matériaux. La jonction n-GaAs/n+-GaAs représente une transition de dopage au sein du GaAs, et non une seconde hétérojonction. L’éditeur de couches (??) et la configuration des contacts (??) sont présentés ensemble ci-dessous.
3.2 Dopage et composition en aluminium
Ouvrez l’éditeur « Doping / Ions » à partir du ruban « Electrical ». Le côté AlGaAs dopé accepteur et le côté GaAs dopé donneur forment la jonction de redressement (??). Le type de dopage change à l’interface des matériaux, ce qui crée un champ électrique interne. La région n+ augmente la concentration d’électrons près de la cathode.
Ensuite, ouvrez l’éditeur de composition/de gradients d’alliage. L’alliage AlGaAs présente une fraction molaire d’aluminium de x = 0,30, tandis que GaAs correspond à x = 0. La variation abrupte de la composition (??) est l’une des principales raisons pour lesquelles le diagramme des bords de bande présente une discontinuité, en plus de la courbure électrostatique habituelle.
3.3 Pourquoi l’interface est-elle importante ?
À température ambiante, le GaAs présente une largeur de bande interdite d’environ 1,42 eV, tandis que l’Al0,30Ga0,70As présente une largeur de bande interdite plus importante, généralement d’environ 1,8 eV, en fonction des paramètres de l’alliage. Dans une hétérojonction, la différence totale de largeur de bande interdite est répartie entre les décalages de bande de conduction et de bande de valence :
\[\Delta E_g = \Delta E_c + \Delta E_v\]Ici, les symboles de décalage indiquent l’amplitude des étapes correspondantes. Leurs valeurs numériques sont déterminées par le modèle de matériau et les hypothèses concernant l’alignement des bandes, et non uniquement par le profil de dopage. Lors de l’analyse ultérieure des limites de bande simulées, il est important de distinguer un décalage net lié au matériau de la courbure douce de la bande produite par la charge d’espace.
4. Paramètres des matériaux électriques, interfaces et recombinaison.
4.1 La base de données de matériaux et les modèles Lua.
Ouvrez l’éditeur des paramètres électriques pour chaque couche. Dans cet exemple, Inorganic/AlGaAs/AlGaAs.lua est utilisé pour l’alliage de type p et Inorganic/GaAs/GaAs.lua pour les deux régions de type n en GaAs. Le modèle Lua AlGaAs (??) et ses paramètres (??) peuvent être comparés au modèle GaAs (??) et à ses paramètres (??).
La simulation est configurée de manière à utiliser les valeurs clés des matériaux – notamment la largeur de bande interdite, l’affinité électronique, la mobilité des porteurs, les densités d’états effectives et les paramètres de recombinaison – à partir de la base de données de matériaux Lua, plutôt que d’imposer des valeurs locales. Ceci est particulièrement important pour l’AlGaAs, car ses paramètres dépendent de la composition de l’alliage. Lisez les valeurs affichées à l’écran et renvoyées par le modèle, plutôt que d’interpréter les valeurs de secours inactives présentes dans le fichier du projet comme les propriétés réelles du matériau.
4.2 Alignement des bandes et paramètres d’interface
Ouvrez l’éditeur d’interfaces, illustré dans ??. Cet exemple représente une interface abrupte. Dans un modèle de base de transport par dérive et diffusion, l’injection de porteurs dépend des discontinuités des bandes de conduction et de valence, des statistiques des porteurs et du traitement de la limite de l’interface. Les échantillons réels peuvent également présenter des états d’interface, un désordre d’alliage, une émission thermoïonique, un effet tunnel ou des imperfections induites par le processus de fabrication ; il ne faut pas supposer que ces effets sont pris en compte simplement parce qu’une hétérojonction est présente.
4.3 Recombinaison radiative et assistée par des défauts.
À cette composition, le GaAs et l’Al0,30Ga0,70As sont des matériaux à bande interdite directe, ce qui signifie que la recombinaison radiative peut se produire lorsque des porteurs de charge sont injectés. Un modèle courant est le suivant :
\[R_{\mathrm{rad}} = B\left(np-n_i^2\right)\]où B est un coefficient qui dépend du matériau, et n et p représentent les densités locales de porteurs libres. La configuration d’entrée permet également une recombinaison de type Shockley–Read–Hall (SRH) à un seul niveau dans les trois régions actives du semi-conducteur. Une expression idéalisée est la suivante :
\[R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}\]La méthode SRH décrit la recombinaison assistée par des pièges et est sensible aux paramètres des pièges définis par les modèles choisis. Dans cet exemple, la recombinaison d’Auger est désactivée. Il ne faut pas conclure qu’un seul processus de recombinaison est dominant sans comparer son ampleur à celle des autres processus activés dans les mêmes conditions de polarisation.
L’anode et la cathode sont soumises à des conditions aux limites ohmiques idéalisées. Ces conditions permettent d’isoler les phénomènes physiques qui se produisent au niveau de l’hétérojonction, mais elles ne permettent pas de prédire le comportement d’un schéma de métallisation particulier ni sa résistance de contact.
5. Exécution de la simulation J–V en mode sombre.
Retournez à la fenêtre principale de simulation et cliquez sur « Lancer la simulation ». Le fichier d’entrée spécifie une variation de tension allant de −0,1 V à +2,5 V, par incréments de 0,025 V, sans illumination optique intentionnelle. OghmaNano résout les équations couplées de Poisson et de continuité des électrons/trous à chaque point de polarisation successif.
Pendant le calcul, passez à l’onglet « Terminal » pour suivre l’avancement du solveur, puis ouvrez la fenêtre « Sortie » une fois le calcul terminé. Le navigateur de sortie (??) et le terminal (??) illustrent ces différentes étapes. Les captures d’écran du terminal sont prises à différents moments, et non à des moments fixes correspondant à des valeurs de polarisation électrique spécifiques.
Double-cliquez sur jv.csv pour ouvrir la courbe simulée de la densité de courant en fonction de la tension ( ?? ). La caractéristique de redressement doit présenter un courant relativement faible sous une faible polarisation inverse et un courant direct de plus en plus important à mesure que l’injection augmente. Pour cette hétérojonction particulière, le début et la forme détaillée dépendent des décalages de bande, des paramètres des matériaux et des taux de recombinaison, et non pas seulement de la bande interdite du GaAs.
L’expression de la diode idéale peut s’avérer utile comme outil d’interprétation locale :
\[J \approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right]\]Cependant, un facteur d’idéalité ajusté, n, ou un courant de saturation, J0, ne devraient être déterminés qu’à partir d’un intervalle de tension approprié. La résistance série, l’injection élevée et le transport à l’hétérojonction peuvent entraîner des écarts par rapport à une simple loi exponentielle. La capture d’écran fournie est à titre d’illustration et ne constitue pas une validation expérimentale indépendante.
6. Analyse des captures d’écran de la simulation : que se passe-t-il à l’intérieur de la diode ?
OghmaNano enregistre les résultats électriques résolus spatialement pendant le balayage J–V. Ouvrez le visualiseur d’instantanés et utilisez son curseur pour comparer les positions de polarisation sélectionnées. Les quatre images de chaque ligne ci-dessous correspondent respectivement au début, à environ un tiers, à environ deux tiers et à la fin du balayage enregistré. Les étiquettes 0 à 3 indiquent les positions des captures d’écran, et non les tensions en volts ; utilisez l’indicateur de polarisation du visualiseur pour obtenir les valeurs exactes pour une exécution spécifique.
6.1 Bande de conduction, bande de valence et niveaux quasi-Fermi.
Superposez les courbes Ec.csv, Ev.csv, Fn.csv et Fp.csv. La première capture d’écran du diagramme de bandes (??) montre le bord de la bande de conduction Ec et le bord de la bande de valence Ev, révélant à la fois le potentiel électrostatique continu et la discontinuité abrupte de la bande à l’interface entre les matériaux AlGaAs et GaAs. Les quasi-niveaux de Fermi des électrons et des trous révèlent les écarts par rapport à l’équilibre et la distribution des porteurs injectés.
Comparez la première et la dernière images. En polarisation directe, la barrière de la jonction se modifie et les niveaux quasi-Fermi se séparent. Ne confondez pas un palier au niveau du bord de bande à l’interface avec du bruit numérique : une véritable hétérojonction peut présenter des discontinuités au niveau des bords de bande, même si le potentiel électrostatique reste continu.
6.2 Densités d’électrons et de trous.
Dans le cadre de ce projet, représentez graphiquement les valeurs de concentration des électrons libres et des trous libres, Q_nfree.csv et Q_pfree.csv. Le côté p-AlGaAs est constitué d’un matériau à prédominance de trous, tandis que le côté GaAs est constitué d’un matériau à prédominance d’électrons. L’augmentation de la polarisation directe injecte des porteurs minoritaires à travers la jonction et modifie la distribution spatiale des porteurs. Comme le montrent les résultats relatifs à la densité des porteurs (??), une limite matérielle abrupte modifie les états électroniques disponibles et l’alignement des bandes, de sorte que le profil de densité des porteurs n’a pas besoin de ressembler à celui d’une jonction PN homogène en GaAs.
6.3 Courants d’électrons et de trous.
Tracez les courbes Jn.csv et Jp.csv afin de séparer les contributions des électrons et des trous au courant total. À l’état stationnaire, les deux composantes peuvent se redistribuer en fonction de la position, les électrons et les trous se recombinant, tandis que la somme des densités de courant électrique, en tenant compte de leur signe, doit rester approximativement constante en l’absence de sources de courant supplémentaires. Les pics locaux ou les variations rapides observées à proximité de l’hétérojonction sur les graphiques du courant porteur (??) doivent être vérifiés par rapport au maillage et au diagramme des bandes.
Les simulations de faibles courants peuvent être délicates sur le plan numérique. Il est préférable d’évaluer la tendance générale et la continuité du courant plutôt que d’attribuer chaque oscillation mineure à un mécanisme de transport physique. Si vous constatez des écarts importants, affinez le maillage et vérifiez la convergence du solveur avant de tirer des conclusions.
6.4 Recombinaison radiative et recombinaison de Shockley-Read-Hall
Analysez R_nfree_to_pfree.csv (recombinaison radiative libre-libre) et R_ss_srh.csv (recombinaison SRH à un seul niveau). Les captures d’écran de la recombinaison (??) affichent ces canaux ensemble. En cas d’injection, le taux de recombinaison absolu peut augmenter considérablement ; son emplacement est déterminé par le chevauchement des porteurs libres, les propriétés du matériau et les paramètres de piège choisis.
Portez une attention particulière aux caractéristiques présentes à l’interface p-AlGaAs/n-GaAs. Une caractéristique locale peut avoir une signification physique, mais sa simple amplitude ne suffit pas à prouver l’existence de pièges à l’interface : cet exercice utilise des modèles de matériaux en volume et les paramètres SRH configurés. Comparez la contribution intégrée et la distribution spatiale des deux processus avant de conclure à un mécanisme de perte dominant.
6.5 Profils de la vitesse des électrons et des trous.
Les résultats distincts de R_n.csv et R_p.csv mettent en évidence les équilibres de taux liés aux électrons et aux trous utilisés dans la simulation. Ces graphiques (??) permettent de vérifier plus précisément où les porteurs sont éliminés ou redistribués pendant le balayage. Interprétez les signes et les éléments représentés sur les graphiques en fonction de leurs en-têtes de données, plutôt que de supposer que toutes les courbes représentent un taux de recombinaison positif unique.
6.6 Mobilité des porteurs de charge
La fonction de mobilité permet d’examiner les paramètres de transport en fonction de la position. On s’attend à ce que des variations se produisent à l’interface, car les modèles AlGaAs et GaAs ne sont pas tenus de renvoyer la même valeur de mobilité. La capture d’écran de l’exemple de mobilité (??) montre une position d’affichage enregistrée ; elle ne constitue pas une preuve d’un effet particulier de la mobilité en fonction du champ.
7. Exercices : établir un lien entre les décalages de bande et le comportement des diodes.
8. Hypothèses et limites de l’exemple.
Il s’agit d’un exemple didactique simplifié en une dimension. Les épaisseurs du semi-conducteur et les concentrations de dopage sont plausibles pour une structure épitaxiale III-V, mais cet exemple, pris seul, ne permet pas de définir un produit réalisable à l’échelle industrielle ou optimisé en termes de performances. Le calcul suppose des contacts ohmiques idéalisés, une structure uniforme en termes de dimensions latérales, une interface nette entre les matériaux, ainsi que les phénomènes de transport et de recombinaison décrits par les modèles sélectionnés.
L’interface AlGaAs/GaAs peut impliquer des phénomènes de transport supplémentaires dans les dispositifs réels. Les contacts pratiques, les défauts, la recombinaison en surface, le désordre de l’alliage, l’auto-échauffement et les effets du substrat peuvent également modifier les caractéristiques mesurées. Dans cette configuration, la recombinaison d’Auger est désactivée, et la simulation est une mesure électrique effectuée dans l’obscurité, et non une prédiction de l’efficacité d’une LED, d’une photodiode ou d’une cellule solaire calibrée. Une validation quantitative nécessite des données expérimentales et des paramètres de matériaux vérifiés de manière indépendante.
9. Prochaines étapes
Après avoir reproduit la courbe J–V et localisé l’hétérojonction dans le diagramme de bandes, le calcul le plus pertinent à effectuer est une comparaison contrôlée de la composition de l’alliage : il faut faire varier la proportion d’aluminium, maintenir le dopage et l’épaisseur constants, et suivre l’évolution du saut de la bande de conduction, du saut de la bande de valence et du courant en polarisation directe. Une simulation d’une homojonction GaAs fournit une référence utile.
Pour obtenir de plus amples informations sur la manière dont OghmaNano enregistre la configuration des appareils et les paramètres des modèles de matériaux, consultez la documentation sur le format de fichier Oghma avancé.