Performance of DFT hybrid functionals in the theoretical treatment of H-bonds: Analysis term-by-term

The performance of DFT to reproduce noncovalent interactions like H-bonds was assessed. Three DFT hybrid functionals (B3LYP, BHandHLYP from Jaguar5.5, and BHandHLYP from G03) were used to calculate the interaction energies and H-bond distances of several host-guest complexes theoretically designed. Two reference systems (whose experimental data were available) were also calculated at the same levels of theory. In all the cases B3LYP and BHandHLYP from G03 gave rather the same results of interaction energy and distances, whereas the functional BHandHLYP from Jaguar5.5 overestimated the interaction energies and underestimated the H-bond distances. Considering the construction of each functional, the terms responsible for such differences are the gradient correction to the exchange functional (Becke88) and the correlation functional (VWN, LYP) and not the Hartree-Fock contribution.