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

Latin American contributions to quantum chemical topology

We survey the contributions from Latin American theoretical chemists to the field of quantum chemical topology (QCT) over nearly the last 30 years with emphasis on the developments and applications of the quantum theory of atoms in molecules (QTAIM). Applications of QCT in the fields of excited states, electron delocalization, chemical bond, aromaticity, conformational analysis, spectroscopic properties, and chemical reactivity are presented. We also consider the coupling of QTAIM with time-dependent density functional theory, the virial theorem in the Kohn-Sham method and the inclusion of electron dynamical correlation in the interacting quantum atoms method using coupled cluster and multi-configurational densities. Additionally, we describe the development of efficient algorithms for the calculation of topological properties derived from the electron density. This review is aimed not only at providing an account of the contributions to QCT in Latin America but also at stimulating guides for further progress in the field.     

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.     

Commentaries on quantum similarity (1): Density gradient quantum similarity

Computation of density gradient quantum similarity integrals is analyzed, while comparing such integrals with overlap density quantum similarity measures. Gradient quantum similarity corresponds to another kind of numerical similarity assessment between a pair of molecular frames, which contrarily to the usual up to date quantum similarity definitions are not measures, that is: strictly positive definite integrals. As the density gradient quantum similarity integrals are defined as scalar products of three real functions, they appear to possess a richer structure than the corresponding positive definite density overlap quantum similarity measures, while preserving the overall similarity trends, when the molecular frames are relatively moved in three‐dimensional space. Similarity indices are also studied when simple cases are analyzed in order to perform more comparisons with density overlap quantum similarity. Multiple gradient quantum similarity integrals are also defined. General GTO formulae are given. Numerical results within the atomic shell approximation (ASA) framework are presented as simple examples showing the new performances of the gradient density quantum similarity. Fortran 90 programs illustrating the proposed theoretical development can be downloaded from appropriate websites. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

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.     

Rhorix: An interface between quantum chemical topology and the 3D graphics program blender

Chemical research is assisted by the creation of visual representations that map concepts (such as atoms and bonds) to 3D objects. These concepts are rooted in chemical theory that predates routine solution of the Schrödinger equation for systems of interesting size. The method of Quantum Chemical Topology (QCT) provides an alternative, parameter‐free means to understand chemical phenomena directly from quantum mechanical principles. Representation of the topological elements of QCT has lagged behind the best tools available. Here, we describe a general abstraction (and corresponding file format) that permits the definition of mappings between topological objects and their 3D representations. Possible mappings are discussed and a canonical example is suggested, which has been implemented as a Python “Add‐On” named Rhorix for the state‐of‐the‐art 3D modeling program Blender. This allows chemists to use modern drawing tools and artists to access QCT data in a familiar context. A number of examples are discussed. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Theory of time‐dependent transport in quantum dot systems

The formulation of the time‐dependent tunnel current in quantum dot (QD) systems in the (generalized) transfer Hamiltonian formalism is reconsidered, taking into account the nonorthogonality between the subsystems. Exploiting the fact that only the total charge in the system is conserved, in general, gives rise to a new formulation of the transport theory. As a result, it is shown necessary to depart from the orthodox picture, in which the current is treated as a local property of the system. A general formula for the current is derived. By expressing properties of the quantum dot in the Hubbard I approximation with the loop correction, thereby including strong correlations of the QD states, numerical results are given for various time‐dependent source‐drain voltages. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Identification of hydrogen bonds using quantum electrodynamics

A method for the identification of hydrogen bonds was investigated from the viewpoint of the stress tensor density proposed by Tachibana and following other works in this field. Hydrogen bonds are known to exhibit common features with ionic and covalent bonds. In quantum electrodynamics, the covalent bond has been demonstrated to display a spindle structure of the stress tensor density. Importantly, this spindle structure is also seen in the hydrogen bond, although the covalency is considerably weaker than in a typical covalent bond. Distinguishing it from the ionic bond is most imperative for the identification of the hydrogen bond. In the present study, the directionality of the hydrogen bond is investigated as the ionic bond is nearly isotropic, while the hydrogen bond exhibits the directionality. It was demonstrated that the hydrogen bond can be distinguished from the ionic bond using the angle dependence of the largest eigenvalue of the stress tensor density.     

Generalized formulation of quantum and classical crystallography using Green's functions

