ABSTRACT: The composition of forming planets is strongly affected by the protoplanetary disc’s thermal structure. This thermal structure is predominantly set by dust radiative transfer and viscous (accretional) heating and can be impacted by gaps – regions of low dust and gas density that can occur when planets form. The effect of variations in dust surface density on disc temperature has been poorly understood to date. In this work, we use the radiative transfer code MCMax to model the 2D dust thermal structure with individual gaps corresponding to planets with masses of 0.1 MJ –5 MJ and orbital radii of 3, 5, and 10 au. Low dust opacity in the gap allows radiation to penetrate deeper and warm the mid-plane by up to 16 K, but only for gaps located in the region of the disc where stellar irradiation is the dominant source of heating. In viscously heated regions, the mid-plane of the gap is relatively cooler by up to 100 K. Outside of the gap, broad radial oscillations in heating and cooling are present due to disc flaring. These thermal features affect local dust–gas segregation of volatile elements (H2O, CH4, CO2, and CO). We find that icelines experience dramatic shifts relative to gapless models: up to 6.5 au (or 71 per cent) closer to the star and 4.3 au (or 100 per cent) closer to the mid-plane. While quantitative predictions of iceline deviations will require more sophisticated models, which include transport, sublimation/condensation kinetics, and gas–dust thermal decoupling in the disc atmosphere, our results suggest that planet-induced iceline variations represent a potential feedback from the planet on to the composition of material it is accreting.

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability1. However, no unambiguous photochemical products have been detected in exoplanet atmospheres so far. Recent observations from the JWST Transiting Exoplanet Community Early Release Science Program2,3 found a spectral absorption feature at 4.05 μm arising from sulfur dioxide (SO2) in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 MJ) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of around 1,100 K (ref. 4). The most plausible way of generating SO2 in such an atmosphere is through photochemical processes5,6. Here we show that the SO2 distribution computed by a suite of photochemical models robustly explains the 4.05-μm spectral feature identified by JWST transmission observations7 with NIRSpec PRISM (2.7σ)8 and G395H (4.5σ)9. SO2 is produced by successive oxidation of sulfur radicals freed when hydrogen sulfide (H2S) is destroyed. The sensitivity of the SO2 feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of about 10× solar. We further point out that SO2 also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations.