A dual-level Shepard interpolation method for generating potential energy surfaces for dynamics calculations

Nguyen, Kiet A.; Rossi, Ivan; Truhlar, Donald G.
October 1995
Journal of Chemical Physics;10/1/1995, Vol. 103 Issue 13, p5522
Academic Journal
We present a new dual-level approach to representing potential energy surfaces in which a very small number of high-level electronic structure calculations are combined with a lower-level global surface, e.g., one defined implicitly by neglect-of-diatomic-differential-overlap calculations with specific reaction parameters, to generate the potential at any geometry where it may be needed. We interpolate the potential energy surface with a small number of accurate data points (the higher level) that are placed along the reaction path by using information on the global shape of the potential from less accurate calculations (the lower level). We confirm the findings of Ischtwan and Collins on the usefulness of single-level schemes including Hessians, and we delineate the regime of usefulness of single-level schemes based on gradients or even single-point energies. Furthermore we find that dual-level interpolation can offer cost savings over single-level schemes, and dual-level methods employing Hessians, gradients, or even only simple energy evaluations can yield reasonable potential energy surfaces with relatively low cost, with the potentials being more accurate along the reaction path. For all methods considered in this paper the accuracy of the interpolation for our test cases is lower when the potentials at points significantly removed from the reaction path are predicted from data that lie entirely on the reaction path. © 1995 American Institute of Physics.


Related Articles

  • Molecular potential energy surfaces by interpolation. Ischtwan, Josef; Collins, Michael A. // Journal of Chemical Physics;6/1/1994, Vol. 100 Issue 11, p8080 

    A moving interpolation technique which provides an accurate representation of potential energy surfaces for polyatomic molecules is presented. The method uses the ab initio energy, energy gradient, and second derivatives calculated at dynamically important configurations. The interpolant of the...

  • Dissociative dynamics of spin-triplet and spin-singlet O2 on Ag(100). Alducin, M.; Busnengo, H. F.; Muiño, R. Díez // Journal of Chemical Physics;12/14/2008, Vol. 129 Issue 22, p224702 

    We study the dissociative dynamics of O2 molecules on the Ag(100) surface. Initially, the impinging molecules are either in the spin-triplet ground state or in the spin-singlet excited state. The molecule-surface interaction is obtained in each case by constructing the six-dimensional potential...

  • Neural network models of potential energy surfaces. Blank, Thomas B.; Brown, Steven D.; Calhoun, August W.; Doren, Douglas J. // Journal of Chemical Physics;9/8/1995, Vol. 103 Issue 10, p4129 

    Neural networks provide an efficient, general interpolation method for nonlinear functions of several variables. This paper describes the use of feed-forward neural networks to model global properties of potential energy surfaces from information available at a limited number of configurations....

  • The utility of higher order derivatives in constructing molecular potential energy surfaces by interpolation. Jordan, Meredith J. T.; Thompson, Keiran C.; Collins, Michael A. // Journal of Chemical Physics;12/8/1995, Vol. 103 Issue 22, p9669 

    In this paper we evaluate the use of higher order derivatives in the construction of an interpolated potential energy surface for the OH+H2→H2O+H reaction. The surface construction involves interpolating between local Taylor expansions about a set of known data points. We examine the use...

  • Potential energy surface for unimolecular dissociation and rearrangement reactions of the ground electronic state of HFCO. Kamiya, Kenshu; Morokuma, Keiji // Journal of Chemical Physics;6/1/1991, Vol. 94 Issue 11, p7287 

    The potential energy surface of the HFCO molecule in its electronic ground state has been investigated with ab initio method, at levels up to MP4(SDTQ)/6-311G**//MP2/6-31G*. At the highest level, the barrier height for molecular dissociation (HFCO→HF+CO) was calculated to be 46.9 kcal/mol...

  • Quasiclassical trajectory study of the dynamics of the H+N[sub 2]O reaction on a new potential energy surface. Castillo, J. F.; Collins, M. A.; Aoiz, F. J.; Bañares, L. // Journal of Chemical Physics;4/22/2003, Vol. 118 Issue 16, p7303 

    A new ab initio potential energy surface (PES) for the H+N[sub 2]O→OH+N[sub 2] reaction has been constructed using the GR OW package of Collins and co-workers. The ab initio calculations have been done using the Becke three-parameter nonlocal exchange functional with the nonlocal...

  • Ab initio ground potential energy surface and quasiclassical trajectory study of the... Gonzalez, Miguel; Hernando, Jordi; Banos, Irene; Sayos, R. // Journal of Chemical Physics;11/15/1999, Vol. 111 Issue 19, p8913 

    Discusses the ab initio study of the ground potential energy surface (PES) and quasiclassical trajectory of the oxygen atom reaction dynamics using the second and fourth order Moller-Plesset methods with a large basis set. Analysis of the microscopic reaction mechanism; Root-mean-square...

  • Efficient potential energy surfaces from partially filled ab initio data over arbitrarily shaped regions. Hollebeek, Timothy; Ho, Tak-San; Rabitz, Herschel // Journal of Chemical Physics;3/1/2001, Vol. 114 Issue 9 

    The reproducing kernel Hilbert space (RKHS) method has been previously shown to be accurate and efficient in the construction of potential energy surfaces (PES) by interpolating fully gridded high level ab initio data. This paper extends the RKHS method to handle partially filled data calculated...

  • Hyperspherical coordinates for molecular dynamics by the method of trees and the mapping of potential energy surfaces for triatomic systems. Aquilanti, Vincenzo; Cavalli, Simonetta; Grossi, Gaia // Journal of Chemical Physics;8/1/1986, Vol. 85 Issue 3, p1362 

    Some results on hyperspherical coordinates and harmonics for the representation of the many-body problem are presented, extensive use being made of the method of trees. Properties of these trees are examined: a lemma on the simplification of trees possessing a particular symmetry is proven, and...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics