Two-electron integrations in the quantum theory of atoms in molecules

A method to compute two-electron integrals over arbitrary regions of space is introduced and particularized to the basins appearing in the quantum theory of atoms in molecules. The procedure generalizes the conventional multipolar approach to account for overlapping densities. We show that the approach is always convergent and computationally efficient, scaling as N(4) in the worst, two-center case. Several numerical results supporting our claims are also presented.     

Inductive effect: a quantum theory of atoms in molecules perspective

Substituent effects are ubiquitous in chemistry and the most fundamental is the inductive effect. In this study, the so-called inductive effect was probed in derivatives of bicyclo[1.1.1]pentane-1-carboxylic acid using the isodesmic reaction energy of the acid-base deprotonation, calculated at the PBE0/6-31++G(d,p) level of theory (used throughout). Although structure, molecular orbitals, and nuclear magnetic shielding parameters are discussed, the main focus of this study is the use of the quantum theory of atoms in molecules to analyze the electron density distribution. It was observed that the effect propagates via the manipulation of atomic dipole moments controlled by that of the substituent. As the dipole moment conforms to the principle of atomic transferability, it is found that the substituent dipole determined in simple systems (e.g., R-H) can be used to describe the effect upon the bicyclo[1.1.1]pentane-1-carboxylic acid system.     

Revisiting the calculation of condensed Fukui functions using the quantum theory of atoms in molecules

The analysis of previously reported shortcomings of the condensed Fukui functions obtained making use of the quantum theory of atoms in molecules indicates these drawbacks are due to the inadequacy of the definition employed to compute them and not to the partitioning. A new procedure, which respects the mathematical definition and solves these problems, is presented for the calculation of condensed Fukui functions for atomic basins defined according to the quantum theory of atoms in molecules. It is tested in a set of 18 molecules, which includes the most controversial reported cases.     

Relativistic (SR‐ZORA) quantum theory of atoms in molecules properties

The Quantum Theory of Atoms in Molecules (QTAIM) is used to elucidate the effects of relativity on chemical systems. To do this, molecules are studied using density‐functional theory at both the nonrelativistic level and using the scalar relativistic zeroth‐order regular approximation. Relativistic effects on the QTAIM properties and topology of the electron density can be significant for chemical systems with heavy atoms. It is important, therefore, to use the appropriate relativistic treatment of QTAIM (Anderson and Ayers, J. Phys. Chem. 2009, 115, 13001) when treating systems with heavy atoms. © 2016 Wiley Periodicals, Inc.     

Toward the multi-component quantum theory of atoms in molecules: a variational derivation

The general formalism of an extended quantum theory of atoms in molecules (QTAIM) dealing with the multi-component quantum systems, composed of various types of quantum particles, is disclosed in this contribution. This novel methodology, termed as the multi-component QTAIM (MC-QTAIM), is able to deal with non-adiabatic ab initio wavefunctions extracting atoms in molecules quantifying their properties. It can also be applied to elucidate the AIM structure of exotic species and bound quantum systems consisting of fundamental elementary particles like positrons and muons. The formalism is based on the previously disclosed density combination idea that is extended to derive the multi-component subsystem hypervirial theorem as well as the extended subsystem energy functional. Through the extended subsystem variational procedure, inspired from Schrödinger’s original variational principle, the surface terms containing the flux of the current property densities are derived. Accordingly, the extended Gamma field is introduced during this variational procedure that is used as the basic scalar field in the topological analysis yielding atoms in molecules and their real space boundaries. The Gamma field is central to the MC-QTAIM, replacing the usual one-electron density employed in the orthodox QTAIM and corresponding topological analysis. Through the multi-component hypervirial theorem, various regional theorems are derived which are then used to quantify the mechanical properties of atoms in molecules; these include the force, virial, torque, power, continuity and current theorems. In order to demonstrate the capability of the formalism, isotopically asymmetric hydrogen molecules, HD, HT and DT as well as YX systems (Y = ⁶Li, ⁷Li; X = H, D, T) composed of electrons and two different nuclei, all treated equally as quantum waves instead of clamped particles, are analyzed within context of the MC-QTAIM. The resulting computational analysis demonstrates that the MC-QTAIM is able to yield reasonable topological structures similar to those observed previously for diatomic species within context of the orthodox QTAIM. The asymmetrical nature of these species, inherent in their non-Born–Oppenhiemer wavefunctions, manifests itself clearly in the MC-QTAIM analysis yielding two distinguishable atomic basins with different properties. These differences are rationalized generally by the observed electron transfer from one basin to the other. Finally, some possible future theoretical extensions are considered briefly.     

