Simulation d’une diode à jonction PN en GaAs : modèles de dérive-diffusion, courbes caractéristiques courant-tension et recombinaison.

2. Création d’une nouvelle simulation.

Pour commencer, créez une nouvelle simulation à partir de la fenêtre principale d’OghmaNano. Cliquez sur le bouton « Nouvelle simulation » dans la barre d’outils. Cela ouvre la boîte de dialogue de sélection du type de simulation (voir ??).

Dans le dialogue de type simulation, double-cliquez sur les exemples GaAs, puis sélectionnez l’exemple de jonction/diode GaAs (voir ??). OghmaNano chargera une structure de jonction GaAs prédéfinie que nous traiterons comme une diode PN.

La structure du dispositif, une fois alimenté, est affichée dans la fenêtre principale de simulation (voir ??). Bien que le problème électrique résolu dans ce tutoriel soit unidimensionnel, la vue en 3D permet de visualiser clairement l’empilement des couches verticales et les zones qui participent au transport et à la recombinaison des porteurs de charge.

La diode est constituée d’une séquence de couches de GaAs empilées verticalement, comprenant une région p+ fortement dopée, une région p moins fortement dopée, une région n faiblement dopée et une région n+ fortement dopée. Cette structure est explicitement répertoriée dans l’éditeur de couches (voir ??), où chaque couche se voit attribuer une épaisseur, un matériau et un rôle électrique.

Les couches centrales de type p et de type n forment la jonction PN active. À l’équilibre, une zone de déplétion se forme à l’interface, ce qui donne lieu à un champ électrique intrinsèque qui contrôle la séparation et le transport des porteurs de charge. Les couches minces et fortement dopées de type p+ et de type n+ servent de zones de contact à faible résistance, ce qui garantit que la tension appliquée se concentre principalement sur la jonction plutôt qu’aux contacts.

Dans les sections suivantes, cette structure sera traitée comme un dispositif unidimensionnel : toutes les variations seront prises en compte dans le sens de la croissance, tandis que les variations latérales seront négligées. Malgré cette simplification, le modèle prend en compte les éléments essentiels de l’électrostatique, du transport des porteurs et de la physique de la recombinaison qui régissent le comportement en mode obscur des diodes PN en GaAs, utilisées dans des dispositifs optoélectroniques et électroniques à haute vitesse.

3. Analyse du profil de dopage.

Le profil de dopage définit la jonction PN du GaAs et, par conséquent, détermine les caractéristiques électrostatiques fondamentales de la diode. Il détermine la position de la jonction, le potentiel intégré, la largeur de la zone de déplétion et le champ électrique interne qui se développe à l’équilibre et sous polarisation.

Pour afficher la configuration de dopage, ouvrez l’éditeur « Dopage / Ions » dans le ruban « Électrique » (voir ??). L’éditeur affiche la distribution spatiale des donneurs et des accepteurs ionisés en fonction de la profondeur (voir ??).

Dans ce tutoriel, la diode est fabriquée en utilisant un profil de dopage GaAs p+/p/n/n+ classique. Les régions centrales p et n sont modérément dopées et forment la jonction PN active, où se développent la zone de déplétion et le champ électrique intrinsèque.

Les couches minces, fortement dopées de type p+ et n+, servent de zones de contact à faible résistance. Leur rôle est de favoriser une bonne injection et extraction des porteurs de charge, tout en veillant à ce que la majeure partie de la tension appliquée se dissipe au niveau de la jonction elle-même plutôt qu’aux contacts.

Dans le cadre de ce tutoriel, l’élément essentiel à vérifier est que le dispositif contienne une région dopée principalement de type accepteur et une région dopée principalement de type donneur, avec une transition nette entre les deux. Les valeurs numériques exactes des concentrations de dopage affectent principalement la largeur de la zone de déplétion et l’intensité du champ intégré, ce qui sera étudié indirectement à travers les caractéristiques courant-tension en mode bloqué de la diode dans les sections suivantes.

