Extension of the B3LYP - Dispersion-Correcting Potential Approach to the Accurate Treatment of both Inter- and Intramolecular Interactions

TL;DR

This paper enhances the B3LYP dispersion-correcting potentials by developing a new carbon DCP, significantly improving the accuracy of non-covalent interaction predictions and reaction barrier calculations in large molecular systems.

Contribution

The authors introduce a new carbon DCP for B3LYP that corrects previous errors and improves accuracy in modeling both inter- and intramolecular interactions.

Findings

New C-DCP improves accuracy of non-covalent interaction predictions.

B3LYP-DCP/6-31+G(2d,2p) achieves low MAEs on multiple databases.

Method accurately predicts reaction barriers within 1.7 kcal/mol of high-level calculations.

Abstract

We recently showed that dispersion-correcting potentials (DCPs), atom-centered Gaussian-type functions developed for use with B3LYP (J. Phys. Chem. Lett. 2012, 3, 1738-1744) greatly improved the ability of the underlying functional to predict non-covalent interactions. However, the application of B3LYP-DCP for the {\beta}-scission of the cumyloxyl radical led a calculated barrier height that was over-estimated by ca. 8 kcal/mol. We show in the present work that the source of this error arises from the previously developed carbon atom DCPs, which erroneously alters the electron density in the C-C covalent-bonding region. In this work, we present a new C-DCP with a form that was expected to influence the electron density farther from the nucleus. Tests of the new C-DCP, with previously published H-, N- and O-DCPs, with B3LYP-DCP/6-31+G(2d,2p) on the S66, S22B, HSG-A, and HC12 databases of non-covalently interacting dimers showed that it is one of the most accurate methods available for treating intermolecular interactions, giving mean absolute errors (MAEs) of 0.19, 0.27, 0.16, and 0.18 kcal/mol, respectively. Additional testing on the S12L database of complexation systems gave an MAE of 2.6 kcal/mol, showing that the B3LYP-DCP/6-31+G(2d,2p) approach is one of the best-performing and feasible methods for treating large systems dominated by non-covalent interactions. Finally, we showed that C-C making/breaking chemistry is well-predicted using the newly developed DCPs. In addition to predicting a barrier height for the {\beta}-scission of the cumyloxyl radical that is within 1.7 kcal/mol of the high-level value, application of B3LYP-DCP/6-31+G(2d,2p) to 10 databases that include reaction barrier heights and energies, isomerization energies and relative conformation energies gives performance that is amongst the best of all available dispersion-corrected density-functional theory approaches.