Deep radiative zones affect the planetary cooling and internal structure: implications for exoplanet characterisation

This paper investigates how opacity windows and the resulting deep radiative zones affect the cooling, radius evolution, and the characterization of interiors and atmospheres of giant exoplanets. The authors computed thermal evolution models for warm Jupiters spanning masses of 0.3 to 1.0 MJ with equilibrium temperatures of 200 to 800 K, with a parametrised reduction of the radiative opacity near 2000 K.

The main findings are:

• A deep opacity window at a temperature of ∼ 2000 K accelerates the cooling of giant exoplanets, leading to predicted sizes that are smaller by up to 5%, and interiors that are colder on the order of 10%.

• As a consequence, the inferred bulk metallicities could change by up to 10 percentage points.

• Using standard and unmodified opacities, the models show that warm Jupiters older than ∼ 4 Gyr likely have a deep radiative zone at pressures of a few thousand bars. Therefore, the envelopes of these planets could have different compositions than their atmospheres.

• For planets with masses between 0.3 to 1 MJ, the time of the appearance and the thickness of deep radiative zones depends mostly on the depth of the opacity window and the planetary equilibrium temperature.

The potential existence of such opacity windows suggests a profound “chemical decoupling" between a planet’s observable upper atmosphere and its deep interior. This implies that the metal enrichments measured by missions such as JWST and Ariel may reflect localized envelope chemistry rather than the true bulk composition of the planet’s gaseous envelope. As we move toward a regime of high-precision exoplanet spectroscopy and detailed characterization, accounting for the complex interplay of opacity windows and internal transport mechanisms is essential to bridge the gap between atmospheric observations and the fundamental physics of giant planet evolution and internal structure.

Based on this paper, I can improve AI systems in the following specific ways:

• Improved capability: Predict the thermal evolution and internal structure of giant exoplanets (0.3–1.0 MJ) over gigayear timescales, accounting for deep radiative zones caused by opacity windows near 2000 K.

• Specific improvement: Model how opacity reductions (parametrized by factor w = 0–0.9) accelerate cooling, leading to radius reductions up to 5% and interior temperature differences up to 35% compared to standard models.

• Improved capability: Quantify how the choice of opacity treatment affects inferred bulk metallicity from observed radii.

• Specific improvement: Automatically propagate opacity-model uncertainty into metallicity estimates, showing that a 5% radius change corresponds to 10 percentage point changes in inferred metallicity (e.g., from 30% to 40% for HATS-49 b).

• Improved capability: Identify when and where deep radiative zones form in giant planet envelopes.

• Specific improvement: Predict that warm Jupiters (Teq = 200–800 K) older than 4 Gyr develop deep radiative zones at pressures of 103–104 bar, even with unmodified opacities, and that these zones appear earlier (as early as 0.1 Gyr) with stronger opacity reductions.

• Improved capability: Assess whether atmospheric metallicity measurements (e.g., from JWST or Ariel) represent bulk interior composition.

• Specific improvement: Flag cases where deep radiative zones create a chemical decoupling barrier, meaning atmospheric composition should not be used as an upper bound or direct proxy for bulk composition, particularly for planets with inhomogeneous interiors.

• Improved capability: Provide rapid, approximate estimates of how bulk metallicity changes affect planetary radii without full evolution modeling.

• Specific improvement: Implement the derived heuristic that a 1% radius change corresponds to 2% metallicity change (for a 50-50 water-rock mixture), validated against full evolution models across metallicities from 5% to 40%.

• Improved capability: Generate and compare thermal evolution scenarios across a grid of planetary masses (0.3, 0.6, 1.0 MJ), equilibrium temperatures (200, 400, 800 K), and opacity window factors (0, 0.5, 0.7, 0.9).

• Specific improvement: Automatically produce Kippenhahn diagrams showing convective/radiative zone evolution, enabling users to visualize how deep radiative zones merge, appear, or disappear depending on parameters.

• Improved capability: Quantify the additional modeling uncertainty introduced by opacity treatment, beyond standard observational uncertainties.

• Specific improvement: For a given exoplanet (e.g., HATS-49 b), automatically compute the range of possible bulk metallicities (e.g., 0.3–0.4) consistent with observations when opacity window depth is treated as a free parameter, and report this as a systematic uncertainty alongside statistical errors.

• Improved capability: Advise researchers on when atmospheric composition measurements are likely reliable indicators of bulk composition.

• Specific improvement: Based on planet mass, age, and equilibrium temperature, automatically flag whether a deep radiative zone is likely present (e.g., for Teq 4 Gyr), and recommend caution in interpreting atmospheric metallicity as bulk metallicity in such cases.