高压高温地质流体中氢氟酸酸性常数(pKa)的机器学习分子动力学研究及其对成矿作用的启示

本文采用自动化机器学习分子动力学（MLMD）方法，结合第一性原理（ab initio）数据，系统计算了氢氟酸（HF(aq)）在高温高压（150-600 °C, 500 bar至20 kbar）条件下的酸性常数（pKa）。研究结果表明，HF(aq)在深部地壳和上地幔环境（尤其是超临界条件）下是比以往估计更强的酸，其解离程度增强，氟离子（F–）的有效浓度提高，这显著增强了稀土元素（REE）、铀（U）等金属的溶解度和迁移能力。所得pKa值被成功拟合到基于密度依赖的Ryzhenko-Bryzgalin热力学模型（MRB）中，为深入理解氟在热液成矿流体中控制元素迁移性和pH的关键作用提供了新的热力学数据。

引言

氟（F）是许多热液成矿流体的关键组分，常见于萤石、磷灰石、云母等矿物中。其在热液中的行为，特别是氟离子（F^–）的可用性，受氢氟酸（HF(aq)）解离平衡的控制。然而，由于实验技术的限制，获取在地球深部（地壳和上地幔）典型的高压高温（P-T）条件下HF(aq)可靠的热力学性质一直充满挑战。本研究旨在利用一种基于第一性原理数据训练的自动化机器学习分子动力学（MLMD）方法，来计算HF(aq)在宽泛温压范围内的pK_a值，并据此建立预测模型，以评估HF(aq)在深部地质流体中对元素迁移性和pH的控制作用。

方法

研究首先进行了基于第一性原理的分子动力学（AIMD）模拟，使用CP2K/QUICKSTEP软件包，采用BLYP交换关联泛函、GTH赝势和DZVP基组，模拟了包含55个水分子和1个HF分子的体系在不同温压条件（100-600 °C, 500 bar至20 kbar）下的行为，获取了径向分布函数（RDF）和配位数（CN）等结构信息。

酸性常数（pK_a）的计算采用了垂直能隙法，通过自由能微扰（FEP）和热力学积分（TI）来计算反应的标准吉布斯自由能变。为了避免质子转移反应带来的非物理构象，在计算中引入了弱谐束缚势。作为参照，同时对水合氢离子（H₃O⁺）进行了相同的去质子化自由能计算。pK_a值最终通过比较HF(aq)和H₃O⁺的去质子化自由能差获得。

为了提升采样效率和精度，研究构建了机器学习势函数（MLP）。利用从AIMD模拟中获得的初始构型及其对应的DFT能量和原子力数据，使用DeepMD-kit包训练了多个深度神经网络势函数。通过模型偏差方法，迭代地筛选出MLMD模拟轨迹中力偏差较大的构型，补充到训练集中重新训练MLP，直至收敛。最终训练得到的MLP用于在LAMMPS软件中进行长时间尺度（0.5纳秒）的NVT系综模拟，以计算各温压条件下的pK_a值。

结果与讨论

AIMD模拟结果显示，在低温（100-200 °C）条件下，F-H键长略长，配位数低于1.0，表明HF存在部分解离。而在高温（>250 °C）条件下，F-H键长缩短至约1.00 Å或更短，配位数接近或达到1.0，表明HF(aq)的缔合程度更高。高压有助于维持HF的缔合状态。

MLMD模拟的验证表明，其预测的能量和原子力与DFT计算结果高度吻合，均方根误差（RMSE）极低，证明了MLP的准确性。研究特别指出，由于采用了弱束缚势，其对溶剂化质子结构动态的影响可忽略不计，因此无需像一些早期AIMD研究那样引入质子平移熵校正项，这保证了pK_a计算在高热条件下的自洽性。

计算得到的pK_a值显示其与水的密度（ρ_H2O）存在明显的依赖关系：pK_a随密度降低（高温或低压）而增大，意味着HF(aq)在低密度超临界流体中酸性减弱（更倾向于以分子形态存在）；反之，在高压高密度条件下，pK_a值较低，HF(aq)解离程度更高。与基于Helgeson-Kirkham-Flowers状态方程（HKF EOS）外推的文献值相比，本研究MLMD计算出的pK_a值在中等至低流体密度条件下要低1.5-2个对数单位，表明HF(aq)在深部高温条件下是比以往认知更强的酸。

为了便于地质应用，研究将MLMD计算出的pK_a值拟合到了修正的Ryzhenko-Bryzgalin（MRB）热力学模型中。拟合得到的模型参数允许对HF(aq)的解离常数在更广泛的P-T范围内进行外推。图示结果表明，在低压下，pK_a随温度升高先增后降，在超临界区附近出现急剧变化，这与流体密度的剧烈变化相关；而在高压下，pK_a随温度的变化相对平缓，且整体数值较低，反映了高压促进酸解离的效应。

地质意义

本研究计算出的HF(aq) pK_a值低于传统外推值，意味着在深部热液条件下，流体中游离F^–的浓度可能显著高于之前的估计。这直接影响到氟与金属离子（如Fe、REE、U等硬阳离子）形成络合物的能力。更高的F^–活度以及由此导致的流体酸性增强，可以显著提高诸如稀土元素和铀等关键金属在热液中的溶解度和迁移效率。这一发现有助于解释在斑岩型铜钼金矿床、铁氧化物铜金（IOCG）矿床等大型成矿系统中普遍存在的富氟现象（如萤石的广泛出现）。新的热力学数据为更精确地模拟氟在控制深部地质流体中元素迁移、沉淀以及流体-岩石相互作用过程中的角色提供了关键参数，对理解这些重要矿床的形成机制具有深远意义。