Structures and Energetics of Elemental Sulfur in Hydrogen Sulfide

Formation pathways of hydrogen polysulfides in sulfur-bearing natural gas reservoirs from density functional theory calculations

CONTEXT

The interaction mechanisms between a sulfur atom (S) and hydrogen sulfide (H2S), as well as the formation and stability of H2Sn (n = 2-9), are fundamental to understanding sulfur chemistry in natural gas reservoirs. Despite their importance, the abiogenic origins and reaction pathways of H2Sn in natural gas fields remain inadequately understood. Clarifying these mechanisms is essential for addressing sulfur deposition challenges, which have direct implications for extraction efficiency, operational safety, and reservoir management.

METHODS

This study utilized quantum chemistry calculations to systematically investigate the reaction mechanisms between sulfur atoms and hydrogen sulfide, with a particular focus on the formation of H2Sn. Transition state (TS) searches were conducted to identify energetically favorable reaction pathways, and intrinsic reaction coordinate (IRC) analyses were performed to validate the reaction trajectories. The kinetics and thermodynamics of H2S2 formation from elemental sulfur and H2S were comprehensively evaluated. Additionally, stability analyses were carried out to assess the relative stability of H2Sn under varying reservoir conditions, offering insights into their decomposition tendencies and subsequent formation of H2S and elemental sulfur (S8).

Molecular mechanism of elemental sulfur dissolution in H2S under stratal conditions

Resulting from the solubility reduction of elemental sulfur during the development of high sulfur gas formations, sulfur deposition often occurs to reduce the gas production and threaten the safety of gas wells. Understanding the dissolution mechanism of elemental sulfur in natural gas is essential to reduce the risk caused by sulfur deposition. Because of the harsh conditions in the high-sulfur formations, it remains challenging to in situ characterize the dissolution-precipitation processes, making deficient the knowledge of sulfur dissolution mechanism. The dissolution of sulfur allotropes (SN, N = 2, 4, 6 and 8) in H2S, the main solvent of sulfur in natural gas, is studied in this work by means of first-principles calculations and molecular dynamics simulations. While S6 and S8 undergo physical interaction with H2S under the conditions corresponding to those at 1-6 km stratigraphic depths, S2 and S4 react with H2S and form stable polysulfides. Unravelling the mechanism would be helpful for understanding and controlling the sulfur deposition in high-sulfur gas development.

Solubility evolution of elemental sulfur in natural gas with a varying H2S content

CONTEXT

The solubility variations of elemental sulfur are of great importance in the prevention of sulfur deposition during the development of high-sulfur gas formations. It has been observed that the solubility varies with H2S content, which is the main solvent of elemental sulfur in natural gas. Moreover, the addition of small amounts of CH4 and/or CO2 in H2S leads to a dramatic solubility reduction of which the mechanism remains unclear. Using a modified direct coexistence method, molecular dynamics simulations are conducted to uncover the molecular mechanism of the solubility reduction. The observed solubility variations with H2S content are reproduced, and the solubility reduction is interpreted by the antisolvent effect of CH4 and CO2. While the H2S content varies in a wide range in the known high-sulfur gas formation, our simulations provide useful information for controlling the sulfur deposition in gas development.

METHODS

Molecular dynamics simulations are carried out using the LAMMPS package. The initial models are constructed with the Packmol software. The Ballone and Jones' potential is used for S8, the Galliero's potential for H2S and CO2, and the transferable potentials for phase equilibria-united atom (TraPPE-UA) force field for CH4. The time step is set to 1 fs, and the molecular trajectories of additional 2 ns after equilibrium are collected for analysis.

DFT Model of Elemental Sulfur in Sulfolane Solutions

The structure and the thermodynamic and optical (UV) properties of elemental sulfur solution in sulfolane (Sl) have been studied using density functional theory methods. The cyclic molecular form of sulfur (S8 "crown") was found using PBE1PBE/6-311+G(d,p) approximation in combination with a polarizable continuum model (the integral equation formalism variant) to exist in sulfolane medium as a Sl-S8-Sl solvate. It has been theoretically established that sulfur can form stable (S8)n clusters in concentrated solutions. An increase in the extent of association (n) of the sulfur cluster leads to a decrease in the extinction coefficient [TD-DFT(TPSSTPSS/6-311+G(d,p))] of the most intense absorption maximum lying at about 50,000 cm-1 while maintaining the shape of the remaining part of the spectrum. The observed pattern qualitatively expresses the spectral regularities of solutions with different concentrations of sulfur in sulfolane. It has been proposed that a model of the absorption spectrum of elemental sulfur suggests a minor contribution of the S12 molecular form (G298°((S12)2) - G298°((S8)3) ≈ -15.5 kJ mol-1). The findings of the study will provide deeper insights into the transformation of molecular forms of sulfur and more precisely analyze processes involving sulfur as an acting species using electronic spectroscopy.