界面拓扑选择将结构基元与甲烷水合物成核中的自由能势垒联系起来

《Journal of Colloid and Interface Science》：Interfacial topological selection links structural motifs to the free-energy barrier in methane hydrate nucleation

　　**摘要****假说**甲烷水合物的成核是一个随机过程，对于能源和流动保障技术至关重要，然而与其自由能势垒跨越相关的结构变化仍未得到充分解析。我们提出假说：界面氢键网络中的瞬态拓扑基元能够提供超越体相核生长的势垒跨越信息。**模拟**为检验该假说，我们将分子动力学（MD）模拟与机器学习（ML）相结合。我们进行了20次独立的长时标甲烷水合物成核MD模拟，并使用分层拓扑环算法识别瞬态结构基元，包括多面体笼和三元水环聚集（TWRA）。应用时滞独立成分分析（tICA）构建自由能景观，同时利用极端梯度提升（XGBoost）配合局部可解释模型无关解释（LIME）量化结构描述符与自由能坐标之间的关联。采用轨迹级数据划分进行模型评估。**发现**在成核轨迹中，势垒攀升伴随着拓扑不兼容的2-1-0型TWRA的耗竭以及水合物-水界面上兼容的0-3-0和0-2-1型基元的富集。XGBoost/LIME分析识别出这些TWRA比例是自由能坐标的信息性结构描述符。在轨迹级留一交叉验证中，将三个关键TWRA比例与瞬时自由能结合，相比于仅使用瞬时自由能，改善了在50-100 ns后自由能演化的预测。这些结果支持界面拓扑选择作为势垒跨越的结构特征，并将局部氢键网络重组与甲烷水合物成核的自由能演化联系起来。**引言**笼形水合物是在高压低温条件下形成的类冰晶体，其中水分子组装成笼状框架，包裹甲烷、二氧化碳、四氢呋喃或环戊烷等客体分子。天然气水合物（NGHs）主要由甲烷水合物组成，具有巨大的潜在能源资源和异常高的能量密度，在标准条件下每单位体积可储存多达160立方米的甲烷[1]。这种高储存容量使NGH储层成为极具前景但技术上具有挑战性的未来能源。相反，水合物在油气管道内自发生成（图1a）会导致严重堵塞，这是一种被称为"流动保障"的主要运营危害，造成巨大的经济损失[2]。除了天然存在外，笼形水合物还可被合成用于工业应用。其笼状结构如同分子筛，通过范德华（vdW）作用选择性地捕获小分子。这一独特性质使其可用于二氧化碳封存[3]、[4]、[5]、[6]，氢气和天然气储存[7]、[8]，气体混合物分离[9]以及海水淡化[4]、[10]。控制水合物形成的能力——无论是在管道中阻止其形成还是促进其储存/分离——关键依赖于对成核分子机制的理解，即该过程难以捉摸的初始步骤。成核是在过冷/过饱和系统中发生的随机过程，其中必须先形成临界核才能进行快速生长。尽管其具有根本重要性，但水合物成核的分子途径仍未被充分理解，并由几个竞争性假说所解释，包括不稳定团簇假说（LCH）[11]、[12]，局部结构假说（LSH）[13]，笼吸附假说[14]和团块假说[15]。Walsh [16]、[17]、[18]等人揭示了早期水合物核是由主要由5126?笼组成的非晶态聚集体，这与LSH提出的晶体结构相反。Jacobson及同事提出了团块假说，认为初始步骤涉及"团块"的形成，即由溶剂分隔的客体分子聚集体，随后周围水分子重组为多面体笼和非晶核。从非晶态到晶态水合物核的多步演化已在多次分子动力学（MD）模拟中被观察到[19]、[20]、[21]。Li等人[22]进一步识别出三元水环聚集（TWRA）是该演化中的关键结构。然而，特定分子水平结构基元、支配成核的自由能景观（FEL）以及随机成核时间之间的定量联系仍然是一个重大知识空白。弥合这一空白对于预测性控制至关重要。在本工作中，我们将MD模拟与机器学习（ML）相结合，以解析局部拓扑如何沿甲烷水合物成核途径演化。我们进行了20次独立的长时标MD模拟，并使用我们的分层拓扑环（HTR）算法[23]识别和量化瞬态结构元素，包括笼类型和拓扑水环聚集（TWRA，图1b）。对局部水序描述符进行时滞独立成分分析（tICA）以构建成核的自由能景观（图1c），随后将笼和TWRA描述符与所得自由能坐标进行分析（图1d）。我们重点关注当系统接近自由能势垒时，含有两个四元环的TWRA的显著耗竭，以及0-3-0和0-2-1型基元的富集。随后使用可解释的XGBoost分析量化这些结构描述符对自由能预测的贡献（图1e），并通过轨迹级验证测试它们是否改善了后续自由能演化的预测（图1f）。我们不是将任何单一TWRA视为成核的直接证实原因，而是评估拓扑选择作为势垒跨越过程中伴随界面重组的集体结构过程。**章节片段****用于甲烷水合物成核的分子动力学模拟**所有分子模拟均使用GROMACS [24]进行。初始模拟盒子尺寸为50 × 50 × 50 ?3，包含3487个H?O分子和550个CH?分子。该分子组成与Li等人[22]用于研究甲烷水合物成核过程中TWRA演变的组成相匹配；Xu等人[25]的纯水参考系统也使用了相同的组成和初始盒子尺寸。TIP4P/Ice [26]和OPLS-AA [27]力场被用于...**成核途径的整合结构与动力学映射**为建立水合物成核的基础结构描述符，首先分析系统构型的时序演化（图2）。笼状结构在200至800 ns之间逐渐形成，表明非晶核的出现（图2a）。当甲烷分子按其四体序参量对F_IS着色时，F_IS升高的区域始终与笼的形成重合，确认F_IS作为成核的可靠结构指标（图2b）。**结论**本工作将MD模拟与ML相结合，解析了甲烷水合物成核过程中水合物-水界面上三阶段拓扑选择模式。第一阶段包含TWRA之间的快速互变，产生局部基元的波动池。在第二阶段，拓扑不兼容的2-1-0型基元的耗竭伴随向自由能势垒的攀升，而0-3-0和0-2-1型基元变得富集。经LIME解释的XGBoost识别出这些TWRA比例...**作者贡献声明****郝永超：** 撰写–原始稿件，可视化，软件，方法论，调查研究，形式化分析，数据整理。**曲永晓：** 可视化，软件，方法论，形式化分析，数据整理。**罗高阳：** 可视化，软件，方法论，形式化分析，数据整理。**张志森：** 撰写–审阅与编辑，监督，项目管理，资金获取，概念化。**吴建阳：** 撰写–审阅与编辑，监督，资源，项目管理...**利益冲突声明**作者声明他们没有已知的竞争性经济利益或个人关系，这些关系可能会被认为影响本论文报告的工作。**致谢**本工作得到**国家自然科学基金委员会**（资助编号：12572127和12172314）、**中国-甘肃省发展项目：产业多元化组成部分**（资助编号：25YFGA076）、**福建省自然科学基金**（资助编号：2025 J01006）以及**"111"项目**（B16029）的资金支持。感谢厦门大学信息网络中心的方少荣和吴天府提供的超级计算机GPU时间...郝永超 | 曲永晓 | 罗高阳 | 张志森 | 吴建阳