Display options
Share it on

J Chem Theory Comput. 2011 Jan 11;7(1):88-96. doi: 10.1021/ct100469b. Epub 2010 Dec 02.

Assessment of the Performance of DFT and DFT-D Methods for Describing Distance Dependence of Hydrogen-Bonded Interactions.

Journal of chemical theory and computation

Kanchana S Thanthiriwatte, Edward G Hohenstein, Lori A Burns, C David Sherrill

Affiliations

  1. Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States.

PMID: 26606221 DOI: 10.1021/ct100469b

Abstract

Noncovalent interactions such as hydrogen bonds, van der Waals forces, and π-π interactions play important roles influencing the structure, stability, and dynamic properties of biomolecules including DNA and RNA base pairs. In an effort to better understand the fundamental physics of hydrogen bonding (H-bonding), we investigate the distance dependence of interaction energies in the prototype bimolecular complexes of formic acid, formamide, and formamidine. Potential energy curves along the H-bonding dissociation coordinate are examined both by establishing reference CCSD(T) interaction energies extrapolated to the complete basis set limit and by assessing the performance of the density functional methods B3LYP, PBE, PBE0, B970, PB86, M05-2X, and M06-2X and empirical dispersion corrected methods B3LYP-D3, PBE-D3, PBE0-D3, B970-D2, BP86-D3, and ωB97X-D, with basis sets 6-311++G(3df,3pd), aug-cc-pVDZ, and aug-cc-pVTZ. Although H-bonding interactions are dominated by electrostatics, it is necessary to properly account for dispersion interactions to obtain accurate energetics. In order to quantitatively probe the nature of hydrogen bonding interactions as a function of distance, we decompose the interaction energy curves into physically meaningful components with symmetry-adapted perturbation theory (SAPT). The SAPT results confirm that the contribution of dispersion and induction are significant at and near equilibrium, although electrostatics dominate. Among the DFT/DFT-D techniques, the best overall results are obtained utilizing counterpoise-corrected ωB97X-D with the aug-cc-pVDZ basis set.

Publication Types