Wang Xiaohong, Carter Stuart, Bowman Joel M
Department of Chemistry, Emory University , Atlanta Georgia 30322, United States ;
J Phys Chem A. 2015 Nov 25;119(47):11632-40. doi: 10.1021/acs.jpca.5b09816. Epub 2015 Nov 16.
We report vibrational self-consistent field/virtual state configuration interaction energies of nitromethane using the code MULTIMODE and a new full-dimensional potential energy surface (PES). The PES is a precise, permutationally invariant linear least-squares fit to 17,049 electronic energies, using the CCSD(T)-F12b method with HaDZ basis (cc-pVDZ basis for H atoms, and aug-cc-pVDZ basis for C, O, N atoms). Nitromethane has 15 vibrational degrees of freedom, including one that is a nearly free internal methyl torsion, which is accurately described by the PES. This torsional mode makes vibrational calculations very challenging and here we present the results of calculations without it. Nevertheless, 14-mode calculations are still challenging and can lead to very large Hamiltonian matrices. To address this issue, we apply a pruning scheme, suggested previously by Handy and Carter, that reduces the size of matrix without sacrificing accuracy in the eigenvalues. The method is briefly described here in the context of partitioning theory. A new and more efficient implementation of it, coded in the latest version of MULTIMODE program, is described. The accuracy and efficiency are demonstrated for 12-mode C2H4 and then applied to CH3NO2. Agreement with available experimental values of the CH3NO2 14 fundamentals is very good. Diffusion Monte Carlo calculations in full dimensionality are done for the zero-point energy and wavefuction. These indicate that the torisonal motion is nearly a free-rotor in this state.
我们报告了使用MULTIMODE代码和一个新的全维势能面(PES)计算得到的硝基甲烷的振动自洽场/虚态组态相互作用能。该势能面是对17049个电子能量进行精确的、具有置换不变性的线性最小二乘拟合得到的,采用了带有HaDZ基组(氢原子采用cc - pVDZ基组,碳、氧、氮原子采用aug - cc - pVDZ基组)的CCSD(T)-F12b方法。硝基甲烷有15个振动自由度,其中包括一个近似自由的内部甲基扭转振动,势能面能够准确描述该振动。这种扭转模式使得振动计算极具挑战性,在此我们展示了不考虑该扭转振动时的计算结果。尽管如此,14模式的计算仍然具有挑战性,并且可能导致非常大的哈密顿矩阵。为了解决这个问题,我们应用了一种由汉迪和卡特之前提出的剪枝方案,该方案在不牺牲本征值精度的情况下减小了矩阵的大小。本文在分区理论的背景下简要描述了该方法。还描述了在最新版本的MULTIMODE程序中编码实现的一种新的、更高效的该方法。通过对12模式的C2H4的计算展示了该方法的准确性和效率,然后将其应用于CH3NO2。与CH3NO2的14个基本振动模式的现有实验值吻合得非常好。对零点能量和波函数进行了全维的扩散蒙特卡罗计算。这些结果表明在该状态下扭转运动近似为自由转子运动。