Grid-based empirical improvement of molecular potential energy surfaces

Tamás Szidarovszky, A. Császár

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

A grid-based method designed to refine adiabatic potential energy surfaces (PES) of molecules via minimizing a suitable objective function is described. The objective function contains deviations from the reference (experimental) (ro)vibrational energy levels and is based on PES correction values determined at the grid points within a discrete-variable-representation nuclear-motion algorithm and first-order perturbation theory (PT). The proposed PES refinement technique is tested on the ground electronic state of the MgH molecule. The large number of numerical test results obtained suggest the following: (1) first-order PT is able to yield accurate correction values at the grid points representing the PES, and for practical cases there seems to be no need to go to higher orders of PT; (2) with the number of grid points greatly exceeding the number of experimental energy levels included in the refinement procedure, terms additional to the obs-calc term, including numerical first and second derivatives of the correction surface, are necessary in the objective function to arrive at a physically meaningful, smooth correction surface; (3) for a given J rotational quantum number, the corrected PES is able to reproduce experimental (ro)vibrational energies to within tenths of cm-1 if they are included in the refinement or interpolated between states that are involved in the optimization, whereas extrapolated states tend to have somewhat larger remaining discrepancies; (4) the PES refined only for the J = 0 states introduces a minor systematic error for J > 0 states, with discrepancies growing with J; (5) when the number of experimental energies included in the refinement greatly exceeds the number of grid points upon which the PES is optimized, the systematic error of treating states with different J rotational quantum numbers can be reduced and an impressive average accuracy can be achieved for all rovibrational states; and (6) in the case of quasibound (also known as resonance) rovibrational states, energies can be computed to accuracies similar to those of the bound states and excellent lifetimes (widths) can also be determined. Changes in thermochemical functions upon inclusion of quasibound states during direct summation is discussed.

Original languageEnglish
Pages (from-to)6256-6265
Number of pages10
JournalJournal of Physical Chemistry A
Volume118
Issue number32
DOIs
Publication statusPublished - Aug 14 2014

Fingerprint

Potential energy surfaces
potential energy
grids
Electron energy levels
Systematic errors
perturbation theory
systematic errors
quantum numbers
Molecules
energy levels
Electronic states
energy
Derivatives
molecules
inclusions
deviation
life (durability)
optimization

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Grid-based empirical improvement of molecular potential energy surfaces. / Szidarovszky, Tamás; Császár, A.

In: Journal of Physical Chemistry A, Vol. 118, No. 32, 14.08.2014, p. 6256-6265.

Research output: Contribution to journalArticle

@article{3e92b2dfa11548428335c05f10901e57,
title = "Grid-based empirical improvement of molecular potential energy surfaces",
abstract = "A grid-based method designed to refine adiabatic potential energy surfaces (PES) of molecules via minimizing a suitable objective function is described. The objective function contains deviations from the reference (experimental) (ro)vibrational energy levels and is based on PES correction values determined at the grid points within a discrete-variable-representation nuclear-motion algorithm and first-order perturbation theory (PT). The proposed PES refinement technique is tested on the ground electronic state of the MgH molecule. The large number of numerical test results obtained suggest the following: (1) first-order PT is able to yield accurate correction values at the grid points representing the PES, and for practical cases there seems to be no need to go to higher orders of PT; (2) with the number of grid points greatly exceeding the number of experimental energy levels included in the refinement procedure, terms additional to the obs-calc term, including numerical first and second derivatives of the correction surface, are necessary in the objective function to arrive at a physically meaningful, smooth correction surface; (3) for a given J rotational quantum number, the corrected PES is able to reproduce experimental (ro)vibrational energies to within tenths of cm-1 if they are included in the refinement or interpolated between states that are involved in the optimization, whereas extrapolated states tend to have somewhat larger remaining discrepancies; (4) the PES refined only for the J = 0 states introduces a minor systematic error for J > 0 states, with discrepancies growing with J; (5) when the number of experimental energies included in the refinement greatly exceeds the number of grid points upon which the PES is optimized, the systematic error of treating states with different J rotational quantum numbers can be reduced and an impressive average accuracy can be achieved for all rovibrational states; and (6) in the case of quasibound (also known as resonance) rovibrational states, energies can be computed to accuracies similar to those of the bound states and excellent lifetimes (widths) can also be determined. Changes in thermochemical functions upon inclusion of quasibound states during direct summation is discussed.",
author = "Tam{\'a}s Szidarovszky and A. Cs{\'a}sz{\'a}r",
year = "2014",
month = "8",
day = "14",
doi = "10.1021/jp504348f",
language = "English",
volume = "118",
pages = "6256--6265",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "32",

}

