Where does the simplified Stellar Contamination Model fail in Exoplanet Transmission Spectroscopy?

TL;DR

This study evaluates the limitations of the simplified Stellar Contamination Model in exoplanet transmission spectroscopy, highlighting the importance of including limb darkening and stellar heterogeneity for accurate atmospheric characterization.

Contribution

The paper introduces ECLIPSE-Xlambda, a pixel-resolved framework that improves stellar contamination modeling by incorporating limb darkening and active-region distribution, surpassing the simplified R-TLSE approach.

Findings

Disc-averaged corrections can differ by up to 400 ppm in the optical.

Limb darkening effects are minimal in the near-infrared.

Extreme faculae are required to explain observed spectral slopes without atmospheric contributions.

Abstract

Stellar photospheric heterogeneities (e.g., starspots, faculae) distort the stellar spectrum in transit and imprint wavelength-dependent biases on the planet-to-star radius ratio (Transit Light Source Effect, TLSE). The Rackham-TLSE (R-TLSE) prescription applies a disc-averaged correction based solely on filling factor and spectral contrast, but transmission spectroscopy also depends on limb darkening, active-region distribution, and transit geometry. We include these in a pixel-resolved framework, ECLIPSE-Xlambda, and run idealised noise-free model-model comparisons to R-TLSE. For LHS 1140 b, K2-18 b, and WASP-69 b, disc-averaged corrections differ from the pixel model by up to about 400 ppm in the optical for active hosts and non-equatorial transits, but stay below about 10 ppm in the near-infrared where limb darkening is weak. We then apply both approaches to the JWST/NIRISS SOSS spectrum of LHS 1140 b. With limb darkening set to zero, ECLIPSE-Xlambda recovers stellar-contamination parameters matching the reference R-TLSE solution, confirming consistency in the disc-averaged limit. With wavelength-dependent limb darkening, reproducing the short-wavelength slope via stellar contamination alone requires hot faculae (delta Tfac about 600 K; ffac about 0.35), equivalent to a circular facular region of radius about 0.6 Rstar (about 60% of the stellar radius) on the disc; such an extended unocculted region is physically unlikely even for an active M dwarf. Purely stellar contamination would therefore require extreme faculae, whereas a genuine atmospheric contribution complementing a more modest facular signal is more plausible. These results delineate the validity regime of R-TLSE and underscore the need for geometry-aware stellar-heterogeneity models including limb darkening in high-precision transmission spectroscopy.