We have obtained uniform frequency scaling factors λ(harm) (for harmonic frequencies), λ(fund) (for fundamentals), and λ(ZPVE) (for zero-point vibrational energies (ZPVEs)) for the Weigend-Ahlrichs and other selected basis sets for MP2, SCS-MP2, and a variety of DFT functionals including double hybrids. For selected levels of theory, we have also obtained scaling factors for true anharmonic fundamentals and ZPVEs obtained from quartic force fields. For harmonic frequencies, the double hybrids B2PLYP, B2GP-PLYP, and DSD-PBEP86 clearly yield the best performance at RMSD = 10-12 cm(-1) for def2-TZVP and larger basis sets, compared to 5 cm(-1) at the CCSD(T) basis set limit. For ZPVEs, again, the double hybrids are the best performers, reaching root-mean-square deviations (RMSDs) as low as 0.05 kcal/mol, but even mainstream functionals like B3LYP can get down to 0.10 kcal/mol. Explicitly anharmonic ZPVEs only are marginally more accurate. For fundamentals, however, simple uniform scaling is clearly inadequate.

The S66x8 dataset for noncovalent interactions of biochemical relevance has been re-examined by means of MP2-F12 and CCSD(F12*)(T) methods. We deem our revised benchmark data to be reliable to about 0.05 kcal mol(-1) RMS. Most levels of DFT perform quite poorly in the absence of dispersion corrections: somewhat surprisingly, that is even the case for the double hybrids and for dRPA75. Analysis of optimized D3BJ parameters reveals that the main benefit of dRPA75 and DSD double hybrids alike is the treatment of midrange dispersion. dRPA75-D3BJ is the best performer overall at RMSD = 0.10 kcal mol(-1). The nonlocal VV10 dispersion functional is especially beneficial for the double hybrids, particularly in DSD-PBEP86-NL (RMSD = 0.12 kcal mol(-1)). Other recommended dispersion-corrected functionals with favorable price/performance ratios are ωB97X-V, and, surprisingly, B3LYP-D3BJ and BLYP-D3BJ (RMSDs of 0.23, 0.20 and 0.23 kcal mol(-1), respectively). Without dispersion correction (but parametrized for midrange interactions) M06-2X has the lead (RMSD = 0.45 kcal mol(-1)). A collection of three energy-based diagnostics yields similar information to an SAPT analysis about the nature of the noncovalent interaction. Two of those are the percentages of Hartree-Fock and of post-MP2 correlation effects in the interaction energy; the third, CSPI = [IE - IE]/[IE + IE] or its derived quantity DEBC = CSPI/(1 + CSPI(2))(1/2), describes the character of the MP2 correlation contribution, ranging from 0 (purely dispersion) to 1 (purely other effects). In addition, we propose an improved, parameter-free scaling for the (T) contribution based on the Ecorr[CCSD-F12b]/Ecorr[CCSD] and Ecorr[CCSD(F12*)]/Ecorr[CCSD] ratios. For Hartree-Fock and conventional DFT calculations, full counterpoise generally yields the fastest basis set convergence, while for double hybrids, half-counterpoise yields faster convergence, as previously established for correlated ab initio methods.

The basis set convergence of explicitly correlated ab initio methods, when applied to noncovalent interactions, has been considered in the presence (and absence) of Boys-Bernardi counterpoise corrections, as well as using "half-counterpoise" (the average of raw and counterpoise-corrected values) as recently advocated in this journal [Burns, L. A.; Marshall, M. S.; Sherrill, C. D. J. Chem. Theory Comput. 2014, 10, 49-57]. Reference results were obtained using basis sets so large that BSSE (basis set superposition error) can be shown to be negligible. For the HF+CABS component, full counterpoise unequivocally exhibits the fastest basis set convergence. However, at the MP2-F12 and CCSD(T*)-F12b levels, surprisingly good uncorrected results can be obtained with small basis sets like cc-pVDZ-F12, owing to error compensation between basis set superposition error (which overbinds) and intrinsic basis set insufficiency (which underbinds). For intermediate sets like cc-pVTZ-F12, "half-half" averages work best, while for large basis sets like cc-pVQZ-F12, full counterpoise may be preferred but BSSE in uncorrected values is tolerably small for most purposes. A composite scheme in which CCSD(T)-MP2 "high level corrections" obtained at the CCSD(T*)-F12b/cc-pVDZ-F12 level are combined with "half-counterpoise" MP2-F12/cc-pVTZ-F12 interaction energies yields surprisingly good performance for standard benchmark sets like S22 and S66.

