Estimation, Computation, and Experimental Correction of Molecular Zero-Point Vibrational Energies

G. Csonka, Adrienn Ruzsinszky, John P. Perdew

Research output: Contribution to journalArticle

32 Citations (Scopus)

Abstract

For accurate thermochemical tests of electronic structure theory, accurate true anharmonic zero-point vibrational energies ZPVE true are needed. We discuss several possibilities to extract this information for molecules from density functional or wave function calculations and/or available experimental data: (1) Empirical universal scaling of density-functional-calculated harmonic ZPVE harms, where we find that polyatomics require smaller scaling factors than diatomics. (2) Direct density-functional calculation by anharmonic second-order perturbation theory PT2. (3) Weighted averages of harmonic ZPVE harm and fundamental ZPVE fund (from fundamental vibrational transition frequencies), with weights (3/4, 1/4) for diatomics and (5/8,3/8) for polyatomics. (4) Experimental correction of the PT2 harmonic contribution, i.e., the estimate ZPVE PT2 true + (ZPVE expt fund -ZPVE PT2 fund) for ZPVE true. The (5/8,3/8) average of method 3 and the additive correction of method 4 have been proposed here. For our database of experimental ZPVE true, consisting of 27 diatomics and 8 polyatomics, we find that methods 1 and 2, applied to the popular B3LYP and the nonempirical PBE and TPSS functional and their one-parameter hybrids, yield polyatomic errors on the order of 0.1 kcal/mol. Larger errors are expected for molecules larger than those in our database. Method 3 yields errors on the order of 0.02 kcal/mol, but requires very accurate (e.g., experimental, coupled cluster, or best-performing density functional) input harmonic ZPVE harm. Method 4 is the best-founded one that meets the requirements of high accuracy and practicality, requiring as experimental input only the highly accurate and widely available ZPVE expt fund and producing errors on the order of 0.05 kcal/mol that are relatively independent of functional and basis set. As a part of our study, we also test the ability of the density functionals to predict accurate equilibrium bond lengths and angles for a data set of 21 mostly polyatomic molecules (since all calculated ZPVEs are evaluated at the correspondingly calculated molecular geometries).

Original languageEnglish
Pages (from-to)6779-6789
Number of pages11
JournalJournal of Physical Chemistry A
Volume109
Issue number30
DOIs
Publication statusPublished - Aug 4 2005

Fingerprint

harmonics
Molecules
energy
Bond length
scaling
Wave functions
Electronic structure
Density functional theory
polyatomic molecules
functionals
molecules
perturbation theory
Geometry
wave functions
electronic structure
requirements
estimates
geometry

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Estimation, Computation, and Experimental Correction of Molecular Zero-Point Vibrational Energies. / Csonka, G.; Ruzsinszky, Adrienn; Perdew, John P.

In: Journal of Physical Chemistry A, Vol. 109, No. 30, 04.08.2005, p. 6779-6789.

Research output: Contribution to journalArticle

@article{71dbd5a9d4a94bd0a584d7ba8e2a1e37,
title = "Estimation, Computation, and Experimental Correction of Molecular Zero-Point Vibrational Energies",
abstract = "For accurate thermochemical tests of electronic structure theory, accurate true anharmonic zero-point vibrational energies ZPVE true are needed. We discuss several possibilities to extract this information for molecules from density functional or wave function calculations and/or available experimental data: (1) Empirical universal scaling of density-functional-calculated harmonic ZPVE harms, where we find that polyatomics require smaller scaling factors than diatomics. (2) Direct density-functional calculation by anharmonic second-order perturbation theory PT2. (3) Weighted averages of harmonic ZPVE harm and fundamental ZPVE fund (from fundamental vibrational transition frequencies), with weights (3/4, 1/4) for diatomics and (5/8,3/8) for polyatomics. (4) Experimental correction of the PT2 harmonic contribution, i.e., the estimate ZPVE PT2 true + (ZPVE expt fund -ZPVE PT2 fund) for ZPVE true. The (5/8,3/8) average of method 3 and the additive correction of method 4 have been proposed here. For our database of experimental ZPVE true, consisting of 27 diatomics and 8 polyatomics, we find that methods 1 and 2, applied to the popular B3LYP and the nonempirical PBE and TPSS functional and their one-parameter hybrids, yield polyatomic errors on the order of 0.1 kcal/mol. Larger errors are expected for molecules larger than those in our database. Method 3 yields errors on the order of 0.02 kcal/mol, but requires very accurate (e.g., experimental, coupled cluster, or best-performing density functional) input harmonic ZPVE harm. Method 4 is the best-founded one that meets the requirements of high accuracy and practicality, requiring as experimental input only the highly accurate and widely available ZPVE expt fund and producing errors on the order of 0.05 kcal/mol that are relatively independent of functional and basis set. As a part of our study, we also test the ability of the density functionals to predict accurate equilibrium bond lengths and angles for a data set of 21 mostly polyatomic molecules (since all calculated ZPVEs are evaluated at the correspondingly calculated molecular geometries).",
author = "G. Csonka and Adrienn Ruzsinszky and Perdew, {John P.}",
year = "2005",
month = "8",
day = "4",
doi = "10.1021/jp0519464",
language = "English",
volume = "109",
pages = "6779--6789",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "30",

}

