Journal
JOURNAL OF COMPUTATIONAL CHEMISTRY
Volume 40, Issue 32, Pages 2868-2881Publisher
WILEY
DOI: 10.1002/jcc.26068
Keywords
hydrogen bond; binding energy; supermolecular approach; symmetry-adapted perturbation theory; correlation analysis; quantum theory of atoms-in-molecules; coupled cluster
Categories
Ask authors/readers for more resources
This work studies the underlying nature of H-bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry-adapted perturbation theory (SAPT). Quantum theory of atoms-in-molecules (QTAIM) descriptors and other HB indices were calculated based on high-quality density functional theory wave functions. We propose a new and rigorous classification of H-bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as very weak HBs with a BE >= -2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as weak-to-medium HBs with a BE varying between -2.5 and -14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into medium HBs with a BE in the range of -11.0 to -15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into strong HBs whose BE is more negative than -15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function-based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; rho(BCP)) of HBs. The fitted equation for neutral complexes is BE/kcal/mol = - 223.08 x rho(BCP)/a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE/kcal/mol = - 332.34 x rho(BCP)/a. u. - 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. (c) 2019 Wiley Periodicals, Inc.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available