Partitioning dynamic electron correlation energy: Viewing Møller-Plesset correlation energies through Interacting Quantum Atom (IQA) energy partitioning

The nature of halogen bonding: insights from interacting quantum atoms and source function studies

A detailed study of the X...N (X = I, Br) halogen bonds in complexes formed by an extended set of substituted pyridines with D-X molecules (D = X, CN) is reported here. The nature of these interactions has been investigated at different (MP2 and DFT) levels of theory through Bader's quantum theory of atoms in molecules (QTAIM) and Pendás' interacting quantum atoms (IQA) scheme, focusing on the role of the local environment (i.e. the substituent on the pyridine ring and the halogenated residue) on the halogen bond features. We found that the exchange-correlation energy represents a substantial contribution to the IQA total energy, in some cases comparable to (I2 complexes) or even dominating (ICN complexes) the electrostatic term. Meaningful information is provided by the source function, indicating that the major contribution to the electron density at the bond critical point of the X...N interaction is derived from the halogen atom, while a much lower contribution comes from the nitrogen atom, which acts as either source or sink for electron density. A relevant contribution from distal atoms, including the various electron-donor and electron-withdrawing substituents in different positions of the pyridine ring, is also determined, highlighting the non-local character of the electron density. The existence of possible relationships between binding energies, interaction energies according to IQA, and QTAIM descriptors such as delocalization indices and source function, has been inspected. In general, good correlations are only found when the local environment, external to the directly involved halogen and nitrogen atoms, plays a minor role in the interaction.

Interacting Quantum Atoms and Multipolar Electrostatic Study of XH···π Interactions

The interaction energies of nine XH···π (X = C, N, and O) benzene-containing van der Waals complexes were analyzed, at the atomic and fragment levels, using QTAIM multipolar electrostatics and the energy partitioning method interacting quantum atoms/fragment (IQA/IQF). These descriptors were paired with the relative energy gradient method, which solidifies the connection between quantum mechanical properties and chemical interpretation. This combination provides a precise understanding, both qualitative and quantitative, of the nature of these interactions, which are ubiquitous in biochemical systems. The formation of the OH···π and NH···π systems is electrostatically driven, with the Qzz component of the quadrupole moment of the benzene carbons interacting with the charges of X and H in XH. There is the unexpectedly intramonomeric role of X-H (X = O, N) where its electrostatic energy helps the formation of the complex and its covalent energy thwarts it. However, the CH···π interaction is governed by exchange-correlation energies, thereby establishing a covalent character, as opposed to the literature's designation as a noncovalent interaction. Moreover, dispersion energy is relevant, statically and in absolute terms, but less relevant compared to other energy components in terms of the formation of the complex. Multipolar electrostatics are similar across all systems.

Non-covalent interactions from a Quantum Chemical Topology perspective

About half a century after its little-known beginnings, the quantum topological approach called QTAIM has grown into a widespread, but still not mainstream, methodology of interpretational quantum chemistry. Although often confused in textbooks with yet another population analysis, be it perhaps an elegant but somewhat esoteric one, QTAIM has been enriched with about a dozen other research areas sharing its main mathematical language, such as Interacting Quantum Atoms (IQA) or Electron Localisation Function (ELF), to form an overarching approach called Quantum Chemical Topology (QCT). Instead of reviewing the latter's role in understanding non-covalent interactions, we propose a number of ideas emerging from the full consequences of the space-filling nature of topological atoms, and discuss how they (will) impact on interatomic interactions, including non-covalent ones. The architecture of a force field called FFLUX, which is based on these ideas, is outlined. A new method called Relative Energy Gradient (REG) is put forward, which is able, by computation, to detect which fragments of a given molecular assembly govern the energetic behaviour of this whole assembly. This method can offer insight into the typical balance of competing atomic energies both in covalent and non-covalent case studies. A brief discussion on so-called bond critical points is given, highlighting concerns about their meaning, mainly in the arena of non-covalent interactions.

