Reconnaissance of the TRAPPIST-1 exoplanet system in the Lyman-\alpha line

Reconnaissance of the TRAPPIST-1 exoplanet system in the Lyman- line

V. Bourrier Observatoire de l’Université de Genève, 51 chemin des Maillettes, 1290 Sauverny, Switzerland    D. Ehrenreich Observatoire de l’Université de Genève, 51 chemin des Maillettes, 1290 Sauverny, Switzerland    P.J. Wheatley Department of Physics, University of Warwick, Coventry CV4 7AL, UK    E. Bolmont Laboratoire AIM Paris-Saclay, CEA/DRF - CNRS - Univ. Paris Diderot - IRFU/SAp, Centre de Saclay, F- 91191 Gif-sur-Yvette Cedex, France    M. Gillon Institut d’Astrophysique et de Géophysique, Université de Liège, Allée du 6 Aout 19C, 4000 Liège, Belgium    J. de Wit Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA    A.J Burgasser Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA 92093, USA    E. Jehin Institut d’Astrophysique et de Géophysique, Université de Liège, Allée du 6 Aout 19C, 4000 Liège, Belgium    D. Queloz Cavendish Laboratory, J J Thomson Avenue, Cambridge, CB3 0HE, UK    Observatoire de l’Université de Genève, 51 chemin des Maillettes, 1290 Sauverny, Switzerland    A.H.M.J. Triaud Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK vincent.bourrier@unige.ch
Key Words.:
planetary systems - stars: individual: TRAPPIST-1 - techniques: spectroscopic
offprints: V.B. ()

The TRAPPIST-1 system offers the opportunity to characterize terrestrial, potentially habitable planets orbiting a nearby ultracool dwarf star. We performed a four-orbit reconnaissance with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope to study the stellar emission at Lyman-, to assess the presence of hydrogen exospheres around the two inner planets, and to determine their UV irradiation. We detect the Lyman- line of TRAPPIST-1, making it the coldest exoplanet host star for which this line has been measured. We reconstruct the intrinsic line profile, showing that it lacks broad wings and is much fainter than expected from the stellar X-ray emission. TRAPPIST-1 has a similar X-ray emission as Proxima Cen but a much lower Ly- emission. This suggests that TRAPPIST-1 chromosphere is only moderately active compared to its transition region and corona. We estimated the atmospheric mass loss rates for all planets, and found that despite a moderate extreme UV emission the total XUV irradiation could be strong enough to strip the atmospheres of the inner planets in a few billions years. We detect marginal flux decreases at the times of TRAPPIST-1b and c transits, which might originate from stellar activity, but could also hint at the presence of extended hydrogen exospheres. Understanding the origin of these Lyman- variations will be crucial in assessing the atmospheric stability and potential habitability of the TRAPPIST-1 planets.

1 Introduction

