Tutorial sul dominio degli eccitoni (parte A): esecuzione di una simulazione 3D del dominio degli eccitoni.Tutorial sul dominio degli eccitoni (parte A): esecuzione di una simulazione 3D del dominio degli eccitoni.
1. Introduzione1. Introduzione
Comprendere come vengono generati, trasportati e dissociati gli eccitoni è fondamentale per il funzionamento delle celle solari organiche a eterogiunzione sfusa (BHJ). In questi sistemi, la luce viene assorbita principalmente nella fase donatrice, generando coppie elettrone-lacuna legate, note come eccitoni. Questi eccitoni devono diffondersi attraverso il materiale donatore e raggiungere un’interfaccia donatore-accettore prima di decadere. All’interfaccia, possono dissociarsi in elettroni e lacune libere, con una probabilità determinata dalla morfologia locale, dalla dimensionalità e dalle costanti di velocità cinetiche rilevanti.Comprendere come vengono generati, trasportati e dissociati gli eccitoni è fondamentale per il funzionamento delle celle solari organiche a eterogiunzione sfusa (BHJ). In questi sistemi, la luce viene assorbita principalmente nella fase donatrice, generando coppie elettrone-lacuna legate, note come eccitoni. Questi eccitoni devono diffondersi attraverso il materiale donatore e raggiungere un’interfaccia donatore-accettore prima di decadere. All’interfaccia, possono dissociarsi in elettroni e lacune libere, con una probabilità determinata dalla morfologia locale, dalla dimensionalità e dalle costanti di velocità cinetiche rilevanti.
L’equilibrio tra la diffusione degli eccitoni, la dissociazione interfacciale e i processi di perdita concorrenti svolge quindi un ruolo determinante nella definizione dell’efficienza effettiva della fotogenerazione. Tuttavia, la modellizzazione esplicita degli eccitoni è stata storicamente evitata non a causa di difficoltà concettuali, ma a causa dell’incertezza dei parametri. Le lunghezze di diffusione degli eccitoni, i tempi di vita, i tassi di dissociazione interfacciale, i canali di decadimento radiativo e non radiativo e i processi di annichilimento erano spesso definiti in modo insufficiente per un determinato sistema di materiali. In molti flussi di lavoro di modellizzazione pratici, questi effetti sono quindi stati inclusi in un singolo fattore scalare di efficienza fotonica, \(\eta_{\mathrm{photon}}\), che rappresenta le perdite nette senza analizzare i meccanismi di trasporto e le dinamiche sottostanti.L’equilibrio tra la diffusione degli eccitoni, la dissociazione interfacciale e i processi di perdita concorrenti svolge quindi un ruolo determinante nella definizione dell’efficienza effettiva della fotogenerazione. Tuttavia, la modellizzazione esplicita degli eccitoni è stata storicamente evitata non a causa di difficoltà concettuali, ma a causa dell’incertezza dei parametri. Le lunghezze di diffusione degli eccitoni, i tempi di vita, i tassi di dissociazione interfacciale, i canali di decadimento radiativo e non radiativo e i processi di annichilimento erano spesso definiti in modo insufficiente per un determinato sistema di materiali. In molti flussi di lavoro di modellizzazione pratici, questi effetti sono quindi stati inclusi in un singolo fattore scalare di efficienza fotonica, \(\eta_{\mathrm{photon}}\), che rappresenta le perdite nette senza analizzare i meccanismi di trasporto e le dinamiche sottostanti.
La situazione sta cambiando. I progressi nelle tecniche sperimentali di caratterizzazione stanno iniziando a fornire misurazioni dirette dei tempi di vita degli eccitoni, delle lunghezze di diffusione e dei canali di perdita nei sistemi simili a BHJ. Man mano che questi parametri vengono definiti in modo più preciso, la modellizzazione esplicita dei domini degli eccitoni diventa sempre più utile, non tanto come strumento per prevedere l’efficienza assoluta dei dispositivi, quanto come metodo per analizzare come le dimensioni del dominio, la dimensionalità e la cinetica dei materiali determinano congiuntamente un’efficace resa nella generazione di carica. Il modello dei domini degli eccitoni utilizzato in questo tutorial è stato sviluppato e pubblicato nell’ambito di questo recente lavoro sperimentale e di modellizzazione (si veda, ad esempio, Nature Materials 21, 55–61 (2022)).La situazione sta cambiando. I progressi nelle tecniche sperimentali di caratterizzazione stanno iniziando a fornire misurazioni dirette dei tempi di vita degli eccitoni, delle lunghezze di diffusione e dei canali di perdita nei sistemi simili a BHJ. Man mano che questi parametri vengono definiti in modo più preciso, la modellizzazione esplicita dei domini degli eccitoni diventa sempre più utile, non tanto come strumento per prevedere l’efficienza assoluta dei dispositivi, quanto come metodo per analizzare come le dimensioni del dominio, la dimensionalità e la cinetica dei materiali determinano congiuntamente un’efficace resa nella generazione di carica. Il modello dei domini degli eccitoni utilizzato in questo tutorial è stato sviluppato e pubblicato nell’ambito di questo recente lavoro sperimentale e di modellizzazione (si veda, ad esempio, Nature Materials 21, 55–61 (2022)).
In questo tutorial, adottiamo quindi una geometria idealizzata, ma completamente tridimensionale, per la cella unitaria: un dominio donatore incorporato in una matrice accettore, inizialmente rappresentato come una sfera donatrice all’interno di un cubo accettore. Sebbene la geometria sia semplificata, il modello considera esplicitamente la generazione, la diffusione, la dissociazione interfacciale e i processi di perdita concorrenti degli eccitoni. Ciò lo rende uno strumento pratico e fisicamente trasparente per esplorare come i parametri misurati sperimentalmente interagiscono tra loro, per verificarne la coerenza interna e per sviluppare una comprensione di come la morfologia e la cinetica, insieme, controllano l’efficienza complessiva di fotogenerazione in tre dimensioni.In questo tutorial, adottiamo quindi una geometria idealizzata, ma completamente tridimensionale, per la cella unitaria: un dominio donatore incorporato in una matrice accettore, inizialmente rappresentato come una sfera donatrice all’interno di un cubo accettore. Sebbene la geometria sia semplificata, il modello considera esplicitamente la generazione, la diffusione, la dissociazione interfacciale e i processi di perdita concorrenti degli eccitoni. Ciò lo rende uno strumento pratico e fisicamente trasparente per esplorare come i parametri misurati sperimentalmente interagiscono tra loro, per verificarne la coerenza interna e per sviluppare una comprensione di come la morfologia e la cinetica, insieme, controllano l’efficienza complessiva di fotogenerazione in tre dimensioni.
2. Equazioni fondamentali2. Equazioni fondamentali
In questo tutorial, il trasporto degli eccitoni viene trattato esplicitamente risolvendo l’equazione di diffusione degli eccitoni in tre dimensioni, nell’intero dominio di simulazione. Quando il modello degli eccitoni è attivo, l’assorbimento ottico influenza direttamente la popolazione degli eccitoni, che quindi evolve secondo l’equazione seguente:In questo tutorial, il trasporto degli eccitoni viene trattato esplicitamente risolvendo l’equazione di diffusione degli eccitoni in tre dimensioni, nell’intero dominio di simulazione. Quando il modello degli eccitoni è attivo, l’assorbimento ottico influenza direttamente la popolazione degli eccitoni, che quindi evolve secondo l’equazione seguente:
\[ \frac{\partial X}{\partial t} = \nabla \cdot \left( D \nabla X \right) + G_{\mathrm{optical}} - k_{\mathrm{dis}} X - k_{\mathrm{FRET}} X - k_{\mathrm{PL}} X - \alpha X^2. \] \[ \frac{\partial X}{\partial t} = \nabla \cdot \left( D \nabla X \right) + G_{\mathrm{optical}} - k_{\mathrm{dis}} X - k_{\mathrm{FRET}} X - k_{\mathrm{PL}} X - \alpha X^2. \]
In questa equazione, \(X(\mathbf{r},t)\) rappresenta la densità degli eccitoni (\(\mathrm{m^{-3}}\)) e \(D\) rappresenta il coefficiente di diffusione degli eccitoni (\(\mathrm{m^2\,s^{-1}}\)). Il termine di sorgente \(G_{\mathrm{optical}}\) (\(\mathrm{m^{-3}\,s^{-1}}\)) indica la generazione di eccitoni dovuta all’assorbimento ottico. I termini rimanenti descrivono i processi concorrenti che eliminano gli eccitoni dalla popolazione:In questa equazione, \(X(\mathbf{r},t)\) rappresenta la densità degli eccitoni (\(\mathrm{m^{-3}}\)) e \(D\) rappresenta il coefficiente di diffusione degli eccitoni (\(\mathrm{m^2\,s^{-1}}\)). Il termine di sorgente \(G_{\mathrm{optical}}\) (\(\mathrm{m^{-3}\,s^{-1}}\)) indica la generazione di eccitoni dovuta all’assorbimento ottico. I termini rimanenti descrivono i processi concorrenti che eliminano gli eccitoni dalla popolazione:
- Dissociazione (\(k_{\mathrm{dis}} X\)): conversione degli eccitoni in elettroni liberi e lacune nelle interfacce donatore-accettore. Questo è il meccanismo principale che contribuisce alla generazione di carica.Dissociazione (\(k_{\mathrm{dis}} X\)): conversione degli eccitoni in elettroni liberi e lacune nelle interfacce donatore-accettore. Questo è il meccanismo principale che contribuisce alla generazione di carica.
- Trasferimento di tipo Förster (\(k_{\mathrm{FRET}} X\)): trasferimento di energia non radiativa degli eccitoni ai siti o ai domini adiacenti, che ridistribuisce l’energia degli eccitoni senza generare direttamente portatori di carica liberi.Trasferimento di tipo Förster (\(k_{\mathrm{FRET}} X\)): trasferimento di energia non radiativa degli eccitoni ai siti o ai domini adiacenti, che ridistribuisce l’energia degli eccitoni senza generare direttamente portatori di carica liberi.
- Perdita per emissione di radiazione (\(k_{\mathrm{PL}} X\)): ricombinazione degli eccitoni tramite emissione di fotoni (fotoluminescenza), che rappresenta un meccanismo di perdita intrinseco.Perdita per emissione di radiazione (\(k_{\mathrm{PL}} X\)): ricombinazione degli eccitoni tramite emissione di fotoni (fotoluminescenza), che rappresenta un meccanismo di perdita intrinseco.
- Annihilazione (\(\alpha X^2\)): annichilazione bimolecolare di eccitoni, che si verifica a elevate concentrazioni e diventa significativa in condizioni di forte eccitazione o in geometrie ristrette.Annihilazione (\(\alpha X^2\)): annichilazione bimolecolare di eccitoni, che si verifica a elevate concentrazioni e diventa significativa in condizioni di forte eccitazione o in geometrie ristrette.
Il coefficiente di diffusione \(D\) viene definito utilizzando la lunghezza di diffusione degli eccitoni \(L\) (m) e il tempo di vita degli eccitoni \(\tau\) (s) tramiteIl coefficiente di diffusione \(D\) viene definito utilizzando la lunghezza di diffusione degli eccitoni \(L\) (m) e il tempo di vita degli eccitoni \(\tau\) (s) tramite
\[ D = \frac{L^2}{\tau}. \] \[ D = \frac{L^2}{\tau}. \]
In questo tutorial eseguiamo una simulazione che considera esclusivamente gli eccitoni: il risolutore calcola la densità degli eccitoni tridimensionale \(X(\mathbf{r})\) e i relativi termini di reazione, ma non risolve il trasporto di elettroni e lacune. Il termine di dissociazione \(k_{\mathrm{dis}} X\) viene quindi interpretato come una velocità di generazione di carica potenziale, risolta spazialmente, che indica dove e con quale efficacia gli eccitoni si trasformerebbero in portatori di carica liberi alle interfacce donatore-accettore.In questo tutorial eseguiamo una simulazione che considera esclusivamente gli eccitoni: il risolutore calcola la densità degli eccitoni tridimensionale \(X(\mathbf{r})\) e i relativi termini di reazione, ma non risolve il trasporto di elettroni e lacune. Il termine di dissociazione \(k_{\mathrm{dis}} X\) viene quindi interpretato come una velocità di generazione di carica potenziale, risolta spazialmente, che indica dove e con quale efficacia gli eccitoni si trasformerebbero in portatori di carica liberi alle interfacce donatore-accettore.
3. Creare la simulazione del dominio degli eccitoni.3. Creare la simulazione del dominio degli eccitoni.
Aprire la finestra «Nuova simulazione» (Figura ??). È possibile accedervi facendo clic su «Nuova simulazione» nella finestra principale. Se si fa doppio clic su «Simulazioni di eccitoni», si visualizzerà la libreria degli eccitoni (Figura ??). Fare doppio clic su «Dominio eccitoni» per aprire il progetto di esempio.Aprire la finestra «Nuova simulazione» (Figura ??). È possibile accedervi facendo clic su «Nuova simulazione» nella finestra principale. Se si fa doppio clic su «Simulazioni di eccitoni», si visualizzerà la libreria degli eccitoni (Figura ??). Fare doppio clic su «Dominio eccitoni» per aprire il progetto di esempio.
4. Verificare la geometria e i parametri.4. Verificare la geometria e i parametri.
Dopo aver aperto l’esempio, la finestra principale visualizza una semplice scena tridimensionale: una sfera donatrice inserita all’interno di un contenitore accettore (Figura ??). I raggi verdi indicano l’illuminazione proveniente dall’alto. Per esaminare i parametri degli eccitoni, fare clic su «Parametri elettrici» nel pannello a sinistra (all’interno della scheda «Struttura del dispositivo»). Questo aprirà l’editor dei parametri elettrici. Se si scorre fino in fondo a questa finestra, si troverà una sezione intitolata «Eccitoni», in cui sono elencati separatamente i parametri specifici per ciascun oggetto nella scena (lo strato circostante e la sfera inserita), come mostrato nelle figure ?? e ??.Dopo aver aperto l’esempio, la finestra principale visualizza una semplice scena tridimensionale: una sfera donatrice inserita all’interno di un contenitore accettore (Figura ??). I raggi verdi indicano l’illuminazione proveniente dall’alto. Per esaminare i parametri degli eccitoni, fare clic su «Parametri elettrici» nel pannello a sinistra (all’interno della scheda «Struttura del dispositivo»). Questo aprirà l’editor dei parametri elettrici. Se si scorre fino in fondo a questa finestra, si troverà una sezione intitolata «Eccitoni», in cui sono elencati separatamente i parametri specifici per ciascun oggetto nella scena (lo strato circostante e la sfera inserita), come mostrato nelle figure ?? e ??.
La sezione relativa agli eccitoni definisce i parametri presenti nell’equazione di diffusione-reazione. Questi includono la lunghezza di diffusione degli eccitoni \(L\), che determina il coefficiente di diffusione tramite \(D = L^{2}/\tau\), e il tempo di vita degli eccitoni \(\tau\), che controlla la scala temporale complessiva per il trasporto in assenza di ulteriori canali di perdita. La ricombinazione radiativa è descritta dalla velocità di fotoluminescenza \(k_{\mathrm{PL}}\), mentre i processi di trasferimento di energia non radiativo sono descritti dalla velocità di tipo Förster \(k_{\mathrm{FRET}}\). A concentrazioni di eccitoni più elevate, l’annichilazione bimolecolare è inclusa tramite il coefficiente \(\alpha\), che corrisponde al termine \(-\alpha X^{2}\) nell’equazione principale. Infine, la velocità di dissociazione \(k_{\mathrm{dis}}\) specifica l’efficienza con cui gli eccitoni si trasformano in stati di trasferimento di carica alle interfacce donatore-accettore ed è il parametro chiave che controlla l’efficienza effettiva di generazione di carica in questo modello.La sezione relativa agli eccitoni definisce i parametri presenti nell’equazione di diffusione-reazione. Questi includono la lunghezza di diffusione degli eccitoni \(L\), che determina il coefficiente di diffusione tramite \(D = L^{2}/\tau\), e il tempo di vita degli eccitoni \(\tau\), che controlla la scala temporale complessiva per il trasporto in assenza di ulteriori canali di perdita. La ricombinazione radiativa è descritta dalla velocità di fotoluminescenza \(k_{\mathrm{PL}}\), mentre i processi di trasferimento di energia non radiativo sono descritti dalla velocità di tipo Förster \(k_{\mathrm{FRET}}\). A concentrazioni di eccitoni più elevate, l’annichilazione bimolecolare è inclusa tramite il coefficiente \(\alpha\), che corrisponde al termine \(-\alpha X^{2}\) nell’equazione principale. Infine, la velocità di dissociazione \(k_{\mathrm{dis}}\) specifica l’efficienza con cui gli eccitoni si trasformano in stati di trasferimento di carica alle interfacce donatore-accettore ed è il parametro chiave che controlla l’efficienza effettiva di generazione di carica in questo modello.
L’elemento concettuale fondamentale di questo esempio è che la sfera (donatore) e lo strato (matrice accettore) sfruttano principi fisici diversi relativi agli eccitoni. Nelle immagini riportate sopra, la sfera donatrice include un valore diverso da zero per \(k_{\mathrm{PL}}\) (un canale di perdita). La regione circostante è configurata con un valore molto elevato per \(k_{\mathrm{dis}}\), quindi, una volta che gli eccitoni raggiungono (o lasciano) tale regione, si dissociano rapidamente. Questo fa sì che il confine tra donatore e accettore si comporti come un efficiente punto di raccolta, producendo il caratteristico profilo con valori elevati al centro e valori ridotti in prossimità del confine.L’elemento concettuale fondamentale di questo esempio è che la sfera (donatore) e lo strato (matrice accettore) sfruttano principi fisici diversi relativi agli eccitoni. Nelle immagini riportate sopra, la sfera donatrice include un valore diverso da zero per \(k_{\mathrm{PL}}\) (un canale di perdita). La regione circostante è configurata con un valore molto elevato per \(k_{\mathrm{dis}}\), quindi, una volta che gli eccitoni raggiungono (o lasciano) tale regione, si dissociano rapidamente. Questo fa sì che il confine tra donatore e accettore si comporti come un efficiente punto di raccolta, producendo il caratteristico profilo con valori elevati al centro e valori ridotti in prossimità del confine.
5. Avvia la simulazione.5. Avvia la simulazione.
Fare clic sul triangolo blu «Avvia simulazione» per avviare il risolutore. L’output del terminale verrà visualizzato nella finestra principale, come mostrato nella Figura ??. Le prime righe forniscono informazioni generali sulla configurazione e sulla geometria della simulazione. Le informazioni diagnostiche principali vengono visualizzate una volta che il risolutore di eccitoni inizia a funzionare e hanno il formato Exciton f() = … , steps = …. Questa riga indica l’avanzamento del risolutore di eccitoni tridimensionale. La quantità f() è un residuo che misura quanto il campo di densità di eccitoni corrente si discosta dall’autocoerenza; man mano che il risolutore esegue iterazioni, questo valore dovrebbe diminuire. Il contatore steps indica il numero di iterazioni. Nell’esempio mostrato, il residuo diminuisce da un valore iniziale a 2.97 × 104 al passaggio 2, dopodiché il risolutore indica che l’equazione è stata risolta e termina.Fare clic sul triangolo blu «Avvia simulazione» per avviare il risolutore. L’output del terminale verrà visualizzato nella finestra principale, come mostrato nella Figura ??. Le prime righe forniscono informazioni generali sulla configurazione e sulla geometria della simulazione. Le informazioni diagnostiche principali vengono visualizzate una volta che il risolutore di eccitoni inizia a funzionare e hanno il formato Exciton f() = … , steps = …. Questa riga indica l’avanzamento del risolutore di eccitoni tridimensionale. La quantità f() è un residuo che misura quanto il campo di densità di eccitoni corrente si discosta dall’autocoerenza; man mano che il risolutore esegue iterazioni, questo valore dovrebbe diminuire. Il contatore steps indica il numero di iterazioni. Nell’esempio mostrato, il residuo diminuisce da un valore iniziale a 2.97 × 104 al passaggio 2, dopodiché il risolutore indica che l’equazione è stata risolta e termina.
Su un tipico computer portatile moderno, questo esempio dovrebbe essere completato in circa 5–10 secondi. Se il tempo di esecuzione si prolunga fino a diversi minuti, ciò indica solitamente che la griglia è eccessivamente dettagliata oppure che la geometria o i parametri sono stati modificati in modo da aumentare notevolmente il numero di punti della griglia. In tal caso, è consigliabile rivedere le impostazioni della griglia o ripetere i passaggi di configurazione precedenti prima di procedere.Su un tipico computer portatile moderno, questo esempio dovrebbe essere completato in circa 5–10 secondi. Se il tempo di esecuzione si prolunga fino a diversi minuti, ciò indica solitamente che la griglia è eccessivamente dettagliata oppure che la geometria o i parametri sono stati modificati in modo da aumentare notevolmente il numero di punti della griglia. In tal caso, è consigliabile rivedere le impostazioni della griglia o ripetere i passaggi di configurazione precedenti prima di procedere.
Exciton f() diminuisce man mano che il risolutore converge.Output del risolutore durante la risoluzione delle equazioni dell’eccitone in 3D. L’indicatore del residuo Exciton f() diminuisce man mano che il risolutore converge.
exciton_output/.La scheda «Output». Il risolutore delle equazioni relative agli eccitoni scrive i risultati in exciton_output/.
exciton_output/. I file sono in formato CSV e possono essere aperti con i visualizzatori di OghmaNano o con altri programmi.Contenuto di exciton_output/. I file sono in formato CSV e possono essere aperti con i visualizzatori di OghmaNano o con altri programmi.6. Visualizzazione dei risultati.6. Visualizzazione dei risultati.
Una volta completata l’esecuzione del programma, passare alla scheda «Output» (Figura ??). Gli elementi più rilevanti per questo tutorial sono la directory exciton_output/, che contiene i risultati relativi agli eccitoni risolti spazialmente, e il file exciton_sim_info.json, che riassume le statistiche globali sulla generazione e la perdita e che verrà utilizzato in seguito nella Parte C. Facendo doppio clic su exciton_output/ si visualizza il contenuto mostrato nella Figura ??.Una volta completata l’esecuzione del programma, passare alla scheda «Output» (Figura ??). Gli elementi più rilevanti per questo tutorial sono la directory exciton_output/, che contiene i risultati relativi agli eccitoni risolti spazialmente, e il file exciton_sim_info.json, che riassume le statistiche globali sulla generazione e la perdita e che verrà utilizzato in seguito nella Parte C. Facendo doppio clic su exciton_output/ si visualizza il contenuto mostrato nella Figura ??.
Aprire G.csv per visualizzare la velocità di generazione degli eccitoni prodotta dal modello ottico (Figura ??). Nella finestra del grafico 3D è possibile ruotare la scena con il mouse e utilizzare i cursori Z e Y nella parte inferiore per sezionare il volume. Questo è spesso il metodo più rapido per verificare che la geometria, la direzione dell’illuminazione e la griglia si comportino come previsto.Aprire G.csv per visualizzare la velocità di generazione degli eccitoni prodotta dal modello ottico (Figura ??). Nella finestra del grafico 3D è possibile ruotare la scena con il mouse e utilizzare i cursori Z e Y nella parte inferiore per sezionare il volume. Questo è spesso il metodo più rapido per verificare che la geometria, la direzione dell’illuminazione e la griglia si comportino come previsto.
Successivamente, aprire exciton.csv (Figura ??), che mostra la densità degli eccitoni allo stato stazionario. Si dovrebbe osservare la caratteristica impronta della diffusione verso un punto di accumulo interfacciale: un’elevata densità di eccitoni all’interno della sfera donatrice, con una diminuzione verso il confine donatore-accettore, dove la dissociazione è più intensa. Utilizzare i cursori di selezione per verificare che questa diminuzione segua la geometria dell’interfaccia e non sia semplicemente un artefatto dell’angolo di visualizzazione.Successivamente, aprire exciton.csv (Figura ??), che mostra la densità degli eccitoni allo stato stazionario. Si dovrebbe osservare la caratteristica impronta della diffusione verso un punto di accumulo interfacciale: un’elevata densità di eccitoni all’interno della sfera donatrice, con una diminuzione verso il confine donatore-accettore, dove la dissociazione è più intensa. Utilizzare i cursori di selezione per verificare che questa diminuzione segua la geometria dell’interfaccia e non sia semplicemente un artefatto dell’angolo di visualizzazione.
Infine, aprire Gn.csv (Figura ??), che mostra la velocità di generazione degli elettroni, espressa in funzione della posizione, derivante dalla dissociazione degli eccitoni. In questo esempio, si dovrebbe osservare un anello (o una zona) ben definito di generazione localizzata in prossimità dell’interfaccia donatore-accettore. Utilizzare i cursori Z e Y per spostare il piano di sezionamento all’interno del volume ed esplorare come questo profilo di dissociazione si modifica in funzione della posizione.Infine, aprire Gn.csv (Figura ??), che mostra la velocità di generazione degli elettroni, espressa in funzione della posizione, derivante dalla dissociazione degli eccitoni. In questo esempio, si dovrebbe osservare un anello (o una zona) ben definito di generazione localizzata in prossimità dell’interfaccia donatore-accettore. Utilizzare i cursori Z e Y per spostare il piano di sezionamento all’interno del volume ed esplorare come questo profilo di dissociazione si modifica in funzione della posizione.
G.csv: frequenza di generazione degli eccitoni \(G_{\mathrm{optical}}\) calcolata tramite il modello ottico (unità di misura \(\mathrm{m^{-3}\,s^{-1}}\)). Ruotare e sezionare (tramite i cursori Z/Y) per analizzare la localizzazione spaziale.G.csv: frequenza di generazione degli eccitoni \(G_{\mathrm{optical}}\) calcolata tramite il modello ottico (unità di misura \(\mathrm{m^{-3}\,s^{-1}}\)). Ruotare e sezionare (tramite i cursori Z/Y) per analizzare la localizzazione spaziale.
exciton.csv: densità degli eccitoni \(X\) (unità di misura \(\mathrm{m^{-3}}\)). Il profilo «alta densità al centro / bassa densità ai bordi» è il risultato previsto della diffusione verso un punto in cui avviene una rapida dissoluzione.exciton.csv: densità degli eccitoni \(X\) (unità di misura \(\mathrm{m^{-3}}\)). Il profilo «alta densità al centro / bassa densità ai bordi» è il risultato previsto della diffusione verso un punto in cui avviene una rapida dissoluzione.
Gn.csv: velocità di generazione di elettroni dovuta alla dissociazione degli eccitoni (unità di misura \(\mathrm{m^{-3}\,s^{-1}}\)). La dissociazione è concentrata all’interfaccia donatore/accettore, creando un caratteristico anello/strato nelle sezioni. Utilizzare i cursori Z/Y per esplorare come questa caratteristica varia in funzione della posizione.Gn.csv: velocità di generazione di elettroni dovuta alla dissociazione degli eccitoni (unità di misura \(\mathrm{m^{-3}\,s^{-1}}\)). La dissociazione è concentrata all’interfaccia donatore/accettore, creando un caratteristico anello/strato nelle sezioni. Utilizzare i cursori Z/Y per esplorare come questa caratteristica varia in funzione della posizione.7. Quali sono i file contenuti in exciton_output/?7. Quali sono i file contenuti in exciton_output/?
La directory exciton_output/ contiene due tipi distinti di file. Il primo gruppo è costituito da campi calcolati generati dal risolutore di eccitoni, come la densità degli eccitoni e i tassi di generazione derivati dalla dissociazione. Il secondo gruppo contiene mappe dei parametri: copie tridimensionali dei parametri degli eccitoni definiti nell’interfaccia grafica. Questi file vengono principalmente salvati per garantire la coerenza e facilitare il debug, consentendo di verificare che i parametri dei materiali previsti siano stati applicati correttamente nello spazio. Le tabelle 1 e 2 riassumono queste due categorie. Tutte le quantità sono espresse in unità del Sistema Internazionale. (Alcuni file potrebbero non essere salvati se i relativi sottomodelli sono disabilitati).La directory exciton_output/ contiene due tipi distinti di file. Il primo gruppo è costituito da campi calcolati generati dal risolutore di eccitoni, come la densità degli eccitoni e i tassi di generazione derivati dalla dissociazione. Il secondo gruppo contiene mappe dei parametri: copie tridimensionali dei parametri degli eccitoni definiti nell’interfaccia grafica. Questi file vengono principalmente salvati per garantire la coerenza e facilitare il debug, consentendo di verificare che i parametri dei materiali previsti siano stati applicati correttamente nello spazio. Le tabelle 1 e 2 riassumono queste due categorie. Tutte le quantità sono espresse in unità del Sistema Internazionale. (Alcuni file potrebbero non essere salvati se i relativi sottomodelli sono disabilitati).
| Nome del fileNome del file | DescrizioneDescrizione | Unità tipicheUnità tipiche |
|---|---|---|
exciton.csvexciton.csv |
Campo di densità degli eccitoni \(X(\mathbf{r})\)Campo di densità degli eccitoni \(X(\mathbf{r})\) | \(\mathrm{m^{-3}}\)\(\mathrm{m^{-3}}\) |
G.csvG.csv |
Frequenza di generazione degli eccitoni \(G_{\mathrm{optical}}(\mathbf{r})\)Frequenza di generazione degli eccitoni \(G_{\mathrm{optical}}(\mathbf{r})\) | \(\mathrm{m^{-3}\,s^{-1}}\)\(\mathrm{m^{-3}\,s^{-1}}\) |
Gn.csvGn.csv |
Frequenza di generazione di elettroni dovuta alla dissociazione degli eccitoniFrequenza di generazione di elettroni dovuta alla dissociazione degli eccitoni | \(\mathrm{m^{-3}\,s^{-1}}\)\(\mathrm{m^{-3}\,s^{-1}}\) |
Gp.csvGp.csv |
Frequenza di generazione di lacune dovuta alla dissociazione degli eccitoniFrequenza di generazione di lacune dovuta alla dissociazione degli eccitoni | \(\mathrm{m^{-3}\,s^{-1}}\)\(\mathrm{m^{-3}\,s^{-1}}\) |
DD |
Coefficiente di diffusione degli eccitoni \(D\)Coefficiente di diffusione degli eccitoni \(D\) | \(\mathrm{m^{2}\,s^{-1}}\)\(\mathrm{m^{2}\,s^{-1}}\) |
alphaalpha |
Contributo dovuto all’annichilazione eccitone-eccitone \(\alpha X^2\)Contributo dovuto all’annichilazione eccitone-eccitone \(\alpha X^2\) | \(\mathrm{m^{3}\,s^{-1}}\)\(\mathrm{m^{3}\,s^{-1}}\) |
| Nome del fileNome del file | DescrizioneDescrizione | Unità tipicheUnità tipiche |
|---|---|---|
k_dis.csvk_dis.csv |
Velocità di dissociazione \(k_{\mathrm{dis}}(\mathbf{r})\)Velocità di dissociazione \(k_{\mathrm{dis}}(\mathbf{r})\) | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) |
k_fretk_fret |
Tasso di trasferimento di tipo Förster \(k_{\mathrm{FRET}}(\mathbf{r})\)Tasso di trasferimento di tipo Förster \(k_{\mathrm{FRET}}(\mathbf{r})\) | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) |
k_plk_pl |
Velocità di perdita per emissione radiativa (fotoluminescenza) \(k_{\mathrm{PL}}(\mathbf{r})\)Velocità di perdita per emissione radiativa (fotoluminescenza) \(k_{\mathrm{PL}}(\mathbf{r})\) | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) |
LL |
Lunghezza di diffusione (dispersione) degli eccitoni \(L(\mathbf{r})\)Lunghezza di diffusione (dispersione) degli eccitoni \(L(\mathbf{r})\) | \(\mathrm{m}\)\(\mathrm{m}\) |
tautau |
Tempo di vita degli eccitoni \(\tau(\mathbf{r})\)Tempo di vita degli eccitoni \(\tau(\mathbf{r})\) | \(\mathrm{s}\)\(\mathrm{s}\) |
I file contenenti i parametri elencati nella Tabella 2 non rappresentano i risultati ottenuti dal risolutore nel senso tradizionale del termine. Essi si limitano a riflettere i parametri degli eccitoni definiti nella sezione «Eccitoni» dell’interfaccia grafica, e vengono scritti sulla stessa griglia tridimensionale utilizzata per i campi di soluzione. Il loro scopo è fornire un metodo trasparente per verificare che i parametri siano assegnati agli oggetti e alle regioni corrette del dominio, in particolare quando si lavora con geometrie complesse.I file contenenti i parametri elencati nella Tabella 2 non rappresentano i risultati ottenuti dal risolutore nel senso tradizionale del termine. Essi si limitano a riflettere i parametri degli eccitoni definiti nella sezione «Eccitoni» dell’interfaccia grafica, e vengono scritti sulla stessa griglia tridimensionale utilizzata per i campi di soluzione. Il loro scopo è fornire un metodo trasparente per verificare che i parametri siano assegnati agli oggetti e alle regioni corrette del dominio, in particolare quando si lavora con geometrie complesse.
8.1 Interpretazione di exciton_sim_info.json (efficienza di dissociazione)8.1 Interpretazione di exciton_sim_info.json (efficienza di dissociazione)
exciton_sim_info.json contiene i dati relativi ai canali di generazione e di perdita, integrati a livello spaziale, per ciascun oggetto, nonché l’efficienza complessiva del processo di dissociazione.Il file di riepilogo exciton_sim_info.json contiene i dati relativi ai canali di generazione e di perdita, integrati a livello spaziale, per ciascun oggetto, nonché l’efficienza complessiva del processo di dissociazione.Sebbene le rappresentazioni grafiche tridimensionali siano estremamente utili per visualizzare dove vengono generati, persi o dissociati gli eccitoni, da sole non forniscono riepiloghi numerici esaustivi. Il file exciton_sim_info.json svolge questa funzione, fornendo dati sull’andamento spaziale dei tassi di generazione e di perdita per ciascun oggetto presente nella simulazione (lo strato circostante, la sfera incorporata e qualsiasi altra forma aggiunta).Sebbene le rappresentazioni grafiche tridimensionali siano estremamente utili per visualizzare dove vengono generati, persi o dissociati gli eccitoni, da sole non forniscono riepiloghi numerici esaustivi. Il file exciton_sim_info.json svolge questa funzione, fornendo dati sull’andamento spaziale dei tassi di generazione e di perdita per ciascun oggetto presente nella simulazione (lo strato circostante, la sfera incorporata e qualsiasi altra forma aggiunta).
Per visualizzare queste informazioni, aprire la scheda principale denominata «Output» e fare doppio clic su exciton_sim_info.json. Il file contiene una struttura JSON nidificata che elenca, per ciascun oggetto, il tasso totale di generazione di eccitoni e i tassi totali associati a ciascun meccanismo di perdita o dissociazione. Questi valori consentono di quantificare come le variazioni dei parametri del materiale, della geometria o della morfologia si traducano in variazioni dell’efficienza di generazione di carica.Per visualizzare queste informazioni, aprire la scheda principale denominata «Output» e fare doppio clic su exciton_sim_info.json. Il file contiene una struttura JSON nidificata che elenca, per ciascun oggetto, il tasso totale di generazione di eccitoni e i tassi totali associati a ciascun meccanismo di perdita o dissociazione. Questi valori consentono di quantificare come le variazioni dei parametri del materiale, della geometria o della morfologia si traducano in variazioni dell’efficienza di generazione di carica.
💡 Unità di misura: le quantità indicate in questo file rappresentano i totali calcolati integrando i valori su tutto il volume. A livello locale, i tassi di generazione e di reazione sono espressi come campi con unità di misura quali \(\mathrm{m^{-3}\,s^{-1}}\) (o \(\mathrm{m^{-3}}\) per le densità). Dopo l’integrazione sul volume, i totali corrispondenti sono espressi con unità di misura \(\mathrm{s^{-1}}\), che rappresentano il numero di eventi al secondo.💡 Unità di misura: le quantità indicate in questo file rappresentano i totali calcolati integrando i valori su tutto il volume. A livello locale, i tassi di generazione e di reazione sono espressi come campi con unità di misura quali \(\mathrm{m^{-3}\,s^{-1}}\) (o \(\mathrm{m^{-3}}\) per le densità). Dopo l’integrazione sul volume, i totali corrispondenti sono espressi con unità di misura \(\mathrm{s^{-1}}\), che rappresentano il numero di eventi al secondo.
La chiave di livello superiore segment0 indica semplicemente l’inizio del rapporto. Al suo interno, il campo name identifica l’oggetto contenitore (in questo caso, il livello), mentre il blocco nidificato shape0 corrisponde all’oggetto incorporato (la sfera). Ogni blocco riporta il contributo totale di tale oggetto ai processi di generazione, dissociazione e perdita.La chiave di livello superiore segment0 indica semplicemente l’inizio del rapporto. Al suo interno, il campo name identifica l’oggetto contenitore (in questo caso, il livello), mentre il blocco nidificato shape0 corrisponde all’oggetto incorporato (la sfera). Ogni blocco riporta il contributo totale di tale oggetto ai processi di generazione, dissociazione e perdita.
Nella tabella 2 sono riassunti i campi che si trovano in questo file e viene spiegato come interpretarli.Nella tabella 2 sono riassunti i campi che si trovano in questo file e viene spiegato come interpretarli.
| Campo JSONCampo JSON | SignificatoSignificato | Unità tipicheUnità tipiche | Dove compare.Dove compare. |
|---|---|---|---|
segment0segment0 |
Nodo contenitore per il report (non si tratta di una grandezza fisica).Nodo contenitore per il report (non si tratta di una grandezza fisica). | —— | Livello superioreLivello superiore |
namename |
Nome dell’oggetto (livello o forma incorporata)Nome dell’oggetto (livello o forma incorporata) | —— | All’interno di segment0 e shape0All’interno di segment0 e shape0 |
tot_Gtot_G |
Tasso totale di generazione di eccitoni, calcolato sull’intero volume dell’oggetto.Tasso totale di generazione di eccitoni, calcolato sull’intero volume dell’oggetto. | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) | Disponi e modella i blocchi.Disponi e modella i blocchi. |
tot_Rk_pltot_Rk_pl |
Tasso totale di perdita per radiazione (fotoluminescenza)Tasso totale di perdita per radiazione (fotoluminescenza) | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) | Disponi e modella i blocchi.Disponi e modella i blocchi. |
tot_Rk_frettot_Rk_fret |
Perdita totale di trasferimento o tasso di interazione di tipo FörsterPerdita totale di trasferimento o tasso di interazione di tipo Förster | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) | Disponi e modella i blocchi.Disponi e modella i blocchi. |
tot_Ralphatot_Ralpha |
Tasso totale di perdita per annichilazione eccitone-eccitoneTasso totale di perdita per annichilazione eccitone-eccitone | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) | Disponi e modella i blocchi.Disponi e modella i blocchi. |
tot_Rk_distot_Rk_dis |
Tasso totale di dissociazione degli eccitoni in portatori di carica liberiTasso totale di dissociazione degli eccitoni in portatori di carica liberi | \(\mathrm{s^{-1}}\)\(\mathrm{s^{-1}}\) | Disponi e modella i blocchi.Disponi e modella i blocchi. |
dis_effdis_eff |
Efficienza complessiva di dissociazione (frazione degli eccitoni generati che si dissociano)Efficienza complessiva di dissociazione (frazione degli eccitoni generati che si dissociano) | \(\%\)\(\%\) | Livello superioreLivello superiore |
Considerati nel loro insieme, i grafici tridimensionali e i riepiloghi numerici contenuti in exciton_sim_info.json offrono una visione complementare della simulazione: i grafici mostrano dove avvengono i processi, mentre il file JSON quantifica il contributo di ciascun processo. Questo lo rende uno strumento utile per analizzare in modo sistematico come le modifiche ai parametri o alla geometria influiscono sulla dissociazione degli eccitoni e sulla generazione effettiva di carica.Considerati nel loro insieme, i grafici tridimensionali e i riepiloghi numerici contenuti in exciton_sim_info.json offrono una visione complementare della simulazione: i grafici mostrano dove avvengono i processi, mentre il file JSON quantifica il contributo di ciascun processo. Questo lo rende uno strumento utile per analizzare in modo sistematico come le modifiche ai parametri o alla geometria influiscono sulla dissociazione degli eccitoni e sulla generazione effettiva di carica.
9. Assemblaggio.9. Assemblaggio.
Considerati nel loro insieme, i risultati tridimensionali descrivono una storia fisica coerente. L’assorbimento ottico genera eccitoni principalmente all’interno del dominio donatore. Questi eccitoni si diffondono quindi attraverso il donatore, si esauriscono in prossimità dell’interfaccia donatore-accettore e si dissociano in modo efficiente una volta raggiunta la regione accettore circostante. Le rappresentazioni spaziali mostrano dove avvengono questi processi; il rapporto riassuntivo exciton_sim_info.json quantifica il contributo di ciascun processo.Considerati nel loro insieme, i risultati tridimensionali descrivono una storia fisica coerente. L’assorbimento ottico genera eccitoni principalmente all’interno del dominio donatore. Questi eccitoni si diffondono quindi attraverso il donatore, si esauriscono in prossimità dell’interfaccia donatore-accettore e si dissociano in modo efficiente una volta raggiunta la regione accettore circostante. Le rappresentazioni spaziali mostrano dove avvengono questi processi; il rapporto riassuntivo exciton_sim_info.json quantifica il contributo di ciascun processo.
In questo esempio, il rapporto riassuntivo suddivide chiaramente i principi fisici di base in base alle diverse aree spaziali, consentendo di interpretare in modo trasparente i processi di generazione, perdita e dissociazione.In questo esempio, il rapporto riassuntivo suddivide chiaramente i principi fisici di base in base alle diverse aree spaziali, consentendo di interpretare in modo trasparente i processi di generazione, perdita e dissociazione.
- Sfera (donatore): la maggior parte dell’assorbimento ottico – e quindi la maggior parte della generazione di eccitoni – avviene all’interno della sfera donatrice, come dimostra l’elevato valore di
tot_Gnel bloccoshape0/sphere. All’interno del donatore, il principale meccanismo di perdita concorrente in questo esempio è la decadenza radiativa, quantificata datot_Rk_pl. Il rapporto tratot_Rk_pletot_Gfornisce una misura diretta dell’intensità con cui le perdite radiative sopprimono la popolazione di eccitoni prima che questa possa raggiungere l’interfaccia donatore-accettore.Sfera (donatore): la maggior parte dell’assorbimento ottico – e quindi la maggior parte della generazione di eccitoni – avviene all’interno della sfera donatrice, come dimostra l’elevato valore ditot_Gnel bloccoshape0/sphere. All’interno del donatore, il principale meccanismo di perdita concorrente in questo esempio è la decadenza radiativa, quantificata datot_Rk_pl. Il rapporto tratot_Rk_pletot_Gfornisce una misura diretta dell’intensità con cui le perdite radiative sopprimono la popolazione di eccitoni prima che questa possa raggiungere l’interfaccia donatore-accettore. - Strato (matrice accettore): nello specifico, lo strato circostante assorbe solo debolmente, quindi il suo contributo alla generazione di eccitoni (
tot_G) è limitato. Il suo ruolo principale è invece la conversione degli eccitoni: lo strato è configurato con un forte canale di dissociazione, il che fa sì chetot_Rk_dissia il termine dominante nel bloccolayer. Questo corrisponde direttamente alla dissociazione interfacciale osservata nelle rappresentazioni spaziali, doveGn.csvmostra una regione a forma di anello o guscio in cui si concentra la generazione di elettroni al confine tra il donatore e l’accettore.Strato (matrice accettore): nello specifico, lo strato circostante assorbe solo debolmente, quindi il suo contributo alla generazione di eccitoni (tot_G) è limitato. Il suo ruolo principale è invece la conversione degli eccitoni: lo strato è configurato con un forte canale di dissociazione, il che fa sì chetot_Rk_dissia il termine dominante nel bloccolayer. Questo corrisponde direttamente alla dissociazione interfacciale osservata nelle rappresentazioni spaziali, doveGn.csvmostra una regione a forma di anello o guscio in cui si concentra la generazione di elettroni al confine tra il donatore e l’accettore.
L’indicatore chiave dis_eff fornisce il valore complessivo dell’efficienza di dissociazione del sistema (in questo caso, approssimativamente \(89\%\)). A livello fisico, ciò indica che la maggior parte degli eccitoni generati nel donatore si diffonde verso l’interfaccia e si dissocia, mentre solo una piccola frazione viene persa a causa della decadenza radiativa. Altri meccanismi, come il trasferimento di Förster o l’annichilazione eccitone-eccitone, sono inattivi in questo specifico esempio.L’indicatore chiave dis_eff fornisce il valore complessivo dell’efficienza di dissociazione del sistema (in questo caso, approssimativamente \(89\%\)). A livello fisico, ciò indica che la maggior parte degli eccitoni generati nel donatore si diffonde verso l’interfaccia e si dissocia, mentre solo una piccola frazione viene persa a causa della decadenza radiativa. Altri meccanismi, come il trasferimento di Förster o l’annichilazione eccitone-eccitone, sono inattivi in questo specifico esempio.
In altre parole, la simulazione si svolge in un regime in cui la dissociazione interfacciale prevale sulla decadenza degli eccitoni, con conseguente elevata resa effettiva di fotoni convertiti in portatori di carica liberi. Questa combinazione di visualizzazione tridimensionale e analisi quantitativa fornisce una rappresentazione concisa ma completa dal punto di vista fisico di come la geometria e i parametri cinetici, interagendo tra loro, determinano l’efficienza di generazione di carica in un modello a dominio di eccitoni.In altre parole, la simulazione si svolge in un regime in cui la dissociazione interfacciale prevale sulla decadenza degli eccitoni, con conseguente elevata resa effettiva di fotoni convertiti in portatori di carica liberi. Questa combinazione di visualizzazione tridimensionale e analisi quantitativa fornisce una rappresentazione concisa ma completa dal punto di vista fisico di come la geometria e i parametri cinetici, interagendo tra loro, determinano l’efficienza di generazione di carica in un modello a dominio di eccitoni.
👉 Fase successiva: proseguire con Parte B, dove modificheremo la geometria del dominio (da sfera a una griglia arbitraria) e regoleremo le impostazioni di generazione ottica per osservare come la geometria e l’assorbimento interagiscono con la diffusione e la dissociazione.👉 Fase successiva: proseguire con Parte B, dove modificheremo la geometria del dominio (da sfera a una griglia arbitraria) e regoleremo le impostazioni di generazione ottica per osservare come la geometria e l’assorbimento interagiscono con la diffusione e la dissociazione.
📘 Torna a introduzione all’elettronica organica per una panoramica più completa dell’argomento. 📘 Torna a introduzione all’elettronica organica per una panoramica più completa dell’argomento.