A JWST transiting survey of FGK stellar limb darkening: empirical evidence for quadratic laws and atmospheric model comparisons

This study presents an analysis of stellar limb-darkening (LD) using high-precision transit observations obtained with the James Webb Space Telescope (JWST). The research focuses on seven exoplanets orbiting FGK host stars, spanning temperatures from 4200 K to 6800 K. The wide wavelength coverage and high signal-to-noise ratio (S/N) provided by NIRISS/SOSS and NIRSpec/PRISM enable precise constraints on the wavelength-dependent limb darkening.

Methodology and Target Selection

The study utilized JWST red-optical to near-IR transit data for these seven exoplanets. The targets were selected based on favorable orbital parameters, ensuring that a large portion of the stellar disk is sampled during transit. The maximum intensity (mu max) probed was used to assess suitability; targets were chosen where mu max was high, with all analyzed targets having mu max > 0.95.

The data reduction and light-curve fitting process involved using the Fast InfraRed Exoplanet Fitting Lyghtcurve (FIREFLy) suite. The white light curves were fitted using four different limb darkening parameterizations: a linear, a quadratic, three-parameter, and four-parameter non-linear law. For model selection, the Bayesian Predictive Information Criterion Simplified (BPICS) statistic was utilized to determine the best predictive power for each fit.

Empirical Findings on Limb Darkening Laws

When comparing the results across individual targets:

• For each individual target the statistics between the higher-order limb darkening laws are not often definitive, with only the linear law clearly disfavored in all cases.

• In most cases, we find the quadratic law usually has a better fit and marginally lower BPICS value... when compared to the 3-parameter and 4-parameter respectively.

When aggregating the statistics across all seven stars:

• Six of the seven stars favor a quadratic LD law.

• The aggregate analysis shows that the quadratic and 3-parameter LD laws offer good fits to the JWST transit light curves while the 4-parameter law is disfavored.

Key Results and Model Bias The study found significant empirical evidence regarding the applicability of standard models:

1. Quadratic Law Preference: The quadratic limb-darkening law is statistically preferred for most FGK stars.

2. Minimal Bias: The quadratic model "recovers consistent intensity profiles for >95% of the stellar disk area, and introduces a minimal bias in the derived transit depths of only about 14 ppm (1 sigma) for 6/7 targets. This finding contrasts with previous studies relying on stellar models."

3. Linear Nature: Empirical observations show that FGK limb darkening is generally more linear than models predict (delta 0.1).

4. Model Discrepancies: The study identified wavelength-independent offsets between data and model quadratic coefficients, which were minimized by adopting mu min = 0.2 in intensity calculations. This was attributed to models overpredict[ing] limb-darkening near the limb where the plane-parallel approximation breaks down.

Addressing Model Uncertainties

To reconcile the discrepancies between empirical data and theoretical models, specific corrections were introduced:

• For spherical PHOENIX models, a mu rescaling method was derived using a tau = 1 photospheric radius. This correction is necessary because "Exotic-LD (Grant & Wakeford 2024) does not apply it natively and will otherwise over-predict limb darkening."

• The study found that the models tend to overpredict the overall limb darkening, as quantified by the linear coefficient, which was observed to be larger for 1-D models compared to 3-D models.

Recommendations for Future Analysis

Based on these findings, the authors provide specific recommendations for researchers using JWST data:

• We generally recommend fitting the LD with a quadratic law, statistically checking the quality of fit and transit depth vs. a 3-parameter law.

• The study recommends first fitting for gray LD offsets [the model differences]... before re-fitting setting priors on the coefficients with offsets applied.

• We provide recommended limb-darkening offset priors for use in JWST transit analyses, enabling more accurate constraints on exoplanet transmission spectra while accounting for residual stellar model uncertainties.

Conclusion

The study concludes that the quadratic law is statistically preferred, offers a minimal bias, and provides a robust empirical benchmark. The authors emphasize that the key result is not that higher-order laws are incorrect, but that they are unnecessary, as the additional parameters in these laws were poorly constrained by the data.

Based on my role as an expert AI researcher, I have analyzed this paper to identify critical methodological improvements that can be implemented in next-generation astronomical and exoplanet analysis AI systems.

The core of the improvement lies in moving beyond simply fitting a model (e.g., quadratic law) to correctly interpreting the resulting parameters by accounting for systematic, model-dependent biases identified in the paper.

