An official website of the United States government
In this work, we present NENCI-2021, a benchmark database of ∼8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (Eint) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The Eint values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from −38.5 to +186.8 kcal/mol, with a mean (median) Eint value of −1.06 kcal/mol (−2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the Eint values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive Eint values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ∼±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.
more »
« less
The influence of a solvent environment on direct non-covalent interactions between two molecules: A symmetry-adapted perturbation theory study of polarization tuning of π – π interactions by water
High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π– π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller–Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π– π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol −1 , while significant changes of up to 3 kcal mol −1 of the electrostatic component are seen for the cationic pyridinium–benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little “tuning” of the direct solute–solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute–solute and solute–solvent interactions.
more »
« less
- Award ID(s):
- 1955940
- PAR ID:
- 10353744
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 19
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 194306
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Both methyl groups and benzene rings are exceedingly common, and they lie near one another in many chemical situations. DFT calculations are used to gauge the strength of the attractive forces between them, and to better understand the phenomena that underlie this attraction. Methane and benzene are taken as the starting point, and substituents of both electron‐withdrawing and donating types are added to each. The interaction energy varies between 1.4 and 5.0 kcal/mol, depending upon the substituents placed on the two groups. The nature of the binding is analyzed via Atoms in Molecules (AIM), Natural Bond Orbital (NBO), Symmetry‐Adapted Perturbation Theory (SAPT), nuclear magnetic resonance (NMR) chemical shifts, and electron density shift diagrams. While there is a sizable electrostatic component, it is dispersion that dominates these interactions, particularly the weaker ones. As such, these interactions cannot be categorized unambiguously as either H‐bonds or tetrel bonds.more » « less
-
The need for a deep understanding of CO2 interactions with other mols. is significant given the importance of supercrit. CO2 (s.c.-CO2) as a green solvent, and interest in design of novel materials for CO2 capture and storage. Soly. of fluorinated compds. in s.c.-CO2 is generally higher than their hydrocarbon analogs and fluorination of hydrocarbon compds. improves "CO2-philicity". Although dissoln. of a compd. in s.c.-CO2 (or any solvent) is a complex process, much can be learned by systematic examn. of solute-solvent interactions as a function of solute mol. properties in dimers and other small clusters. In the present study, Fourier-transform microwave spectroscopy and Symmetry Adapted Perturbation Theory (SAPT) calcns. have been used to examine intermol. interactions between CO2 and a series of fluorinated ethylene mols. with varying degree and position of fluorination (1-3 F atoms). While 1-fluoro-, 1,1-difluoro- and 1,1,2-trifluoroethylene...CO2 complexes are planar (with two isomers obsd. for 1-fluoroethylene), cis-1,2-difluoroethylene...CO2 is nonplanar. Although lowest energy structures predicted by MP2/6-311++G(2d,2p) calcns. are not always in agreement with obsd. configurations, SAPT calcns. provide binding energies consistent with the obsd. structures. We are now extending these studies to trimers contg. one fluorinated ethylene with two CO2 "solvent" mols. in order to gain further insight into structural and energetic changes as a solvation shell is formed.more » « less
-
A quantum chemical study of the interactions between N‐methylacetamide (NMA), which serves as a peptide model, and poloxamer building blocks made up of polyethylene oxide (PEO) and polypropylene oxide (PPO) dimers shows a range of noncovalent interactions. Eight stable NMA···PEO complexes were found, with interaction energies between 2.24 and 6.73 kcal/mol. NMA···PPO complexes show slightly stronger interactions, reaching up to 7.39 kcal/mol. Besides NH···O and CH···O hydrogen bonding, several complexes also exhibit notable nonclassical C═O···O interactions. These are stabilized by bidirectional charge transfer (CT) that includes C═O π → O–C σ* and O(lp) → C═O π* donations. Such interactions improve binding strengths to levels similar to those of typical hydrogen bonds. The findings highlight that hydrogen bonding and CT work together in peptide–poloxamer interactions.more » « less
-
Hypohalous acids (HOX) are a class of molecules that play a key role in the atmospheric seasonal depletion of ozone and have the ability to form both hydrogen and halogen bonds. The interactions between the HOX monomers (X = F, Cl, Br) and water have been studied at the CCSD(T)/aug-cc-pVTZ level of theory with the spin free X2C-1e method to account for scalar relativistic effects. Focal point analysis was used to determine CCSDT(Q)/CBS dissociation energies. The anti hydrogen bonded dimers were found with interaction energies of −5.62 kcal mol −1 , −5.56 kcal mol −1 , and −4.97 kcal mol −1 for X = F, Cl, and Br, respectively. The weaker halogen bonded dimers were found to have interaction energies of −1.71 kcal mol −1 and −3.03 kcal mol −1 for X = Cl and Br, respectively. Natural bond orbital analysis and symmetry adapted perturbation theory were used to discern the nature of the halogen and hydrogen bonds and trends due to halogen substitution. The halogen bonds were determined to be weaker than the analogous hydrogen bonds in all cases but close enough in energy to be relevant, significantly more so with increasing halogen size.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1063/5.0087302
