Performance of the density matrix functional theory in the quantum theory of atoms in molecules

Electron density-based protocol to recover the interacting quantum atoms components of intermolecular binding energy

A new approach for obtaining interacting quantum atoms-defined components of binding energy of intermolecular interactions, which bypasses the use of standard six-dimensional integrals and two-particle reduced density matrix (2-RDM) reconstruction, is proposed. To examine this approach, three datasets calculated within the density functional theory framework using the def2-TZVP basis have been explored. The first two, containing 53 weakly bound bimolecular associates and 13 molecular clusters taken from the crystal, were used in protocol refinement, and the third one containing other 20 bimolecular and three cluster systems served as a validation reference. In addition, to verify the performance of the proposed approach on an exact 2-RDM, calculations within the coupled cluster formalism were performed for part of the first set systems using the cc-pVTZ basis set. The process of optimization of the proposed parametric model is considered, and the role of various energy contributions in the formation of non-covalent interactions is discussed with regard to the obtained trends.

A fast approximate extension of the interacting quantum atoms energy decomposition to excited states

An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact V_xc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.

Collective interactions among organometallics are exotic bonds hidden on lab shelves

Recent discovery of an unusual bond between Na and B in NaBH3- motivated us to look for potentially similar bonds, which remained unnoticed among systems isoelectronic with NaBH3-. Here, we report a novel family of collective interactions and a measure called exchange-correlation interaction collectivity index (ICIXC; [Formula: see text]) to characterize the extent of collective versus pairwise bonding. Unlike conventional bonds in which ICIXC remains close to one, in collective interactions ICIXC may approach zero. We show that collective interactions are commonplace among widely used organometallics, as well as among boron and aluminum complexes with the general formula [Ma+AR3]b- (A: C, B or Al). In these species, the metal atom interacts more efficiently with the substituents (R) on the central atoms than the central atoms (A) upon forming efficient collective interactions. Furthermore, collective interactions were also found among fluorine atoms of XFn systems (X: B or C). Some of organolithium and organomagnesium species have the lowest ICIXC among the more than 100 studied systems revealing the fact that collective interactions are rather a rule than an exception among organometallic species.

Planar Hypercoordinate Carbons in Alkali Metal Decorated CE3 2- and CE2 2- Dianions

After exploring the potential energy surfaces of M_m CE₂ ^p (E=S-Te, M=Li-Cs, m=2, 3 and p=m-2) and M_n CE₃ ^q (E=S-Te, M=Li-Cs, n=1, 2, q=n-2) combinations, we introduce 38 new global minima containing a planar hypercoordinate carbon atom (24 with a planar tetracoordinate carbon and 14 with a planar pentacoordinate carbon). These exotic clusters result from the decoration of V-shaped CE₂ ^2- and Y-shaped CE₃ ^2- dianions, respectively, with alkali counterions. All these 38 systems fulfill the geometrical and electronic criteria to be considered as true planar hypercoordinate carbon systems. Chemical bonding analyses indicate that carbon is covalently bonded to chalcogens and ionically connected to alkali metals.

Neither too Classic nor too Exotic: One-Electron Na⋅B Bond in NaBH3 - Cluster

It is here reported that the NaBH₃ ^- cluster exhibits a Na⋅B one-electron bond, a well-established type of electron-deficient bonding in the literature. The topological analysis of the electron localization function, at the correlated level, reveals that Na^- , when approaching the bonding distance, fairly distributes its valence electron pair between two lobes. One of these electrons is used to bond with BH₃ , which participates through its boron empty p-orbital. Furthermore, the bonding situation of LiBH₃ ^- , KBH₃ ^- , MgBH₃ , and CaBH₃ global minima structures are similar to that of NaBH₃ ^- , extending the family of these new one-electron bond systems with biradicaloid character.

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.

Substituent effects in the so-called cationπ interaction of benzene and its boron-nitrogen doped analogues: overlooked role of σ-skeleton

