扩散量子蒙特卡洛方法计算氢提取反应的能垒
详细信息    查看官网全文
摘要
量子蒙特卡洛方法具有很高的并行效率,计算量随体系大小的标度较低,而且具有很高的计算精度,能用于较大体系的高精度计算。本文中,我们采用固定节面的扩散量子蒙特卡洛方法(FN-DMC)计算燃烧化学中的一类重要反应:氢提取反应的能垒高度。我们用FN-DMC计算了Trular等发展交换相关泛函时所用HTBH38/04的19个氢提取反应的正逆向势垒高度。在FN-DMC计算中采用Jastrow因子与单行列式波函数乘积形式的引导波函数。其中的单行列式波函数中的轨道分别采用全电子基组6-311+G(2d,p)以及两个不同的赝势及相应基函数,在LDA近似所得到。我们用FN-DMC计算的结果与B3LYP、M06-2X、CCSD、CCSD(T)和显相关CCSD(T)计算结果以及参考结果做了对比。计算结果显示CCSD(T)方法和显相关CCSD(T)方法得到的结果与参考结果吻合最好,而FN-DMC对这些反应的能垒高度也给出了高精度的结果:与最优势垒估计值的平均偏差为0.98kcal/mol,在化学精度以内。但是CCSD(T)由于计算量非常大,难以用于较大的分子体系,而FN-DMC可用于较大分子的氢提取反应能垒的高精度计算。
Quantum monte-carlo method can be applied to calculations of relatively large systems with high accuracy due to its high parallelization efficiency and low computational scaling.In this work, the fixed-node diffuse quantum monte-carlo(FN-DMC) is adopted to calculate barrier heights of hydrogen abstraction reactions.Forward and backward barrier heights of 19 H-abstraction reactions employed in Trular's group for developing XC functional are calculated with FN-DMC.Single Slater determinant with Jastrow factor is used as guide functions in calculations with the Slater determinant obtained from LDA calculations using all-electron 6-311+G(2d,p)or pseudopotentials and their corresponding basis functions.The calculated results are compared with those from B3 LYP, M06-2X, CCSD(T) as well as reference values.CCSD(T) results are in best agreement with reference values, while FN-DMC also affords barrier heights with an average error of 0.98kcal/mol.This result indicates that FN-DMC can be employed to calculate barrier heights of H-abstraction reactions of large molecules with high accuracy.
引文
[1]Manolo C.per;Kelly A.Walker;Salvy P.Russo.J.Chem.Theory Comput.2012,8,2255-2259
    [2]Zhao Y,Lynch BJ,Truhlar DG.Phys.Chem.Chem.Phys.2005,7,43-52.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700