Long range behavior of high-rank topological multipole moments

The construction of a high-rank multipolar force field (for peptides) is a complex task, leading to several intermediate questions in need of a clear answer. Here we focus on the convergence of the (electrostatic) multipolar expansion at medium and long range. Using molecular electron densities, quantum chemical topology (QCT) defines the atoms as finite volumes, each endowed with multipole moments. The terms in the multipole expansion are grouped according to powers of the internuclear distance, R(-L). Given two atom types at a given distance, we determine which rank (L) is necessary for the electrostatic energy to converge to the exact interaction energy within a certain error. With this information, the rank of the expansion for each interaction can be adapted to the required accuracy and the available computing power.     

Realistic sampling of amino acid geometries for a multipolar polarizable force field

The Quantum Chemical Topological Force Field (QCTFF) uses the machine learning method kriging to map atomic multipole moments to the coordinates of all atoms in the molecular system. It is important that kriging operates on relevant and realistic training sets of molecular geometries. Therefore, we sampled single amino acid geometries directly from protein crystal structures stored in the Protein Databank (PDB). This sampling enhances the conformational realism (in terms of dihedral angles) of the training geometries. However, these geometries can be fraught with inaccurate bond lengths and valence angles due to artefacts of the refinement process of the X‐ray diffraction patterns, combined with experimentally invisible hydrogen atoms. This is why we developed a hybrid PDB/nonstationary normal modes (NM) sampling approach called PDB/NM. This method is superior over standard NM sampling, which captures only geometries optimized from the stationary points of single amino acids in the gas phase. Indeed, PDB/NM combines the sampling of relevant dihedral angles with chemically correct local geometries. Geometries sampled using PDB/NM were used to build kriging models for alanine and lysine, and their prediction accuracy was compared to models built from geometries sampled from three other sampling approaches. Bond length variation, as opposed to variation in dihedral angles, puts pressure on prediction accuracy, potentially lowering it. Hence, the larger coverage of dihedral angles of the PDB/NM method does not deteriorate the predictive accuracy of kriging models, compared to the NM sampling around local energetic minima used so far in the development of QCTFF. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

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.     

Incorporation of local structure into kriging models for the prediction of atomistic properties in the water decamer

Machine learning algorithms have been demonstrated to predict atomistic properties approaching the accuracy of quantum chemical calculations at significantly less computational cost. Difficulties arise, however, when attempting to apply these techniques to large systems, or systems possessing excessive conformational freedom. In this article, the machine learning method kriging is applied to predict both the intra‐atomic and interatomic energies, as well as the electrostatic multipole moments, of the atoms of a water molecule at the center of a 10 water molecule (decamer) cluster. Unlike previous work, where the properties of small water clusters were predicted using a molecular local frame, and where training set inputs (features) were based on atomic index, a variety of feature definitions and coordinate frames are considered here to increase prediction accuracy. It is shown that, for a water molecule at the center of a decamer, no single method of defining features or coordinate schemes is optimal for every property. However, explicitly accounting for the structure of the first solvation shell in the definition of the features of the kriging training set, and centring the coordinate frame on the atom‐of‐interest will, in general, return better predictions than models that apply the standard methods of feature definition, or a molecular coordinate frame. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

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.     

The topology of the Ehrenfest force density revisited. A different perspective based on Slater‐type orbitals

The topology of the Ehrenfest force density was studied with Slater‐type orbitals (STO). At larger distances from the nuclei, STOs generate similar artefacts as noticed before with Gaussian‐type orbitals. The topology of the Ehrenfest force density was found to be mainly homeomorphic with the topology of the electron density. For the first time, reliable integrations of several properties over force density atomic basins were performed successfully. Integration of the electron density of a number of hydrides, fluorides, and chlorides of first row elements over force density basins indicate substantial differences between the partial charges of the atoms as compared with those obtained from electron density basins. Calculations on saturated hydrocarbons confirm that the electronegativity of carbon atoms increases with increasing geometrical strain. Atomic interaction lines are observed to exist in the Ehrenfest force density between the hydrogen atoms of several so‐called “congested” molecules, and also in some inclusion complexes of alkanes with helium. However, interaction lines are lacking in several other controversial cases. © 2015 Wiley Periodicals, Inc.     

Point atomic multipole moments for simulation of electrostatic potential and field in all-siliceous zeolites

Calibration method of atomic multipole moments (AMMs) is presented with respect to geometries of all-siliceous zeolite models obtained with X-ray diffraction (XRD) methods. Mulliken atomic charges and AMMs are calculated for all-siliceous types possessing small size elementary unit cells at the hybrid density functional theory (DFT) (B3LYP) and general gradient approximation (GGA) Perdew-Burke-Ernzerhof (PBE) levels and then used to fit the dependences versus geometry variables for the Mulliken charges and versus special coordinate for the AMMs. Fitted and exact charges and AMMs are used to compute electrostatic potential (EP) and electric field (EF) for all-siliceous zeolites with CRYSTAL. A possibility of application of the point AMMs to quantum mechanical/molecular mechanics computations or classic simulation of physical adsorption is evaluated. The considered models expand over wide range of structural parameters and could be applied even to amorphous all-siliceous systems.     

