Estimations of Fe–N2 Intrinsic Interaction Energies of Iron–Sulfur/Nitrogen–Carbon Sites: A Deeper Bonding Insight by EDA-NOCV Analysis of a Model Complex of the Nitrogenase Cofactor

The intrinsic (gas-phase) stacking energies of natural and artificial nucleobases were explored using density functional theory (DFT) and correlated ab initio methods. Ranking the stacking strength of natural nucleobase dimers revealed a preference in binding partner similar to that seen from experiments, namely G > C > A > T > U. Decomposition of these interaction energies using symmetry-adapted perturbation theory (SAPT) showed that these dispersion dominated interactions are modulated by electrostatics. Artificial nucleobases showed a similar stacking preference for natural nucleobases and were also modulated by electrostatic interactions. A robust predictive multivariate model was developed that quantitively predicts the maximum stacking interaction between natural and a wide range of artificial nucleobases using molecular descriptors based on computed electrostatic potentials (ESPs) and the number of heavy atoms. This model should find utility in designing artificial nucleobase analogs that exhibit stacking interactions comparable to those of natural nucleobases. Further analysis of the descriptors in this model unveil the origin of superior stacking abilities of certain nucleobases, including cytosine and guanine.

Download Full-text

On the Nature of Halide-Water Interactions: Insights from Many-Body Representations and Density Functional Theory

10.26434/chemrxiv.7613042.v2 ◽

2019 ◽

Author(s):

Brandon B. Bizzarro ◽

Colin K. Egan ◽

Francesco Paesani

Keyword(s):

Density Functional Theory ◽

Charge Transfer ◽

Molecular Orbital ◽

Density Functional ◽

Decomposition Analysis ◽

Orbital Energy ◽

Energy Decomposition Analysis ◽

Interaction Energies ◽

Many Body ◽

Functional Theory

<div> <div> <div> <p>Interaction energies of halide-water dimers, X<sup>-</sup>(H<sub>2</sub>O), and trimers, X<sup>-</sup>(H<sub>2</sub>O)<sub>2</sub>, with X = F, Cl, Br, and I, are investigated using various many-body models and exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. Analysis of the results obtained with the many-body models demonstrates the need to capture important short-range interactions in the regime of large inter-molecular orbital overlap, such as charge transfer and charge penetration. Failure to reproduce these effects can lead to large deviations relative to reference data calculated at the coupled cluster level of theory. Decompositions of interaction energies carried out with the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA) method demonstrate that permanent and inductive electrostatic energies are accurately reproduced by all classes of XC functionals (from generalized gradient corrected (GGA) to hybrid and range-separated functionals), while significant variance is found for charge transfer energies predicted by different XC functionals. Since GGA and hybrid XC functionals predict the most and least attractive charge transfer energies, respectively, the large variance is likely due to the delocalization error. In this scenario, the hybrid XC functionals are then expected to provide the most accurate charge transfer energies. The sum of Pauli repulsion and dispersion energies are the most varied among the XC functionals, but it is found that a correspondence between the interaction energy and the ALMO EDA total frozen energy may be used to determine accurate estimates for these contributions. </p> </div> </div> </div>

Download Full-text