Revisiting the variational nature of the quantum theory of atoms in molecules

A simplified derivation of the variational nature of the real space basins used in the quantum theory of atoms in molecules (QTAM) is presented. We focus on pointing out the non-standard characteristics of the variational problem that is solved, and on clarifying some points that tend to be misinterpreted. An explicit discussion of the meaning of the functional minimized is also presented and used to derive both a new form for the atomic virial theorem and to show how basin chemical potentials may be defined.     

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

The generalization to arbitrary molecular geometries of the energetic partitioning provided by the atomic virial theorem of the quantum theory of atoms in molecules (QTAIM) leads to an exact and chemically intuitive energy partitioning scheme, the interacting quantum atoms (IQA) approach, that depends on the availability of second-order reduced density matrices (2-RDMs). This work explores the performance of this approach in particular and of the QTAIM in general with approximate 2-RDMs obtained from the density matrix functional theory (DMFT), which rests on the natural expansion (natural orbitals and their corresponding occupation numbers) of the first-order reduced density matrix (1-RDM). A number of these functionals have been implemented in the promolden code and used to perform QTAIM and IQA analyses on several representative molecules and model chemical reactions. Total energies, covalent intra- and interbasin exchange-correlation interactions, as well as localization and delocalization indices have been determined with these functionals from 1-RDMs obtained at different levels of theory. Results are compared to the values computed from the exact 2-RDMs, whenever possible.     

The two-component quantum theory of atoms in molecules (TC-QTAIM): foundations

The foundations of the two-component quantum theory of atoms in molecules (TC-QTAIM) are addressed in this contribution. In this regard, the theory is presented in an axiomatic manner and the main theorems describing regional properties of atoms in molecules are considered in detail. This is an extension of the orthodox quantum theory of atoms in molecules (QTAIM) for dealing with non-adiabatic wavefunctions of usual molecules as well as extracting the regional quantum structure of exotic species from the corresponding wavefunctions. The best examples of the latter are positronic and muonic species. The computational study of a model system consisting of a clamped lithium nucleus, four electrons, and a positively charged quantum particle carrying a unit of positive charge with a variable mass, m = 200–10¹³ m _e, supplements the theoretical argument demonstrating unambiguously that the TC-QTAIM analysis yields reasonable results. It reveals that the contribution of the positively charged particle in the topological analysis and basin properties is non-negligible. Most importantly, it is demonstrated that by increasing the mass of the positive particle, the TC-QTAIM analysis tends toward the QTAIM analysis of the lithium hydride system considered within the clamped nucleus paradigm. This result seems to indicate that the orthodox QTAIM is just the asymptote of the TC-QTAIM, the latter encompasses the former. Thus, one may claim that the TC-QTAIM is a unified framework for the AIM analysis of vast variety of quantum systems.     

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: 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     

A quantum theory atoms in molecules investigation of Lewis base protonation

This investigation uses atomic properties derived from the quantum theory of atoms in molecules formalism to rationalize the infrared intensity of the stretching vibration that arises as a Lewis base (B) is protonated (B‐H mode). Moreover, the interacting quantum atom (IQA) partition is employed to evaluate the energetics of protonation. All calculations are performed at the CCSD/cc‐pVQZ level except by the IQA analysis, which is carried out by means of the B3LYP/cc‐pVQZ//CCSD/cc‐pVQZ treatment. First, an efficiency scale is established for Lewis bases in terms of the electronic charge transfer potential. Next, this study shows that the intensity of the B‐H stretching depends mostly on the electronic charge amount transferred to the proton. Thus, intensity data provide empirical assessment of Lewis base charge transfer efficiency. Finally, the group separation observed during correlation of proton affinities and electronic charge transfer potential is explained by the interaction energy between fragments of the protonated system.     

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     

Stress tensor eigenvector following with next-generation quantum theory of atoms in molecules

The eigenvectors of the electronic stress tensor have been identified as useful for the prediction of chemical reactivity because they determine the most preferred directions to move the bonds. A new 3–D vector based interpretation of the chemical bond that we refer to as the bond-path framework set B provides a version of the quantum theory of atoms in molecules (QTAIM) beyond the minimum definition for bonding that is particularly suitable for understanding changes in molecular electronic structure that occur during reactions. We demonstrate that the most preferred direction for bond motion using the stress tensor corresponds to the most compressible direction and not to the least compressible direction as previously reported. We show the necessity for a directional approach constructed using the eigenvectors along the entire bond-length and demonstrate the insufficiency of scalar measures for capturing the nature of the stress tensor within the QTAIM partitioning.     

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.     