For close-in planets, heating by the stellar energy can lead to the expansion of the upper atmospheric layers and their eventual escape. Transit spectroscopy in the UV is a powerful way to probe these extended atmospheres, in particular in the stellar Lyman- (Ly-) line of neutral hydrogen (e.g., VM2003, Lecav2012). To date, the warm Neptune GJ436b is the lowest-mass planet found evaporating (Ehrenreich2015), with transit absorption depths up to 60% in the Ly- line. The formation of such an extended exosphere is due in great part to the low mass of GJ436b and the gentle irradiation from its M-dwarf host (Bourrier2015_GJ436, Bourrier2016), suggesting that small planets around cold stars can produce very deep UV transit signatures. For telluric planets, large amounts of hydrogen in their upper atmosphere could indicate the presence of a steam envelope being photo-dissociated, and replenished by evaporating water oceans (Jura2004). Few attempts have been made, however, to detect atmospheric escape from Earth-sized planets, and Ly- transit observations of the super Earth 55 Cnc e (Ehrenreich2012) and HD 97658 b (Bourrier2016_HD976) showed no evidence for hydrogen exospheres.
The nearby ultracool dwarf TRAPPIST-1 was found to harbor seven short-period, transiting, Earth-sized planets (Gillon2016; Gillon2017). TRAPPIST-1 e, f, and g are within the habitable zone, and all planets have equilibrium temperatures low enough to make liquid water possible on their surface. The combined transmission spectrum of the inner planets, b and c, in the near-infrared (deWit2016) revealed no significant features, which rules out cloud-free hydrogen-dominated atmospheres but still allows for denser envelopes, such as a cloud-free water-vapor atmosphere. The high level of X-rays to ultraviolet (XUV, in 5–912 Å) irradiation (Wheatley2016) and the sustained activity of M dwarfs over approximately one billion years could however impact the putative atmospheres of the TRAPPIST-1 planets and their habitability (Miguel2015, Bolmont2016).
The system proximity (12 pc), its large planet-to-star radii ratios, and high systemic velocity (-56 km s) make it an interesting target for Ly- transit spectroscopy. Little is known about the FUV emission of such cold stars, and we present in this Letter a reconnaissance study of TRAPPIST-1 with the Hubble Space Telescope (HST). Our objectives were to try and measure the stellar Ly- emission, to estimate the extreme UV (EUV, in 100–912 Å) irradiative environment of the planets, and to search for signatures of planetary hydrogen escape.

2 The Ly- line of TRAPPIST-1

We used the Space Telescope Imaging Spectrograph (STIS) instrument onboard the HST to observe TRAPPIST-1 at FUV wavelengths (G140M grating at 1222 Å). Four HST orbits were obtained in 2016 in the frame of the Mid-Cycle Program 14493 (PI: Bourrier), with exposure time of approximately 35 min. We observed TRAPPIST-1 during two consecutive HST orbits on 26 September, outside of the two inner planets’ transits, to search for the stellar Ly- emission and obtain a reference for the out-of-transit flux. We observed the star during the transit of TRAPPIST-1b on September 30th (one HST orbit) and approximately 1.7 h after the transit of TRAPPIST-1c on November 23rd (one HST orbit). The data were reduced with the CALSTIS pipeline. In the region of the Ly- line, the background is dominated by Earth geocoronal airglow emission. The airglow varies in strength and position with the epoch of observation, and its correction by the pipeline can yield spurious flux values when the airglow is much stronger than the stellar line. We thus excluded from our analysis the contaminated ranges [-60 ; 74] km s (out-of-transits), [-60 ; 49] km s (TRAPPIST-1b), and [-60 ; 54] km s (TRAPPIST-1c), defined in the heliocentric rest frame (Fig. 1). No breathing effect was detected in the time-tagged data (Bourrier2016_HD976), most likely because it is dominated by the photon noise from the very faint Ly- line.

Stellar Ly- lines are difficult to observe because of absorption by the interstellar medium (ISM) between the star and the Earth. However, the proximity of the TRAPPIST-1 system (12 pc) limits the depth of this absorption, and the high systemic velocity of the star (-56.3 km s; Reiners2015) shifts the stellar line blueward of the spectrum. This allowed us to measure the blue wing and the core of TRAPPIST-1 Ly- line, freed up from both ISM and airglow contamination (Fig. 1). The stellar line is spectrally resolved with STIS, and measured with similar flux levels in our four epochs of observations. To our knowledge, TRAPPIST-1 is the coldest exoplanet host star from which UV emission has been detected, and the coldest star with a resolved observation of the Lyman- line. We discuss the importance of this detection in Sect. 3.

