Interpolating moving least-squares methods for fitting potential energy surfaces: applications to classical dynamics calculations.

As a continuation of our efforts to develop efficient and accurate interpolating moving least-squares (IMLS) methods for generating potential energy surfaces, we carry out classical trajectories and compute kinetics properties on higher degree IMLS surfaces. In this study, we have investigated the choice of coordinate system, the range of points (i.e., the cutoff radius) used in fitting, and strategies for selections of data points and basis elements. We illustrate and test the method by applying it to hydrogen peroxide (HOOH). In particular, reaction rates for the O-O bond breaking in HOOH are calculated on fitted surfaces using the classical trajectory approach to test the accuracy of the IMLS method for providing potentials for dynamics calculations.

[1]  Akio Kawano,et al.  Improving the accuracy of interpolated potential energy surfaces by using an analytical zeroth-order potential function. , 2004, The Journal of chemical physics.

[2]  D. Thompson,et al.  A classical trajectory study of bond dissociation in HONO and HOOH , 2003 .

[3]  Donald L. Thompson,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: Detailed analysis of one-dimensional applications , 2003 .

[4]  George C. Schatz,et al.  A local interpolation scheme using no derivatives in potential sampling: Application to O(1D) + H2 system , 2003, J. Comput. Chem..

[5]  G. G. Maisuradze,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: Illustrative approaches and applications , 2003 .

[6]  B. Kuhn,et al.  State-to-state studies of intramolecular energy transfer in highly excited HOOH(D): Dependencies on vibrational and rotational excitation , 2000 .

[7]  George C. Schatz,et al.  A local interpolation scheme using no derivatives in quantum-chemical calculations , 1999 .

[8]  B. Kuhn,et al.  A new six-dimensional analytical potential up to chemically significant energies for the electronic ground state of hydrogen peroxide , 1999 .

[9]  Herschel Rabitz,et al.  A global A-state potential surface for H2O: Influence of excited states on the O(1D)+H2 reaction , 1997 .

[10]  Michael A. Collins,et al.  Molecular Potential Energy Surfaces by Interpolation , 1994, International Conference on Computational Science.

[11]  G. Schatz,et al.  A coupled channel study of HN2 unimolecular decay based on a global ab initio potential surface , 1991 .

[12]  H. W. Schranz,et al.  An efficient microcanonical sampling procedure for molecular systems , 1991 .

[13]  G. Schatz,et al.  State‐selective studies of T→R, V energy transfer: The H+CO system , 1988 .

[14]  J. Murrell,et al.  Molecular Potential Energy Functions , 1985 .

[15]  Anders Öberg,et al.  Algorithms for approximation of invariant measures for IFS , 2005 .

[16]  A. Laganá,et al.  Reaction and molecular dynamics : proceedings of the European School on Computational Chemistry, Perugia, Italy, July (1999) , 2000 .

[17]  P. Knowles,et al.  Coupled ab initio potential energy surfaces for the reaction Cl(2P)+HCl→ClH+Cl(2P) , 1999 .

[18]  Michael A. Collins,et al.  Molecular potential-energy surfaces by interpolation: Furtherrefinements , 1997 .