---
title: Sulphur Atom Popping Model Mechanisms
url: https://www.emergentmind.com/topics/sulphur-atom-popping-model
type: topic
---

# Sulphur Atom Popping Model Mechanisms

Searching arXiv for the cited papers and related usage of the term.
arXiv search query: 2111.08109 OR 1406.2278 OR 2509.20179 OR 2212.13538
“Sulphur atom popping model” denotes several sulfur-centered mechanisms in which sulfur or sulfur-bearing species undergo a discrete release, relocation, or state-selection event that then controls a larger chemical or electronic response. In the Orion KL astrochemical literature, the term refers to a two-phase UCL_CHEM scenario in which sulfur is depleted onto grains during collapse, converted predominantly into H\(_2\)S on ice mantles, and then returned to the gas during protostellar heating or C-type shocks [1406.2278]. In monolayer MoS\(_2\) memristors, it refers to an intrinsic electric-field-driven process in which a sulfur atom adjacent to a vacancy is displaced from the chalcogen plane into the Mo plane, generating a metastable, charge-trapping defect associated with nonvolatile resistive switching [2509.20179]. In dissociative recombination of SH\(^+\), an analogous “popping” picture describes the state-selective appearance of atomic sulfur in either the ground \(^{3}P\) or the lowest excited \(^{1}D\) state, with the branching controlled by dissociative pathways and Landau–Zener flux redistribution [2212.13538]. A separate sulfur literature models liquid polyamorphism through maximum valence and polymerization-driven phase separation rather than through a popping mechanism [2111.08109].

## 1. Terminological scope and domains of use

In the cited literature, the expression is not attached to a single universal formalism. Instead, it labels distinct sulfur-centered events at different scales, from grain-surface desorption chemistry to defect-mediated electronic switching and state-resolved atomic production.

| Context | What “pops” | Governing picture |
|---|---|---|
| Orion KL astrochemistry | Sulfur-bearing mantle species return to gas | Phase I freeze-out, Phase II sublimation during warm-up or shock |
| MoS\(_2\) memristor | A sulfur atom moves into the Mo plane at a vacancy | Field-driven reactive dynamics plus charge trapping |
| SH\(^+\) dissociative recombination | Atomic sulfur emerges in a specific electronic state | MQDT branching and Landau–Zener redistribution |
| Liquid sulfur polyamorphism | No popping event is invoked | Maximum-valence polymerization and LLPT |

This multiplicity matters because superficially similar language can conceal very different ontologies. In Orion KL, the relevant event is thermal or shock-induced release of mantle material. In MoS\(_2\), it is a local lattice rearrangement under electric field. In SH\(^+\) recombination, it is a state-selective dissociation product channel. The sulfur polyamorphism model, by contrast, concerns reversible covalent bonding and an additional attraction among fully bonded atoms, not a sulfur-pop process [2111.08109].

## 2. Orion KL: freeze-out, mantle storage, and thermal or shock release

The Orion KL formulation is a two-stage gas-grain model implemented with UCL_CHEM. Phase I is a modified collapse with freeze-out of gas-phase material onto grain mantles; Phase II begins either when a central protostar forms in the hot core or when a non-dissociative C-type shock passes through the plateau [1406.2278]. The model adopts initial abundances relative to H nuclei of H \(=1.0\), He \(=0.075\), O \(=4.45\times10^{-4}\), C \(=1.79\times10^{-4}\), N \(=8.52\times10^{-5}\), and S \(=1.43\times10^{-5}\), and also explores sulfur reduced to \(0.1\,S_\odot\) and \(0.01\,S_\odot\).

A central control parameter in Phase I is the accretion efficiency,
\[
f_r = \frac{mCO}{mCO+CO_{\mathrm{gas}}},
\]
where \(mCO\) is mantle CO at the end of Phase I. The explored ranges are \(f_r=0.30\) to \(0.85\) for the hot core and \(f_r=0.85\) for the plateau. Non-thermal desorption is switched off.

The model’s key sulfur assumption is that most sulfur is ionised during Phase I. This is tied to an inhomogeneous, clumpy cloud in which UV photons penetrate and ionise sulfur; the paper emphasizes sulfur’s ionisation potential of \(10.36\) eV and states that the collision timescale of S\(^+\) with negatively charged grains is about an order of magnitude smaller than for neutral carbon-, oxygen-, or nitrogen-bearing species. The consequence is efficient freeze-out of S\(^+\), followed by grain-surface conversion mainly into H\(_2\)S and secondarily into OCS.

