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.
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Page Range or eLocation-ID:
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.
  2. 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, themore »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.« less
  3. 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.
  4. The high-pressure structure and stability of the calcic amphibole tremolite (Ca2Mg5Si8O22(OH)2) was investigated to ~40 GPa at 300 K by single-crystal X-ray diffraction using synchrotron radiation. C2/m symmetry tremolite displays a broader metastability range than previously studied clinoamphiboles, exhibiting no first-order phase transition up to 40 GPa. Axial parameter ratios a/b and a/c, in conjunction with finite strain versus normalized pressure trends, indicate that changes in compressional behavior occur at pressures of ~5 and ~20 GPa. An analysis of the finite strain trends, using third-order Birch-Murnaghan equations of state, resulted in bulk moduli (𝐾) of 72(7), 77(2), and 61(1) GPa for the compressional regimes from 0-5 GPa (regime I), 5-20 GPa (II), and above 20 GPa (III), respectively, and accompanying pressure-derivatives of the bulk moduli (𝐾′) of 8.6(42), 6.0(3), and 10.0(2). The results are consistent with first-principle theoretical calculations of tremolite elasticity. The axial compressibility ratios of tremolite, determined as 𝛽a : 𝛽b : 𝛽c = 2.22:1.0:0.78 (regime I), 2.12:1.0:0.96 (II), and 1.03:1.0:0.75 (III), demonstrate a substantial reduction of the compressional anisotropy of tremolite at high pressures, which is a notable contrast with the increasingly anisotropic compressibility observed in the high-pressure polymorphs of the clinoamphibole grunerite. The shift in compression-regimemore »at 5 GPa (I-II) transition is ascribed to stiffening along the crystallographic a-axis corresponding to closure of the vacant A-site in the structure, and a shift in the topology of the a-oriented surfaces of the structural I-beam from concave to convex. The II-III regime shift at 20 GPa corresponds to an increasing rate of compaction of the Ca-polyhedra and increased distortion of the Mg-octahedral sites, processes which dictate compaction in both high-pressure compression-regimes. Bond-valence analyses of the tremolite structure under pressure show dramatic overbonding of the Ca-cations (75% at 30 GPa), with significant Mg-cation overbonding as well (40%). These imply that tremolite’s notable metastability range hinges on the calcium cation’s bonding environment. The 8-fold coordinated Ca-polyhedron accommodates significant compaction under pressure, while the geometry of the Ca-O polyhedron becomes increasingly regular and inhibits the reorientation of the tetrahedral chains that generate phase transitions observed in other clinoamphiboles. Peak/background ratio of diffraction data collected above 40 GPa and our equation of state determination of bulk moduli and compressibilities of tremolite in regime III, in concert with the results of our previous Raman study, suggest that C2/m tremolite may be approaching the limit of its metastability above 40 GPa. Our results have relevance for both the metastable compaction of tremolite during impact events, and for possible metastable persistence of tremolite within cold subduction zones within the Earth.« less
  5. Abstract Studies of molecular mixtures containing hydrogen sulfide (H 2 S) could open up new routes towards hydrogen-rich high-temperature superconductors under pressure. H 2 S and ammonia (NH 3 ) form hydrogen-bonded molecular mixtures at ambient conditions, but their phase behavior and propensity towards mixing under pressure is not well understood. Here, we show stable phases in the H 2 S–NH 3 system under extreme pressure conditions to 4 Mbar from first-principles crystal structure prediction methods. We identify four stable compositions, two of which, (H 2 S) (NH 3 ) and (H 2 S) (NH 3 ) 4 , are stable in a sequence of structures to the Mbar regime. A re-entrant stabilization of (H 2 S) (NH 3 ) 4 above 300 GPa is driven by a marked reversal of sulfur-hydrogen chemistry. Several stable phases exhibit metallic character. Electron–phonon coupling calculations predict superconducting temperatures up to 50 K, in the Cmma phase of (H 2 S) (NH 3 ) at 150 GPa. The present findings shed light on how sulfur hydride bonding and superconductivity are affected in molecular mixtures. They also suggest a reservoir for hydrogen sulfide in the upper mantle regions of icy planets in a potentially metallic mixture, which couldmore »have implications for their magnetic field formation.« less