Implementation of the interacting quantum atom energy decomposition using the CASPT2 method

We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.

The nature of the intermolecular interaction in (H2X)2 (X = O, S, Se)

Hydrogen bonds (HBs) are crucial non-covalent interactions in chemistry. Recently, the occurrence of an HB in (H2S)2 has been reported (Arunan et al., Angew. Chem., Int. Ed., 2018, 57, 15199), challenging the textbook view of H2S dimers as mere van der Waals clusters. We herein try to shed light on the nature of the intermolecular interactions in the H2O, H2S, and H2Se dimers via correlated electronic structure calculations, Symmetry Adapted Perturbation Theory (SAPT) and Quantum Chemical Topology (QCT). Although (H2S)2 and (H2Se)2 meet some of the criteria for the occurrence of an HB, potential energy curves as well as SAPT and QCT analyses indicate that the nature of the interaction in (H2O)2 is substantially different (e.g. more anisotropic) from that in (H2S)2 and (H2Se)2. QCT reveals that the HB in (H2O)2 includes substantial covalent, dispersion and electrostatic contributions, while the last-mentioned component plays only a minor role in (H2S)2 and (H2Se)2. The major contributions to the interactions of the dimers of H2S and H2Se are covalency and dispersion as revealed by the exchange-correlation components of QCT energy partitions. The picture yielded by SAPT is somewhat different but compatible with that offered by QCT. Overall, our results indicate that neither (H2S)2 nor (H2Se)2 are hydrogen-bonded systems, showing how the nature of intermolecular contacts involving hydrogen atoms evolves in a group down the periodic table.

TWOE Code: An Efficient Tool for Explicit Partition of Coupled Cluster and Configuration Interaction Energies into Atomic and Diatomic Contributions

The efficient implementation of the TWOE program for evaluating the atomic and interatomic energy components at post-HF level was developed. The systematic convergence of these terms up to a near full-CI limit was performed for the first time for a series of coupled cluster methods: CCSD → CCSDT → CCSDTQ → CCSDTQP. A comparison with corresponding CI approaches (up to fifth excitation level) is additionally discussed. For a set of diatomic systems, it was demonstrated that, along with a full molecular energy convergence, all its components are also converged but with different patterns. It was found that not all components are decreased in their values at increasing computational rank. For instance, atomic energy parts are decreased while interatomic (interaction) energies are increased as the limiting level is approached. Two schemes were employed for atomic partition of molecules: the Baders approach and planes dissection. Influence of dynamical correlation effects on atomic energy components was analyzed in detail. Current TWOE implementation allows one, in principle, to work with any ab initio method providing the two-particle density matrix. It is believed that the developed program will be a useful tool for a real space energy decomposition that helps to reveal the most peculiar points in the structure of the total and correlation energies of a molecule.

Origin-Independent Decomposition of the Static Polarizability

Real-space analysis tools afford additive and transferable contributions of atoms to molecular properties. In the case of the molecular (hyper)polarizabilities, the atomic contributions that have been derived so far include a charge-transfer term that is origin-dependent. In this letter, we present the first genuinely origin-independent energy-based (OIEB) methodology for the decomposition of the static (hyper)polarizabilities that benefits from real-space molecular energy decomposition schemes, focusing on the static polarizability and showing that extension to static hyperpolarizabilities is straightforward. The numerical realization of the OIEB method shows the expected origin independence, atomic additivity, and transferability of atomic and functional group polarizability tensors. Furthermore, the OIEB atomic (fragment) polarizability tensors are symmetric by definition.

On the many-body nature of intramolecular forces in FFLUX and its implications

