This work addresses the pathological behavior of the energetics of dimethyl sulfoxide and related sulfur-containing compounds by providing the computational benchmark energetics of R2E2 species, where R = H/CH3 and E = O/S, with bent and pyramidal geometries using state-of-the-art methodologies. These 22 geometries were fully characterized with coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)], second-order Møller–Plesset perturbation theory (MP2), and 22 density functional theory (DFT) methods with 8, 12, and 12, respectively, correlation consistent basis sets of double-, triple-, or quadruple-ζ quality. The relative energetics were determined at the MP2 and CCSD(T) complete basis set (CBS) limits using 17 basis sets up to sextuple-ζ and include augmented, tight-d, and core–valence correlation consistent basis sets. The relative energies of oxygen-/sulfur-containing compounds exhibit exceptionally slow convergence to the CBS limit with canonical methods as well as significant basis set dependence. CCSD(T) with quadruple-ζ basis sets can give qualitatively incorrect relative energies. Explicitly correlated MP2-F12 and CCSD(T)-F12 methods dramatically accelerate the convergence of the relative energies to the CBS limit for these problematic compounds. The F12 methods with a triple-ζ quality basis set give relative energies that deviate no more than 0.41 kcal mol−1 from the benchmark CBS limit. The correlation consistent Composite Approach (ccCA), ccCA-TM (TM for transition metals), and G3B3 deviated by no more than 2 kcal mol−1 from the benchmark CBS limits. Relative energies for oxygen-/sulfur-containing systems fully characterized with DFT are quite unreliable even with triple-ζ quality basis sets, and 13 out of 45 combinations fortuitously give a relative energy that is within 1 kcal mol−1 on average from the benchmark CCSD(T) CBS limit for these systems.
REFERENCES
“The G3large basis set can be downloaded from the website http://chemistry.anl.gov/compmat/g3theory.htm. The G3 total energies of the molecules in the G2/97 test set are also available from this website.
“The 6-31G(d), 6-311G(d), and 6-31G(2d f,p) basis sets use six Cartesian d-functions (6d), while the G3large basis set uses five “pure” d-functions (5d) as in 6-311G(d). Both the 6-31G(2d f,p) and G3large basis sets use a “pure” 7 f set.”