Guide pratique sur la simulation d’excitons.
1. Introduction
Dans ce tutoriel, nous allons étudier comment simuler la dynamique des excitons dans les dispositifs photovoltaïques organiques (DPO). Lorsqu’un photon est absorbé dans un semi-conducteur organique, il génère une paire électron-trou liée, appelée exciton. excitonDans les dispositifs photovoltaïques organiques (OPV), ces excitons doivent se déplacer ou diffuser jusqu’à une interface donneur-accepteur, où ils peuvent se séparer en porteurs de charge libres. Ce processus peut être modélisé à l’aide d’un… dissociation d’excitons équation, qui décrit à la fois le transport des excitons et leur transformation finale en électrons et en trous. Dans OghmaNano, ces dynamiques des excitons sont intégrées aux équations standard de dérive-diffusion afin de prendre en compte l’ensemble du phénomène photophysique. Une discussion plus détaillée sur les moments où la modélisation des excitons est nécessaire (et sur les moments où elle peut être omise) est présentée dans Excitons et recombinaison geminate.
2. Premiers pas
À partir de Nouvelle simulation. dans le ruban du fichier, cliquez pour ouvrir Nouvelle simulation. fenêtre (voir) Figure 1a). Dans cette fenêtre, vous trouverez de nombreuses catégories d’exemples préconfigurés. Double-cliquez sur Simulations d’excitons dossier pour ouvrir le sous-menu des exemples disponibles relatifs aux excitons (voir Figure 1b). Pour ce tutoriel, nous utiliserons Dispositif à base d’excitons modèle, qui constitue un point de départ simple pour simuler la génération, le transport et la dissociation des excitons dans les structures OPV.
3. Lancez la simulation et examinez les résultats concernant les excitons.
Une fois l’appareil ouvert, vous pourrez voir la pile dans la fenêtre principale. Avant de lancer la simulation, allez dans l’onglet « Électrique » et vérifiez que le bouton « Solveur d’excitons » est enfoncé (il l’est généralement par défaut, mais il est préférable de le confirmer ??). Cela garantit que la dynamique des excitons est activée dans le modèle. Une fois cette opération effectuée, revenez au menu « Fichier » et appuyez sur le bouton « Lancer la simulation » pour exécuter le modèle ??.
Une fois la simulation terminée, vous pouvez examiner les résultats dans l’onglet « Résultats » (??). Celui-ci affiche la liste de tous les fichiers enregistrés sur le disque, y compris les résultats liés aux excitons. En double-cliquant sur jv.csv, vous pouvez afficher la courbe courant-tension (??). Cette courbe ressemble à une courbe courant-tension standard que vous obtiendriez à partir de n’importe quelle simulation de dispositif, mais dans ce cas, le solveur d’excitons a été activé, ce qui signifie que la dynamique des excitons est prise en compte dans la physique sous-jacente.
jv.csv. Bien qu’elle ressemble à une courbe J-V classique, le solveur d’excitons est activé.Si vous double-cliquez sur le dossier exciton_output situé dans le répertoire de sortie, vous pouvez accéder aux résultats détaillés du solveur d’excitons (voir ??). Ce répertoire contient toutes les données de sortie du solveur, y compris les constantes en fonction de la position et les grandeurs calculées dérivées de la dynamique des excitons. Par exemple, en double-cliquant sur exciton.csv, vous obtenez un graphique de la distribution des excitons dans l’épaisseur du dispositif (voir ??). De même, Gn.csv indique le taux de génération d’électrons en fonction de la position, tandis que Gp.csv affiche le taux de génération de trous (voir ??). Dans le cas d’un dispositif 1D simple, ces deux fichiers de taux de génération sont en réalité identiques.
exciton_output dans le répertoire de sortie de la simulation.
4. Équation et paramètres du transport des excitons
La distribution des excitons à l’intérieur du dispositif est régie par l’équation de transport des excitons :
\[ \frac{\partial X}{\partial t} = \nabla \!\cdot \!\big(D\,\nabla X\big) + G_{\mathrm{optical}} - k_{\mathrm{dis}}\,X - k_{\mathrm{FRET}}\,X - k_{\mathrm{PL}}\,X - \alpha\,X^{2} \]
Ici, \(X(\mathbf{r},t)\) représente la densité d’excitons (m\(^{-3}\)) ; \(D\) représente le coefficient de diffusion des excitons (m\(^2\)s\(^{-1}\)) ; \(G_{\mathrm{optical}}\) représente le taux de génération locale d’excitons, fourni par le modèle optique (proportionnel au nombre de photons absorbés) ; \(k_{\mathrm{dis}}\) représente le taux de dissociation en charges libres ; \(k_{\mathrm{FRET}}\) représente le taux de transfert d’énergie par résonance de Förster ; \(k_{\mathrm{PL}}\) représente le taux de désintégration radiative ; et \(\alpha\) représente le coefficient d’annihilation exciton-exciton (m\(^3\)s\(^{-1}\)). Lorsque ce modèle est activé, les équations de dérive-diffusion électroniques utilisent un terme de génération défini comme \(G = k_{\mathrm{dis}}\,X\), qui peut être limité aux régions interfaciales, si cela est configuré.
Tous ces paramètres peuvent être consultés et modifiés dans l’éditeur de paramètres électriques, accessible depuis le menu principal (voir ??). Les champs spécifiques aux excitons sont regroupés sous la rubrique « Excitons ». Ils sont résumés dans le tableau ci-dessous.
| Paramètre | Signification | Unités |
|---|---|---|
| Longueur de diffusion. | Longueur de diffusion efficace des excitons avant la diffusion. | m |
| Durée de vie | Durée moyenne de la persistance d’un exciton avant sa désintégration ou sa dissociation. | s |
| kPL | Taux de désintégration radiative (photoluminescence). | s⁻¹ |
| kFRET | Taux de transfert d’énergie par résonance de Förster. | s⁻¹ |
| kα | Coefficient d’annihilation exciton-exciton. | m³/s |
| kdis | Constante de vitesse de dissociation, qui permet de transformer les excitons en charges libres. | s⁻¹ |
5. Comment le solveur s’intègre-t-il au processus de simulation ?
Lorsque le solveur d’excitons est désactivé (en haut de ??), l’optique basée sur la matrice de transfert calcule le taux de génération de porteurs et le transmet directement au solveur de dérive-diffusion. Lorsque le solveur d’excitons est activé (en bas de ??), l’optique fournit plutôt un taux de génération d’excitons au solveur d’excitons. Ce dernier fait ensuite évoluer cette population, en simulant la diffusion, le transfert (FRET), la dissociation, la désintégration radiative et l’annihilation, en utilisant les paramètres définis dans l’éditeur de paramètres électriques, et fournit le taux final de génération de porteurs pour le solveur de dérive-diffusion. En résumé, le solveur d’excitons s’insère entre l’optique et le transport électrique, ce qui vous permet de modéliser la physique des excitons sans modifier la configuration optique ou les équations de dérive-diffusion.
La séquence des opérations dans la simulation basée sur les excitons est illustrée dans ??. Tout d’abord, le solveur optique est exécuté, calculant le profil d’absorption des photons, section par section, dans tout le dispositif. Ensuite, le solveur d’excitons est lancé, au cours duquel les excitons générés sont propagés, transférés ou dissociés jusqu’à ce que le solveur converge, généralement en quelques dizaines d’étapes. Enfin, le taux de génération des excitons est examiné en fonction de la profondeur du dispositif et de la longueur d’onde, tel que calculé par le modèle de la matrice de transfert. Ensemble, ces résultats constituent la chaîne de traitement qui relie le modèle optique au solveur d’excitons, et, en fin de compte, aux équations de dérive-diffusion.
💡 Tâches : Essayez ces modifications simples pour explorer le modèle d’exciton (indice : modifiez-les d’un ou deux ordres de grandeur afin d’observer un effet).
- Modifiez la durée de vie de l’exciton dans l’éditeur des paramètres électriques, puis relancez la simulation.
- Réduisez le taux de dissociation (kdis) et observez son effet sur la génération de porteurs.
- Augmentez le coefficient d’annihilation exciton-exciton (α) afin de simuler des interactions exciton-exciton plus intenses.
- Comparez la courbe J-V lorsque le solveur d’excitons est activé et lorsqu’il est désactivé.
✅ Résultats attendus
- Une durée de vie plus longue des excitons entraîne une densité d’excitons plus élevée et une augmentation de la génération de porteurs de charge.
- Une diminution des taux de dissociation entraîne une réduction du photCourant, car un plus petit nombre d’excitons sont convertis en porteurs libres.
- Des taux d’annihilation élevés élargissent ou réduisent le profil de densité des excitons, ce qui diminue l’efficacité.
- Les courbes caractéristiques courant-tension (J-V) obtenues avec le solveur d’excitons activé peuvent présenter une valeur de courant de court-circuit (JSC) plus faible que celle obtenue avec le modèle standard de dérive-diffusion, en particulier si les phénomènes de piégeage ou d’annihilation sont importants.
📘 Retour à introduction à l’électronique organique pour une vue d’ensemble plus large du sujet.