Red- and blue-shifted hydrogen bonds: the Bent rule from quantum theory of atoms in molecules perspective

MP2/6-311++G(3pd,3df) calculations were performed on complexes of acetylene and fluoroform acting as the proton donating systems and different Lewis bases being the proton acceptors since these complexes are linked through C-H···Y hydrogen bonds. Quantum Theory of Atoms in Molecules (QTAIM) is applied to explain the nature of these interactions. The characteristics of bond critical points are presented for these complexes. The inter-relations between energetic and geometrical parameters as well as the parameters derived from the Natural Bond Orbital (NBO) theory are analyzed here. Red- and blue-shifted hydrogen bonds are detected for the complexes investigated and the differences between those interactions are analyzed from the QTAIM perspective. It is shown that such differences are in agreement with the Bent rule. The position of the bond critical point of the proton donating bond is connected with the nature of hydrogen bonding, that is, if it is blue- or red-shifted.     

Next-generation quantum theory of atoms in molecules for the photochemical ring-opening reactions of oxirane

The conical intersections corresponding to the C─O and C─C ring opening were optimized and the reaction paths traversing these intersections were obtained. Investigation of the C─O ring opening revealed that when traversing the lowest energy conical intersection, the reaction path returns to the closed ring geometry. The C─O path traversing the intersection featuring torsion of terminal CH₂ group however, led to a ring-opened geometry, an H-shift and the formation of acetaldehyde that can undergo further dissociation. The observation of different reaction paths was explained by the 3-D paths from quantum theory of atoms in molecules (QTAIM) that defined the most preferred direction of electronic motion that precisely tracked the mechanisms of bond breaking and formation throughout the photo-reactions. The size, orientation, and location of these most preferred 3-D paths indicated the extent and direction of motion of atoms, bonds, and the degree of torsion or planarity of a bond indicating a predictive ability.     

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     

Simulant molecules with trivalent or pentavalent phosphorus atoms: bond dissociation energies and other thermodynamic and structural properties from quantum chemical models

The CBS-QB3 and G4 thermochemical models have been used to generate energetic, structural, and spectroscopic data on a set of molecules with trivalent or pentavalent phosphorus atoms that can serve as simulants of chemical warfare agents. Based on structural data, the conformational stabilities of these molecules are explained in terms of the anomeric interaction within the OPOC and OPSC fragments. For those cases where experimental data are available, comparisons have been made between calculated and previously reported vibrational frequencies. All varieties of bond dissociation energies have been examined except those for C-H and P═O bonds. In trivalent phosphorus molecules, the O-C and S-C bonds have the lowest dissociation energies. In the pentavalent phosphorus set, the S-C bonds, followed by P-S bonds, have the lowest dissociation energies. In the fluorinated simulant molecules, the P-F bond is strongest, and the P-C or O-C bonds are weakest.     

Interpretation of anomeric effect in the N-C-N unit with the quantum theory of atoms in molecules

The conformational preferences of six model compounds for the N-C-N anomeric unit (methanediamine, 2,2-propanediamine, N,N,N',N'-tetramethyl-methanediamine, 1,3-dizacyclohexane, 1,3,5-triazacyclohexane, and 2-aminopiperidine) were analyzed within the framework of the Quantum Theory of Atoms in Molecules. The relative stabilization of the conformers is related to two factors: (i) the reduction of the electron population experienced by the hydrogens of the central methylene when they display more gauche arrangements to lone pairs (lp) and (ii) the reduction of the electron population of aminic hydrogens when the corresponding N-H bond is in a parallel arrangement to the lone pair of another N. The former depletion takes place in lp-N-C-N antiperiplanar dispositions, whereas the latter is shown in lp-N-C-N gauche arrangements. Therefore, we can say that the electron density removed from the central hydrogens is moved to an aminic one on going from an antiperiplanar to a gauche disposition of a lp-N-C-N unit. The relative energies of aminic and central hydrogens in the conformer series is the main factor determining the conformational preference. In contrast to what happens in O-C-O containing compounds (where both N(H) depletions take place in the O-C-O-H gauche dispositions), the stabilization gained by N and C atoms plays a secondary role. This is in line with a general trend exhibited by hydrogens as the most available (less energy cost) atomic basins for receiving or providing electron density along a chemical change. It also explains why the anomeric conformational stabilization due to the N-C-N units is significantly less than that of the O-C-O- units. Moreover, the variations of electron population due to conformational changes are not in keeping with the stereoelectronic model of the anomeric effect, as was previously found for diverse molecules containing the O-C-O anomeric unit.