Automatic potential energy surface generation directly from ab initio calculations using Shepard interpolation: A test calculation for the H2+H system

Toshimasa Ishida, George C Schatz

Research output: Contribution to journalArticle

56 Citations (Scopus)

Abstract

A modified version of the Shepard interpolation scheme proposed recently by Collins et al. is used to determine a potential energy surface for the H2-H reaction. Our modifications are based on preliminary calculations in which the Liu-Siegbahn-Truhlar-Horowitz surface for H3 is used to study convergence of the Shepard procedure. Included in the modifications are changes to the form of the weight functions, the coordinate sets used in the zeroth-order surface, and the transformation of the Cartesian first and second derivatives to internal coordinates. The new ab initio potential energy surface is based on calculations using the basis set of Siegbahn and Liu and second order Møller-Plesset (MP2) perturbation calculations. We have not experienced convergence problems with the self-consistent-field iteration while "growing" the surface. Furthermore, a converged fit can be obtained using only 85 ab initio calculations. This demonstrates that the Shepard interpolation scheme is a powerful candidate for automatic determination of medium quality potential surfaces for dynamical calculations.

Original languageEnglish
Pages (from-to)3558-3568
Number of pages11
JournalJournal of Chemical Physics
Volume107
Issue number9
Publication statusPublished - Sep 1 1997

Fingerprint

Potential energy surfaces
interpolation
Interpolation
potential energy
iteration
self consistent fields
Derivatives
perturbation

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

@article{b4bde46c74254970818f079a7b51d86f,
title = "Automatic potential energy surface generation directly from ab initio calculations using Shepard interpolation: A test calculation for the H2+H system",
abstract = "A modified version of the Shepard interpolation scheme proposed recently by Collins et al. is used to determine a potential energy surface for the H2-H reaction. Our modifications are based on preliminary calculations in which the Liu-Siegbahn-Truhlar-Horowitz surface for H3 is used to study convergence of the Shepard procedure. Included in the modifications are changes to the form of the weight functions, the coordinate sets used in the zeroth-order surface, and the transformation of the Cartesian first and second derivatives to internal coordinates. The new ab initio potential energy surface is based on calculations using the basis set of Siegbahn and Liu and second order M{\o}ller-Plesset (MP2) perturbation calculations. We have not experienced convergence problems with the self-consistent-field iteration while {"}growing{"} the surface. Furthermore, a converged fit can be obtained using only 85 ab initio calculations. This demonstrates that the Shepard interpolation scheme is a powerful candidate for automatic determination of medium quality potential surfaces for dynamical calculations.",
author = "Toshimasa Ishida and Schatz, {George C}",
year = "1997",
month = "9",
day = "1",
language = "English",
volume = "107",
pages = "3558--3568",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "9",

}

TY - JOUR

T1 - Automatic potential energy surface generation directly from ab initio calculations using Shepard interpolation

T2 - A test calculation for the H2+H system

AU - Ishida, Toshimasa

AU - Schatz, George C

PY - 1997/9/1

Y1 - 1997/9/1

N2 - A modified version of the Shepard interpolation scheme proposed recently by Collins et al. is used to determine a potential energy surface for the H2-H reaction. Our modifications are based on preliminary calculations in which the Liu-Siegbahn-Truhlar-Horowitz surface for H3 is used to study convergence of the Shepard procedure. Included in the modifications are changes to the form of the weight functions, the coordinate sets used in the zeroth-order surface, and the transformation of the Cartesian first and second derivatives to internal coordinates. The new ab initio potential energy surface is based on calculations using the basis set of Siegbahn and Liu and second order Møller-Plesset (MP2) perturbation calculations. We have not experienced convergence problems with the self-consistent-field iteration while "growing" the surface. Furthermore, a converged fit can be obtained using only 85 ab initio calculations. This demonstrates that the Shepard interpolation scheme is a powerful candidate for automatic determination of medium quality potential surfaces for dynamical calculations.

AB - A modified version of the Shepard interpolation scheme proposed recently by Collins et al. is used to determine a potential energy surface for the H2-H reaction. Our modifications are based on preliminary calculations in which the Liu-Siegbahn-Truhlar-Horowitz surface for H3 is used to study convergence of the Shepard procedure. Included in the modifications are changes to the form of the weight functions, the coordinate sets used in the zeroth-order surface, and the transformation of the Cartesian first and second derivatives to internal coordinates. The new ab initio potential energy surface is based on calculations using the basis set of Siegbahn and Liu and second order Møller-Plesset (MP2) perturbation calculations. We have not experienced convergence problems with the self-consistent-field iteration while "growing" the surface. Furthermore, a converged fit can be obtained using only 85 ab initio calculations. This demonstrates that the Shepard interpolation scheme is a powerful candidate for automatic determination of medium quality potential surfaces for dynamical calculations.

UR - http://www.scopus.com/inward/record.url?scp=0001201569&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001201569&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0001201569

VL - 107

SP - 3558

EP - 3568

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 9

ER -