No Arabic abstract
Warm Neptune GJ3470b has been recently observed in 23S-23P transition of metastable helium, yielding absorption of about 1% in Doppler velocity range of [-40; 10] km/s. Along with previous detection of absorption in Ly{alpha} with depth of 20-40% in the blue and red wings of the line, it offers a complex target for simulation and testing of the current models. Obtained results suggest that absorption in both these lines comes from interaction of expanding upper planetary atmosphere with stellar plasma wind, allowing to constrain the stellar plasma parameters and the helium abundance in planet atmosphere.
Using a global 3D, fully self-consistent, multi-fluid hydrodynamic model, we simulate the escaping upper atmosphere of the warm Neptune GJ436b, driven by the stellar XUV radiation impact and gravitational forces and interacting with the stellar wind. Under the typical parameters of XUV flux and stellar wind plasma expected for GJ436, we calculate in-transit absorption in Ly{alpha} and find that it is produced mostly by Energetic Neutral Atoms outside of the planetary Roche lobe, due to the resonant thermal line broadening. At the same time, the influence of radiation pressure has been shown to be insignificant. The modelled absorption is in good agreement with the observations and reveals such features as strong asymmetry between blue and red wings of the absorbed Ly{alpha} line profile, deep transit depth in the high velocity blue part of the line reaching more than 70%, and the timing of early ingress. On the other hand, the model produces significantly deeper and longer egress than in observations, indicating that there might be other processes and factors, still not accounted, that affect the interaction between the planetary escaping material and the stellar wind. At the same time, it is possible that the observational data, collected in different measurement campaigns, are affected by strong variations of the stellar wind parameters between the visits, and therefore, they cannot be reproduced altogether with the single set of model parameters.
Transmission spectroscopy of WASP-107b revealed 7-8% absorption at the position of metastable HeI triplet at 10830 {AA} in Doppler velocity range of [-20; 10] km/s, which is stronger than that measured in other exoplanets. With a dedicated 3D self-consistent hydrodynamic multi-fluid model we calculated the expanding upper atmosphere of WASP-107b and reproduced within the observations accuracy the measured HeI absorption profiles, constraining the stellar XUV flux to (6-10) erg cm-2 s-1 at 1 a.u., and the upper atmosphere helium abundance He/H to 0.075-0.15. The radiation pressure acting on the metastable HeI atoms was shown to be an important factor affecting the shape of the absorption profiles. Its effect is counterbalanced by the processes of collisional depopulation of the HeI metastable state. Altogether, the observed HeI absorption in WASP-107b can be interpreted with the expected reasonable parameters of the stellar-planetary system and appropriate account of the electron and atom impact processes.
To shed more light on the nature of the observed Ly{alpha} absorption during transits of HD 209458b and to quantify the major mechanisms responsible for the production of fast hydrogen atoms (the so called energetic neutral atoms, ENAs) around the planet, 2D hydrodynamic multifluid modeling of the expanding planetary upper atmosphere, which is driven by stellar XUV, and its interaction with the stellar wind has been performed. The model selfconsistently describes the escaping planetary wind, taking into account the generation of ENAs due to particle acceleration by the radiation pressure and by the charge exchange between the stellar wind protons and planetary atoms. The calculations in a wide range of stellar wind parameters and XUV flux values showed that under typical Sun-like star conditions, the amount of generated ENAs is too small, and the observed absorption at the level of 6-8 percent can be attributed only to the non-resonant natural line broadening. For lower XUV fluxes, e.g., during the activity minima, the number of planetary atoms that survive photoionization and give rise to ENAs increases, resulting in up to 10-15 percent absorption at the blue wing of the Lya line, caused by resonant thermal line broadening. A similar asymmetric absorption can be seen under the conditions realized during coronal mass ejections, when sufficiently high stellar wind pressure confines the escaping planetary material within a kind of bowshock around the planet. It was found that the radiation pressure in all considered cases has a negligible contribution to the production of ENAs and the corresponding absorption.
Characterising the atmospheres of exoplanets is key to understanding their nature and provides hints about their formation and evolution. High-resolution measurements of the helium triplet, He(2$^{3}$S), absorption of highly irradiated planets have been recently reported, which provide a new mean to study their atmospheric escape. In this work, we study the escape of the upper atmospheres of HD 189733 b and GJ 3470 b by analysing high-resolution He(2$^{3}$S) absorption measurements and using a 1D hydrodynamic model coupled with a non-LTE model for the He(2$^{3}$S) state. We also use the H density derived from Ly$alpha$ observations to further constrain their temperatures, T, mass-loss rates,$dot M$, and H/He ratios. We have significantly improved our knowledge of the upper atmospheres of these planets. While HD 189733 b has a rather compressed atmosphere and small gas radial velocities, GJ 3470 b, with a gravitational potential ten times smaller, exhibits a very extended atmosphere and large radial outflow velocities. Hence, although GJ 3470 b is much less irradiated in the XUV, and its upper atmosphere is much cooler, it evaporates at a comparable rate. In particular, we find that the upper atmosphere of HD 189733 b is compact and hot, with a maximum T of 12400$^{+400}_{-300}$ K, with very low mean molecular mass (H/He=(99.2/0.8)$pm0.1$), almost fully ionised above 1.1 R$_p$, and with $dot M$=(1.1$pm0.1$)$times$10$^{11}$ g/s. In contrast, the upper atmosphere of GJ 3470 b is highly extended and relatively cold, with a maximum T of 5100$pm900$ K, also with very low mean molecular mass (H/He=(98.5/1.5)$^{+1.0}_{-1.5}$), not strongly ionised and with $dot M$=(1.9$pm1.1$)$times$10$^{11}$ g/s. Furthermore, our results suggest that the upper atmospheres of giant planets undergoing hydrodynamic escape tend to have very low mean molecular mass (H/He$gtrsim$97/3).
The GJ 436 planetary system is an extraordinary system. The Neptune-size planet that orbits the M3 dwarf revealed in the Ly$alpha$ line an extended neutral hydrogen atmosphere. This material fills a comet-like tail that obscures the stellar disc for more than 10 hours after the planetary transit. Here, we carry out a series of 3D radiation hydrodynamic simulations to model the interaction of the stellar wind with the escaping planetary atmosphere. With these models, we seek to reproduce the $sim56%$ absorption found in Ly$alpha$ transits, simultaneously with the lack of absorption in H$alpha$ transit. Varying the stellar wind strength and the EUV stellar luminosity, we search for a set of parameters that best fit the observational data. Based on Ly$alpha$ observations, we found a stellar wind velocity at the position of the planet to be around [250-460] km s$^{-1}$ with a temperature of $[3-4]times10^5$ K. The stellar and planetary mass loss rates are found to be $2times 10^{-15}$ M$_odot$ yr$^{-1}$ and $sim[6-10]times10^9$ g s$^{-1}$, respectively, for a stellar EUV luminosity of $[0.8-1.6]times10^{27}$ erg s$^{-1}$. For the parameters explored in our simulations, none of our models present any significant absorption in the H$alpha$ line in agreement with the observations.