Theoretical Chemistry, Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands.
J Chem Phys. 2017 Aug 28;147(8):084306. doi: 10.1063/1.4990661.
We compute four-dimensional diabatic potential energy surfaces and transition dipole moment surfaces of O-O, relevant for the theoretical description of collision-induced absorption in the forbidden XΣ → aΔ and XΣ → bΣ bands at 7883 cm and 13 122 cm, respectively. We compute potentials at the multi-reference configuration interaction (MRCI) level and dipole surfaces at the MRCI and complete active space self-consistent field (CASSCF) levels of theory. Potentials and dipole surfaces are transformed to a diabatic basis using a recent multiple-property-based diabatization algorithm. We discuss the angular expansion of these surfaces, derive the symmetry constraints on the expansion coefficients, and present working equations for determining the expansion coefficients by numerical integration over the angles. We also present an interpolation scheme with exponential extrapolation to both short and large separations, which is used for representing the O-O distance dependence of the angular expansion coefficients. For the triplet ground state of the complex, the potential energy surface is in reasonable agreement with previous calculations, whereas global excited state potentials are reported here for the first time. The transition dipole moment surfaces are strongly dependent on the level of theory at which they are calculated, as is also shown here by benchmark calculations at high symmetry geometries. Therefore, ab initio calculations of the collision-induced absorption spectra cannot become quantitatively predictive unless more accurate transition dipole surfaces can be computed. This is left as an open question for method development in electronic structure theory. The calculated potential energy and transition dipole moment surfaces are employed in quantum dynamical calculations of collision-induced absorption spectra reported in Paper II [T. Karman et al., J. Chem. Phys. 147, 084307 (2017)].
我们计算了 O-O 的四元非绝热势能面和跃迁偶极矩面,这对于理论描述分别在 7883 cm 和 13122 cm 处的禁戒 XΣ→aΔ 和 XΣ→bΣ 带中碰撞诱导吸收是至关重要的。我们在多参考组态相互作用(MRCI)水平上计算势能,在 MRCI 和完全活性空间自洽场(CASSCF)水平上计算偶极子表面。使用最近的基于多属性的二分法算法,将势能和偶极子表面转换为非绝热基。我们讨论了这些表面的角扩展,推导了扩展系数的对称约束,并提出了通过在角度上进行数值积分来确定扩展系数的工作方程。我们还提出了一种具有指数外推的插值方案,用于表示角扩展系数对 O-O 距离的依赖性。对于复合物的三重基态,势能面与以前的计算结果吻合良好,而全局激发态势能则是首次报道。跃迁偶极矩面强烈依赖于计算它们的理论水平,这也通过在高对称几何形状下的基准计算得到了证明。因此,除非能够计算出更准确的跃迁偶极子表面,否则碰撞诱导吸收光谱的从头算计算不能成为定量预测。这在电子结构理论的方法发展中留下了一个悬而未决的问题。计算得到的势能和跃迁偶极矩面被用于 Paper II [T. Karman 等人,J. Chem. Phys. 147, 084307(2017)]中报道的量子动力学碰撞诱导吸收光谱的计算。