[论文解读] Simulation of the carbon dioxide hydrate-water interfacial energy

Jesús Algabaa Esteban Acuña, José Manuel Míguez+3 more|arXiv (Cornell University)|Jan 23, 2026

Methane Hydrates and Related Phenomena被引用 0

一句话总结

论文在共存条件下使用 Mold Integration，结合 TIP4P/Ice 水和 TraPPE CO2，计算 CO2 水合物–水界面自由能，并与热力学定义及实验数据一致性进行验证。

ABSTRACT

Carbon dioxide hydrates are ice-like nonstoichiometric inclusion solid compounds with importance to global climate change, and gas transportation and storage. The thermodynamic and kinetic mechanisms that control carbon dioxide nucleation critically depend on hydrate-water interfacial free energy. Interfacial energies show large uncertainties due to the conditions at which experiments are performed. Under these circumstances, we hypothesize that accurate molecular models for water and carbon dioxide combined with computer simulation tools can offer an alternative but complementary way to estimate interfacial energies at coexistence conditions from a molecular perspective. We have evaluated the interfacial free energy of carbon dioxide hydrates at coexistence conditions (three-phase equilibrium or dissociation line) implementing advanced computational methodologies, including the novel Mold Integration methodology. Our calculations are based on the definition of the interfacial free energy, standard statistical thermodynamic techniques, and the use of the most reliable and used molecular models for water (TIP4P/Ice) and carbon dioxide (TraPPE) available in the literature. We find that simulations provide an interfacial energy value, at coexistence conditions, consistent with the experiments from its thermodynamic definition. Our calculations are reliable since are based on the use of two molecular models that accurately predict: (1) The ice-water interfacial free energy; and (2) the dissociation line of carbon dioxide hydrates. Computer simulation predictions provide alternative but reliable estimates of the carbon dioxide interfacial energy. Our pioneering work demonstrates that is possible to predict interfacial energies of hydrates from a truly computational molecular perspective and opens a new door to the determination of free energies of hydrates.

研究动机与目标

• 解释水合物–水界面自由能对 CO2 水合物在气候、存储与运输环境中的重要性。
• 提出通过分子模拟的方法在共存条件下估算界面自由能。
• 将 Mold Integration 方法应用于使用可靠水和 CO2 模型的 CO2 水合物–水体系。
• 评估基于模拟的界面能量是否与实验推断一致，并为水合物热力学提供计算路径。

提出的方法

• 采用 Mold Integration (MI) 来计算水合物–水界面自由能 gamma_hw，从形成水合物板条的可逆功中得到。
• 在 CO2–水混合物中使用锚定在水氧晶格位置的魅力位点模具来成核一个水合物板。
• 通过热力学积分在开启水–模具相互作用时计算可逆功 Delta G_hw（Delta G_hw = N_w*epsilon_m + Delta G_m）。
• 确定使可逆性达到最佳的模具半径 r_w^0，并将 gamma_hw 外推至 r_w^0。
• 在 coexistence 条件下的 NP_zAT 集合中进行模拟（P=40 MPa, T=287 K），使用 TIP4P/Ice 水和 TraPPE CO2，并对交互作用采用修正的 Lorentz–Berthelot 规则。
• 考虑到三相水合物–水–CO2 体系的有限尺寸及共存性影响。

Figure 1: Pressure-temperature projection of the three-phase coexistence line ( \ch CO2 hydrate – \ch H2O – \ch CO2) of the \ch CO2 hydrate. Red circles represent the experimental data taken from the literature Nakano1998a , the blue curve the experimental vapor pressure of pure \ch CO2 Lemmon2019a

实验结果

研究问题

• RQ1 在水合物–水共存条件下，CO2 水合物与水之间的界面自由能是多少？
• RQ2 Mold Integration 能否为采用真实水/CO2 模型的 CO2 水合物提供可靠的第一性原理界面能估计？
• RQ3 选择的模具半径 r_w 如何影响计算得到的界面能量，如何将 gamma_hw 外推至最佳 r_w^0？
• RQ4 模拟得到的界面能是否与现有的间接实验测量与热力学定义一致？

主要发现

• 基于 MI 的计算给出在共存时的 CO2 水合物–水界面自由能。
• 使用经过调参的 TIP4P/Ice 水和 TraPPE CO2 及交叉相互作用，可再现已知的解离/三相线，从而实现共存模拟。
• 该方法需要确定一个最佳模具半径 r_w^0，并将 gamma_hw 外推至该值。
• 使用固定模具深度 epsilon_m = 8 k_B T 以确保模具位点的完全占据，同时在路径上避免跨越一级相变。

Figure 2: Snapshots showing the crystallization of the \ch CO2 hydrate phase from the water- \ch CO2 two-phase coexistence at $40\,\text{MPa}$ and $287\,\text{K}$ . Different parts show the simulation box for (top) $t=0$ , (middle) $t=50\,\text{ns}$ , and (bottom) $t=200\,\text{ns}$ . Red and white