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

Planteamiento del Problema
Los hidratos de hidrógeno (H₂) representan una alternativa prometedora para el transporte y almacenamiento de gas, ofreciendo potencialmente una reducción en los costos de materia prima en comparación con los hidruros metálicos, manteniendo una capacidad de almacenamiento volumétrico comparable. Sin embargo, la aplicación práctica de los hidratos de H₂ se ve dificultada por la falta de una comprensión exhaustiva de sus equilibrios de fase y de los factores específicos que rigen la ocupación de las jaulas. Debido a que la determinación experimental de las propiedades de los hidratos es un desafío debido a las condiciones extremas de alta presión y baja temperatura requeridas para su formación, los enfoques teóricos y de simulación ofrecen una perspectiva molecular. Un debate significativo persiste con respecto a la ocupación óptima de las jaulas pequeñas (D) y grandes (H) en la estructura II (sII) del hidrato de H₂. Estudios previos han arrojado resultados contradictorios, donde algunos sugieren que la doble ocupación de las jaulas D es posible y otros indican que esto desestabiliza la estructura. Además, los modelos de simulación existentes a menudo luchan por predecir simultáneamente la solubilidad de H₂ en agua y la temperatura de coexistencia de tres fases ($T_3$) a través de un rango de presiones.

Metodología
Los autores emplearon simulaciones de dinámica molecular (MD) utilizando el paquete de software GROMACS para determinar la línea de disociación de los hidratos de H₂. Las moléculas de agua fueron modeladas utilizando el potencial TIP4P/Ice, mientras que las moléculas de H₂ fueron descritas mediante una versión modificada del modelo de Silvera y Goldman (SG). Un aspecto crítico de la metodología involucró la modificación de la regla de combinación de Berthelot para las interacciones cruzadas entre el oxígeno del agua y las moléculas de H₂. Los autores introdujeron un factor de corrección dependiente de la temperatura, $\chi$, y propusieron un nuevo parámetro, $b_0$, para refinar el potencial de interacción.

Se utilizaron dos métodos de simulación distintos para determinar la temperatura de coexistencia de tres fases ($T_3$), donde coexisten la fase de hidrato (H), la fase líquida acuosa (Lw) y la fase de líquido de H₂ puro (LH₂):

1. Método de Solubilidad: Este enfoque involucró la simulación de dos equilibrios separados: la fase acuosa en contacto con H₂ puro (Lw–LH₂) y la fase acuosa en contacto con el hidrato de H₂ (H–Lw). Al graficar la solubilidad de H₂ en la fase acuosa en función de la temperatura para ambos sistemas a una presión constante (185 MPa), la intersección de las dos curvas de solubilidad proporcionó el valor de $T_3$. Este método se aplicó a diversas ocupaciones de hidratos, incluyendo ocupaciones simples y múltiples de las jaulas D y H (específicamente las configuraciones 1-1, 1-2, 2-2, 1-3 y 1-4).
2. Técnica de Coexistencia Directa: En este método, un único cuadro de simulación contenía las tres fases (H, Lw, LH₂) en contacto. Al variar la temperatura a presiones fijas (100, 185 y 300 MPa) y monitorear la energía potencial del sistema, los autores determinaron si la fase de hidrato crecía o se derretía. El $T_3$ se identificó como el rango de temperatura entre la temperatura más alta donde el hidrato crece y la más baja donde se derrite.

Además, la fuerza impulsora para la nucleación ($\Delta\mu_N^{EC}$) se calculó como una función de la sobrefusión y la ocupación del hidrato utilizando aproximaciones de la Teoría de Nucleación Clásica (CNT).

Resultados Clave

• Efecto de la Ocupación en $T_3$: Utilizando el método de solubilidad a 185 MPa, el estudio encontró que la ocupación del hidrato tiene un efecto insignificante en la temperatura de disociación. Todas las ocupaciones probadas (1-1, 1-2, 2-2, 1-3 y 1-4) arrojaron valores de $T_3$ dentro de las barras de error de aproximadamente 264–266 K.
• Estabilidad Termodinámica y Nucleación: A pesar de los valores similares de $T_3$, el análisis de la fuerza impulsora para la nucleación reveló diferencias termodinámicas significativas. La configuración con 1 molécula de H₂ en las jaulas pequeñas D y 3 moléculas de H₂ en las jaulas grandes H (ocupación 1-3) exhibió el $\Delta\mu_N^{EC}$ más negativo, lo que indica que es la configuración termodinámicamente más favorecida. Por el contrario, la doble ocupación de las jaulas pequeñas D (2-2) resultó ser la menos favorecida debido a las altas fuerzas de repulsión, resultando en los valores de $\Delta\mu_N^{EC}$ más positivos.
• Mejora del Modelo: La modificación original de Berthelot propuesta por Michalis et al. (con $b_0=0$) predijo con precisión la solubilidad de H₂, pero subestimó ligeramente los valores experimentales de $T_3$ en un 5–7 K. Los autores propusieron una nueva modificación introduciendo un parámetro $b_0 = 0.1$ en la expresión del factor $\chi$. Este ajuste mejoró significativamente la predicción de los valores de $T_3$, logrando una excelente concordancia con los datos experimentales a 100, 185 y 300 MPa, aunque resultó en una ligera sobreestimación de la solubilidad de H₂.
• Constantes de Red: El estudio observó que la constante de red ($a$) del hidrato sII aumenta tanto con la temperatura como con la ocupación. Notablemente, la ocupación 2-2 resultó en una constante de red mayor que la ocupación 1-4, atribuido a la fuerte repulsión entre dos moléculas de H₂ confinadas en las pequeñas jaulas D.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera determinación mediante simulación computacional del efecto de la ocupación del hidrato de H₂ tanto en la temperatura de disociación como en la fuerza impulsora para la nucleación. Los autores enfatizan que su trabajo resuelve la ambigüedad respecto a la estequiometría más estable, identificando la ocupación 1-3 como el estado termodinámicamente preferido bajo las condiciones estudiadas. Este hallazgo sugiere una capacidad de almacenamiento teórica de 3.2 % en peso para los hidratos de H₂, la cual es inferior a la capacidad del 5.0 % asociada con la doble ocupación de todas las jaulas, pero representa la configuración más estable.

El estudio también destaca la importancia de la parametrización del campo de fuerza, demostrando que una modificación simple a la regla de combinación de Berthelot puede reconciliar las discrepancias entre la solubilidad simulada y las temperaturas de equilibrio de fase. Los autores concluyen que, si bien la ocupación 1-3 es la más estable, la elección de los parámetros de simulación (específicamente el parámetro $b_0$) permite a los investigadores priorizar ya sea la predicción precisa de la solubilidad o la predicción precisa de la línea de disociación, dependiendo de los requisitos específicos de su investigación.