{"id":"d2833254-9b18-4c8c-986d-c1f7d4eca769","arxiv_id":"2509.08349","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"Multi-epoch Na, K, and SO2 variations around WASP-39b are attributed to a hypothetical tidally heated volcanic exomoon rather than to the planet's atmosphere alone.","lead":"This paper compares sodium, potassium, and sulfur dioxide absorption seen around WASP-39b at different times and reports large changes between ground-based and JWST observations. It argues that a small, tidally heated volcanic moon venting gas into a cloud or torus could explain those changes, though the case depends on many unconstrained assumptions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Cross-instrument systematics may fully explain the claimed Na/K variability; a common-pipeline re-analysis is needed before the exomoon interpretation can be supported.","rationale":"The reader's weakest assumption correctly identifies the cross-instrument systematics as the critical premise. My stress-test agrees and specifies a concrete test: re-analysis with a common pipeline. The internal inconsistency in SO2 mass-loss rates is real but secondary; the variability claim is the foundation of the paper, and if that variability is not established, the exomoon scenario loses its observational motivation. Since the reader's CONDITIONAL verdict already incorporates this uncertainty, I do not propose changing the verdict. The paper should explicitly state that the variability is contingent on a common systematics analysis, and the mass-loss inconsistency in Table 2 vs §4.1 should be resolved before the tidal-consistency claim is made.","tokens_in":14124,"tokens_out":9875,"duration_ms":97777,"concrete_test":"Download the publicly available reduced spectra for WASP-39 b from Nikolov et al. (2016, VLT/FORS2), Fischer et al. (2016, HST/STIS), Rustamkulov et al. (2023, JWST/NIRSpec), and Feinstein et al. (2023, JWST/NIRISS). Re-extract Na I and K I equivalent widths using one pipeline: (1) fit a common stellar+planet model to line-free continuum regions within ±50 Å of each D line, (2) convolve the same intrinsic line model with each instrument's LSF, (3) compute Wλ with the same integration limits, and (4) include instrument-specific systematics (wavelength solution, flat-field, normalization) in the error budget. If the resulting logN values still differ by >5x between epochs with non-overlapping uncertainties, the variability claim is robust; if they converge or the errors expand to overlap, the exomoon interpretation is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Na, K, and SO2 column densities vary by >order of magnitude across epochs rests on comparing VLT/FORS2, HST/STIS, JWST/NIRSpec, and JWST/NIRISS spectra without a common systematics budget (§2, Tables 1–2, Fig. 1). The paper normalizes all transit depths to a constant c=0.0210 and allows a free continuum normalization per spectrum in the MCMC retrievals (§2.2). Because the alkali lines are unresolved at all resolutions, the measured equivalent widths (Eq. 4) are highly sensitive to continuum placement, instrument line-spread function, and wavelength calibration. For K, the detections in VLT (logN=11.42) and NIRISS (10.61) are separated by ~6x, while HST and NIRSpec yield only upper limits (<10.43, <10.47) – differences that could easily arise from normalization or telluric/systematic offsets. No cross-instrument calibration check is provided to show these variations are astrophysical; indeed, §4.2 admits that 'characterizing possible instrumental and stellar variability would benefit from dedicated repeat observations.' If the apparent variability is an artifact, the exomoon hypothesis loses its observational motivation. Separately, the SO2 mass-loss rates in Table 2 (e.g., (0.6–40)×10^11 kg/s for G395H) are inconsistent with the torus-limited range quoted in §4.1 (10^8.2–9.8 kg/s), which also undermines the tidal-consistency claim, but the cross-instrument issue is more fundamental because it attacks the reality of the variability itself.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes archival transmission spectra of the hot Saturn WASP-39b in the optical (Na I, K I) and infrared (SO2, CO2), reporting epoch-to-epoch variability of more than an order of magnitude in Na, K, and SO2 column densities. The authors interpret this variability as evidence for a tidally heated volcanic exomoon ('WASP-39b I') venting gas into a cloud or torus around the planet. They use analytical equivalent-width estimates, MCMC retrievals of molecular opacities, and the open-source Monte-Carlo codes prometheus and serpens to simulate sputtered clouds and tori. They conclude that the simulations can reproduce the observed line-of-sight column-density variations and that the inferred SO2 mass-loss rate is consistent with tidal-heating predictions. They also note that high-resolution alkali Doppler observations are needed to constrain a putative orbit.","tokens_in":14651,"tokens_out":5680,"duration_ms":63067,"significance":"If the variability is real and the exomoon interpretation is correct, this would be one of the first indirect detections of an exomoon through its volcanic gas output, with major implications for satellite formation and tidal physics. The paper builds on a series of prior works by the same group and makes use of public codes and JWST archival data. It explicitly acknowledges some limitations, such as the need for repeat observations. However, the current evidence is not yet convincing because the central variability claim rests on heterogeneous datasets and the model 'reproduction' is partly normalized to the observations. The work is important as a hypothesis-generating study but needs substantial additional analysis before the exomoon claim can be supported.","major_comments":[{"comment":"The central variability claim rests on comparing VLT/FORS2, HST/STIS, JWST/NIRISS, JWST/NIRSpec, G395H, and MIRI data without a common systematics budget. The MCMC retrieval in §2.2 treats the continuum normalization as a free parameter per dataset, and Eqn. (4) for unresolved alkali equivalent widths is highly sensitive to continuum placement, line-spread function, and wavelength calibration. The K detections/non-detections differ by factors of several between instruments (Table 2: VLT logN=11.42, NIRISS 10.61, HST <10.43, NIRSpec <10.47), and these differences could plausibly be instrumental. A homogeneous re-analysis with a single pipeline, or at least a quantitative cross-instrument systematics budget, is required before the astrophysical variability that motivates the exomoon hypothesis can be considered established.","section":"§2.2"},{"comment":"The simulations in Figure 4 are normalized by maintaining a constant total number of atoms N, where N is set from the observed column densities via Eqn. (2). Thus the models are not independent predictions of the absolute density; they are scaled to the data. The resulting phase curves can match the observed points with a suitable choice of geometry and satellite phase at each epoch. To claim that the simulations 'reproduce' the observed variations, the paper should either (a) fix the mass-loss rate from tidal theory (e.g., ~10^7.5 kg/s for SO2) and predict the observed column density, or (b) clearly state that N is an input and provide a quantitative goodness-of-fit metric (e.g., chi-square or likelihood) that accounts for the free normalization, phase, and geometry choices.","section":"§3"},{"comment":"There is an internal inconsistency in the reported SO2 mass-loss rates. Table 2 lists SO2 Mdot values of (0.6–40)×10^11 kg/s (G395H) and (0.15–10)×10^10 kg/s (MIRI), whereas §4.1 quotes a range of 10^8.2–9.8 kg/s for SO2 and Figure 4iii shows a tidal prediction of 10^7.5±1 kg/s. These differ by orders of magnitude. Because the abstract and conclusions claim that the 'estimated SO2 flux is consistent with tidal gravitation predictions,' the definition of which Mdot is tabulated (gamma-limited vs torus-limited) and the numerical reconciliation must be provided. As written, the reader cannot assess the consistency claim.","section":"Table 2"},{"comment":"The geometry is selected per species after the fact: Na is claimed to be a cloud, K a cloud, and SO2 a torus. The satellite orbital phase at each epoch appears to be free, and the satellite mass can be varied (Io-mass to Earth-mass, Fig. 4iii). With such flexibility, the model can accommodate nearly any observed sequence of column densities. Please provide a self-consistent model in which a single satellite mass, orbit, and geometry (or a small discrete set) is fit simultaneously to all species and epochs, reporting the best-fit parameters and uncertainties. Alternatively, state explicitly how many degrees of freedom the model has and why the current post-hoc choices are not overfitting.","section":"§4"}],"minor_comments":[{"comment":"The phrase 'Roche limit interior to the planetary photosphere' is confusing. Does it mean the Roche limit lies at an orbital period shorter than 8 h, and thus inside the planet? Clarify the statement; Figure 5 seems to show a Roche limit outside the photosphere.","section":"Abstract"},{"comment":"Upper limits are written as '10 10.43', which is easily misread. Use standard notation such as logN < 10.43 or an arrow symbol.","section":"Table 2"},{"comment":"In Figure 4iii, the label 'M x 100' is unclear; specify which model is multiplied by 100 and whether it is a mass-loss rate or column density. Also define M_tor, M_cloud, and M_tides in the caption.","section":"Fig. 4"},{"comment":"The symbol N is used both for the total number of particles (Eqn. 2) and for column density (Eqn. 3). This is a common source of confusion; please use N_total or similar in Eqn. (2).","section":"Eqn. (2)"},{"comment":"The statement that the Na/SO2 ratio is 'far smaller than Io's' would benefit from a quantitative value. Currently the reader cannot assess the claimed factor.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper presents an intriguing hypothesis and uses appropriate public data and tools, but the evidence as it stands is too fragile: the variability could be instrumental, the model normalization makes the 'reproduction' partly circular, and the SO2 flux values are internally inconsistent. A major revision that (i) re-analyzes the data homogeneously or provides a systematics budget, (ii) turns the model into a forward prediction with fixed tidal input, and (iii) reconciles the SO2 mass-loss rates would be needed before this can be considered for publication. The topic is well within the scope of MNRAS and the paper could become a useful contribution after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper argues that order-of-magnitude variability in Na, K, and SO2 around WASP-39b across six years is caused by a tidally heated volcanic exomoon, using Prometheus and Serpens simulations. The claim is plausible but not yet supported. The variability is the load-bearing observation, and it rests entirely on comparing VLT/FORS2, HST/STIS, and two JWST instruments without a common systematics budget. The alkali lines are unresolved and the measured equivalent widths depend on continuum placement and line-spread function; the paper does not show those differences are astrophysical. The paper itself admits that \"characterizing possible instrumental and stellar variability would benefit from dedicated repeat observations,\" which undercuts the central inference.\n\nWhat is genuinely new is the multi-epoch compilation itself and the target-specific phase-curve simulations. The optically-thin retrievals are clearly described, the CO2 stability across epochs is a useful internal control, and the code is public. The paper is honest about calling the exomoon \"putative\" and pointing to the Doppler observations that would actually test it.\n\nThe soft spots are real. The simulations are normalized to the observed column densities (Eqn. 2), so reproducing those values is partly by construction. Satellite mass, radius, orbit, phase, and geometry are free, so the fit is not constrained enough to be evidential. And there is an internal inconsistency: Table 2 reports SO2 mass-loss rates around 10^11 kg/s for the G395H epoch, while Section 4.1 quotes a torus-limited range of 10^8.2–9.8 kg/s. That needs to be reconciled before the \"consistent with tidal predictions\" claim can be taken seriously.\n\nIf the variability evaporates under a common-pipeline re-analysis, the exomoon hypothesis is unmotivated. But the paper does not pretend to have a detection; it offers a scenario and a path to test it. I would send this to peer review, with a referee asked specifically to check the cross-instrument normalization and the mass-loss arithmetic. It belongs in the literature as a hypothesis-generating paper, and it will be a useful reference for the Doppler follow-up that is now the obvious next step.","headline":"A speculative but well-framed exomoon scenario that hinges on cross-instrument variability the paper does not yet prove; worth refereeing mainly to force a common-pipeline check.","tokens_in":15218,"tokens_out":2228,"would_cite":true,"duration_ms":25233,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A tidally heated volcanic satellite orbiting WASP-39 b can reproduce the observed order-of-magnitude swings in neutral sodium, potassium, and sulfur dioxide column densities.","keywords":["exomoons","tidal heating","volcanic satellites","WASP-39 b","transmission spectroscopy","sodium and potassium lines","sulfur dioxide","hot Saturn atmospheres"],"falsifier":"Take a single high-resolution spectrograph and observe WASP-39 b over many consecutive transits at R around 100,000. If the Na and K lines always remain at the planetary rest frame with constant equivalent widths, or if the apparent column-density swings disappear when all data are re-reduced with a common continuum normalization, the volcanic-satellite explanation fails. A positive test would be a periodic Doppler shift and phase-locked column-density modulation with a period shorter than 15.3 hours.","tokens_in":14063,"feed_emoji":"🌋","tokens_out":5942,"duration_ms":64334,"temperature":0.7,"pith_summary":"This paper argues that the irregular, order-of-magnitude swings in neutral sodium, potassium, and sulfur dioxide seen in WASP-39 b's transit spectra are not instrumental artifacts or a static planetary atmosphere, but the signature of an unseen, tidally heated volcanic moon venting gas into a cloud or torus around the planet. Combining ground-based VLT, HST, and JWST observations from several epochs, the authors measure line-of-sight column densities of Na I, K I, SO2, and CO2 in the optically thin limit and find variability far exceeding CO2's stability. Monte Carlo simulations that place an outgassing satellite in orbit reproduce both the magnitude and the phase-dependent behavior of the observed column densities, with implied mass-loss rates consistent with three-body tidal-heating predictions. If correct, the result turns apparent inconsistencies between instruments into evidence for an exomoon, and makes WASP-39 b a test bed for detecting volcanic moons around hot Jupiters.","feed_headline":"Volcanic moon drives WASP-39 b's shifting sodium and sulfur gas","feed_subtitle":"Gas vented by a tidally heated satellite matches the decade of sodium, potassium, and SO2 swings seen in transit spectra.","key_machinery":"The argument rests on the optically thin absorption relation transit depth ≈ N σ(λ), which converts measured line depths into line-of-sight column densities, plus the mass-loss relation N = (Mdot/m) τ, where the lifetime τ is set either by photoionization for a cloud geometry or by charge-exchange in a plasma torus. On top of this framework, the paper runs its own three-dimensional Monte Carlo and radiative-transfer simulations of gas evaporating from a satellite under stellar radiation pressure; the key comparison is between predicted phase-dependent column-density curves and the epoch-stamped data points. The central identity is that column density scales directly with mass-loss rate, so t","core_discovery":"The paper's central claim is that a tidally heated volcanic satellite, an exo-Io, orbiting WASP-39 b can account for the more-than-order-of-magnitude variability in Na I, K I, and SO2 columns measured across 2013-2023. Using optically thin transmission spectroscopy, the authors compute Na and K column densities from alkali equivalent widths and SO2/CO2 columns from NIR and MIR spectra, then compare these with three-dimensional test-particle simulations of gas sputtered from a satellite into either a localized cloud or a toroidal structure. The simulations reproduce the observed line-of-sight column density variations, and the estimated SO2 flux is consistent with tidal-gravitation prediction","pith_inferences":["A natural next test is to stack multiple transits at the same satellite phase: if column densities repeat with orbital phase, that would be a direct orbital fingerprint separate from any planetary atmosphere signal.","The torus-versus-cloud distinction suggests a population-level ordering: short-lived species like potassium should vanish when the moon is occulted, while longer-lived species like SO2 should vary more smoothly; this prediction could be checked in other systems.","If confirmed, WASP-39 b would offer a gas-based exomoon detection path that is cheaper and faster than light-curve searches, since it uses existing transmission spectra rather than dedicated moon transits.","The inferred mass loss places the satellite near the Roche limit, so monitoring over years could reveal orbital decay or the early stages of ring formation."],"forward_implications":["If the variability is real, single-epoch retrievals of SO2 in hot-Saturn atmospheres may be contaminated by an exogenic source, so abundance estimates should account for possible satellite venting.","High-resolution alkali observations of WASP-39 b across multiple transits should reveal a Doppler-shifting, phase-dependent signal that pinpoints the satellite's orbit, with a period shorter than about 15.3 hours.","The Na/SO2 ratio becomes a diagnostic: values far below Io's point to a hotter, more efficiently stripped volcanic source, distinguishing satellite venting from photochemical SO2 production.","The inferred mass-loss rates imply a dusty, volcanically sourced component that JWST mid-infrared observations could detect directly.","Other hot Jupiters with unexplained alkali variability could show similar cloud or torus signatures if volcanic moons are common around such planets."],"supporting_citations":[{"why":"Supplies the three-body tidal-heating formalism and predicted mass-loss rates that the paper compares with its inferred SO2 flux.","marker":"Oza et al. (2019)"},{"why":"Provides the evaporative transmission spectroscopy forward model used to convert transit depths into column densities and to simulate cloud and torus geometries.","marker":"Gebek & Oza (2020)"},{"why":"Supplies the three-dimensional Monte Carlo code and stellar radiation-pressure treatment used for the phase-dependent density predictions.","marker":"Meyer zu Westram et al. (2024)"},{"why":"Supplies VLT/FORS2 visible alkali transit spectra that anchor the high column-density epochs.","marker":"Nikolov et al. (2016)"},{"why":"Supplies HST/STIS alkali spectra used for a second optical epoch.","marker":"Fischer et al. (2016)"},{"why":"Supplies JWST/NIRSpec spectra from which Na, K, SO2, and CO2 column densities and their epoch variability are retrieved.","marker":"Rustamkulov et al. (2023)"},{"why":"Supplies the JWST/NIRISS potassium epoch used to show potassium transitioning from detection to non-detection.","marker":"Feinstein et al. (2023)"},{"why":"Supplies JWST/MIRI SO2 and CO2 data and the atmospheric retrieval baseline that the volcanic hypothesis challenges.","marker":"Powell et al. (2024)"},{"why":"Gives the dynamical stability boundary used to bound the putative satellite's orbit at 0.41 Hill radius and under 15.3 hours.","marker":"Kisare & Fabrycky (2024)"}],"fun_headline_variants":["Volcanic moon likely drives WASP-39 b's gas variability","Exo-Io volcano may explain shifting Na, K, SO2","Tidally heated volcano moon fits WASP-39b gas swings","Exo-Io's vented gas matches hot Saturn's sodium shifts","Volcanic satellite model mimics WASP-39b's alkali swings"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that the epoch-to-epoch differences in Na, K, and SO2 column densities are real astrophysical variations rather than artifacts of comparing different instruments, normalizations, spectral ranges, and retrieval assumptions.","fun_headline_variants_meta":{"raw":{"variants":["Volcanic moon likely drives WASP-39 b's gas variability","Exo-Io volcano may explain shifting Na, K, SO2","Tidally heated volcano moon fits WASP-39b gas swings","Exo-Io's vented gas matches hot Saturn's sodium shifts","Volcanic satellite model mimics WASP-39b's alkali swings"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000423,"raw_usage":{"total_tokens":2013,"prompt_tokens":755,"completion_tokens":1258,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":499,"completion_tokens_details":{"reasoning_tokens":1174}},"tokens_in":499,"tokens_out":1258,"duration_ms":10713,"temperature":1.0,"reasoning_tokens":1174,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T20:40:57.840393+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a single high-resolution spectrograph and observe WASP-39 b over many consecutive transits at R around 100,000. If the Na and K lines always remain at the planetary rest frame with constant equivalent widths, or if the apparent column-density swings disappear when all data are re-reduced with a common continuum normalization, the volcanic-satellite explanation fails. A positive test would be a periodic Doppler shift and phase-locked column-density modulation with a period shorter than 15.3 hours.","supporting_citations":[],"review_version":1}