First-Principles Modeling of Non-Covalent Interactions in Supramolecular Systems: The Role of Many-Body Effects

Journal Article · Thu Oct 18 04:00:00 EDT 2012 · Journal of Chemical Theory and Computation

DOI:https://doi.org/10.1021/ct300711r· OSTI ID:1564977

Tkatchenko, Alexandre ^[1]; Alfè, Dario ^[2]; Kim, Kwang S. ^[3]

---

1. Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany; Department of Chemistry, Pohang University of Science and Technology, Pohang 790−784, South Korea
2. London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom
3. Department of Chemistry, Pohang University of Science and Technology, Pohang 790−784, South Korea

Supramolecular host–guest systems play an important role for a wide range of applications in chemistry and biology. The prediction of the stability of host–guest complexes represents a great challenge to first-principles calculations due to an interplay of a wide variety of covalent and noncovalent interactions in these systems. In particular, van der Waals (vdW) dispersion interactions frequently play a prominent role in determining the structure, stability, and function of supramolecular systems. On the basis of the widely used benchmark case of the buckyball catcher complex (C60@C60H28), we assess the feasibility of computing the binding energy of supramolecular host–guest complexes from first principles. Large-scale diffusion Monte Carlo (DMC) calculations are carried out to accurately determine the binding energy for the C60@C60H28 complex (26 ± 2 kcal/mol). On the basis of the DMC reference, we assess the accuracy of widely used and efficient density-functional theory (DFT) methods with dispersion interactions. The inclusion of vdW dispersion interactions in DFT leads to a large stabilization of the C60@C60H28 complex. However, DFT methods including pairwise vdW interactions overestimate the stability of this complex by 9–17 kcal/mol compared to the DMC reference and the extrapolated experimental data. A significant part of this overestimation (9 kcal/mol) stems from the lack of dynamical dielectric screening effects in the description of the molecular polarizability in pairwise dispersion energy approaches. The remaining overstabilization arises from the isotropic treatment of atomic polarizability tensors and the lack of many-body dispersion interactions. A further assessment of a different supramolecular system – glycine anhydride interacting with an amide macrocycle – demonstrates that both the dynamical screening and the many-body dispersion energy are complex contributions that are very sensitive to the underlying molecular geometry and type of bonding. We discuss the required improvements in theoretical methods for achieving “chemical accuracy” in the first-principles modeling of supramolecular systems.

• Research Organization: Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); UT-Battelle LLC/ORNL, Oak Ridge, TN (Unted States)
• Sponsoring Organization: USDOE Office of Science
• DOE Contract Number: AC05-00OR22725
• OSTI ID: 1564977
• Journal Information: Journal of Chemical Theory and Computation, Journal Name: Journal of Chemical Theory and Computation Journal Issue: 11 Vol. 8; ISSN 1549-9618
• Publisher: American Chemical Society
• Country of Publication: United States
• Language: English

Related Subjects

Chemistry
Physics