The absorption of stellar radiation observed by the HD209458b in resonant lines of OI and CII has not yet been satisfactorily modeled. In our previous 2D simulations we have shown that the hydrogen-dominated upper atmosphere of HD209458b, heated by XUV radiation, expands supersonically beyond the Roche lobe and drags the heavier species along with it. Assuming solar abundances, OI and CII particles accelerated by tidal forces to velocities up to 50 km/s should produce the absorption due to Doppler resonance mechanism at the level of 6-10%, consistent with the observations. Since the 2D geometry does not take into account the Coriolis force in the planet reference frame, the question remained to which extent the spiraling of the escaping planetary material and its actually achieved velocity may influence the conclusions made on the basis of 2D modeling. In the present paper we apply for the first time in the study of HD209458b a global 3D hydrodynamic multi-fluid model that self-consistently describes the formation and expansion of the escaping planetary wind, affected by the tidal and Coriolis forces, as well as by the surrounding stellar wind. The modeling results confirm our previous findings that the velocity and density of the planetary flow are sufficiently high to produce the absorption in HI, OI, and CII resonant lines at the level close to the in-transit observed values. The novel finding is that the matching of the absorption measured in MgII and SiIII lines requires at least 10 times lower abundances of these elements than the Solar system values.

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.

We present two-dimensional multi-fluid numerical modelling of the upper atmosphere of the hot Jupiter WASP-12b. The model includes hydrogen chemistry, and self-consistently describes the expansion of the planetary upper atmosphere and mass loss due to intensive stellar irradiation, assuming a weakly magnetized planet. We simulate the planetary upper atmosphere and its interaction with the stellar wind (SW) with and without the inclusion of tidal force and consider different XUV irradiation conditions and SW parameters. With the inclusion of tidal force, even for a fast SW, the escaping planetary material forms two streams, propagating towards and away from the star. The atmospheric escape and related mass loss rate reaching the value of 10^12 gs^-1 appear to be mostly controlled by the stellar gravitational pull. We computed the column density and dynamics of MgII ions considering three different sets of SW parameters and XUV fluxes. The simulations enable to compute the absorption at the position of the Mg h line and to reproduce the times of ingress and egress. In case of a slow SW and without accounting for tidal force, the high orbital velocity leads to the formation of a shock approximately in the direction of the planetary orbital motion. In this case, mass loss is proportional to the stellar XUV flux. At the same time, ignoring of tidal effects for WASP-12b is a strong simplification, so the scenario with a shock, altogether is an unrealistic one.

We present a 3D fully selfconsistent multi-fluid hydrodynamic aeronomy model to study the structure of a hydrogen dominated expanding upper atmosphere around the hot Jupiter HD 209458b and the warm Neptune GJ 436b. In comparison to previous studies with 1D and 2D models, the present work finds such 3D features as zonal flows in upper atmosphere reaching up to 1 km/s, the tilting of the planetary outflow by Coriolis force by up to 45 degrees and its compression around equatorial plane by tidal forces. We also investigated in details the influence of Helium (He) on the structure of the thermosphere. It is found that by decrease of the barometric scale-height, the He presence in the atmosphere strongly affects the H2 dissociation front and the temperature maximum.

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.