Evaluation of molecular moments by three methods

The three evaluative methods of the moment by means of the coefficient a_k of the HMO characteristic polynomial, the tree graph, and walk graph have been studied. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 237–244, 2000     

Molecular QTAIM Topology Is Sensitive to Relativistic Corrections

The topology of the molecular electron density of benzene dithiol gold cluster complex Au₄−S−C₆H₄−S′−Au′₄ changed when relativistic corrections were made and the structure was close to a minimum of the Born–Oppenheimer energy surface. Specifically, new bond paths between hydrogen atoms on the benzene ring and gold atoms appeared, indicating that there is a favorable interaction between these atoms at the relativistic level. This is consistent with the observation that gold becomes a better electron acceptor when relativistic corrections are applied. In addition to relativistic effects, here, we establish the sensitivity of molecular topology to basis sets and convergence thresholds for geometry optimization.     

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.     

Accuracy and tractability of a kriging model of intramolecular polarizable multipolar electrostatics and its application to histidine

We propose a generic method to model polarization in the context of high‐rank multipolar electrostatics. This method involves the machine learning technique kriging, here used to capture the response of an atomic multipole moment of a given atom to a change in the positions of the atoms surrounding this atom. The atoms are malleable boxes with sharp boundaries, they do not overlap and exhaust space. The method is applied to histidine where it is able to predict atomic multipole moments (up to hexadecapole) for unseen configurations, after training on 600 geometries distorted using normal modes of each of its 24 local energy minima at B3LYP/apc‐1 level. The quality of the predictions is assessed by calculating the Coulomb energy between an atom for which the moments have been predicted and the surrounding atoms (having exact moments). Only interactions between atoms separated by three or more bonds (“1, 4 and higher” interactions) are included in this energy error. This energy is compared with that of a central atom with exact multipole moments interacting with the same environment. The resulting energy discrepancies are summed for 328 atom–atom interactions, for each of the 29 atoms of histidine being a central atom in turn. For 80% of the 539 test configurations (outside the training set), this summed energy deviates by less than 1 kcal mol^?1. © 2013 Wiley Periodicals, Inc.     

Testing theory beyond molecular structure: Electron density distributions of complex molecules

The electron density distribution (EDD) of a molecular system can be determined experimentally from elaborate X‐ray diffraction measurements or calculated with quantum mechanical methods: This provides a unique opportunity for mutual validation of the experimental and theoretical methods—a validation that goes far beyond comparison of molecular structures. Two examples of complex molecular systems of biologic relevance are presented. The first is the cocrystallized complex of betaine, imidazole, and picric acid, 1, which is a 75‐atom molecular complex serving as a model for the active site in the serine proteases class of enzymes, the so‐called catalytic triad. For 1 the experimental charge density was determined by combined modeling of single crystal synchrotron X‐ray and neutron diffraction data measured at 28(1) K, and it is compared with ab initio theoretical calculations at the B3LYP/6‐311G(d,p) level of theory. Overall, the agreement is good, but in one strong N? H? O hydrogen bond clear differences are observed. The second example concerns the EDD of the mixed valence trinuclear oxo‐centered iron carboxylate, [Fe₃O(OOCC(CH₃)₃)₆(NC₅H₅)₃], 2. This molecule contains 133 atoms (542 electrons) including three open‐shell iron atoms, and the experimental investigation is based on synchrotron X‐ray diffraction data. Calculations in the experimental geometry at the commonly used UB3LYP/LanL2DZ level of theory are not able to reproduce a number of experimentally observed electron density features. In particular, the sp³‐like distribution on the central oxygen atom and the electron deformations on the iron centers are at variance with experiment. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

The localized electrons detector as an ab initio representation of molecular structures

The local single particle momentum is proposed as a localized‐electrons detector (LED) that provides a direct three‐dimensional representation of bonding interactions in molecules. It is given exclusively in terms of the electron density and its gradient. We show that the graphical representation of bonding interactions given by LED is consistent with the local curvatures of the electron density as given by the eigenvalues of the Hessian matrix, according to a local symmetry classification of the critical points here introduced. LED consistently complements the topological analysis of the electron density given by the quantum theory of atoms in molecules, by providing a graphical representation of the symmetry of the bonding interactions in molecular systems. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2418–2425, 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     

Atomic charges,dipole moments,and Fukui functions using the Hirshfeld partitioning of the electron density

In the Hirshfeld partitioning of the electron density, the molecular electron density is decomposed in atomic contributions, proportional to the weight of the isolated atom density in the promolecule density, constructed by superimposing the isolated atom electron densities placed on the positions the atoms have in the molecule. A maximal conservation of the information of the isolated atoms in the atoms-in-molecules is thereby secured. Atomic charges, atomic dipole moments, and Fukui functions resulting from the Hirshfeld partitioning of the electron density are computed for a large series of molecules. In a representative set of organic and hypervalent molecules, they are compared with other commonly used population analysis methods. The expected bond polarities are recovered, but the charges are much smaller compared to other methods. Condensed Fukui functions for a large number of molecules, undergoing an electrophilic or a nucleophilic attack, are computed and compared with the HOMO and LUMO densities, integrated over the Hirshfeld atoms in molecules.     

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.