skip to main content

Title: First principles study of dense and metallic nitric sulfur hydrides
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 could more » have implications for their magnetic field formation. « less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Communications Chemistry
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.
  2. 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 Hmore »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.« less
  3. 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
  4. Formation of vitreous ice during rapid compression of water at room temperature is important for biology and the study of biological systems. Here, we show that Raman spectra of rapidly compressed water at greater than 1 GPa at room temperature exhibits the signature of high-density amorphous ice, whereas the X-ray diffraction (XRD) pattern is dominated by crystalline ice VI. To resolve this apparent contradiction, we used molecular dynamics simulations to calculate full vibrational spectra and diffraction patterns of mixtures of vitreous ice and ice VI, including embedded interfaces between the two phases. We show quantitatively that Raman spectra, which probe the local polarizability with respect to atomic displacements, are dominated by the vitreous phase, whereas a small amount of the crystalline component is readily apparent by XRD. The results of our combined experimental and theoretical studies have implications for detecting vitreous phases of water, survival of biological systems under extreme conditions, and biological imaging. The results provide additional insight into the stable and metastable phases of H 2 O as a function of pressure and temperature, as well as of other materials undergoing pressure-induced amorphization and other metastable transitions.
  5. Sub-Neptunes are common among the discovered exoplanets. However, lack of knowledge on the state of matter inH2O-rich setting at high pressures and temperatures (PT) places important limitations on our understanding of this planet type. We have conducted experiments for reactions betweenSiO2andH2O as archetypal materials for rock and ice, respectively, at highPT. We found anomalously expanded volumes of dense silica (up to 4%) recovered from hydrothermal synthesis above ∼24 GPa where theCaCl2-type (Ct) structure appears at lower pressures than in the anhydrous system. Infrared spectroscopy identified strong OH modes from the dense silica samples. Both previous experiments and our density functional theory calculations support up to 0.48 hydrogen atoms per formula unit of (Si1xH4x)O2(x=0.12). At pressures above 60 GPa,H2O further changes the structural behavior of silica, stabilizing a niccolite-type structure, which is unquenchable. From unit-cell volume and phase equilibrium considerations, we infer that the niccolite-type phase may contain H with an amount at least comparable with or higher than that of the Ct phase. Our results suggest that the phases containing both hydrogen and lithophile elements could bemore »the dominant materials in the interiors of water-rich planets. Even for fully layered cases, the large mutual solubility could make the boundary between rock and ice layers fuzzy. Therefore, the physical properties of the new phases that we report here would be important for understanding dynamics, geochemical cycle, and dynamo generation in water-rich planets.

    « less