4. Analyse des paramètres électriques et des mécanismes de recombinaison.

Les paramètres des matériaux électriques sont définis pour chaque région et contrôlent le transport des porteurs de charge, la recombinaison et l’électrostatique dans toute la diode. Ouvrez l’éditeur de paramètres électriques à partir de la fenêtre principale en sélectionnant Structure du dispositif → Paramètres électriques. Chaque couche de la structure du dispositif possède son propre onglet de paramètres. Dans ce tutoriel, nous utilisons le même modèle de matériau GaAs dans les quatre régions, mais nous les interprétons différemment : les régions p+ et n+ servent de régions de contact à faible résistance, tandis que les régions p et n forment la jonction PN active.

Les figures ?? à ?? présentent l’éditeur des paramètres électriques pour chaque région (p+, p, n, n+). Pour l’arséniure de gallium (GaAs), les paramètres les plus remarquables sont les mobilités des porteurs : la mobilité des électrons est élevée par rapport au silicium, tandis que la mobilité des trous est plus modeste. Cette asymétrie est une caractéristique essentielle du GaAs et explique son utilisation généralisée dans l’électronique à haute fréquence.

Étant donné que le transport des porteurs dans le GaAs est généralement rapide, en particulier pour les électrons, le comportement en courant continu (I–V) d’une jonction PN en GaAs n’est souvent pas limité par le transport dans le volume du matériau. Au contraire, l’équilibre entre l’injection, l’électrostatique et la recombinaison détermine la réponse du dispositif. Les paramètres présentés dans l’éditeur contrôlent donc principalement l’efficacité avec laquelle les porteurs injectés sont éliminés par les processus de recombinaison une fois que la barrière de la jonction est réduite.

Les densités d’états effectives affichées dans l’éditeur (par exemple, ??) déterminent les statistiques des porteurs de charge et les concentrations d’équilibre des porteurs de charge pour le GaAs. Ces valeurs diffèrent de celles du silicium en raison de la structure de bande et des masses effectives différentes du GaAs, et elles influencent directement les niveaux d’injection des porteurs de charge et les taux de recombinaison par le biais des densités d’électrons et de trous \(n\) et \(p\).

Recombinaison radiative libre-libre

Étant donné que le GaAs est un semi-conducteur à bande interdite directe, la recombinaison entre les électrons libres et les trous libres peut se produire efficacement par le biais de transitions radiatives. Dans le GaAs de haute qualité, ce canal de recombinaison entre électrons et trous libres est souvent le principal mécanisme de perte dans la jonction active sous polarisation directe. Chez OghmaNano, ce processus est contrôlé par la constante de vitesse de recombinaison entre les électrons et les trous libres, qui est visible dans l’éditeur des paramètres électriques (voir, par exemple, ??).

Le taux de recombinaison correspondant est de la forme suivante :

où \(B\) est le coefficient de recombinaison radiative libre-libre. Dans le GaAs, ce terme peut dominer l’équilibre de la recombinaison une fois que les porteurs de charge sont injectés dans la jonction, en particulier lorsque la densité des défauts est faible. Par conséquent, la caractéristique courant-tension dans l’obscurité d’une diode GaAs de haute qualité reflète souvent la physique de la recombinaison radiative plutôt qu’un comportement limité par les défauts.

Recombinaison de Shockley-Read-Hall (SRH)

La recombinaison de Shockley–Read–Hall (SRH) prend en compte la recombinaison médiatisée par les défauts à travers les états électroniques dans la bande interdite. Dans OghmaNano, ce processus est contrôlé à l’aide des paramètres d’équilibre des pièges SRH, qui sont affichés dans l’éditeur (voir ??) : une énergie de piège \(E_t\) (par rapport au milieu de la bande interdite), une densité de piège \(N_t\), et les sections efficaces de capture des électrons et des trous \(\sigma_n\) et \(\sigma_p\).

