精确计算氢键键能的Gaussian理论.pdf

第 20 卷 第 1 期2003 年?2 月贵州大学学报(自然科学版)Journal of Guizhou University(Natural Sciences)Vol.20 No.1Feb.2003文章编号?1000-5269(2003)01-0036-06精确计算氢键键能的 Gaussian 理论*?梁?雪,王一波(贵州大学化学系,贵州 贵阳?550025)摘?要?将研究单分子能量非常成功的 Gaussian 理论用于氢键键能的精确计算,选择计算精度高、所需计算资源较少的 Gaussian-3(G3)理论研究了一系列典型氢键体系,与 CCSD(T)结果比较,误差均在?2kJ/mol 范围内。同时还能大量地节省计算机时和磁盘空间。关键词?Gaussian 理论;氢键;键能中图分类号?O641?13?文献标识码?A0?引?言氢键是一类既复杂又重要的相互作用,尽管这种作用比共价键弱得多,却对描述分子许多行为以及决定物质的物理性质和化学性质非常重要。尤其在生物分子体系中,许多生理过程都是通过氢键作用体现出来的。Pauling L.在经典著作?化学键的本质?中曾指出:?当结构化学方法进一步被应用到生理问题上时,人们将会发现氢键在生理学上的意义比任何一个结构特点都大?。可见氢键因其与生命过程密切相关,早就受到人们的高度重视。由于实验研究的困难性,从理论上精确地计算氢键键能引起了人们的广泛兴趣。它需要采用高精度的量子化学计算方法,有效地计算电子相关能,提高基函数效率和消除基函数重叠误差 BSSE,比单分子问题的计算困难得多。近年来,Pople J.A.及其合作者提出了一系列用于精确计算分子能量的 Gaussian-1,Gaussian-2 和Gaussian-3理论 1 6,该方法无论在结果的精度上还是计算资源消耗上,与其它许多高精度的传统量子化学计算方法相比都有明显的优势。它主要用于计算分子在其平衡构型下的能量性质,包括离解能、电离能、质子亲和能和电子亲合能等,其结果与精确实验值比较的误差在?2kcal/mol 范围内,对于尚无实验结果的分子体系能给出相当准确的预测。此外,还可用于热化学 7 12和化学反应途径 13 17的研究。Gaussian 理论的创立者把研究的目标放在了单分子上,在其几百个测试算例中,没有一个涉及到氢键体系。有关 Gaussian 理论应用的文献很多,但研究的对象几乎仅局限于单分子,即使有见其用于氢键体系的研究报道 18,也是着重于对体系本身的探讨上,而没有对方法进行研究。Gaussian 理论用于氢键领域的研究几乎还是一片空白。Gaussian 理论采用了高级相关方法用于计算分子的能量,基函数中包含了高角量子数的极化函数和弥散函数,这些与研究氢键所用的方法同源,也是能够把它用于氢键键能研究的首要和必要条件。本文就Gaussian 理论对于氢键键能计算进行了系统的研究。取得令人满意的结果。1?方?法Gaussian 理论基于 ab initio 分子轨道理论,它经历了 Gaussian-1(简称 G1)、Gaussian-2(简称 G2)和Gaussian-3(简称 G3)几个主要阶段,此外,在这些方法的基础上稍作改进用来处理具体问题的研究还有很多。19 32我们选取了计算精度较高,计算量较小的 G3理论,将其用于氢键键能的计算。Gaussian-3 理论的基本思想如下:(1)在 HF/6-31G(d)水平下获得初步的平衡几何结构。该水平下获得的几何构型将用于振动频率的计算获得零点振动能。?收稿日期:2002-12-30项目来源:国家自然科学基金资助项目(2976001,29992590-1)。作者简介:梁?雪(1977-),女,硕士研究生。研究方向:量子化学计算;王一波(1962-),男,教授。研究方向:量子化学从头计算。(2)用 MP2(full)/6-31G(d)确定几何结构,它将作为最终的平衡结构用于后面一系列单点能量的计算。(3)进行更高级别的一系列单点能量的计算。首先用四级微扰获得基本能量,即完成 MP4/6-31G(d)的计算,该结果还必须加入一系列修正。(a)获得弥散函数的修正?E(+)?E(+)=E MP4/6-31+G(d)-E MP4/6-31G(d)(b)对非氢原子加入高极化函数和对氢原子加入 p 极化函数的修正?E(2df,p).?E(2df,p)=E MP4/6-31G(2df,p)-E MP4/6-31G(d)(c)二次组态相互作用获得相关能的修正?E(QCI)?E(QCI)=E QCISD(T)/6-31G(d)-E MP4/6-31G(d)(d)对于大基函数和由于弥散函数和极化函数各自基函数扩大产生的非加和性的修正?E(G3Large)?E(G3Large)=E MP2(full)/G3Large-E MP2/6-31G(2df,p)-E MP2/6-31+G(d)+E MP2/6-31G(d)(4)MP4/6-31 G(d)获得的基本能量加上上述四个修正后,再加上自旋-轨道修正?E(SO)(仅对原子),就可以得到一个组合的能量 E组合E组合=E MP4/6-31G(d)+?E(+)+?E(2df,p)+?E(QCI)+?E(G3Large)+?E(SO)(5)高级相关修正(HLC)用来说明能量计算中的残余误差?Ee(G3)=E组合+E(HLC)这里,HLC 包括两个部分,-An?-B(n?-n?)用于分子,-Cn?-D(n?-n?)用于原子。n?和 n?分别是?和?的价电子数,并且 n?n?。在 G3 理论中,A=6.386 millihartrees,B=2.977 millihartrees,C=6.219millihartrees,D=1.185 millihartrees。(6)最后,加上第一步频率计算时获得的零点能,就获得了 0K 时 G3 的总能量?E0(G3)=Ee(G3)+E(ZPE)=E MP4/6-31G(d)+?E(+)+?E(2df,p)+?E(QCI)+?E(G3Large)+?E(SO)+E(HLC)+E(ZPE)它是 QCISD(T)(full)/G3Large 的一个很好的近似。2?结果与讨论2?1?参考值?EHF、?EB3LYP、?ECCSD(T)的计算由于从实验上获取氢键的结合能非常困难,所以大多数体系尚未有较好的实验数据参考。为此,我们对选取的体系分别采用从头算 Hartree-Fock、密度泛函 B3LYP 和簇偶合 CCSD(T)方法作了一个参考值,将 G3 的计算结果与其比较。几何构型均采用 G3 中获得的结构,同时还采用了 Boys 和 Bernardi 提出的均衡校正法 CP(Counterpoise Procedure)33进行BSSE 校正,获得的结合能?EHF、?EB3LYP和?ECCSD(T)值列于表 2。所有计算均采用 Gaussian98 程序 48,全部计算工作在贵州省高性能计算化学重点实验室完成。2?2?Gaussian-3(G3)理论用于研究氢键键能G3 获得的结合能以及各相关修正项的值列于表 1。表 1?Gaussian-3 计算的能量和修正(kJ?mol-1)体?系MP4/6-31G(d)?E(+)?E(2df,p)?E(QCI)?E(G3Large)?E(HLC)?E(ZPE)E0(H2O)2-29.642.261.630.343.880.008.14-13.38(HF)2-32.8517.830.420.740.540.006.93-6.39(HCl)2-7.32-0.830.180.100.290.002.76-5.41NH3?NH3(CS)-19.072.412.040.121.320.00-5.89-7.29NH3?NH3(C2h)-19.082.472.080.121.270.00-5.02-8.12H2O?HF-44.740.082.820.444.010.00-9.66-27.77H2O?HCl-30.593.02-0.490.684.830.007.01-15.55H2O?H2S(A)-11.83-2.011.690.15-0.100.00-4.82-7.28?37?第 1 期?梁?雪,王一波:精确计算氢键键能的 Gaussian 理论H2O?H2S(B)-15.910.99-1.020.434.160.00-4.74-6.61H2O?NH3-33.74-1.014.600.232.370.008.32-19.22H2O?N2-6.83-0.520.380.300.760.002.40-3.51NH3?HF-58.52-4.437.420.422.010.0011.27-41.83NH3?HCl-43.674.652.670.682.060.009.02-24.59H2S?HF-19.34-3.972.190.27-0.170.006.58-14.43H2S?HCl-13.04-0.830.780.23-0.530.004.85-8.55C2H2?HF-19.78-0.35-0.490.550.460.004.77-14.86C2H2?HCl-12.440.51-1.360.340.950.003.39-8.61C2H2?H2O-18.742.03-0.150.473.990.004.17-8.22HF?CO-17.351.530.791.22-2.900.004.50-12.22Cl-?H2O-68.985.503.770.09-2.61-1.755.30-58.69?注:以上各项的值均为 EAB-(EA+EB)由于参考值中没有考虑零点能的计算,所以我们将 E0的结果中去掉第二步频率计算时获得的零点能值(所得结果为 EG3),将其与参考值进行了比较,结果列于表 2。表 2?Gaussian-3 获得的?EG3与?EHF、?EB3LY P、?EC CSD(T)结果比较(kJ?mol-1)体?系?EG3?EHF?EB3LYP?ECCSD(T)Expt.(H2O)2-21.53-14.28-18.79-20.00-22.30?2.09(HF)2-13.32-8.35-12.53-18.23-18.65(HCl)2-8.17-2.18-4.91-7.17NH3?NH3(CS)-13.14-5.98-13.38-12.58NH3?NH3(C2h)-13.18-5.98-13.42-12.56H2O?HF-37.43-32.51-39.02-35.97-30.1?7.1H2O?HCl-22.56-14.20-22.87-21.25H2O?H2S(A)-12.10-5.81-8.96-11.51H2O?H2S(B)-11.35-3.52-7.35-10.49H2O?NH3-27.55-19.73-27.79-26.29H2O?N2-5.91-0.40-1.68-4.68NH3?HF-53.10-49.01-57.71-54.11NH3?HCl-33.62-24.25-41.85-35.62H2S?HF-21.01-14.66-21.09-20.59H2S?HCl-13.39-5.62-12.18-12.81C2H2?H2O-12.39-7.66-9.10-11.60C2H2?HF-19.62-13.72-17.56-17.97C2H2?HCl-12.00-5.06-8.32-11.09HF?CO-16.83-7.54-15.55-14.22Cl-?H2O-62.24-53.85-61.31-62.89?从表 2 的计算结果显示,?EG3与?EHF相差很大,这主要是因为Hartree-Fock 方法方法没有考虑电子相关,而电子相关对弱相互作用来说不能忽略,有的体系的电子相关能在总相互作用能中所占的比例几乎达到了一半34,所以该方法在用于研究氢键键能时所得结果不可靠。密度泛函 B3LYP 方法考虑了电子相关,可以用于研究氢键,但是我们可以发现,对于相互作用较强的体系,B3LYP 的结果与 CCSD(T)结果相差不大,与?EG3也能较好地吻合,但对于较弱的相互作用体系,B3LYP 获得的结合能?EB3LYP与 CCSD(T)以及G3的结果有了一定的差距,可以说,密度泛函方法只适用于研究较强相互作用体系。而用CCSD(T)方法获得的结合能?ECCSD(T)与 G3 获得的?EG3结果较为接近。CCSD(T)结果与 G3 结果比较除(HF)2体系误差稍大(为?4.91kJ/mol),同时还有 HF?CO 体系的误差为 2.61kJ/mol 外,其余体系?EG3与 CCSD(T)结果的误差均在?2kJ/mol 范围内,这已大大超过了实验所能及的精度。同时,对于有较好实验数据的体系,G3 获得的计算值与实验值也能很好地吻合。对于 H2O?H2S 体系,有两种不同的结构,分别代表H2O 作为质子授体的结构 A 和 H2S 作为质子授体的结构 B。从表 2 的计算结果显示,结构 A 较结构 B 稳定。这与我们先前的工作获得的结论一致 35。我们计算了?EG3与?ECCSD(T)的总的平均绝对偏差,对所选取的体系,?EG3与?ECCSD(T)的平均绝对偏差为 1.30 kJ/mol(即 0.31kcal/mol)。2?3?计算机时和所需磁盘空间的比较这里,我们选取了几个体系,分别用 G3 和 CCSD(T)方法进行了计算机时和所需磁盘空间的比较,结?38?贵州大学学报(自然科学版)?第 20 卷果列于表 3。很明显,在计算结果精度相当的同时,CCSD(T)方法所需的计算机时和磁盘空间都要大得多。普遍来看,CCSD(T)所需的 CPU 时间是 G3 的 20 40 倍,所需的外存是 G3 的 10 50 倍。如 H2O?H2O 的CPU 时间 CCSD(T)约是 G3 的 45 倍,H2O?HF 体系 CCSD(T)所需的 CPU 时间约是 G3 的 36 倍,所需硬盘空间 CCSD(T)约是 G3 的 50 倍。表 3?G3 与 CCSD(T)方法 CPU时间(分)和磁盘空间(MB)的比较体系方法CPU timeDisk storageH2O?H2OG315100CCSD(T)6862015H2O?HFG316116CCSD(T)5785798N2?H2OG337325CCSD(T)7662770?注:a.在 CCSD(T)计算中,基函数采用的是 aug-cc-pVTZ 基组b.计算在 P?800,SuSE Linux 环境下完成。cpu时间为复合物与单分子的 cpu时间的总和(几何优化和零点能计算的时间不包括在内)。2?4?关于零点能的讨论零点能的获得是在第二步 HF/6-31G(d)频率的计算中,并乘上校正因子 0.8929。由于从头算法获得的谐振动频率在与 IR 实验的直接观测频率相比较时都存在一定的误差,为了使理论计算能很好地重现观测频率,一般都需要乘上一个校正因子 36。该校正因子的取值一般由相应方法所产生的误差大小决定。由于氢键的研究必须采用高精度的量子化学方法,我们认为采用 HF/6-31G(d)方法,对于氢键体系零点能的精确计算是远远不够的。从前人发表的研究氢键的工作来看,考虑零点能计算的并不多见。由于零点能的计算比较困难,并受计算条件的限制,目前想要进一步提高零点能的计算精度并非易事。想要让结果更加精确,尤其是能更好地处理一些热化学问题,对零点能修正方法的改进是一项重要而迫切的工作。2?5?对特殊体系的讨论从表 2 中我们发现,G3 对该体系的处理不太不理想。对于(HF)2已有大量的实验 37 40和理论研究41 47。和二聚水一样,它也是强氢键的一个典型例子。为此,我们将该体系作为特例对 G3 进行了改进。由于 F 原子具有很强的电负性,使得它对 H 原子的吸引增强,为此我们在几何构型优化时加入了一个sp型弥散函数,并对氢原子加入了一个 p 极化函数,即初步几何构型优化在 HF/6-31+G(d,p)计算下完成,平衡结构的确定由原来的 MP2(full)/6-31G(d)变为 MP2(full)/6-31+G(d,p),将该方法用于计算(HF)2的结合能?E修。为了验证该方法对其他体系是否有效,我们同时选取了具有代表性的(H2O)2进行研究,其结果列于表 4。表 4?HF?HF体系结合能体系?Ee?E修?ECCSD(T)Expt.(HF)2?13.32?20.53?18.23-18.65(H2O)2-21.53(26)-22.33(31)-20.00-22.30?2.09?注:括符内的数值为复合物时间与单分子的 cpu时间总和,其单位为分。结合能的单位为 kJ?mol-1。从表 4 显示的结果可见,适当考虑弥散和极化对 HF?HF 体系的影响是非常有效的。对于(H2O)2结果也较好,但所需的计算资源有所增加。我们认为如果不是特殊体系,原始的 G3 方法获得的结合能的精度已经足够了,但是针对HF?HF 体系对 G3进行的修改对一些特殊体系的研究提供了一种思路。3?结论本文将研究单分子能量非常成功的 Gaussian 理论用于处理氢键键能的计算,其结果也同样成功。用Gaussian 理论计算氢键键能,不需考虑 BSSE 校正,其结果的精度与 CCSD(T)方法相当,对所选取的 20 个对象,其计算值与 CCSD(T)结果比较,平均绝对偏差分别为 1.30kJ/mol(即 0.31kcal/mol)。这不仅达到了化学精度,而且与原有 Gaussian 理论对单分子处理结果精度比较,也要好得多。原有 Gaussian 理论对单分子的处理,即使是比较好的 G3 的理论值与实验值比较,平价绝对偏差也有 1.01 kcal/mol。同时,该?39?第 1 期?梁?雪,王一波:精确计算氢键键能的 Gaussian 理论方法能节约大量的计算机时和磁盘空间。与实验值比较,G3的计算结果也能很好地重现实验数据。可以说,Gaussian 理论用于研究氢键键能,不失为一个好方法,其结果的精度足以为实验工作提供有价值的参考,我们也期待着更精确的实验和理论计算的验证。参考文献 1 Pople J A,Head-Gordon M,Fox D J,Raghavachari K and A.Curtiss L.Gaussian-1 theory:A general procedure for prediction of molecu?lar energies J.J Chem Phys,1989,0:5622.2 Curtiss L A,Raghavachari K,T rucks G W and Pople J A.Gaussian-2 theory formolecular energies of first-and second-row compounds J.J Chem Phys,1991,94:7221.3 Curtiss L A,Raghavachari,K and Pople J A.Gaussian-2 theory using reduced M?ller-Plesset orders J.J Chem Phys,1993,98:1293.4 Curtiss L A,Raghavachari K,Redfern P C,Rassolov V,and Pople J A.Gaussian-3(G3)theory for molecular containing first and second-row compounds J.J Chem Phys,1998,109:7764.5 Curtiss L A,Redfern P C,K.Raghavachari and Pople J A.Gaussian-3 theory using reduced M?ller-Plesset order J.J ChemPhys,1999,110:4703.6 Curtiss L A,Redfern P C,Raghavachari K and Pople J A.Gaussian-3 theory:a variation based on third-order perturbation theory and anassessment of the contribution of core-related correlation J.Chem Phys Lett,1999,313:600.7 Curtiss L A,Raghavachari K,Redfern P C and Pople J A.Assessment of Gaussian-2 and density functional theories for the computation ofenthalpies of formation J.J Chem Phys,1997,106:1063.8 Curtiss L A,Redfern P C,Raghavachari K and.Pople J A.Assessment of Gaussian-2 and density functional theories for the computation ofionization potentials and electron affinities J.J Chem Phys,1998,109:42.9 Curtiss L A,Raghavachari K,Redfern P C and Pople J A.Assessment of Gaussian-3 and density functional theories for a larger experimen?tal test set J.J Chem Phys,2000,112:7374.10 Raghavachari K,Stefanov B BandCurtiss L A.Accurate thermochemistry for larger molecules:Gaussian-2 theory with bond separation en?ergies J.J Chem Phys,1997,106:6764.11 Notario R,Castano O,Gomperts R,M.Frutos L and Palmeiro R.Organic thermo-chemistry at high ab initio Levels.3.A G3 study ofcyclic saturated and unsaturated hydrocarbons J.J Org Chem,2000,65:4298.12 Curtiss L A,Raghavachari K,C.Redfern P,Kedziora G S andPople J A.On comparison of experimental thermochemical data with G3 theo?ry J.J Phys Chem A,2001,105:227.13 Hrusak J,Schwarz H.A Gaussian-1 study of the CHOF+potential energy surface.T heoretical evidence for the existence of three experi?mentally distinguishable isomers in the gas phase J.J Phys Chem,1993,97:4659.14 Ma N L,Smith B J,Pople J A and Radorn L.Rearrangement and dissociative reactions of the methanol radical cation(CH3OHo+):A com?parison of theory and ex perimentJ.J A m Chem Soc,1991,113:7903.15 Wong M W,Radom L.Thermochemistry and ion-molecule reactions of isomeric C3H2o+cations J.J A m Chem Soc,1993,115:1507.16 Luna A,Yanez M.G2 ab initio calculations on the thermochemistry of P,N,Hn(n=0,2)and P,N,Hn+n=0,3)species and on thepotential energy surfaces of P,N,H3+singlet and triplet state cations J.J Phys Chem,1993,97:10659.17 Mebel A M,K.Morokuma and Lin M C.Ab initio molecular orbital study of potential energy surface for the NH+NO2reaction J.J ChemPhys,1994,101:3916.18 David W,Ball.G2 calculation of weak molecular complexes:H2O-NO J.J Phys Chem A,1997,101:4835.19 Baboul A G,Curtiss L A,Redfern P C and Raghavachari K.Gaussian-3 theory using density functional geometries and zero-point ener?gies J.J Chem Phys,1999,110:7650.20 Curtiss L A,Raghavachari K,Redfern P C and Pople J A.Gaussian-3 theory using scaled energies J.J Chem Phys,1997,112:1125.21 Curtiss L A,Redfern P C,Raghavachari K and Pople J A.Gaussian-3X(G3X)theory:Use of improved geometries,zero-point ener?gies,and Hartree-Fock basis sets J.J Chem Phys,2001,114:108.22 Curtiss L A,Carpenter J E,Raghavachari K and Pople J A.Validity of additivity approximations used in Gaussian-2 theory J.J ChemPhys,1992,96:9030.23 Curtiss L A,Raghavachari K and Pople J A.Gaussian-2 theory:Use of higher level correlation methods,quadratic configuration interac?tion geometries,and second-order M?ller-Plesset zero-point energies J.J Chem Phys,1995,103:4192.24 Smith B J and Radom L.Calculation of proton affinities using the G2(M P2,SVP)procedure J.J Phys Chem,1995,99:6468.25 Curtiss L A,Redfern P C,Smith B J and Radom L.Gaussian-2(G2)theory:Reduced basis set requirements J.J Chem Phys,1996,104:5148.26 Bauschlicher C W and Partridge H.A modification of the Gaussian-2 approach using density functional theory J.J Chem Phys,1995,103:1788.27 M ebel A M,Morokuma K and Lin M C.Modification of the GAUSSIAN-2 theoretical model:T he use of coupled-cluster energies,densi?ty-functional geometries,and frequencies J.J Chem Phys,1995,103:7414.28 Curtissad L A,Raghavacharib K,Redferna P C,Baboula A G andPople J A.Gaussian-3 theory using coupled cluster energies J.ChemPhys Lett,1999,314:101.29 Fast P L,S?nchez M L and Truhlara D G.Multi-coefficient Gaussian-3 method for calculating potential energy syrfaces J.Chem PhysLett,1999,306:407.30 Fast,Patton L,Sanchez Maria Luz,Corchado Jose C T ruhlar Donald G.The Gaussian-2 method with proper dissociation,improved accu?40?贵州大学学报(自然科学版)?第 20 卷racy,and less cost J.J Chem Phys,1999,110:11679.31 Redferna P C,Blaudeau J P and Curtissad L A.Assessment of modified Gaussian-2(G2)and density functional theories for molecules con?taining third-row atoms Ga-Kr+J.J Phys Chem A,1997,101:8701.32 Curtiss L A,Redfern P C,Rassolov V,Kedziora G and Pople J A.J Chem Phys,2001,114:9287.33 Boys S F Bernardi F.T he calculation of small molecular interactions by the differences of separate total energies.Some procedures with re?duced error J.Mol Phys,1970,19:553.34 王一波,陶福明,潘毓刚.氢键的精确从头计算方法及其用于 NH3,H2O 和 HF 分子间氢键的研究 J.中国科学(B 辑)1995,25(10):1016-1025.35 王一波.高级量子化学从头计算法研究 H2S 和 H2O 分子间相互作用 J.中国科学(B辑).1995,25:673.36 Roger S Grev,Curtiss L Janssen,and Henry Schaefer F.?.Concerning zero-point vibrational energy corrections to electronic energies J.J Chem Phys,1991,95:5128.37 Dyke T R,Howard B J and Klemperer W.Radiofrequency and microwave spectrum of the hydrogen fluoride dimer;a nonrigid molecule J.J Chem Phys,1972,56:2442.38 Howard B J,Dyke T R and Klemperer W.The molecular beam spectrum and the structure of the hydrogen fluoride dimer J.J Chem Phys,1984,81:5417.39 Pine A S,Lafferty W J and Howard B J.Vibrational predissociation,tunneling,and rotational saturation in the HF and DF dimers J.JChem Phys,1984,81:2939.40 Wu X T,Hayes E F,Mccoy A B.Rotation vibration interactions in(HF)2:2.rotation-vibration interactions in low-lying vibrationalstates J?。J Chem Phys,1999,110:2365.41 Schwenke D W and Truhlar D G.Systematic study of basis set superposition errors in the calculated interaction energy of two HF molecules J.J Chem Phys,1985,82:2418.42 Schwenke D W and T ruhlar D G.A new potential energy surface for vibration-vibration coupling in HF-HF collisions:formulation andquantal scattering calculations J.J Chem Phys,1988,88:4800.43 Rijks W andWormer P E S.Correlated Van der Waals coefficients.II.Dimers consisting of CO,HF,H2Oand NH3 J.J Chem Phys,1989,90:6507.44 Quack M and Suhm,M A.Potenial energy surfaces,quasiadiabatic channels,rovibrational spectra,and intramolecular dynamics of(HF)2and its isotopomers from quantum Monte Carlo calculations J.J Chem Phys,1991,95:28.45 Racine S C and Davidson E R.Electron correlation contribution to the hydrogen bond in(HF)2 J.J Phys Chem,1993,97:6367.46 Bohac E J,Marshall M D and Miller R E.Initial State Effects in the Vibrational Predissociation Dynamics of Hydrogen Fluoride Dimer J.J Chem Phys,1992,96:6681.47 Peterson K A,Dunning T H.Benchmark calculations with correlated molecular wave functions.VII.Binding energy and