Density functional benchmark for quadruple hydrogen bonds

Hydrogen bonding is an important non-covalent interaction that plays a major role in molecular self-organization and supramolecular structures. It can be described accurately with ab initio quantum chemical wave function methods, which become computationally expensive for large molecular assemblies. Density functional theory (DFT) offers a better balance between accuracy and computational cost, and can be routinely applied to large systems. A large number of density functional approximations (DFAs) has been developed, but their accuracy depend on the application, necessitating benchmark studies to guide their selection for use in applications. Some of us have recently determined highly accurate hydrogen bonding energies of 14 quadruply hydrogen-bonded dimers by extrapolating coupled-cluster energies to the complete basis set limit as well as extrapolating electron correlation contributions with a continued-fraction approach [U. Ahmed et al., Phys. Chem. Chem. Phys., 2024, 26, 24470-24476]. In this work, we study the reproduction of these bonding energies at the DFT level using 152 DFAs. The top ten density functional approximations are composed of eight variants of the Berkeley functionals both with and without dispersion corrections, and two Minnesota 2011 functionals augmented with a further dispersion correction. We find the B97M-V functional with the non-local correlation functional replaced by an empirical D3BJ dispersion correction to be the best DFA, while changes to the dispersion part in other Berkeley functionals lead to poorer performance in our study.