skip to main content


Title: A First-Principles Exploration of NaxSy Binary Phases at 1 atm and Under Pressure
Interest in Na-S compounds stems from their use in battery materials at 1 atm, as well as the potential for superconductivity under pressure. Evolutionary structure searches coupled with Density Functional Theory calculations were employed to predict stable and low-lying metastable phases of sodium poor and sodium rich sulfides at 1 atm and within 100–200 GPa. At ambient pressures, four new stable or metastable phases with unbranched sulfur motifs were predicted: Na2S3 with C 2 / c and Imm2 symmetry, C 2 -Na2S5 and C 2 -Na2S8. Van der Waals interactions were shown to affect the energy ordering of various polymorphs. At high pressure, several novel phases that contained a wide variety of zero-, one-, and two-dimensional sulfur motifs were predicted, and their electronic structures and bonding were analyzed. At 200 GPa, P 4 / m m m -Na2S8 was predicted to become superconducting below 15.5 K, which is close to results previously obtained for the β -Po phase of elemental sulfur. The structures of the most stable M3S and M4S, M = Na, phases differed from those previously reported for compounds with M = H, Li, K.  more » « less
Award ID(s):
1827815
NSF-PAR ID:
10120800
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Crystals
Volume:
9
Issue:
9
ISSN:
2073-4352
Page Range / eLocation ID:
441
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. With the motivation of searching for new superconductors in the Mg–B system, we performed ab initio evolutionary searches for all the stable compounds in this binary system in the pressure range of 0–200 GPa. We found previously unknown, yet thermodynamically stable, compositions MgB 3 and Mg 3 B 10 . Experimentally known MgB 2 is stable in the entire pressure range 0–200 GPa, while MgB 7 and MgB 12 are stable at pressures below 90 GPa and 35 GPa, respectively. We predict a reentrant behavior for MgB 4 , which becomes unstable against decomposition into MgB 2 and MgB 7 at 4 GPa and then becomes stable above 61 GPa. We find ubiquity of phases with boron sandwich structures analogous to the AlB 2 -type structure. However, with the exception of MgB 2 , all other magnesium borides have low electron–phonon coupling constants λ of 0.32–0.39 and are predicted to have T c below 3 K. 
    more » « less
  2. Abstract

    Here, the combination of theoretical computations followed by rapid experimental screening and in situ diffraction studies is demonstrated as a powerful strategy for novel compounds discovery. When applied for the previously “empty” Na−Zn−Bi system, such an approach led to four novel phases. The compositional space of this system was rapidly screened via the hydride route method and the theoretically predicted NaZnBi (PbClF type,P4/nmm) and Na11Zn2Bi5(Na11Cd2Sb5type,P) phases were successfully synthesized, while other computationally generated compounds on the list were rejected. In addition, single crystal X‐ray diffraction studies of NaZnBi indicate minor deviations from the stoichiometric 1 : 1 : 1 molar ratio. As a result, two isostructural (PbClF type,P4/nmm) Zn‐deficient phases with similar compositions, but distinctly different unit cell parameters were discovered. The vacancies on Zn sites and unit cell expansion were rationalized from bonding analysis using electronic structure calculations on stoichiometric “NaZnBi”.In‐situsynchrotron powder X‐ray diffraction studies shed light on complex equilibria in the Na−Zn−Bi system at elevated temperatures. In particular, the high‐temperature polymorphHT‐Na3Bi (BiF3type,Fmm) was obtained as a product of Na11Zn2Bi5decomposition above 611 K.HT‐Na3Bi cannot be stabilized at room temperature by quenching, and this type of structure was earlier observed in the high‐pressure polymorphHP‐Na3Bi above 0.5 GPa. The aforementioned approach of predictive synthesis can be extended to other multinary systems.

     
    more » « less
  3. null (Ed.)
    X-ray diffraction indicates that the structure of the recently discovered carbonaceous sulfur hydride (C-S-H) room temperature superconductor is derived from previously established van der Waals compounds found in the H2S-H2 and CH4-H2 systems. Crystals of the superconducting phase were produced by a photochemical synthesis technique leading to the superconducting critical temperature Tc of 288 K at 267 GPa. X-ray diffraction patterns measured from 124 to 178 GPa, within the pressure range of the superconducting phase, are consistent with an orthorhombic structure derived from the Al2Cu-type determined for (H2S)2H2 and (CH4)2H2 that differs from those predicted and observed for the S-H system to these pressures. The formation and stability of the C-S-H compound can be understood in terms of the close similarity in effective volumes of the H2S and CH4 components, and denser carbon-bearing S-H phases may form at higher pressures. The results are crucial for understanding the very high superconducting Tc found in the C-S-H system at megabar pressures. 
    more » « less
  4. Abstract

    The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X-ray absorption near-edge structure spectroscopy (S-XANES) where S 6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms.

    One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel.

    The thermodynamics of incorporation of disulfide and bisulfide into apatite is evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S is also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS(aq)(2−n) and HnS(aq)(2−n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite, and decreases for hydroxylapatite, as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.

     
    more » « less
  5. The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X‑ray absorption near-edge structure spectroscopy (S-XANES) where S6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms. One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel. The thermodynamics of incorporation of disulfide and bisulfide into apatite are evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S are also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS2(aq)(2–n) and HnS(aq)(2–n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite and decreases for hydroxylapatite as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended also to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures. 
    more » « less