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

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.

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.

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.