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     

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.     

Quantum topology phase diagrams for the cis‐ and trans‐isomers of the cyclic contryphan‐Sm peptide

Within the quantum theory of atoms in molecules (QTAIM) framework we present a quantum topology phase diagram (QTPD) using the Poincaré–Hopf relation of a total of 17 all new QTAIM topologies of the cis‐ and trans‐isomers of the cyclic contryphan‐Sm peptide. The resultant QTPD consists of separate regions for the cis‐ and trans‐isomers that only overlap for topologies associated with the lowest energy minima of the cis‐ and trans‐isomers. We determine the QTAIM topologies of 29 “missing” isomers. A new, contracted formulation of the QTPD is presented, this contracted formulation includes the interamino acid bond critical points (BCPs) that link together the amino acid units, the disulphide bridge “pivot” BCP and side chain bonding interactions. The seven interamino acid BCPs linking the amino acid units coincide with the so‐called peptide backbone, the conventional qualitative approach to reduce the complexity of the peptide. We expand the interpretation of ellipticity to include the associated eigenvectors and find that higher values of the ellipticity ? are associated with a greater preference to conserve folding states. We quantify previous qualitative findings that suggested the disulfide bond is central to the folding behavior of the cyclic contryphan‐Sm peptide and why the cis‐isomer is the major form of the cyclic contryphan‐Sm peptide. © 2014 Wiley Periodicals, Inc.     

Revisiting the foundations of the quantum theory of atoms in molecules: The subsystem variational procedure and the finite nuclear models

The role of finite nuclear models (FNMs) is scrutinized within the context of the quantum theory of atoms in molecules (QTAIMs). It is demonstrated that the newly proposed analytic‐algebraic definition of the topological atoms is consistently extendable to the cases where a FNM is employed to construct the molecular hamiltonian. The whole variational procedure is reconsidered, and the insensitivity of final results relative to the employed FNMs is explicitly demonstrated. The analysis once again clearly demonstrates that the analytic‐algebraic condition is an independent axiom that must be added to the subsystem variational procedure to construct the QTAIMs. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

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     

Chemical bonding in view of electron charge density and kinetic energy density descriptors

Stalke's dilemma, stating that different chemical interpretations are obtained when one and the same density is interpreted either by means of natural bond orbital (NBO) and subsequent natural resonance theory (NRT) application or by the quantum theory of atoms in molecules (QTAIM), is reinvestigated. It is shown that within the framework of QTAIM, the question as to whether for a given molecule two atoms are bonded or not is only meaningful in the context of a well‐defined reference geometry. The localized‐orbital‐locator (LOL) is applied to map out patterns in covalent bonding interaction, and produces results that are consistent for a variety of reference geometries. Furthermore, LOL interpretations are in accord with NBO/NRT, and assist in an interpretation in terms of covalent bonding. © 2008 Wiley Periodicals, Inc.J Comput Chem, 2009.     

Extending the Domain‐Averaged Exchange‐Correlation Energies Within the Context of the MC‐QTAIM: Tracing Subtle Variations Induced by Isotope Substitution

Recently, it has been demonstrated that the domain‐averaged exchange‐correlation energies, V^xc, are capable of tracing the covalent character of atom–atom interactions unequivocally and thus pave the way for detailed bonding analysis within the context of the quantum theory of atoms in molecules (QTAIM) [M. García‐Revilla, E. Francisco, P. L. Popelier, A. Martín Pendás, ChemPhysChem 2013 , 14, 1211–1218]. Herein, the concept of V^xc is extended within the context of the newly developed multicomponent QTAIM (MC‐QTAIM). The extended version, , is capable of analyzing nonadiabatic wavefunctions and thus is sensitive to the mass of nuclei and can trace “locally” the subtle electronic variations induced by isotope substitution. To demonstrate this capability in practice, ab initio nonadiabatic wavefunctions for three isotopically substituted hydrogen cyanide molecules, in which the hydrogen nucleus was assumed to be a proton, deuterium, or tritium, were derived. The resulting wavefunctions were then used to compute and it emerged that for the hydrogen–carbon bond, the was distinct for each isotopic composition and varied in line with chemical expectations. Indeed, the introduction of paves the way for the investigation of vast numbers of structural and kinetic isotope effects within the context of the MC‐QTAIM.     

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.     