TY - JOUR

T1 - Grid-based empirical improvement of molecular potential energy surfaces

AU - Szidarovszky, Tamás

AU - Császár, A.

PY - 2014/8/14

Y1 - 2014/8/14

N2 - A grid-based method designed to refine adiabatic potential energy surfaces (PES) of molecules via minimizing a suitable objective function is described. The objective function contains deviations from the reference (experimental) (ro)vibrational energy levels and is based on PES correction values determined at the grid points within a discrete-variable-representation nuclear-motion algorithm and first-order perturbation theory (PT). The proposed PES refinement technique is tested on the ground electronic state of the MgH molecule. The large number of numerical test results obtained suggest the following: (1) first-order PT is able to yield accurate correction values at the grid points representing the PES, and for practical cases there seems to be no need to go to higher orders of PT; (2) with the number of grid points greatly exceeding the number of experimental energy levels included in the refinement procedure, terms additional to the obs-calc term, including numerical first and second derivatives of the correction surface, are necessary in the objective function to arrive at a physically meaningful, smooth correction surface; (3) for a given J rotational quantum number, the corrected PES is able to reproduce experimental (ro)vibrational energies to within tenths of cm-1 if they are included in the refinement or interpolated between states that are involved in the optimization, whereas extrapolated states tend to have somewhat larger remaining discrepancies; (4) the PES refined only for the J = 0 states introduces a minor systematic error for J > 0 states, with discrepancies growing with J; (5) when the number of experimental energies included in the refinement greatly exceeds the number of grid points upon which the PES is optimized, the systematic error of treating states with different J rotational quantum numbers can be reduced and an impressive average accuracy can be achieved for all rovibrational states; and (6) in the case of quasibound (also known as resonance) rovibrational states, energies can be computed to accuracies similar to those of the bound states and excellent lifetimes (widths) can also be determined. Changes in thermochemical functions upon inclusion of quasibound states during direct summation is discussed.

AB - A grid-based method designed to refine adiabatic potential energy surfaces (PES) of molecules via minimizing a suitable objective function is described. The objective function contains deviations from the reference (experimental) (ro)vibrational energy levels and is based on PES correction values determined at the grid points within a discrete-variable-representation nuclear-motion algorithm and first-order perturbation theory (PT). The proposed PES refinement technique is tested on the ground electronic state of the MgH molecule. The large number of numerical test results obtained suggest the following: (1) first-order PT is able to yield accurate correction values at the grid points representing the PES, and for practical cases there seems to be no need to go to higher orders of PT; (2) with the number of grid points greatly exceeding the number of experimental energy levels included in the refinement procedure, terms additional to the obs-calc term, including numerical first and second derivatives of the correction surface, are necessary in the objective function to arrive at a physically meaningful, smooth correction surface; (3) for a given J rotational quantum number, the corrected PES is able to reproduce experimental (ro)vibrational energies to within tenths of cm-1 if they are included in the refinement or interpolated between states that are involved in the optimization, whereas extrapolated states tend to have somewhat larger remaining discrepancies; (4) the PES refined only for the J = 0 states introduces a minor systematic error for J > 0 states, with discrepancies growing with J; (5) when the number of experimental energies included in the refinement greatly exceeds the number of grid points upon which the PES is optimized, the systematic error of treating states with different J rotational quantum numbers can be reduced and an impressive average accuracy can be achieved for all rovibrational states; and (6) in the case of quasibound (also known as resonance) rovibrational states, energies can be computed to accuracies similar to those of the bound states and excellent lifetimes (widths) can also be determined. Changes in thermochemical functions upon inclusion of quasibound states during direct summation is discussed.

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

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

U2 - 10.1021/jp504348f

DO - 10.1021/jp504348f

M3 - Article

AN - SCOPUS:84906229987

VL - 118

SP - 6256

EP - 6265

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 32

ER -