Search

Ca-Al-rich inclusions (CAIs) are the oldest dated solid materials in the solar system, found as light-coloured crystalline ingredients in meteorites. Their formation time is commonly associated with age zero of the Solar System. Yet, the physical and chemical processes that once led to the formation of these sub-millimetre to centimetre-sized mineral particles in the early solar nebula are still a matter of debate. This paper proposes a pathway to form such inclusions during the earliest phases of disc evolution. We combine 1D viscous disc evolutionary models with 2D radiative transfer, equilibrium condensation, and new dust opacity calculations. We show that the viscous heating associated with the high accretion rates in the earliest evolutionary phases causes the midplane inside of about 0.5 au to heat up to limiting temperatures of about 1500-1700 K, but no further. These high temperatures force all refractory material components of the inherited interstellar dust grains to sublimate - except for a few Al-Ca-Ti oxides such as Al2O3, Ca2Al2SiO7, and CaTiO3. Once the Mg-Fe silicates are gone, the dust becomes more transparent and the heat is more efficiently transported to the disc surface, which prevents any further warm-up. This thermostat mechanism keeps these minerals above their annealing temperature for hundreds of thousands of years, which creates large, pure and crystalline particles. These particles are dragged out by the viscously spreading disc. Beyond about 0.5 au, the silicates re-condense on the Ca-Al-rich particles, adding an amorphous silicate matrix. We estimate that this mechanism to produce CAIs works during the first 50000 years of disc evolution. These particles then continue to move outward and populate the entire disc up to radii of about 50 au, before, eventually, the accretion rate subsides, the disc cools, and the particles start to drift inwards., Comment: accepted by Astronomy & Astrophysics, 17 pages, 10 figures, 3 tables

The detection of radio emission from an exoplanet would constitute the best way to determine its magnetic field. Indeed, the presence of a planetary magnetic field is a necessary condition for radio emission via the Cyclotron Maser Instability. The presence of a magnetic field is, however, not sufficient. At the emission site, the local cyclotron frequency has to be sufficiently high compared to the local plasma frequency. As strong stellar insolation on a low-mass planet can lead to an extended planetary atmosphere, the magnetospheric plasma frequency depends on the planetary mass, its orbital distance, and its host star. We show that an extended planetary atmosphere can quench the radio emission. This seems to be true, in particular, for an important fraction of the planets less massive than approximately two Jupiter masses and with orbital distances below $\sim$0.2 AU. Most of the best candidates suggested by radio scaling laws lie in this parameter space. Taking this effect quenching into account will have important implications for the target selection of observation campaigns. At the same time, this effect will have consequences for the interpretation of observational data., Comment: 6 pages, 1 figure. Proceedings of Planetary, Solar and Heliospheric Radio Emissions IX (2023)

The initial abundance of radioactive heat producing isotopes in the interior of a terrestrial planet are important drivers of its thermal evolution and the related tectonics and possible evolution to an Earth-like habitat. The moderately volatile element K can be outgassed from a magma ocean into H$_2$-dominated primordial atmospheres of protoplanets with assumed masses between 0.55-1.0$ M_{\rm Earth}$ at the time when the gas disk evaporated. We estimate this outgassing and let these planets grow through impacts of depleted and non-depleted material that resembles the same $^{40}$K abundance of average carbonaceous chondrites until the growing protoplanets reach 1.0 $M_{\rm Earth}$. We examine different atmospheric compositions and, as a function of pressure and temperature, calculate the proportion of K by Gibbs Free Energy minimisation using the GGChem code. We find that for H$_2$-envelopes and for magma ocean surface temperatures that are $\ge$ 2500 K, no K condensates are thermally stable, so that outgassed $^{40}$K can populate the atmosphere to a great extent. However, due to magma ocean turn-over time and the limited diffusion of $^{40}$K into the upper atmosphere, from the entire $^{40}$K in the magma ocean only a fraction may be available for escaping into space. The escape rates of the primordial atmospheres and the dragged $^{40}$K are further simulated for different stellar EUV-activities with a multispecies hydrodynamic upper atmosphere evolution model. Our results lead to different abundances of heat producing elements within the fully grown planets which may give rise to different thermal and tectonic histories of terrestrial planets and their habitability conditions., Comment: 22 pages, 11 figures. This is a preprint of a 2nd revision submitted to MNRAS

A 3D fully self-consistent multi-fluid hydrodynamic aeronomy model is applied to simulate the hydrogen-helium expanding upper atmosphere of the hot Jupiter HD189733b, and related absorption in the Lya line and the 10830 A line of metastable helium. We studied the influence of a high-energy stellar flux, stellar wind, and Lya cooling to reproduce the available observations. We found that to fit the width of the absorption profile in 10830 A line the escaping upper atmosphere of planet should be close to the energy limited escape achieved with a significantly reduced Lya cooling at the altitudes with HI density higher than 3*10^6 cm^-3. Based on the preformed simulations, we constrain the helium abundance in the upper atmosphere of HD189733b by a rather low value of He/H~0.005. We show that under conditions of a moderate stellar wind similar to that of the Sun the absorption of Lya line takes place mostly within the Roche lobe due to thermal broadening at a level of about 7%. At an order of magnitude stronger wind, a significant absorption of about 15% at high blue shifted velocities of up to 100 km/s is generated in the bowshock region, due to Doppler broadening. These blue shifted velocities are still lower than those (~200 km/s) detected in one of the observations. We explain the differences between performed observations, though not in all the details, by the stellar activity and the related fluctuations of the ionizing radiation (in case of 10830 A line), and stellar wind (in case of Lya line).

