Electron correlation in the interacting quantum atoms partition via coupled-cluster lagrangian densities

Exact decompositions of the total KS-DFT exchange-correlation energy into one- and two-center terms

In the so-called Interacting Quantum Atoms (IQA) approach, the molecular energy is numerically decomposed as a sum of atomic and diatomic contributions. While proper formulations have been put forward for both Hartree-Fock and post-Hartree-Fock wavefunctions, this is not the case for the Kohn-Sham density functional theory (KS-DFT). In this work, we critically analyze the performance of two fully additive approaches for the IQA decomposition of the KS-DFT energy, namely, the one from Francisco et al., which uses atomic scaling factors, and that from Salvador and Mayer based upon the bond order density (SM-IQA). Atomic and diatomic exchange-correlation (xc) energy components are obtained for a molecular test set comprising different bond types and multiplicities and along the reaction coordinate of a Diels-Alder reaction. Both methodologies behave similarly for all systems considered. In general, the SM-IQA diatomic xc components are less negative than the Hartree-Fock ones, which is in good agreement with the known effect of electron correlation upon (most) covalent bonds. In addition, a new general scheme to minimize the numerical error of the sum of two-electron energy contributions (i.e., Coulomb and exact exchange) in the framework of overlapping atoms is described in detail.

The IQA Energy Partition in a Drug Design Setting: A Hepatitis C Virus RNA-Dependent RNA Polymerase (NS5B) Case Study

The interaction of the thumb site II of the NS5B protein of hepatitis C virus and a pair of drug candidates was studied using a topological energy decomposition method called interacting quantum atoms (IQA). The atomic energies were then processed by the relative energy gradient (REG) method, which extracts chemical insight by computation based on minimal assumptions. REG reveals the most important IQA energy contributions, by atom and energy type (electrostatics, sterics, and exchange–correlation), that are responsible for the behaviour of the whole system, systematically from a short-range ligand–pocket interaction until a distance of approximately 22 Å. The degree of covalency in various key interatomic interactions can be quantified. No exchange–correlation contribution is responsible for the changes in the energy profile of both pocket–ligand systems investigated in the ligand–pocket distances equal to or greater than that of the global minimum. Regarding the hydrogen bonds in the system, a “neighbour effect” was observed thanks to the REG method, which states that a carbon atom would rather not have its covalent neighbour oxygen form a hydrogen bond. The combination of IQA and REG enables the automatic identification of the pharmacophore in the ligands. The coarser Interacting Quantum Fragments (IQF) enables the determination of which amino acids of the pocket contribute most to the binding and the type of energy of said binding. This work is an example of the contribution topological energy decomposition methods can make to fragment-based drug design.

Binding Energy Partition of Promising IRAK-4 Inhibitor (Zimlovisertib) for the Treatment of COVID-19 Pneumonia

The technique of Fragment-Based Drug Design (FBDD) considers the interactions of different moieties of molecules with biological targets for the rational construction of potential drugs. One basic assumption of FBDD is that the different functional groups of a ligand interact with a biological target in an approximately additive, that is, independent manner. We investigated the interactions of different fragments of ligands and Interleukin-1 Receptor-Associated Kinase 4 (IRAK-4) throughout the FBDD design of Zimlovisertib, a promising anti-inflammatory, currently in trials to be used for the treatment of COVID-19 pneumonia. We utilised state-of-the-art methods of wave function analyses mainly the Interacting Quantum Atoms (IQA) energy partition for this purpose. By means of IQA, we assessed the suitability of every change to the ligand in the five stages of FBDD which led to Zimlovisertib on a quantitative basis. We determined the energetics of the interaction of different functional groups in the ligands with the IRAK-4 protein target and thereby demonstrated the adequacy (or lack thereof) of the changes made across the design of this drug. This analysis permits to verify whether a given alteration of a prospective drug leads to the intended tuning of non-covalent interactions with its protein objective. Overall, we expect that the methods exploited in this paper will prove valuable in the understanding and control of chemical modifications across FBDD processes.

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.

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.

Photochemistry in Real Space: Batho- and Hypsochromism in the Water Dimer

The development of chemical intuition in photochemistry faces several difficulties that result from the inadequacy of the one-particle picture, the Born-Oppenheimer approximation, and other basic ideas used to build models. It is shown herein how real-space approaches can be efficiently used to gain valuable insights in photochemistry through a simple example of red and blue shift effects: the double hypso- and bathochromic shifts in the low-lying valence excited states of (H₂ O)₂ . It is demonstrated that 1) the use of these techniques allows the perturbative language used in the theory of intermolecular interactions, even in the strongly interacting short-range regime, to be maintained; 2) one and only one molecule is photoexcited in each of the addressed excited states and 3) the electrostatic interaction between the in-the-cluster molecular dipoles provides a fairly intuitive rationalisation of the observed batho- and hypsochromism. The methods exploited and illustrated herein are able to maintain the individuality and properties of the interacting entities in a molecular aggregate, and thereby they allow chemical intuition in general states, at any geometry and using a broad variety of electronic structure methods to be kept and built.

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.

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.

Tetrel Interactions from an Interacting Quantum Atoms Perspective