Quantum crystallography is a methodology by which structural information about a crystalline material obtained from X‐ray crystallography is combined with quantum mechanical methods. The objective is to enhance the data obtained from the X‐ray diffraction experiment, which are related to the atomic structure of the crystal, and to predict the properties and efficacy of those chemical compounds from which the crystals are derived. One approach in quantum crystallography is to use a projector matrix with a normalized trace. In this approach, quantum mechanical parameters in the projector matrix are fit into crystallographic data. During this fitting, the properties of the projector matrix called idempotency and normalization are used. Throughout this implementation procedure, Clinton's iteration scheme has been used in addition to the least‐squares technique. The purpose of the present study is to generalize Clinton's iterative equations in quantum crystallography by means of single‐particle Green's functions with the aid of the equal atoms model in the theory of direct methods. Convergency characters of the novel iterative equations are discussed by the steepest descent procedure. Furthermore, whether the calculations are valid in nonorthogonal bases was also examined. The iteration schemes widely used in quantum crystallography have been generalized but, in addition, the generalized expressions relating to the phase determination procedure and the probabilities of the sign relations between the structure factors are obtained and discussed comprehensively. The phrase order of crystallography has been put forward as a new concept. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

An attempt at quantum thermal physics

Summary This paper outlines an alternative exposition of the structure of quantum thermodynamics which is essentially based on Carnot's theory where fluxes of caloric are identified with negative information fluxes. It is further assumed that the thermal energy evolved by thermal processes is identical with the electromagnetic zero-point background energy evolved by the destruction of information inscribed in a structural unit (qubit). Theoretical arguments on an elementary level are accompanied by illustrative examples.     

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     

Density dynamics in some quantum systems

The quantum domain behavior of classical nonintegrable systems is well‐understood by the implementation of quantum fluid dynamics and quantum theory of motion. These approaches properly explain the quantum analogs of the classical Kolmogorov–Arnold–Moser type transitions from regular to chaotic domain in different anharmonic oscillators. Field‐induced tunneling and chaotic ionization in Rydberg atoms are also analyzed with the help of these theories. Quantum fluid density functional theory may be used to understand different time‐dependent processes like ion‐atom/molecule collisions, atom‐field interactions, and so forth. Regioselectivity as well as confined atomic/molecular systems and their reactivity dynamics have also been explained. © 2013 Wiley Periodicals, Inc.     

Visualization and integration of quantum topological atoms by spatial discretization into finite elements

We present a novel algorithm to integrate property densities over the volume of a quantum topological atom. Atoms are grown outward, starting from a sphere centered on the nucleus, by means of a finite element meshing algorithm. Bond critical points and ring critical points require special treatment. The overall philosophy as well as intricate features of this meshing algorithm are given, followed by details of the quadrature over the finite elements. An effort has been made to design a streamlined and compact algorithm, focusing on the core of challenges arising in tracing the electron density's gradient vector field. The current algorithm also generates a new type of pictures that can be a Graphical User Interface. Excellent integration errors, L(Omega), are obtained, even for atoms with (narrow) tails or sharp corners.     

Calculation of wave‐functions with frozen orbitals in mixed quantum mechanics/molecular mechanics methods. II. Application of the local basis equation

The application of the local basis equation (Ferenczy and Adams, J. Chem. Phys. 2009 , 130, 134108) in mixed quantum mechanics/molecular mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) methods is investigated. This equation is suitable to derive local basis nonorthogonal orbitals that minimize the energy of the system and it exhibits good convergence properties in a self‐consistent field solution. These features make the equation appropriate to be used in mixed QM/MM and QM/QM methods to optimize orbitals in the field of frozen localized orbitals connecting the subsystems. Calculations performed for several properties in divers systems show that the method is robust with various choices of the frozen orbitals and frontier atom properties. With appropriate basis set assignment, it gives results equivalent with those of a related approach [G. G. Ferenczy previous paper in this issue] using the Huzinaga equation. Thus, the local basis equation can be used in mixed QM/MM methods with small size quantum subsystems to calculate properties in good agreement with reference Hartree–Fock–Roothaan results. It is shown that bond charges are not necessary when the local basis equation is applied, although they are required for the self‐consistent field solution of the Huzinaga equation based method. Conversely, the deformation of the wave‐function near to the boundary is observed without bond charges and this has a significant effect on deprotonation energies but a less pronounced effect when the total charge of the system is conserved. The local basis equation can also be used to define a two layer quantum system with nonorthogonal localized orbitals surrounding the central delocalized quantum subsystem. © 2013 Wiley Periodicals, Inc.     

Local physical quantities for spin based on the relativistic quantum field theory in molecular systems

Local physical quantities for spin are investigated on the basis of the four‐ and two‐component relativistic quantum theory. In the quantum field theory, local physical quantities for spin such as the spin angular momentum density, spin torque density, zeta force density, and zeta potential play important roles in spin dynamics. We discuss how to calculate these local physical quantities based on the two‐component relativistic quantum theory. Some different types of relativistic numerical calculations of local physical quantities in Li atom and C₆H₆ are demonstrated and compared. Local physical quantities for each orbital are also discussed, and it is seen that a total local zeta potential is given as a result of some cancellation of large contributions from each orbital. © 2016 Wiley Periodicals, Inc.