In this review we discuss all the relevant solar/stellar radiation and plasma parameters and processes that act together in the formation and modification of atmospheres and exospheres that consist of surface-related minerals. Magma ocean degassed silicate atmospheres or thin gaseous envelopes from planetary building blocks, airless bodies in the inner Solar System, and close-in magmatic rocky exoplanets such as CoRot-7b, HD219134b and 55 Cnc e are addressed. The depletion and fractionation of elements from planetary embryos, which act as the building blocks for protoplanets are also discussed. In this context the formation processes of the Moon and Mercury are briefly reviewed. The Lunar surface modification since its origin by micrometeoroids, plasma sputtering, plasma impingement as well as chemical surface alteration and the search of particles from the early Earth's atmosphere that were collected by the Moon on its surface are also discussed. Finally, we address important questions on what can be learned from the study of Mercury's environment and its solar wind interaction by MESSENGER and BepiColombo in comparison with the expected observations at exo-Mercurys by future space-observatories such as the JWST or ARIEL and ground-based telescopes and instruments like SPHERE and ESPRESSO on the VLT, and vice versa., Comment: 105 pages, 26 figures. This is a preprint of an article that is accepted for publication in Space Science Reviews

The signs of an expanding atmosphere of HD209458b have been observed with far-ultraviolet transmission spectroscopy and in the measurements of transit absorption by metastable HeI. These observations are interpreted using the hydrodynamic and Monte-Carlo numerical simulations of various degree of complexity and consistency. At the same time, no attempt has been made to model atmospheric escape of a magnetized HD209458b, to see how the planetary magnetic field might affect the measured transit absorption lines. This paper presents the global 3D MHD self-consistent simulations of the expanding upper atmosphere of HD209458b interacting with the stellar wind, and models the observed HI (Lya), OI (1306 A), CII (1337 A), and HeI (10830 A) transit absorption features. We find that the planetary dipole magnetic field with the equatorial surface value of Bp = 1 G profoundly changes the character of atmospheric material outflow and the related absorption. We also investigate the formation of planetary magnetosphere in the stellar wind and show that its size is more determined by the escaping atmosphere flow rather than by the strength of magnetic field. Fitting of the simulation results to observations enables constraining the stellar XUV flux and He abundance at ~10 erg cm2/s (at 1 a.u.) and He/H=0.02, respectively, as well as setting an upper limit for the dipole magnetic field of Bp<0.1 G on the planetary surface at the equator. This implies that the magnetic dipole moment of HD209458b should be less than 0.06 of the Jovian value.

Despite their importance for determining the evolution of the Earth's atmosphere and surface conditions, the evolutionary histories of the Earth's atmospheric CO$_2$ abundance during the Archean eon and the Sun's activity are poorly constrained. In this study, we apply a state-of-the-art physical model for the upper atmosphere of the Archean Earth to study the effects of different atmospheric CO$_2$/N$_2$ mixing ratios and solar activity levels on the escape of the atmosphere to space. We find that unless CO$_2$ was a major constituent of the atmosphere during the Archean eon, enhanced heating of the thermosphere by the Sun's strong X-ray and ultraviolet radiation would have caused rapid escape to space. We derive lower limits on the atmospheric CO$_2$ abundance of approximately 40\% at 3.8~billion years ago, which is likely enough to counteract the faint young Sun and keep the Earth from being completely frozen. Furthermore, our results indicate that the Sun was most likely born as a slow to moderate {rotating young G-star} to prevent rapid escape, putting essential constraints on the Sun's activity evolution throughout the solar system's history. In case that there were yet unknown cooling mechanisms present in the Archean atmosphere, this could reduce our CO$_2$ stability limits, and it would allow a more active Sun., Comment: 40 pages, 6 figures. This is a preprint accepted to be published in Earth and Planetary Science Letters

The ion escape of Mars' CO$_2$ atmosphere caused by its dissociation products C and O atoms is {simulated} from present time to $\sim 4.1$ billion years ago (Ga) by {numerical models of the upper atmosphere and its interaction with the solar wind}. The planetward-scattered pick-up ions are used for sputtering estimates of exospheric particles including $^{36}$Ar and $^{38}$Ar isotopes. Total ion escape, sputtering and photochemical escape rates are compared. For solar EUV fluxes $\geq$\,3 times that of today's Sun (earlier than $\sim 2.6$ Ga) ion escape becomes the dominant atmospheric non-thermal loss process until thermal escape takes over during the pre-Noachian eon (earlier than $\sim 4.0\,-\,4.1$ Ga). If we extrapolate the total escape of CO$_2$-related dissociation products back in time until $\sim$ 4.1 Ga we obtain a {maximum} theoretical equivalent to CO$_2$ partial pressure of more than {$\sim 0.4$ bar through non-thermal escape during quiet solar wind conditions}. {However, surface-atmosphere interaction and/or extreme solar events such as frequent CMEs could have increased this value even further. By including the surface as a sink, up to 0.9\,bar, or even up to 1.8\,bar in case of hidden carbonate reservoirs, could have been present at 4.1\,Ga} The fractionation of $^{36}$Ar/$^{38}$Ar isotopes through sputtering and volcanic outgassing from its initial chondritic value of 5.3, as measured in the 4.1 billion years old Mars meteorite ALH 84001, until the present day, however, can be reproduced for assumed CO$_2$ partial pressures of {$\sim0.01 -- 0.4$\,bar without, and $\sim0.4 -- 1.8$\,bar including surface sinks, and} depending on the cessation time of the Martian dynamo (assumed between 3.6\,--\,4.0 Ga) - if atmospheric sputtering of Ar started afterwards., Comment: 63 pages (including 8 pages of supplement), 21 figures and 20 tables (including 12 figures and 15 tables in the supplement), resubmitted to Icarus

