LAUSR

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.1MJ - 5MJ and orbital radii of 3, 5, 10 au. Low dust opacity in the gap allows radiation to penetrate deeper and warm the midplane 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 midplane 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, CO). We find that icelines experience dramatic shifts relative to gapless models: up to 6.5 au (or 71%) closer to the star and 4.3 au (or 100%) closer to the midplane. 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.