DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States

Abstract

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV–vis radiation.

Picking up Good Vibrations through Quartic Force Fields and Vibrational Perturbation Theory

“…A similar order-of-magnitude savings as the hybrid F12-TcCR+DZ QFF has been shown for combining F12-TcCR harmonics with density functional theory (DFT) cubic and quatic terms . Most DFT functionals behave the same, but PBE0 typically provides the most consistent results .…”

“…A similar order-of-magnitude savings as the hybrid F12-TcCR+DZ QFF has been shown for combining F12-TcCR harmonics with density functional theory (DFT) cubic and quatic terms . Most DFT functionals behave the same, but PBE0 typically provides the most consistent results . However, as the size of the molecule in question increases, the number of harmonic points in the QFF quickly decreases relative to the whole total number of points. , Hence, these hybrid methods will provide better scaling as the size of (number of atoms within) the molecule increases.…”