We discuss the current state of knowledge of terrestrial planet formation from the aspects of different planet formation models and isotopic data from 182Hf-182W, U-Pb, lithophile-siderophile elements, 48Ca/44Ca isotope samples from planetary building blocks, 36Ar/38Ar, 20Ne/22Ne, 36Ar/22Ne isotope ratios in Venus' and Earth's atmospheres, the expected solar 3He abundance in Earth's deep mantle and Earth's D/H sea water ratios that shed light on the accretion time of the early protoplanets. Accretion scenarios that can explain the different isotope ratios, including a Moon-forming event after ca. 50 Myr, support the theory that the bulk of Earth's mass (>80%) most likely accreted within 10-30 Myr. From a combined analysis of the before mentioned isotopes, one finds that proto-Earth accreted 0.5-0.6 MEarth within the first ~4-5 Myr, the approximate lifetime of the protoplanetary disk. For Venus, the available atmospheric noble gas data are too uncertain for constraining the planet's accretion scenario accurately. However, from the available Ar and Ne isotope measurements, one finds that proto-Venus could have grown to 0.85-1.0 MVenus before the disk dissipated. Classical terrestrial planet formation models have struggled to grow large planetary embryos quickly from the tiniest materials within the typical lifetime of protoplanetary disks. Pebble accretion could solve this long-standing time scale controversy. Pebble accretion and streaming instabilities produce large planetesimals that grow into Mars-sized and larger planetary embryos during this early accretion phase. The later stage of accretion can be explained well with the Grand-Tack, annulus or depleted disk models. The relative roles of pebble accretion and planetesimal accretion/giant impacts are poorly understood and should be investigated with N-body simulations that include pebbles and multiple protoplanets., Comment: 49 pages, 7 figures. This is a preprint of an article published in Space Science Reviews. The final authenticated version can be found online at: https://doi.org/10.1007/s11214-020-00778-4

It is not yet entirely clear whether Mars began as a warm and wet planet that evolved towards the present-day cold and dry body or if it always was cold and dry with just some sporadic episodes of liquid water on its surface. An important clue into this question can be gained by studying the earliest evolution of the Martian atmosphere and whether it was dense and stable to maintain a warm and wet climate or tenuous and susceptible to strong atmospheric escape. We discuss relevant aspects for the evolution and stability of a potential early Martian atmosphere. This contains the solar EUV flux evolution, the formation timescale and volatile inventory of the planet including volcanic degassing, impact delivery and removal, the loss of a catastrophically outgassed steam atmosphere, atmosphere-surface interactions, and thermal and non-thermal escape processes affecting any secondary atmosphere. While early non-thermal escape at Mars before 4 billion years ago (Ga) is poorly understood, particularly in view of its ancient intrinsic magnetic field, research on thermal escape processes indicate that volatile delivery and volcanic degassing cannot counterbalance the strong thermal escape. Therefore, a catastrophically outgassed steam atmosphere of several bars of CO2 and H2O, or CO and H2 for reduced conditions, could have been lost within just a few million years (Myr). Thereafter, Mars likely could not build up a dense secondary atmosphere during its first ~400 Myr but might only have possessed an atmosphere sporadically during events of strong volcanic degassing, potentially also including SO2. This indicates that before ~4.1 Ga Mars indeed might have been cold and dry. A denser CO2- or CO-dominated atmosphere, however, might have built up afterwards but must have been lost later-on due to non-thermal escape processes and sequestration into the ground., Comment: 120 pages, 12 figures. This is a preprint of an article published in Space Science Reviews. The final authenticated version can be found online at https://doi.org/10.1007/s11214-020-00779-3

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., Comment: to be published

Here we discuss the current state of knowledge on how atmospheric escape processes can fractionate noble gas isotopes and moderately volatile rock-forming elements that populate primordial atmospheres, magma ocean related environments, and catastrophically outgassed steam atmospheres. Variations of isotopes and volatile elements in different planetary reservoirs keep information about atmospheric escape, composition and even the source of accreting material. We summarize our knowledge on atmospheric isotope ratios and discuss the latest evidence that proto-Venus and Earth captured small H$_2$-dominated primordial atmospheres that were lost by EUV-driven hydrodynamic escape after the disk dispersed. All relevant thermal and non-thermal atmospheric escape processes that can fractionate various isotopes and volatile elements are discussed. Erosion of early atmospheres, crust and mantle by large planetary impactors are also addressed. Further, we discuss how moderately volatile elements such as the radioactive heat producing element $^{40}$K and other rock-forming elements such as Mg can also be outgassed and lost from magma oceans that originate on large planetary embryos and accreting planets. Outgassed elements escape from planetary embryos with masses that are $\geq$\,M$_{\rm Moon}$ directly, or due to hydrodynamic drag of escaping H atoms originating from primordial- or steam atmospheres at more massive embryos. We discuss how these processes affect the final elemental composition and ratios such as K/U, Fe/Mg of early planets and their building blocks. Finally, we review modeling efforts that constrain the early evolution of Venus, Earth and Mars by reproducing their measured present day atmospheric $^{36}$Ar/$^{38}$Ar, $^{20}$Ne/$^{22}$Ne noble gas isotope ratios and the role of isotopes on the loss of water and its connection to the redox state on early Mars., Comment: 65 pages, 17 figures. This is a preprint of an article published in Space Science Reviews. The final authenticated version can be found online at: https://doi.org/10.1007/s11214-020-00701-x

This brief review will discuss the current knowledge on the origin and evolution of the nitrogen atmospheres of the icy bodies in the solar system, particularly of Titan, Triton and Pluto. An important tool to analyse and understand the origin and evolution of these atmospheres can be found in the different isotopic signatures of their atmospheric constituents. The $^{14}$N/$^{15}$N ratio of the N$_2$-dominated atmospheres of these bodies serve as a footprint of the building blocks from which Titan, Triton and Pluto originated and of the diverse fractionation processes that shaped these atmospheres over their entire evolution. Together with other measured isotopic and elemental ratios such as $^{12}$C/$^{13}$C or Ar/N these atmospheres can give important insights into the history of the icy bodies in the solar system, the diverse processes that affect their N$_2$-dominated atmospheres, and the therewith connected solar activity evolution. Titan's gaseous envelope most likely originated from ammonia ices with possible contributions from refractory organics. Its isotopic signatures can yet be seen in the - compared to Earth - comparatively heavy $^{14}$N/$^{15}$N ratio of 167.7, even though this value slightly evolved over its history due to atmospheric escape and photodissociation of N$_2$. The origin and evolution of Pluto's and Triton's tenuous nitrogen atmospheres remain unclear, even though it might be likely that their atmospheres originated from the protosolar nebula or from comets. An in-situ space mission to Triton such as the recently proposed Trident mission, and/or to the ice giants would be a crucial cornerstone for a better understanding of the origin and evolution of the icy bodies in the outer solar system and their atmospheres in general., Comment: 60 pages, 6 figures. This is a preprint of an article published in Space Science Reviews. The final authenticated version can be found online at : https://doi.org/10.1007/s11214-020-00752-0