Dans le GaAs, la recombinaison de Shockley-Read-Hall n’est généralement pas considérée comme la limite intrinsèque de recombinaison du matériau. Elle permet plutôt d’évaluer la qualité du matériau et la densité des défauts. Dans le GaAs épitaxial de haute qualité, la recombinaison de Shockley-Read-Hall est faible et la recombinaison radiative est prédominante ; dans les matériaux de qualité inférieure ou près des interfaces et des zones endommagées par le processus de fabrication, la recombinaison de Shockley-Read-Hall peut devenir significative.

La durée de vie des SRH est définie par

et le taux de recombinaison qui en résulte est de

Dans ce tutoriel, la recombinaison SRH est intentionnellement maintenue. Cela vous permet d’étudier comment l’augmentation de la densité de défauts ou la réduction de la durée de vie des porteurs modifie le comportement de la diode, la faisant passer d’un fonctionnement dominé par l’émission radiative à un fonctionnement limité par la recombinaison, ce qui est particulièrement pertinent pour comprendre le fonctionnement des dispositifs réels, où les processus de fabrication et les interfaces jouent un rôle essentiel.

Électrostatique et paramètres de bande.

Enfin, la structure de bande et les paramètres électrostatiques utilisés pour définir le GaAs sont visibles dans chaque onglet de la région (par exemple, ??) : l’affinité électronique, la largeur de bande (\(E_g \approx 1.42\,\mathrm{eV}\) à température ambiante) et la permittivité relative (\(\varepsilon_r \approx 12.9\)). Ces paramètres définissent le potentiel intégré de la jonction PN et déterminent l’échelle de tension sur laquelle se produit une injection significative de porteurs.

Dans le cas du GaAs, la combinaison d’une bande interdite directe et d’une mobilité élevée des porteurs signifie qu’une fois la barrière de jonction réduite, les porteurs sont injectés efficacement et la physique de la recombinaison devient le facteur dominant qui façonne la courbe courant-tension en l’absence de lumière. Il est donc essentiel de comprendre comment ces paramètres électriques interagissent afin d’interpréter à la fois les résultats de la simulation et le comportement réel des dispositifs en GaAs.

5. Exécution de la simulation, obtention de courbes caractéristiques courant-tension (I-V) en mode sombre, et extraction des paramètres.

Une fois que la structure du dispositif, le profil de dopage et les paramètres électriques sont définis, la simulation de la diode peut être lancée directement à partir de la fenêtre principale. Cliquez sur « Lancer la simulation » pour démarrer le solveur. Pendant l’exécution, les informations de convergence pour chaque point de polarisation sont affichées dans la fenêtre de sortie, ce qui vous permet de surveiller la stabilité et l’avancement du solveur (voir ??).

Pour vérifier les caractéristiques de la diode, ouvrez l’onglet « Sortie » et double-cliquez sur jv.csv (voir ??). Pour une diode au GaAs correctement configurée, la courbe I-V doit être régulière et monotone. En polarisation inverse, le courant reste faible et dépend faiblement de la tension, ce qui indique une saturation limitée par la recombinaison. En polarisation directe, le courant augmente rapidement avec la tension appliquée, ce qui correspond à l’injection de porteurs à travers la jonction PN.

La forme de la zone de polarisation directe contient des informations physiques utiles. Sur un graphique semi-logarithmique, la pente de la zone exponentielle peut être utilisée pour déterminer un facteur d’idéalité, ce qui permet de savoir si le courant est principalement limité par la diffusion (\(n \approx 1\)) ou par des processus limités par la recombinaison (\(n \approx 2\)). L’ordonnée à l’origine extrapolée de cette zone fournit une estimation du courant de saturation inverse, qui est directement liée aux paramètres de recombinaison SRH et Auger mentionnés dans la section 4. Dans le cas du GaAs, la bande interdite directe et les mobilités généralement élevées signifient que l’injection peut être très efficace, de sorte que les paramètres de recombinaison présentent souvent une signature particulièrement nette dans la pente et l’ordonnée à l’origine de la zone de polarisation directe.

