Aims: We aim at constraining the conditions of the wind and high-energy emission of the host star reproducing the non-detection of Ly$\alpha$ planetary absorption. Methods: We model the escaping planetary atmosphere, the stellar wind, and their interaction employing a multi-fluid, three-dimensional hydrodynamic code. We assume a planetary atmosphere composed of hydrogen and helium. We run models varying the stellar high-energy emission and stellar mass-loss rate, further computing for each case the Ly$\alpha$ synthetic planetary atmospheric absorption and comparing it with the observations. Results: We find that a non-detection of Ly$\alpha$ in absorption employing the stellar high-energy emission estimated from far-ultraviolet and X-ray data requires a stellar wind with a stellar mass-loss rate about six times lower than solar. This result is a consequence of the fact that, for $\pi$ Men c, detectable Ly$\alpha$ absorption can be caused exclusively by energetic neutral atoms, which become more abundant with increasing the velocity and/or the density of the stellar wind. By considering, instead, that the star has a solar-like wind, the non-detection requires a stellar ionising radiation about four times higher than estimated. This is because, despite the fact that a stronger stellar high-energy emission ionises hydrogen more rapidly, it also increases the upper atmosphere heating and expansion, pushing the interaction region with the stellar wind farther away from the planet, where the planet atmospheric density that remains neutral becomes smaller and the production of energetic neutral atoms less efficient. Conclusions: Comparing the results of our grid of models with what is expected and estimated for the stellar wind and high-energy emission, respectively, we support the idea that the atmosphere of $\pi$ Men c is likely not hydrogen-dominated.

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 modeled the transit signatures in the Lya line of a putative Earth-sized planet orbiting in the HZ of the M dwarf GJ436. We estimated the transit depth in the Lya line for an exo-Earth with three types of atmospheres: a hydrogen-dominated atmosphere, a nitrogen-dominated atmosphere, and a nitrogen-dominated atmosphere with an amount of hydrogen equal to that of the Earth. We calculated the in-transit absorption they would produce in the Lya line. We applied it to the out-of-transit Lya observations of GJ 436 obtained by the HST and compared the calculated in-transit absorption with observational uncertainties to determine if it would be detectable. To validate the model, we also used our method to simulate the deep absorption signature observed during the transit of GJ 436b and showed that our model is capable of reproducing the observations. We used a DSMC code to model the planetary exospheres. The code includes several species and traces neutral particles and ions. At the lower boundary of the DSMC model we assumed an atmosphere density, temperature, and velocity obtained with a hydrodynamic model for the lower atmosphere. We showed that for a small rocky Earth-like planet orbiting in the HZ of GJ436 only the hydrogen-dominated atmosphere is marginally detectable with the STIS/HST. Neither a pure nitrogen atmosphere nor a nitrogen-dominated atmosphere with an Earth-like hydrogen concentration in the upper atmosphere are detectable. We also showed that the Lya observations of GJ436b can be reproduced reasonably well assuming a hydrogen-dominated atmosphere, both in the blue and red wings of the Lya line, which indicates that warm Neptune-like planets are a suitable target for Lya observations. Terrestrial planets can be observed in the Lya line if they orbit very nearby stars, or if several observational visits are available.

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.

This is the second paper in a series where we build a self-consistent model to simulate the mass-loss process of a close-orbit magnetized giant exoplanet, so-called hot Jupiter (HJ). In this paper we generalize the hydrodynamic (HD) model of an HJ expanding hydrogen atmosphere, proposed in the first paper, to include the effects of intrinsic planetary magnetic field. The proposed self-consistent axisymmetric 2D magnetohydrodynamics model incorporates radiative heating and ionization of the atmospheric gas, basic hydrogen chemistry for the appropriate account of major species composing HJ's upper atmosphere and related radiative energy deposition, and H3+ and Ly{\alpha} cooling processes. The model also takes into account a realistic solar-type X-ray/EUV spectrum for calculation of intensity and column density distribution of the radiative energy input, as well as gravitational and rotational forces acting in a tidally locked planet-star system. An interaction between the expanding atmospheric plasma and an intrinsic planetary magnetic dipole field leads to the formation of a current-carrying magnetodisk that plays an important role for topology and scaling of the planetary magnetosphere. A cyclic character of the magnetodisk behavior, composed of consequent phases of the disk formation followed by the magnetic reconnection with the ejection of a ring-type plasmoid, has been discovered and investigated. We found that the mass-loss rate of an HD 209458b analog planet is weakly affected by the equatorial surface field <0.3 G, but is suppressed by an order of magnitude at the field of 1 G.