The discovery of thousands of highly irradiated, low-mass, exoplanets has led to the idea that atmospheric escape is an important process that can drive their evolution. Of particular interest is the inference from recent exoplanet detections that there is a large population of low mass planets possessing significant, hydrogen dominated atmospheres, even at masses as low as $\sim 2$~M$_\oplus$. The size of these hydrogen dominated atmospheres indicates the the envelopes must have been accreted from the natal protoplanetary disc. This inference is in contradiction with the Solar System terrestrial planets, that did not reach their final masses before disc dispersal, and only accreted thin hydrogen dominated atmospheres. In this review, we discuss the evidence for hydrogen dominated atmospheres on terrestrial mass ($\lesssim$ 2~M$_\oplus$) planets. We then discuss the possible origins and evolution of these atmospheres with a focus on the role played by hydrodynamic atmospheric escape driven by the stellar high-energy emission (X-ray and EUV; XUV)., Comment: Accepted for publication in Space Science Reviews. Part of ISSI special collection on 'Understanding the Diversity of Planetary Atmospheres'

We apply a 1D upper atmosphere model to study thermal escape of nitrogen over Titan's history. Significant thermal escape should have occurred very early for solar EUV fluxes 100 to 400 times higher than today with escape rates as high as $\approx 1.5\times 10^{28}$ s$^{-1}$ and $\approx 4.5\times 10^{29}$ s$^{-1}$, respectively, while today it is $\approx 7.5\times 10^{17}$ s$^{-1}$. Depending on whether the Sun originated as a slow, moderate or fast rotator, thermal escape was the dominant escape process for the first 100 to 1000 Myr after the formation of the solar system. If Titan's atmosphere originated that early, it could have lost between $\approx 0.5 - 16$ times its present atmospheric mass depending on the Sun's rotational evolution. We also investigated the mass-balance parameter space for an outgassing of Titan's nitrogen through decomposition of NH$_3$-ices in its deep interior. Our study indicates that, if Titan's atmosphere originated at the beginning, it could have only survived until today if the Sun was a slow rotator. In other cases, the escape would have been too strong for the degassed nitrogen to survive until present-day, implying later outgassing or an additional nitrogen source. An endogenic origin of Titan's nitrogen partially through NH$_3$-ices is consistent with its initial fractionation of $^{14}$N/$^{15}$N $\approx$ 166 - 172, or lower if photochemical removal was relevant for longer than the last $\approx$ 1,000 Myr. Since this ratio is slightly above the ratio of cometary ammonia, some of Titan's nitrogen might have originated from refractory organics., Comment: 16 pages, 6 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS) following peer review. The version of record [https://doi.org/10.1093/mnras/staa3151] is available online at: https://academic.oup.com/mnras/advance-article-abstract/doi/10.1093/mnras/staa3151/5930166

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.

The development of habitable conditions on Earth is tightly connected to the evolution of its atmosphere which is strongly influenced by atmospheric escape. We investigate the evolution of the polar ion outflow from the open field line bundle which is the dominant escape mechanism for the modern Earth. We perform Direct Simulation Monte Carlo (DSMC) simulations and estimate the upper limits on escape rates from the Earth's open field line bundle starting from three gigayears ago (Ga) to present assuming the present-day composition of the atmosphere. We perform two additional simulations with lower mixing ratios of oxygen of 1% and 15% to account for the conditions shortly after the Great Oxydation Event (GOE). We estimate the maximum loss rates due to polar outflow three gigayears ago of $3.3 \times10^{27}$ s$^{-1}$ and $2.4 \times 10^{27}$ s$^{-1}$ for oxygen and nitrogen, respectively. The total integrated mass loss equals to 39% and 10% of the modern atmosphere's mass, for oxygen and nitrogen, respectively. According to our results, the main factors that governed the polar outflow in the considered time period are the evolution of the XUV radiation of the Sun and the atmosphere's composition. The evolution of the Earth's magnetic field plays a less important role. We conclude that although the atmosphere with the present-day composition can survive the escape due to polar outflow, a higher level of CO$_2$ between 3.0 and 2.0~Ga is likely necessary to reduce the escape., Comment: 14 pages, 6 figures, published by JGR Space Physics

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., Comment: Accepted for publication in A&A. The abstract has been shortened to fit the arXiv form

In the present chapter we present the results of evolutionary studies of exoplanetary atmospheres. We mostly focus on the sub- to super-Earth domain, although these methods are applicable to all types of exoplanets. We consider both thermal and nonthermal loss processes. The type of thermal loss mechanism depends on so-called escape parameter $\beta$, which is the ratio of the gravitational energy of a particle to its thermal energy. While $\beta$ is decreasing, an exoplanet switches from classical Jeans to modified Jeans and finally to blow-off escape mechanisms. During blow-off the majority of the atmospheric particles dispose of enough energy to escape the planet's gravity field. This leads to extreme gas losses. Although nonthermal losses never exceed blow-off escape, they are of significant importance for planets with relatively weak Jeans-type escape. From the diversity of nonthermal escape mechanisms, in the present chapter we focus on ion pickup and discuss the importance of other loss mechanisms. The general conclusion of the chapter is, that escape processes strongly shape the evolution of the exoplanets and determine, if the planet loses its atmosphere due to erosion processes or, on the contrary, stays as mini-Neptune type body, which can probably not be considered as a potential habitat as we know it., Comment: Originally published by: Springer International Publishing Switzerland 2015 H. Lammer, M. Khodachenko (eds.), Characterizing Stellar and Exoplanetary Environments, Astrophysics and Space Science Library 411

