Precise photoionisation treatment and hydrodynamic effects in atmospheric modelling of warm and hot Neptunes

Observational breakthroughs in the exoplanet field of the last decade motivated the development of numerous theoretical models describing atmospheres and mass loss, which is believed to be one of the main drivers of planetary evolution. We aim to outline for which types of close-in planets in the Neptune-mass range the accurate treatment of photoionisation effects is most relevant concerning atmospheric escape and the parameters relevant for interpreting observations. We developed the CHAIN (Cloudy e Hydro Ancora INsieme) model combining 1D hydrodynamic upper atmosphere model with the non-LTE photoionisation and radiative transfer code Cloudy accounting for photochemistry, detailed atomic level populations, and chemical reactions for all elements up to zinc. We apply CHAIN to model the upper atmospheres of a range of Neptune-like planets with masses between 1 and 50 M$_{\oplus}$, varying also the orbital parameters. For the majority of warm and hot Neptunes, we find slower and denser outflows, with lower ion fractions, compared to the predictions of the hydrodynamic model alone. Furthermore, we find significantly different temperature profiles between CHAIN and the hydrodynamic model alone, though the peak values are similar for similar atmospheric compositions. The mass-loss rates predicted by CHAIN are higher for hot, strongly irradiated planets and lower for more moderate planets. All differences between the two models are strongly correlated with the amount of high-energy irradiation. Finally, we find that the hydrodynamic effects impact significantly ionisation and heating. The impact of the precise photoionisation treatment provided by Cloudy strongly depends on the system parameters. This suggests that some of the simplifications typically employed in hydrodynamic modelling might lead to systematic errors when studying planetary atmospheres, even at a population-wide level.