Fermi and Coulomb correlation effects upon the interacting quantum atoms energy partition

How Many Electrons Does a Molecular Electride Hold?

Electrides are very peculiar ionic compounds where electrons occupy the anionic positions. In a crystal lattice, these isolated electrons often form channels or surfaces, furnishing electrides with many traits with promising technological applications. Despite their huge potential, thus far, only a few stable electrides have been produced because of the intricate synthesis they entail. Due to the difficulty in assessing the presence of isolated electrons, the characterization of electrides also poses some serious challenges. In fact, their properties are expected to depend on the arrangement of these electrons in the molecule. Among the criteria that we can use to characterize electrides, the presence of a non-nuclear attractor (NNA) of the electron density is both the rarest and the most salient feature. Therefore, a correct description of the NNA is crucial to determine the properties of electrides. In this paper, we analyze the NNA and the surrounding region of nine molecular electrides to determine the number of isolated electrons held in the electride. We have seen that the correct description of a molecular electride hinges on the electronic structure method employed for the analyses. In particular, one should employ a basis set with sufficient flexibility to describe the region close to the NNA and a density functional approximation that does not suffer from large delocalization errors. Finally, we have classified these nine molecular electrides according to the most likely number of electrons that we can find in the NNA. We believe this classification highlights the strength of the electride character and will prove useful in designing new electrides.

From the Linnett-Gillespie model to the polarization of the spin valence shells of metals in complexes

In this paper, we present a novel approach to track the origin of the metal complex structure from the topology of the α and β spin densities as an extension of the Linnett-Gillespie model. Usually, the theories that explain the metal-ligand interactions consider the disposition and the relative energies of the empty or occupied set of d orbitals, ignoring the spin contribution explicitly. Our quantum topological approach considers the spatial distribution of the α and β spin valence shells, and the energy interaction between them. We used the properties of the atomic graph, a topological object that summarises the charge concentrations and depletions on the valence shell of an atom in a molecule, and the interacting quantum atoms (IQA) energy partition scheme. Unlike the Linnett-Gillespie model, which is based on electron-electron repulsion, our approach states that the ligands provoke a redistribution of the electron density to maximize the nuclear-electron interactions in each spin valence shell to bypass the concentration of electron-electron interactions, resulting in a polarization pattern which determines the position of the ligands.

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.

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.

The Coulomb Hole of the Ne Atom

We analyze the Coulomb hole of Ne from highly-accurate CISD wave functions obtained from optimized even-tempered basis sets. Using a two-fold extrapolation procedure we obtain highly accurate results that recover 97 % of the correlation energy. We confirm the existence of a shoulder in the short-range region of the Coulomb hole of the Ne atom, which is due to an internal reorganization of the K-shell caused by electron correlation of the core electrons. The feature is very sensitive to the quality of the basis set in the core region and it is not exclusive to Ne, being also present in most of second-row atoms, thus confirming that it is due to K-shell correlation effects.

A relative energy gradient (REG) study of the planar and perpendicular torsional energy barriers in biphenyl

Biphenyl is a prototype molecule, the study of which is important for a proper understanding of stereo-electronic effects. In the gas phase it has an equilibrium central torsion angle of ~ 45° and shows both a planar (0°) and a perpendicular (90°) torsional energy barrier. The latter is analysed for the first time. We use the newly proposed REG method, which is an exhaustive procedure that automatically ranks atomic energy contributions according to their importance in explaining the energy profile of a total system. Here, the REG method operates on energy contributions computed by the interacting quantum atoms method. This method is minimal in architecture and provides a crisp picture of well-defined and well-separated electrostatic, steric and exchange (covalent) energies at atomistic level. It is shown that the bond critical point occurring between the ortho-hydrogens in the planar geometry has been wrongly interpreted as a sign of repulsive interaction. A convenient metaphor of analysing football matches is introduced to clarify the role of a REG analysis.

Revitalizing the concept of bond order through delocalization measures in real space

Ab initio quantum chemistry is an independent source of information supplying an ever widening group of experimental chemists. However, bridging the gap between these ab initio data and chemical insight remains a challenge. In particular, there is a need for a bond order index that characterizes novel bonding patterns in a reliable manner, while recovering the familiar effects occurring in well-known bonds. In this article, through a large body of calculations, we show how the delocalization index derived from Quantum Chemical Topology (QCT) serves as such a bond order. This index is defined in a parameter-free, intuitive and consistent manner, and with little qualitative dependency on the level of theory used. The delocalization index is also able to detect the subtler bonding effects that underpin most practical organic and inorganic chemistry. We explore and connect the properties of this index and open the door for its extensive usage in the understanding and discovery of novel chemistry.

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.

Comprehensive benchmarking of density matrix functional approximations

The energy usually serves as a yardstick in assessing the performance of approximate methods in computational chemistry. After all, these methods are mostly used for the calculation of the electronic energy of chemical systems. However, computational methods should be also aimed at reproducing other properties, such strategy leading to more robust approximations with a wider range of applicability. In this study, we suggest a battery of ten tests with the aim to analyze density matrix functional approximations (DMFAs), including several properties that the exact functional should satisfy. The tests are performed on a model system with varying electron correlation, carrying a very small computational effort. Our results not only put forward a complete and exhaustive benchmark test for DMFAs, currently lacking, but also reveal serious deficiencies of existing approximations that lead to important clues in the construction of more robust DMFAs.