Constraining planetary mass-loss rates by simulating Parker wind profiles with Cloudy

TL;DR

This paper introduces a framework combining Parker wind models with Cloudy simulations to better constrain exoplanet atmospheric mass-loss rates and temperatures, reducing degeneracies in spectroscopic analysis.

Contribution

It presents an iterative method integrating Cloudy with Parker wind models to constrain thermospheric temperature and mass-loss rates from observational data.

Findings

Ruled out large portions of parameter space for HD 209458 b.

Estimated a thermospheric temperature of 8200 K for HD 209458 b.

Constrained mass-loss rates for multiple exoplanets.

Abstract

Models of exoplanet atmospheres based on Parker-wind density and velocity profiles are a common choice in fitting spectroscopic observations tracing planetary atmospheric escape. Inferring atmospheric properties using these models often results in a degeneracy between the temperature and the mass-loss rate, and thus provides weak constraints on either parameter. We present a framework that can partially resolve this degeneracy by placing more stringent constraints on the expected thermospheric temperature. We use the photoionization code Cloudy within an iterative scheme to compute the temperature structure of a grid of 1D Parker wind models, including the effects of radiative heating/cooling, as well as the hydrodynamic effects (expansion cooling and heat advection). We constrain the parameter space by identifying models that are not self-consistent through a comparison of the simulated temperature in the He 10830 $\r{A}$ line-forming region to the temperature assumed in creating the models. We demonstrate this procedure on models based on HD 209458 b. By investigating the Parker wind models with an assumed temperature between 4000 and 12000 K, and a mass-loss rate between $10^{8}$ and $10^{11}$ g s$^{-1}$, we are able to rule out a large portion of this parameter space. Furthermore, we fit the models to previous observational data and combine both constraints to find a preferred thermospheric temperature of $T=8200^{+1200}_{-1100}$ K and a mass-loss rate of $\dot{M}=10^{9.84^{+0.24}_{-0.27}}$ g s$^{-1}$ assuming a fixed atmospheric composition and no gas pressure confinement by the stellar wind. Using the same procedure, we constrain the temperatures and mass-loss rates of WASP-69 b, WASP-52 b, HAT-P-11 b, HAT-P-18 b and WASP-107 b.