Coronal Mass Ejections (CMEs) may have major importance for planetary and stellar evolution. Stellar CME parameters, such as mass and velocity, have yet not been determined statistically. So far only a handful of stellar CMEs has been detected mainly on dMe stars using spectroscopic observations. We therefore aim for a statistical determination of CMEs of solar-like stars by using spectroscopic data from the ESO phase 3 and Polarbase archives. To identify stellar CMEs we use the Doppler signal in optical spectral lines being a signature of erupting filaments which are closely correlated to CMEs. We investigate more than 3700 hours of on-source time of in total 425 dF-dK stars. We find no signatures of CMEs and only few flares. To explain this low level of activity we derive upper limits for the non detections of CMEs and compare those with empirically modelled CME rates. To explain the low number of detected flares we adapt a flare power law derived from EUV data to the H{\alpha} regime, yielding more realistic results for H{\alpha} observations. In addition we examine the detectability of flares from the stars by extracting Sun-as-a-star H{\alpha} light curves. The extrapolated maximum numbers of observable CMEs are below the observationally determined upper limits, which indicates that the on-source times were mostly too short to detect stellar CMEs in H{\alpha}. We conclude that these non detections are related to observational biases in conjunction with a low level of activity of the investigated dF-dK stars., Comment: MNRAS accepted

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., Comment: Monthly Notices of the Royal Astronomical Society (2019)

With a Monte-Carlo model we investigate the escape of hot oxygen and carbon from the martian atmosphere for four points in time in its history corresponding to 1, 3, 10, and 20 times the present solar EUV flux. We study and discuss different sources of hot oxygen and carbon atoms in the thermosphere and their changing importance with the EUV flux. The increase of the production rates due to higher densities resulting from the higher EUV flux competes against the expansion of the thermosphere and corresponding increase in collisions. We find that the escape due to photodissociation increases with increasing EUV level. However, for the escape via some other reactions, e.g.~dissociative recombination of O$_2^+$, this is only true until the EUV level reaches 10 times the present EUV flux, and then the rates start to decrease. Furthermore, our results show that Mars could not have had a dense atmosphere at the end of the Noachian epoch, since such an atmosphere would not have been able to escape until today. In the pre-Noachian era, most of a magma ocean and volcanic activity related outgassed CO$_2$ atmosphere could have been lost thermally until the Noachian epoch, when non-thermal loss processes such as suprathermal atom escape became dominant. Thus, early Mars could have been hot and wet during the pre-Noachian era with surface CO$_2$ pressures larger than 1 bar during the first 300 Myr after the planet's origin.

The influence of the hydrogen hydrodynamic upper atmosphere escape, driven by the solar soft X-ray and extreme ultraviolet radiation (XUV) flux, on an expected magma ocean outgassed steam atmosphere of early Venus is studied. By assuming that the young Sun was either a weak or moderate active young G star, we estimated the water loss from a hydrogen dominated thermosphere due to the absorption of the solar XUV flux and the precipitation of solar wind produced energetic hydrogen atoms (ENAs). The production of ENAs and their interaction with the hydrodynamic extended upper atmosphere, including collision-related feedback processes, have been calculated by means of Monte Carlo models. ENAs that collide in the upper atmosphere deposit their energy and heat the surrounding gas mainly above the main XUV energy deposition layer. It is shown that precipitating ENAs modify the thermal structure of the upper atmosphere, but the enhancement of the thermal escape rates caused by these energetic hydrogen atoms is negligible. Our results also indicate that the majority of oxygen arising from dissociated H$_2$O molecules is left behind during the first 100 Myr. It is thus suggested that the main part of the remaining oxygen has been absorbed by crustal oxidation., Comment: 6 figures and 2 tables

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., Comment: 34 pages and 18 figures

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters., Comment: 8 pages, 3 figures

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., Comment: 23 pages, 8 figures

Aims. In this letter, we calculate for the first time the full transonic hydrodynamic escape of an Earth-like atmosphere. We consider the case of an Earth-mass planet with an atmospheric composition identical to that of the current Earth orbiting at 1 AU around a young and very active solar mass star. Methods. To model the upper atmosphere, we used the Kompot Code, which is a first-principles model that calculates the physical structures of the upper atmospheres of planets, taking into account hydrodynamics and the main chemical and thermal processes taking place in the upper atmosphere of a planet. This model enabled us to calculate the 1D vertical structure of the atmosphere using as input the high-energy spectrum of a young and active Sun. Results. The atmosphere has the form of a transonic hydrodynamic Parker wind, which has an outflow velocity at the upper boundary of our computational domain that exceeds the escape velocity. The outflowing gas is dominated by atomic nitrogen and oxygen and their ion equivalents and has a maximum ionization fraction of 20%. The mass outflow rate is found to be 1.8x10^9 g s^-1, which would erode the modern Earth's atmosphere in less than 0.1 Myr. Conclusions. This extreme mass loss rate suggests that an Earth-like atmosphere cannot form when the planet is orbiting within the habitable zone of a very active star. Instead, such an atmosphere can only form after the activity of the star has decreased to a much lower level. This happened in the early atmosphere of the Earth, which was likely dominated by other gases such as CO2. Since the time it takes for the activity of a star to decay is highly dependent on its mass, this is important for understanding possible formation timescales for planets orbiting low-mass stars., Comment: Accepted for publication by A&A Letters

