基于B3LYP/LanL2DZ和B3LYP/aug‐cc‐pVTZ方法计算铜络合物的铜同位素分馏

详细信息   

• 责任者：杨军丽^c中国科学院计算地球动力学重点实验室^d1987^e硕士研究生 ; 李永兵 ; 田会全 ; 刘善琪 ; 孙玥^c西安科技大学 ; 任东方^c中央民族大学 ; 刘建明^c中国科学院地质与地球物理研究所 ; 石耀霖
• 刊名：地学前缘
• 出版年：2014
• 卷期：21^b5总109
• 页码：116～127
• 索书号：P206.6/208-8
• 图表：^c表10参44插图
• 文摘期号：2014/11

摘要

铜的氯络合物和硫氢络合物是铜在成矿流体中稳定存在并参与地球化学迁移的重要形式;文中利用量子化学方法B3LYP/LanL2DZ和B3LYP/aug‐cc‐pVTZ研究计算了Cu⁺络合物(CuCl、[CuCl₂]^-、[CuCl₃]^2-、CuCl(H₂O)(真空)、CuHS，0～600 ℃)和Cu²⁺络合物(CuCl⁺、CuCl₂(真空)、CuOH⁺、CuHS⁺、Cu(HS)₂，0～300 ℃)的简约配分函数比。然后对比前人的理论与计算成果，讨论了不同基组方法计算铜络合物的铜同位素分馏的适应性。计算结果表明：利用LanL2DZ得到的Cu⁺氯络合物的分馏结果103lnβ65‐63比应用其他方法计算的结果偏大；而基于LanL2DZ基组计算Cu²⁺络合物的结果比aug‐cc‐pVTZ基组偏小。利用赝势基组LanL2DZ计算铜同位素分馏和实验结果偏差比较大，可能的原因是基于赝势基组LanL2DZ对上述络合物开展结构优化时，键长值比实验值偏大所致。因此，从理论计算上看，利用6‐311+G(d，p)基组可能更适合铜络合物的铜同位素分馏计算。虽然这些不同基组计算的结果存在差异，但与前人的实验结果相比，各种理论计算结果都在误差允许范围之内。鉴于此，在利用第一性原理计算同位素分馏系数时，如果计算条件允许，最好利用多种基组计算并作对比分析。