We reconstructed the intrinsic Ly- profile using the method detailed in Bourrier2015_GJ436. A model profile of the stellar Ly- line is absorbed by the interstellar hydrogen and deuterium, and convolved by STIS line spread function (LSF). A Metropolis-Hasting Markov chain Monte Carlo algorithm is used to find the best fit between the model spectrum and the observed spectrum, which are compared over the range shown in Fig. 2. The observed spectrum was calculated as the mean of all spectra, excluding the ranges contaminated by airglow emission or displaying possible flux variations (Sec. 5). The heliocentric velocity of the star was fixed to its aforementioned value. The LISM Kinematic Calculator111http://sredfield.web.wesleyan.edu/ predicts that the line of sight (LOS) toward TRAPPIST-1 crosses the Local Interstellar Cloud (LIC). Since the fit is mostly sensitive to the hydrogen column density in the cloud, we fixed its heliocentric radial velocity, temperature, and turbulent velocity to the values of the LIC (-1.25 km s, 7500 K, 1.62 km s, respectively; Redfield_Linsky2008). Trying out different profiles for the intrinsic stellar Ly- line, we found it was best represented as a single-peaked Gaussian profile, displayed in Fig. 2. The best-fit hydrogen column density for the ISM is log(H i) = 18.30.2, similar to other LOS at a distance of 12 pc (Fig. 14 in Wood2005).

Figure 1: Ly- line spectrum of TRAPPIST-1, plotted as a function of Doppler velocity in the heliocentric rest frame. The black line has been averaged over the out-of-transit observations, while the blue and orange spectra were obtained at the times of TRAPPIST-1b and c transits, respectively. The dashed black line indicates the velocity of the star, and the dashed red line the velocity of the ISM/LIC cloud. The red wing of the Ly- line cannot be observed from Earth because of ISM absorption along the line of sight, and contamination from geocoronal emission (hatched region).
Figure 2: Intrinsic Ly- line profile of TRAPPIST-1 (orange line). The solid blue line shows the Ly- line profile after absorption by the ISM hydrogen and deuterium (dotted black line). The solid black line shows the line profile convolved with STIS LSF and compared to the observations shown as a red histogram. Hatched regions were excluded from the fit. The right axis shows the ratio between radiation pressure and stellar gravity.

3 High-energy emission from an ultracool dwarf

The X-ray emission of TRAPPIST-1 was recently studied by Wheatley2016, who found that the star is a relatively strong and variable coronal X-ray source with an X-ray luminosity similar to that of the quiet Sun, despite a much lower bolometric luminosity. We report the weighted average of their APEC- and CEMEKL-derived fluxes in 5–100 Å  in Table 1. We searched the GALEX data222http://galex.stsci.edu/galexview/ for near-UV (NUV) or far-UV (FUV) observations of TRAPPIST-1, but the closest star with a detection is 57” away. The stellar EUV emission would, in any case, be absorbed by the ISM, but it can be estimated through indirect methods. We used the Linsky2014 relation for M dwarfs, based on our estimate for the intrinsic Ly- flux, to calculate the flux in 100–912 Å  (Table 1). We find a much lower EUV flux than derived by Wheatley2016 from the X-ray flux and Chadney2015 empirical relation, possibly due to the very late type of TRAPPIST-1.

Little is known about the UV emission of low-mass, cool stars such as TRAPPIST-1 (Table 2). Because of their low temperature, the photosphere of M dwarfs contributes little to the UV continuum, compared to bright chromospheric and transition region emission lines. Ly- in particular is the brightest UV line of low-mass stars (e.g., France2013). Assuming the Linsky2014 relations are relevant to such a cold star, the Ly- flux of TRAPPIST-1 is indeed as large as its entire EUV flux. In Table 2, we compare TRAPPIST-1 with the coldest M dwarf hosts with a measured Ly- flux. Ly- emission decreases with stellar temperature (Linsky2013). TRAPPIST-1 is more than 400 K colder than those stars, hence it is not surprising that its Ly- flux is more than twice as small. Interestingly, TRAPPIST-1 has a similar X-ray emission as Proxima Cen (0.136 and 0.142 erg s cm at 1 au, respectively) but a much lower Ly- emission. This might be an indication of a lower chromospheric activity for TRAPPIST-1 compared to its corona, and is consistent with the decrease in the ratio of Ly- to X-ray emission noted for later-type stars by Linsky2013. Furthermore, XUV activity as a function of age has been shown to remain constant for M stars younger than approximately 300 Myr (Shkolnik2014). After this saturation phase, chromospheric activity and coronal emission steadily decrease due to the spin down of the star, but the decrease in X-ray emission is steeper than in the NUV and FUV (Stelzer2013). The fact that TRAPPIST-1 emits nearly three times less flux at Ly- than in the X-ray would thus suggest it is still relatively young (age is constrained to be more than 500 Myr, Filippazzo2015).