En règle générale, il est toujours conseillé d’examiner la courbe I–V avant d’interpréter les valeurs dérivées. Les discontinuités, les conventions de signe inattendues ou les variations de courant non physiques indiquent généralement des problèmes liés aux conditions aux limites, aux étapes de polarisation, aux paramètres de recombinaison ou à la convergence du solveur. Pour une simple diode PN au GaAs comme celle-ci, la courbe I–V dans l’obscurité doit être facilement compréhensible et intuitive.

6. Analyse des captures d’écran de la simulation : bandes, recombinaison et flux de courant.

Lors d’une mesure de la courbe I–V, OghmaNano enregistre la solution interne des équations de dérive-diffusion à chaque point de polarisation dans le répertoire des instantanés. Ces fichiers permettent de visualiser les prédictions du solveur concernant ce qui se passe à l’intérieur de la diode : courbure des bandes d’énergie, séparation du niveau quasi-Fermi, activité de recombinaison et transport du courant. L’analyse de ces paramètres est essentielle pour comprendre pourquoi une courbe I–V particulière se manifeste.

Dans cette section, nous examinons trois points représentatifs de polarisation : une polarisation inverse proche de l’équilibre (−0,1 V), une polarisation directe modérée proche du seuil de fonctionnement (≈ 0,45 V) et une polarisation directe élevée (0,8 V). Ensemble, ces exemples illustrent la transition de l’équilibre, en passant par un transport limité par l’injection, jusqu’à un fonctionnement à injection élevée. Pour le GaAs, la même transition qualitative s’applique, mais le principal canal de recombinaison dans la région active est généralement une recombinaison radiative plutôt qu’une recombinaison limitée par les défauts (SRH).

6.1 Bords de bande et niveaux quasi-Fermi

Pour reproduire les diagrammes de bandes, ouvrez le visualiseur d’instantanés et ajoutez les fichiers Ec.csv, Ev.csv, Fn.csv et Fp.csv. Ils correspondent respectivement au bord de la bande de conduction, au bord de la bande de valence, au niveau quasi-Fermi des électrons et au niveau quasi-Fermi des trous.

À −0,1 V (figure ??), la diode est presque à l’équilibre. La courbure des bandes reflète le potentiel intégré imposé par le profil de dopage, et les niveaux quasi-Fermi sont presque plats et coïncident, ce qui indique un faible flux de courant net. La zone de déplétion est clairement visible en tant que zone de forte courbure des bandes au niveau de la jonction. À environ 0,45 V (figure ??), la polarisation directe réduit la barrière de la jonction. Les niveaux quasi-Fermi des électrons et des trous se séparent à travers la zone de déplétion, ce qui constitue la signature interne de l’injection de porteurs. Cette séparation des niveaux quasi-Fermi est directement responsable de l’augmentation exponentielle du courant observée dans la courbe I–V. À 0,8 V (figure ??), la jonction est fortement polarisée en direct. La barrière est fortement réduite, les niveaux quasi-Fermi sont largement séparés et le dispositif fonctionne dans un régime d’injection élevée où les densités de porteurs sont importantes dans une grande partie de la structure. Dans le GaAs, cette séparation des niveaux quasi-Fermi est particulièrement révélatrice car elle est directement liée à une forte injection de porteurs et à une recombinaison efficace des porteurs libres dans la région active.

6.4 Recombinaison radiative libre-libre dans la région active.

La recombinaison libre-libre peut être étudiée en traçant R_nfree_to_pfree.csv, qui représente le taux de recombinaison radiative local en fonction de la position. Dans le GaAs, ce processus est généralement le principal mécanisme de recombinaison dans la jonction active sous polarisation directe, car le GaAs est un semi-conducteur à bande interdite directe et les transitions radiatives sont efficaces. Lors de l’interprétation de ces graphiques, l’élément le plus important n’est pas la structure spatiale fine, mais plutôt la forme générale du profil de recombinaison et la manière dont son amplitude évolue en fonction de la polarisation.