Toward a fuzzy atom view within the context of the quantum theory of atoms in molecules: quasi-atoms

The notion of quasi-atoms is introduced within the context of the quantum theory of atoms in molecules. Being a subset of the quantum divided basins that were introduced previously, quasi-atoms are the quantum subsystems that are practically indistinguishable from the topological atoms; thus, revealing the continuous evolution of quantum divided basins into topological atoms. This indistinguishability is rooted in the limited accuracy of chemical observations; they are not sensitive to discriminate a topological atom from its associated quasi-atoms. In this regard, it is disclosed that the set of quantum atoms is in a wide-range including members other than topological atoms; the quasi-atoms are concrete examples. Finally, the idea of the fuzzy set of atoms that is foreign to the disjoint partitioning schemes for which the orthodox QTAIM is a classic example is extended employing the set of quasi-atoms.     

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.     

On the connections between the quantum theory of atoms in molecules (QTAIM) and density functional theory (DFT): a letter from Richard F. W. Bader to Lou Massa

A 31-year-old letter from Professor Richard F. W. Bader to Professor Lou Massa outlining the connections between the quantum theory of atoms in molecules (QTAIM) and density functional theory (DFT) especially with regard to the first Hohenberg-Kohn theorem is brought to light. This connection has not often been the topic of such a focused review by Bader and is presented here for the first time. The scientific importance of this letter is, in the opinion of the presenter, as timely today as it was back then in 1986. In Bader’s own opening words: “... that if I sent you a summary of what I think are the important connections between our work and density functional theory, ...”. He then takes us in a grand tour of the foundations of QTAIM culminating into the antecedents of a paper he later published with Professor Pierre Becker, whereby the Hohenberg-Kohn theorem is shown to operate at the level of an atom-in-a-molecule. Bader closes his letter by suggesting to Massa: “Study these two charge distributions – they are proof of the theorem of Hohenberg and Kohn”. By that Bader meant that when the charge distributions of two atoms or groups are identical within a given precision, then the kinetic and total energy contributions of these atoms to the corresponding molecular quantities are also identical. It is revealing to follow the intellectual threads weaved by Bader which provides us with a glimpse of his thought processes and intuition that guided him to some of his key discoveries. The lucidity, rigor, and clarity characteristic of Bader and the informality of style of a letter makes it of pedagogic and historic interest.

     

A vector‐based representation of the chemical bond for the normal modes of benzene

We introduce a vector‐based interpretation of the chemical bond within the quantum theory of atoms in molecules (QTAIM), the bond‐path framework set B = {p, q, r}, to follow variations in the 3D morphology of all bonds for the four infrared active normal modes of benzene. The bond‐path framework set B comprises three unique paths p, q, and r where r is the familiar QTAIM bond concept of bond‐path (r) while the two new paths p and q are formulated from the least and most preferred directions of electron density accumulation, respectively. We find 3D distortions including bond stretching/compression, torsion, and curving. We introduce two fractional measures to quantify these variations away from linearity of the bond.     

Toward a Consistent Interpretation of the QTAIM: Tortuous Link between Chemical Bonds,Interactions, and Bond/Line Paths

Currently, bonding analysis of molecules based on the Quantum Theory of Atoms in Molecules (QTAIM) is popular; however, “misinterpretations” of the QTAIM analysis are also very frequent. In this contribution the chemical relevance of the bond path as one of the key topological entities emerging from the QTAIM’s topological analysis of the one‐electron density is reconsidered. The role of nuclear vibrations on the topological analysis is investigated demonstrating that the bond paths are not indicators of chemical bonds. Also, it is argued that the detection of the bond paths is not necessary for the “interaction” to be present between two atoms in a molecule. The conceptual disentanglement of chemical bonds/interactions from the bonds paths, which are alternatively termed “line paths” in this contribution, dismisses many superficial inconsistencies. Such inconsistencies emerge from the presence/absence of the line paths in places of a molecule in which chemical intuition or alternative bonding analysis does not support the presence/absence of a chemical bond. Moreover, computational QTAIM studies have been performed on some “problematic” molecules, which were considered previously by other authors, and the role of nuclear vibrations on presence/absence of the line paths is studied demonstrating that a bonding pattern consistent with other theoretical schemes appears after a careful QTAIM analysis and a new “interpretation” of data is performed.     

