Richard Chatterjee

Abstract:

Following seminal studies such as Muñoz’s 2007 work on HD 209458b, which simulated heavy element escape beyond the Roche lobe, one-dimensional hydrocodes have flourished, routinely solving the Euler equations to model transonic outflows across an increasingly diverse population of exoplanets. However, the modelling frontier of escape is often shaped by the hand-off from continuum to rarefied flow (Kn ≳ 0.1) and non-equilibrium processes. Molecular-kinetic techniques, long the workhorse of Solar-System aeronomy, naturally bridge this gap, providing a self-consistent description of collisional, transitional and free-molecular regimes in a single framework. Here we make the case for a concerted push toward large-scale molecular-kinetic simulations of exoplanet outflows, highlighting two end-member scenarios along the escape spectrum where forthcoming observations may allow the theory to be tested and refined.Cosmic Shoreline. Characterising the transition from Jeans (particle-by-particle) escape to subsonic and ultimately transonic bulk outflow remains an open problem in escape theory. The onset of rapid escape (~1 bar Myr⁻¹) as ionising irradiation increases is a key parameter defining the phase boundary between airless and airy rocky worlds—the “Cosmic Shoreline” (Zahnle & Catling 2017; Ji et al. 2025). Johnson et al. (2013) combined an analytic treatment with Direct Simulation Monte Carlo (DSMC; Bird 1994) to derive a critical heating rate for triggering transonic flow, working with the ansatz that the scaling of this transition extends smoothly from Pluto- to Earth-sized bodies. We will present new DSMC simulations that probe this transition for high-molecular-weight atmospheres on Earth-mass and super-Earth planets, refining the dynamics of rapid escape across this regime.Helium triplet and fractionation. Fractionation may help explain some of the non-detections of the neutral-helium triplet (1083 nm) in giant-planet outflows (Schulik & Owen 2024). Multi-fluid hydrodynamics simulations have found that the neutral helium can actually be accelerated by gravity to accrete out of the flow at a downward velocity of ~1 km s⁻¹ (Xing et al. 2023; Schulik & Owen 2024). We note that the ratio of the slip velocity to the thermal speed of the outflow scales with the Knudsen number for collisionality, ΔU/ Vth~ KnHe . Thus, we will discuss how a significant slip velocity may require Kn ≳ 0.1, a regime in which the fractionation process may be better described with molecular-kinetics, possibly with implications for predictions of the transit depth of the helium triplet.Moreover, the Direct Simulation Monte Carlo (DSMC) method offers some desirable properties over hydrocodes: it scales naturally to fully three-dimensional geometries, albeit at significant computational cost, and naturally treats non-equilibrium phenomena such as photoelectron heating and excited-state populations.

Abstract:

The “cosmic shoreline”, a semi-empirical relation that separates airless worlds from worlds with atmospheres as proposed by Zahnle & Catling (2017), is now guiding large-scale JWST surveys aimed at detecting rocky exoplanet atmospheres. We expand upon this framework by revisiting the shorelines using existing hydrodynamic escape models applied to Earth-like, Venus-like, and steam atmospheres for rocky exoplanets, and we estimate energy-limited escape rates for CH4 atmospheres. We determine the critical instellation required for atmospheric retention by calculating time-integrated atmospheric mass loss. Our analysis introduces a new metric for target selection in the Rocky Worlds DDT and refines expectations for rocky planet atmosphere searches in Cycle 4. Exploring initial volatile inventory ranging from 0.01% to 1% of planetary mass, we find that its variation prevents the definition of a unique clear-cut shoreline, though non-linear escape physics can reduce this sensitivity to initial conditions. Additionally, uncertain distributions of high-energy stellar evolution and planet age further blur the critical instellations for atmospheric retention, yielding broad shorelines. Hydrodynamic escape models find atmospheric retention is markedly more favorable for higher-mass planets orbiting higher-mass stars, with carbon-rich atmospheres remaining plausible for 55 Cancri e despite its extreme instellation. Dedicated modelling efforts are needed to better constrain the escape dynamics of secondary atmospheres, such as the role of atomic line cooling, especially for Earth-sized planets. Finally, we illustrate how density measurements can be used to statistically test the existence of the cosmic shorelines, emphasizing the need for more precise mass and radius measurements.

Abstract:

We present a JWST Near Infrared Imager and Slitless Spectrograph/Single Object Slitless Spectroscopy transmission spectrum of the super-Earth GJ 357 b: the first atmospheric observation of this exoplanet. Despite missing the first 40 per cent of the transit due to using an out-of-date ephemeris, we still recover a transmission spectrum that does not display any clear signs of atmospheric features. We perform a search for Gaussian-shaped absorption features within the data but find that this analysis yields comparable fits to the observations as a flat line. We compare the transmission spectrum to a grid of atmosphere models and reject, to 3 confidence, atmospheres with metallicities solar (4 g mol−1) with clouds at pressures down to 0.01 bar. We analyse how the retention of a secondary atmosphere on GJ 357 b may be possible due to its higher escape velocity compared to an Earth-sized planet and the exceptional inactivity of its host star relative to other M2.5V stars. The star’s XUV luminosity decays below the threshold for rapid atmospheric escape early enough that the volcanic revival of an atmosphere of several bars of CO is plausible, though subject to considerable uncertainty. Finally, we model the feasibility of detecting an atmosphere on GJ 357 b with MIRI/LRS, MIRI photometry, and NIRSpec/G395H. We find that, with two eclipses, it would be possible to detect features indicative of an atmosphere or surface. Further to this, with three to four transits, it would be possible to detect a 1 bar nitrogen-rich atmosphere with 1000 ppm of CO.