Energetic neutral atoms as the explanation for the high velocity hydrogen around HD 209458b

Energetic neutral atoms as the explanation for the high velocity hydrogen around HD 209458b

M. Holmström, A. Ekenbäck11footnotemark: 1, F. Selsis , T. Penz, H. Lammer, and P. Wurz Swedish Institute of Space Physics, PO Box 812, SE-98128 Kiruna, Sweden.Centre de Recherche Astrophysique de Lyon (CNRS, Université de Lyon, Ecole Normale Supérieure de Lyon), 46 Allée d’Italie, F-69007, Lyon, France.Laboratoire d’Astrophysique de Bordeaux (CNRS, Université Bordeaux 1), BP 89, F-33270, Floirac, France.INAF – Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, I-90134 Palermo, Italy.Space Research Institute, Austrian Academy of Sciences, Schmiedlstr. 6, A-8042, Graz, Austria.Physikalisches Institut, University of Bern, Sidlerstr. 5, CH-3012 Bern, Switzerland.
December 7, 2007

Absorption in the stellar Lyman- (Ly-) line observed during the transit of the extrasolar planet HD 209458b reveals high velocity atomic hydrogen at great distances from the planet[1, 2]. This has been interpreted as hydrogen atoms escaping from the exosphere of the planet[1, 3], possibly undergoing hydrodynamic blow-off[4], being accelerated by stellar radiation pressure. However, around solar system planets the production of energetic neutral atoms from charge exchange between solar wind protons and neutral hydrogen from the exospheres has been observed[5, 6, 7], and should also occur at extrasolar planets. Here we show that the measured transit-associated Ly- absorption can be explained by the interaction between the exosphere of HD 209458b and the stellar wind, and that radiation pressure alone cannot explain the observation. This is the first observation of energetic neutral atoms outside the solar system. Since the stellar wind protons are the source of the observed energetic neutral atoms, this provides a completely new method of probing stellar wind conditions, and our model suggests a slow and hot stellar wind near HD 209458b at the time of the observation.

Energetic neutral atoms (ENAs) are produced wherever energetic ions meet a neutral atmosphere, and solar wind ENAs have been observed at every planet in the solar system where ENA instrumentation has been available — at Earth[5], at Mars[6], and at Venus[7].

By energetic we mean that the ions have a much greater velocity than the thermal velocities of the exospheric neutrals. During the charge exchange process, an electron is transferred from the neutral to the ion, resulting in a neutral atom and an ionized neutral. Due to the large relative velocities of the ions and the exospheric neutrals, the momenta of the individual atoms are preserved to a good approximation. Thus, the produced ENAs will have the same velocity distribution as the source population of ions.

When first observed (also by their Ly- signature[8, 9]), the extended hydrogen coronae of Mars and Venus were assumed to constitute the uppermost layers of an escaping exosphere. The observed densities were used to infer exospheric scale heights and temperatures, which proved to be extremely high compared to theoretical predictions (up to 700 K). In situ spacecraft observations later found exospheric temperatures of 210 and 270 K[10]. The discrepancy was eventually explained by photochemically-produced energetic particles, and by ENAs, produced by charge exchange between energetic solar wind protons and the planetary exosphere. This mechanism, well-known in the solar system, has however not been considered as a possible origin of the atomic hydrogen corona revealed by HST observations of HD 209458b.

HD 209458b is a Jupiter-type gas giant with a mass of 0.65 M and a size of 1.32 R that orbits around its host star HD 209458 at 0.045 AU[12], which is a solar-like G-type star with an age of about 4 Gyr. The activity of the star can be estimated from its X-ray luminosity measured by the XMM-Newton space observatory, and is comparable to that of the present Sun during a moderately quiet phase[11]. Because of its Sun-like stellar type and average activity, it is justified to use the energy environment observed at the Sun as inputs for our model.

