A Closed-Form Analytical Theory of Non-Isobaric Transmission Spectroscopy for Exoplanet Atmospheres

TL;DR

This work addresses the limitation of isobaric opacity in exoplanet transmission spectroscopy by introducing a general power-law opacity $\kappa \propto P^{n}$ with an explicit reference opacity $\kappa_{0}$ at $P_{0}$. Using an Abel transform of the radial absorption coefficient in an isothermal, hydrostatic atmosphere, it derives a closed-form expression for the effective transit radius $R_p = R_0 + \frac{H}{1+n} (\gamma + \ln \tau_0)$, where $\tau_0$ is linked to $\kappa_{0}$, $P_{0}$, gravity $g$, and scale height $H$. The method also prescribes constructing the wavelength-dependent opacity exponent $n(\lambda)$ from opacity grids to form an effective $n_{\rm mix}(\lambda)$ for the mixed atmosphere. Application to Earth and WASP-39b shows significant improvements over the isobaric model, with large reductions in residuals and information criteria, demonstrating that pressure-dependent opacity captures essential vertical structure in transmission spectra and enabling efficient semi-analytical retrievals for upcoming JWST and ARIEL data.

Abstract

Analytical models are essential for building physical intuition and guiding the interpretation of exoplanet observations by clarifying the dependencies that shape atmospheric signatures. We present a generalization of the classical isothermal, isobaric transmission model by allowing the opacity to vary with pressure as a power law, $κ\propto P^{n}$, and explicitly defining the reference opacity $κ_{0}$ at a chosen pressure $P_{0}$. By treating the slant optical depth as an Abel transform of the radial absorption coefficient, we derive a closed-form expression for the effective transit radius in a hydrostatic, isothermal atmosphere with pressure-dependent opacity. The solution provides a compact framework for exploring non-isobaric effects and explicitly links the vertical opacity gradient to observable spectral features. We benchmark the model against empirical transmission spectra of Earth and the hot Jupiter WASP-39b, finding a significantly improved fit relative to the isobaric formula. This generalized expression offers a physically interpretable foundation for analyzing high-precision spectra from JWST and upcoming ARIEL observations, and can serve as a basis for semi-analytical retrieval approaches optimized for computational efficiency.

Figures (2)

• Figure 1: Comparison between the analytical transmission models and the empirical spectra. Top: Earth's transmission spectrum from Macdonald & Cowan (2019, 2023) derived from ACE-FTS solar-occultation measurements. The empirical spectrum is shown as a solid black line, the isobaric analytical model as a blue curve, and the generalized non-isobaric model as an orange curve. Bottom: JWST transmission spectrum of WASP-39b from the ERS program (Rustamkulov et al. 2023; Tsai et al. 2023), shown as black points with $1\sigma$ error bars, compared with the same analytical models (blue: isobaric; orange: non-isobaric). In both cases the non-isobaric formulation reproduces the relative depths and shapes of the major molecular bands, notably the SO$_2$ and CO$_2$ features, more faithfully than the classical isobaric approximation.
• Figure 2: Pressure- and wavelength-dependent opacity behavior for the main atmospheric absorbers considered in the WASP--39b analysis. Top: Monochromatic opacities $\kappa_i(\lambda,P_0)$ at $P_0=1$ bar and $T=1100$ K for the dominant species (H$_2$O, CO$_2$, CO, SO$_2$, and CH$_4$), computed at $R=1000$ using the ExoMol/HITEMP/HITRAN line lists. Bottom: Corresponding pressure-dependence exponents $n_i(\lambda)=\partial\ln\kappa_i/\partial\ln P|_{T,\lambda}$ and the abundance-weighted mixture value $n_{\mathrm{mix}}(\lambda)$ for WASP--39b, smoothed over 0.3 $\mu$m. The results illustrate that $n(\lambda)$ varies strongly with wavelength and species, reflecting the dominance of different opacity sources, while the weighted $n_{\mathrm{mix}}(\lambda)$ provides an effective parameter entering the generalized analytical transmission formula (Eq. 21).