Atoms in molecules in real space: a fertile field for chemical bonding.

In this perspective, we review some recent advances in the concept of atoms-in-molecules from a real space perspective. We first introduce the general formalism of atomic weight factors that allows unifying the treatment of fuzzy and non-fuzzy decompositions under a common algebraic umbrella. We then show how the use of reduced density matrices and their cumulants allows partitioning any quantum mechanical observable into atomic or group contributions. This circumstance provides access to electron counting as well as energy partitioning, on the same footing. We focus on how the fluctuations of atomic populations, as measured by the statistical cumulants of the electron distribution functions, are related to general multi-center bonding descriptors. Then we turn our attention to the interacting quantum atom energy partitioning, which is briefly reviewed since several general accounts on it have already appeared in the literature. More attention is paid to recent applications to large systems. Finally, we consider how a common formalism to extract electron counts and energies can be used to establish an algebraic justification for the extensively used bond order-bond energy relationships. We also briefly review a path to recover one-electron functions from real space partitions. Although most of the applications considered will be restricted to real space atoms taken from the quantum theory of atoms in molecules, arguably the most successful of all the atomic partitions devised so far, all the take-home messages from this perspective are generalizable to any real space decompositions.

The Ehrenfest force field: A perspective based on electron density functions.

The topology of the Ehrenfest force field (EhF) is investigated as a tool for describing local interactions in molecules and intermolecular complexes. The EhF is obtained by integrating the electronic force operator over the coordinates of all but one electron, which requires knowledge of both the electron density and the reduced pair density. For stationary states, the EhF can also be obtained as minus the divergence of the kinetic stress tensor, although this approach leads to well-documented erroneous asymptotic behavior at large distances from the nuclei. It is shown that these pathologies disappear using the electron density functions and that the EhF thus obtained displays the correct behavior in real space, with no spurious critical points or attractors. Therefore, its critical points can be unambiguously obtained and classified. Test cases, including strained molecules, isomerization reactions, and intermolecular interactions, were analyzed. Various chemically relevant facts are highlighted: for example, non-nuclear attractors are generally absent, potential hydrogen-hydrogen interactions are detected in crowded systems, and a bifurcation mechanism is observed in the isomerization of HCN. Moreover, the EhF atomic basins are less charged than those of the electron density. Although integration of the EhF over regions of real space can also be performed to yield the corresponding atomic forces, several numerical drawbacks still need to be solved if electron density functions are to be used for that purpose. Overall, the results obtained support the Ehrenfest force field as a reliable descriptor for the definition of atomic basins and molecular structure.

QM/MM Energy Decomposition Using the Interacting Quantum Atoms Approach.

The interacting quantum atoms (IQA) method decomposes the quantum mechanical (QM) energy of a molecular system in terms of one- and two-center (atomic) contributions within the context of the quantum theory of atoms in molecules. Here, we demonstrate that IQA, enhanced with molecular mechanics (MM) and Poisson-Boltzmann surface-area (PBSA) solvation methods, is naturally extended to the realm of hybrid QM/MM methodologies, yielding intra- and inter-residue energy terms that characterize all kinds of covalent and noncovalent bonding interactions. To test the robustness of this approach, both metal-water interactions and QM/MM boundary artifacts are characterized in terms of the IQA descriptors derived from QM regions of varying size in Zn(II)- and Mg(II)-water clusters. In addition, we analyze a homologous series of inhibitors in complex with a matrix metalloproteinase (MMP-12) by carrying out QM/MM-PBSA calculations on their crystallographic structures followed by IQA energy decomposition. Overall, these applications not only show the advantages of the IQA QM/MM approach but also address some of the challenges lying ahead for expanding the QM/MM methodology.

Partition of the electronic energy of the PM7 method via the interacting quantum atoms approach.

Partitions of the electronic energy such as that provided by the Interacting Quantum Atoms (IQA) approach have given valuable insights for numerous chemical systems and processes. Unfortunately, this kind of analysis may involve the integration of scalar fields over very irregular volumes, a condition which leads to a large and often prohibitive computational effort. These circumstances have limited the use of these energy partitions to systems comprising a few tens of atoms at most. On the other hand, semiempirical methods have proved useful in the study of systems of several thousands of atoms. Therefore, the goal of this work is to carry out partitions of the semiempirical method PM7 in compliance with the IQA approach. For this purpose, we computed one- and two-atomic energetic contributions whose sum equals the PM7 electronic energy. We illustrate how one might exploit the partition of electronic energies computed via the PM7 method by considering small organic and inorganic molecules and the energetics of individual hydrogen bond interactions within several water clusters which include (H2O)30, (H2O)50 and (H2O)100. We also considered the solvation of the amphiphilic caprylate anion to exemplify how to exploit the energy partition proposed in this paper. Overall, this investigation shows how the approach put forward herein might give further insights of the interactions occurring within complex systems in physical and biological chemistry.

