Wang Lu, Zhang Xiao-Long, Zhai Yu, Nooijen Marcel, Li Hui
Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, China.
Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada.
J Chem Phys. 2020 Aug 7;153(5):054303. doi: 10.1063/5.0009098.
An ab initio intermolecular potential energy surface (PES) for the van der Waals complex of HO-N that explicitly incorporates the intramolecular Q bending normal mode of the HO monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster theory [CCSD(T)-F12] with an augmented correlation-consistent triple zeta basis set and an additional bond function. Analytic five-dimensional intermolecular PESs for ν(HO) = 0 and 1 are obtained by fitting to the multi-dimensional Morse/long-range potential function form. These fits to 40 890 points have the root-mean-square (rms) discrepancy of 0.88 cm for interaction energies less than 2000.0 cm. The resulting vibrationally averaged PESs provide good representations of the experimental microwave and infrared data: for microwave transitions of HO-N, the rms discrepancy is only 0.0003 cm, and for infrared transitions of the A symmetry of the HO(ν = 1 ← 0)-N, the rms discrepancy is 0.001 cm. The calculated infrared band origin shifts associated with the ν bending vibration of water are 2.210 cm and 1.323 cm for HO-N and DO-N, respectively, in good agreement with the experimental values of 2.254 cm and 1.266 cm. The benchmark tests and comparisons of the predicted spectral properties are carried out between CCSD(T)-F12a and CCSD(T)-F12b approaches.
本文提出了一种用于HO-N范德华复合物的从头算分子间势能面(PES),该势能面明确纳入了HO单体的分子内Q弯曲简正模式。电子结构计算是在显式相关耦合簇理论[CCSD(T)-F12]下进行的,采用了增强的相关一致三重ζ基组和一个额外的键函数。通过拟合多维莫尔斯/长程势函数形式,获得了ν(HO)=0和1时的解析五维分子间PES。对于相互作用能小于2000.0cm的情况,对40890个点的这些拟合的均方根(rms)偏差为0.88cm。所得的振动平均PES很好地表示了实验微波和红外数据:对于HO-N的微波跃迁,rms偏差仅为0.0003cm,对于HO(ν = 1 ← 0)-N的A对称性的红外跃迁,rms偏差为0.001cm。对于HO-N和DO-N,与水的ν弯曲振动相关的计算红外带起源位移分别为2.210cm和1.323cm,与实验值2.254cm和1.266cm吻合良好。在CCSD(T)-F12a和CCSD(T)-F12b方法之间进行了预测光谱性质的基准测试和比较。