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

Enunciado do Problema
Os hidratos de H₂ representam uma alternativa promissora para o transporte e armazenamento de gás, oferecendo potencialmente uma redução nos custos de matéria-prima em comparação com os hidretos metálicos, mantendo uma capacidade de armazenamento volumétrico comparável. No entanto, a aplicação prática dos hidratos de H₂ é dificultada pela falta de uma compreensão abrangente de seus equilíbrios de fase e dos fatores específicos que regem a ocupação das gaiolas. A determinação experimental das propriedades dos hidratos é desafiadora devido às condições extremas de alta pressão e baixa temperatura necessárias para a formação, enquanto abordagens teóricas e de simulação oferecem uma perspectiva molecular. Um debate significativo persiste em relação à ocupação ideal das gaiolas pequenas (D) e grandes (H) na estrutura II (sII) do hidrato de H₂. Estudos anteriores produziram resultados conflitantes, com alguns sugerindo que a ocupação dupla das gaiolas D é possível e outros indicando que ela desestabiliza a estrutura. Além disso, os modelos de simulação existentes frequentemente têm dificuldade em prever simultaneamente a solubilidade precisa de H₂ em água e a temperatura de coexistência de três fases ($T_3$) através de uma gama de pressões.

Metodologia
Os autores empregaram simulações de dinâmica molecular (MD) utilizando o pacote de software GROMACS para determinar a linha de dissociação dos hidratos de H₂. As moléculas de água foram modeladas usando o potencial TIP4P/Ice, enquanto as moléculas de H₂ foram descritas usando uma versão modificada do modelo Silvera e Goldman (SG). Um aspecto crítico da metodologia envolveu a modificação da regra de combinação de Berthelot para as interações cruzadas entre o oxigênio da água e as moléculas de H₂. Os autores introduziram um fator de correção dependente da temperatura, $\chi$, e propuseram um novo parâmetro, $b_0$, para refinar o potencial de interação.

Dois métodos distintos de simulação foram utilizados para determinar a temperatura de coexistência de três fases ($T_3$), onde as fases do hidrato (H), líquido aquoso (Lw) e H₂ líquido (LH₂) coexistem:

1. Método de Solubilidade: Esta abordagem envolveu a simulação de dois equilíbrios separados: a fase aquosa em contato com H₂ puro (Lw–LH₂) e a fase aquosa em contato com o hidrato de H₂ (H–Lw). Ao plotar a solubilidade de H₂ na fase aquosa como uma função da temperatura para ambos os sistemas sob uma pressão constante (185 MPa), a interseção das duas curvas de solubilidade forneceu o valor de $T_3$. Este método foi aplicado a várias ocupações de hidrato, incluindo ocupações simples e múltiplas das gaiolas D e H (especificamente as configurações de ocupação 1-1, 1-2, 2-2, 1-3 e 1-4).
2. Técnica de Coexistência Direta: Neste método, uma única caixa de simulação continha todas as três fases (H, Lw, LH₂) em contato. Ao variar a temperatura em pressões fixas (100, 185 e 300 MPa) e monitorar a energia potencial do sistema, os autores determinaram se a fase de hidrato crescia ou derretia. O $T_3$ foi identificado como a faixa de temperatura entre a temperatura mais alta onde o hidrato cresce e a mais baixa onde ele derrete.

Adicionalmente, a força motriz para a nucleação ($\Delta\mu_N^{EC}$) foi calculada como uma função do superresfriamento e da ocupação do hidrato usando aproximações da Teoria da Nucleação Clássica (CNT).

Principais Resultados

• Efeito da Ocupação em $T_3$: Usando o método de solubilidade a 185 MPa, o estudo descobriu que a ocupação do hidrato tem um efeito negligenciável na temperatura de dissociação. Todas as ocupações testadas (1-1, 1-2, 2-2, 1-3 e 1-4) forneceram valores de $T_3$ dentro das margens de erro de aproximadamente 264–266 K.
• Estabilidade Termodinâmica e Nucleação: Apesar dos valores de $T_3$ semelhantes, a análise da força motriz para a nucleação revelou diferenças termodinâmicas significativas. A configuração com 1 molécula de H₂ nas gaiolas pequenas D e 3 moléculas de H₂ nas gaiolas grandes H (ocupação 1-3) exibiu o $\Delta\mu_N^{EC}$ mais negativo, indicando ser a configuração termodinamicamente mais favorável. Por outro lado, a ocupação dupla das gaiolas pequenas D (2-2) foi considerada a menos favorável devido às altas forças de repulsão, resultando nos valores de $\Delta\mu_N^{EC}$ mais positivos.
• Melhoria do Modelo: A modificação original de Berthelot proposta por Michalis et al. (com $b_0=0$) previu com precisão a solubilidade de H₂, mas subestimou ligeiramente os valores experimentais de $T_3$ em 5–7 K. Os autores propuseram uma nova modificação introduzindo um parâmetro $b_0 = 0,1$ na expressão do fator $\chi$. Este ajuste melhorou significativamente a previsão dos valores de $T_3$, trazendo-os para um excelente acordo com os dados experimentais em 100, 185 e 300 MPa, embora tenha resultado em uma ligeira superestimação da solubilidade de H₂.
• Constantes de Rede: O estudo observou que a constante de rede ($a$) do hidrato sII aumenta com a temperatura e a ocupação. Notavelmente, a ocupação 2-2 resultou em uma constante de rede maior do que a ocupação 1-4, atribuída à forte repulsão entre duas moléculas de H₂ confinadas nas pequenas gaiolas D.

Significância e Alegações
O artigo afirma fornecer a primeira determinação por simulação computacional do efeito da ocupação do hidrato de H₂ tanto na temperatura de dissociação quanto na força motriz para a nucleação. Os autores enfatizam que seu trabalho resolve a ambiguidade em relação à estequiometria mais estável, identificando a ocupação 1-3 como o estado termodinamicamente preferido sob as condições estudadas. Esta descoberta sugere uma capacidade de armazenamento teórica de 3,2 % em massa para hidratos de H₂, que é inferior à capacidade de 5,0 % em massa associada à ocupação dupla de todas as gaiolas, mas representa a configuração mais estável.

O estudo também destaca a importância da parametrização do campo de força, demonstrando que uma modificação simples na regra de combinação de Berthelot pode reconciliar as discrepâncias entre a solubilidade simulada e as temperaturas de equilíbrio de fase. Os autores concluem que, embora a ocupação 1-3 seja a mais estável, a escolha dos parâmetros de simulação (especificamente o parâmetro $b_0$) permite que os pesquisadores priorizem tanto a previsão precisa da solubilidade quanto a previsão precisa da linha de dissociação, dependendo dos requisitos específicos de sua investigação.