Atomic shell structure from Born probabilities: Comparison to other shell descriptors and persistence in molecules.

Real space chemical bonding descriptors, such as the electron localization function or the Laplacian of the electron density, have been widely used in electronic structure theory thanks to their power to provide chemically intuitive spatial images of bonded and non-bonded interactions. This capacity stems from their ability to display the shell structure of atoms and its distortion upon molecular formation. Here, we examine the spatial position of the N electrons of an atom at the maximum of the square of the wavefunction, the so-called Born maximum, as a shell structure descriptor for ground state atoms with Z = 1-36, comparing it to other available indices. The maximization is performed with the help of variational quantum Monte Carlo calculations. We show that many electron effects (mainly Pauli driven) are non-negligible, that Born shells are closer to the nucleus than any other of the examined descriptors, and that these shells are very well preserved in simple molecules.

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.

Lewis Structures from Open Quantum Systems Natural Orbitals: Real Space Adaptive Natural Density Partitioning.

Building chemical models from state-of-the-art electronic structure calculations is not an easy task, since the high-dimensional information contained in the wave function needs to be compressed and read in terms of the accepted chemical language. We have already shown ( Phys. Chem. Chem. Phys. 2018, 20, 21368) how to access Lewis structures from general wave functions in real space by reformulating the adaptive natural density partitioning (AdNDP) method proposed by Zubarev and Boldyrev ( Phys. Chem. Chem. Phys. 2008, 10, 5207). This provides intuitive Lewis descriptions from fully orbital invariant position space descriptors but depends on not immediately accessible higher order cumulant density matrices. By using an open quantum systems (OQS) perspective, we here show that the rigorously defined OQS fragment natural orbitals can be used to build a consistent real space adaptive natural density partitioning based only on spatial information and the system's one-particle density matrix. We show that this rs-AdNDP approach is a cheap, efficient, and robust technique that immerses electron counting arguments fully in the real space realm.

Interacting Quantum Atoms Method for Crystalline Solids.

An implementation of the Interacting Quantum Atoms method for crystals is presented. It provides a real space energy decomposition of the energy of crystals in which all energy components are physically meaningful. The new package ChemInt enables one to compute intra-atomic and inter-atomic energies, as well as electron population measures used for quantitative description of chemical bonds in crystals. The implementation is tested and applied to characteristic molecular and crystalline systems with different types of bonding.

Efficient implementation of the interacting quantum atoms energy partition of the second-order Møller-Plesset energy.

We describe an efficient implementation of the partition of the second-order Møller-Plesset (MP2) correlation energy within the interacting quantum atoms (IQA) energy decomposition. We simplify the IQA integration bottleneck by considering only the occupied to virtual elements of the second order reduced density matrix, a procedure that reduces substantially the size of the two-electron matrix, which has to be addressed. The algorithmic improvements described herein allow to perform the decomposition of the MP2 correlation energy for medium size molecular systems using moderate computational resources. We expect that the methods developed in this investigation will prove useful to understand electron correlation effects through a real space perspective.

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 (H2 O)2 . 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.

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

Invited for the cover of this issue is Alberto Fernández-Alarcón and co-workers at The Institute of Chemistry of the National Autonomous University of Mexico and The School of Chemistry of the University of Oviedo. The image depicts the real space analysis of the excitation energies in the double blue and red shift of the water dimer. Read the full text of the article at 10.1002/chem.202002854.

Bond Order Densities in Real Space.

In this contribution we introduce the concept of bond order density (BOD) on the basis of a previous work on natural adaptive orbitals. We show that BODs may be used to visualize both the global spatial distribution of the covalent bond order and its eigencomponents, which we call bond(ing) channels. BODs can be equally computed at correlated and noncorrelated levels of theory and in ground or excited states, thus offering an appealing description of bond-forming, bond-breaking, and bond-evolution processes. We show the power of the approach by examining a number of homo- and heterodiatomics, including the controversial existence of a fourth bonding component in dicarbon, by analyzing a few interesting bonding situations in polyatomics and chemical transformations, and by exemplifying exotic bonding behaviors in simple excited electronic states.

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.

