skip to main content


Title: Introducing DDEC6 atomic population analysis: part 3. Comprehensive method to compute bond orders
Developing a comprehensive method to compute bond orders is a problem that has eluded chemists since Lewis's pioneering work on chemical bonding a century ago. Here, a computationally efficient method solving this problem is introduced and demonstrated for diverse materials including elements from each chemical group and period. The method is applied to non-magnetic, collinear magnetic, and non-collinear magnetic materials with localized or delocalized bonding electrons. Examples studied include the stretched O 2 molecule, 26 diatomic molecules, 3d and 5d transition metal solids, periodic materials with 1 to 8748 atoms per unit cell, a biomolecule, a hypercoordinate molecule, an electron deficient molecule, hydrogen bound systems, transition states, Lewis acid–base complexes, aromatic compounds, magnetic systems, ionic materials, dispersion bound systems, nanostructures, and other materials. From near-zero to high-order bonds were studied. Both the bond orders and the sum of bond orders for each atom are accurate across various bonding types: metallic, covalent, polar-covalent, ionic, aromatic, dative, hypercoordinate, electron deficient multi-centered, agostic, and hydrogen bonding. The method yields similar results for correlated wavefunction and density functional theory inputs and for different S Z values of a spin multiplet. The method requires only the electron and spin magnetization density distributions as input and has a computational cost scaling linearly with increasing number of atoms in the unit cell. No prior approach is as general. The method does not apply to electrides, highly time-dependent states, some extremely high-energy excited states, and nuclear reactions.  more » « less
Award ID(s):
1555376
NSF-PAR ID:
10042797
Author(s) / Creator(s):
Date Published:
Journal Name:
RSC Adv.
Volume:
7
Issue:
72
ISSN:
2046-2069
Page Range / eLocation ID:
45552 to 45581
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Molecules with unstable isotopes often contain heavy and deformed nuclei and thus possess a high sensitivity to parity-violating effects, such as the Schiff moments. Currently the best limits on Schiff moments are set with diamagnetic atoms. Polar molecules with quantum-enhanced sensing capabilities, however, can offer better sensitivity. In this work, we consider the prototypical 223 Fr 107 Ag molecule, as the octupole deformation of the unstable 223 Fr francium nucleus amplifies the nuclear Schiff moment of the molecule by two orders of magnitude relative to that of spherical nuclei and as the silver atom has a large electron affinity. To develop a competitive experimental platform based on molecular quantum systems, 223 Fr atoms and 107 Ag atoms have to be brought together at ultracold temperatures. That is, we explore the prospects of forming 223 Fr 107 Ag from laser-cooled Fr and Ag atoms. We have performed fully relativistic electronic-structure calculations of ground and excited states of FrAg that account for the strong spin-dependent relativistic effects of Fr and the strong ionic bond to Ag. In addition, we predict the nearest-neighbor densities of magnetic-field Feshbach resonances in ultracold 223 Fr + 107 Ag collisions with coupled-channel calculations. These resonances can be used for magneto-association into ultracold, weakly-bound FrAg. We also determine the conditions for creating 223 Fr 107 Ag molecules in their absolute ground state from these weakly-bound dimers via stimulated Raman adiabatic passage using our calculations of the relativistic transition electric dipole moments. 
    more » « less
  2. The DDEC6 method is one of the most accurate and broadly applicable atomic population analysis methods. It works for a broad range of periodic and non-periodic materials with no magnetism, collinear magnetism, and non-collinear magnetism irrespective of the basis set type. First, we show DDEC6 charge partitioning to assign net atomic charges corresponds to solving a series of 14 Lagrangians in order. Then, we provide flow diagrams for overall DDEC6 analysis, spin partitioning, and bond order calculations. We wrote an OpenMP parallelized Fortran code to provide efficient computations. We show that by storing large arrays as shared variables in cache line friendly order, memory requirements are independent of the number of parallel computing cores and false sharing is minimized. We show that both total memory required and the computational time scale linearly with increasing numbers of atoms in the unit cell. Using the presently chosen uniform grids, computational times of ∼9 to 94 seconds per atom were required to perform DDEC6 analysis on a single computing core in an Intel Xeon E5 multi-processor unit. Parallelization efficiencies were usually >50% for computations performed on 2 to 16 cores of a cache coherent node. As examples we study a B-DNA decamer, nickel metal, supercells of hexagonal ice crystals, six X@C 60 endohedral fullerene complexes, a water dimer, a Mn 12 -acetate single molecule magnet exhibiting collinear magnetism, a Fe 4 O 12 N 4 C 40 H 52 single molecule magnet exhibiting non-collinear magnetism, and several spin states of an ozone molecule. Efficient parallel computation was achieved for systems containing as few as one and as many as >8000 atoms in a unit cell. We varied many calculation factors ( e.g. , grid spacing, code design, thread arrangement, etc. ) and report their effects on calculation speed and precision. We make recommendations for excellent performance. 
    more » « less
  3. null (Ed.)
    Spin crossover (SCO) is a phenomenon observed for certain transition metal complexes with electronic configuration 3d4-3d7. The conversion between the low-spin (LS) and high-spin (HS) states is usually driven by a variety of external perturbations, such as temperature, pressure, or light. The switching between the enthalpically preferred LS state and entropically favorable HS state is accompanied by dramatic changes in the metal-ligand bond lengths, unit cell volume, optical absorption spectrum, and magnetic susceptibility.1 These changes make SCO materials suitable for applications in sensors, memory, and display devices. One of the central challenges in the SCO research is to initiate strongly cooperative interactions known to lead to abrupt spin transitions and thermal hysteresis that can be harvested as a memory effect. One of the strategies to enhance the cooperativity is to design SCO complexes with supramolecular interactions such as π-stacking of aromatic fragments or hydrogen bonding.2 In this work, we report syntheses and characterization of heteroleptic complexes of [Fe(tpma)(L)](ClO4)2 (tpma = tris(pyridin-2-ylmethyl)amine) with novel π-extended biimidazole-type ligands (L) bearing 2,3-dimethyl-naphthalene-, 6,7-dimethyl-2,3-diphenyl-quinoxaline, and 2,3-dimethyl-anthracene pendant fragments. Solvent-free naphthalene-functionalized complex [Fe(tpma)(xnap-bim)](ClO4)2 exhibits abrupt spin transition at T1/2 = 127K with a narrow 1 K hysteresis loop. In contrast, polymorph of this complex that contains one interstitial molecules of pyridine exhibits gradual SCO. Anthracene-functionalized complex [Fe(tpma)(anthra-bim)](ClO4)2 also crystallizes as two polymorphs. Structural studies at 100, 230, and 300 K revealed dramatic changes in the N-Fe-N biting angles and Fe-N distances, indicating the occurrence of temperature-induced SCO. Complex [Fe(tpma)(quin-bim)](ClO4)2 (quin-bim = 6,7-dimethyl-2,3-diphenyl-quinoxaline-2,2’-biimidazole) showed only HS state at 100 and 230 K. In the crystal packing the mononuclear cations form stacks along b axis. We discuss how the observed magnetic behavior correlates with changes in the crystal packing and interactions between the pendant aromatic substituents on the aforementioned complexes. 
    more » « less
  4. Clathrate phases with crystal structures exhibiting complex disorder have been the subject of many prior studies. Here we report syntheses, crystal and electronic structure, and chemical bonding analysis of a Li-substituted Ge-based clathrate phase with the refined chemical formula Ba8Li5.0(1)Ge41.0, which is a rare example of ternary clathrate-I where alkali metal atoms substitute framework Ge atoms. Two different synthesis methods to grow single crystals of the new clathrate phase are presented, in addition to the classical approach towards polycrystalline materials by combining pure elements in desired stoichiometric ratios. Structure elucidations for samples from different batches were carried out by single-crystal and powder X-ray diffraction methods. The ternary Ba8Li5.0(1)Ge41.0 phase crystallizes in the cubic type-I clathrate structure (space group no. 223, a  10.80 Å), with the unit cell being substantially larger compared to the binary phase Ba8Ge43 (Ba8□3Ge43, a  10.63 Å). The expansion of the unit cell is the result of the Li atoms filling vacancies and substituting atoms in the Ge framework, with Li and Ge co-occupying one crystallographic (6c) site. As such, the Li atoms are situated in four-fold coordination environment surrounded by equidistant Ge atoms. Analysis of chemical bonding applying the electron density/ electron localizability approach reveals ionic interaction of barium with the Li–Ge framework, while the lithium-germanium bonds are strongly polar covalent. 
    more » « less
  5. Bond order quantifies the number of electrons dressed-exchanged between two atoms in a material and is important for understanding many chemical properties. Diatomic molecules are the smallest molecules possessing chemical bonds and play key roles in atmospheric chemistry, biochemistry, lab chemistry, and chemical manufacturing. Here we quantum-mechanically calculate bond orders for 288 diatomic molecules and ions. For homodiatomics, we show bond orders correlate to bond energies for elements within the same chemical group. We quantify and discuss how semicore electrons weaken bond orders for elements having diffuse semicore electrons. Lots of chemistry is effected by this. We introduce a first-principles method to represent orbital-independent bond order as a sum of orbital-dependent bond order components. This bond order component analysis (BOCA) applies to any spin-orbitals that are unitary transformations of the natural spin-orbitals, with or without periodic boundary conditions, and to non-magnetic and (collinear or non-collinear) magnetic materials. We use this BOCA to study all period 2 homodiatomics plus Mo 2 , Cr 2 , ClO, ClO − , and Mo 2 (acetate) 4 . Using Manz's bond order equation with DDEC6 partitioning, the Mo–Mo bond order was 4.12 in Mo 2 and 1.46 in Mo 2 (acetate) 4 with a sum of bond orders for each Mo atom of ∼4. Our study informs both chemistry research and education. As a learning aid, we introduce an analogy between bond orders in materials and message transmission in computer networks. We also introduce the first working quantitative heuristic model for all period 2 homodiatomic bond orders. This heuristic model incorporates s–p mixing to give heuristic bond orders of ¾ (Be 2 ), 1¾ (B 2 ), 2¾ (C 2 ), and whole number bond orders for the remaining period 2 homodiatomics. 
    more » « less