Simulation d’une diode à jonction PN en InP : dérive-diffusion, rectification et recombinaison des porteurs de charge.

1. Introduction : une jonction redresseuse en InP.

Le phosphure d’indium (InP) est un semi-conducteur de type III-V à bande interdite directe, utilisé comme base pour les dispositifs électroniques à haute vitesse et les technologies photoniques intégrées. Une diode à jonction PN est l’une de ses structures de dispositifs les plus simples et les plus instructives : une région de type p rencontre une région de type n, ce qui crée un champ électrique interne et une barrière de potentiel qui contrôle le mouvement des électrons et des trous. La diode conduit fortement en polarisation directe et supprime l’injection de porteurs en polarisation inverse, ce qui produit une rectification électrique.

L’InP présente, à température ambiante, une largeur de bande interdite d’environ 1,34 eV et une mobilité électronique relativement élevée. Ces propriétés font de la conduction dans les jonctions InP un élément important dans l’électronique des semi-conducteurs composés et permettent de mieux comprendre les interfaces semi-conductrices utilisées dans les plateformes d’appareils à base d’InP. L’exemple présenté dans ce tutoriel illustre le fonctionnement électrique d’une jonction homojonction compacte p-InP/n-InP dotée d’une couche de contact n⁺.

Au lieu d’ajuster la courbe J–V à une simple équation de diode idéale, OghmaNano résout les équations couplées de Poisson et de continuité de dérive-diffusion. Les données obtenues révèlent à la fois le courant aux bornes et la distribution interne des énergies des bandes de conduction et de valence, les densités d’électrons et de trous, les niveaux quasi-Fermi et la recombinaison. Cela nous permet d’expliquer pourquoi la diode effectue une rectification, plutôt que de simplement constater qu’elle le fait.

La fenêtre principale de simulation, illustrée dans ??, affiche la structure InP stratifiée verticalement. Bien que la géométrie soit représentée en trois dimensions, le problème électrique résolu ici est unidimensionnel, le long de la direction de croissance.

2. Créer une nouvelle simulation.

Dans la fenêtre principale d’OghmaNano, sélectionnez « Nouvelle simulation ». L’écran de sélection des simulations s’ouvre, affichant les différentes catégories d’exemples (??). Ouvrez la section « Démonstrations III–V/InP » et choisissez l’exemple de diode InP PN (??). L’emplacement exact de l’exemple peut varier selon les versions d’OghmaNano.

Une fois l’exemple ouvert, sélectionnez l’onglet « Structure du dispositif » et ouvrez l’éditeur de couches. Vérifiez les noms des couches, leur épaisseur et les modèles de matériaux électriques associés avant de lancer la simulation. Les sections suivantes expliquent la contribution de chaque région à la réponse globale de la diode.

3. Analyse de la structure du dispositif, des contacts et du dopage.

3.1 Structure de la couche semi-conductrice

Le dispositif électrique contient trois couches d’InP, pour une épaisseur totale de matériau semi-conducteur de 700 nm. Les couches p-InP et n-InP forment la jonction homojonction à effet redresseur, tandis que la région n⁺-InP constitue une région cathodique plus fortement dopée. La disposition des couches est visible dans l’éditeur (??).

La zone d’air de 200 nm affichée dans l’éditeur de structure fait partie de l’affichage et de la géométrie optique environnante, et non de la pile électrique du semi-conducteur. Par conséquent, son épaisseur n’est pas prise en compte dans l’épaisseur totale de 700 nm de l’appareil électrique.

3.2 Dopage et effets électrostatiques inhérents

Ouvrez « Electrical » puis « Doping / Ions » pour afficher les concentrations de donneurs et d’accepteurs (??). La variation nette du dopage change de signe à l’interface entre le semi-conducteur de type p-InP et le semi-conducteur de type n-InP, à environ 250 nm de la partie supérieure de la pile de semi-conducteurs électriquement actifs. La deuxième variation, entre le semi-conducteur de type n-InP et le semi-conducteur de type n⁺-InP, modifie la densité d’électrons en direction de la cathode, mais ne crée pas une autre jonction PN.

