Inclusion of sulfur chemistry in a validated C/H/O/N chemical network: identification of key C/S coupling pathways

This study presents a validated C/H/O/N/S chemical network that integrates sulfur kinetics with carbon and nitrogen chemistry through combustion data and ab initio calculations, revealing that C/S coupling significantly alters exoplanet atmospheric abundances and observables, particularly by increasing predicted CS2 levels.

The Gist

Imagine the atmosphere of a distant planet as a giant, chaotic kitchen where chemical ingredients are constantly being chopped, mixed, heated, and burned. For years, astronomers have been trying to write the "recipe book" (a chemical network) to predict what these atmospheres look like.

However, until now, this recipe book had a major blind spot. It was great at describing the main ingredients like Carbon, Hydrogen, Oxygen, and Nitrogen (the "CHON" team), but it completely ignored Sulfur.

This paper is like a chef realizing, "Wait a minute! We forgot the sulfur!" and then rushing to the library to find the missing recipes, testing them in a lab, and finally adding them to the master cookbook.

Here is the breakdown of what they did, using some everyday analogies:

1. The Problem: The Missing Ingredient

Recently, the James Webb Space Telescope (JWST) looked at two planets, WASP-39 b and WASP-107 b, and found Sulfur Dioxide (SO₂). This was a big deal because it was the first time sulfur was spotted in an exoplanet's sky.

The problem? The existing computer models used to predict these atmospheres didn't know how to handle sulfur properly. They treated sulfur like a guest who just sits in the corner, ignoring how sulfur interacts with the other ingredients (like Carbon). It's like trying to bake a cake but ignoring how the eggs interact with the flour; the result won't taste right.

2. The Solution: Borrowing from Firefighters

The authors needed a better sulfur recipe, but there wasn't one for space. So, they looked at combustion science (the study of fire and burning).

• The Analogy: Think of a burning forest fire or a car engine. These environments are hot and messy, just like the atmospheres of "Hot Jupiters" (giant, hot planets). Scientists who study fires have spent decades figuring out exactly how sulfur behaves when it burns.
• The Move: The team took these "fire recipes" and adapted them for space. They didn't just copy-paste; they tested them against 1,600 different experiments (like burning hydrogen sulfide or methanethiol in a lab) to make sure the math worked.

3. The Big Discovery: The "Sulfur Bridge"

When they added sulfur to their model, they found something surprising: Sulfur doesn't just sit there; it shakes hands with Carbon.

• The Key Player: They discovered a tiny, unstable molecule called CH₂S (Methylenethio). Think of CH₂S as a bridge or a connector.
• The Effect: In the old models, Carbon and Sulfur were strangers. In the new model, CH₂S acts like a matchmaker, allowing them to react.
• The Result: This connection creates a massive amount of Carbon Disulfide (CS₂).
  • Imagine: You thought there were only a few cars in a parking lot. But once you realize there's a secret tunnel connecting the lot to a highway, you suddenly find thousands of cars. That's what happened to CS₂. The new model predicts there is way more CS₂ in these atmospheres than we thought—sometimes 100 to 1,000 times more!

4. Why Does This Matter? (The Spectral Fingerprint)

Astronomers don't see planets with their eyes; they look at the light passing through the planet's atmosphere (a "transit spectrum"). Every gas leaves a unique "fingerprint" on that light.

• The Old View: The models predicted a certain pattern of light based on the old, incomplete recipes.
• The New View: Because there is so much more CS₂ and other sulfur-carbon mixtures, the "fingerprint" changes.
  • For Hot Planets: The changes are subtle, mostly affecting the amount of Methane (CH₄).
  • For Cooler Planets: The changes are huge! The new model predicts a strong signal from CS₂ at specific wavelengths (like a loud note in a song). This matches recent, real-life observations of a planet called TOI-270 d, where astronomers actually did find hints of CS₂.

5. The Takeaway

