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 more » 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. « less
Award ID(s):
Publication Date:
Journal Name:
RSC Adv.
Page Range or eLocation-ID:
45552 to 45581
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 canmore »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.« 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, nickelmore »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.« less
  3. 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 thatmore »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.« less
  4. 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 amore »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.« less
  5. The Kohn–Sham potentialveff(r)is the effective multiplicative operator in a noninteracting Schrödinger equation that reproduces the ground-state density of a real (interacting) system. The sizes and shapes of atoms, molecules, and solids can be defined in terms of Kohn–Sham potentials in a nonarbitrary way that accords with chemical intuition and can be implemented efficiently, permitting a natural pictorial representation for chemistry and condensed-matter physics. Letϵmaxbe the maximum occupied orbital energy of the noninteracting electrons. Then the equationveff(r)=ϵmaxdefines the surface at which classical electrons with energyϵϵmaxwould be turned back and thus determines the surface of any electronic object. Atomic and ionic radii defined in this manner agree well with empirical estimates, show regular chemical trends, and allow one to identify the type of chemical bonding between two given atoms by comparing the actual internuclear distance to the sum of atomic radii. The molecular surfaces can be fused (for a covalent bond), seamed (ionic bond), necked (hydrogen bond), or divided (van der Waals bond). This contribution extends the pioneering work of Z.-Z. Yang et al. [Yang ZZ, Davidson ER (1997)Int J Quantum Chem62:47–53; Zhao DX, et al. (2018)Mol Phys116:969–977] by our consideration of the Kohn–Sham potential, protomolecules,more »doubly negative atomic ions, a bond-type parameter, seamed and necked molecular surfaces, and a more extensive table of atomic and ionic radii that are fully consistent with expected periodic trends.

    « less