À l’équilibre, les électrons ont tendance à diffuser de la région de type n vers la région de type p, tandis que les trous diffusent dans la direction opposée. Les impuretés ionisées non compensées génèrent un champ qui s’oppose à une diffusion supplémentaire. La zone de déplétion et la barrière de potentiel interne qui en résultent constituent la base de la réponse de redressement.

Ici, \(\phi\) représente le potentiel électrostatique, \(\varepsilon\) la permittivité diélectrique et \(\rho\) la densité de charge locale. Les dopants ionisés et les porteurs de charge mobiles déterminent le champ électrique local, tandis que la tension appliquée modifie la barrière effective pour l’injection des porteurs de charge.

3.3 Contacts électriques

Le modèle utilise une barrière ohmique où les trous sont majoritaires au niveau de l’anode et une barrière ohmique où les électrons sont majoritaires au niveau de la cathode. Ces conditions aux limites sont appropriées pour étudier la réponse de la jonction sans introduire de barrière d’injection supplémentaire au niveau des contacts. Dans les dispositifs fabriqués, la métallisation et le traitement des contacts déterminent à quel point les contacts se rapprochent d’un comportement ohmique idéal.

4. Modèles de matériaux et paramètres de transport

4.1 Modèles de matériaux Lua : origine des paramètres de l’InP.

L’une des caractéristiques importantes de cet exemple de diode InP est que les paramètres des matériaux électriques sont définis à l’aide de modèles de matériaux Lua. Cela établit un lien clair entre les propriétés physiques d’un semi-conducteur et la solution de diffusion-dérive. Les régions p-InP, n-InP et n⁺-InP font toutes référence à la même définition de matériau, Inorganic/InP/InP.lua. Leur comportement électrique différent découle principalement du profil de dopage, plutôt que de l’attribution de trois ensembles distincts de constantes de matériau.

Pour examiner le modèle, ouvrez la section « Structure du dispositif » → « Paramètres électriques » (l’éditeur de matériaux électriques) et sélectionnez l’une des couches d’InP. L’onglet « Éditeur de code » affiche le script Lua du matériau associé à cette région ; l’onglet « Paramètres » affiche les propriétés du matériau calculées. Ces deux vues doivent être utilisées ensemble. Le code définit la manière dont les propriétés sont obtenues, tandis que les valeurs calculées révèlent ce que le solveur électrique reçoit réellement.

Les images présentées dans ?? et ?? illustrent la manière d’accéder à l’éditeur de code. La configuration de la couche de la simulation est la référence pour le modèle de matériau attribué : dans cet exemple, chaque région active est associée au fichier Lua InP. Les noms des fichiers image sont hérités de la séquence de captures d’écran et ne déterminent pas le type de semi-conducteur simulé.

Dans la simulation fournie, les commutateurs db_overrides sont réglés sur false pour les principaux paramètres de transport et de bande. Par conséquent, OghmaNano est configuré pour obtenir les valeurs correspondantes à partir de la base de données de matériaux/du modèle Lua, au lieu d’utiliser les valeurs par défaut enregistrées ailleurs dans le fichier JSON. Cela est important lors de l’examen de la configuration : les valeurs effectives affichées dans l’onglet Paramètres sont plus informatives que les valeurs par défaut stockées dans un champ de remplacement inactif.

4.2 Structure de bandes, densité d’états et mobilité.

