Revisiting the foundations of the quantum theory of atoms in molecules: Toward a rigorous definition of topological atoms

An explicit classification of consistent variational constraints within the context of the “quantum theory of proper open subsystems” as well as the “quantum theory of atoms in molecules” (QTAIM) it presented. It is demonstrated that the general variational procedure is not sensitive enough to discriminate between different mathematically consistent variational conditions. The uniqueness of the regional kinetic energy is employed to derive the net zero‐flux condition and the regions satisfying this condition are named as quantum divided basins. A modified form of the local zero‐flux is proposed in order to define topological atoms within the context of the orthodox QTAIM. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

The quantum theory of atoms in positronic molecules: The subsystem variational procedure

This contribution deals with the subsystem variational procedure within the context of the quantum theory of atoms in positronic molecules (QTAIPM). Before introducing the subsystem energy functional termed as joint subsystem energy functional, a novel notation and the combination strategy are disclosed in detail by restating the positronic subsystem hypervirial theorem. They are employed in proposing the proper subsystem energy functional, the validity of which is checked by various criteria. The zero flux surfaces of the joint density distribution are used to define the topological atoms in the positronic molecules, and they are incorporated into the subsystem variational procedure as proper real space boundary conditions. The variational procedure finally yields the flux of the joint current property density that also appears in the positronic subsystem hypervirial theorem. At every stage, the corresponding equations for the purely electronic systems within the context of the quantum theory of atoms in molecules (QTAIM) are presented to clearly reveal the analogy between these two formalisms and to emphasize the importance of combining the property density distributions in the QTAIPM. The presented material demonstrates the internal consistency of the whole framework and discloses the fact that the QTAIM must be regarded as a variant of the QTAIPM. Furthermore, this formalism promises an extended QTAIM, which is hoped to resolve the issue of molecular structure beyond the clamp nuclei approximation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Revisiting the foundations of the quantum theory of atoms in molecules: Some open problems

In a series of papers in the last 10 years, various aspects of the mathematical foundations of the quantum theory of atoms in molecules have been considered by this author and his coworkers in some details. Although these considerations answered part of the questions raised by some critics on the mathematical foundations of the quantum theory of atoms in molecules, however, new mathematical problems also emerged during these studies that were reviewed elsewhere [Sh. Shahbazian Int. J. Quantum Chem. 2011 , 111, 4497.]. Beyond mathematical subtleties of the formalism that were the original motivation for initial exchanges and disputes, the questions raised by critics had a constructive effect and prompted the author to propose a novel extension of the theory, now called the multi‐component quantum theory of atoms in molecules [M. Goli, Sh. Shahbazian Theor. Chem. Acc. 2013 , 132, 1365.]. Taking this background into account, in this paper a new set of open problems is put forward that the author believes proper answers to these questions, may open new doors for future theoretical developments of the quantum theory of atoms in molecules. Accordingly, rather than emphasizing on the rigorous mathematical formulation, the practical motivations behind proposing these questions are discussed in detail and the relevant literature are reviewed while when possible, evidence and routes to answers are also provided. The author hopes that proposing these open questions as a compact package may motivate more mathematically oriented people to participate in future developments of the quantum theory of atoms in molecules and its multi‐component version.     

The quantum theory of atoms in positronic molecules: A case study on diatomic species

This article presents the first systematic study of a series of diatomic positronic species using the recently proposed regional approach: the quantum theory of atoms in positronic molecules (QTAIPM). This survey includes the LiH,e⁺, NaH,e⁺, LiF,e⁺, NaF,e⁺, BeO,e⁺, MgO,e⁺, CN^?,e⁺, and OH^?,e⁺ species as typical examples. The computational algorithm of the whole analysis is communicated and reviewed in detail. The topological analysis of the joint density distribution reveals topological structures similar to those observed for the purely electronic systems; that is, each system decomposes into two quantum atoms. By considering some of the regional properties of these quantum atoms, it is demonstrated that the positron affects them seriously through two different mechanisms: direct and indirect contributions, the latter refers to electronic and geometric relaxations. The computational results clearly reveal the fact that the regional properties of the quantum atoms of positronic molecules are not deducible from their purely electronic counterparts; thus, an independent analysis is required for each positronic molecule. The positronic population is considered as a typical regional property showing that the attachment of a positron to a purely electronic system enhances the polarization of the electronic distribution. The concept of regional positron affinities is also introduced and discussed as a nonroutine application of the QTAIPM. The results of this article set the stage for further study on the quantum atoms of polyatomic positronic species. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Biphenyl: A stress tensor and vector‐based perspective explored within the quantum theory of atoms in molecules

