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.     

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.     

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.     

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.     

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.     

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     

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.     

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.     

Congested molecules. Where is the steric repulsion? An analysis of the electron density by the method of interacting quantum atoms

The computed electron density of several congested saturated hydrocarbons and halogenated derivatives has been analyzed by the method of interacting quantum atoms (IQA). For all the molecules studied, the calculations show the existence of a bond path between the congested atoms and which, according to the Quantum Theory of Atoms in Molecules, indicates that there is a stabilizing interaction between these atoms. The bond path is found to exist up to interatomic distances well‐beyond the sum of the van der Waals radii. The IQA results indicate that steric hindrance is not a repulsive force between the congested atoms but that is the result of an increase in the intra‐atomic or self‐energy of the congested atoms. This increase in self‐energy is caused by the deformation of the atomic basin of the congested atoms. © 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.     

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.     

The effects of higher orders of perturbation theory on the correlation energy of atoms and bonds in molecules

We examine, for the first time, the effects of higher orders of Møller–Plesset perturbation theory on the individual atoms within a molecule and the bonds between them, via the topological energy partitioning method of interacting quantum atoms. In real terms (i.e., not by absolute value) MP3 decreases the correlation energy of a bond, and MP4SDQ also decreases the energy of the atoms at either end of the bond. In addition, we investigated long‐range through‐space dispersive effects on a H₂ oligomer. Overall, MP3 is the largest correction to the correlation energy, and most of that energy is allocated to chemical bonds, reducing their values in actual terms. The MP4SDQ bond correlation correction, despite being relatively small, tends to have two effects: (i) for small or negative correlation energies MP4SDQ tends to decrease the bond correlation values even more, and (ii) for large (positive) bond correlation energies MP4SDQ tends to restore the bond correlation energies from the MP3 back toward the MP2 values. Furthermore, each individual part of a molecule or complex (atom or bond) has a specific convergence pattern for the MPn series: through‐space interactions converge at MP2 but bonds converge at MP3 level. The atomic correlation energy appears to head toward convergence at the MP4 level.     

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     

Prediction of conformationally dependent atomic multipole moments in carbohydrates

The conformational flexibility of carbohydrates is challenging within the field of computational chemistry. This flexibility causes the electron density to change, which leads to fluctuating atomic multipole moments. Quantum Chemical Topology (QCT) allows for the partitioning of an “atom in a molecule,” thus localizing electron density to finite atomic domains, which permits the unambiguous evaluation of atomic multipole moments. By selecting an ensemble of physically realistic conformers of a chemical system, one evaluates the various multipole moments at defined points in configuration space. The subsequent implementation of the machine learning method kriging delivers the evaluation of an analytical function, which smoothly interpolates between these points. This allows for the prediction of atomic multipole moments at new points in conformational space, not trained for but within prediction range. In this work, we demonstrate that the carbohydrates erythrose and threose are amenable to the above methodology. We investigate how kriging models respond when the training ensemble incorporating multiple energy minima and their environment in conformational space. Additionally, we evaluate the gains in predictive capacity of our models as the size of the training ensemble increases. We believe this approach to be entirely novel within the field of carbohydrates. For a modest training set size of 600, more than 90% of the external test configurations have an error in the total (predicted) electrostatic energy (relative to ab initio) of maximum 1 kJ mol^?1 for open chains and just over 90% an error of maximum 4 kJ mol^?1 for rings. © 2015 Wiley Periodicals, Inc.