Hydrogen-dominated atmospheres of hot exoplanets expand and escape due to the intense heating by the X-ray and extreme ultraviolet (XUV) irradiation of their host stars. Excess absorption of neutral hydrogen has been observed in the Ly$\alpha$ line during transits of several close-in exoplanets, indicating such extended atmospheres. For the hot Jupiter HD 189733b, this absorption shows temporal variability. Variations in stellar XUV emission and/or stellar wind conditions have been invoked to explain this effect. We apply a 1D hydrodynamic upper atmosphere model and a 3D MHD stellar wind flow model to study the effect of variations of the stellar XUV and wind conditions on the neutral hydrogen distribution, including the production of energetic neutral atoms (ENAs), and the related Ly$\alpha$ transit signature. We obtain comparable, albeit slightly higher Ly$\alpha$ absorption as observed in 2011 with a stellar XUV flux of $1.8\times10^4$ erg cm$^{-2}$ s$^{-1}$, rather typical activity conditions for this star. Flares similar to the one observed 8 h before the transit are unlikely to have caused a significant modulation of the transit signature. The resulting Ly$\alpha$ absorption is dominated by atmospheric broadening, whereas the contribution of ENAs is negligible, as they are formed inside the bow shock from decelerated wind ions that are heated to high temperatures. Thus, within our modeling framework and assumptions, we find an insignificant dependence on the stellar wind parameters. Since the transit absorption can be modeled with typical stellar XUV and wind conditions, it is possible that the non-detection of the absorption in 2010 was affected by less typical stellar activity conditions, such as a very different magnitude and/or shape of the star's spectral XUV emission, or temporal/spatial variations in Ly$\alpha$ affecting the determination of the transit absorption., Comment: 22 pages, 19 figures, 4 tables; A&A, published

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., Comment: 17 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

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.

We investigate under which conditions supermassive hot Jupiters can sustain source regions for radio emission, and whether this emission could propagate to an observer outside the system. We study Tau Bootis b-like planets (a supermassive hot Jupiter with 5.84 Jupiter masses and 1.06 Jupiter radii), but located at different orbital distances (between its actual orbit of 0.046 AU and 0.2 AU). Due to the strong gravity of such planets and efficient radiative cooling, the upper atmosphere is (almost) hydrostatic and the exobase remains very close to the planet, which makes it a good candidate for radio observations. We expect similar conditions as for Jupiter, i.e. a region between the exobase and the magnetopause that is filled with a depleted plasma density compared with cases where the whole magnetosphere cavity is filled with hydrodynamically outward flowing ionospheric plasma. Thus, unlike classical hot Jupiters like the previously studied planets HD 209458b and HD 189733b, supermassive hot Jupiters should be in general better targets for radio observations.

Planets that are embedded in the changing magnetic fields of their host stars can experience significant induction heating in their interiors caused by the planet's orbital motion. For induction heating to be substantial, the planetary orbit has to be inclined with respect to the stellar rotation and dipole axes. Using WX~UMa, for which the rotation and magnetic axes are aligned, as an example, we show that for close-in planets on inclined orbits, induction heating can be stronger than the tidal heating occurring inside Jupiter's satellite Io; namely, it can generate a surface heat flux exceeding 2\,W\,m$^{-2}$. An internal heating source of such magnitude can lead to extreme volcanic activity on the planet's surface, possibly also to internal local magma oceans, and to the formation of a plasma torus around the star aligned with the planetary orbit. A strongly volcanically active planet would eject into space mostly SO$_2$, which would then dissociate into oxygen and sulphur atoms. Young planets would also eject CO$_2$. Oxygen would therefore be the major component of the torus. If the O{\sc i} column density of the torus exceeds $\approx$10$^{12}$\,cm$^{-2}$, the torus could be revealed by detecting absorption signatures at the position of the strong far-ultraviolet O{\sc i} triplet at about 1304\,\AA. We estimate that this condition is satisfied if the O{\sc i} atoms in the torus escape the system at a velocity smaller than 1--10\,km\,s$^{-1}$. These estimates are valid also for a tidally heated planet., Comment: 8 pages, 6 figures, accepted for publication in ApJ

The K2-33 planetary system hosts one transiting ~5 R_E planet orbiting the young M-type host star. The planet's mass is still unknown, with an estimated upper limit of 5.4 M_J. The extreme youth of the system (<20 Myr) gives the unprecedented opportunity to study the earliest phases of planetary evolution, at a stage when the planet is exposed to an extremely high level of high-energy radiation emitted by the host star. We perform a series of 1D hydrodynamic simulations of the planet's upper atmosphere considering a range of possible planetary masses, from 2 to 40 M_E, and equilibrium temperatures, from 850 to 1300 K, to account for internal heating as a result of contraction. We obtain temperature profiles mostly controlled by the planet's mass, while the equilibrium temperature has a secondary effect. For planetary masses below 7-10 M_E, the atmosphere is subject to extremely high escape rates, driven by the planet's weak gravity and high thermal energy, which increase with decreasing mass and/or increasing temperature. For higher masses, the escape is instead driven by the absorption of the high-energy stellar radiation. A rough comparison of the timescales for complete atmospheric escape and age of the system indicates that the planet is more massive than 10 M_E., Comment: 11 pages, 7 figures

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., Comment: open access

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 present a study of an impacting descent probe that increases the science return of spacecraft orbiting or passing an atmosphere-less planetary body of the solar system, such as the Galilean moons of Jupiter. The descent probe is a carry-on small spacecraft (< 100 kg), to be deployed by the mother spacecraft, that brings itself onto a collisional trajectory with the targeted planetary body in a simple manner. A possible science payload includes instruments for surface imaging, characterisation of the neutral exosphere, and magnetic field and plasma measurement near the target body down to very low-altitudes (~1 km), during the probe's fast (~km/s) descent to the surface until impact. The science goals and the concept of operation are discussed with particular reference to Europa, including options for flying through water plumes and after-impact retrieval of very-low altitude science data. All in all, it is demonstrated how the descent probe has the potential to provide a high science return to a mission at a low extra level of complexity, engineering effort, and risk. This study builds upon earlier studies for a Callisto Descent Probe (CDP) for the former Europa-Jupiter System Mission (EJSM) of ESA and NASA, and extends them with a detailed assessment of a descent probe designed to be an additional science payload for the NASA Europa Mission., Comment: 34 pages, 11 figures

