Jin Dongliang, Zhong Jing
School of Petrochemical Engineering, Changzhou University, 213164 Changzhou, P.R. China.
J Phys Chem B. 2023 Jun 15;127(23):5295-5307. doi: 10.1021/acs.jpcb.3c01768. Epub 2023 Jun 5.
The Gibbs-Thomson (GT) equation finds that the shift in the freezing/melting temperature under confinement with respect to its bulk counterpart is inversely proportional to the pore size. This century old relation successfully elaborates the freezing experiments of many fluids (e.g., water, molten salt), while it fails in quantitatively predicting the phase stability of the nonstoichoimetric crystals (e.g., gas hydrates). Based only on the crystal/liquid coexistence, we here revisit the GT equation to treat the multicomponent compounds within a slit confined geometry. In addition to the interfacial energy contribution, the extended GT equation accounts for the excess free energies associated with the composition variations upon the freezing/melting transition. Using the direct coexisting method (DCM), we first probe the melting temperatures of a face-centered cubic (fcc) crystal confined in slit pores as well as its bulk counterpart. The melting temperature under confinement is shown to be depressed compared to the bulk. We then turn to estimate the parameters entering the GT equation using several independent molecular simulations. The melting temperature depression observed in the DCM simulations is found to be well described by the GT equation if used with accurate estimates of the pore/crystal and pore/liquid interfacial tensions. Finally, using the above molecular modeling strategies, we show that the GT equation with the composition correction successfully predicts the shifted melting temperature of methane hydrate confined in porous solids. For such nonstoichoimetric compounds under confinement, accounting for the composition effects is of utmost importance as it exhibits a non-negligible contribution to the GT description. The extended GT equation can be expected to investigate the capillary freezing of the nonstoichoimetric compound in nanopores and to provide a better understanding of the pore body.
吉布斯-汤姆逊(GT)方程表明,受限条件下的凝固/熔化温度相对于其体相的变化与孔径成反比。这个有百年历史的关系成功地阐述了许多流体(如水、熔盐)的凝固实验,但在定量预测非化学计量晶体(如气体水合物)的相稳定性方面却失败了。仅基于晶体/液体共存,我们在此重新审视GT方程,以处理狭缝受限几何结构中的多组分化合物。除了界面能贡献外,扩展的GT方程还考虑了凝固/熔化转变时与成分变化相关的过量自由能。使用直接共存方法(DCM),我们首先探究了狭缝孔中受限的面心立方(fcc)晶体及其体相的熔化温度。结果表明,受限条件下的熔化温度低于体相。然后,我们使用几个独立的分子模拟来估计进入GT方程的参数。如果使用准确估计的孔/晶体和孔/液体界面张力,发现DCM模拟中观察到的熔化温度降低可以用GT方程很好地描述。最后,使用上述分子建模策略,我们表明具有成分校正的GT方程成功地预测了受限在多孔固体中的甲烷水合物的移动熔化温度。对于这种受限条件下的非化学计量化合物,考虑成分效应至关重要,因为它对GT描述有不可忽略的贡献。扩展的GT方程有望用于研究纳米孔中非化学计量化合物的毛细管凝固,并更好地理解孔体。