Pith. sign in

REVIEW 3 major objections 4 minor 10 references

STCTM: a forward modeling and retrieval framework for stellar contamination and stellar spectra

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper presents STCTM, a Bayesian framework that models how star spots and faculae contaminate exoplanet transmission spectra, letting observers test whether an apparent atmospheric signal could actually be stellar contamination.

desk verdict Useful and already-used tool paper that combines TLSE and stellar retrievals, but it lacks an in-paper validation benchmark and slightly overstates the gap against POSEIDON and StellarFit. read the letter →

arxiv 2508.19297 v1 pith:NCLT3WRW submitted 2025-08-25 astro-ph.IM astro-ph.EPastro-ph.SR

Caroline Piaulet-Ghorayeb This is my paper · ORCID
classification astro-ph.IMastro-ph.EPastro-ph.SR
keywords transitlightsourceeffectstellarcontaminationexoplanettransmissionspectroscopyMdwarfstarsstarspotsfaculaeBayesianretrievalMCMC
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper introduces STCTM, a flexible, open-source Bayesian retrieval framework for the transit light source effect (TLSE). The TLSE is the wavelength-dependent imprint that unocculted stellar spots and faculae leave on a transmission spectrum, and it can mimic or mask planetary atmospheric features, especially for small planets around M dwarfs. The paper argues that the community needs a tool that can model this contamination generically, fit it to any transmission spectrum, and infer the stellar surface parameters that could explain the data without invoking a planetary atmosphere. STCTM also includes an exotune sub-module that retrieves stellar parameters directly from out-of-transit spectra, providing data-driven priors on how much contamination a host star can inject. If the framework works as claimed, it gives observers a practical way to test the robustness of exoplanet atmosphere detections.

What carries the argument

The load-bearing object is the TLSE forward model: a parameterized stellar surface made of spots and faculae with covering fractions and temperature contrasts, combined with precomputed synthetic spectra from model grids such as PHOENIX or SPHINX. The model produces a wavelength-dependent contamination spectrum that MCMC compares to the observed transmission spectrum. Exotune applies the same surface parameterization to out-of-transit spectra. The framework's genericity comes from allowing any user-specified stellar model grid and input data format.

What would settle it

Take a star with an independently measured surface map, for example from Doppler imaging, and a JWST transmission spectrum of one of its planets; run STCTM and compare the retrieved spot covering fractions and temperature contrasts to the mapped values. A systematic mismatch would show that the forward model's simplified surface representation is insufficient for real stellar surfaces.

Watch

Extended reading notes

Core claim

The central claim is that the transit light source effect can be forward-modeled and retrieved with a general, user-tunable Bayesian framework. STCTM parameterizes the stellar surface using covering fractions and temperature contrasts for spots and faculae, computes the resulting contamination spectrum from user-supplied synthetic stellar model grids, and fits this model to observed transmission spectra with MCMC. In exotune mode, the same machinery is applied to out-of-transit stellar spectra, allowing independent constraints on the contamination a host star can produce. The paper positions this as filling a gap: existing tools are either not public, not built for inference, serial-only, or

Load-bearing premise

The load-bearing premise is that a few uniform spots and faculae with fixed temperature contrasts, rendered with precomputed stellar model grids, faithfully represent the real heterogeneous surface of the star; if the grids are biased or the parameterization misses key effects, the retrieved contamination spectrum and any conclusion that no planetary atmosphere is required will be biased.

Editorial extensions

If this is right

  • Transmission spectra of small planets around M dwarfs can be checked against a TLSE-only model; if that model fits, an apparent atmospheric detection is not secure without further evidence.
  • Out-of-transit stellar spectra can be used to place data-driven priors on the amplitude and wavelength shape of stellar contamination for planets around the same host star.
  • Because stellar model grids are plug-in, results can be cross-checked across different synthetic spectral libraries without rewriting the retrieval code.
  • Fully parallelized MCMC makes the inference fast enough for routine application to JWST and HST transmission spectra.
  • The framework outputs model comparison metrics, sample spectra, and parameter samples, making contamination fits reproducible across different studies.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same spot-and-faculae retrieval applied to out-of-transit spectra could be time-resolved to track stellar activity evolution, an extension the paper does not explicitly develop.
  • A natural next step is adding wavelength-dependent limb darkening or temperature gradients across spots and faculae; the current simplified surface model may absorb such effects into broader parameter uncertainties.
  • If STCTM's TLSE-only fits are statistically preferred for a sizable sample of M-dwarf planets, estimates of how common featureless or atmosphere-free rocky planets are may need to be revisited; the paper enables such an analysis but does not quantify it.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. STCTM is an open-source Bayesian retrieval framework designed to model the transit light source effect (TLSE) in exoplanet transmission spectra and to infer the stellar surface parameters that are consistent with observations when no planetary signal is assumed. The paper describes the high-level workflow (configuration via input.toml, MCMC fitting with emcee, post-processing and publication-ready output), highlights the exotune sub-module for retrieving stellar parameters from out-of-transit spectra, and lists several recent peer-reviewed applications (Lim et al. 2023; Radica et al. 2025; Piaulet-Ghorayeb et al. 2024, 2025; Ahrer et al. 2025; Roy et al. 2023). The paper does not, however, provide an equation-level specification of the forward model, any validation benchmarks, or quantitative tests of the retrieval pipeline.