Tetrel bonds, the purportedly non-covalent interaction between a molecule that contains an atom of group 14 and an anion or (more generally) an atom or molecule with lone electron pairs, are under intense scrutiny. In this work, we perform an interacting quantum atoms (IQA) analysis of several simple complexes formed between an electrophilic fragment (A) (CH₃F, CH₄, CO₂, CS₂, SiO₂, SiH₃F, SiH₄, GeH₃F, GeO₂, and GeH₄) and an electron-pair-rich system (B) (NCH, NCO-, OCN-, F-, Br-, CN-, CO, CS, Kr, NC-, NH₃, OC, OH₂, SH-, and N₃-) at the aug-cc-pvtz coupled cluster singles and doubles (CCSD) level of calculation. The binding energy ( E bind AB ) is separated into intrafragment and inter-fragment components, and the latter in turn split into classical and covalent contributions. It is shown that the three terms are important in determining E bind AB , with absolute values that increase in passing from electrophilic fragments containing C, Ge, and Si. The degree of covalency between A and B is measured through the real space bond order known as the delocalization index ( δ AB ). Finally, a good linear correlation is found between δ AB and E xc AB , the exchange correlation (xc) or covalent contribution to E bind AB .

Reliability of HF/IQA, B3LYP/IQA, and MP2/IQA data in interpreting the nature and strength of interactions

Using CCSD/BBC1/IQA as a reference to describe interactions in glycol conformers, including O4⋯H6 and O3⋯O4 interactions, B3LYP/IQA data produced an exact qualitative picture and exceptionally reliable IQA energy changes, resulting in a perfectly comparable FAMSEC-based quantitative interpretation.

Partition of electronic excitation energies: the IQA/EOM-CCSD method

Different developments in chemistry and emerging technologies have generated a renewed interest in the properties of molecular excited states. We present herein the partition of black-box, size-consistent equation-of-motion coupled cluster singles and doubles (EOM-CCSD) excitation energies within the framework of the interacting quantum atoms (IQA) formalism. We denote this method as IQA/EOM-CCSD. We illustrate this approach by considering small molecules used often in the study of excited states. This investigation shows how the combination of IQA and EOM-CCSD may provide valuable insights into the molecular changes induced by electron excitation via the real space distribution of the energy of an absorbed photon in a molecular system. Our results reveal (i) the most energetically deformed atomic basins and (ii) the most affected covalent and non-covalent interactions within a molecule due to a given photoexcitation. In other words, this kind of analysis provides insights into the spatial energetic redistribution accompanying an electronic excitation, with interesting foreseeable applications in the rational design of photoexcitations with tailored chemical effects. Altogether, we expect that the IQA/EOM-CCSD excitation energy partition will prove useful in the understanding of systems and processes of interest in photophysics and photochemistry.

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.

Perspective: Found in translation: Quantum chemical tools for grasping non-covalent interactions

Today's quantum chemistry methods are extremely powerful but rely upon complex quantities such as the massively multidimensional wavefunction or even the simpler electron density. Consequently, chemical insight and a chemist's intuition are often lost in this complexity leaving the results obtained difficult to rationalize. To handle this overabundance of information, computational chemists have developed tools and methodologies that assist in composing a more intuitive picture that permits better understanding of the intricacies of chemical behavior. In particular, the fundamental comprehension of phenomena governed by non-covalent interactions is not easily achieved in terms of either the total wavefunction or the total electron density, but can be accomplished using more informative quantities. This perspective provides an overview of these tools and methods that have been specifically developed or used to analyze, identify, quantify, and visualize non-covalent interactions. These include the quantitative energy decomposition analysis schemes and the more qualitative class of approaches such as the Non-covalent Interaction index, the Density Overlap Region Indicator, or quantum theory of atoms in molecules. Aside from the enhanced knowledge gained from these schemes, their strengths, limitations, as well as a roadmap for expanding their capabilities are emphasized.

Electron-Pair Distribution in Chemical Bond Formation

The chemical formation process has been studied from relaxation holes, Δh(u), resulting from the difference between the radial intracule density and the nonrelaxed counterpart, which is obtained from atomic radial intracule densities and the pair density constructed from the overlap of the atomic densities. Δh(u) plots show that the internal reorganization of electron pairs prior to bond formation and the covalent bond formation from electrons in separate atoms are completely recognizable processes from the shape of the relaxation hole, Δh(u). The magnitude of Δh(u), the shape of Δh(u) ∀ u < R_eq, and the distance between the minimum and the maximum in Δh(u) provide further information about the nature of the chemical bond formed. A computational affordable approach to calculate the radial intracule density from approximate pair densities has been also suggested, paving the way to study electron-pair distributions in larger systems.

Cooperative and anticooperative effects in resonance assisted hydrogen bonds in merged structures of malondialdehyde

We analyzed non-additive effects in resonance assisted hydrogen bonds (RAHBs) in different β-enolones, which are archetypal compounds of these types of interactions. For this purpose, we used (i) potential energy curves to compute the formation energy, ΔE, of the RAHBs of interest in different circumstances along with (ii) tools offered by quantum chemical topology, namely, the Quantum Theory of Atoms In Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) electronic energy partition. We established the effect that a given H-bond exerts over ΔE associated with another RAHB, determining in this way the cooperativity or the anticooperativity of these interactions. The mesomeric structures and the QTAIM delocalisation indices are consistent with the determined cooperative or anticooperative character of two given RAHBs. The HB cooperativity and anticooperativity studied herein are directly reflected in the IQA interaction energy E, but they are modulated by the surrounding hydrocarbon chain. The IQA decomposition of ΔE_coop, a measure of the cooperativity between a pair of interacting RAHBs, indicates that the analyzed H-bond cooperative/anticooperative effects are associated with greater/smaller (i) strengthening of the pseudo-bicyclic structure of the compounds of interest and (ii) electron localisations with its corresponding changes in the intra and intermolecular exchange-correlation contributions to ΔE. Overall, we expect that this investigation will provide valuable insights into the interplay among hydrogen bonded atoms and the π system in RAHBs contributing in this way to the understanding of the general features of H-bonds.

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.