The preferred mantle composition at the end of Phase I is reported as \(m(\mathrm{H_2S})=1.357\times10^{-6}\), \(m(\mathrm{OCS})=7.300\times10^{-8}\), and \(m(\mathrm{H_2O})=4.278\times10^{-4}\), giving \(m(\mathrm{H_2S})/m(\mathrm{H_2O})=3\times10^{-3}\) and \(m(\mathrm{OCS})/m(\mathrm{H_2O})=2\times10^{-4}\). The paper states that this corresponds to “\(\sim95\%\)” of frozen-out sulfur ending up as H\(_2\)S when S\(^+\) is the accreting ion. A recurrent clarification follows from this: the model is not a literal release of atomic S from grains, but a release of sulfur primarily as H\(_2\)S, with OCS as a secondary reservoir.

Phase II assumes \(100\%\) sublimation efficiency from the grains as the region warms. In the hot core, a central protostar forms, the temperature rises to \(300\) K, and evaporation is time dependent rather than instantaneous; the warm-up timescale depends on stellar mass, with \(5\), \(10\), and \(15\,M_\odot\) cases modeled. In the plateau, the driver is a non-dissociative C-type shock with velocities of \(\sim20\) km s\(^{-1}\) for \(1000\) K and \(\sim30\) km s\(^{-1}\) for \(2000\) K; the density is held fixed during Phase II at either \(5\times10^5\) cm\(^{-3}\) or \(5\times10^6\) cm\(^{-3}\). The shock cools on the order of hundreds of years, and the dynamical timescale is \(\sim1000\) years.

## 3. Orion KL after release: gas-phase sulfur chemistry and observational constraints

Once mantle species desorb, the post-release sulfur chemistry diverges between hot core and plateau conditions. In the hot core, the dominant routes reported are
\[
\mathrm{O + SH \rightarrow SO + H}
\]
and
\[
\mathrm{OH + SO \rightarrow SO_2 + H},
\]
with a later-time recycling channel
\[
\mathrm{SO_2 + C \rightarrow SO + CO}.
\]
The stated consequence is that SO forms mainly from SH, SO is converted into SO\(_2\) by OH, and the SO\(_2\)/SO ratio increases with time [1406.2278].

In the plateau, the chemistry is explicitly time-stratified. Early after the shock, while the gas remains hot (\(t\lesssim10^3\) yr), the dominant channels are
\[
\mathrm{O_2 + S \rightarrow SO + O}
\]
and
\[
\mathrm{O_2 + SO \rightarrow SO_2 + O}.
\]
After the gas cools (\(t\gtrsim10^3\) yr), the chemistry shifts to
\[
\mathrm{OH + S \rightarrow SO + H}
\]
and
\[
\mathrm{OH + SO \rightarrow SO_2 + H}.
\]
The paper identifies this as a chemical transition from O\(_2\)-dominated to OH-dominated SO/SO\(_2\) formation and states that this transition is currently taking place in Orion KL’s plateau.

The best-fitting hot-core conditions are initial sulfur abundance \(=0.1\,S_\odot\), density \(n_H=10^7\) cm\(^{-3}\), and stellar mass \(10\)–\(15\,M_\odot\). The preferred hot-core realization is Model 5, with \(M_\star=10\,M_\odot\), \(f_r=0.85\), initial sulfur \(=0.1\,S_\odot\), and \(n_H=10^7\) cm\(^{-3}\). Under these conditions the model reproduces all six species well at about \(t\sim6\times10^4\) yr. The paper also notes that \(0.01\,S_\odot\) can work only if the density is increased to \(10^8\) cm\(^{-3}\), whereas solar sulfur produces H\(_2\)S column densities that are too high.