Significance. If the framework performs as claimed, it fills a clear community need: a flexible, modular, open-source tool for TLSE modeling and inference that supports arbitrary stellar model grids and is parallelized for tractable MCMC retrievals. The existing literature applications are a notable strength, providing indirect evidence of practical utility and community adoption. The emphasis on user-configurable inputs, reproducible outputs, and integration with widely-used stellar model grids (PHOENIX, SPHINX, MSG) is also valuable. The central scientific claim, however, rests on the correctness of the underlying forward model and the inference procedure. Because the paper itself contains no derivation or test of that model, the evidence currently available to a reader is mostly circumstantial; the paper's credibility would be substantially improved by explicitly presenting the forward-model equations and a validation demonstration.

major comments (3)
  1. [Summary / Main features of the code] The paper does not specify the forward model used to compute the TLSE contrast spectrum. It never defines the relationship between spot/faculae covering fractions, temperature contrasts, and the observed transmission spectrum. Without this specification, readers cannot assess whether the MCMC samples the intended parameter space or whether the implementation is correct. Please add the relevant equations (e.g., the standard contamination factor expression relating observed and true transit depth to spot/faculae properties), or at minimum give an explicit pointer to the equations used in Lim et al. 2023 / Radica et al. 2025 and the online documentation, and state how the stellar model grids are interpolated (e.g., via MSG) and how synthetic spectra for spot/faculae regions are combined.
  2. [Documentation / Uses of STCTM in the literature] No validation benchmarks or synthetic tests are presented in the paper. The 'Uses of STCTM in the literature' section provides examples of application, but a software paper whose central claim is a reliable inference framework needs more direct evidence. Please include at least one injection-recovery test showing that STCTM recovers input spot/faculae parameters and correctly identifies 'no planetary contribution' when none is present. Alternatively, include a representative retrieval on a published dataset and compare with the published posteriors. The test should be described in the paper itself (not only linked in the documentation) so that the reader can verify the pipeline's behavior.
  3. [Main features / exotune] The retrieved parameters are sensitive to systematic errors in the adopted stellar model grids and to the simplified parameterization of stellar heterogeneities as single-temperature spot/faculae regions with covering fractions. The paper mentions an error-inflation parameter only for exotune retrievals, not for the transmission-spectrum TLSE retrievals. Since the central claim is that STCTM can infer the range of stellar parameters 'compatible with the observations in the absence of any planetary contribution,' systematic model error could directly bias this range. Please add a discussion of this limitation and, ideally, demonstrate robustness by running a retrieval with two different model grids (e.g., PHOENIX vs. SPHINX) or by including an error-inflation / model-marginalization term in the transmission-spectrum case as well.
minor comments (4)
  1. [References] The Piaulet-Ghorayeb et al. 2025 reference has a URL pointing to a 2024 ApJ paper (2024ApJ...974L..10P). Please correct or clarify whether this is the intended citation.
  2. [Summary] The phrase 'spectral contrasts between bright and dark spots' is a little ambiguous: spots are usually the dark regions and faculae the bright ones. Consider rephrasing to 'spots and faculae' for clarity.
  3. [Similar Tools] For POSEIDON, the description 'including TLSE-only retrievals on transmission spectra' would benefit from a precise reference to the relevant POSEIDON documentation or paper, since the current sentence is vague.
  4. [Future Developments] The version number (v2.1.1) is mentioned, but no repository URL appears in the printed text. For reproducibility, please include the repository and documentation URLs explicitly in the manuscript body or a footnote.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: STCTM is a software/retrieval framework paper; its inference procedure is not presented as a prediction derived from its inputs.

full rationale

