A recently suggested constrained quasiclassical trajectory (c-QCT) method for avoiding the zero-point leak in the water dimer  is applied to the water trimer, employing an ab initio full-dimensional potential energy surface. We demonstrate the failure of the standard/unconstrained QCT method for (H 2O)3 and show the utility of c-QCT dynamics. In addition, standard classical molecular dynamics and c-QCT dynamics are contrasted for the time-dependence of isomerizations between the multiple global and local minima as well as radial distribution functions are obtained at low temperature and at 300 K. Results from these calculations are compared with rigorous quantum path integral Monte Carlo calculations.
ASJC Scopus subject areas
- Physics and Astronomy(all)
- Physical and Theoretical Chemistry