For the plateau, the best fit is Model 4: initial sulfur abundance \(=0.1\,S_\odot\), final density \(=5\times10^6\) cm\(^{-3}\), and maximum shock temperature closer to \(2000\) K than \(1000\) K. Lowering the density to \(5\times10^5\) cm\(^{-3}\) worsens SO and often OCS, and \(0.01\,S_\odot\) underproduces SO and OCS by up to an order of magnitude. The observed source-averaged column densities used for comparison are, for the plateau, SO \(=18\pm10\), SO\(_2\) \(=10\pm4\), CS \(=2.4\pm0.5\), OCS \(=12\pm3\), H\(_2\)CS \(=1.0\pm0.3\), and no H\(_2\)S measurement; for the hot core, SO \(=9\pm4\), SO\(_2\) \(=65\pm20\), CS \(=1.7\pm0.9\), OCS \(=15\pm4\), H\(_2\)CS \(=1.6\pm0.5\), and H\(_2\)S \(=950\pm200\), all in units of \(\times10^{15}\,\mathrm{cm^{-2}}\).

A common misunderstanding is therefore excluded by the model itself: the “sulphur atom popping” picture in Orion KL does not mean direct popping of neutral sulfur atoms from grains. It means efficient freeze-out of ionised sulfur, conversion mainly to H\(_2\)S on mantles, and thermal desorption followed by rapid gas-phase reprocessing.

## 4. MoS\(_2\) memristors: field-driven sulfur displacement and nonvolatile switching

In monolayer MoS\(_2\) memristors, the sulphur atom popping model is an intrinsic switching mechanism in the active layer rather than an electrode-filament model. The proposed event is that a sulfur atom adjacent to a vacancy is pulled by a sufficiently strong electric field out of its original chalcogen plane into the molybdenum plane, producing a metastable popped configuration [2509.20179]. The physical consequences listed in the paper are the appearance of additional defect states near the Fermi level, stronger local bonding that can resemble a “virtual filament,” and charge trapping under bias. The nonvolatile state is therefore a coupled structure-plus-charge configuration.

The mechanism was first established in reactive-force-field molecular dynamics using ReaxFF. The follow-up paper stresses that ReaxFF is essential because it allows bond breaking and bond formation together with dynamic charge redistribution through charge equilibration. Partial charges are updated every timestep by a QEq scheme, so the force on an atom under bias follows the conceptual relation
\[
F = q\cdot E.
\]
Under field, the simulations show sulfur atoms moving from the opposite sulfur plane into the Mo plane at vacancy sites, while neighboring atoms also rearrange in-plane during prolonged field exposure. This dynamic evolution is central to the authors’ interpretation: the popped state is not treated as a geometry imposed a priori, but as a field-induced non-equilibrium product of time-dependent atomic motion and charge response.

The density-functional-theory analysis is used to reinterpret the energetics. In a relaxed \(8\times8\) monolayer MoS\(_2\) supercell with a sulfur vacancy and a manually placed sulfur atom in the Mo plane, the survival of the popped atom after relaxation depends strongly on how the structure is prepared. If the vacancy structure is first relaxed and the sulfur atom is inserted afterward, the atom can revert; if the atom is moved directly into the Mo plane in the vacancy-containing structure and then relaxed, it can remain there. The paper argues from this that earlier claims against the model were inconclusive because the local neighborhood and force environment depend on the relaxation history.

At zero field or without extra charge, the popped configuration is reported to be higher in energy than the vacancy state by about \(2.74\) eV. Under increasing electric field, the popped-state energy drops while the vacancy-state energy changes little, and the CINEB calculations display an energy valley only when the field is present. Charge-density-difference plots show charge accumulation around the popped atom, supporting the interpretation that the defect becomes a favorable electron trap. Because sulfur vacancies in MoS\(_2\) already generate mid-gap states, the popped state is described as creating or enhancing a second vacancy-like trapping environment.

To mimic field-induced electron accumulation after the field is removed, the DFT calculations add electrons manually. With two extra electrons, the popped sulfur atom remains in the Mo plane after relaxation and the added charge localizes around the popped site rather than spreading uniformly. With one extra electron the popped state is unstable, whereas with two or more electrons it becomes stable. The AIMD results further report that at \(300\) K with two extra electrons the popped state survives for about \(0.5\) ps, with three extra electrons it remains stable for about \(1\) ps, and at \(500\) K with three extra electrons it persists for about \(3\) ps. The stated interpretation is that trapped charge stabilizes the popped defect and accounts for persistence after the applied field is removed.

The paper also delineates methodological limits. It argues that ground-state DFT is not well suited to validate a non-equilibrium field-induced state unless the charge state is treated appropriately, that NEB is conceptually questionable when used between a stable vacancy state and an allegedly unstable popped state, and that universal machine-learning interatomic potentials such as M3GNet are unreliable for this problem because they are not charge-aware and are trained mainly on bulk materials rather than on surfaces, interfaces, vacancies, and electric-field-driven reactive processes. These caveats do not eliminate the proposed mechanism; they specify the conditions under which the model is claimed to be meaningful.