We use quantum theory of atoms in molecules (QTAIM) and the stress tensor topological approaches to explain the effects of the torsion φ of the C‐C bond linking the two phenyl rings of the biphenyl molecule on a bond‐by‐bond basis using both a scalar and vector‐based analysis. Using the total local energy density H( r _b), we show the favorable conditions for the formation of the controversial H–H bonding interactions for a planar biphenyl geometry. This bond‐by‐bond QTAIM analysis is found to be agreement with an earlier alternative QTAIM atom‐by‐atom approach that indicated that the H–H bonding interaction provided a locally stabilizing effect that is overwhelmed by the destabilizing role of the C‐C bond. This leads to a global destabilization of the planar biphenyl conformation compared with the twisted global minimum. In addition, the H( r _b) analysis showed that only the central torsional C‐C bond indicated a minimum for a torsion φ value coinciding with that of the conventional global energy minimum. The H–H bonding interactions are found to be topologically unstable for any torsion of the central C‐C bond away from the planar biphenyl geometry. Conversely, we demonstrate that for 0.0° < φ < 39.95° there is a resultant increase in the topological stability of the C nuclei comprising the central torsional C‐C bond. Evidence is found of the effect of the H–H bonding interactions on the torsion φ of the central C‐C bond of the biphenyl molecule in the form of the QTAIM response β of the total electronic charge density ρ( r _b). Using a vector‐based treatment of QTAIM we confirm the presence of the sharing of chemical character between adjacent bonds. In addition, we present a QTAIM interpretation of hyperconjugation and conjugation effects, the former was quantified as larger in agreement with molecular orbital (MO) theory. The stress tensor and the QTAIM H atomic basin path set areas are independently found to be new tools relevant for the incommensurate gas to solid phase transition occurring in biphenyl for a value of the torsion reaction coordinate φ ≈ 5°. © 2015 Wiley Periodicals, Inc.     

The general setting for the zero‐flux condition: The lagrangian and zero‐flux conditions that give the heisenberg equation of motion

Generalizing our recent work on relativistic generalizations of the quantum theory of atoms in molecules, we present the general setting under which the principle of stationary action for a region leads to open quantum subsystems. The approach presented here is general and works for any Hamiltonian, and when a reasonable Lagrangian is selected, it often leads to the integral of the Laplacian of the electron density on the region vanishing as a necessary condition for the zero‐flux surface. Alternatively, with this method, one can design a Lagrangian that leads to a surface of interest (though this Lagrangian may not be, and indeed probably will not be, “reasonable”). For any reasonable Lagrangian for the electronic wave function and any two‐component method (related by integration by parts to the Hamiltonian) considered, the Bader definition of an atom is recaptured. © 2018 Wiley Periodicals, Inc.     

The quantum divided basins: A new class of quantum subsystems

The rigorous theory of the quantum divided basins (QDB), the quantum subsystems emerging from the net zero‐flux equation, is considered in this article. This framework, the quantum theory of proper open subsystems, is derived from the extension of the quantum theory of atoms in molecules to encompass the new class of quantum subsystems. It is demonstrated that the regional hypervirial theorem and the associated regional observables as well as the subsystem variational procedure are all expressible for the QDB. The history of QDB is briefly reviewed and the bundles, which were introduced by other researchers, are offered as typical examples whereas new examples of QDB (not yet mentioned in literature) are also presented. Based on some model systems as well as the nitrogen molecule, the regional properties and the morphologies of typical QDB are scrutinized. It is also demonstrated that the QDB may be used to study the fine structure of the electron localization and delocalization. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms

Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra‐atomic and inter‐atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra‐atomic and inter‐atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange‐correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange‐correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc.     

Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions

Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺ together with (ii) the anions F⁻, Cl⁻, Br⁻, NO₃⁻ and SO₄²⁻ with water clusters (H₂O)_n, n=1–8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.     

Two-electron integrations in the quantum theory of atoms in molecules with correlated wave functions