Dans l’onglet « Paramètres évalués », les premières valeurs à vérifier sont l’énergie de bande interdite \(E_g\) et l’affinité électronique \(\chi\). L’énergie de bande interdite définit la séparation entre les bandes de conduction et de valence, tandis que l’affinité électronique indique la position du bord de la bande de conduction par rapport à une référence d’énergie du vide. Pour l’InP, l’énergie de bande interdite à température ambiante est d’environ \(1.34\,\mathrm{eV}\). Une vérification initiale utile consiste à déterminer si la séparation entre \(E_c\) et \(E_v\) dans le diagramme de bandes obtenu correspond à l’énergie de bande interdite évaluée.

Les densités d’états efficaces des bandes de conduction et de valence, \(N_c\) et \(N_v\), influencent la relation entre la concentration de porteurs et la position du niveau quasi-Fermi. Elles contribuent donc aux populations d’équilibre des électrons et des trous, ainsi qu’à la réponse à l’injection sous polarisation directe. La permittivité relative \(\varepsilon_r\) intervient dans l’équation de Poisson et influence le potentiel électrostatique et la largeur de la zone de déplétion. Ces quantités doivent être considérées comme un ensemble cohérent, et non comme des paramètres de réglage numériques indépendants.

Les mobilités des électrons et des trous, \(\mu_n\) et \(\mu_p\), déterminent la réponse au déplacement des porteurs de charge dans un champ électrique ; la diffusion dépend également des gradients de densité des porteurs de charge. Dans le modèle unidimensionnel, les densités de courant sont exprimées comme suit :

Ici, \(n\) et \(p\) représentent les densités locales d’électrons et de trous, \(E\) représente le champ électrique, et \(D_n\) et \(D_p\) représentent les coefficients de diffusion. Les deux composantes du courant peuvent varier considérablement dans la diode, même lorsque le courant total en régime permanent est approximativement indépendant de la position. C’est pourquoi les sections suivantes examinent \(J_n\) et \(J_p\) en même temps que les profils de densité de porteurs.

4.3 Comment vérifier et comparer les paramètres des matériaux.

Dans l’éditeur électrique, sélectionnez le matériau p-InP et examinez son code Lua ainsi que ses paramètres calculés. Ensuite, passez au matériau n-InP, puis enfin au matériau n⁺-InP. Comparez le modèle de matériau, la largeur de bande interdite, l’affinité électronique, les mobilités, les densités d’états effectives et les coefficients de recombinaison utilisés. Pour cette hétérojonction, le même modèle de base InP est utilisé dans les trois régions ; la principale modification intentionnelle entre les régions est leur concentration d’accepteurs ou de donneurs.

Il est utile de distinguer deux concepts : le modèle des matériaux décrit les propriétés de l’InP, tandis que le profil de dopage détermine quels types de porteurs de charge prédominent localement et où se trouve la jonction. Cette distinction devient particulièrement pertinente lorsqu’on passe d’une jonction homojonction PN à une hétérojonction à base d’InP, où les différences entre les modèles des matériaux peuvent elles-mêmes générer des décalages de bande. Dans cette diode, la courbure des bandes de conduction et de valence, visible sur les schémas, reflète principalement les effets électrostatiques à travers la transition p/n.

Pour obtenir des résultats reproductibles, conservez le modèle de matériau Lua constant lorsque vous modifiez les concentrations de dopage. Inversement, si vous étudiez la sensibilité du résultat à la mobilité ou au coefficient de recombinaison, modifiez délibérément ce paramètre du matériau et comparez la courbe J–V, les densités de porteurs et les profils de courant interne obtenus. L’éditeur de paramètres et les vues du code Lua établissent un lien entre les hypothèses physiques utilisées comme données d’entrée et les résultats simulés.

4.4 Paramètres de recombinaison

Les régions d’InP comprennent la recombinaison de Shockley-Read-Hall (SRH) et un canal de recombinaison des porteurs libres, tel que décrit dans le modèle de matériau utilisé. La recombinaison SRH décrit la capture des électrons et des trous facilitée par les défauts ; la recombinaison porteur-porteur est particulièrement importante dans l’InP à bande interdite directe, car les transitions électron-trou peuvent se produire de manière radiative. Dans le cadre de cet exercice, les deux processus sont utiles en tant que mécanismes de perte de porteurs qui modifient le courant électrique.