TY - JOUR

T1 - Estimation, Computation, and Experimental Correction of Molecular Zero-Point Vibrational Energies

AU - Csonka, G.

AU - Ruzsinszky, Adrienn

AU - Perdew, John P.

PY - 2005/8/4

Y1 - 2005/8/4

N2 - For accurate thermochemical tests of electronic structure theory, accurate true anharmonic zero-point vibrational energies ZPVE true are needed. We discuss several possibilities to extract this information for molecules from density functional or wave function calculations and/or available experimental data: (1) Empirical universal scaling of density-functional-calculated harmonic ZPVE harms, where we find that polyatomics require smaller scaling factors than diatomics. (2) Direct density-functional calculation by anharmonic second-order perturbation theory PT2. (3) Weighted averages of harmonic ZPVE harm and fundamental ZPVE fund (from fundamental vibrational transition frequencies), with weights (3/4, 1/4) for diatomics and (5/8,3/8) for polyatomics. (4) Experimental correction of the PT2 harmonic contribution, i.e., the estimate ZPVE PT2 true + (ZPVE expt fund -ZPVE PT2 fund) for ZPVE true. The (5/8,3/8) average of method 3 and the additive correction of method 4 have been proposed here. For our database of experimental ZPVE true, consisting of 27 diatomics and 8 polyatomics, we find that methods 1 and 2, applied to the popular B3LYP and the nonempirical PBE and TPSS functional and their one-parameter hybrids, yield polyatomic errors on the order of 0.1 kcal/mol. Larger errors are expected for molecules larger than those in our database. Method 3 yields errors on the order of 0.02 kcal/mol, but requires very accurate (e.g., experimental, coupled cluster, or best-performing density functional) input harmonic ZPVE harm. Method 4 is the best-founded one that meets the requirements of high accuracy and practicality, requiring as experimental input only the highly accurate and widely available ZPVE expt fund and producing errors on the order of 0.05 kcal/mol that are relatively independent of functional and basis set. As a part of our study, we also test the ability of the density functionals to predict accurate equilibrium bond lengths and angles for a data set of 21 mostly polyatomic molecules (since all calculated ZPVEs are evaluated at the correspondingly calculated molecular geometries).

AB - For accurate thermochemical tests of electronic structure theory, accurate true anharmonic zero-point vibrational energies ZPVE true are needed. We discuss several possibilities to extract this information for molecules from density functional or wave function calculations and/or available experimental data: (1) Empirical universal scaling of density-functional-calculated harmonic ZPVE harms, where we find that polyatomics require smaller scaling factors than diatomics. (2) Direct density-functional calculation by anharmonic second-order perturbation theory PT2. (3) Weighted averages of harmonic ZPVE harm and fundamental ZPVE fund (from fundamental vibrational transition frequencies), with weights (3/4, 1/4) for diatomics and (5/8,3/8) for polyatomics. (4) Experimental correction of the PT2 harmonic contribution, i.e., the estimate ZPVE PT2 true + (ZPVE expt fund -ZPVE PT2 fund) for ZPVE true. The (5/8,3/8) average of method 3 and the additive correction of method 4 have been proposed here. For our database of experimental ZPVE true, consisting of 27 diatomics and 8 polyatomics, we find that methods 1 and 2, applied to the popular B3LYP and the nonempirical PBE and TPSS functional and their one-parameter hybrids, yield polyatomic errors on the order of 0.1 kcal/mol. Larger errors are expected for molecules larger than those in our database. Method 3 yields errors on the order of 0.02 kcal/mol, but requires very accurate (e.g., experimental, coupled cluster, or best-performing density functional) input harmonic ZPVE harm. Method 4 is the best-founded one that meets the requirements of high accuracy and practicality, requiring as experimental input only the highly accurate and widely available ZPVE expt fund and producing errors on the order of 0.05 kcal/mol that are relatively independent of functional and basis set. As a part of our study, we also test the ability of the density functionals to predict accurate equilibrium bond lengths and angles for a data set of 21 mostly polyatomic molecules (since all calculated ZPVEs are evaluated at the correspondingly calculated molecular geometries).

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

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

U2 - 10.1021/jp0519464

DO - 10.1021/jp0519464

M3 - Article

VL - 109

SP - 6779

EP - 6789

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 30

ER -