A new look into the atmospheric composition of WASP-39 b
Pith reviewed 2026-05-22 20:17 UTC · model grok-4.3
The pith
Silicon monoxide rather than sulfur dioxide can explain the 4.1 micron absorption in WASP-39 b while keeping the atmosphere in chemical equilibrium.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By iterating between equilibrium chemistry models and blind retrievals the analysis identifies a silicon-dominated chemistry in which SiO produces the 4.1 μm feature instead of SO2. The retrieved gases remain consistent with an atmosphere in thermochemical equilibrium whose temperature-pressure structure is fixed by the H2O and CO2 absorption bands, while silicate clouds and hazes account for the mid-infrared continuum observed by MIRI.
What carries the argument
Iterative loop of ab initio equilibrium chemistry models and blind retrievals that pre-selects species and post-reconstructs the solution to enforce physical consistency.
If this is right
- The atmosphere remains consistent with chemical equilibrium once the temperature-pressure profile is anchored by H2O and CO2.
- Silicate clouds and hazes reproduce the 5-12 μm MIRI window without additional ad-hoc opacity sources.
- Na and K are required alongside the molecular species to match the optical and near-infrared data.
- The same iterative retrieval strategy can be applied to other JWST targets to test whether silicon chemistry has been overlooked.
Where Pith is reading between the lines
- If confirmed, the result would imply that silicon-bearing gases deserve systematic inclusion in retrieval grids for hot Jupiters observed at similar wavelengths.
- The method highlights how pre-selection of chemically plausible species before retrieval can reduce the number of free parameters and limit unphysical solutions.
- Future instruments with resolving power high enough to separate SiO and SO2 line lists could directly test the assignment without relying on the iterative framework.
Load-bearing premise
The iterative cycle between equilibrium chemistry and blind retrievals converges on a physically correct solution without introducing biases from the order of steps or from later adjustments to fit the data.
What would settle it
A higher-resolution spectrum around 4.1 μm that shows the line positions and widths expected for SO2 rather than SiO, or that detects strong SO2 bands at other wavelengths where the current model predicts none.
read the original abstract
Being one of the first exoplanets observed by the James Webb Space Telescope, WASP-39 b has become an iconic target and many transit spectra recorded with different instruments (NIRISS, NIRCAM, NIRSpec G395H, NIRSpec PRISM and MIRI) are currently available, allowing in-depth studies of its atmosphere. We present here a novel approach to interpret WASP-39 b's transit spectroscopic data, consisting of a multi-step process where ab initio equilibrium chemistry models and blind retrievals are used iteratively to find physically robust, optimal solutions. Following this approach, we have identified a new scenario to explain WASP-39 b's atmospheric composition, in which silicon-based chemistry plays a major role. In this scenario, SiO may explain the spectral absorption at 4.1 $\mu$m, currently interpreted as being due to SO$_2$. SiO and the other gas species identified by the retrieval models, i.e. H$_2$O, CO$_2$, Na and K, are consistent with an atmosphere in chemical equilibrium with a temperature-pressure profile constrained by H$_2$O and CO$_2$ absorption bands. In addition, silicate clouds and hazes can produce the spectral features observed by MIRI in the spectral window 5-12 $\mu$m. While we advocate the need for more data, possibly at higher spectral resolution, to confirm our results for WASP-39 b's atmospheric composition, we highlight a refined atmospheric retrieval strategy with pre-selection and post-reconstruction to guide the next generation of transit spectroscopy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a multi-step iterative procedure that alternates ab initio equilibrium chemistry calculations with blind atmospheric retrievals to re-analyze the combined JWST transmission spectrum of WASP-39 b (NIRISS, NIRCam, NIRSpec G395H/PRISM, MIRI). The central claim is that SiO, rather than SO2, produces the 4.1 μm absorption feature; the retrieved mixing ratios of H2O, CO2, Na and K are consistent with thermochemical equilibrium along a T–P profile anchored by the H2O and CO2 bands; and silicate clouds plus hazes reproduce the 5–12 μm MIRI continuum and features.
Significance. If the iterative scheme can be shown to be free of data-driven bias, the result would provide a concrete alternative interpretation of an iconic JWST target and illustrate how equilibrium-chemistry priors can be folded into retrievals without circularity. The work also supplies a practical template for the next generation of retrieval pipelines that combine forward-model grids with data-driven exploration.
major comments (2)
- [§3.2] §3.2 (Iterative retrieval procedure): The manuscript asserts that the multi-step loop of pre-selection, equilibrium-chemistry reconstruction and post-reconstruction yields “physically robust, optimal solutions,” yet no quantitative control is presented (e.g., evidence ratio between the final SiO model and an otherwise identical SO2 model, or a cross-validation test that holds out the 4.1 μm region). Without such a metric it remains unclear whether the preference for SiO is an independent outcome or an artifact of the iteration itself.
- [§4.1, Figure 5] §4.1 and Figure 5: The claim that the retrieved T–P profile is “constrained by H2O and CO2 absorption bands” and therefore independent of the SiO choice is not demonstrated by a direct comparison; the posterior T–P contours shown appear to be derived from the full data set that already includes the 4.1 μm window. A dedicated retrieval run that masks the 4.1 μm region and then predicts the SiO abundance would test whether the equilibrium consistency is truly prior-driven.
minor comments (2)
- [Figure 3] The caption of Figure 3 should explicitly state the wavelength range and instrument combination used for each panel to avoid ambiguity when comparing to earlier SO2-only models.
- [Table 2] A short table listing the prior ranges adopted for the SiO, SO2 and cloud parameters in the final retrieval stage would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive report. We address each of the major comments below. We believe the suggested additions will strengthen the paper and have incorporated revisions to provide the requested quantitative controls and tests.
read point-by-point responses
-
Referee: [§3.2] §3.2 (Iterative retrieval procedure): The manuscript asserts that the multi-step loop of pre-selection, equilibrium-chemistry reconstruction and post-reconstruction yields “physically robust, optimal solutions,” yet no quantitative control is presented (e.g., evidence ratio between the final SiO model and an otherwise identical SO2 model, or a cross-validation test that holds out the 4.1 μm region). Without such a metric it remains unclear whether the preference for SiO is an independent outcome or an artifact of the iteration itself.
Authors: We agree that providing a quantitative metric would strengthen the claim of robustness. The iterative approach was intended to use equilibrium chemistry to guide but not dictate the retrieval, with blind retrievals allowing data to drive the initial selection. However, we recognize the value of an explicit comparison. In the revised manuscript, we will add a calculation of the Bayesian evidence ratio between the SiO-inclusive model and an SO2 model, as well as discuss the results of a cross-validation test holding out the 4.1 μm region. This will demonstrate that the preference is not an artifact. revision: yes
-
Referee: [§4.1, Figure 5] §4.1 and Figure 5: The claim that the retrieved T–P profile is “constrained by H2O and CO2 absorption bands” and therefore independent of the SiO choice is not demonstrated by a direct comparison; the posterior T–P contours shown appear to be derived from the full data set that already includes the 4.1 μm window. A dedicated retrieval run that masks the 4.1 μm region and then predicts the SiO abundance would test whether the equilibrium consistency is truly prior-driven.
Authors: The referee correctly notes that the presented T-P posterior uses the full dataset. While the H2O and CO2 bands are the dominant constraints on temperature across a broad wavelength range, a masked retrieval would provide clearer evidence. We will conduct and include results from a retrieval excluding the 4.1 μm region to show that the T-P profile remains consistent and that the equilibrium chemistry still supports the SiO abundance, confirming the independence from that specific feature. revision: yes
Circularity Check
Iterative equilibrium-chemistry + retrieval loop makes consistency claim non-independent
specific steps
-
fitted input called prediction
[Abstract]
"a multi-step process where ab initio equilibrium chemistry models and blind retrievals are used iteratively to find physically robust, optimal solutions. Following this approach, we have identified a new scenario to explain WASP-39 b's atmospheric composition, in which silicon-based chemistry plays a major role. In this scenario, SiO may explain the spectral absorption at 4.1 μm, currently interpreted as being due to SO2. SiO and the other gas species identified by the retrieval models, i.e. H2O, CO2, Na and K, are consistent with an atmosphere in chemical equilibrium with a temperature-press"
The retrievals fit abundances and T-P parameters directly to the spectral data; the equilibrium chemistry models are incorporated inside the same iterative loop. Declaring that the fitted species 'are consistent with an atmosphere in chemical equilibrium' therefore follows by construction from the procedure that already enforces or selects for equilibrium solutions, rather than constituting an independent prediction or test.
full rationale
The paper describes a multi-step iterative procedure that interleaves ab initio equilibrium chemistry models with blind retrievals on the same transit spectra to arrive at 'physically robust, optimal solutions.' The central claims—that SiO (not SO2) accounts for the 4.1 μm feature and that the retrieved gases plus T-P profile are consistent with chemical equilibrium—are presented as outputs of this loop. Because the equilibrium models are used to guide or constrain the retrievals that are then checked for consistency with equilibrium, the consistency statement reduces to a property of the fitting strategy rather than an independent confirmation. No quantitative control (evidence ratio, cross-validation, or hold-out test) is quoted to show the iteration does not preferentially favor SiO. This matches the 'fitted input called prediction' pattern at the level of the method itself. The remainder of the analysis (cloud/haze modeling, call for higher-resolution data) is independent and does not increase the score.
Axiom & Free-Parameter Ledger
free parameters (3)
- molecular abundances
- T-P profile parameters
- cloud/haze parameters
axioms (2)
- domain assumption The atmosphere is in chemical equilibrium
- ad hoc to paper SiO produces the 4.1 μm absorption feature
Forward citations
Cited by 2 Pith papers
-
Magnetic field strengths of hot giant exoplanets consistent with Solar System values
Wind speed measurements in seven ultra-hot Jupiters decrease with temperature, consistent with magnetic drag and implying magnetic field strengths of a few gauss.
-
Effect of tidal gravity and planetary rotation on the retrieved atmospheric abundances of close-in exoplanets
Including tidal and centrifugal gravity corrections increases retrieved molecular abundances in transmission spectra retrievals for WASP-12b (HST) and WASP-39b (JWST).
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.