For a first estimate of the ENA production near HD 209458b, we assume that the charge exchange takes place in an undisturbed stellar wind that is flowing radially away from the star. At 0.045 AU from HD 209458, the stellar wind is most likely subsonic[13] and does not produce a planetary bow shock. Simulations indicate a subsolar magnetopause distance of about four planetary radii if the planet is magnetized[14]. If the planet is not magnetized, we would expect the undisturbed stellar wind to get even closer to the planet. Here we model the ENA production by a particle model that includes stellar wind protons and atomic hydrogen. Charge exchange between stellar wind protons and exospheric hydrogen atoms takes place outside a conic obstacle that represents the magnetosphere of the planet (Supplementary Fig. 1). The resulting exospheric cloud, along with the produced ENAs, covers a region larger than the stellar disc, as seen from Earth (Fig. 1). The cloud is shaped like a comet tail due to the stellar radiation pressure, curved by the Coriolis force, as predicted[15] and seen in earlier numerical simulations[1, 3]. There is a population of atoms with high velocity — these are the stellar wind protons that have charge exchanged, becoming ENAs. In the velocity spectrum along the -axis (the planet–star line), the ENAs are clearly visible as a distribution that is separate from the main exospheric hydrogen component, due to the different bulk velocities and temperatures (Fig. 2).

Now we estimate how the ENA cloud would affect the observed Ly- absorption spectrum of HD 209458b[1]. The line profile was observed outside and during transit, and the difference between the two profiles correspond to the attenuation by hydrogen atoms (Fig. 3).

There are several features of the transit spectrum that any proposed source of the observed hydrogen atoms need to account for. First, hydrogen atoms with velocities of up to -130 km/s (away from the star). Second, a fairly uniform absorption over the whole velocity range -130 to -45 km/s. Third, absorption in the velocity range between 30 and 105 km/s (toward the star).

The current explanation of the observation is that hydrogen atoms in the exosphere are undergoing hydrodynamic escape, and are then further accelerated by the stellar radiation pressure[1, 4]. There are however some difficulties in explaining the observations by this process, as can be seen by examining the three features listed above.

First of all, a large radiation pressure on the hydrogen atoms is needed to accelerate them to a velocity of 130 km/s before they are photoionized. The acceleration also has to occur before they move out from the region in front of the star, due to the orbital motion of the planet. The second feature is difficult to explain, since if hydrogen atoms were driven to speeds of up to 130 km/s, we would expect the velocity spectrum to have an exponential decay for higher velocities, due to the finite lifetime of hydrogen atoms because of photoionization (four hours on average). This drop-off for high velocities is independent of the details of the model, e.g., the values of radiation pressure and photoionization lifetime used. This would lead to a decay in the absorption spectrum, inconsistent with the observed fairly uniform absorption over the whole velocity range -130 to -45 km/s. Finally, an exosphere driven by radiation pressure cannot explain hydrogen atoms moving toward the star with speeds between 30 and 105 km/s. However, this feature of the observation is not completely certain, and it has been stated[1] that more observations are needed to clarify whether an absorption is present in the red part of the line (toward the star).

Our model shows that all observed features listed above can be explained by ENAs. If we turn off the ENA production in the model, none of these features are explained. When we compare the modeled Ly- profile with the observed ones, we find that the modeled spectrum leads to attenuation over the whole velocity range from -130 to -45 km/s, as is seen in the observation. The model also shows some absorption in the red part of the velocity spectrum, i.e. hydrogen atoms moving at high velocities toward the star, since for this stellar wind (50 km/s and  K), some part of the proton velocity distribution will have positive velocities along the -axis (toward the star), resulting in an ENA flux toward the star. This slow and hot stellar wind is not unrealistic at such small orbital distances[13].

If ENAs from charge exchange are responsible for the observed attenuation, we have on one hand, much less information on the main exospheric component than suggested by previous explanations. Since what we observe are mainly ENAs, a range of exospheric conditions and atmospheric loss rates can be consistent with the observation. Thus, the observation only constrains the radiation pressure driven atmospheric escape insofar that the exosphere has to be extended enough to reach the stellar wind outside the magnetopause of the planet. The atmospheric escape through ENA production is small. For the model parameters used here, the loss is  kg/s, which is more than an order of magnitude smaller than the estimated thermal loss of about  kg/s for similar exospheric conditions[16, 17, 18].

On the other hand, we gain information on the underlying plasma flows, and if it is the undisturbed stellar wind, we have a way of observing stellar wind properties such as temperature and velocity around other stars, at the location of extrasolar planets. By varying the stellar wind temperature, the stellar wind velocity and the radiation pressure in the model, we find a best fit of the modeled Ly- absorption to the observation for a stellar wind velocity of 50 km/s, and a temperature of 10 K (Fig. 3).

