Guide sur le GaAs (Partie A) : Exécution d’une simulation 3D de dérive-diffusion avec un défaut vertical.

1. Introduction

Dans ce tutoriel, vous effectuerez une simulation 3D de dérive-diffusion d’une diode pn en GaAs, contenant un défaut vertical intentionnellement évident (un court-circuit de haut en bas à travers la structure). Il s’agit d’un problème délibérément « inéquitable » pour les modèles 1D : le défaut force un flux de courant latéral, une concentration de courant et une recombinaison localisée qui ne peuvent exister dans une simulation purement dans le sens de l’épaisseur. Pourquoi cela est important : dans les dispositifs réels, les performances sont souvent limitées non pas par la structure idéale des couches, mais par la manière dont le courant circule réellement. Les défauts locaux, les micro-trous, les fuites sur les bords, les contacts de taille finie et les matériaux non uniformes peuvent avoir une incidence importante sur les mesures de la courbe courant-tension (JV). Un modèle de dérive-diffusion 3D vous permet de relier directement les caractéristiques de la courbe courant-tension à des régions spécifiques du dispositif, ce qui permet de répondre à des questions auxquelles les modèles 1D ne peuvent pas répondre.

Comment utiliser cette approche dans votre propre travail : considérez le défaut présenté ici comme un modèle pour toute non-idéalité localisée, par exemple un filament de court-circuit, un trou dans une couche de blocage, une limite de grain, une zone de sonde endommagée ou une électrode de surface finie. Le processus est le même : créez une géométrie 3D minimale, maintenez la taille du maillage latéral aussi petite que possible, exécutez le solveur, puis utilisez les courbes de courant-tension résolues par contact et les instantanés (jn.csv, \(\phi\), recombinaison) pour établir un lien entre le comportement électrique et l’emplacement physique.

L’objectif est double : (1) vous fournir un point de départ pratique pour une modélisation électrique véritablement tridimensionnelle dans OghmaNano, et (2) vous enseigner la prise de décision en matière de modélisation : quand la 3D est essentielle (défauts, contacts finis, variation latérale) et quand elle peut être réduite à une dimension (1D). Dans la partie A, vous exécutez la simulation du dispositif défectueux, dans la partie B, vous supprimez le défaut et vous montrez que les résultats en 3D, 2D et 1D concordent une fois la symétrie rétablie, et dans la partie C, vous activez l’éclairage et vous observez où apparaît la sensibilité numérique à faible courant.

2. Premiers pas

Dans la fenêtre principale d’OghmaNano, cliquez sur « Nouvelle simulation » pour ouvrir la bibliothèque d’appareils (voir ??). Double-cliquez sur « GaAs demos ». Cela ouvre un dossier contenant des exemples d’arséniure de gallium (voir ??). Enfin, double-cliquez sur « GaAs - 3D defect » et enregistrez la simulation dans un dossier auquel vous avez accès en écriture.

Une fois l’exemple ouvert, la fenêtre principale affiche une vue en 3D du dispositif (voir ??). La structure est une diode pn en GaAs, contenant deux blocs visibles qui représentent une zone de défaut vertical intentionnel (par exemple, un filament d’impureté ou un trajet de perte similaire à un court-circuit). Le défaut s’étend de haut en bas à travers le dispositif, ce qui le rend intrinsèquement tridimensionnel et lui permet de créer une concentration de courant et une recombinaison spatialement localisée, phénomènes qui seraient absents dans un modèle purement unidimensionnel.

Pour vérifier les caractéristiques du matériau et les paramètres relatifs aux défauts, accédez à l’onglet « Structure du dispositif » et cliquez sur « Paramètres électriques » pour ouvrir l’éditeur de paramètres (voir ??). Vous y trouverez des ensembles de paramètres distincts pour :

• GaAs de type p+ et GaAs de type n+ (les couches de la diode), et
• défaut0 et défaut1 (les zones présentant des défauts).

L’exemple utilise des valeurs de mobilité et de densité de porteurs réalistes et typiques du GaAs. La recombinaison inclut un canal bimoléculaire simple de type « libre-libre » (souvent désigné dans le contexte des dispositifs organiques par n3–p3) avec un coefficient représentatif de 1 × 10−15 m3s−1.

Vous remarquerez également que les pièges d’équilibre de Shockley-Read-Hall (SRH) sont activés. Dans de nombreuses simulations de dispositifs organiques, les pièges dynamiques de SRH sont nécessaires, car l’occupation des pièges emmagasine une quantité importante de charge et influence les propriétés électrostatiques. Cependant, pour le GaAs, les densités de pièges sont généralement suffisamment faibles pour que les pièges puissent souvent être traités comme un simple mécanisme de perte par recombinaison. L’utilisation des pièges d’équilibre de SRH permet de prendre en compte cette physique sans introduire de variables d’état supplémentaires pour la charge emmagasinée dans les pièges, ce qui permet de maintenir un modèle à la fois stable et efficace.

3. Exploration de la simulation : dopage et maillage.

La barre d’outils Électrique permet d’accéder aux principaux outils d’édition électrique utilisés pour définir et analyser la simulation (voir ??). Dans cette section, nous nous concentrerons sur deux outils : Doping/Ions, qui définit la jonction pn, et Electrical mesh, qui contrôle la résolution spatiale du solveur de dérive-diffusion.

Cliquez sur « Doping/Ions » pour ouvrir l’éditeur de dopage (voir ??). Cet exemple présente une diode pn volontairement simple, avec des régions de type p et de type n clairement séparées. En maintenant un profil de dopage simple ici, il est plus facile d’attribuer tout comportement inhabituel ultérieur à la géométrie (le défaut) plutôt qu’à l’électrostatique.

