Simulation de cellules solaires organiques à grande surface (PM6:Y6) : des réseaux de contacts tridimensionnels aux dispositifs éclairés.
1. Introduction
Dans ce tutoriel, nous poursuivons la série sur... simulation de dispositifs sur de grandes surfacesen mettant l’accent sur les dispositifs qui ont dépassé le stade des petites cellules de laboratoire pour atteindre des géométries adaptées à des applications à grande échelle. Il se situe entre deux processus existants décrits dans le manuel : conception de contacts de grande surface (sans dispositifs actifs)et modules pérovskites tridimensionnels complexesCet exemple fait le pont entre les deux concepts, en introduisant les diodes et l’absorption de la lumière dans un système conçu précédemment. structure de contactmais sans aller jusqu’à une solution complète. module interconnecté.
Dans cet exemple, nous simulons une seule cellule solaire organique de grande surface, basée sur une architecture de type PM6:Y6, avec des contacts métalliques hexagonaux imprimés. L’objectif est de comprendre comment un dispositif qui fonctionne bien sur une petite surface se comporte lorsqu’il est agrandi, et comment le flux de courant latéral, la chute de tension et la résistance de contact commencent à prédominer. Bien que ce tutoriel utilise des matériaux OPV, le processus est général : en remplaçant les données optiques n/k et en ajustant les valeurs des composants, le même modèle peut représenter votre dispositif de recherche.
Pour appliquer ce processus à votre propre dispositif, vous devez simplement fournir deux ensembles de données : les résistivités des couches (généralement issues de la littérature ou de structures de test simples), et un petit nombre de paramètres de diode, tels que le facteur d’idéalité, le photocourant et le courant de saturation. Une fois ces paramètres définis, le modèle vous permet de simuler l’agrandissement d’un dispositif à l’échelle du laboratoire et d’identifier les mécanismes de perte qui limiteront son efficacité lorsqu’il sera utilisé dans des géométries réelles de grande taille.
Dans ce cas, nous n’utilisons pas le solveur 3D complet de diffusion-dérive, bien qu’OghmaNano en propose un. Pour les dispositifs de grande surface, la résolution complète en 3D du modèle de diffusion-dérive devient rapidement trop coûteuse en termes de calcul et est souvent inutile, car les principales limitations proviennent plutôt du transport latéral et de la résistance de contact que du mouvement des porteurs à travers chaque couche. Au lieu de cela, OghmaNano construit un modèle de circuit 3D tenant compte de la géométrie, dans lequel le transport latéral et les contacts sont représentés par des éléments résistifs, et le comportement photovoltaïque est modélisé à l’aide de diodes intégrées.
Bien que le problème électrique soit résolu à l’aide d’une formulation de circuit, l’optique est traitée de manière rigoureuse. La photogénération est calculée à l’aide d’un solveur optique basé sur une matrice de transfert, évalué en coupes unidimensionnelles sur la géométrie tridimensionnelle, ce qui permet de conserver les effets d’interférence et d’absorption des couches minces, tandis que le solveur de circuit gère la collecte de courant sur une grande surface.
2. Créez la simulation.
Dans la fenêtre principale, cliquez sur « Nouvelle simulation ». Vous verrez une liste des catégories de simulation (voir ??). Ensuite, double-cliquez sur « Cellule solaire PM6:Y6 à grande surface » (voir ??).
Une fois la simulation créée, la fenêtre principale s’ouvre et affiche la géométrie du dispositif en trois dimensions (voir ??). La couche de contact supérieure est une grille métallique hexagonale (en forme de nid d’abeilles) étiquetée Ag. Cette couche d’argent agit comme une grille de collecte de courant hautement conductrice. Étant donné qu’elle est réalisée sous forme de grille, la majeure partie de la surface reste ouverte, ce qui permet à la lumière de traverser la couche de contact. Sous la grille d’argent se trouve le matériau PEDOT:PSS, qui sert de couche de diffusion latérale du courant. Comme dans le tutoriel sur la couche de contact de grande surface, la charge est d’abord collectée latéralement dans la couche de polymère, puis transférée dans la grille d’argent à faible résistance pour une extraction efficace.
La couche photoactive est constituée de PM6:Y6. Sur le plan électrique, cette couche est représentée par un élément de diode avec génération optique, tandis que des composants résistifs supplémentaires permettent de prendre en compte sa conductivité finie dans les calculs. Au fond de la pile se trouve un contact continu en argent qui couvre toute la surface du dispositif et assure la connexion électrique.
La barre verte visible dans la fenêtre principale peut sembler flotter, mais il s’agit simplement d’un indicateur visuel de la condition limite d’extraction. Elle délimite la zone dans laquelle le solveur applique le contact électrique utilisé pour extraire le courant du réseau de circuits.
3. Vérifiez l’état de la superposition des couches épitaxiales et des contacts.
Cette simulation (et l’exemple précédent concernant une grande zone de contact) est définie à l’aide d’une approche par couches. épitaxial structure. L’appareil est constitué d’une pile de couches minces, réalisée à l’aide de l’éditeur de couches, ce qui garantit que chaque couche est parfaitement superposée à la couche suivante, sans espaces ni variations d’épaisseur. Cette approche est particulièrement adaptée aux dispositifs planaires de grande surface, où la pile verticale est fixe et le comportement latéral est introduit par le biais du modèle de circuit. Elle diffère des simulations ultérieures, qui utilisent un modèle modulaire (voir la suite). tutoriel sur les modules interconnectés), où les contacts, les barres omnibus et les configurations plus complexes sont définies de manière explicite à l’aide d’une géométrie 3D complète.
Pour examiner la structure, ouvrez l’éditeur de calques afin d’afficher la pile de calques (voir ??), puis ouvrez le panneau Contacts dans l’éditeur de structure du dispositif (voir ??). Sur le plan géométrique, la configuration est simple, mais importante : la contact inférieur couvre toute la surface du dispositif, tandis que le contact supérieur est beaucoup plus petit et définit la région où le courant est extrait, ce qui correspond au marqueur d’extraction vert affiché dans la vue 3D. Il est essentiel de comprendre la géométrie de ces contacts, car elle détermine les trajets du courant latéral qui influencent le comportement des dispositifs de grande surface.
4. Vérification de la structure en treillis.
Passez à l’onglet « Schéma de circuit », puis cliquez sur l’icône de rafraîchissement située en bas à gauche (le bouton en forme de symbole de recyclage). Cette action permet de créer la représentation du circuit à partir de la géométrie et de l’agencement des couches définies. La grille de circuit résultante est affichée dans ??. Chaque élément de la grille correspond à un composant du circuit, tel qu’une résistance ou une diode, et la couleur indique la couche ou la région du matériau dont il provient.
En faisant pivoter la vue, vous pouvez clairement identifier les zones de contact supérieure et inférieure. Dans les figures ?? et ??, les points bleus indiquent les nœuds à partir desquels le courant est extrait du maillage du circuit. Ces nœuds définissent les conditions aux limites de contact et représentent les points où le solveur de circuit collecte le courant du dispositif.
5. Paramètres électriques et génération optique.
Ouvrez les paramètres électriques à partir du panneau de structure du dispositif. Cet éditeur contrôle les propriétés de la résistance et de la diode utilisées par le modèle de circuit à grande surface. Des exemples de paramètres sont présentés dans ?? et ??. Pour la couche PM6:Y6, vous verrez les paramètres spécifiques à un élément de diode, notamment le facteur d’idéalité n et le courant de polarisation inverse (de saturation) I0. Si vous remplacez ces valeurs par des valeurs issues de vos propres courbes expérimentales, vous pourrez constater comment votre dispositif se comportera à plus grande échelle. Vous pouvez également activer la génération optique de porteurs de charge. Lorsque la génération optique est activée, le solveur optique à matrice de transfert fournit un terme de photogénération qui est intégré au modèle électrique sous la forme d’un photocourant Iph dans l’équation de la diode éclairée :
I(V) = I0 [exp(qV / (n k T)) − 1] − Iph
Dans ce tutoriel, l’interrupteur permet de contrôler efficacement le terme Iph. Lorsque la génération optique est désactivée, le circuit se comporte comme un réseau de diodes en mode bloqué. Lorsque la génération optique est activée, l’illumination injecte un courant dans le circuit par le biais de l’élément diode, et la valeur de Iph est déterminée par le flux de photons absorbés par l’optique à matrice de transfert.
💡 Détail de la mise en œuvre : la couche active peut être divisée en plusieurs éléments résistifs à travers son épaisseur afin de prendre en compte le transport volumique et la résistance en série. Toutefois, le modèle de circuit contient toujours une seule couche de diode, généralement placée au fond de la région active. Cela reflète le fait que la diode représente la réponse globale courant-tension et la photogénération du dispositif, tandis que les éléments résistifs tiennent compte de la chute de tension et de la répartition du courant à l’intérieur du matériau.
6. Lancer la simulation.
Pour lancer la simulation, cliquez sur « Exécuter la simulation » (le triangle bleu dans la barre d’outils principale) ou appuyez sur la touche F9. Le solveur commencera immédiatement à s’exécuter et la fenêtre de sortie sera mise à jour en temps réel (voir ??).
Les lignes colorées qui apparaissent au début de la simulation correspondent au calcul optique. À cette étape, le solveur de matrice de transfert évalue la photogénération dans tout le dispositif, section par section, afin de créer le profil de génération résolu spatialement. Une fois cette opération terminée, le solveur passe à la résolution du problème électrique et commence à analyser le réseau de circuits.
Il est important de surveiller les données affichées à la fin de la simulation. Ces informations ne sont pas du bruit parasite, mais elles fournissent des indications précieuses sur le comportement physique de la simulation. Dans la capture d’écran présentée, on peut constater qu’environ −9,48 × 102 A m−2 s’échappe d’un des contacts, tandis que −9,48 × 101 A m−2 s’échappe de l’autre. Au même moment, le résidu du solveur indiqué, F = 1,7 × 10−14, révèle que la continuité du courant est respectée avec une précision extrêmement élevée, et que le premier point de la courbe J-V a été calculé en environ 17 ms.
La différence de densité de courant entre les deux contacts est tout à fait normale. Le contact inférieur couvre toute la surface du dispositif, tandis que le contact supérieur n’occupe qu’une petite zone. Par conséquent, le même courant total est extrait sur une surface beaucoup plus petite au niveau du contact supérieur, ce qui entraîne une densité de courant nettement plus élevée à cet endroit. L’analyse des résultats de simulation de cette manière permet de vérifier rapidement et efficacement que la géométrie et les conditions aux limites de la simulation sont physiquement cohérentes. Si l’erreur résiduelle devait augmenter jusqu’à atteindre une valeur proche de l’unité, ou si le temps d’exécution devenait anormalement long, cela indiquerait généralement un problème, tel qu’un maillage défectueux ou un contact déconnecté. Une fois la simulation terminée, ouvrez le fichier. Résultat tab (voir) ??). Les différents résultats possibles sont décrits plus en détail dans le tutoriel précédent et approfondi dans prochain tutorielIci, nous nous concentrons uniquement sur les courbes caractéristiques tension-courant obtenues en résolvant les contacts individuellement.
Chaque contact électrique possède son propre fichier JV, jv_contact0.csv et jv_contact1.csv, qui enregistre le courant qui s’écoule de ce contact spécifique. Bien que les densités de courant soient différentes, le courant total extrait est le même une fois que les surfaces de contact sont prises en compte. L’ouverture du fichier jv_contact0.csv permet d’afficher la courbe JV présentée ci-dessous.
7. Modification des paramètres du matériau.
Il est possible de modifier directement les propriétés des matériaux à l’aide de l’éditeur d’objets. En faisant un clic droit sur n’importe quelle couche dans la vue 3D et en sélectionnant « Modifier l’objet », l’éditeur d’objets s’ouvre (voir ?? et ??).
À partir de cet écran, vous pouvez modifier le matériau optique attribué à l’objet, modifier sa forme géométrique (par exemple, lorsque vous utilisez des lentilles à motifs, comme des structures en nid d’abeille), ou ajuster les attributs visuels tels que la couleur. Dans cette simulation de type épitaxial, la position absolue de l’objet est largement déterminée par l’empilement des couches, de sorte que les paramètres tels que la translation et la rotation ont généralement peu d’effet. Ce qui compte le plus, ce sont l’attribution du matériau optique et la géométrie de l’objet.
💡 Essayez ceci : faites un clic droit sur la couche PM6:Y6 et remplacez son matériau optique par un pérovskite. Relancez ensuite la simulation et observez comment la courbe courant-tension (JV) évolue. Si vous le souhaitez, vous pouvez également mettre à jour les paramètres de la diode (par exemple, le facteur d’idéalité, le courant de saturation et le photCourant) en utilisant les valeurs indiquées dans la littérature sur les pérovskites.
📘 Retour à introduction à l’électronique organique pour une vue d’ensemble plus large du sujet.