The results of harmonic and anharmonic frequency calculations on a guanine-cytosine complex with an enolic structure (a tautomeric form with cytosine in the enol form and with a hydrogen at the 7-position on guanine) are presented and compared to gas-phase IR-UV double resonance spectral data. Harmonic frequencies were obtained at the RI-MP2/cc-pVDZ, RI-MP2/TZVPP, and semiempirical PM3 levels of electronic structure theory. Anharmonic frequencies were obtained by the CC-VSCF method with improved PM3 potential surfaces; the improved PM3 potential surfaces are obtained from standard PM3 theory by coordinate scaling such that the improved PM3 harmonic frequencies are the same as those computed at the RI-MP2/cc-pVDZ level. Comparison of the data with experimental results indicates that the average absolute percentage deviation for the methods is 2.6% for harmonic RI-MP2/cc-pVDZ (3.0% with the inclusion of a 0.956 scaling factor that compensates for anharmonicity), 2.5% for harmonic RI-MP2/TZVPP (2.9% with a 0.956 anharmonicity factor included), and 2.3% for adapted PM3 CC-VSCF; the empirical scaling factor for the ab initio harmonic calculations improves the stretching frequencies but decreases the accuracy of the other mode frequencies. The agreement with experiment supports the adequacy of the improved PM3 potentials for describing the anharmonic force field of the G...C base pair in the spectroscopically probed region. These results may be useful for the prediction of the pathways of vibrational energy flow upon excitation of this system. The anharmonic calculations indicate that anharmonicity along single mode coordinates can be significant for simple stretching modes. For several other cases, coupling between different vibrational modes provides the main contribution to anharmonicity. Examples of strongly anharmonically coupled modes are the symmetric stretch and group torsion of the hydrogen-bonded NH2 group on guanine, the OH stretch and torsion of the enol group on cytosine, and the NH stretch and NH out-of-plane bend of the non-hydrogen-bonded NH group on guanine.

No abstract

The anharmonic vibrational spectra of α-D-glucose, β-D-glucose, and sucrose are computed by the vibrational self-consistent field (VSCF) method, using potential energy surfaces from electronic structure theory, for the lowest energy conformers that correspond to the gas phase and to the crystalline phase, respectively. The results are compared with ultraviolet-infrared (UV-IR) spectra of phenyl β-D-glucopyranoside in a molecular beam, with literature results for sugars in matrices and with new experimental data for the crystalline state. Car-Parrinello dynamics simulations are also used to study temperature effects on the spectra of α-D-glucose and β-D-glucose and to predict their vibrational spectra at 50, 150, and 300 K. The effects of temperature on the spectral features are analyzed and compared with results of the VSCF calculations conducted at 0 K. The main results include: (i) new potential surfaces, constructed from Hartree-Fock, adjusted to fit harmonic frequencies from Møller-Plesset (MP2) calculations, that give very good agreement with gas phase, matrix, and solid state spectra; (ii) computed infrared spectra of the crystalline solid of α-glucose, which are substantially improved by including mimic groups that represent the effect of the solid environment on the sugar; and (iii) identification of a small number of combination-mode transitions, which are predicted to be strong enough for experimental observation. The results are used to assess the role of anharmonic effects in the spectra of the sugars in isolation and in the solid state and to discuss the spectroscopic accuracy of potentials from different electronic structure methods.

The experimental mid-and far-IR spectra of six conformers of phenylalanine in the gas phase are presented. The experimental spectra are compared to spectra calculated at the B3LYP and at the MP2 level. The differences between B3LYP and MP2 IR spectra are found to be small. The agreement between experiment and theory is generally found to be very good, however strong discrepancies exist when -NH 2 out-of-plane vibrations are involved. The relative energies of the minima as well as of some transition states connecting the minima are explored at the CCSD(T) level. Most transition states are found to be less than 2000 cm À1 above the lowest energy structure. A simple model to describe the observed conformer abundances based on quasiequilibria near the barriers is presented and it appears to describe the experimental observation reasonably well. In addition, the vibrations of one of the conformers are investigated using the correlation-corrected vibrational self-consistent field method.

No abstract