Despite massive efforts to pinpoint the substituent effects in the so-called cationπ systems, no consensus has been yet reached on how substituents exercise their effects in the interaction of the aromatic molecule with the metal ion. The π-polarization (the Hunter model) and the direct local effect (the Wheeler-Houk model) are two lines of thought applied to this problem, but the justification of both approaches is based on insufficiently proven assumptions and approximations. In order to shed more light on this issue we propose a new approach which enables us to gauge directly the energetic trends resulting from the interaction of the ring with the cation. In our method we add one more partitioning level to the interacting quantum atoms (IQA) scheme and decompose the IQA interaction energies into contributions resulting from σ and π electron densities of the aromatic ring. The new approach, which is named partitioned-IQA, abbreviated as p-IQA, has been applied to complexes of derivatives of benzene or azaborine interacting with a sodium cation. The p-IQA approach reveals that in these systems both σ and π electronic moieties are polarized. Interestingly, for the majority of cases the σ-polarization outweighs the π one, contrary to the Hunter model. However, the Wheeler-Houk model is not precise, either, since the σ-polarization shows some degree of non-locality. In addition, the substituents are found to have a negligible influence on the ring orbital-overlapping capability, i.e. the covalency. Therefore, the substituent effect in the cationπ interaction is a nonlocal classical effect, indicating that neither Hunter model nor Wheeler-Houk model is able to fully describe all the aspects of the substituent effects. The p-IQA conclusions for the considered systems have been compared with the results from the functional-group SAPT (F-SAPT) method. We believe that the presented partitioning in the IQA framework will provide a deeper insight into the substituent effects in the cationπ interactions, which is beyond the σ-π atomic charge population separation.

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.

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.

Reliability of interacting quantum atoms (IQA) data computed from post-HF densities: impact of the approximation used

The BBC1 approximation is recommended for IQA calculations; MP2/BBC1 and CCSD/BBC1 produced highly comparable FAMSEC-based interpretations of intramolecular interactions.

Natural occupation numbers in two-electron quantum rings

Natural orbitals (NOs) are central constituents for evaluating correlation energies through efficient approximations. Here, we report the closed-form expression of the NOs of two-electron quantum rings, which are prototypical finite-extension systems and new starting points for the development of exchange-correlation functionals in density functional theory. We also show that the natural occupation numbers for these two-electron paradigms are in general non-vanishing and follow the same power law decay as atomic and molecular two-electron systems.

Forced bonding and QTAIM deficiencies: a case study of the nature of interactions in He@adamantane and the origin of the high metastability

Calculations within the framework of the interacting quantum atoms (IQA) approach have shown that the interactions of the helium atom with both tertiary, tC, and secondary, sC, carbon atoms in the metastable He@adamantane (He@adam) endohedral complex are bonding in nature, whereas the earlier study performed within the framework of Bader's quantum theory of atoms in molecules (QTAIM) revealed that only He---tC interactions are bonding. The He---tC and He---sC bonding interactions are shown to be forced by the high pressure that the helium and carbon atoms exert upon each other in He@adam. The occurrence of a bonding interaction between the helium and sC atoms, which are not linked by a bond path, clearly shows that the lack of a bond path between two atoms does not necessarily indicate the lack of a bonding interaction, as is asserted by QTAIM. IQA calculations showed that not only the destabilization of the adamantane cage, but also a huge internal destabilization of the helium atom, contribute to the metastability of He@adam, these contributions being roughly equal. This result disproves previous opinions based on QTAIM analysis that only the destabilization of the adamantane cage accounts for the endothermicity of He@adam. Also, it was found that there is no homeomorphism of the ρ(r) and -v(r) fields of He@adam. Comparison of the IQA and QTAIM results on the interactions in He@adam exposes other deficiencies of the QTAIM approach. The reasons for the deficiencies in the QTAIM approach are analyzed.

Electron delocalization and aromaticity in low-lying excited states of archetypal organic compounds

Aromaticity is a property usually linked to the ground state of stable molecules. Although it is well-known that certain excited states are unquestionably aromatic, the aromaticity of excited states remains rather unexplored. To move one step forward in the comprehension of aromaticity in excited states, in this work we analyze the electron delocalization and aromaticity of a series of low-lying excited states of cyclobutadiene, benzene, and cyclooctatetraene with different multiplicities at the CASSCF level by means of electron delocalization measures. While our results are in agreement with Baird's rule for the aromaticity of the lowest-lying triplet excited state in annulenes having 4nπ-electrons, they do not support Soncini and Fowler's generalization of Baird's rule pointing out that the lowest-lying quintet state of benzene and septet state of cyclooctatetraene are not aromatic.