地质流体中络合物同位素平衡分馏系数的第一性原理计算  

First principles calculations of equilibrium isotope fractionations for complexes in geological fluid

作　　者：孙明光 王佳新[1] 袁顺达[1] SUN MingGuang;WANG JiaXin;YUAN ShunDa(MLR Key Laboratory of Metallogeny and Mineral Assessment,Institute of Mineral Resources,Chinese Academy of Geological Sciences,Beijing 100037,China;MOE Key Laboratory of Orogenic Belts and Crustal Evolution,Peking University,Beijing 100871,China)

机构地区：[1]中国地质科学院矿产资源研究所自然资源部成矿作用与资源评价重点实验室,北京100037 [2]北京大学造山带与地壳演化教育部重点实验室,北京100871

出　　处：《矿床地质》2022年第2期372-382,共11页Mineral Deposits

基　　金：国家自然科学基金重大研究计划重点支持项目(编号:92062218)资助。

摘　　要：地质流体中不同络合物之间的稳定同位素平衡分馏系数是利用稳定同位素示踪地质流体演化过程的必要前提。但由于受实验条件的限制,流体中部分元素不同络合物之间的同位素分馏系数难以通过实验精确测定。基于第一性原理的理论计算是获得溶液中不同络合物之间稳定同位素平衡分馏系数的一种重要手段,能够获取各种极端条件下同位素分馏系数。在应用第一性原理计算溶液中不同络合物间同位素平衡分馏系数时,需要对溶液的动态特性和水分子之间的相互作用进行模拟。目前,显式溶剂模型和隐式溶剂模型是模拟地质流体中水溶液环境的2种主要模型,均已成功应用于同位素分馏系数的计算。文章系统介绍了这2种模型的计算方法,并且对比研究了这2种模型计算约化配分函数比(β)的优势及不足,结果显示,2种模型的主要差异表现在3个方面:(1)计算耗时的不同,隐式溶剂模型计算的原子数目少,结构单一,因此耗时少;显式溶剂模型计算的原子数目多,结构复杂多样,因此在同一计算和理论水平条件下,耗时比隐式溶剂模型多;(2)水分子数和构型的不同,隐式溶剂模型只有一个构型,因而不能评估络合物构型对β值的影响;而显式溶剂模型能够通过添加不同数量的水分子及建立不同的初始构型有效评估水分子数量及构型对计算结果的影响;(3)隐式溶剂模型将溶质周围的水环境定义为具有均匀介电常数的连续介质,没有充分考虑水的动态特性和水分子之间的相互作用,导致计算结果误差较大;而显式溶剂模型能更好的处理水分子的动态特征及分子间形成的氢键,计算结果更为准确。通过2种模型的对比研究表明,显式溶剂模型能够更好的模拟水分子和溶质间的相互作用,计算获得的各络合物之间的同位素平衡分馏系数更为准确。本研究结果可以为今后开展相关同�The equilibrium isotope fractionation factors between complexes in nature fluids play a key role in using sta-ble isotope to trace the evolution process of nature fluids.However,due to the limitations of experiments,it is dif-ficult to accurately determine the isotope fractionation factors between complexes.Based on the first principle,quantum calculation is an important method to obtain the stable isotope equilibrium fractionation factor between different complexes,which cannot be obtained by experimental method for various extreme conditions.It is nec-essary to simulate the dynamic characteristics of solution and the interaction between water molecules when ap-plying the first principle to calculate the equilibrium isotope fractionation factors between complexes.At present,there are two methods to simulate the environment of aqueous solution in geological fluid,explicit solvent model and implicit solvent model,which have been successfully applied to the calculation of equilibrium isotope frac-tionation factors of many elements.In this study,we systematically introduce the two models,and compare the advantages and disadvantages of the two models in calculating reduced isotope partition function ratio(RPFR orβ)value.The results show that the differences between the two methods are mainly reflected in three aspects:①for computing time,due to the small number of atoms and a single structure of implicit solvent model,the cal-culation time is less;however,due to the large number of atoms and complex structure of the explicit solvent model,the calculation takes a long time under the same calculation and theoretical level with the implicit solvent model;②For different configurations,the implicit solvent model has only one configuration,which cannot evalu-ate the effect of configuration onβvalues,while the explicit solvent model can effectively evaluate the influence of configuration on the results by establishing different initial configurations;③The implicit solvent model de-fines the aqueous solution environm

关 键 词：地球化学 第一性原理 同位素分馏 络合物 含矿流体 显式溶剂模型和隐式溶剂模型 

分 类 号：P597[天文地球—地球化学]