Molecular structure does not easily identify the intricate non-covalent interactions that govern many areas of biology and chemistry, including design of new materials and drugs. We develop an approach to detect non-covalent interactions in real space, based on the electron density and its derivatives. Our approach reveals underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of non-bonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.
Non-covalent interactions hold the key to understanding many chemical, biological, and technological problems. Describing these non-covalent interactions accurately, including their positions in real space, constitutes a first step in the process of decoupling the complex balance of forces that define non-covalent interactions. Because of the size of macromolecules, the most common approach has been to assign van der Waals interactions (vdW), steric clashes (SC), and hydrogen bonds (HBs) based on pairwise distances between atoms according to their van der Waals radii. We recently developed an alternative perspective, derived from the electronic density: the Non-Covalent Interactions (NCI) index [J. Am. Chem. Soc. 2010, 132, 6498]. This index has the dual advantages of being generally transferable to diverse chemical applications and being very fast to compute, since it can be calculated from promolecular densities. Thus, NCI analysis is applicable to large systems, including proteins and DNA, where analysis of non-covalent interactions is of great potential value. Here, we describe the NCI computational algorithms and their implementation for the analysis and visualization of weak interactions, using both selfconsistent fully quantum-mechanical, as well as promolecular, densities. A wide range of options for tuning the range of interactions to be plotted is also presented. To demonstrate the capabilities of our approach, several examples are given from organic, inorganic, solid state, and macromolecular chemistry, including cases where NCI analysis gives insight into unconventional chemical bonding. The NCI code and its manual are available for download at
A procedure based on density functional theory is used for the calculation of the gas-phase bond dissociation enthalpy (BDE) and ionization potential for molecules belonging to the class of phenolic antioxidants. We show that use of locally dense basis sets (LDBS) vs full basis sets gives very similar results for monosubstituted phenols, and that the LDBS procedure gives good agreement with the change in experimental BDE values for highly substituted phenols in benzene solvent. Procedures for estimating the O--H BDE based on group additivity rules are given and tested. Several interesting classes of phenolic antioxidants are studied with these methods, including commercial antioxidants used as food additives, compounds related to Vitamin E, flavonoids in tea, aminophenols, stilbenes related to resveratrol, and sterically hindered phenols. On the basis of these results we are able to interpret relative rates for the reaction of antioxidants with free radicals, including a comparison of both H-atom-transfer and single-electron-transfer mechanisms, and conclude that in most cases H-atom transfer will be dominant.
We have recently introduced [J. Chem. Phys. 122, 154104 (2005)] a simple parameter-free model of the dispersion interaction based on the instantaneous in space, dipole moment of the exchange hole. The model generates remarkably accurate interatomic and intermolecular C6 dispersion coefficients, and geometries and binding energies of intermolecular complexes. The model involves, in its original form, occupied Hartree-Fock or Kohn-Sham orbitals. Here we present a density-functional reformulation depending only on total density, the gradient and Laplacian of the density, and the kinetic-energy density. This density-functional model performs as well as the explicitly orbital-dependent model, yet offers obvious computational advantages.
We have previously demonstrated that the dipole moment of the exchange hole can be used to derive intermolecular C(6) dispersion coefficients [J. Chem. Phys. 122, 154104 (2005)]. This was subsequently the basis for a novel post-Hartree-Fock model of intermolecular interactions [J. Chem. Phys. 123, 024101 (2005)]. In the present work, the model is extended to include higher-order dispersion coefficients C(8) and C(10). The extended model performs very well for prediction of intermonomer separations and binding energies of 45 van der Waals complexes. In particular, it performs twice as well as basis-set extrapolated MP2 theory for dispersion-bound complexes, with minimal computational cost.
The optimized effective potential (OEP) for exchange was introduced some time ago by Sharp and Horton and by Talman and Shadwick. The integral equation for the OEP is difficult to solve, however, and a variety of approximations have therefore been proposed. These are explicitly orbital dependent and require the same two-electron integrals as Hartree-Fock theory. We have found a remarkably simple approximate effective potential that closely resembles the Talman-Shadwick potential in atoms. It depends only on total densities and requires no two-electron integrals.
Intermolecular interactions are of great importance in chemistry but are difficult to model accurately with computational methods. In particular, Hartree-Fock and standard density-functional approximations do not include the physics necessary to properly describe dispersion. These methods are sometimes corrected to account for dispersion by adding a pairwise C6R6 term, with C6 dispersion coefficients dependent on the atoms involved. We present a post-Hartree-Fock model in which C6 coefficients are generated by the instantaneous dipole moment of the exchange hole. This model relies on occupied orbitals only, and involves only one, universal, empirical parameter to limit the dispersion energy at small interatomic separations. The model is extensively tested on isotropic C6 coefficients of 178 intermolecular pairs. It is also applied to the calculation of the geometries and binding energies of 20 intermolecular complexes involving dispersion, dipole-induced dipole, dipole-dipole, and hydrogen-bonding interactions, with remarkably good results.
A benchmark for non-covalent interactions in solids (C21) based on the experimental sublimation enthalpies and geometries of 21 molecular crystals is presented. Thermal and zero-point effects are carefully accounted for and reference lattice energies and thermal pressures are provided, which allow dispersion-corrected density functionals to be assessed in a straightforward way. Other thermal corrections to the sublimation enthalpy (the 2RT term) are reexamined. We compare the recently implemented exchange-hole dipole moment (XDM) model with other approaches in the literature to find that XDM roughly doubles the accuracy of DFT-D2 and non-local functionals in computed lattice energies (4.8 kJ/mol mean absolute error) while, at the same time, predicting cell geometries within less than 2% of the experimental result on average. The XDM model of dispersion interactions is confirmed as a very promising approach in solid-state applications.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.