p-winds: an open-source Python code to model planetary outflows and upper atmospheres
arXiv:2111.11370 · doi:10.1051/0004-6361/202142038
Abstract
Atmospheric escape is considered to be one of the main channels for evolution in sub-Jovian planets, particularly in their early lives. While there are several hypotheses proposed to explain escape in exoplanets, testing them with atmospheric observations remains a challenge. In this context, high-resolution transmission spectroscopy of transiting exoplanets for the metastable helium triplet (He 2S) at nm has emerged as a reliable technique to observe and measure escape. To aid in the prediction and interpretation of metastable He transmission spectroscopy observations, we developed the code p-winds. This is an open-source, fully documented, scalable Python implementation of the one-dimensional, purely H+He Parker wind model for upper atmospheres coupled with ionization balance, ray-tracing, and radiative transfer routines. We demonstrate an atmospheric retrieval by fitting p-winds models to the observed metastable He transmission spectrum of the warm Neptune HAT-P-11 b, and take into account the variation of the in-transit absorption caused by transit geometry. For this planet, our best fit yields a total atmospheric escape rate of approximately g s and wind temperature of K. The range of retrieved mass loss rates increases significantly when we let the H atom fraction be a free parameter, but the posterior distribution of the latter remains unconstrained by He observations alone. The stellar host limb darkening does not have a significant impact in the retrieved escape rate or outflow temperature for HAT-P-11 b. Based on the non-detection of escaping He for GJ 436 b, we are able to rule out total escape rates higher than g s at 99.7% (3) confidence.
13 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics. We encourage the community to send comments and test the code, which is openly available at https://github.com/ladsantos/p-winds
References in corpus (24)
- Array Programming with NumPy
- A giant comet-like cloud of hydrogen escaping the warm Neptune-mass exoplanet GJ 436b
- Spectrally resolved detection of sodium in the atmosphere of HD189733b with the HARPS spectrograph
- Using Stellar Limb-Darkening to Refine the Properties of HD 209458b
- Lyman- Transit Spectroscopy and the Neutral Hydrogen Tail of the Hot Neptune GJ436b
- Hubble PanCET: An extended upper atmosphere of neutral hydrogen around the warm Neptune GJ 3470 b
- High-resolution confirmation of an extended helium atmosphere around WASP-107b
- The long egress of GJ~436b's giant exosphere
- Reconnaissance of the TRAPPIST-1 exoplanet system in the Lyman- line
- Mass loss rate and local thermodynamic state of KELT-9 b thermosphere from the hydrogen Balmer series
- Confirmation of WASP-107b's extended Helium atmosphere with Keck II/NIRSPEC
- Wind of Change: retrieving exoplanet atmospheric winds from high-resolution spectroscopy
- Constraints on Metastable Helium in the Atmospheres of WASP-69b and WASP-52b with Ultra-Narrowband Photometry
- An Overabundance of Low-density Neptune-like Planets
- A Heavy Molecular Weight Atmosphere for the Super-Earth π Men c
- Modelling the He I triplet absorption at 10830 Angstroms in the atmospheres of HD 189733 b and GJ 3470 b
- Metastable Helium Reveals an Extended Atmosphere for the Gas Giant HAT-P-18b
- The dichotomy of atmospheric escape in AU Mic b
- The high-energy environment and atmospheric escape of the mini-Neptune K2-18 b
- GJ 436b and the stellar wind interaction: simulations constraints using Ly and H transits
- Search for helium in the upper atmosphere of the hot Jupiter WASP-127 b using Gemini/Phoenix
- Zodiacal Exoplanets in Time. XI. The Orbit and Radiation Environment of the Young M Dwarf-Hosted Planet K2-25b
- Lyman-alpha in the GJ 1132 System: Stellar Emission and Planetary Atmospheric Evolution
- The XUV irradiation and likely atmospheric escape of the super-Earth Men c