Le coefficient \(B\), les durées de vie des pièges \(\tau_n\) et \(\tau_p\), ainsi que les densités auxiliaires \(n_1\) et \(p_1\) déterminent l’importance de chaque processus. L’élément clé de cet exercice pratique est de comprendre comment la distribution spatiale de la recombinaison évolue lorsque la polarisation entraîne un nombre croissant d’électrons et de trous dans la même région.

4.5 Maillage et solution numérique

Le calcul électrique utilise un maillage unidimensionnel de 150 points répartis sur 700 nm, ce qui correspond à un espacement moyen d’environ 4,7 nm. Les équations non linéaires couplées sont résolues à l’aide de l’itération de Newton et du solveur linéaire creux UMFPACK. Ce maillage constitue un point de départ utile pour l’étude de la jonction et des régions quasi-neutres.

Les résultats comprennent des informations sur la convergence et des profils de dispositif en fonction des décalages. Une première vérification consiste à examiner la courbe J–V pour s’assurer qu’elle est régulière, puis à confirmer que les courants d’électrons et de trous respectent le principe général de conservation du courant en régime permanent. L’éditeur de maillage est présenté dans ??.

5. Exécution de la simulation et interprétation de la courbe courant-tension.

Retournez à la fenêtre principale et sélectionnez « Exécuter la simulation ». OghmaNano effectue une variation de la tension appliquée et enregistre la solution convergée à chaque point. Le fichier d’exemple enregistré utilise une variation de tension de 0 à 3 V, avec des incréments de 25 mV. La fenêtre de terminal enregistre l’avancement du solveur pendant la variation, tandis que l’onglet « Sortie » permet d’accéder aux données obtenues.

Une fois l’exécution terminée, les résultats sont affichés dans l’onglet « Résultats » (??). Ouvrez jv.csv pour afficher un graphique de la densité de courant en fonction de la tension appliquée. Les captures d’écran comprennent à la fois une représentation linéaire de la courbe J–V (??) et une représentation logarithmique du courant (??).

5.1 Comprendre les caractéristiques.

À basse tension appliquée, la barrière de jonction intégrée limite fortement l’injection de porteurs. L’augmentation de la tension de polarisation directe réduit cette barrière, ce qui facilite l’entrée des électrons et des trous dans la zone opposée, et le courant total augmente. Dans le régime intermédiaire, le courant est souvent approximativement proportionnel à la tension ; pour des courants plus importants, les pertes dues au transport et à la résistance peuvent modifier la pente de la courbe.

Cette équation de la diode idéale constitue une référence utile pour l’analyse d’un intervalle à décroissance exponentielle approximative, mais elle n’est pas imposée par le modèle de dérive-diffusion. Le facteur d’idéalité \(n\) et le coefficient \(J_0\) ne doivent être déterminés que sur un intervalle de polarisation approprié ; ils peuvent varier en fonction de la recombinaison, de la résistance série et du comportement du contact.

Le graphique logarithmique permet de déterminer le début d’un courant significatif et les variations de pente. Le graphique linéaire facilite l’interprétation de l’amplitude et de la forme du courant élevé. Pour obtenir une représentation claire du fonctionnement de la diode, examinez les deux graphiques en même temps que les profils de la bande interne et de la densité des porteurs.

6. Fonctionnement interne : comment fonctionne la diode.

Le répertoire « snapshots » contient des profils spatiaux enregistrés aux points de polarisation sélectionnés pendant le balayage J–V. Les captures d’écran fournies illustrent quatre positions au cours de ce balayage. Elles montrent une progression allant d’une faible injection vers une polarisation directe plus forte. Les quatre images de chaque groupe correspondent aux mêmes positions de curseur sélectionnées, ce qui permet de comparer directement les différentes grandeurs physiques.