Several studies have shown that stellar activity features, such as occulted and non-occulted starspots, can affect the measurement of transit parameters biasing studies of transit timing variations and transmission spectra. We present $\texttt{PyTranSpot}$, which we designed to model multiband transit light curves showing starspot anomalies, inferring both transit and spot parameters. The code follows a pixellation approach to model the star with its corresponding limb darkening, spots, and transiting planet on a two dimensional Cartesian coordinate grid. We combine $\texttt{PyTranSpot}$ with an MCMC framework to study and derive exoplanet transmission spectra, which provides statistically robust values for the physical properties and uncertainties of a transiting star-planet system. We validate $\texttt{PyTranSpot}$'s performance by analyzing eleven synthetic light curves of four different star-planet systems and 20 transit light curves of the well-studied WASP-41b system. We also investigate the impact of starspots on transit parameters and derive wavelength dependent transit depth values for WASP-41b covering a range of 6200-9200 $\AA$, indicating a flat transmission spectrum., Comment: 17 pages, 22 figures; accepted for publication in Astronomy & Astrophysics

Low-mass M stars are plentiful in the Universe and often host small, rocky planets detectable with the current instrumentation. Recently, seven small planets have been discovered orbiting the ultracool dwarf TRAPPIST-1\cite{Gillon16,Gillon17}. We examine the role of electromagnetic induction heating of these planets, caused by the star's rotation and the planet's orbital motion. If the stellar rotation and magnetic dipole axes are inclined with respect to each other, induction heating can melt the upper mantle and enormously increase volcanic activity, sometimes producing a magma ocean below the planetary surface. We show that induction heating leads the three innermost planets, one of which is in the habitable zone, to either evolve towards a molten mantle planet, or to experience increased outgassing and volcanic activity, while the four outermost planets remain mostly unaffected., Comment: Published in Nature Astronomy; https://www.nature.com/articles/s41550-017-0284-0

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., Comment: 8 pages, 6 figures, 2 tables, accepted to MNRAS

Stellar coronal mass ejections (CMEs) may play an important role in mass- and angular momentum loss of young Sun-like stars. If occurring frequently, they may also have a strong effect on planetary evolution by increasing atmospheric erosion. So far it has not been possible to infer the occurrence frequency of stellar CMEs from observations. Based on their close relation with flares on the Sun, we develop an empirical model combining solar flare-CME relationships with stellar flare rates to estimate the CME activity of young Sun-like and late-type main-sequence stars. By comparison of the obtained CME mass-loss rates with observations of total mass-loss rates, we find that our modeled rates may exceed those from observations by orders of magnitude for the most active stars. This reveals a possible limit to the extrapolation of such models to the youngest stars. We find that the most uncertain component in the model is the flare-CME association rate adopted from the Sun, which does not properly account for the likely stronger coronal confinement in active stars. Simple estimates of this effect reveal a possible suppression of CME rates by several orders of magnitude for young stars, indicating that this issue should be addressed in more detail in the future., Comment: 18 pages, 8 figures, 1 table, accepted to MNRAS

Planetary embryos form protoplanets via mutual collisions, which can lead to the development of magma oceans. During their solidification, large amounts of the mantles' volatile contents may be outgassed. The resulting H$_2$O/CO$_2$ dominated steam atmospheres may be lost efficiently via hydrodynamic escape due to the low gravity and the high stellar EUV luminosities. Protoplanets forming later from such degassed building blocks could therefore be drier than previously expected. We model the outgassing and subsequent hydrodynamic escape of steam atmospheres from such embryos. The efficient outflow of H drags along heavier species (O, CO$_2$, noble gases). The full range of possible EUV evolution tracks of a solar-mass star is taken into account to investigate the escape from Mars-sized embryos at different orbital distances. The envelopes are typically lost within a few to a few tens of Myr. Furthermore, we study the influence on protoplanetary evolution, exemplified by Venus. We investigate different early evolution scenarios and constrain realistic cases by comparing modeled noble gas isotope ratios with observations. Starting from solar values, consistent isotope ratios (Ne, Ar) can be found for different solar EUV histories, as well as assumptions about the initial atmosphere (either pure steam or a mixture with accreted H). Our results generally favor an early accretion scenario with a small amount of accreted H and a low-activity Sun, because in other cases too much CO$_2$ is lost during evolution, which is inconsistent with Venus' present atmosphere. Important issues are likely the time at which the initial steam atmosphere is outgassed and/or the amount of CO$_2$ which may still be delivered at later evolutionary stages. A late accretion scenario can only reproduce present isotope ratios for a highly active young Sun, but then very massive steam atmospheres would be required., Comment: 61 pages, 7 figures, 3 tables, accepted to Icarus

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., Comment: Published on ApJ