On Electrostatics, Covalency, and Chemical Dashes: Physical Interactions versus Chemical Bonds.

The increasing availability of real-space interaction energies between quantum atoms or fragments that provide a chemically intuitive decomposition of intrinsic bond energies into electrostatic and covalent terms [see, for instance, Chem. Eur. J. 2018, 24, 9101] provides evidence for differences between the physicist's concept of interaction and the chemist's concept of a bond. Herein, it is argued that, for the former, all types of interactions are treated equally, whereas, for the latter, only the covalent short-range interactions have actually been used to build intuition about chemical graphs and chemical bonds. This has led to the bonding role of long-range Coulombic terms in molecular chemistry being overlooked. Simultaneously, blind consideration of electrostatic terms in chemical bonding parlance may lead to confusion. The relationship between these concepts is examined herein, and some notes of caution on how to merge them are proposed.

Exotic Bonding Regimes Uncovered in Excited States.

Real-space tools were employed to show that the chemical bonding scenario used routinely to understand ground states lacks the necessary flexibility in excited states. It is shown that, even for two-center, two-electron bonds, the real-space bond orders have exotic values that have never been reported. The nature of these situations was uncovered by using electron-counting techniques that provide an appealing statistical interpretation of bonding descriptors, together with simple physical models. Bond orders greater than one as well as negative bond orders for a single bonding electron pair emerge in situations in which the electrons in the pair show a gregarious (bosonic) instead of the usual lonely (fermionic) behavior. In the first case the gregarious pair is intra-atomic, whereas the coupling is interatomic in the second. A number of examples are used to substantiate these claims.

Chemical Bonding from the Statistics of the Electron Distribution.

An introduction to the theory of chemical bonding from the point of view of the statistics of the electron distribution is presented. When atoms bind to form a molecule, their originally fixed number of electrons ceases to be a well-defined observable, and this implies that their in-the-molecule electron populations fluctuate. If a chemically meaningful definition of an atom in a molecule is assumed, the probabilities of finding a given number of electrons in each of the atoms comprising the molecule can be computed. We show in this review how the complete electron distribution function (EDF) can be used to reconstruct the basic concepts and quantities used in chemical bonding without recourse to the orbital paradigm. From the statistical point of view, which inherits Born's probabilistic interpretation of quantum mechanics, a set of atoms are bonded when their electron populations are mutually dependent. We quantify this statistical dependence by the cumulant moments of the EDF, which provide a consistent description of both two- and multi-center bonding. Particular attention is paid to building EDFs from model wavefunctions. With this, a simple bridge with orbital thinking is built. The statistical interpretation allows to easily classify all possible bonds of a given kind. We show that there are vast unexplored territories that should receive due consideration. Although building EDFs from models is easy and very instructive, the contrary is considerably more difficult. Recipes to extract chemical information from computed EDFs are also reviewed and, in all the cases, simple toy systems are used to show how the methodology works, allowing non-experts to follow easily the presentation.

Fluorine conformational effects characterized by energy decomposition analysis.

Electrostatic and stereoelectronic effects associated with fluorine atoms can be exploited as conformational tools for the design of shape-controlled functional molecules. To gain further insight into the nature and strength of these effects, we use the Interacting Quantum Atoms (IQA) method augmented with the semiclassical pairwise dispersion potential to decompose the conformational energies of fluoro-substituted molecules into fragment-based energy contributions, which include deformation/distortion terms and the electrostatic, exchange-correlation and dispersion interactions. The studied molecules comprise various F-CH2-CH2-X and F-CH2-CO-X systems, as well as selected conformers of an α,ß-difluoro-γ-amino-acid derivative that is potentially useful for the design of shape-controlled bioactive amino acids and peptides. We identify the most relevant exchange-correlation and/or electrostatic interaction terms contributing to the stability of the various conformers, and we show that IQA can be used to assess the gauche/anti or trans/cis preferences in molecules with two or more rotatable bonds as well as to study the roles played by other concomitant effects (e.g., CH/OH/NHF contacts). For the α,ß-difluoro-γ-amino acid derivatives, our theoretical analysis indicates that the gauche/anti and trans/cis effects associated with fluorine bonds can be significantly attenuated by other specific intra-molecular contacts.

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.

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) (CH3F, CH4, CO2, CS2, SiO2, SiH3F, SiH4, GeH3F, GeO2, and GeH4) and an electron-pair-rich system (B) (NCH, NCO-, OCN-, F-, Br-, CN-, CO, CS, Kr, NC-, NH3, OC, OH2, SH-, and N3-) 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 .