An official website of the United States government
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
Introducing DDEC6 atomic population analysis: part 4. Efficient parallel computation of net atomic charges, atomic spin moments, bond orders, and more
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
- Award ID(s):
- 1555376
- PAR ID:
- 10055288
- Date Published:
- Journal Name:
- RSC Advances
- Volume:
- 8
- Issue:
- 5
- ISSN:
- 2046-2069
- Page Range / eLocation ID:
- 2678 to 2707
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
We present a new implementation for computing spin–orbit couplings (SOCs) within a time-dependent density-functional theory (TD-DFT) framework in the standard spin-conserving formulation as well in the spin–flip variant (SF-TD-DFT). This approach employs the Breit–Pauli Hamiltonian and Wigner–Eckart’s theorem applied to the reduced one-particle transition density matrices, together with the spin–orbit mean-field treatment of the two-electron contributions. We use a state-interaction procedure and compute the SOC matrix elements using zero-order non-relativistic states. Benchmark calculations using several closed-shell organic molecules, diradicals, and a single-molecule magnet illustrate the efficiency of the SOC protocol. The results for organic molecules (described by standard TD-DFT) show that SOCs are insensitive to the choice of the functional or basis sets, as long as the states of the same characters are compared. In contrast, the SF-TD-DFT results for small diradicals (CH 2 , [Formula: see text], SiH 2 , and [Formula: see text]) show strong functional dependence. The spin-reversal energy barrier in a Fe(III) single-molecule magnet computed using non-collinear SF-TD-DFT (PBE0, ωPBEh/cc-pVDZ) agrees well with the experimental estimate.more » « less
-
Abstract Rare-earth monopnictides are a family of materials simultaneously displaying complex magnetism, strong electronic correlation, and topological band structure. The recently discovered emergent arc-like surface states in these materials have been attributed to the multi-wave-vector antiferromagnetic order, yet the direct experimental evidence has been elusive. Here we report observation of non-collinear antiferromagnetic order with multiple modulations using spin-polarized scanning tunneling microscopy. Moreover, we discover a hidden spin-rotation transition of single-to-multiple modulations 2 K below the Néel temperature. The hidden transition coincides with the onset of the surface states splitting observed by our angle-resolved photoemission spectroscopy measurements. Single modulation gives rise to a band inversion with induced topological surface states in a local momentum region while the full Brillouin zone carries trivial topological indices, and multiple modulation further splits the surface bands via non-collinear spin tilting, as revealed by our calculations. The direct evidence of the non-collinear spin order in NdSb not only clarifies the mechanism of the emergent topological surface states, but also opens up a new paradigm of control and manipulation of band topology with magnetism.more » « less
-
Abstract Emerging from the intricate interplay of topology and magnetism, the giant anomalous Hall effect (AHE) is the most known topological property of the recently discovered kagomé ferromagnetic Weyl semimetal Co3Sn2S2with the magnetic Co atoms arranged on a kagomé lattice. Here it is reported that the AHE in Co3Sn2S2can be fine‐tuned by an applied magnetic field orientated within ≈2° of the kagomé plane, while beyond this regime, it stays unchanged. Particularly, it can vanish in magnetic fields parallel to the kagomé plane and even decrease in magnetic fields collinear with the spin direction. This tunable AHE can be attributed to local spin switching enabled by the geometrical frustration of the magnetic kagomé lattice, revealing that spins in a kagomé ferromagnet change their switching behavior as the magnetic field approaches the kagomé plane. These results also suggest a versatile way to tune the properties of a kagomé magnet.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/c7ra11829e