Stimulated by the discovery of a number of close-in low-density planets, we generalise the Jeans escape parameter taking hydrodynamic and Roche lobe effects into account. We furthermore define $\Lambda$ as the value of the Jeans escape parameter calculated at the observed planetary radius and mass for the planet's equilibrium temperature and considering atomic hydrogen, independently of the atmospheric temperature profile. We consider 5 and 10 $M_{\oplus}$ planets with an equilibrium temperature of 500 and 1000 K, orbiting early G-, K-, and M-type stars. Assuming a clear atmosphere and by comparing escape rates obtained from the energy-limited formula, which only accounts for the heating induced by the absorption of the high-energy stellar radiation, and from a hydrodynamic atmosphere code, which also accounts for the bolometric heating, we find that planets whose $\Lambda$ is smaller than 15-35 lie in the "boil-off" regime, where the escape is driven by the atmospheric thermal energy and low planetary gravity. We find that the atmosphere of hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) low-mass ($M_{\rm pl}\lessapprox$ 5 $M_{\oplus}$) planets with $\Lambda$ < 15-35 shrinks to smaller radii so that their $\Lambda$ evolves to values higher than 15-35, hence out of the boil-off regime, in less than $\approx$500 Myr. Because of their small Roche lobe radius, we find the same result also for hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) higher mass ($M_{\rm pl}\lessapprox$ 10 $M_{\oplus}$) planets with $\Lambda$ < 15-35, when they orbit M-dwarfs. For old, hydrogen-dominated planets in this range of parameters, $\Lambda$ should therefore be $\geq$15-35, which provides a strong constraint on the planetary minimum mass and maximum radius and can be used to predict the presence of aerosols and/or constrain planetary masses, for example., Comment: Accepted for publication by A&A

We present the analysis of six nights of spectroscopic monitoring of two young and fast rotating late-type stars, namely the dMe star HK Aqr and the dG/dK star PZ Tel. On both stars we detect absorption features reminiscent of signatures of co-rotating cool clouds or prominences visible in H$\alpha$. Several prominences on HK Aqr show periodic variability in the prominence tracks which follow a sinusoidal motion (indication of prominence oscillations). On PZ Tel we could not find any periodic variability in the prominence tracks. By fitting sinusoidal functions to the prominence tracks we derive amplitudes and periods which are similar to those of large amplitude oscillations seen in solar prominences. In one specific event we also derive a periodic variation of the prominence track in the H$\beta$ spectral line which shows an anti-phase variation with the one derived for the H$\alpha$ spectral line. Using these parameters and estimated mass density of a prominence on HK Aqr we derive a minimum magnetic field strength of $\sim$2G. The relatively low strength of the magnetic field is explained by the large height of this stellar prominence ($\ge$ 0.67 stellar radii above the surface).

For the hot exoplanets CoRoT-24b and CoRoT-24c, observations have provided transit radii R$_{\rm T}$ of 3.7$\pm$0.4 R$_{\oplus}$ and 4.9$\pm$0.5 R$_{\oplus}$, and masses of $\le$5.7 M$_{\oplus}$ and 28$\pm$11 M$_{\oplus}$, respectively. We study their upper atmosphere structure and escape applying an hydrodynamic model. Assuming R$_{\rm T} \approx$ R$_{\rm PL}$, where R$_{\rm PL}$ is the planetary radius at the pressure of 100 mbar, we obtained for CoRoT-24b unrealistically high thermally-driven hydrodynamic escape rates. This is due to the planet's high temperature and low gravity, independent of the stellar EUV flux. Such high escape rates could last only for $<$100 Myr, while R$_{\rm PL}$ shrinks till the escape rate becomes less than or equal to the maximum possible EUV-driven escape rate. For CoRoT-24b, R$_{\rm PL}$ must be therefore located at $\approx 1.9-2.2$ R$_{\oplus}$ and high altitude hazes/clouds possibly extinct the light at R$_{\rm T}$. Our analysis constraints also the planet's mass to be 5$-$5.7 M$_{\oplus}$. For CoRoT-24c, R$_{\rm PL}$ and R$_{\rm T}$ lie too close together to be distinguished in the same way. Similar differences between R$_{\rm PL}$ and R$_{\rm T}$ may be present also for other hot, low-density sub-Neptunes., Comment: Accepted for publication by MNRAS

We suggest an additional possible plasma source to explain part of the phenomena observed for the transiting hot Jupiters WASP-12b and HD 189733b in their ultraviolet (UV) light curves. In the proposed scenario, material outgasses from the molten surface of Trojan satellites on tadpole orbits near the Lagrange points L$_4$ and L$_5$. We show that the temperature at the orbital location of WASP-12b is high enough to melt the surface of rocky bodies and to form shallow lava oceans on them. In case of WASP-12b, this leads to the release of elements such as Mg and Ca, which are expected to surround the system. The predicted Mg and Ca outgassing rates from two Io-sized WASP-12b Trojans are $\approx 2.2 \times 10^{27}$ s$^{-1}$ and $\approx 2.2 \times 10^{26}$ s$^{-1}$, respectively. Trojan outgassing can lead to the apparent lack of emission in Mg{\sc ii}\,h\&k and Ca{\sc ii}\,H\&K line cores of WASP-12. For HD 189733b, the mechanism is only marginally possible due to the lower temperature. This may be one of the reasons that couldn't explain the early ingress of HD 189733b observed in the far-UV (FUV) C{\sc ii} doublet due to absence of carbon within elements outgassed by molten lava. We investigate the long-term stability region of WASP-12b and HD 189733b in case of planar and inclined motion of these satellites and show that unlike the classical exomoons orbiting the planet, Io-sized Trojans can be stable for the whole systems life time., Comment: 14 pages, 7 figures, accepted to MNRAS

We investigate the loss rates of the hydrogen atmospheres of terrestrial planets with a range of masses and orbital distances by assuming a stellar extreme ultraviolet (EUV) luminosity that is 100 times stronger than that of the current Sun. We apply a 1D upper atmosphere radiation absorption and hydrodynamic escape model that takes into account ionization, dissociation and recombination to calculate hydrogen mass loss rates. We study the effects of the ionization, dissociation and recombination on the thermal mass loss rates of hydrogen-dominated super-Earths and compare the results to those obtained by the energy-limited escape formula which is widely used for mass loss evolution studies. Our results indicate that the energy-limited formula can to a great extent over- or underestimate the hydrogen mass loss rates by amounts that depend on the stellar EUV flux and planetary parameters such as mass, size, effective temperature, and EUV absorption radius., Comment: 8 pages, 2 figures, 1 table, submitted to MNRAS