skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The location of the chemical bond. Application of long covalent bond theory to the structure of silica
Oxygen is the most abundant terrestrial element and is found in a variety of materials, but still wanting is a universal theory for the stability and structural organization it confers. Herein, a computational molecular orbital analysis elucidates the structure, stability, and cooperative bonding of α-quartz silica (SiO 2 ). Despite geminal oxygen-oxygen distances of 2.61–2.64 Å, silica model complexes exhibit anomalously large O-O bond orders (Mulliken, Wiberg, Mayer) that increase with increasing cluster size—as the silicon-oxygen bond orders decrease. The average O-O bond order in bulk silica computes to 0.47 while that for Si-O computes to 0.64. Thereby, for each silicate tetrahedron, the six O-O bonds employ 52% (5.61 electrons) of the valence electrons, while the four Si-O bonds employ 48% (5.12 electrons), rendering the O-O bond the most abundant bond in the Earth’s crust. The isodesmic deconstruction of silica clusters reveals cooperative O-O bonding with an O-O bond dissociation energy of 4.4 kcal/mol. These unorthodox, long covalent bonds are rationalized by an excess of O 2 p –O 2 p bonding versus anti-bonding interactions within the valence molecular orbitals of the SiO 4 unit (48 vs. 24) and the Si 6 O 6 ring (90 vs. 18). Within quartz silica, oxygen 2 p orbitals contort and organize to avoid molecular orbital nodes, inducing the chirality of silica and resulting in Möbius aromatic Si 6 O 6 rings, the most prevalent form of aromaticity on Earth. This long covalent bond theory (LCBT) relocates one-third of Earth’s valence electrons and indicates that non-canonical O-O bonds play a subtle, but crucial role in the structure and stability of Earth’s most abundant material.  more » « less
Award ID(s):
1904768
PAR ID:
10451370
Author(s) / Creator(s):
Date Published:
Journal Name:
Frontiers in Chemistry
Volume:
11
ISSN:
2296-2646
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The origin in deshielding of 29 Si NMR chemical shifts in R 3 Si–X, where X = H, OMe, Cl, OTf, [CH 6 B 11 X 6 ], toluene, and O X (O X = surface oxygen), as well as i Pr 3 Si + and Mes 3 Si + were studied using DFT methods. At the M06-L/6-31G(d,p) level of theory the geometry optimized structures agree well with those obtained experimentally. The trends in 29 Si NMR chemical shift also reproduce experimental trends; i Pr 3 Si–H has the most shielded 29 Si NMR chemical shift and free i Pr 3 Si + or isolable Mes 3 Si + have the most deshielded 29 Si NMR chemical shift. Natural localized molecular orbital (NLMO) analysis of the contributions to paramagnetic shielding ( σ p ) in these compounds shows that Si–R (R = alkyl, H) bonding orbitals are the major contributors to deshielding in this series. The Si–R bonding orbitals are coupled to the empty p-orbital in i Pr 3 Si + or Mes 3 Si + , or to the orbital in R 3 Si–X. This trend also applies to surface bound R 3 Si–O X . This model also explains chemical shift trends in recently isolated t Bu 2 SiH 2 + , t BuSiH 2 + , and SiH 3 + that show more shielded 29 Si NMR signals than R 3 Si + species. There is no correlation between isotropic 29 Si NMR chemical shift and charge at silicon. 
    more » « less
  2. Abstract Several renewable energy schemes aim to use the chemical bonds in abundant molecules like water and ammonia as energy reservoirs. Because the O-H and N-H bonds are quite strong (>100 kcal/mol), it is necessary to identify substances that dramatically weaken these bonds to facilitate proton-coupled electron transfer processes required for energy conversion. Usually this is accomplished through coordination-induced bond weakening by redox-active metals. However, coordination-induced bond weakening is difficult with earth’s most abundant metal, aluminum, because of its redox inertness under mild conditions. Here, we report a system that uses aluminum with a redox non-innocent ligand to achieve significant levels of coordination-induced bond weakening of O-H and N-H bonds. The multisite proton-coupled electron transfer manifold described here points to redox non-innocent ligands as a design element to open coordination-induced bond weakening chemistry to more elements in the periodic table. 
    more » « less
  3. Abstract What happens when a C−H bond is forced to interact with unpaired pairs of electrons at a positively charged metal? Such interactions can be considered as “contra‐electrostatic” H‐bonds, which combine the familiar orbital interaction pattern characteristic for the covalent contribution to the conventional H‐bonding with an unusual contra‐electrostatic component. While electrostatics is strongly stabilizing component in the conventional C−H⋅⋅⋅X bonds where X is an electronegative main group element, it is destabilizing in the C−H⋅⋅⋅M contacts when M is Au(I), Ag(I), or Cu(I) of NHC−M−Cl systems. Such remarkable C−H⋅⋅⋅M interaction became experimentally accessible within (α‐ICyDMe)MCl, NHC‐Metal complexes embedded into cyclodextrins. Computational analysis of the model systems suggests that the overall interaction energies are relatively insensitive to moderate variations in the directionality of interaction between a C−H bond and the metal center, indicating stereoelectronic promiscuity of fully filled set ofd‐orbitals. A combination of experimental and computational data demonstrates that metal encapsulation inside the cyclodextrin cavity forces the C−H bond to point toward the metal, and reveals a still attractive “contra‐electrostatic” H‐bonding interaction. 
    more » « 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 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
  5. null (Ed.)
    The degree to which the Earth’s mantle stores and cycles water in excess of the storage capacity of nominally anhydrous minerals is dependent upon the stability of hydrous phases under mantle-relevant pressures, temperatures, and compositions. Two hydrous phases, phase D and phase H, are stable to the pressures and temperatures of the Earth’s lower mantle, suggesting that the Earth’s lower mantle may participate in the cycling of water. We build on our prior work of density functional theory calculations on phase H with the stability, structure, and bonding of hydrous phases D, and we predict the aluminum partitioning with H in the Al 2 O 3 -SiO 2 -MgO-H 2 O system. We address the solid solutions through a statistical sampling of site occupancy and calculation of the partition function from the grand canonical ensemble. We show that each phase has a wide solid solution series between MgSi 2 O 6 H 2 -Al 2 SiO 6 H 2 and MgSiO 4 H 2 -2 δ AlOOH + SiO 2 , in which phase H is more aluminum rich than phase D at a given bulk composition. We predict that the addition of Al to both phases D and H stabilizes each phase to higher temperatures through additional configurational entropy. While we have shown that phase H does not exhibit symmetric hydrogen bonding at high pressure, we report here that phase D undergoes a gradual increase in the number of symmetric H-bonds beginning at ∼30 GPa, and it is only ∼50% complete at 60 GPa. 
    more » « less