This paper is a reminder that in the era of the James Webb Space Telescope, we can't just guess. We need validated, comprehensive recipes.

• The Lesson: If you want to understand the atmosphere of a planet, you can't ignore the "weird" elements like sulfur. They might be small, but they change the whole flavor of the dish.
• The Future: By using data from fire labs (combustion) to fix our space models, we can finally stop guessing and start knowing exactly what these alien worlds are made of.

In short: The authors took the "fire science" of sulfur, mixed it with their "space science" of carbon, and discovered that the two are best friends. This friendship creates a new chemical party (CS₂) that changes how we see and understand the atmospheres of distant worlds.

Technical Summary

1. Problem Statement

Recent detections of sulfur dioxide ($SO_2$) in the atmospheres of exoplanets WASP-39 b and WASP-107 b by the James Webb Space Telescope (JWST) have highlighted a critical gap in current atmospheric modeling.

• The Gap: Existing chemical networks used for exoplanets typically treat sulfur chemistry as a separate sub-mechanism added to a Carbon/Hydrogen/Oxygen/Nitrogen (C/H/O/N) network. This approach neglects the complex kinetic coupling between sulfur and carbon/nitrogen species.
• Data Scarcity: Reliable kinetic data for sulfur reactions, particularly those involving coupling with carbon (C/S) and nitrogen (N/S), is scarce in the literature. Most combustion networks for hetero-atomic compounds are incomplete or lack validation for the specific conditions (high temperature, hydrogen-dominated atmospheres) found in exoplanets.
• Consequence: Current models may produce inaccurate abundance profiles and synthetic spectra, leading to potential misinterpretations of observational data (e.g., metallicity, formation history).

2. Methodology

The authors developed a comprehensive, fully coupled C/H/O/N/S kinetic network using a dual approach combining combustion literature and ab initio calculations, followed by rigorous experimental validation.

• Network Construction:

  • Base: Built upon a previously validated C0-C2 C/H/O/N mechanism (Veillet et al. 2024).
  • Sulfur Integration: Merged the Stagni et al. (2022) H$_2$S mechanism with the C/H/O/N core.
  • C/S Coupling: Integrated mechanisms for Methanethiol (CH$_3$SH) and Carbon Disulfide (CS$_2$). Since existing models (e.g., Colom-Díaz et al. 2021) showed numerical instabilities or lacked key pathways, the authors developed a new detailed sub-mechanism for CH$_3$SH.
  • Theoretical Calculations: Used ab initio calculations (CCSD(T)-F12/CBS level) to determine rate coefficients for sensitive reactions and unknown thermochemistry, particularly for the CH$_3$SH pyrolysis pathways.
  • N/S Coupling: Included the N/S reaction subset from Cooper et al. (2022) to account for nitrogen-sulfur interactions.

• Validation Strategy:

  • The network was validated against 1,606 experimental data points from the literature covering combustion and pyrolysis of H$_2$S, CH$_3$SH, CS$_2$, and OCS.
  • Data sources included shock tubes, flow reactors, and jet-stirred reactors across a wide range of temperatures (500–2500 K) and pressures.
  • Thermal Decomposition of H$_2$S: The authors addressed a known conflict between theory and experiment regarding the $H_2S \rightarrow S + H_2$ pathway (which is spin-forbidden). They introduced an intermediate excited state species, $S(^1D)$, to reconcile the observed activation energy with the correct reverse reaction pathway.

• Exoplanet Application:

  • Models: Used the 1D kinetic model FRECKLL to simulate six exoplanets: GJ 436 b, GJ 1214 b, HD 189733 b, HD 209458 b, WASP-39 b, and WASP-107 b.
  • Comparison: Results were compared against:
    1. The VULCAN network (Tsai et al. 2022) used in previous JWST studies.
    2. The authors' previous C/H/O/N network (Veillet et al. 2024) without sulfur.
  • Observables: Synthetic transmission spectra were generated using TauREx 3.1 to assess the impact on observables.
  • Pathway Analysis: Used graph theory (NetworkX) to solve the minimum flow problem and identify dominant formation and destruction pathways.

