Quantum and classical studies of vibrational motion of CH5+ on a global potential energy surface obtained from a novel ab initio direct dynamics approach

We report a full dimensional, ab initio based potential energy surface for CH(5) (+). The ab initio electronic energies and gradients are obtained in direct-dynamics calculations using second-order M?ller-Plesset perturbation theory with the correlation consistent polarized valence triple zeta basis. The potential energy and the dipole moment surfaces are fit using novel procedures that ensure the full permutational symmetry of the system. The fitted potential energy surface is tested by comparing it against additional electronic energy calculations and by comparing normal mode frequencies at the three lowest-lying stationary points obtained from the fit against ab initio ones. Well-converged diffusion Monte Carlo zero-point energies, rotational constants, and projections along the CH and HH bond lengths and the tunneling coordinates are presented and compared with the corresponding harmonic oscillator and standard classical molecular dynamics ones. The delocalization of the wave function is analyzed through comparison of the CH(5) (+) distributions with those obtained when all of the hydrogen atoms are replaced by (2)H and (3)H. The classical dipole correlation function is examined as a function of the total energy. This provides a further probe of the delocalization of CH(5) (+).