Parameterization of a B3LYP Specific Correction for Noncovalent Interactions and Basis Set Superposition Error on a Gigantic Data Set of CCSD(T) Quality Noncovalent Interaction Energies

A vast number of noncovalent interaction energies at the counterpoise corrected CCSD(T) level have been collected from the literature to build a diverse new data set. The whole data set, which consists of 2027 CCSD(T) energies, includes most of the published data at this level. A large subset of the data was then used to train a novel, B3LYP specific, empirical correction scheme for noncovalent interactions and basis set superposition error (abbreviated as B3LYP-MM). Results obtained with our new correction scheme were directly compared to benchmark results obtained with B3LYP-D3 and M06-2X (two popular density functionals designed specifically to accurately model noncovalent interactions). For noncovalent complexes dominated by dispersion or dipole-dipole interactions, all three tested methods give accurate results with the medium-sized aug-cc-pVDZ basis set with MUEs of 0.27 (B3LYP-MM), 0.32 (B3LYP-D3), and 0.47 kcal/mol (M06-2X) (with explicit counterpoise corrections). These results validate both B3LYP-D3 and M06-2X for interactions of this type using a much larger data set than was presented in prior work. However, our new dispersion correction scheme shows some clear advantages for dispersion and dipole-dipole dominated complexes with the small LACVP* basis set, which is very popular in use due to its low associated computational cost: The MUE for B3LYP-MM with the LACVP* basis set for this subset of complexes (without explicit counterpoise corrections) is only 0.28 kcal/mol, compared to 0.65 kcal/mol for M06-2X or 1.16 kcal/mol for B3LYP-D3. Additionally, our new correction scheme also shows major improvements in accuracy for hydrogen-bonded systems and for systems involving ionic interactions, for example, cation-pi interactions. Compared to B3LYP-D3 and M06-2X, we also find that our new B3LYP-MM correction scheme gives results of higher or equal accuracy for a large data set of conformer energies of di- and tripeptides, sugars, and cysteine.