À très faible polarisation (figure ??), la recombinaison libre-libre est faible, car le dispositif est proche de l’équilibre et le produit de concentration des porteurs \(np\) reste proche de \(n_i^2\). Bien que l’arséniure de gallium (GaAs) permette des transitions radiatives efficaces, les concentrations de porteurs nécessaires pour générer un taux de transition radiative élevé ne sont pas encore atteintes.

Lorsque la diode est soumise à une polarisation directe (figure ??), les électrons et les trous injectés coexistent dans la majeure partie de la structure. Le profil de recombinaison résultant, où les porteurs se recombinent librement, prend une forme de « calotte » large et légèrement incurvée qui s’étend sur presque tout le dispositif. Cette forme reflète le fait que la recombinaison radiative dépend principalement du produit local des porteurs \(np\), qui est relativement uniforme dans les régions quasi-neutres une fois l’injection établie. Contrairement à la recombinaison SRH, il n’y a pas de localisation marquée dans la zone de déplétion.

À forte polarisation directe (figure ??), la forme du profil de recombinaison libre-libre reste globalement similaire, tandis que son amplitude augmente considérablement. Il s’agit d’une caractéristique diagnostique essentielle de la recombinaison à prédominance radiative dans le GaAs : l’augmentation de la polarisation entraîne principalement une augmentation de la densité des porteurs dans tout le dispositif, plutôt qu’une modification de la zone où se produit la recombinaison. Le profil uniforme qui s’étend sur tout le dispositif indique que la recombinaison radiative agit comme un processus volumique dans toutes les régions où des porteurs sont injectés.

Ce comportement contraste fortement avec la recombinaison médiée par les défauts. La recombinaison de type SRH produit généralement des pics étroits, centrés sur la jonction, qui ne s’élargissent qu’à des niveaux d’injection élevés, tandis que la recombinaison libre-libre dans le GaAs produit un seul « dôme » de recombinaison étendu, dont la hauteur – et non son étendue spatiale – indique le niveau d’injection. Dans les dispositifs à bande interdite directe, ce profil radiatif large est donc un indicateur de la haute qualité du matériau et de l’efficacité de l’injection des porteurs, et il contribue souvent de manière prépondérante à l’équilibre de recombinaison qui détermine la caractéristique I–V.

6.2 Recombinaison de Shockley-Read-Hall

Afin d’étudier la recombinaison médiée par les défauts, représentez R_ss_srh.csv, qui indique le taux de recombinaison de Shockley-Read-Hall résolu spatialement à l’intérieur de la diode. Les trois graphiques ci-dessous correspondent aux mêmes points de polarisation utilisés dans l’analyse du diagramme de bandes : −0,1 V, ≈0,45 V et 0,8 V. Pour le GaAs, ces graphiques s’interprètent le mieux comme un indicateur des pertes assistées par des pièges plutôt que comme le principal canal de recombinaison dans une jonction de haute qualité. Le point essentiel à considérer est la manière dont la localisation spatiale de la recombinaison SRH se compare au profil plus large de recombinaison libre-libre présenté ci-dessus.

À −0,1 V (figure ??), la diode est proche de l’équilibre. Les électrons prédominent dans la partie de type n et les trous dans la partie de type p, de sorte qu’une recombinaison SRH importante ne peut se produire que dans la zone étroite autour de la jonction, où les deux types de porteurs sont présents simultanément. Par conséquent, le taux de recombinaison SRH est fortement localisé au centre du dispositif, ce qui correspond à la zone de déplétion. À environ 0,45 V (figure ??), la polarisation directe injecte des porteurs à travers la jonction et augmente le produit local des densités d’électrons et de trous. Le pic SRH augmente considérablement en amplitude, mais il reste spatialement confiné à la région centrale du dispositif. Cela indique que la recombinaison SRH reste un canal de perte centré sur la jonction, contrôlé par le chevauchement des porteurs et l’activité des pièges dans la zone de déplétion et à proximité de celle-ci. À 0,8 V (figure ??), le comportement s’étend à mesure que l’injection de porteurs augmente dans toute la structure. Cependant, même lorsque la recombinaison SRH se propage, elle reste généralement plus fortement concentrée dans la région de la jonction que la recombinaison libre-libre, qui tend à se répartir dans les régions quasi-neutres injectées dans le GaAs.