6.1 Bandes de conduction et de valence, et niveaux quasi-Fermi.

Ouvrez le visualiseur d’instantanés et sélectionnez Ec.csv, Ev.csv, Fn.csv et Fp.csv. Il s’agit respectivement du bord de la bande de conduction, du bord de la bande de valence, du niveau quasi-Fermi des électrons et du niveau quasi-Fermi des trous. Leur séparation spatiale et leurs pentes révèlent la barrière de potentiel, la réponse à la tension appliquée et l’écart par rapport à l’équilibre.

À faible courant d’injection, les bords de la bande d’énergie se courbent de part et d’autre de la jonction en raison du champ électrostatique inhérent. Lorsque la polarisation directe augmente, la barrière effective diminue et les quasi-niveaux de Fermi des électrons et des trous se séparent à mesure que les porteurs traversent le dispositif. L’élément essentiel est que l’InP reste le même semi-conducteur des deux côtés de la jonction : les profils qui évoluent reflètent ses caractéristiques électrostatiques et l’injection de porteurs, qui sont déterminées par le dopage.

La comparaison entre ?? et ?? illustre la manière dont le biais terminal modifie la structure de bande interne. Plus précisément, les niveaux quasi-Fermi permettent de déterminer directement l’état hors équilibre établi pendant la conduction en régime permanent.

6.2 Densités d’électrons et de trous.

Ensuite, examinez les graphiques Q_nfree.csv et Q_pfree.csv, qui illustrent la répartition spatiale des électrons libres et des trous. À faible polarisation, les électrons prédominent sur le côté de type n et les trous sur le côté de type p. Leur densité diminue considérablement autour de la jonction, ce qui reflète la zone de déplétion.

À mesure que la polarisation directe augmente, les porteurs minoritaires sont injectés dans les régions opposées. Cela accroît le chevauchement des populations d’électrons et de trous et modifie la forme des deux profils de densité. Pour des niveaux d’injection plus importants, les populations de porteurs libres de part et d’autre de la jonction diffèrent considérablement de leurs valeurs à faible polarisation.

Comparez la variation de la densité au point de transition avec la position de la transition de dopage, puis avec les zones où la courbure de la bande est la plus forte. Cela permet de relier les porteurs fournis par le contact à la barrière électrostatique créée par la jonction PN.

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

Représentez les courbes Jn.csv et Jp.csv. Ces profils illustrent la contribution des électrons et des trous au courant électrique qui traverse le dispositif. Les différentes composantes peuvent varier en fonction de la position, car la recombinaison transfère les porteurs entre les populations d’électrons et de trous, mais le courant total à l’état stationnaire devrait rester approximativement constant en l’absence d’accumulation nette de charges.

À des courants très faibles, la compensation entre les termes de dérive et de diffusion peut donner l’impression que les profils de courant individuels sont bruyants ou présentent une structure très marquée. Les captures d’écran obtenues à faible polarisation mettent davantage en évidence ce comportement que les captures d’écran obtenues à des courants plus élevés. Il est utile d’examiner les profils à une échelle appropriée et de comparer le courant total à la valeur J–V au niveau de la borne avant d’interpréter les détails numériques.

6.4 Profil de recombinaison

La recombinaison détermine l’endroit où les électrons et les trous injectés sont éliminés de leurs populations de porteurs respectives. Ouvrez le visualiseur d’instantanés avec R_nfree_to_pfree.csv et R_ss_srh.csv pour comparer les canaux de recombinaison directe et les canaux assistés par défaut. Les deux processus dépendent de la disponibilité des électrons et des trous, de sorte que leur distribution spatiale évolue en fonction des changements dans les populations de porteurs injectées.