We investigate the interaction between the magnetized stellar wind plasma and the partially ionized hydrodynamic hydrogen outflow from the escaping upper atmosphere of non- or weakly magnetized hot Jupiters. We use the well-studied hot Jupiter HD 209458b as an example for similar exoplanets, assuming a negligible intrinsic magnetic moment. For this planet, the stellar wind plasma interaction forms an obstacle in the planet's upper atmosphere, in which the position of the magnetopause is determined by the condition of pressure balance between the stellar wind and the expanded atmosphere, heated by the stellar extreme ultraviolet (EUV) radiation. We show that the neutral atmospheric atoms penetrate into the region dominated by the stellar wind, where they are ionized by photo-ionization and charge exchange, and then mixed with the stellar wind flow. Using a 3D magnetohydrodynamic (MHD) model, we show that an induced magnetic field forms in front of the planetary obstacle, which appears to be much stronger compared to those produced by the solar wind interaction with Venus and Mars. Depending on the stellar wind parameters, because of the induced magnetic field, the planetary obstacle can move up to ~0.5-1 planetary radii closer to the planet. Finally, we discuss how estimations of the intrinsic magnetic moment of hot Jupiters can be inferred by coupling hydrodynamic upper planetary atmosphere and MHD stellar wind interaction models together with UV observations. In particular, we find that HD 209458b should likely have an intrinsic magnetic moment of 10-20% that of Jupiter.

Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress.

We present a study of plasma conditions in the atmospheres of the Hot Jupiters HD 209458b and HD 189733b and for an HD 209458b-like planet at orbit locations between 0.2-1 AU around a Sun-like star. We discuss how these conditions influence the radio emission we expect from their magnetospheres. We find that the environmental conditions are such that the cyclotron maser instability (CMI), the process responsible for the generation of radio waves at magnetic planets in the solar system, most likely will not operate at Hot Jupiters. Hydrodynamically expanding atmospheres possess extended ionospheres whose plasma densities within the magnetosphere are so large that the plasma frequency is much higher than the cyclotron frequency, which contradicts the condition for the production of radio emission and prevents the escape of radio waves from close-in exoplanets at distances <0.05 AU from a Sun-like host star. The upper atmosphere structure of gas giants around stars similar to the Sun changes between 0.2 and 0.5 AU from the hydrodynamic to a hydrostatic regime and this results in conditions similar to solar system planets with a region of depleted plasma between the exobase and the magnetopause where the plasma frequency can be lower than the cyclotron frequency. In such an environment, a beam of highly energetic electrons accelerated along the field lines towards the planet can produce radio emission. However, even if the CMI could operate the extended ionospheres of Hot Jupiters are too dense to let the radio emission escape from the planets.

The interaction of escaping upper atmosphere of a hydrogen rich non-magnetized analog of HD209458b with a stellar wind of its host G-type star at different orbital distances is simulated with a 2D axisymmetric multi-fluid hydrodynamic model. A realistic sun-like spectrum of XUV radiation which ionizes and heats the planetary atmosphere, hydrogen photo-chemistry, as well as stellar-planetary tidal interaction are taken into account to generate self-consistently an atmospheric hydrodynamic outflow. Two different regimes of the planetary and stellar winds interaction have been modelled. These are: 1) the "captured by the star" regime, when the tidal force and pressure gradient drive the planetary material beyond the Roche lobe towards the star, and 2) the "blown by the wind" regime, when sufficiently strong stellar wind confines the escaping planetary atmosphere and channels it into the tail. The model simulates in details the hydrodynamic interaction between the planetary atoms, protons and the stellar wind, as well as the production of energetic neutral atoms (ENAs) around the planet due to charge-exchange between planetary atoms and stellar protons. The revealed location and shape of the ENA cloud either as a paraboloid shell between ionopause and bowshock (for the "blown by the wind" regime), or a turbulent layer at the contact boundary between the planetary stream and stellar wind (for the "captured by the star" regime) are of importance for the interpretation of Ly{\alpha} absorption features in exoplanetary transit spectra and characterization of the plasma environments.

In the present series of papers we propose a consistent description of the mass loss process. To study the effects of intrinsic magnetic field of a close-orbit giant exoplanet (so-called Hot Jupiter) on the atmospheric material escape and formation of planetary inner magnetosphere in a comprehensive way, we start with a hydrodynamic model of an upper atmosphere expansion presented in this paper. While considering a simple hydrogen atmosphere model, we focus on selfconsistent inclusion of the effects of radiative heating and ionization of the atmospheric gas with its consequent expansion in the outer space. Primary attention is paid to investigation of the role of specific conditions at the inner and outer boundaries of the simulation domain, under which different regimes of material escape (free- and restricted- flow) are formed. Comparative study of different processes, such as XUV heating, material ionization and recombination, H3+ cooling, adiabatic and Lyman-alpha cooling, Lyman-alpha reabsorption is performed. We confirm basic consistence of the outcomes of our modeling with the results of other hydrodynamic models of expanding planetary atmospheres. In particular, we obtain that under the typical conditions of an orbital distance 0.05 AU around a Sun-type star a Hot Jupiter plasma envelope may reach maximum temperatures up to ~9000K with a hydrodynamic escape speed ~9 km/s resulting in the mass loss rates ~(4-7)*10^10 g*s . In the range of considered stellar-planetary parameters and XUV fluxes that is close to mass loss in the energy limited case. The inclusion of planetary intrinsic magnetic fields in the model is a subject of the following up paper (Paper II).