3. Key Contributions

• First Validated C/H/O/N/S Network: Creation of the first sulfur kinetic network for exoplanets that is extensively validated against combustion and pyrolysis experiments, ensuring reliability.
• Identification of Key Coupling Species:
  • CH$_2$S (Thioformaldehyde): Identified as the critical "bridge" species coupling carbon and sulfur chemistry. It is missing in many current exoplanet networks but is essential for accurate modeling.
  • CS$_2$ (Carbon Disulfide): Found to be a major product of C/S coupling, with abundances significantly higher than previously predicted.
• New Reaction Pathways:
  • Discovered a dominant pathway for CS$_2$ formation via $CH_3 \rightarrow CH_2S \rightarrow CHS \rightarrow CS \rightarrow CS_2$.
  • Identified a bypass mechanism for ethylene ($C_2H_4$) formation: $CH_2S + CH_3 \rightarrow C_2H_4 + SH$, which skips the rate-limiting $CH_3 + CH_3$ recombination step.
  • Mapped complex N/S coupling pathways involving species like NS, SNO, and NO.

4. Key Results

• Abundance Profiles:

  • CS$_2$ Abundance: The new network predicts CS$_2$ abundances 2 to 3 orders of magnitude higher than the Tsai 2022 network for warm Neptunes (GJ 436 b, GJ 1214 b) and hot Jupiters. In some cases, CS$_2$ reaches ~1000 ppm.
  • Methane ($CH_4$): Sulfur coupling significantly reduces $CH_4$ abundance in certain layers (e.g., HD 189733 b) due to oxidation pathways involving CS and SO.
  • Hydrogen Cyanide ($HCN$) & Ammonia ($NH_3$): The coupling alters $HCN$ and $NH_3$ abundances. For HD 189733 b, $HCN$ abundance drops by ~2 orders of magnitude due to the depletion of its precursor $CH_4$.
  • Acetylene ($C_2H_2$) & Ethylene ($C_2H_4$): Abundances of these species increase significantly (up to 2 orders of magnitude) due to the new sulfur-mediated formation pathways.

• Spectral Impact:

  • CS$_2$ Features: Strong spectral features of CS$_2$ appear at 11 $\mu$m and 25 $\mu$m. These are particularly prominent in cooler planets like GJ 1214 b and WASP-107 b.
  • WASP-107 b: The new model predicts a massive reduction in $CH_4$ features (amplitude change of ~2000 ppm at 3.3 $\mu$m) and a significant shift in $NH_3$ features (disappearance of the 20 $\mu$m feature).
  • WASP-39 b: While $SO_2$ profiles remain similar to previous models, the inclusion of C/S coupling reveals subtle differences in $CH_4$ and $NH_3$ that affect the continuum.

• Validation of Theory: The new network successfully reproduces experimental CS$_2$ yields from CH$_3$SH pyrolysis, a feat the previous Tsai 2022 network could not achieve, confirming the necessity of the CH$_2S$ intermediate.

5. Significance

• JWST Era Relevance: As JWST detects more sulfur-bearing species (SO$_2$, and potentially CS$_2$ in TOI-270 d), accurate kinetic networks are essential to avoid "blind spots" in atmospheric retrieval.
• Chemical Markers: The study suggests CS$_2$ could serve as a photochemical marker for cooler exoplanets, similar to how SO$_2$ marks photochemistry in hot Jupiters.
• Methodological Shift: The paper demonstrates that combustion and pyrolysis data are invaluable resources for exoplanet atmospheric modeling. By leveraging these terrestrial datasets, researchers can build more robust, validated networks that capture complex hetero-atomic couplings previously ignored.
• Future Directions: The authors call for further experimental and theoretical work on CS$_2$ oxidation in hydrogen-rich environments and N/S coupling to fully resolve remaining uncertainties in exoplanet chemistry.