Ensuite, cliquez sur « Grille électrique » pour ouvrir l’éditeur de grille (voir ??). Vous constaterez que les axes x, y et z sont tous activés, ce qui signifie que la simulation est véritablement tridimensionnelle. Dans ce système de coordonnées :

• y représente la direction principale du transport (à travers l’épaisseur de la diode).
• Les axes x et z représentent les directions latérales utilisées pour corriger le défaut et contrôler le flux de courant latéral.

Le maillage dans la direction y est délibérément divisé entre les couches de type n et de type p, de manière à ce que chaque région conserve une résolution appropriée. Cela devient important si vous modifiez ultérieurement l’épaisseur des couches ou si vous introduisez une asymétrie. En revanche, les maillages latéraux (x et z) utilisent relativement peu de points (par exemple, environ 10 dans la direction x et environ 5 dans la direction z). C’est intentionnel : en 3D, le coût de calcul augmente rapidement. Si vous augmentez la densité du maillage dans chaque dimension d’un facteur \(k\), le nombre total d’inconnues augmente approximativement selon le facteur \(k^3\). En termes pratiques, les simulations 3D deviennent beaucoup plus lentes à mesure que vous affinez le maillage. Il est donc préférable de maintenir la résolution latérale aussi faible que possible, tout en conservant une résolution suffisante pour les phénomènes physiques qui vous intéressent.

4. Lancer la simulation 3D et afficher les résultats.

Pour lancer la simulation, cliquez sur « Lancer la simulation » (le triangle bleu). Étant donné qu’il s’agit d’un problème de dérive-diffusion véritablement tridimensionnel, le processus prendra plus de temps qu’une simulation 1D. Sur un ordinateur portable relativement moderne, l’exemple devrait se terminer en environ 30 secondes. Si le temps d’exécution est nettement plus long, vérifiez d’abord que la simulation est enregistrée sur un disque local ; les dossiers montés sur un réseau ou synchronisés avec le cloud peuvent ralentir considérablement l’exécution en raison d’un volume important d’opérations d’entrée/sortie de fichiers. Une fois l’exécution terminée, ouvrez l’onglet « Sortie » (voir ??). Ce répertoire contient les fichiers de sortie standard générés par la résolution du problème de dérive-diffusion, notamment des instantanés des variables internes et des données de courant-tension. À ce stade, il est important de faire la distinction entre les données agrégées et les données résolues au niveau des contacts.

Étant donné que cet appareil est simulé en 3D avec des contacts de surface finie, il est généralement préférable d’ignorer le fichier global jv.csv et d’utiliser plutôt les courbes courant-tension (JV) résolues pour chaque contact. Les fichiers jv_contact0.csv et jv_contact1.csv indiquent respectivement le courant associé aux contacts supérieur et inférieur. Dans les simulations multidimensionnelles, ces courbes constituent la méthode la plus pertinente sur le plan physique pour interpréter le comportement électrique et diagnostiquer les pertes spécifiques à chaque contact. Ouvrez jv_contact0.csv pour afficher la courbe courant-tension du contact supérieur (voir ??). Étant donné que l’appareil contient un défaut vertical, cette courbe ne doit pas ressembler à la caractéristique idéale d’une diode pn unidimensionnelle. Dans la section suivante (Partie B), vous désactiverez le défaut et vous observerez comment la réponse courant-tension change, et comment la même structure revient à une configuration bidimensionnelle et unidimensionnelle une fois que l’asymétrie tridimensionnelle inhérente est supprimée.

5. Analyse des variables internes à l’aide de la fonction d’affichage des captures d’écran.

Pour examiner le fonctionnement interne du solveur, ouvrez le répertoire snapshots/ dans le dossier de sortie. Vous aurez ainsi accès aux variables du simulateur, présentées en fonction de leur position et de l’étape de simulation (généralement la tension appliquée). En double-cliquant sur le répertoire, vous ouvrirez le visualiseur d’instantanés présenté dans ?? et ??.

Dans la fenêtre des captures d’écran, cliquez sur le bouton bleu « + » pour ajouter un graphique et sélectionnez jn.csv dans le menu déroulant « Fichier à utiliser pour le graphique ». Le fichier jn.csv correspond à la densité de courant électronique verticale. À l’aide du curseur principal, vous pouvez observer l’évolution de la distribution du courant en fonction de la tension appliquée, tandis que les curseurs y et z vous permettent de réaliser des coupes dans le volume 3D.

Cet exemple illustre une diode en l’absence de lumière, de sorte que tout le courant est dû à des effets électriques plutôt qu’à un photocourant. L’observation essentielle est que le défaut crée un chemin de conduction très localisé, qui peut être visualisé directement en 3D et associé de manière non équivoque au comportement non idéal de la courbe courant-tension (JV). Le visualiseur d’instantanés peut être utilisé pour explorer de nombreuses autres grandeurs internes, notamment Ec.csv et Ev.csv (bords de bande), les niveaux quasi-Fermi, les taux de recombinaison et les densités de porteurs. L’étude de ces paramètres en parallèle avec la courbe courant-tension est l’une des méthodes les plus efficaces pour mieux comprendre le fonctionnement des solveurs de dérive-diffusion dans les dispositifs multidimensionnels.

Il est important de noter que, bien que le potentiel électrostatique \(\phi\) évolue en fonction de la tension appliquée, la distribution d’impuretés (Nad.csv) est fixée par l’éditeur d’impuretés et ne change pas avec la tension. Néanmoins, la visualisation de la distribution d’impuretés en 3D est utile, en particulier pour les tutoriels ultérieurs qui traitent des variations latérales de la distribution d’impuretés ou des défauts localisés dans l’espace.