A multiple coherent states implementation of the semiclassical approximation is introduced and employed to obtain the power spectra with a few classical trajectories. The method is integrated with the time-averaging semiclassical initial value representation to successfully reproduce anharmonicity and Fermi resonance splittings at a level of accuracy comparable to semiclassical simulations of thousands of trajectories. The method is tested on two different model systems with analytical potentials and implemented in conjunction with the first-principles molecular dynamics scheme to obtain the power spectrum for the carbon dioxide molecule.
View Article and Find Full Text PDFIn this work, we explore the use of the semiclassical initial value representation (SC-IVR) method with first-principles electronic structure approaches to carry out classical molecular dynamics. The proposed approach can extract the vibrational power spectrum of carbon dioxide from a single trajectory providing numerical results that agree with experiment and quantum calculations. The computational demands of the method are comparable to those of classical single-trajectory calculations, while describing uniquely quantum features such as the zero-point energy and Fermi resonances.
View Article and Find Full Text PDFPhys Chem Chem Phys
October 2006
We report an ab initio study of the van der Waals region of the O(3P)-H2 potential energy surface based on RCCSD(T) calculations with an aug-cc-pVQZ basis supplemented by bond functions. In addition, an open-shell implementation of symmetry-adapted perturbation theory (SAPT) is used to corroborate the RCCSD(T) calculations and to investigate the relative magnitudes of the various contributions to the van der Waals interaction. We also investigate the effect of the spin-orbit coupling on the position and depth of the van der Waals well.
View Article and Find Full Text PDFWe report state-to-state cross sections and thermal rate constants for vibrational and rotational relaxation of OH(2pi) by collision with H atoms. The cross sections are calculated by the coupled-states (CS) statistical method including the full open-shell character of the OH + H system. Four potential energy surfaces (PESs) ((1,3)A' and (1,3)A'') describe the interaction of OH(X2pi) with H atoms.
View Article and Find Full Text PDFWe report state-to-state and overall thermal rate constants for the isotope exchange reaction D((2)S)+OH((2)Pi)-->OD((2)Pi)+H((2)S) for 0 K