Accuracy of finite-difference harmonic frequencies in density functional theory
Abstract
Analytic Hessians are often viewed as essential for the calculation of accurate harmonic frequencies, but the implementation of analytic second derivatives is nontrivial and solution of the requisite coupled-perturbed equations engenders a sizable memory footprint for large systems, given that these equations are not required for energy and gradient calculations in density functional theory. Here, we benchmark the alternative approach to harmonic frequencies based on finite differences of analytic first derivatives, a procedure that is amenable to large-scale parallelization. Not only for absolute frequencies but also for isotopic and conformer-dependent frequency shifts in flexible molecules, we find that the finite-difference approach exhibits mean errors < 0.1 cm-1 as compared to results based on an analytic Hessian. For very small frequencies corresponding to nonbonded vibrations in noncovalent complexes (for which the harmonic approximation is questionable anyway), the finite-difference error can be larger, but even in these cases the errors can be reduced below 0.1 cm-1 by judicious choice of the displacement step size and a higher-order finite-difference approach. Finally, the surprising accuracy and robustness of the finite-difference results suggests that availability of the analytic Hessian is not so important in today's era of commodity processors that are readily available in largemore »
- Authors:
-
- The Ohio State Univ., Columbus, OH (United States)
- The Ohio State Univ., Columbus, OH (United States); Max-Planck-Inst. fur Kohlenforschung, Mulheim (Germany)
- Publication Date:
- Research Org.:
- The Ohio State Univ., Columbus, OH (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Basic Energy Sciences (BES). Chemical Sciences, Geosciences, and Biosciences Division; National Science Foundation (NSF)
- OSTI Identifier:
- 1604465
- Alternate Identifier(s):
- OSTI ID: 1401300
- Grant/Contract Number:
- SC0008850; CHE-1300603
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Computational Chemistry
- Additional Journal Information:
- Journal Volume: 38; Journal Issue: 19; Journal ID: ISSN 0192-8651
- Publisher:
- Wiley
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 74 ATOMIC AND MOLECULAR PHYSICS; finite difference; harmonic frequencies; isotopic shifts; noncovalent interactions; density functional theory
Citation Formats
Liu, Kuan-Yu, Liu, Jie, and Herbert, John M. Accuracy of finite-difference harmonic frequencies in density functional theory. United States: N. p., 2017.
Web. doi:10.1002/jcc.24811.
Liu, Kuan-Yu, Liu, Jie, & Herbert, John M. Accuracy of finite-difference harmonic frequencies in density functional theory. United States. https://doi.org/10.1002/jcc.24811
Liu, Kuan-Yu, Liu, Jie, and Herbert, John M. Wed .
"Accuracy of finite-difference harmonic frequencies in density functional theory". United States. https://doi.org/10.1002/jcc.24811. https://www.osti.gov/servlets/purl/1604465.
@article{osti_1604465,
title = {Accuracy of finite-difference harmonic frequencies in density functional theory},
author = {Liu, Kuan-Yu and Liu, Jie and Herbert, John M.},
abstractNote = {Analytic Hessians are often viewed as essential for the calculation of accurate harmonic frequencies, but the implementation of analytic second derivatives is nontrivial and solution of the requisite coupled-perturbed equations engenders a sizable memory footprint for large systems, given that these equations are not required for energy and gradient calculations in density functional theory. Here, we benchmark the alternative approach to harmonic frequencies based on finite differences of analytic first derivatives, a procedure that is amenable to large-scale parallelization. Not only for absolute frequencies but also for isotopic and conformer-dependent frequency shifts in flexible molecules, we find that the finite-difference approach exhibits mean errors < 0.1 cm-1 as compared to results based on an analytic Hessian. For very small frequencies corresponding to nonbonded vibrations in noncovalent complexes (for which the harmonic approximation is questionable anyway), the finite-difference error can be larger, but even in these cases the errors can be reduced below 0.1 cm-1 by judicious choice of the displacement step size and a higher-order finite-difference approach. Finally, the surprising accuracy and robustness of the finite-difference results suggests that availability of the analytic Hessian is not so important in today's era of commodity processors that are readily available in large numbers.},
doi = {10.1002/jcc.24811},
journal = {Journal of Computational Chemistry},
number = 19,
volume = 38,
place = {United States},
year = {Wed May 10 00:00:00 EDT 2017},
month = {Wed May 10 00:00:00 EDT 2017}
}
Web of Science
Works referenced in this record:
Many-Body Expansion with Overlapping Fragments: Analysis of Two Approaches
journal, February 2013
- Richard, Ryan M.; Herbert, John M.
- Journal of Chemical Theory and Computation, Vol. 9, Issue 3
Systematic optimization of long-range corrected hybrid density functionals
journal, February 2008
- Chai, Jeng-Da; Head-Gordon, Martin
- The Journal of Chemical Physics, Vol. 128, Issue 8
Understanding the many-body expansion for large systems. II. Accuracy considerations
journal, April 2016
- Lao, Ka Un; Liu, Kuan-Yu; Richard, Ryan M.
- The Journal of Chemical Physics, Vol. 144, Issue 16
Nonlocal van der Waals density functional: The simpler the better
journal, December 2010
- Vydrov, Oleg A.; Van Voorhis, Troy
- The Journal of Chemical Physics, Vol. 133, Issue 24
Density Functional Theory Calculations on the Mononuclear Non-Heme Iron Active Site of Hmd Hydrogenase: Role of the Internal Ligands in Tuning External Ligand Binding and Driving H 2 Heterolysis
journal, October 2010
- Dey, Abhishek
- Journal of the American Chemical Society, Vol. 132, Issue 39
ωB97X-V: A 10-parameter, range-separated hybrid, generalized gradient approximation density functional with nonlocal correlation, designed by a survival-of-the-fittest strategy
journal, January 2014
- Mardirossian, Narbe; Head-Gordon, Martin
- Physical Chemistry Chemical Physics, Vol. 16, Issue 21
General atomic and molecular electronic structure system
journal, November 1993
- Schmidt, Michael W.; Baldridge, Kim K.; Boatz, Jerry A.
- Journal of Computational Chemistry, Vol. 14, Issue 11, p. 1347-1363
Pair–Pair Approximation to the Generalized Many-Body Expansion: An Alternative to the Four-Body Expansion for ab Initio Prediction of Protein Energetics via Molecular Fragmentation
journal, January 2016
- Liu, Jie; Herbert, John M.
- Journal of Chemical Theory and Computation, Vol. 12, Issue 2
Infrared photodissociation of a water molecule from a flexible molecule-H2O complex: Rates and conformational product yields following XH stretch excitation
journal, April 2007
- Clarkson, Jasper R.; Herbert, John M.; Zwier, Timothy S.
- The Journal of Chemical Physics, Vol. 126, Issue 13
Conformation dependence of the SH and CS stretch frequencies of the cysteine residue
journal, November 1992
- Qian, Weili; Krimm, Samuel
- Biopolymers, Vol. 32, Issue 11
Identification of Isomers from Calculated Vibrational Spectra. A Density Functional Study of Tetrachlorinated Dibenzodioxins
journal, April 1995
- Rauhut, Guntram; Pulay, Peter
- Journal of the American Chemical Society, Vol. 117, Issue 14
A generalized many-body expansion and a unified view of fragment-based methods in electronic structure theory
journal, August 2012
- Richard, Ryan M.; Herbert, John M.
- The Journal of Chemical Physics, Vol. 137, Issue 6
Parallelization of analytical Hartree—Fock and density functional theory Hessian calculations. Part I: parallelization of coupled-perturbed Hartree—Fock equations
journal, June 2002
- Korambath, Prakashan P.; Kong, Jing; Furlani, Thomas R.
- Molecular Physics, Vol. 100, Issue 11
Mapping the genome of meta-generalized gradient approximation density functionals: The search for B97M-V
journal, February 2015
- Mardirossian, Narbe; Head-Gordon, Martin
- The Journal of Chemical Physics, Vol. 142, Issue 7
ω B97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation
journal, June 2016
- Mardirossian, Narbe; Head-Gordon, Martin
- The Journal of Chemical Physics, Vol. 144, Issue 21
Accuracy of Quantum Chemical Methods for Large Noncovalent Complexes
journal, July 2013
- Sedlak, Robert; Janowski, Tomasz; Pitoňák, Michal
- Journal of Chemical Theory and Computation, Vol. 9, Issue 8
The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals
journal, July 2007
- Zhao, Yan; Truhlar, Donald G.
- Theoretical Chemistry Accounts, Vol. 120, Issue 1-3
Understanding the many-body expansion for large systems. I. Precision considerations
journal, July 2014
- Richard, Ryan M.; Lao, Ka Un; Herbert, John M.
- The Journal of Chemical Physics, Vol. 141, Issue 1
Psi4: an open-source ab initio electronic structure program : Psi4: an electronic structure program
journal, October 2011
- Turney, Justin M.; Simmonett, Andrew C.; Parrish, Robert M.
- Wiley Interdisciplinary Reviews: Computational Molecular Science, Vol. 2, Issue 4
Advances in molecular quantum chemistry contained in the Q-Chem 4 program package
journal, September 2014
- Shao, Yihan; Gan, Zhengting; Epifanovsky, Evgeny
- Molecular Physics, Vol. 113, Issue 2
Isotope effects in the vibrational spectrum of the SF6 molecule
journal, April 2002
- Kolomiĭtsova, T. D.; Kondaurov, V. A.; Sedelkova, E. V.
- Optics and Spectroscopy, Vol. 92, Issue 4
NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations
journal, September 2010
- Valiev, M.; Bylaska, E. J.; Govind, N.
- Computer Physics Communications, Vol. 181, Issue 9, p. 1477-1489
Benchmark Structures and Harmonic Vibrational Frequencies Near the CCSD(T) Complete Basis Set Limit for Small Water Clusters: (H 2 O) n = 2, 3, 4, 5, 6
journal, April 2015
- Howard, J. Coleman; Tschumper, Gregory S.
- Journal of Chemical Theory and Computation, Vol. 11, Issue 5
A standard grid for density functional calculations
journal, July 1993
- Gill, Peter M. W.; Johnson, Benny G.; Pople, John A.
- Chemical Physics Letters, Vol. 209, Issue 5-6
Ground State Conformational Preferences and CH Stretch–Bend Coupling in a Model Alkoxy Chain: 1,2-Diphenoxyethane
journal, February 2013
- Buchanan, Evan G.; Sibert, Edwin L.; Zwier, Timothy S.
- The Journal of Physical Chemistry A, Vol. 117, Issue 13
A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu
journal, April 2010
- Grimme, Stefan; Antony, Jens; Ehrlich, Stephan
- The Journal of Chemical Physics, Vol. 132, Issue 15
Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections
journal, January 2008
- Chai, Jeng-Da; Head-Gordon, Martin
- Physical Chemistry Chemical Physics, Vol. 10, Issue 44
Works referencing / citing this record:
Fantasy versus reality in fragment-based quantum chemistry
journal, November 2019
- Herbert, John M.
- The Journal of Chemical Physics, Vol. 151, Issue 17
Efficient Calculation of the Negative Thermal Expansion in ZrW2O8
journal, July 2018
- Vila, Fernando D.; Hayashi, Scott T.; Rehr, John J.
- Frontiers in Chemistry, Vol. 6