Reactive Force Field for P/Sn/I System: Atomistic Insight into the Early Stage of Black Phosphorus and Phosphorene Synthesis Process

Sintesi Tecnica: Campo di Forza Reattivo per il Sistema P/Sn/I e Approccio Atomistico alla Sintesi del Fosforo Nero

Problematica
Il fosforo nero (BP) e il suo controparte bidimensionale, il fosforene, sono materiali promettenti per l'elettronica di prossima generazione grazie all'elevata mobilità dei portatori e ai band gap sintonizzabili. Tuttavia, la sintesi di film di BP di alta qualità e di grande area tramite Trasporto Chimico di Vapore (CVT) rimane una sfida. Il processo, che tipicamente utilizza fosforo rosso con additivi di stagno (Sn) e iodio (I₂), non è pienamente compreso a livello atomistico. Persistono incertezze chiave riguardanti i ruoli specifici degli additivi, l'identità delle specie stabili in fase gassosa, i meccanismi di trasporto del fosforo e la natura degli intermedi di pre-nucleazione. L'osservazione sperimentale di queste dinamiche nelle prime fasi è difficile a causa delle alte temperature e degli ambienti reattivi, mentre la teoria del funzionale della densità (DFT) standard è limitata a sistemi piccoli e configurazioni statiche a 0K, fallendo nel catturare le reti di reazione dinamiche a temperatura finita necessarie per modellare i meccanismi di crescita.

Metodologia
Per affrontare queste limitazioni, gli autori hanno sviluppato un nuovo campo di forza reattivo (ReaxFF) specifico per il sistema P/Sn/I, denominato P/Sn/I-2512. Questo campo di forza è stato addestrato contro un ampio dataset di meccanica quantistica (QM) generato tramite DFT (funzionale PBE con base TZ2P e correzione di dispersione D3 di Grimme). Il set di addestramento includeva:

• Energie di dissociazione dei legami Sn–Sn, Sn–I, Sn–P, I–P e I–I.
• Profili dell'angolo di valenza e della torsione per varie configurazioni molecolari.
• Calori di formazione per specie chiave (SnI₂, SnI₄, PI₃).
• Equazioni di stato per lo stagno in fase condensata ($\alpha$-Sn e $\beta$-Sn) ed energie di formazione di cluster di stagno (Sn$_x$, $x=2-6$).

I parametri del campo di forza sono stati ottimizzati utilizzando la Strategia Evolutiva di Adattamento della Matrice di Covarianza (CMA-ES) all'interno dell'Amsterdam Modeling Suite. La validazione è stata eseguita confrontando i risultati di ReaxFF con i dati DFT per la dissociazione dei legami, le equazioni di stato e le energie di reazione, mostrando un forte accordo.

Successivamente, sono state condotte simulazioni di Dinamica Molecolare Reattiva (ReaxFF-MD) utilizzando LAMMPS per investigare i processi di pre-nucleazione nelle prime fasi a 923 K. Lo studio ha impiegato l'ensemble NVT con un termostato Nosé-Hoover. Le simulazioni sono state strutturate per isolare variabili specifiche:

1. Validazione del Sistema: Simulazioni MD di sistemi Sn/I₂ e P/I₂ a diverse densità (pressioni efficaci di 0,01, 0,1 e 1 atm) per verificare la formazione di specie chiave come SnI₂, SnI₄, PI₂ e PI₃.
2. Ruolo dello Iodio: Tre sistemi sono stati confrontati (P solo; P + Sn; P + Sn + I₂) per isolare l'effetto dello iodio sulla ricombinazione del fosforo.
3. Effetti di Densità: Il sistema P/Sn/I₂ è stato simulato sotto tre regimi di densità (pressioni efficaci di 5, 10 e 100 atm) per modellare la transizione da vapore diluito a fasi condensate, imitando le condizioni in cui la formazione di BP si verifica sperimentalmente.

Risultati Chiave

• Validazione del Campo di Forza: Il campo di forza P/Sn/I-2512 ha riprodotto accuratamente i dati di riferimento DFT per la dissociazione dei legami, gli angoli di valenza, le barriere torsionali e le equazioni di stato per i cristalli e i cluster di stagno. Ha predetto con successo la formazione di SnI₂, SnI₄ e PI₃ ad alte temperature.
• Dipendenza dalla Densità in Sn/I₂ e P/I₂: Nel sistema Sn/I₂, l'aumento della densità ha spostato la composizione dalle specie elementari ad ambienti di coordinazione complessi, promuovendo la formazione di SnI₄ e Sn₂I₄ (dimeri di SnI₂) insieme a SnI₂. Nel sistema P/I₂, densità più elevate hanno facilitato la formazione di PI₃ e specie polimeriche PI$_x$, in linea con le difficoltà sperimentali nell'isolare il PI₃ puro.
• Ruolo dello Iodio nella Ricombinazione: Le simulazioni hanno rivelato che, in assenza di iodio, gli atomi di fosforo rimanevano largamente atomici o formavano cluster minimi. L'aggiunta di I₂ ha aumentato significativamente il tasso di ricombinazione P–P, portando alla formazione di motivi binari P$_x$I$_y$ e composti ternari transitori Sn$_x$P$_y$I$_z$. Ciò conferma il ruolo dello iodio nel solubilizzare il fosforo e facilitarne il trasporto.
• Percorso di Pre-nucleazione ad Alta Densità: A basse densità (5 atm), il sistema consisteva principalmente di fosforo atomico e piccoli cluster (P₂, P₃). A densità medie (10 atm), iniziavano a formarsi cluster più grandi (P₄, P₅). Ad alta densità (100 atm), il sistema si è condensato rapidamente in cluster di ordine superiore (P$_x$ con x = 8–12) e grandi composti ternari Sn$_x$P$_y$I$_z$.
• Identificazione degli Intermedi: Alla densità più elevata, le simulazioni hanno prodotto strutture di fosforo tipo Hittorf ai bordi dei cluster Sn$_x$P$_y$I$_z$. Inoltre, sono stati osservati motivi esagonali simili al BP stabilizzati dallo iodio, identificati come potenziali semi di nucleazione.

Significato e Rivendicazioni
L'articolo sostiene di fornire un modello atomistico pratico per studiare le condizioni di crescita e gli effetti degli additivi nella sintesi del vapore di BP/fosforene. I contributi primari sono:

1. Sviluppo di un Campo di Forza Validato: Il campo di forza P/Sn/I-2512 offre uno strumento affidabile per simulare sistemi reattivi complessi coinvolgenti fosforo, stagno e iodio, superando i limiti di dimensione e scala temporale della DFT.
2. Intuizione Meccanicistica: Lo studio elucida che lo iodio e la densità del sistema sono i fattori controllanti per la ricombinazione del fosforo. Dimostra che lo iodio non è meramente un agente di trasporto, ma promuove attivamente la ricombinazione P–P attraverso la formazione di intermedi P$_x$I$_y$ e composti ternari Sn$_x$P$_y$I$_z$.
3. Modello di Pre-nucleazione: Gli autori propongono un percorso specifico di pre-nucleazione in cui le condizioni di alta densità, combinate con lo iodio, facilitano la condensazione del fosforo in intermedi di tipo Hittorf e esagoni simili al BP. Questi elementi sono suggeriti come precursori per la successiva crescita del fosforo nero durante la fase di raffreddamento.

Il lavoro non propone nuovi protocolli sperimentali, ma offre piuttosto un quadro teorico per interpretare le osservazioni sperimentali esistenti, specificamente la necessità di alta pressione (densità) e il ruolo critico dello iodio nella sintesi CVT del fosforo nero.