Impact of dynamic electrostatic disorder on hole mobility in rubrene: a nonadiabatic molecular dynamics investigation

Sintesi Tecnica: Impatto del Disordine Elettrostatico Dinamico sulla Mobilità delle Lacune nel Rubrene

Problematica
I cristalli molecolari organici ad alta mobilità, come il rubrene, sono fondamentali per l'elettronica organica, eppure ottenere una descrizione quantitativamente predittiva del loro trasporto di carica rimane una sfida. Il trasporto di carica in questi materiali avviene in un regime intermedio tra l'hopping di piccolo polarone e il trasporto a banda, caratterizzato da una delocalizzazione transitoria in cui la funzione d'onda del portatore di carica è modulata dalle fluttuazioni termiche. Sebbene la dinamica molecolare non adiabatica quantistico-classica mista (specificamente il Fragment Orbital-based Surface Hopping, o FOB-SH) offra una via promettente per simulare questo regime senza assumere un meccanismo di trasporto specifico, le precedenti simulazioni su larga scala di cristalli apolari come il rubrene hanno tipicamente trascurato il disordine elettrostatico dinamico. Questa omissione è stata guidata dal costo computazionale della valutazione delle interazioni coulombiche a lungo raggio (che tipicamente richiedono la sommatoria di Ewald) e dall'assunzione che tali interazioni siano trascurabili nei sistemi apolari. Tuttavia, trascurare queste interazioni rimuove i contributi elettrostatici all'energia di riorganizzazione e sopprime le fluttuazioni dell'energia di sito, portando potenzialmente a previsioni imprecise della delocalizzazione del portatore e della mobilità.

Metodologia
Questo studio integra il disordine elettrostatico dinamico nelle simulazioni FOB-SH del trasporto di lacune a temperatura ambiente nel rubrene. Gli autori impiegano il metodo di sommatoria elettrostatica nello spazio reale Damped Shifted-Force (DSF), combinato con uno schema di addizione-sottrazione efficiente, per superare i colli di bottiglia computazionali associati alla valutazione dei contributi elettrostatici per molteplici stati localizzati di carica.

I componenti metodologici chiave includono:

• Trattamento Elettrostatico: Il metodo DSF approssima le interazioni Coulomb a lungo raggio utilizzando un'interazione nello spazio reale smorzata (damped) e traslata a zero in un cutoff finito ($R_{cut} = 16$ Å). Uno schema di addizione-sottrazione decompone l'energia elettrostatica di uno stato di carica localizzata nell'energia del sistema neutro più termini di correzione, riducendo la scalabilità computazionale da $O(N_{mol})$ a quasi $O(1)$ rispetto a un singolo calcolo.
• Effetti di Polarizzazione: Per tenere conto della polarizzazione elettronica senza il costo di campi di forza polarizzabili espliciti, è stato applicato uno schema implicito basato sulla scala di carica (fattore $\gamma = 0,80$) alle cariche parziali atomiche.
• Protocollo di Simulazione: Le simulazioni FOB-SH sono state eseguite a 300 K utilizzando una supercella di $50 \times 13 \times 1$ unità cellulari (1300 molecole). Le simulazioni hanno propagato la funzione d'onda della lacuna in modo quantistico mentre evolvevano i nuclei in modo classico. I risultati sono stati confrontati con un modello di base "zero-carica" (dove le cariche parziali sono impostate a zero) e validati rispetto ai calcoli Smooth Particle-Mesh Ewald (SPME) per l'energia di riorganizzazione.
• Analisi: Lo studio ha quantificato l'energia di riorganizzazione ($\lambda$), il disordine dell'energia di sito ($\sigma_E$), le fluttuazioni dell'accoppiamento elettronico, il Rapporto di Partecipazione Inversa (IPR) per misurare la delocalizzazione della funzione d'onda e la mobilità del portatore di carica ($\mu$) derivata dallo spostamento quadratico medio (MSD).

Risultati Chiave

1. Energia di Riorganizzazione: L'inclusione delle interazioni elettrostatiche (tramite DSF o SPME) aumenta l'energia di riorganizzazione per il trasferimento di lacune di circa 20–40 meV (da una linea di base di 152 meV a ~175–191 meV) rispetto al modello a zero carica. L'analisi spettrale rivela che questo aumento deriva da un potenziamento dell'accoppiamento sia attraverso i moti reticolari a bassa frequenza (<200 cm$^{-1}$) che attraverso i modi intramolecolari ad alta frequenza. Il metodo DSF cattura la maggior parte di questo contributo elettrostatico, mostrando un buon accordo con i benchmark SPME.
2. Disordine dell'Energia di Sito: L'inclusione dell'elettrostatica aumenta significativamente le fluttuazioni della radice quadrata media delle energie di sito ($\sigma_E$), passando da 63,5 meV (zero-carica) a 69,2 meV (DSF). Questo aumento è quantitativamente coerente con l'aumento dell'energia di riorganizzazione, confermando che l'elettrostatica modifica principalmente il disordine energetico diagonale lasciando pressoché invariati gli accoppiamenti elettronici off-diagonal.
3. Delocalizzazione della Funzione d'Onda: Il potenziamento del disordine dell'energia di sito riduce l'estensione spaziale della funzione d'onda della lacuna. Il Rapporto di Partecipazione Inversa (IPR) medio diminuisce da 13 (zero-carica) a 9 (DSF), indicando un portatore di carica più localizzato.
4. Mobilità di Carica: La riduzione della delocalizzazione porta a una significativa diminuzione della mobilità delle lacune prevista. Lungo la direzione ad alta mobilità $a$, la mobilità scende da 34,8 cm$^2$V$^{-1}$s$^{-1}$ (zero-carica) a 21,1 cm$^2$V$^{-1}$s$^{-1}$ (DSF). Questo valore è in stretto accordo con il riferimento sperimentale di 20,2 cm$^2$V$^{-1}$s$^{-1}$. Inoltre, l'inclusione dell'elettrostatica migliora il rapporto di anisotropia previsto ($\mu_a/\mu_b$), avvicinandolo alle osservazioni sperimentali.

Significatività e Rivendicazioni
Il documento dimostra che il disordine elettrostatico dinamico è una componente essenziale per le simulazioni quantitativamente predittive del trasporto di carica nei semiconduttori organici molecolari, anche quelli composti da molecole apolari come il rubrene. Includendo esplicitamente queste interazioni tramite il metodo computazionalmente efficiente DSF, gli autori mostrano che i modelli precedenti che trascuravano l'elettrostatica sovrastimavano la delocalizzazione e la mobilità dei portatori. Lo studio stabilisce che le fluttuazioni dell'energia di sito elettrostatiche sono il meccanismo primario responsabile della localizzazione dinamica osservata negli alti cristalli organici ad alta mobilità, e la loro inclusione è necessaria per colmare il divario tra simulazioni teoriche e realtà sperimentale. Il lavoro valida il metodo DSF come un'alternativa scalabile e valida alla sommatoria di Ewald per la dinamica non adiabatica su larga scala in sistemi periodici.