Variational vibrational calculations using high-order anharmonic force fields

Research output: Contribution to journalArticle

55 Citations (Scopus)

Abstract

A simple variational procedure, termed DOPI for discrete variable representation - Hamiltonian in orthogonal coordinates - direct product basis-iterative diagonalization, is described and applied to compute low-lying vibrational band origins (VBOs) of the triatomic systems H2O, CO 2, and N2O, employing published empirical and theoretical sextic force fields. While in these cases no difficulties arise when quartic potentials are used, the limited range of applicability of 6th-order potentials presents difficulties for the variational determination of VBOs, in particular for the higher-lying bending states. For H2O, transformation of quadratic and quartic force fields from simple bond stretching to Simons-Parr-Finlan (SPF) coordinates results in computed VBOs deviating less from experiment. This, however, does not hold for the VBOs computed from the transformed sextic force fields where the two representations provide highly similar results. While use of empirical quartic and sextic force fields result in a much better reproduction of experimental VBOs than that of ab initio force fields, especially at higher (fifth- and sixth-) order the empirical force constants, obtained through different refinement procedures, do not correspond to the associated derivatives of the potential energy surface (PES). Rotational constants characterizing low-lying vibrational states have been evaluated as expectation values using inertia tensor formulas in the Eckart and principal axis frames. Only the Eckart axes should be used for these computations and they yield accurate vibrationally averaged rotational constants.

Original languageEnglish
Pages (from-to)2411-2423
Number of pages13
JournalMolecular Physics
Volume102
Issue number23-24
DOIs
Publication statusPublished - Dec 10 2004

Fingerprint

Hamiltonians
Potential energy surfaces
Carbon Monoxide
field theory (physics)
Stretching
Tensors
Derivatives
Reproduction
Experiments
vibrational states
inertia
potential energy
tensors
products

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Variational vibrational calculations using high-order anharmonic force fields. / Czakó, G.; Furtenbacher, T.; Császár, A.; Szalay, V.

In: Molecular Physics, Vol. 102, No. 23-24, 10.12.2004, p. 2411-2423.

Research output: Contribution to journalArticle

@article{266a1f6738904c56a2b5971425bb930d,
title = "Variational vibrational calculations using high-order anharmonic force fields",
abstract = "A simple variational procedure, termed DOPI for discrete variable representation - Hamiltonian in orthogonal coordinates - direct product basis-iterative diagonalization, is described and applied to compute low-lying vibrational band origins (VBOs) of the triatomic systems H2O, CO 2, and N2O, employing published empirical and theoretical sextic force fields. While in these cases no difficulties arise when quartic potentials are used, the limited range of applicability of 6th-order potentials presents difficulties for the variational determination of VBOs, in particular for the higher-lying bending states. For H2O, transformation of quadratic and quartic force fields from simple bond stretching to Simons-Parr-Finlan (SPF) coordinates results in computed VBOs deviating less from experiment. This, however, does not hold for the VBOs computed from the transformed sextic force fields where the two representations provide highly similar results. While use of empirical quartic and sextic force fields result in a much better reproduction of experimental VBOs than that of ab initio force fields, especially at higher (fifth- and sixth-) order the empirical force constants, obtained through different refinement procedures, do not correspond to the associated derivatives of the potential energy surface (PES). Rotational constants characterizing low-lying vibrational states have been evaluated as expectation values using inertia tensor formulas in the Eckart and principal axis frames. Only the Eckart axes should be used for these computations and they yield accurate vibrationally averaged rotational constants.",
author = "G. Czak{\'o} and T. Furtenbacher and A. Cs{\'a}sz{\'a}r and V. Szalay",
year = "2004",
month = "12",
day = "10",
doi = "10.1080/0026897042000274991",
language = "English",
volume = "102",
pages = "2411--2423",
journal = "Molecular Physics",
issn = "0026-8976",
publisher = "Taylor and Francis Ltd.",
number = "23-24",

}

TY - JOUR

T1 - Variational vibrational calculations using high-order anharmonic force fields

AU - Czakó, G.

AU - Furtenbacher, T.

AU - Császár, A.

AU - Szalay, V.

PY - 2004/12/10

Y1 - 2004/12/10

N2 - A simple variational procedure, termed DOPI for discrete variable representation - Hamiltonian in orthogonal coordinates - direct product basis-iterative diagonalization, is described and applied to compute low-lying vibrational band origins (VBOs) of the triatomic systems H2O, CO 2, and N2O, employing published empirical and theoretical sextic force fields. While in these cases no difficulties arise when quartic potentials are used, the limited range of applicability of 6th-order potentials presents difficulties for the variational determination of VBOs, in particular for the higher-lying bending states. For H2O, transformation of quadratic and quartic force fields from simple bond stretching to Simons-Parr-Finlan (SPF) coordinates results in computed VBOs deviating less from experiment. This, however, does not hold for the VBOs computed from the transformed sextic force fields where the two representations provide highly similar results. While use of empirical quartic and sextic force fields result in a much better reproduction of experimental VBOs than that of ab initio force fields, especially at higher (fifth- and sixth-) order the empirical force constants, obtained through different refinement procedures, do not correspond to the associated derivatives of the potential energy surface (PES). Rotational constants characterizing low-lying vibrational states have been evaluated as expectation values using inertia tensor formulas in the Eckart and principal axis frames. Only the Eckart axes should be used for these computations and they yield accurate vibrationally averaged rotational constants.

AB - A simple variational procedure, termed DOPI for discrete variable representation - Hamiltonian in orthogonal coordinates - direct product basis-iterative diagonalization, is described and applied to compute low-lying vibrational band origins (VBOs) of the triatomic systems H2O, CO 2, and N2O, employing published empirical and theoretical sextic force fields. While in these cases no difficulties arise when quartic potentials are used, the limited range of applicability of 6th-order potentials presents difficulties for the variational determination of VBOs, in particular for the higher-lying bending states. For H2O, transformation of quadratic and quartic force fields from simple bond stretching to Simons-Parr-Finlan (SPF) coordinates results in computed VBOs deviating less from experiment. This, however, does not hold for the VBOs computed from the transformed sextic force fields where the two representations provide highly similar results. While use of empirical quartic and sextic force fields result in a much better reproduction of experimental VBOs than that of ab initio force fields, especially at higher (fifth- and sixth-) order the empirical force constants, obtained through different refinement procedures, do not correspond to the associated derivatives of the potential energy surface (PES). Rotational constants characterizing low-lying vibrational states have been evaluated as expectation values using inertia tensor formulas in the Eckart and principal axis frames. Only the Eckart axes should be used for these computations and they yield accurate vibrationally averaged rotational constants.

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

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

U2 - 10.1080/0026897042000274991

DO - 10.1080/0026897042000274991

M3 - Article

VL - 102

SP - 2411

EP - 2423

JO - Molecular Physics

JF - Molecular Physics

SN - 0026-8976

IS - 23-24

ER -