Dissociation Line and Driving Force for Nucleation of the Multiple Occupied Hydrogen Hydrate from Computer Simulation

Problematica
Gli idrati di idrogeno (H₂) rappresentano una promettente alternativa per il trasporto e l'accumulo di gas, offrendo potenzialmente una riduzione dei costi delle materie prime rispetto agli idruri metallici, pur mantenendo una capacità volumetrica di stoccaggio comparabile. Tuttavia, l'applicazione pratica degli idrati di H₂ è ostacolata dalla mancanza di una comprensione esaustiva dei loro equilibri di fase e dei fattori specifici che governano l'occupazione delle gabbie. Poiché la determinazione sperimentale delle proprietà degli idrati è complessa a causa delle condizioni estreme di alta pressione e bassa temperatura richieste per la formazione, gli approcci teorici e di simulazione offrono una prospettiva molecolare. Un dibattito significativo persiste riguardo all'occupazione ottimale delle gabbie piccole (D) e grandi (H) nella struttura II (sII) dell'idrato di H₂. Studi precedenti hanno prodotto risultati contrastanti, con alcuni che suggeriscono che la doppia occupazione delle gabbie D sia possibile e altri che indicano come questa destabilizzi la struttura. Inoltre, i modelli di simulazione esistenti spesso faticano a prevedere simultaneamente l'accurata solubilità di H₂ in acqua e la temperatura di coesistenza a tre fasi ($T_3$) attraverso un intervallo di pressioni.

Metodologia
Gli autori hanno impiegato simulazioni di dinamica molecolare (MD) utilizzando il pacchetto software GROMACS per determinare la linea di dissociazione degli idrati di H₂. Le molecole d'acqua sono state modellate utilizzando il potenziale TIP4P/Ice, mentre le molecole di H₂ sono state descritte utilizzando una versione modificata del modello Silvera e Goldman (SG). Un aspetto critico della metodologia ha riguardato la modifica della regola di combinazione di Berthelot per le interazioni incrociate tra l'ossigeno dell'acqua e le molecole di H₂. Gli autori hanno introdotto un fattore di correzione dipendente dalla temperatura, $\chi$, e hanno proposto un nuovo parametro, $b_0$, per raffinare il potenziale di interazione.

Sono stati utilizzati due distinti metodi di simulazione per determinare la temperatura di coesistenza a tre fasi ($T_3$), dove coesistono l'idrato (H), la fase liquida acquosa (Lw) e la fase liquida di H₂ pura (LH₂):

1. Metodo della Solubilità: Questo approccio ha previsto la simulazione di due equilibri separati: la fase acquosa in contatto con H₂ pura (Lw–LH₂) e la fase acquosa in contatto con l'idrato di H₂ (H–Lw). Tracciando la solubilità di H₂ nella fase acquosa in funzione della temperatura per entrambi i sistemi a una pressione costante (185 MPa), l'intersezione delle due curve di solubilità ha fornito il valore di $T_3$. Questo metodo è stato applicato a varie occupazioni di idrato, inclusi le occupazioni singole e multiple delle gabbie D e H (specificamente le configurazioni 1-1, 1-2, 2-2, 1-3 e 1-4).
2. Tecnica della Coesistenza Diretta: In questo metodo, un singolo box di simulazione conteneva tutte e tre le fasi (H, Lw, LH₂) in contatto. Variando la temperatura a pressioni fisse (100, 185, 300 MPa) e monitorando l'energia potenziale del sistema, gli autori hanno determinato se la fase idrato cresceva o si scioglieva. La $T_3$ è stata identificata come l'intervallo di temperatura tra la temperatura più alta in cui l'idrato cresce e quella più bassa in cui si scioglie.

Inoltre, la forza motrice per la nucleazione ($\Delta\mu_N^{EC}$) è stata calcolata in funzione del sottoraffreddamento e dell'occupazione dell'idrato utilizzando approssimazioni della Teoria della Nucleazione Classica (CNT).

Risultati Chiave

• Effetto dell'Occupazione su $T_3$: Utilizzando il metodo della solubilità a 185 MPa, lo studio ha rilevato che l'occupazione dell'idrato ha un effetto trascurabile sulla temperatura di dissociazione. Tutte le occupazioni testate (1-1, 1-2, 2-2, 1-3 e 1-4) hanno prodotto valori di $T_3$ entro i margini di errore di circa 264–266 K.
• Stabilità Termodinamica e Nucleazione: Nonostante i simili valori di $T_3$, l'analisi della forza motrice per la nucleazione ha rivelato significative differenze termodinamiche. La configurazione con 1 molecola di H₂ nelle gabbie piccole D e 3 molecole di H₂ nelle gabbie grandi H (occupazione 1-3) ha esibito la $\Delta\mu_N^{EC}$ più negativa, indicando che è la configurazione termodinamicamente favorita. Al contrario, la doppia occupazione delle gabbie piccole D (2-2) è risultata la meno favorita a causa delle elevate forze repulsive, portando ai valori di $\Delta\mu_N^{EC}$ più positivi.
• Miglioramento del Modello: La modifica originale di Berthelot proposta da Michalis et al. (con $b_0=0$) ha predetto accuratamente la solubilità di H₂, ma ha leggermente sottostimato i valori sperimentali di $T_3$ di 5–7 K. Gli autori hanno proposto una nuova modifica introducendo un parametro $b_0 = 0.1$ nell'espressione del fattore $\chi$. Questa regolazione ha migliorato significativamente la previsione dei valori di $T_3$, portandoli in eccellente accordo con i dati sperimentali a 100, 185 e 300 MPa, sebbene abbia comportato una leggera sovrastima della solubilità di H₂.
• Costanti Reticolari: Lo studio ha osservato che la costante reticolare ($a$) dell'idrato sII aumenta con la temperatura e con l'occupazione. Notevolmente, l'occupazione 2-2 ha prodotto una costante reticolare maggiore rispetto all'occupazione 1-4, ciò che è attribuito alla forte repulsione tra due molecole di H₂ confinate nelle piccole gabbie D.

Significatività e Rivendicazioni
L'articolo sostiene di fornire la prima determinazione tramite simulazione al computer dell'effetto dell'occupazione dell'idrato di H₂ sia sulla temperatura di dissociazione che sulla forza motrice per la nucleazione. Gli autori sottolineano che il loro lavoro risolve l'ambiguità riguardante la stechiometria più stabile, identificando l'occupazione 1-3 come lo stato termodinamicamente preferito nelle condizioni studiate. Questa scoperta suggerisce una capacità di stoccaggio teorica del 3,2% in peso per gli idrati di H₂, che è inferiore alla capacità del 5,0% in peso associata alla doppia occupazione di tutte le gabbie, ma rappresenta la configurazione più stabile.

Lo studio evidenzia inoltre l'importanza della parametrizzazione del campo di forza, dimostrando che una semplice modifica alla regola di combinazione di Berthelot può riconciliare le discrepanze tra la solubilità simulata e le temperature di equilibrio di fase. Gli autori concludono che, sebbene l'occupazione 1-3 sia la più stabile, la scelta dei parametri di simulazione (specificamente il parametro $b_0$) consente ai ricercatori di dare priorità alla previsione accurata della solubilità o alla previsione accurata della linea di dissociazione, a seconda dei requisiti specifici della loro investigazione.