Here are the specific improvements and what they will enable an improved AI system can do:

The Improvement: Integrate a derived method for calculating the effective limb angle (mu eff) into spherical atmosphere models (specifically PHOENIX). Instead of using the standard mu=0 definition, the AI system must calculate mu eff —the angle corresponding to where the slant optical depth (tau slant) reaches unity.

What it enables: The AI system can accurately rescale and interpolate intensity profiles from spherical models to match a consistent plane-parallel grid, eliminating the systematic overestimation of limb darkening inherent in 1D spherical geometry. This allows for a direct, physics-based comparison between theoretical predictions and observed JWST data without relying on ad-hoc numerical fixes.

The Improvement: Incorporate the calculated gray offset matrices (, delta) from Table 3 into the fitting pipeline. These offsets are not fixed; they are functions of the stellar parameters (T eff, g) and must be applied based on a specified mu min (e.g., 0.2).

What it enables: The AI system can perform an offset-corrected Bayesian fit. Instead of simply flagging model discrepancies, the the system can automatically adjust its expected prior parameters (and delta) based on the star's characteristics, ensuring that a high-quality fit is achieved even when using suboptimal stellar models (e.g., achieving chi 2 v about 1). This moves AI analysis from finding the best fit to finding the most physically accurate parameter set.

The Improvement: Implement a dynamic selection mechanism that prioritizes using mu min = 0.2 as the default cutoff for intensity calculations, rather than arbitrarily excluding small mu values. The AI system must also be able to calculate and apply the associated offset priors derived from this specific cutoff.

What it enables: The AI system can automatically recognize and compensate for the fact that models tend to overpredict strength near the extreme limb (mu < 0.2). By prioritizing this physically motivated cut, the the AI ensures that its parameter estimates are based on the region of maximum flux contribution, thereby minimizing bias in derived transit depths.

The Improvement: Integrate the Bayesian Predictive Information Criterion Simplified (BPICS) statistic as a primary feature in model selection algorithms. The system must be trained to interpret the resulting BPICS values and associated odds ratios (e.g., 20:1 or 98.8%) rather than relying solely on simple likelihood comparisons.

What it enables: The AI can perform rigorous, statistically informed model comparison across different LD laws (Linear vs. Quadratic vs. 3-Parameter vs. 4-Parameter) for every single target star, automatically selecting the most predictive model and providing a probabilistic justification for its choice, significantly reducing the chance of misclassification based on simple fit metrics.

The improved AI system will be capable of performing Bias-Corrected Chromatic Transit Analysis. It will not just report that a quadratic law is preferred; it will:

1. Physically correct the theoretical model predictions (PHOENIX, etc.) using mu eff and mu min.

2. Automatically apply the necessary parameter offsets (, delta) based on the star's T eff.

3. Select a model based on its predictive power (BPICS) while maintaining high fidelity to the empirical JWST data.4

We present a study of stellar limb-darkening using JWST transit observations of seven exoplanets orbiting FGK host stars, spanning 4200-6800 K. The wide wavelength coverage and high S/N of NIRISS/SOSS and NIRSpec/PRISM enable precise constraints on the wavelength-dependent limb-darkening. Using Bayesian model selection, we find that the quadratic limb-darkening law is statistically preferred over higher-order laws for most FGK stars, recovers consistent intensity profiles for > 95% of the stellar disk area, and introduces a minimal bias in the derived transit depths of only about 14 ppm (1 σ) for 6/7 targets. This contrasts with previous studies relying on stellar models. We compare the empirically derived quadratic coefficients to predictions from the PHOENIX, MPS-ATLAS, MURaM, and Stagger stellar atmosphere grids. We introduce a quadratic limb-darkening parameterization in terms of limb intensity and curvature at mid- μ (δ), finding that empirical FGK limb darkening is generally more linear than models predict (δ 0.1). We identify wavelength-independent offsets between data and model quadratic coefficients, minimized by adopting μ min = 0.2 in intensity calculations; we attribute this in part to models overpredicting limb-darkening near the limb where the plane-parallel approximation breaks down. For spherical PHOENIX models, we derive a μ rescaling method using a τ= 1 photospheric radius. With these corrections, residual offsets are minimized and all stellar models achieve statistically acceptable fits. We provide recommended limb-darkening offset priors for use in JWST transit analyses, enabling more accurate constraints on exoplanet transmission spectra while accounting for residual stellar model uncertainties.