STCTM is described as a Bayesian retrieval framework for modeling the transit light source effect (TLSE). The paper makes no first-principles derivation of a physical result; its central claim is that the code can fit stellar surface parameters to transmission or out-of-transit spectra. That is an inversion/retrieval task, not a circular 'prediction' of the fitted quantity. The forward model is explicitly inherited from external, established sources (Rackham et al. 2018 for the TLSE formalism; PHOENIX/SPHINX model grids; MSG for interpolation), and the fitted parameters—spot/faculae covering fractions and temperatures—are standard retrieval parameters, not constants fitted and then renamed as predictions. The 'Uses of STCTM in the literature' section lists applications, several from the author's own prior work, but the paper does not lean on those citations as proof of the core inference claim; they are contextual citations of where the code was used. The absence of equation-level details or injection-recovery validation in this JOSS paper is a legitimate correctness/robustness concern, but it is not circularity: no step is shown to reduce to its own input by construction. Therefore the appropriate score is 0.

Assumptions & free parameters 4 free parameters · 3 assumptions · 0 invented entities

No new physical entities or effects are introduced. The framework's parameters are standard TLSE fitting parameters that users supply to the code; the paper itself makes no new physical derivation. The central claim depends on the validity of established TLSE theory and stellar model grids.

free parameters (4)
  • spot temperature contrast (T_spot/T_phot)
    Temperature difference between unocculted starspots and quiescent photosphere; retrieved by the code from transmission spectra or stellar spectra.
  • spot covering fraction (f_spot)
    Fractional area of the projected stellar disk covered by spots; primary fitted parameter in TLSE retrievals.
  • faculae covering fraction (f_fac)
    Coverage of bright facular regions if included; fitted to the observed spectra.
  • error inflation parameter (exotune)
    Optional scaling of the data uncertainties to absorb model mismatch in out-of-transit spectral retrievals.
assumptions (3)
  • domain assumption The TLSE can be modeled as the spectral contrast of unocculted spots and faculae against the stellar photosphere (Rackham et al. 2018).
    This is the physical basis of the framework, referenced in the Statement of Need and implemented as a forward model.
  • domain assumption PHOENIX and SPHINX stellar model grids produce reliable wavelength-dependent synthetic spectra for M dwarfs and their heterogeneities.
    The code supports any user-specified model grid; the fidelity of the retrieved contamination depends on these grids.
  • standard math Markov Chain Monte Carlo sampling (emcee) yields a statistically sound posterior when the likelihood is correctly specified.
    The framework performs Bayesian inference via MCMC; this is a standard statistical assumption.

how reviews work

0 comments
Cite this review

Pith. "Pith review of STCTM: a forward modeling and retrieval framework for stellar contamination and stellar spectra." pith.science (2026). https://pith.science/paper/NCLT3WRW

@misc{pith2026250819297,
  author       = {Pith},
  title        = {Pith review of: STCTM: a forward modeling and retrieval framework for stellar contamination and stellar spectra},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NCLT3WRW}},
  note         = {Machine review of arXiv:2508.19297}
}
read the original abstract

Transmission spectroscopy is a key avenue for the near-term study of small-planet atmospheres and the most promising method when it comes to searching for atmospheres on temperate rocky worlds, which are often too cold for planetary emission to be detectable. At the same time, the small planets that are most amenable for such atmospheric probes orbit small and cool M dwarf stars. As the field becomes increasingly ambitious in the search for signs of even thin atmospheres on small exoplanets, the transit light source effect (TLSE), caused by unocculted stellar surface heterogeneities, is becoming a limiting factor: it is imperative to develop robust inference methods to disentangle planetary and stellar contributions to the observed spectra. Here, I present STCTM, the STellar ConTamination Modeling framework, a flexible Bayesian retrieval framework to model the impact of the TLSE on any exoplanet transmission spectrum, and infer the range of stellar surface parameters that are compatible with the observations in the absence of any planetary contribution. With the "exotune" sub-module, users can also perform retrievals directly on out-of-transit stellar spectra in order to place data-driven priors on the extent to which the TLSE can impact any planet's transmission spectrum. The input data formats, stellar models, and fitted parameters are easily tunable using human-readable files and the code is fully parallelized to enable fast inferences. [shortened for arxiv; see full summary in the PDF]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

