A Unified microscopic picture of cation and anion migration in MAPbI$_3$

Utilizzando simulazioni di dinamica molecolare con potenziali di rete neurale, questo studio rivela che la rapida migrazione ionica in MAPbI$_3$ è guidata da un moto collettivo concertato di interstiziali MA e difetti legati allo I dipendenti dalla carica, mentre le vacanze MA rimangono immobili, rivedendo così la comprensione convenzionale del trasporto ionico nelle perovskiti ibride.

L'essenziale

Immagina una cella solare realizzata con un cristallo speciale chiamato MAPbI3. Non pensare a questo cristallo come a un blocco rigido di pietra, ma come a una spugna morbida e comprimibile composta da minuscoli mattoncini. All'interno di questa spugna, ci sono due tipi principali di mattoncini: blocchi di metalli pesanti (Piombo e Iodio) e blocchi organici più leggeri, definiti "molecole" (chiamati MA, che sono come piccole molecole di metilammonio).

Il problema è che questa spugna non è perfetta. A volte mancano dei mattoncini (creando vacanze), e a volte mattoncini extra vengono schiacciati dove non dovrebbero esserci (creando interstiziali). Quando questi "difetti" iniziano a muoversi, possono causare il degrado della cella solare nel tempo.

Da molto tempo, gli scienziati hanno cercato di capire esattamente come si muovono questi difetti e a quale velocità. I numeri ottenuti dagli esperimenti erano estremamente variabili, come un gruppo di persone che indovina la velocità di un'auto e ottiene risposte che vanno dal "passo dell'uomo" al "supersonico".