A recent method proposed to compute two-electron integrals over arbitrary regions of space [Martin Pendas, A. et al., J Chem Phys 2004, 120, 4581] is extended to deal with correlated wave functions. To that end, we use a monadic factorization of the second-order reduced density matrix originally proposed by E. R. Davidson [Chem Phys Lett 1995, 246, 209] that achieves a full separation of the interelectronic components into one-electron terms. The final computational effort is equivalent to that found in the integration of a one determinant wave function with as many orbitals as occupied functions in the correlated expansion. Similar strategies to extract the exchange and self-interaction contributions from the two-electron repulsion are also discussed, and several numerical results obtained in a few test systems are summarized.     

Quasi‐random integration in quantum chemistry: Efficiency,stability, and application to the study of small atoms and molecules constrained in spherical boxes

This project consists of two parts. In the first part, a series of test calculations is performed to verify that the integrals involved in the determination of atomic and molecular properties by standard self‐consistent field (SCF) methods can be obtained through Halton, Korobov, or Hammersley quasi‐random integration procedures. Through these calculations, we confirm that all three methods lead to results that meet the levels of precision required for their use in the calculation of properties of small atoms or molecules at least at a Hartree–Fock level. Moreover, we have ensured that the efficiency of quasi‐random integration methods that we have tested is Halton=Korobov>Hammersley?pseudo‐random. We also find that these results are comparable to those yielded by ordinary Monte Carlo (pseudo‐random) integration, with a calculation effort of two orders of smaller magnitude. The second part, which would not have been possible without the integration method previously analyzed, contains a first study of atoms constrained in spherical boxes through SCF calculations with basis functions adapted to the features of the problem: Slater‐type orbitals (STOs) trimmed by multiplying them by a function that yields 1 for 0 < r < (R‐δ), polynomial values for (R‐δ) < r < R and null for r > R, R being the radius of the box and δ a variationally determined interval. As a result, we obtain a equation of state for electrons of small systems, valid just in the limit of low temperatures, but fairly simple. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A QTAIM-based energy partitioning for understanding the physical origin of conformational preferences: Application to the Z effect in O=C-X-R and related units

A quantum theory of atoms in molecules-based energy partitioning was carried out for Z and E conformers of a series of O=C-X-R containing compounds. The results obtained for the simplest compound (formic acid) indicate that the attraction of the electron density within carbonyl oxygen by the nucleus of the acid hydrogen is the most important energy term for Z preference. This conclusion can be extended (mutatis mutandis) to larger carboxylic acids, esters, sulfur derivatives, secondary amides, and carbonyl isocyanates, and even explains the sequence of relative conformational energies in the HCXOH series (X = O, S, Se). In contrast, although the hyperconjugative model has been traditionally employed to explain this preference, we observe it is incompatible with: (i) relative values of diverse QTAIM atomic populations for the Z / E conformational equilibrium; (ii) conformational energies in the HCXOH series. © 2012 Wiley Periodicals, Inc.     

Energy decomposition in the topological theory of atoms in molecules and in the linear combination of atomic orbitals formalism: a note

 It is shown that the molecular energy calculated at the self-consistent-field level can be strictly expressed as a sum of one- and two-atom energy components in the framework of Bader's topological theory of atoms in molecules (AIM). The expressions of our recent “chemical energy component analysis” can be obtained from the AIM ones as some linear combination of atomic orbitals mappings of the integrations over the atomic basins. Received: 15 June 2000 / Accepted: 4 October 2000 / Published online: 19 January 2001     

The Analysis of Os≡C Bond and Electric Field Influence on the Properties in the Osmium Carbyne Complex OsCl3 (≡CCH2CMe3 )(PH3 )2: A Theoretical Insight

In this paper, we study the effect of electric field on the dipole moment, electronic structure, and frontier orbital energy in the osmium carbyne complex OsCl₃ (≡CCH₂CMe₃ )(PH₃ )₂ using MPW1PW91 quantum chemical computations. We demonstrate the nature of the chemical bond between the [OsCl₃ (PH₃ )₂]⁻ and [CCH₂CMe₃ ]⁺ fragments through energy and charge decomposition analyses. We also estimate the percentage composition in terms of the specified groups of frontier orbitals for this complex to investigate the feature in the metal–ligand bonds. Quantum theory of atoms in molecules (QTAIM ) is applied to elucidate the Os≡C bond in the complex. Also, the influence of external electric field on the energy, frontier orbital energies, and HOMO–LUMO gap values is analyzed.