10 extracted references · 1 canonical work pages

  1. [1]

    S., Welbanks, L., Acuna, L., Allart, R., Coulombe, L.-P., Louca, A

    Ahrer, E.-M., Radica, M., Piaulet-Ghorayeb, C., Raul, E., Wiser, L. S., Welbanks, L., Acuna, L., Allart, R., Coulombe, L.-P., Louca, A. J., MacDonald, R. J., Saidel, M., Evans-Soma, T. M., Benneke, B., Christie, D., Beatty, T. G., Cadieux, C., Cloutier, R., Doyon, R., ... Schlichting, H. E. (2025). Escaping helium and a highly muted spectrum suggest a met...

    arXiv doi:10.48550/arxiv.2504.20428 2025

  2. [9]

    B., Doyon, R., Lafrenière, D., L’Heureux, A., & Lim, O

    https://doi.org/DOIunavailable 4 Radica, M., Piaulet-Ghorayeb, C., Taylor, J., Coulombe, L.-P., Benneke, B., Albert, L., Artigau, É., Cowan, N. B., Doyon, R., Lafrenière, D., L’Heureux, A., & Lim, O. (2025). Promise and peril: Stellar contamination and strict limits on the atmosphere composition of TRAPPIST- 1 c from JWST NIRISS transmission spectra.The A...

    doi:10.3847/2041-8213/ada381 2025

  3. [24]

    W., Lang, D., & Goodman, J

    https://doi.org/10.21105/joss.00024 Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. (2013). Emcee: The MCMC hammer. Publications of the Astronomical Society of the Pacific, 125(925),

    doi:10.21105/joss.00024 2013

  4. [41]

    steam world

    https://doi.org/ 10.3847/1538-4357/acabc2 Lim, O., Benneke, B., Doyon, R., MacDonald, R. J., Piaulet, C., Artigau, É., Coulombe, L.-P., Radica, M., L’Heureux, A., Albert, L., Rackham, B. V., de Wit, J., Salhi, S., Roy, P.-A., Flagg, L., Fournier-Tondreau, M., Taylor, J., Cook, N. J., Lafrenière, D., ... Darveau- Bernier, A. (2023). Atmospheric Reconnaissa...

    doi:10.3847/1538-4357/acabc2 2023

  5. [122]

    STCTM: a forward modeling and retrieval framework for stellar contamination and stellar spectra

    https://doi.org/10.3847/1538-4357/aaa08c Piaulet-Ghorayeb, (2025). STCTM: a forward modeling and retrieval framework for stellar contamination and stellar spectra. Journal of Open Source Software, 0(0),

    doi:10.3847/1538-4357/aaa08c 2025

  6. [123]

    P., Tollerud, E

    https://doi.org/10.3847/1538-3881/aabc4f Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., Greenfield, P., Droettboom, M., Bray, E., Aldcroft, T., Davis, M., Ginsburg, A., Price-Whelan, A. M., Kerzendorf, W. E., Conley, A., Crighton, N., Barbary, K., Muna, D., Ferguson, H., Grollier, F., Parikh, M. M., Nair, P. H., ... Streicher, O. (2013). Astro...

    doi:10.3847/1538-3881/aabc4f 2013

  7. [306]

    R., Millman, K

    https: //doi.org/10.1086/670067 Harris, C. R., Millman, K. J., Walt, S. J. van der, Gommers, R., Virtanen, P., Cournapeau, D., Wieser, E., Taylor, J., Berg, S., Smith, N. J., Kern, R., Picus, M., Hoyer, S., Kerkwijk, M. H. van, Brett, M., Haldane, A., Río, J. F. del, Wiebe, M., Peterson, P., ... Oliphant, T. E. (2020). Array programming with NumPy. Nature...

    doi:10.1086/670067 2020

  8. [4503]

    https://doi.org/10.21105/joss.04503 Foreman-Mackey, D. (2016). Corner.py: Scatterplot matrices in python.The Journal of Open Source Software, 1(2),

    doi:10.21105/joss.04503 2016

Show all 10 references
  1. [4602]

    de, Doyon, R., Rackham, B

    https://doi.org/10.21105/joss.04602 TRAPPIST-1 JWST Community Initiative, Wit, J. de, Doyon, R., Rackham, B. V., Lim, O., Ducrot, E., Kreidberg, L., Benneke, B., Ribas, I., Berardo, D., Niraula, P., Iyer, A., Shapiro, A., Kostogryz, N., Witzke, V., Gillon, M., Agol, E., Meadow...

    2024 doi

  2. [8269]

    D., Caloca, G., Mansfield, M., Zieba, S., Piaulet, C., Benneke, B., Filippazzo, J., May, E., Roy, P.-A., Kreidberg, L., & Stevenson, K

    https://doi.org/DOIunavailable 3 Bell, T., Ahrer, E.-M., Brande, J., Carter, A., Feinstein, A. D., Caloca, G., Mansfield, M., Zieba, S., Piaulet, C., Benneke, B., Filippazzo, J., May, E., Roy, P.-A., Kreidberg, L., & Stevenson, K. (2022). Eureka!: An End-to-End Pipeline for JW...

    2022

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.