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.