Although, for HD 209458b, the available Ly- data are affected by large uncertainties, more accurate observations would improve the derived stellar wind estimates. Also, the depletion of the stellar wind proton flow by charge exchange will change the character of the planet–stellar wind interaction. Present models of the stellar wind interaction with HD 209458b have not taken this process into account[14]. More observations of the Ly- absorption by HD 209458b, and its variation over time, could be used to confirm the origin of the extended hot hydrogen cloud. The variability of the Ly- by ENAs should be larger on short time scales than for other explanations, e.g., hydrodynamic escape, since the solar wind parameters can change significantly on a timescale of hours, as seen near Earth.


  • [1] Vidal-Madjar, A., et al., An extended upper atmosphere around the extrasolar planet HD209458b, Nature, 422, 143-146 (2003)
  • [2] Ben-Jaffel, L., Exoplanet HD 209458b: Inflated hydrogen atmosphere but no sign of evaporation, The Astrophysical Journal, 671, L61-L64 (2007)
  • [3] Vidal-Madjar, A. & Lecavelier des Etangs, A., ”Osiris”(HD209458b), an evaporating planet, in Extrasolar Planets: Today and Tomorrow, ASP Conference Proceedings, Vol. 321. Edited by J.-P. Beaulieu, A. Lecavelier des Etangs, & C. Terquem, 152-159 (2004). ISBN: 1-58381-183-4
  • [4] Vidal-Madjar, et al., Detection of Oxygen and Carbon in the Hydrodynamically Escaping Atmosphere of the Extrasolar Planet HD 209458b, The Astrophysical Journal, 604(1), L69-L72 (2004)
  • [5] Collier, M. R., et al., Observations of neutral atoms from the solar wind, Journal of Geophysical Research, 106(A11), 24893-24906 (2001)
  • [6] Futaana, Y., et al., First ENA observations at Mars: Subsolar ENA jet, Icarus, 182(2), 413–423 (2006)
  • [7] Galli, A., et al., First observation of energetic neutral atoms in the Venus environment, Planetary and Space Science, in press (2008) http://dx.doi.org/10.1016/j.pss.2007.12.011
  • [8] Barth, C. A., Interpretation of the Mariner 5 Lyman alpha measurements, J. Atmos. Sci., 25, 564-567 (1968)
  • [9] Barth, C. A., et al., Marine 6: Ultraviolet Spectrum of Mars Upper Atmosphere, Science, 165, 1004-1005 (1969)
  • [10] Lichtenegger, H.I.M., et al., Effects of low energetic neutral atoms on martian and venusian dayside exospheric temperature estimations, Space Science Reviews, 126, 469­-501 (2006)
  • [11] Penz, T., Micela, G., & Lammer, H., Influence of the evolving stellar X–ray luminosity distribution on exoplanetary mass loss, Astron. Astrophys. 477, 309-314 (2008)
  • [12] Knutson, H., Charbonneau, D., Noyes, R. W., Brown, T.M. & Gilliland, R.L., Using stellar limb-darkening to refine the properties of HD 209458b, Astrophys. J., 655, 564-575 (2007)
  • [13] Preusse, S., Kopp, A., Büchner, J. & Motschmann, U., Stellar wind regimes of close-in extrasolar planets, Astronomy and Astrophysics, 434(3), 1191-1200 (2005)
  • [14] Preusse, S., Kopp, A., Büchner, J. & Motschmann, U., MHD simulation scenarios of the stellar wind interaction with Hot Jupiter magnetospheres, Planetary and Space Science, 55(5), 589-597 (2007)
  • [15] Schneider, J., Rauer, H., Lasota, J. P., Bonazzola, S., & Chassefiere, E., Brown Dwarfs and Extrasolar Planets, ASP Conference Series 134, 241-244 (1998)
  • [16] Yelle, Roger V., Aeronomy of extra-solar giant planets at small orbital distances, Icarus, 170(1), 167-179 (2004)
  • [17] Yelle, Roger V., Corrigendum to Aeronomy of extra-solar giant planets at small orbital distances, Icarus, 183(2), 508 (2006)
  • [18] García Muñoz, A., Physical and chemical aeronomy of HD 209458b, Planetary and Space Science, 55(10), 1426–1455 (2007)