Au début de la conduction, la recombinaison peut être concentrée près de la jonction, là où les populations d’électrons et de trous se chevauchent. En cas d’injection directe plus importante, ce chevauchement s’étend sur une plus grande partie du semi-conducteur, et les profils de recombinaison peuvent devenir plus larges. Comparez la position et la largeur des pics de recombinaison avec les densités de porteurs pour la même polarisation.

Lors de la comparaison des deux canaux, il faut tenir compte à la fois de la valeur maximale et de l’intégrale spatiale. Un pic local marqué n’a pas nécessairement besoin de dominer la recombinaison intégrée si un autre canal est actif sur un volume beaucoup plus grand. Cet aspect est particulièrement important lorsqu’il s’agit de relier les profils de recombinaison au courant électrique total.

6.5 Profils des taux nets de création d’électrons et de trous.

Les captures d’écran étiquetées R_n.csv et R_p.csv fournissent des informations supplémentaires sur les taux utilisés dans les équations de continuité des électrons et des trous. Elles permettent de déterminer où l’équilibre des termes locaux de génération, de capture et de recombinaison est le plus marqué. Interprétez leurs signes en fonction de la convention de représentation graphique utilisée par OghmaNano.

Ensemble, la structure de bande, les densités de porteurs, les courants et les profils de variation permettent d’expliquer complètement les caractéristiques de la diode : la barrière interne contrôle l’injection, la tension appliquée modifie l’électrostatique interne, et le transport et la recombinaison déterminent le courant à l’état stationnaire qui en résulte.

7. Exercices et interprétation.

8. Interprétation du modèle et analyses complémentaires.

Cet exemple utilise un modèle unidimensionnel et à l’état stationnaire de dérive-diffusion pour l’InP, une structure compacte de semi-conducteur p/n/n⁺ et des contacts ohmiques idéalisés. Il est particulièrement adapté pour illustrer comment la rectification électrique émerge de la barrière de jonction et comment ce comportement est lié à la distribution spatiale des porteurs de charge. Des études plus approfondies sur les dispositifs peuvent intégrer des caractéristiques de contact mesurées, des paramètres de matériaux affinés, des variations de température et des détails géométriques supplémentaires.

Les propriétés du matériau InP sont obtenues à partir du modèle Lua sélectionné. Pour les études quantitatives, vérifiez les valeurs calculées de la bande interdite, de l’affinité électronique, de la mobilité et des coefficients de recombinaison dans l’éditeur de matériaux ou dans les données de sortie. Vérifiez la conservation du courant et les résidus du solveur, puis répétez les points de polarisation représentatifs en utilisant un maillage plus fin afin d’établir la stabilité numérique.

Dans l’exemple fourni, la plage de balayage standard s’étend de 0 à 3 V. Une expérience intéressante à mener ensuite consisterait à inclure une plage de polarisation inverse modérée, ce qui permettrait de comparer directement la conduction directe et le courant de fuite inverse. La plage de polarisation inverse et les caractéristiques de fuite pertinentes sur le plan physique dépendent du dispositif et des mécanismes pris en compte dans la simulation.

Pour approfondir l’analyse, faites varier les concentrations de dopage une par une. En général, l’augmentation de la concentration nominale de dopage réduit la zone de déplétion et modifie l’équilibre électrostatique ; la modification des coefficients de recombinaison influe sur le courant, qui dépend de l’injection. Ces comparaisons contrôlées sont utiles, car les matériaux et la géométrie de la couche peuvent rester constants.

9. Prochaines étapes

La diode InP PN constitue un point de départ utile pour la conception de structures semi-conductrices plus avancées à base d’InP. Une évolution naturelle consiste à introduire une région faiblement dopée et à comparer le comportement électrostatique et le comportement en matière de transport d’une jonction PN et d’une jonction PIN. Une autre approche consiste à étudier les hétérostructures impliquant l’InGaAs/InP, où les différences de composition entraînent des décalages de bande supplémentaires et de nouveaux comportements d’injection de porteurs.