Explicitly correlated double hybrid DFT: a comprehensive analysis of the basis set convergence on the GMTKN55 database

Abstract

Double-hybrid density functional theory (DHDFT) offers a pathway to accuracies approaching composite wavefunction approaches like G4 theory. However, the GLPT2 (Görling 2nd order perturbation theory) term causes them to partially inherit the slow ∝ L −3 (with L the maximum angular momentum) basis set convergence of correlated wavefunction methods. This could potentially be remedied by introducing F12 explicit correlation: we investigate the basis set convergence of both DHDFT and DHDFT-F12 for the large and chemically diverse GMTKN55 (general main-group thermochemistry, kinetics, and noncovalent interactions) benchmark suite. The B2GP-PLYP-D3(BJ) and revDSD-PBEP86-D4 double hybrid density functionals (DHDFs) are investigated as test cases, together with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets as large as cc-pV(Q+d)Z-F12. We show that F12 greatly accelerates basis set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12 basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPP in orbital-based approaches, and in fact