F F F
(5-100 Å) (100-912 Å)
T-1b 1102 400 414 7.6 46
T-1c 588 213 221 4.1 14
T-1d 296 107 111 2.0 9
T-1e 171 62 64 1.2 6
T-1f 99 36 37 0.7 5
T-1g 67 24 25 0.5 2
T-1h 34 12 13 0.2 -
  • Note: Fluxes in erg s cm, hydrogen photoionization rate in x10 s, mass loss rate in 10 g s. X-ray flux is derived from Wheatley2016. The mass of TRAPPIST-1 h is unconstrained, and we did not estimate its mass-loss rate. Uncertainties on the EUV flux were estimated by accounting for both the uncertainty on the Ly- flux and the Linsky2014 relation.

Table 1: XUV irradiation of TRAPPIST-1 planets.
Star Type T F
(K) (erg s cm)
TRAPPIST-1 M8 2550 0.05
GJ 1214 M4.5 2935 0.12
Proxima Cen M5.5 3050 0.30
GJ 876 M5 3062 0.37
  • Note: Fluxes are given at 1 au from the stars. References: Gillon2016; Loyd2016; Youngblood2016; Anglada2016; Wood2005.

Table 2: Coldest M dwarf hosts with measured Ly- flux

4 Mass loss from the TRAPPIST-1 planets

While the TRAPPIST-1 planets are unlikely to be surrounded by hydrogen-dominated atmospheres (deWit2016), they may still harbor a substantial amount of water (Gillon2017). The stellar XUV irradiation and internal tidal heating could be sufficient for water oceans to form a supercritical steam atmosphere that could be photodissociated and sustain a significant hydrogen escape (Jura2004). Using the derived XUV flux, we calculated the mass-loss rate from all TRAPPIST-1 planets in the energy-limited regime (Lecav2007), updating the estimates made by Bolmont2016 and Wheatley2016:

(1)

with representing the heating efficiency set to 1% (Ehrenreich2015, Salz2016a). accounts for the increased cross-sectional area of planets to EUV radiation, while accounts for the contribution of tidal forces to the potential energy (Erkaev2007). Both are set to unity for these cool and small planets. The mean density of the planets was calculated using preliminary masses derived from transit timing variations (Gillon2017). The derived escape rates are given in Table 1. Earth-like atmospheres (510 g) and oceans (1.410 g) would take between 5 and 22 Gy to be stripped from TRAPPIST-1d to g, but only approximately 1 and 3 Gy for TRAPPIST-1b and c. While this result strongly depends on the heating efficiency and neglects the complex physics behind atmospheric escape (Bolmont et al. 2017, in prep), it suggests nonetheless that Earth-like TRAPPIST-1b and c would have to be still young - or to have an initial composition richer in water than Earth - to have retained a significant water content.

5 Hydrogen exospheres around TRAPPIST-1b/c