The Role of Molecular Electrostatic Potentials in the Formation of a Halogen Bond in Furan⋅⋅⋅XY and Thiophene⋅⋅⋅XY Complexes

The halogen bonding of furan???XY and thiophene???XY (X=Cl, Br; Y=F, Cl, Br), involving σ‐ and π‐type interactions, was studied by using MP2 calculations and quantum theory of “atoms in molecules” (QTAIM) studies. The negative electrostatic potentials of furan and thiophene, as well as the most positive electrostatic potential (V_S,max) on the surface of the interacting X atom determined the geometries of the complexes. Linear relationships were found between interaction energy and V_S,max of the X atom, indicating that electrostatic interactions play an important role in these halogen‐bonding interactions. The halogen‐bonding interactions in furan???XY and thiophene???XY are weak, “closed‐shell” noncovalent interactions. The linear relationship of topological properties, energy properties, and the integration of interatomic surfaces versus V_S,max of atom X demonstrate the importance of the positive σ hole, as reflected by the computed V_S,max of atom X, in determining the topological properties of the halogen bonds.     

Fatigue and photochromism S1 excited state reactivity of diarylethenes from QTAIM and the stress tensor

The realization of technologically relevant functional systems from idealized photochromic compounds remains elusive due to the double requirement that such switches must possess both highly efficient photo‐isomerization reactivity and extremely low fatigue over a large number of switching cycles. Nowadays, improvements of the switching properties in complex diarylethene structures are mainly attained on a “trial and error” basis through chemical substitutions aimed at tuning the chemical properties of the core of the diarylethene. Therefore, we present new guiding principles to analyze the first excited state reactivity of diarylethenes based on the quantum theory of atoms in molecules (QTAIM) including the stress tensor. This approach straightforwardly provides consistent theoretical justification to partner the already successful symmetric substitution patterns obtained from experiments. The guiding principles provided by QTAIM and stress tensor suggest more complex asymmetric patterns should be included for the systematic design of new technologically relevant functional compounds. The stress tensor trajectory T_σ(s) analysis is used to characterize the photochromism reaction as reusable and the fatigue reaction as irreversible and find candidate sites for alteration by future experiment.     

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.     

Atoms in molecules: beyond Born–Oppenheimer paradigm

This contribution presents the first atoms in molecules study that goes beyond the Born–Oppenheimer paradigm employing the newly developed two-component quantum theory of atoms in molecules (TC-QTAIM). The LiH, LiD, and LiT systems containing quantum instead of clamped hydrogen nuclei are used as typical examples. The computational analysis that is done on non-adiabatic wavefunctions derived from the fully variational multicomponent molecular orbital approach (FV-MC-MO) results in hydrogen atomic basins without any clamped nucleus. The topological analysis of the Γ-field, the field that replaces the usual one-electron density used in the orthodox topological analysis, reveals delicate differences among the considered systems. The calculation of basin properties also demonstrates that the TC-QTAIM differentiates among atomic basins containing isotopes. Since the nuclear dynamics is contained intrinsically in non-adiabatic wavefunctions, the nuclear contribution to both topological analysis and basin properties naturally emerges from the TC-QTAIM analysis resolving the long-standing obstacle of consistent incorporation of nuclear dynamics within the context of the orthodox QTAIM. Also, a similar analysis is done on non-adiabatic wavefunctions describing excited instead of ground nuclear vibrations of the considered systems demonstrating the fact that TC-QTAIM is capable of being employed for both ground and excited nuclear vibrational states.