Interpolating moving least-squares methods for fitting potential energy surfaces: a strategy for efficient automatic data point placement in high dimensions.

An accurate and efficient method for automated molecular global potential energy surface (PES) construction and fitting is demonstrated. An interpolating moving least-squares (IMLS) method is developed with the flexibility to fit various ab initio data: (1) energies, (2) energies and gradients, or (3) energies, gradients, and Hessian data. The method is automated and flexible so that a PES can be optimally generated for trajectories, spectroscopy, or other applications. High efficiency is achieved by employing local IMLS in which fitting coefficients are stored at a limited number of expansion points, thus eliminating the need to perform weighted least-squares fits each time the potential is evaluated. An automatic point selection scheme based on the difference in two successive orders of IMLS fits is used to determine where new ab initio data need to be calculated for the most efficient fitting of the PES. A simple scan of the coordinate is shown to work well to identify these maxima in one dimension, but this search strategy scales poorly with dimension. We demonstrate the efficacy of using conjugate gradient minimizations on the difference surface to locate optimal data point placement in high dimensions. Results that are indicative of the accuracy, efficiency, and scalability are presented for a one-dimensional model potential (Morse) as well as for three-dimensional (HCN), six-dimensional (HOOH), and nine-dimensional (CH4) molecular PESs.

[1]  H. Rabitz,et al.  Reproducing kernel Hilbert space interpolation methods as a paradigm of high dimensional model representations: Application to multidimensional potential energy surface construction , 2003 .

[2]  Akio Kawano,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: applications to classical dynamics calculations. , 2004, The Journal of chemical physics.

[3]  Richard Dawes,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: computing high-density potential energy surface data from low-density ab initio data points. , 2007, The Journal of chemical physics.

[4]  Donald L. Thompson,et al.  Gradient incorporation in one-dimensional applications of interpolating moving least-squares methods for fitting potential energy surfaces , 2007 .

[5]  Joel M. Bowman,et al.  RESONANCES : BRIDGE BETWEEN SPECTROSCOPY AND DYNAMICS , 1998 .

[6]  Lawrence B. Harding,et al.  Ab initio calculations of electronic and vibrational energies of HCO and HOC , 1986 .

[7]  H. Rabitz,et al.  High Dimensional Model Representations , 2001 .

[8]  Herschel Rabitz,et al.  Efficient Implementation of High Dimensional Model Representations , 2001 .

[9]  Mark N. Gibbs,et al.  Combining ab initio computations, neural networks, and diffusion Monte Carlo: An efficient method to treat weakly bound molecules , 1996 .

[10]  H. Schlegel,et al.  Ab initio classical trajectory study of H2CO+H2 + CO dissociation , 1994 .

[11]  J. J. Soares Neto,et al.  The fitting of potential energy surfaces using neural networks. Application to the study of the photodissociation processes , 1998 .

[12]  Kersti Hermansson,et al.  Representation of Intermolecular Potential Functions by Neural Networks , 1998 .

[13]  Jonathan Tennyson,et al.  Ab initio global potential, dipole, adiabatic, and relativistic correction surfaces for the HCN-HNC system , 2001 .

[14]  Steven D. Brown,et al.  Neural network models of potential energy surfaces , 1995 .

[15]  M. A. Collins,et al.  Molecular potential energy surfaces by interpolation , 1994 .

[16]  Herschel Rabitz,et al.  Multicut‐HDMR with an application to an ionospheric model , 2004, J. Comput. Chem..

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

[18]  H. Rabitz,et al.  Random sampling-high dimensional model representation (RS-HDMR) and orthogonality of its different order component functions. , 2006, The journal of physical chemistry. A.

[19]  S. Singh,et al.  Approximation theory and spline functions , 1984 .

[20]  Michel Dupuis,et al.  Theoretical three-dimensional potential-energy surface for the reaction of Be with HF , 1983 .

[21]  T. Carrington,et al.  A nested molecule-independent neural network approach for high-quality potential fits. , 2006, The journal of physical chemistry. A.

[22]  G. Alexits.  Approximation theory , 1983 .

[23]  H. Rabitz,et al.  Practical Approaches To Construct RS-HDMR Component Functions , 2002 .

[24]  H. Rabitz,et al.  High Dimensional Model Representations Generated from Low Dimensional Data Samples. I. mp-Cut-HDMR , 2001 .

[25]  John D. Watts,et al.  Ab initio direct dynamics study of OH+HCl→Cl+H2O , 1997 .

[26]  William L. Hase,et al.  A direct dynamics study of the F + C2H4 → C2H3F + H product energy distributions , 1998 .

[27]  D. Schwenke,et al.  Vibrational energy levels for CH4 from an ab initio potential. , 2001, Spectrochimica acta. Part A, Molecular and biomolecular spectroscopy.

[28]  Donald L Thompson,et al.  Interpolating moving least-squares methods for fitting potential energy surfaces: Improving efficiency via local approximants. , 2007, The Journal of chemical physics.

[29]  David Levin,et al.  The approximation power of moving least-squares , 1998, Math. Comput..

[30]  R Komanduri,et al.  Ab initio potential-energy surfaces for complex, multichannel systems using modified novelty sampling and feedforward neural networks. , 2005, The Journal of chemical physics.

[31]  J. J. Soares Neto,et al.  The fitting of potential energy surfaces using neural networks: Application to the study of vibrational levels of H3+ , 1998 .

[32]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[33]  William L. Hase,et al.  DIRECT DYNAMICS SIMULATION OF THE LIFETIME OF TRIMETHYLENE , 1996 .

[34]  Joel M. Bowman,et al.  The adiabatic rotation approximation for rovibrational energies of many-mode systems: Description and tests of the method , 1998 .

[35]  R. Wyatt,et al.  Dynamics of molecules and chemical reactions , 1996 .

[36]  T Hollebeek,et al.  Constructing multidimensional molecular potential energy surfaces from ab initio data. , 2003, Annual review of physical chemistry.

[37]  D. W. Noid,et al.  Potential energy surfaces for macromolecules. A neural network technique , 1992 .

[38]  A. Gross,et al.  Representing high-dimensional potential-energy surfaces for reactions at surfaces by neural networks , 2004 .

[39]  Akio Kawano,et al.  Interpolating moving least-squares methods for fitting potential-energy surfaces: further improvement of efficiency via cutoff strategies. , 2006, The Journal of chemical physics.

[40]  Anthony Skjellum,et al.  A High-Performance, Portable Implementation of the MPI Message Passing Interface Standard , 1996, Parallel Comput..

[41]  Michael A. Collins,et al.  POLYATOMIC MOLECULAR POTENTIAL ENERGY SURFACES BY INTERPOLATION IN LOCAL INTERNAL COORDINATES , 1998 .

[42]  M. A. Collins,et al.  Convergence of molecular potential energy surfaces by interpolation: Application to the OH+H2→H2O+H reaction , 1995 .

[43]  Sergei Manzhos,et al.  A random-sampling high dimensional model representation neural network for building potential energy surfaces. , 2006, The Journal of chemical physics.

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