La comparaison des deux mécanismes de recombinaison permet d’obtenir une représentation interne utile du fonctionnement des diodes au GaAs. La recombinaison directe crée un puits volumique une fois l’injection établie, et elle domine généralement l’équilibre global de la recombinaison dans le volume actif du dispositif. En revanche, la recombinaison de type SRH met en évidence les zones où les pertes assistées par des pièges sont concentrées, souvent dans la zone de déplétion ou à proximité, ce qui en fait un indicateur sensible de la qualité du matériau et des dommages à l’interface. Cette distinction permet d’expliquer pourquoi les diodes au GaAs peuvent présenter une forte injection et des courants importants tout en restant très sensibles aux défauts au niveau de la jonction.

6.3 Densités de courant des électrons et des trous.

Enfin, examinez les courants porteurs en traçant les graphiques de Jn.csv et Jp.csv, qui représentent respectivement les densités de courant des électrons et des trous en fonction de la position. Ces graphiques permettent de visualiser directement la manière dont la charge est transportée à travers la diode dans différentes conditions de polarisation et comment le dispositif passe de l’état d’équilibre à un état de conduction continue en polarisation directe.

À -0,1 V (figure ??), la diode est proche de l’équilibre et le courant physique réel est extrêmement faible. Les flux d’électrons et de trous sont presque équilibrés dans tout le dispositif, de sorte que le courant net résulte de la différence entre deux quantités presque égales. Dans ce régime, le problème numérique est intrinsèquement mal conditionné, et de petites oscillations ou un bruit apparent dans les profils de courant sont à prévoir. Dans cette démonstration avec du GaAs, l’effet peut sembler plus prononcé : la combinaison d’un transport de porteurs très rapide (en particulier pour les électrons) et d’un courant net extrêmement faible rend le problème d’annulation plus marqué, de sorte que les tracés bruts Jn/Jp peuvent sembler visiblement bruités, même si la physique sous-jacente est simplement celle d’un courant « presque nul ». À environ 0,45 V (figure ??), la polarisation directe entraîne l’injection de porteurs à travers la jonction. Le courant d’électrons domine sur le côté n et le courant de trous domine sur le côté p, mais les deux courants sont continus dans tout le dispositif, ce qui reflète la conservation de la charge à l’état stationnaire. La densité de courant augmente rapidement par rapport au cas proche de l’équilibre, mais les profils peuvent encore présenter une petite structure numérique près de la jonction, où les gradients sont importants. Dans le GaAs, où l’injection peut s’établir rapidement une fois la barrière réduite, il est courant de constater que les courants ne se stabilisent qu’après que le dispositif s’est clairement éloigné du régime de faible courant dominé par l’annulation. À 0,8 V (figure ??), la diode fonctionne en polarisation directe importante. Les densités de porteurs sont élevées dans toute la structure, et les courants d’électrons et de trous deviennent importants, plus lisses et plus uniformes dans les régions quasi neutres. Même ici, les simulations avec du GaAs peuvent conserver une certaine ondulation résiduelle dans les profils de courant bruts, en particulier si le solveur affiche les courants ponctuels non filtrés ; l’élément de diagnostic essentiel est que la tendance à grande échelle est physiquement cohérente, avec un flux de courant continu et sans inversions de signe parasites dans le dispositif actif.

Ensemble, ces graphiques de densité de courant présentent une description cohérente du fonctionnement interne de la diode : ils montrent une annulation quasi parfaite des flux d’électrons et de trous à l’équilibre, puis une conduction directe limitée par l’injection, et enfin un transport en régime permanent à courant élevé sous une forte polarisation directe.