## 5. SH\(^+\) dissociative recombination: state-selective sulfur production

A different sulphur-popping usage appears in the dissociative recombination of SH\(^+\) in the diffuse interstellar medium. The reaction considered is
\[
\mathrm{SH}^+(v_i^+) + e^-(\varepsilon) \rightarrow \mathrm{S} + \mathrm{H},
\]
with competing vibrational excitation
\[
\mathrm{SH}^+(v_i^+) + e^-(\varepsilon) \rightarrow \mathrm{SH}^+(v_f^+) + e^-(\varepsilon').
\]
The focus is on vibrationally relaxed SH\(^+\), \(v_i^+=0\), and the core issue is that the product sulfur atom is not produced in a single undifferentiated state: the branching into S(\(^3P\)) and S(\(^1D\)) is channel dependent [2212.13538].

The calculation is performed in multichannel quantum defect theory. The model includes both the previously studied \(^{2}\Pi\) manifold and the newly added \(^{4}\Pi\) manifold, uses 33 ionization channels in total—18 vibrational levels on the SH\(^+\) \(X\,^3\Sigma^-\) core and 15 on the excited \(A\,^3\Pi\) core—retains only the incident \(p\)-wave electron, and treats the reaction matrix to second order. The DR cross section contains both direct capture into dissociative resonant states and indirect capture into Rydberg states followed by predissociation. Rate coefficients are obtained from the Maxwellian average
\[
k(T)=\frac{8\pi}{\sqrt{m}(2\pi k_B T)^{3/2}} \int_0^\infty \sigma(\varepsilon)\,\varepsilon\,e^{-\varepsilon/k_B T}\,d\varepsilon.
\]

The \(^{4}\Pi\) contribution is reported to be much smaller than the \(^{2}\Pi\) one: at low collision energies it is about \(3\)–\(4\) orders of magnitude below the \(^{2}\Pi\) contribution, and around \(100\) K the \(^{2}\Pi\) rate is near \(10^{-7}\,\mathrm{cm^3\,s^{-1}}\) whereas the \(^{4}\Pi\) rate is only a few \(\times10^{-11}\,\mathrm{cm^3\,s^{-1}}\). The paper attributes this to a less favorable crossing with the ion and relatively weak valence–Rydberg couplings in the Franck–Condon region.

The state-selective sulfur yield follows from the dissociative branches. In the initial “no-flux-redistribution” picture, S(\(^3P\)) comes from capture into \(D_2\,^2\Pi\) and \(D'_1\,^4\Pi\), while S(\(^1D\)) comes from the \(D_1\,^2\Pi\) path. Averaged over \(10\)–\(1000\) K, the raw yields are about \(24\%\) into S(\(^3P\)) and \(76\%\) into S(\(^1D\)). The thermal rate coefficient for S(\(^3P\)) production is about \(2\times10^{-8}\,\mathrm{cm^3\,s^{-1}}\) on average between \(10\) and \(200\) K, with a maximum of \(\sim4\times10^{-8}\,\mathrm{cm^3\,s^{-1}}\) near \(40\) K. For S(\(^1D\)) production, the rate decreases monotonically from about \(2\times10^{-7}\,\mathrm{cm^3\,s^{-1}}\) at \(10\) K to \(3\times10^{-8}\,\mathrm{cm^3\,s^{-1}}\) at \(1000\) K.

The decisive correction arises from the crossing of the two dissociative \(^{2}\Pi\) curves around
\[
R_0\sim8\,a_0.
\]
The crossing is treated with a Landau–Zener model,
\[
P_{12}(\varepsilon)=1-\exp\!\left(-\frac{2\pi V_{12}^2(R_0)}{\hbar a v}\right),
\]
with \(V_{12}(R_0)=0.0314\) Hartree. The redistributed cross sections are
\[
\sigma'_1(\varepsilon)=(1-P_{12})\sigma_1(\varepsilon)+P_{12}\sigma_2(\varepsilon),
\]
\[
\sigma'_2(\varepsilon)=(1-P_{12})\sigma_2(\varepsilon)+P_{12}\sigma_1(\varepsilon).
\]
The quoted probabilities are \(P_{12}=0.99403\) at \(0.001\) eV, \(0.99400\) at \(0.01\) eV, \(0.99373\) at \(0.1\) eV, and \(0.99246\) at \(0.5\) eV. Because \(P_{12}\approx0.994\), the branching ratios invert: for example, at \(0.001\) eV the yields change from \(n(^3P)=0.274\), \(n(^1D)=0.726\) without Landau–Zener redistribution to \(n(^3P)=0.724\), \(n(^1D)=0.276\) with it. After this treatment, the paper states that more than two-thirds of recombination events produce ground-state sulfur at very low collision energy, in good agreement with TSR storage-ring measurements.

For astrochemical use, the DR rate of \(v_i^+=0\) is fitted as
\[
k(T)=a_0\left(\frac{T}{300}\right)^{a_1}e^{-a_2/T},
\]
with \(a_0=6.53104\times10^{-8}\), \(a_1=-0.348405\), \(a_2=-2.59354\), and RMS \(0.00855\), where \(k\) is in \(\mathrm{cm^3\,s^{-1}}\). In this usage, “popping” thus refers neither to desorption nor to lattice displacement, but to the state-resolved emergence of sulfur atoms from a recombination event.

## 6. Relation to sulfur polyamorphism models and broader conceptual distinctions

A separate sulfur modeling framework addresses liquid polyamorphism through a maximum-valence approach rather than through any popping event. In this model, each particle has a limited number of covalent bonds, and the sulfur-like case is \(z=2\), with states \(S_0\), \(S_1\), and \(S_2\) corresponding to zero, one, and two bonds [2111.08109]. The interactions are: a van der Waals core-core attraction that generates the liquid–gas transition, a narrow shell-shell attractive well that represents reversible covalent bonding, and an additional attraction between maximally bonded \(S_2\) atoms when they are not chemically bonded to one another. The baseline parameters are \(w=0.4\sigma\), \(w_b=0.02\sigma\), \(\epsilon_b=\epsilon\), and for the extra sulfur-like term \(\epsilon_{22}=0.5\epsilon\), \(w_{22}=0.3\sigma\).

For \(z=2\), the reversible transitions are
\[
S_0 \leftrightarrow S_1 \leftrightarrow S_2,
\]
with explicit reactions \(S_0+S_0\to S_1+S_1\), \(S_0+S_1\to S_1+S_2\), and \(S_1+S_1\to S_2+S_2\), plus the reverse bond-breaking processes. The model yields both a low-pressure liquid–gas phase transition and a high-pressure liquid–liquid phase transition. The reported critical values are
\[
T_c^{\mathrm{LG}}=1.023,\quad P_c^{\mathrm{LG}}=0.0922,\quad \rho_c^{\mathrm{LG}}=0.35,
\]
and
\[
T_c^{\mathrm{LL}}=1.187,\quad P_c^{\mathrm{LL}}=2.28,\quad \rho_c^{\mathrm{LL}}=0.81.
\]
The LL coexistence line has a positive slope in the \(P\)–\(T\) plane, and the LLPT separates a monomer-rich LDL from a polymer-rich HDL. The paper states explicitly that without the \(S_2\)-\(S_2\) attraction no LLPT appears.

Its order-parameter analysis uses \(\phi_0\), \(\phi_1\), and \(\phi_2\), with \(\phi_0+\phi_1+\phi_2=1\), and identifies \(1-\phi_2\) as an appropriate order parameter for the LLPT coupled to polymerization because \(\phi_2\) is symmetric across the coexistence curve and centered near \(\phi_2=0.5\). This liquid-sulfur framework is mechanistically different from all three popping usages above. It concerns polymerization-driven phase separation, not grain desorption, defect relocation, or state-selective atomic branching.

A plausible implication is that sulfur-focused models often derive their explanatory power from a discrete local transformation that is tightly coupled to a macroscopic observable: mantle conversion and desorption in Orion KL, a defect-plus-charge configuration in MoS\(_2\), state-selective dissociation in SH\(^+\) recombination, and bond-counting-driven polymerization in liquid sulfur. That common structural motif does not make the mechanisms interchangeable. The phrase “sulphur atom popping model” therefore requires domain-specific interpretation, with the operative physics supplied by the surrounding literature rather than by the phrase itself.

Source: https://www.emergentmind.com/topics/sulphur-atom-popping-model