“…We model the lower and middle atmosphere with the photochemical-thermal structure model of Lavvas & Arfaux (2021), assuming global redistribution of heat and constituents. At pressures lower than 1 μbar, we use the multispecies hydrodynamic escape model of Koskinen et al (2022), updated to include the elements H, He, Mg, Si, Fe, O, C, N, Na, K, S, Ca, and relevant ions (Huang et al 2023). For the stellar spectrum, ranging from the FUV to the infrared, we use the spectral energy distribution (SED) computed with the PHOE-NIX stellar atmosphere code (Husser et al 2013) based on the stellar parameters from Lendl et al (2020).…”