FFLUX is a biomolecular force field under construction, based on Quantum Chemical Topology (QCT) and machine learning (kriging), with a minimalistic and physically motivated design. A detailed analysis of the forces within the kriging models as treated in FFLUX is presented, taking as a test example a liquid water model. The energies of topological atoms are modeled as 3N_atoms -6 dimensional potential energy surfaces, using atomic local frames to represent the internal degrees of freedom. As a result, the forces within the kriging models in FFLUX are inherently N-body in nature where N refers to N_atoms . This provides a fuller picture that is closer to a true quantum mechanical representation of interactions between atoms. The presented computational example quantitatively showcases the non-negligible (as much as 9%) three-body nature of bonded forces and angular forces in a water molecule. We discuss the practical impact on the pressure calculation with N-body forces and periodic boundary conditions (PBC) in molecular dynamics, as opposed to classical force fields with two-body forces. The equivalence between the PBC-related correction terms in the general virial equation is shown mathematically.

Interacting Quantum Atoms-A Review

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study of the Halogen Bond with Explicit Analysis of Electron Correlation

Energy profiles of seven halogen-bonded complexes were analysed with the topological energy partitioning called Interacting Quantum Atoms (IQA) at MP4(SDQ)/6-31+G(2d,2p) level of theory. Explicit interatomic electron correlation energies are included in the analysis. Four complexes combine X₂ (X = Cl or F) with HCN or NH₃, while the remaining three combine ClF with HCN, NH₃ or N₂. Each complex was systematically deformed by translating the constituent molecules along its central axis linking X and N, and reoptimising its remaining geometry. The Relative Energy Gradient (REG) method (Theor. Chem. Acc. 2017, 136, 86) then computes which IQA energies most correlate with the total energy during the process of complex formation and further compression beyond the respective equilibrium geometries. It turns out that the covalent energy (i.e., exchange) of the halogen bond, X…N, itself drives the complex formation. When the complexes are compressed from their equilibrium to shorter X…N distance then the intra-atomic energy of N is in charge. When the REG analysis is restricted to electron correlation then the interatomic correlation energy between X and N again drives the complex formation, and the complex compression is best described by the destabilisation of the through-space correlation energy between N and the "outer" halogen.

Contributions of IQA electron correlation in understanding the chemical bond and non-covalent interactions

The quantum topological energy partitioning method Interacting Quantum Atoms (IQA) has been applied for over a decade resulting in an enlightening analysis of a variety of systems. In the last three years we have enriched this analysis by incorporating into IQA the two-particle density matrix obtained from Møller–Plesset (MP) perturbation theory. This work led to a new computational and interpretational tool to generate atomistic electron correlation and thus topologically based dispersion energies. Such an analysis determines the effects of electron correlation within atoms and between atoms, which covers both bonded and non-bonded “through -space” atom–atom interactions within a molecule or molecular complex. A series of papers published by us and other groups shows that the behavior of electron correlation is deeply ingrained in structural chemistry. Some concepts that were shown to be connected to bond correlation are bond order, multiplicity, aromaticity, and hydrogen bonding. Moreover, the concepts of covalency and ionicity were shown not to be mutually excluding but to both contribute to the stability of polar bonds. The correlation energy is considerably easier to predict by machine learning (kriging) than other IQA terms. Regarding the nature of the hydrogen bond, correlation energy presents itself in an almost contradicting way: there is much localized correlation energy in a hydrogen bond system, but its overall effect is null due to internal cancelation. Furthermore, the QTAIM delocalization index has a connection with correlation energy. We also explore the role of electron correlation in protobranching, which provides an explanation for the extra stabilization present in branched alkanes compared to their linear counterparts. We hope to show the importance of understanding the true nature of the correlation energy as the foundation of a modern representation of dispersion forces for ab initio, DFT, and force field calculations.

Atomic Partitioning of the MPn (n = 2, 3, 4) Dynamic Electron Correlation Energy by the Interacting Quantum Atoms Method: A Fast and Accurate Electrostatic Potential Integral Approach