Using the out-of-transit Ly- line as a reference, we identified marginal flux decreases (Fig. 3) during the transit of TRAPPIST-1b (4021% in [-95 ; -55] km s) and after the transit of TRAPPIST-1c (4118% in [-135 ; -40] km s). Since the star has a variable corona (Wheatley2016), this might be an indication of a similarly variable chromosphere. Alternatively, and given that the Ly- line is stable over time outside of the above ranges (Fig. 3), this might hint at the presence of hydrogen exospheres around the two inner planets. The measured variations are too large to be caused by the planetary disks (Llama2016) and would correspond to opaque exospheric disks extending up to approximately seven times the planets’ radii. We note that the putative exosphere of TRAPPIST-1c would absorb the stellar Ly- line up to two hours after the planet transit, at larger velocities than the signal from TRAPPIST-1b. This suggests that the exosphere of planet c would extend into a cometary tail, with the coma around the planet likely producing a much deeper absorption at the time of the transit.
The velocity field of these exospheres would depend on the planet orbital velocity and the relative strength of radiation pressure and stellar gravity (Bourrier2015_GJ436). The faint Ly- line of TRAPPIST-1 yields a radiation pressure that can only compensate for up to 20% of the star gravity (Fig. 2). This is reminiscent of the giant hydrogen exosphere of GJ 436 b (Ehrenreich2015), slowly diffusing around the planet owing to the low radiation pressure from its M dwarf host (70%). However, using the EVE code (Bourrier_lecav2013) to simulate atmospheric escape from TRAPPIST-1b and c, we found that the orbital velocities of the planets (80 and 70 km s, respectively) and the low radiation pressure could impart maximum Doppler velocities below 40 km s to hydrogen atoms at the time of the observations. If the absorption signatures from TRAPPIST1-b and -c are confirmed, this means that their exospheres are subjected to additional dynamical mechanisms (e.g., stellar wind interactions or a high-speed planetary outflow). Finally, we used the XUV spectrum derived in Sect. 2 to calculate the photoionization rate of neutral hydrogen at the semi-major axis of TRAPPIST-1b and c. We found low values that correspond for each planet to lifetimes of approximately 40 and 70 hours (Table 1). This ionization field, lower than for GJ 436 b (Bourrier2016), would further help the formation of giant hydrogen exospheres around the inner TRAPPIST-1 planets.

Figure 3: Ly- line flux integrated over [-95 ; -55] km s (blue points) and the complementary domain in [-210 ; 210] km s, airglow excluded (black points). The flux is normalized by the reference measurements in the left panel, and plotted as a function of time relative to the transit of TRAPPIST-1b (left and middle panels) and TRAPPIST-1c (right panel).

6 Conclusion

The TRAPPIST-1 system is a compelling target for atmospheric characterization, with seven temperate terrestrial planets. We detect the Ly- of its ultracool dwarf star, and derive the EUV irradiation of the planets. The Ly- line contributes as much as the entire EUV spectrum, and combined with an X-ray irradiation more than twice as strong, it could significantly impact the chemical and physical balance of the planets’ atmospheres. In particular, water vapor in the upper atmosphere could photo-dissociate and sustain an outflow of escaping hydrogen. The stellar Ly- line is bright enough to perform transit spectroscopy, and we detect marginal flux decreases in localized, high-velocity ranges during the transit of planet b, and shortly after the transit of planet c. This could hint at the presence of extended hydrogen exospheres around the two inner planets, and suggest that atmospheric escape might play an active role in the evolution of all TRAPPIST-1 planets. Even if the measured variations do not have a planetary origin and arise from activity in the stellar Ly- line, this would have important consequences for the stability of the planetary atmospheres and the production of biomarkers. More observations will be required to characterize the high-energy environment of the TRAPPIST-1 system, and its effect on the planets’ atmospheres.

Acknowledgements.
We thank the referee for their efficient reading of our paper. This work is based on observations made with the NASA/ESA Hubble Space Telescope, for which support was provided by NASA through a grant from the Space Telescope Science Institute. This work has been carried out in the frame of the National Centre for Competence in Research “PlanetS” supported by the Swiss National Science Foundation (SNSF), and was funded in part by the European Research Council (ERC) under the FP/2007-2013 ERC grant 336480. V.B. and D.E. acknowledges the financial support of the SNSF. M.G. and E.J. are Research Associate at the Belgian Fund for Scientific Research (F.R.S-FNRS). E. B. acknowledges funding by the European Research Council through ERC grant SPIRE 647383. PW is supported by a Science and Technology Facilities Council (STFC) consolidated grant (ST/L000733/1).

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
375325
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description