Reduced-dimensional quantum computations for the rotational-vibrational dynamics of F(-)-CH4 and F(-)-CH2D2

J Phys Chem A. 2013 Aug 15;117(32):6975-83. doi: 10.1021/jp312160n. Epub 2013 Mar 5.

Abstract

Variational rotational-vibrational quantum chemical computations are performed for the F(-)-CH4 and F(-)-CH2D2 anion complexes using several reduced-dimensional models in a curvilinear polyspherical coordinate system and utilizing an accurate ab initio potential energy surface (PES). The implementation of the models is made practical by using the general rovibrational code GENIUSH, which constructs the complicated form of the exact rovibrational kinetic energy operator in reduced and full dimensions in any user-specified coordinates and body-fixed frames. A one-dimensional CF stretch, 1D(RCF), a two-dimensional intermolecular bend, 2D(θ,φ), and a three-dimensional intermolecular, 3D(RCF,θ,φ), rigid methane model provide vibrational energies for the low-frequency, large-amplitude modes in good agreement with full-dimensional MCTDH results for F(-)-CH4. The 2D(θ,φ) and 3D(RCF,θ,φ) four-well computations, describing equally the four possible CH-F(-) bonds, show that the ground-state tunneling splitting is less than 0.01 cm(-1). For the hydrogen-bonded CH stretching fundamental a local-mode model is found to have almost spectroscopic accuracy, whereas a harmonic frequency analysis performs poorly. The 2D(θ,φ) and 3D(RCF,θ,φ) rotational-vibrational computations on the Td-symmetric four-well PES reveal that in most cases F(-)-CH4 behaves as a semirigid C3v symmetric top. For the degenerate intermolecular bending vibrational states substantial splittings of the rigid rotor levels are observed. For F(-)-CH2D2 the rotational levels guide the assignment of the vibrational states to either F(-)-H or F(-)-D connectivity.