Lu Shih-I
Department of Applied Chemistry, Fooyin University, 151 Chinhsueh Road, Ta-Liao Hsiang, Kaohsiung Hsien, 831 Taiwan.
J Chem Phys. 2005 May 15;122(19):194323. doi: 10.1063/1.1899125.
Ab initio calculations of transition state structure and reaction enthalpy of the F + H2-->HF + H reaction has been carried out by the fixed-node diffusion quantum Monte Carlo method in this study. The Monte Carlo sampling is based on the Ornstein-Uhlenbeck random walks guided by a trial wave function constructed from the floating spherical Gaussian orbitals and spherical Gaussian geminals. The Monte Carlo calculated barrier height of 1.09(16) kcal/mol is consistent with the experimental values, 0.86(10)/1.18(10) kcal/mol, and the calculated value from the multireference-type coupled-cluster (MRCC) calculation with the aug-cc-pVQZ(F)/cc-pVQZ(H) basis set, 1.11 kcal/mol. The Monte Carlo-based calculation also gives a similar value of the reaction enthalpy, -32.00(4) kcal/mol, compared with the experimental value, -32.06(17) kcal/mol, and the calculated value from a MRCC/aug-cc-pVQZ(F)/cc-pVQZ(H) calculation, -31.94 kcal/mol. This study clearly indicates a further application of the random-walk-based approach in the field of quantum chemical calculation.
本研究采用固定节点扩散量子蒙特卡罗方法对F + H2→HF + H反应的过渡态结构和反应焓进行了从头算计算。蒙特卡罗采样基于由浮动球形高斯轨道和球形高斯双原子构建的试探波函数引导的奥恩斯坦-乌伦贝克随机游走。蒙特卡罗计算得到的势垒高度为1.09(16) kcal/mol,与实验值0.86(10)/1.18(10) kcal/mol以及使用aug-cc-pVQZ(F)/cc-pVQZ(H)基组的多参考型耦合簇(MRCC)计算得到的值1.11 kcal/mol一致。基于蒙特卡罗的计算还给出了与实验值-32.06(17) kcal/mol和MRCC/aug-cc-pVQZ(F)/cc-pVQZ(H)计算得到的值-31.94 kcal/mol相似的反应焓值,即-32.00(4) kcal/mol。本研究清楚地表明了基于随机游走的方法在量子化学计算领域的进一步应用。