Recently, the quantum topological energy partitioning method called interacting quantum atoms (IQA) has been extended to MPn (n = 2, 3, 4) wave functions. This enables the extraction of chemical insight related to dynamic electron correlation. The large computational expense of the IQA-MPn approach is compensated by the advantages that IQA offers compared to older nontopological energy decomposition schemes. This expense is problematic in the construction of a machine learning training set to create kriging models for topological atoms. However, the algorithm presented here markedly accelerates the calculation of atomically partitioned electron correlation energies. Then again, the algorithm cannot calculate pairwise interatomic energies because it applies analytical integrals over whole space (rather than over atomic volumes). However, these pairwise energies are not needed in the quantum topological force field FFLUX, which only uses the energy of an atom interacting with all remaining atoms of the system that it is part of. Thus, it is now feasible to generate accurate and sizeable training sets at MPn level of theory. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Using the Relative Energy Gradient Method with Interacting Quantum Atoms to Determine the Reaction Mechanism and Catalytic Effects in the Peptide Hydrolysis in HIV-1 Protease

The reaction mechanism in an active site is of the utmost importance when trying to understand the role that an enzyme plays in biological processes. In a recently published paper [Theor. Chem. Acc. 2017, 136, 86], we formalised the Relative Energy Gradient (REG) method for automating an Interacting Quantum Atoms (IQA) analysis. Here, the REG method is utilised to determine the mechanism of peptide hydrolysis in the aspartic active site of the enzyme HIV-1 Protease. Using the REG method along with the IQA approach we determine the mechanism of peptide hydrolysis without employing any arbitrary parameters and with remarkable ease (albeit at large computational cost: the system contains 133 atoms, which means that there are 17 689 individual IQA terms to be calculated). When REG and IQA work together it is possible to determine a reaction mechanism at atomistic resolution from data directly derived from quantum calculations, without arbitrary parameters. Moreover, the mechanism determined by this novel method gives concrete insight into how the active site residues catalyse peptide hydrolysis.

MP2-IQA: upscaling the analysis of topologically partitioned electron correlation

When electronic correlation energy is partitioned topologically, a detailed picture of its distribution emerges, both within atoms and between any two atoms. This methodology allows one to study dispersion beyond its more narrow definition in long-range Rayleigh-Schrödinger perturbation theory. The interacting quantum atoms (IQA) method was applied to MP2/6-31G(d,p) (uncontracted) wave functions of a wide variety of systems: glycine…water (hydration), the ethene dimer (π-π interactions), benzene (aromaticity), cyclobutadiene (antiaromaticity), and NH₃BH₃ (dative bond). Through the study of molecular complexes it turns out that dispersion energy is either important to a system's stabilization (for the C₂H₄ dimer) or not important (for Gly…H₂O). We have also discovered that the delocalization in benzene lowers the strength of Coulomb repulsion in the bonds, which has been quantified for the first time through IQA. Finally, we showed that the nature of the dative bond is much different from that of a regular covalent bond as it is not destabilized by electronic correlation. Finally, the conclusions obtained for these archetypical systems have implications for the future of the quantum topological force field FFLUX in the simulation of larger systems. Graphical abstract Atomic and bond dynamic correlation energies are now available thanks to IQA. Larger molecules can now be accessed to include resonance and solvation of FFLUX force field.

Fluorine Gauche Effect Explained by Electrostatic Polarization Instead of Hyperconjugation: An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study

We present an interacting quantum atoms (IQA) study of the gauche effect by comparing 1,2-difluoroethane, 1,2-dichloroethane, and three conformers of 1,2,3,4,5,6-hexafluorocyclohexane. In the 1,2-difluoroethane, the gauche effect is observed in that the gauche conformation is more stable than the anti, whereas in 1,2-dichloroethane the opposite is true. The analysis performed here is exhaustive and unbiased thanks to using the recently introduced relative energy gradient (REG) method [ Thacker , J. C. R. ; Popelier , P. L. A. Theor. Chem. Acc . 2017 , 136 , 86 ], as implemented in the in-house program ANANKE. We challenge the common explanation that hyperconjugation is responsible for the gauche stability in 1,2-difluoroethane and instead present electrostatics as the cause of gauche stability. Our explanation of the gauche effect is also is seen in other molecules displaying local gauche conformations, such as the recently synthesized "all-cis" hexafluorocyclohexane and its conformers where all the fluorine atoms are in the equatorial positions. Using our extension of the traditional IQA methodology that allows for the partitioning of electrostatic terms into polarization and charge transfer, we propose that the cause of gauche stability is 1,3 C···F electrostatic polarization interactions. In other words, if a number of fluorine atoms are aligned, then the stability due to polarization of nearby carbon atoms is increased.