Supplementary Information is linked to the online version of the paper at www.nature.com/nature

M.H., H.L., F.S. and T.P. thanks the International Space Science Institute (ISSI), as this study was carried out within the framework of the ISSI team “Evolution of Exoplanet Atmospheres and their Characterization”. T.P. is supported by the Marie Curie Fellowship project ISHERPA, and the host institution INAF-Osservatorio Astronomico di Palermo. H.L. thanks ASA for funding the CoRoT project. This research was conducted using the resources of the High Performance Computing Center North (HPC2N), Umeå University, Sweden, and the Center for Scientific and Technical Computing (LUNARC), Lund University, Sweden. The software used in this work was in part developed by the DOE-supported ASC / Alliance Center for Astrophysical Thermonuclear Flashes at the University of Chicago.

Author Contributions
A.E. wrote an initial version of the simulation code. F.S. helped in modeling the observation. T.P. provided knowledge on the atmospheres of extrasolar planets. H.L. and F.S. suggested that the HST observation could be due to ENAs. P.W. contributed with expertise in ENA processes.

Author Information
Reprints and permissions information is available at www.nature.com/reprints. The authors declare no competing financial interests. Correspondence and requests for materials should be addressed to M.H. (matsh@irf.se).

Figure 1: The hydrogen cloud around the planet. Shown from (a) above, perpendicular to the planet’s orbital plane, and (b) from Earth, along the -axis direction. Each point corresponds to a hydrogen meta particle. The color of the points shows the velocity of the particles along the -axis. Particles with velocity magnitude smaller than 50 km/s are red, and those with higher velocity are black. The small circles show the planet size. The large circle in (b) shows the star’s position at mid transit. During transit the star moves from left to right in (b). At the outer boundaries of the simulation domain, stellar wind protons are injected with a  cm number density, 50 km/s velocity, and  K temperature. The planet’s interaction with the stellar wind is modeled by removing all stellar wind protons inside a conic obstacle at a sub-stellar distance of about (where the radius of the planet  m). Hydrogen atoms are launched from an inner boundary (a sphere of radius ) assuming a number density of  cm and a temperature of  K, consistent with atmospheric models[18]. The trajectory of each proton and hydrogen atom is followed in time. The forces on a hydrogen atom are the gravity of the planet, the Coriolis force due to the rotating coordinate system, and radiation pressure. After each time step a hydrogen atom can undergo photoionization, elastic collision with another hydrogen atom or charge exchange with a proton. The photoionization time assumed is 4 hours which is a scaled Earth value. The radiation pressure corresponds to a photon-hydrogen collision rate of 0.35 s and is chosen to improve the model fit. It is lower than a scaled Earth value of  s over a solar cycle. The coordinate system used is centered at the planet with its -axis toward the star, and the -axis opposite to the planet’s velocity. Further details of the simulations can be found in the SI.
Figure 2: Velocities of the hydrogen atoms. The modeled -axis (planet–star) velocity spectrum of hydrogen atoms in front of the star at the moment of mid transit, not including atoms in front or behind the planet. The part of the distribution that is due to ENAs is shaded. Varying the stellar wind temperature and velocity in the model confirms that the width of this part of the distribution is proportional to the temperature of the stellar wind, with a larger width for larger temperatures, and the center of the distribution follows the stellar wind velocity. The un-shaded part of the spectrum is due to the exospheric hydrogen atoms.
Figure 3: Comparison of the modeled Ly- profile with the observed ones. In blue is the observed profile before transit. In green is the observed profile during transit. In red is the modeled profile, constructed by applying the attenuations computed from the simulations to the observed profile before transit. The abscissa is the hydrogen velocity along the -axis (away from Earth — toward the star). The regions where there is a significant difference between the profiles are denoted ’In’ and ’Geo’, the latter being the region of geocoronal emission at low velocities that should be excluded. The modeled profile is computed at the instant of mid-transit. The details of computing the Ly- attenuation from the hydrogen cloud are given in the SI. The modeled Ly- absorption shown here is for a stellar wind velocity of 50 km/s and a temperature of 10 K. The fit is worse for stellar wind velocities of 0 or 100 km/s, or stellar wind temperatures of  K or  K as shown in the SI.

Figure 1

Figure 2

Figure 3

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
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

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 description