Questo articolo utilizza una simulazione computerizzata super-intelligente (alimentata dall'Intelligenza Artificiale) per osservare questi difetti muoversi in tempo reale, come una telecamera ad alta velocità che riprende una pista da ballo. Ecco cosa hanno scoperto, spiegato in modo semplice:

1. Il "Fantasma" contro il "Carico Pesante"

In questo cristallo, i difetti di Iodio (gli ioni alogeni) sono come fantasmi. Sono leggeri e scattanti. Che manchi un atomo di Iodio (una vacanza) o che ce ne sia uno extra schiacciato dentro (un interstiziale), esso si muove con grande facilità. L'energia necessaria per metterli in movimento è molto bassa, come spingere un carrello della spesa su un pavimento liscio.

2. Il Ballerino Sorprendente (La Molecola MA)

La grande sorpresa di questo articolo riguarda le molecole MA. Queste sono molto più grandi e pesanti degli atomi di Iodio. Potresti aspettarti che siano lente, goffe e difficili da muovere, come cercare di spingere un pianoforte a coda attraverso una stanza.

• La Vecchia Credenza: Gli scienziati pensavano che queste grandi molecole fossero bloccate o si muovessero molto lentamente.
• La Nuova Scoperta: La simulazione ha mostrato che gli interstiziali MA (le molecole extra) sono in realtà veloci quanto i fantasmi di Iodio!

Com'è possibile?
L'articolo spiega che queste grandi molecole non si muovono da sole. Si muovono in un abbraccio di gruppo. Immagina tre persone su una pista da ballo. Invece che una persona cerchi di schiacciarsi passando, ruotano e si spostano tutte insieme in un movimento coordinato e "concertato". Uno fa un passo avanti, gli altri ruotano per fare spazio, e improvvisamente l'intero gruppo si è spostato. Questo lavoro di squadra permette alle pesanti molecole MA di muoversi quasi alla stessa velocità dei minuscoli atomi di Iodio.

3. Colui che Rimane Fermo

C'è un'eccezione: le Vacanze MA (buchi dove manca una molecola MA). La simulazione ha mostrato che questi buchi sono essenzialmente immobili. Anche quando la temperatura è stata alzata molto nella simulazione, questi buchi non si sono mossi. È come se il buco fosse incollato al pavimento. Questo suggerisce che se vedi il MA muoversi in una cella solare, è probabile che siano le molecole extra a muoversi, non gli spazi vuoti.

4. Perché i Numeri Erano Confusi

L'articolo suggerisce che il motivo per cui gli esperimenti passati hanno dato risposte così diverse (alcuni dicendo che è lento, altri che è veloce) è perché stavano misurando cose diverse.

• Il movimento veloce (barriera energetica di 0,15–0,20 eV) è ciò che accade all'interno profondo del cristallo (diffusione nel volume), ed è su questo che si è concentrato questo studio.
• Il movimento più lento riportato in altri studi potrebbe avvenire ai bordi dei grani cristallini o alle interfacce tra di essi, dove le cose si bloccano e si muovono in modo diverso.

Il Quadro Generale

Questo studio riscrive il manuale su come comprendiamo questi materiali. Ci dice che:

1. Il lavoro di squadra conta: Anche le grandi molecole pesanti possono muoversi velocemente se si muovono insieme in una danza coordinata.
2. La carica non conta molto: A differenza dei difetti di Iodio, la cui velocità cambia in base alla loro carica elettrica, le molecole MA si muovono alla stessa velocità sia che siano cariche che neutre.
3. Il "lento" MA è un mito: L'idea che la parte organica del cristallo sia un collo di bottiglia lento e pigro è sbagliata; è in realtà piuttosto agile quando si muove come squadra.

Comprendendo che questi difetti sono così mobili e si muovono in modi specifici, gli scienziati possono ora progettare meglio modi per "passivare" (tappare) questi difetti o impedirne il movimento, il che dovrebbe aiutare a rendere le celle solari e le luci molto più durature.

Sommario tecnico

Sintesi Tecnica: Un'immagine microscopica unificata della migrazione di cationi e anioni in MAPbI3

Enunciato del Problema
La passivazione dei difetti e la limitazione della mobilità dei difetti sono fondamentali per prolungare la vita utile dei dispositivi in perovskite di alogenuro. Sebbene studi sperimentali abbiano misurato estesamente il trasporto ionico nel materiale prototipale MAPbI$_3$ (ioduro di piombo e metilammonio), i risultati mostrano una significativa dispersione nelle barriere di migrazione riportate (da 0,10 a 0,94 eV). Inoltre, la natura specifica e la composizione dei difetti migranti — se siano vacanze, interstiziali o difetti estesi — rimangono difficili da identificare sperimentalmente. Gli approcci computazionali convenzionali, basati sul metodo della banda elastica spinta (NEB), incontrano difficoltà in questi materiali soffici a causa dell'esistenza di molteplici percorsi a bassa energia in competizione e della bassa simmetria del catione metilammonio (MA). Inoltre, le simulazioni standard di dinamica molecolare (MD) spesso mancano della precisione necessaria per un'estrazione affidabile delle barriere, a causa delle limitazioni dei campi di forza classici.

Metodologia
Per affrontare queste sfide, gli autori hanno impiegato simulazioni di dinamica molecolare (MD) guidate da campi di forza appresi tramite machine learning (MLFF), addestrati su dati della teoria del funzionale densità (DFT).

• Sviluppo del Campo di Forza: Lo studio ha utilizzato il modello Allegro, un potenziale di rete neurale (NNP) locale equivariante E(3) basato su grafi, implementato tramite il pacchetto NequIP.
• Dati di Addestramento: Gli NNP sono stati addestrati su dati DFT generati utilizzando il funzionale SCAN+rVV10 (inclusi gli interazioni di van der Waals) all'interno del pacchetto VASP. I dataset sono stati costruiti utilizzando un approccio di campionamento "on-the-fly" da brevi run di MD, garantendo l'inclusione di strutture diverse.
• Modellazione dei Difetti: Le simulazioni hanno modellato difetti puntuali in supercelle 2$\times$2$\times$2 (96 atomi) e run di produzione in supercelle 6$\times$6$\times$6 (2592 atomi). Lo studio ha considerato vacanze e interstiziali di ioduro (I) e MA nei loro stati di carica più prominenti in condizioni intrinseche: $I^-_I$, $V^+_I$, $I^0_{MA}$, $I^+_{MA}$, $V^0_{MA}$ e $V^-_{MA}$. NNP separati sono stati addestrati per difetti MA neutri e carichi selezionando dati DFT con livelli di Fermi vicini alle bande di valenza o di conduzione.
• Validazione: L'accuratezza degli NNP è stata validata confrontando le barriere di migrazione calcolate tramite NEB a immagine arrampicata (CI-NEB) con i benchmark DFT, mostrando deviazioni inferiori a 0,11 eV. I confronti delle forze hanno prodotto un coefficiente di determinazione ($R^2$) di 0,99.
• Protocollo di Simulazione: Sono state eseguite tre simulazioni MD indipendenti di 2 nanosecondi a cinque temperature comprese tra 500 K e 600 K. I coefficienti di diffusione sono stati estratti utilizzando la relazione di Einstein, trattando tutti gli atomi di una specifica specie per tenere conto dei processi di espulsione ("kick-out") piuttosto che della diffusione per tracciante.

Risultati Chiave

• Barriere di Migrazione e Mobilità: Le simulazioni hanno rivelato che gli interstiziali di ioduro ($I^-_I$) e le vacanze ($V^+_I$) sono altamente mobili, con energie di attivazione ($E_a$) di 0,18 eV e 0,15 eV, rispettivamente. Remarkabilmente, gli interstiziali MA ($I^0_{MA}$ e $I^+_{MA}$) mostrano una mobilità comparabile con energie di attivazione di 0,18 eV e 0,20 eV.
• Vacanze MA Immobile: In contrasto con gli interstiziali, le vacanze MA ($V^0_{MA}$ e $V^-_{MA}$) sono risultate immobili entro le scale temporali e di temperatura simulate (fino a 600 K). Ciò suggerisce un limite inferiore per la loro energia di attivazione di circa 0,4 eV.
• Meccanismi di Migrazione:
  • Difetti di Ioduro: La migrazione avviene tramite un meccanismo di "espulsione" ("kick-out") in cui un atomo di ioduro vicino ruota per riempire una vacanza, oppure un interstiziale salta tra i ponti Pb-I-Pb.
  • Interstiziali MA: Nonostante la loro natura molecolare, gli interstiziali MA migrano tramite un movimento concertato che coinvolge tre ioni MA vicini. Due ioni MA in una gabbia Pb-I e un terzo vicino ruotano per allineare i loro assi C-N, permettendo allo ione centrale di migrare verso una gabbia adiacente.
• Dipendenza dallo Stato di Carica: La diffusione dei difetti legati allo ioduro mostra una dipendenza dallo stato di carica, mentre la diffusione dei difetti legati al MA è virtualmente indipendente dallo stato di carica. Ciò è attribuito al fatto che gli interstiziali e le vacanze MA creano livelli di impurezza molto superficiali che sono ionizzati a temperatura ambiente, portando a una carica delocalizzata che non ostacola significativamente il movimento.
• Coerenza con gli Esperimenti: Le barriere di attivazione calcolate (0,15–0,20 eV) si allineano con l'estremità inferiore dell'ampio intervallo di valori riportati sperimentalmente.

Significato e Affermazioni
Il documento afferma di fornire un'immagine microscopica unificata della migrazione ionica in MAPbI$_3$ che riconcilia le discrepanze nelle energie di attivazione sperimentali. Gli autori sostengono che i valori più bassi osservati negli esperimenti (0,15–0,20 eV) corrispondono alla diffusione intrinseca di massa dominata dai difetti di ioduro e dagli interstiziali MA. I valori più alti riportati in altri studi riflettono probabilmente contributi aggiuntivi, come il trasporto tra grani o la formazione di difetti complessi, piuttosto che barriere di migrazione intrinseche di massa.

Un contributo primario è la revisione della comprensione convenzionale secondo cui la migrazione del MA è intrinsecamente lenta a causa delle sue dimensioni molecolari. Lo studio dimostra che gli interstiziali MA sono quasi tanto mobili quanto i difetti di ioduro grazie a un meccanismo di migrazione collettivo e concertato. Inoltre, il lavoro evidenzia che la specie del catione A (MA) non influisce significativamente sul movimento degli alogenuri, poiché la mobilità dei difetti di alogenuro in MAPbI$_3$ è comparabile a quella nelle perovskiti totalmente inorganiche CsPb(I$_x$Br$_{1-x}$)$_3$. Infine, lo studio sottolinea l'importanza del movimento molecolare collettivo nel permettere una rapida migrazione ionica nelle perovskiti ibride.