Machine Learning of Dynamic Electron Correlation Energies from Topological Atoms

We present an innovative method for predicting the dynamic electron correlation energy of an atom or a bond in a molecule utilizing topological atoms. Our approach uses the machine learning method Kriging (Gaussian Process Regression with a non-zero mean function) to predict these dynamic electron correlation energy contributions. The true energy values are calculated by partitioning the MP2 two-particle density-matrix via the Interacting Quantum Atoms (IQA) procedure. To our knowledge, this is the first time such energies have been predicted by a machine learning technique. We present here three important proof-of-concept cases: the water monomer, the water dimer, and the van der Waals complex H₂···He. These cases represent the final step toward the design of a full IQA potential for molecular simulation. This final piece will enable us to consider situations in which dispersion is the dominant intermolecular interaction. The results from these examples suggest a new method by which dispersion potentials for molecular simulation can be generated.

Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective

Accurate description of the intrinsic preferences of amino acids is important to consider when developing a biomolecular force field. In this study, we use a modern energy partitioning approach called Interacting Quantum Atoms to inspect the cause of the φ and ψ torsional preferences of three dipeptides (Gly, Val, and Ile). Repeating energy trends at each of the molecular, functional group, and atomic levels are observed across both (1) the three amino acids and (2) the φ/ψ scans in Ramachandran plots. At the molecular level, it is surprisingly electrostatic destabilization that causes the high-energy regions in the Ramachandran plot, not molecular steric hindrance (related to the intra-atomic energy). At the functional group and atomic levels, the importance of key peptide atoms (Oi-1 , Ci , Ni , Ni+1 ) and some sidechain hydrogen atoms (Hγ ) are identified as responsible for the destabilization seen in the energetically disfavored Ramachandran regions. Consistently, the Oi-1 atoms are particularly important for the explanation of dipeptide intrinsic behavior, where electrostatic and steric destabilization unusually complement one another. The findings suggest that, at least for these dipeptides, it is the peptide group atoms that dominate the intrinsic behavior, more so than the sidechain atoms. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

The Transferability of Topologically Partitioned Electron Correlation Energies in Water Clusters

The electronic effects that govern the cohesion of water clusters are complex, demanding the inclusion of N-body, Coulomb, exchange and correlation effects. Here we present a much needed quantitative study of the effect of correlation (and hence dispersion) energy on the stabilization of water clusters. For this purpose we used a topological energy partitioning method called Interacting Quantum Atoms (IQA) to partition water clusters into topological atoms, based on a MP2/6-31G(d,p) wave function, and modified versions of GAUSSIAN09 and the Quantum Chemical Topology (QCT) program MORFI. Most of the cohesion in the water clusters provided by electron correlation comes from intramolecular energy stabilization. Hydrogen bond-related interactions tend to largely cancel each other. Electron correlation energies are transferable in almost all instances within 1 kcal mol^-1 . This observed transferability is very important to the further development of the QCT force field FFLUX, especially to the future modelling of liquid water.

Quantifying Electron Correlation of the Chemical Bond

The Interacting Quantum Atoms (IQA) method is used to analyze the correlated part of the Møller-Plesset (MP) perturbation theory two-particle density matrix. Such an analysis determines the effects of electron correlation within atoms and between atoms, which covers both bonds and nonbonded through-space atom-atom interactions within a molecule or molecular complex. Electron correlation lowers the energy of the atoms at either end of a bond, but for the bond itself, it can be stabilizing or destabilizing. Bonds are described in a two-dimensional world of exchange and charge transfer, where covalency is not the opposite of ionicity.