We present an implicit numerical method to solve the time-dependent equations of radiation hydrodynamics (RHD) in axial symmetry assuming hydrostatic equilibrium perpendicular to the equatorial plane (1+1D) of a gaseous disk. The equations are formulated in conservative form on an adaptive grid and the corresponding fluxes are calculated by a spacial second order advection scheme. Self-gravity of the disk is included by solving the Possion equation. We test the resulting numerical method through comparison with a simplified analytical solution as well as through the long term viscous evolution of protoplanetary disk when due to viscosity matter is transported towards the central host star and the disk depletes. The importance of the inner boundary conditions on the structural behaviour of disks is demonstrated with several examples.

Alfv\'en-wave-driven 3D magnetohydrodynamic (MHD) models, which are increasingly used to predict stellar wind properties, contain unconstrained parameters and rely on low-resolution stellar magnetograms. We explore the effects of the input Alfv\'en wave energy flux and the surface magnetogram on the wind properties predicted by the Alfv\'en Wave Solar Model (AWSoM). We lowered the resolution of two solar magnetograms during solar cycle maximum and minimum using spherical harmonic decomposition. The Alfv\'en wave energy was altered based on non-thermal velocities determined from a far ultraviolet (FUV) spectrum of the solar twin 18 Sco. Additionally, low-resolution magnetograms of three solar analogues were obtained using Zeeman Doppler imaging (ZDI). Finally, the simulated wind properties were compared to Advanced Composition Explorer (ACE) observations. AWSoM simulations using well constrained input parameters taken from solar observations can reproduce the observed solar wind mass and angular momentum loss rates. The resolution of the magnetogram has a small impact on the wind properties and only during cycle maximum. However, variation in Alfv\'en wave energy influences the wind properties irrespective of the solar cycle activity level. Furthermore, solar wind simulations carried out using the low-resolution magnetogram of the three stars instead of the solar magnetogram could lead to an order of a magnitude difference in the simulated wind properties. The choice in Alfv\'en energy has a stronger influence on the wind output compared to the magnetogram resolution. The influence could be even stronger for stars whose input boundary conditions are not as well constrained as those of the Sun. Unsurprisingly, replacing the solar magnetogram with a stellar magnetogram could lead to completely inaccurate solar wind properties, and should be avoided in solar and stellar wind simulations.

When exposed to the high energy X-ray and ultraviolet radiation of a very active star, water vapor in the upper atmospheres of planets can be photodissociated and rapidly lost to space. In this paper, I study the chemical, thermal, and hydrodynamic processes in the upper atmospheres of terrestrial planets, concentrating on water vapor dominated atmospheres orbiting in the habitable zones of active stars. I consider different stellar activity levels and find very high levels of atmospheric escape in all cases, with the outflowing gas being dominated by atomic hydrogen and oxygen in both their neutral and ion forms. In the lower activity cases, I find that the accumulation of O$_2$ and increases in the D/H ratios in the atmospheres due to mass fractionation are possible, but in the higher activity cases no mass fractionation takes place. Connecting these results to stellar activity evolution tracks for solar mass stars, I show that huge amounts of water vapor can be lost, and both the losses and the amount of O$_2$ that can be accumulated in the atmosphere depend sensitively on the star's initial rotation rate. For an Earth-mass planet in the habitable zone of a low-mass M-dwarf, my results suggest that the accumulation of atmospheric O$_2$ is unlikely unless water loss can take place after the star's most active phase.

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.

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.

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 suggest to use the Hubble Space Telescople (HST) follow-up observations of the TESS targets for detecting possible plasma tori along the orbits of exoplanets orbiting M dwarfs. The source of the torus could be planetary volcanic activity due to tidal or electromagnetic induction heating. Fast losses to space for planets orbiting these active stars can lead to the lost material forming a torus along the planetary orbit, similar to the Io plasma torus. We show that such torus would be potentially detectable by the HST in the UV.

Planet atmospheric escape induced by high-energy stellar irradiation is a key phenomenon shaping the structure and evolution of planetary atmospheres. Therefore, the present-day properties of a planetary atmosphere are intimately connected with the amount of stellar flux received by a planet during its lifetime, thus with the evolutionary path of its host star. Using a recently developed analytic approximation based on hydrodynamic simulations for atmospheric escape rates, we track within a Bayesian framework the evolution of a planet as a function of stellar flux evolution history, constrained by the measured planetary radius, with the other system parameters as priors. We find that the ideal objects for this type of study are close-in sub-Neptune-like planets, as they are highly affected by atmospheric escape, and yet retain a significant fraction of their primordial hydrogen-dominated atmospheres. Furthermore, we apply this analysis to the HD3167 and K2-32 planetary systems. For HD3167, we find that the most probable irradiation level at 150 Myr was between 40 and 130 times solar, corresponding to a rotation period of 1.78^{+2.69}_{-1.23} days. For K2-32, we find a surprisingly low irradiation level ranging between half and four times solar at 150 Myr. Finally, we show that for multi-planet systems, our framework enables one to constrain poorly known properties of individual planets.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

Context. In tight binary star systems, tidal interactions can significantly influence the rotational and orbital evolution of both stars, and therefore their activity evolution. This can have strong effects on the atmospheric evolution of planets that are orbiting the two stars. Aims. In this paper, we aim to study the evolution of stellar rotation and of X-ray and ultraviolet (XUV) radiation in tight binary systems consisting of two solar mass stars and use our results to study planetary atmosphere evolution in the habitable zones of these systems. Methods. We have applied a rotation model developed for single stars to binary systems, taking into account the effects of tidal interactions on the rotational and orbital evolution of both stars. We used empirical rotation-activity relations to predict XUV evolution tracks for the stars, which we used to model hydrodynamic escape of hydrogen dominated atmospheres. Results. When significant, tidal interactions increase the total amount of XUV energy emitted, and in the most extreme cases by up to factor of $\sim$50. We find that in the systems that we study, habitable zone planets with masses of 1~M$_\oplus$ can lose huge hydrogen atmospheres due to the extended high levels of XUV emission, and the time that is needed to lose these atmospheres depends on the binary orbital separation.For some orbital separations, and when the stars are born as rapid rotators, it is also possible for tidal interactions to protect atmospheres from erosion by quickly spinning down the stars. For very small orbital separations, the loss of orbital angular momentum by stellar winds causes the two stars to merge. We suggest that the merging of the two stars could cause previously frozen planets to become habitable due to the habitable zone boundaries moving outwards.

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.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

The primary objectives of the ExoplANETS-A project are to: establish new knowledge on exoplanet atmospheres; establish new insight on influence of the host star on the planet atmosphere; disseminate knowledge, using online, web-based platforms. The project, funded under the EU's Horizon-2020 programme, started in January 2018 and has a duration ~3 years. We present an overview of the project, the activities concerning the host stars and some early results on the host stars.

Although the present-day Sun rotates too slowly to exhibit centrifugally-supported 'slingshot prominences', at some time during its past it may have formed these clouds of cool gas and ejected them into the interplanetary medium. We determine the time period for this behaviour by using a rotation evolution code to derive the properties of the formation and ejection of slingshot prominences during the lifetime of a star similar to our Sun. The mass, mass loss rate and rate of ejection of these prominences are calculated using the analytical expression derived in our previous work. We find that for stars with an initial rotation rate larger than $4.6\, \Omega_\odot$, about half of all solar mass stars, slingshot prominences will be present even after the star reaches the main sequence phase. In a fast rotator, this means that prominences can form until the star reaches $\sim 800$ Myr old. Our results also indicate that the mass and lifetime of this type of prominence have maximum values when the star reaches the ZAMS at an age of $\sim 40$ Myr for a solar mass star.

Studies of planetary atmospheric composition, variability, and evolution require appropriate theoretical and numerical tools to estimate key atmospheric parameters, among which the mass-loss rate is often the most important. In evolutionary studies, it is common to use the energy-limited formula, which is attractive for its simplicity but ignores important physical effects and can be inaccurate in many cases. To overcome this problem, we consider a recently developed grid of about 7000 one-dimensional upper-atmosphere hydrodynamic models computed for a wide range of planets with hydrogen-dominated atmospheres from which we extract the mass-loss rates. The grid boundaries are [1:39] MEARTH in planetary mass, [1:10] REARTH in planetary radius, [300:2000] K in equilibrium temperature, [0.4:1.3] MSUN in host star's mass, [0.002:1.3] au in orbital separation, and about [10^{26}:5*10^{30}] erg/s in stellar X-ray and extreme ultraviolet luminosity. We then derive an analytical expression for the atmospheric mass-loss rates based on a fit to the values obtained from the grid. The expression provides the mass-loss rates as a function of planetary mass, planetary radius, orbital separation, and incident stellar high-energy flux. We show that this expression is a significant improvement to the energy-limited approximation for a wide range of planets. The analytical expression presented here enables significantly more accurate planetary evolution computations without increasing computing time.

There is growing observational and theoretical evidence suggesting that atmospheric escape is a key driver of planetary evolution. Commonly, planetary evolution models employ simple analytic formulae (e.g., energy limited escape) that are often inaccurate, and more detailed physical models of atmospheric loss usually only give snapshots of an atmosphere's structure and are difficult to use for evolutionary studies. To overcome this problem, we upgrade and employ an already existing upper atmosphere hydrodynamic code to produce a large grid of about 7000 models covering planets with masses 1 - 39 Earth mass with hydrogen-dominated atmospheres and orbiting late-type stars. The modeled planets have equilibrium temperatures ranging between 300 and 2000 K. For each considered stellar mass, we account for three different values of the high-energy stellar flux (i.e., low, moderate, and high activity). For each computed model, we derive the atmospheric temperature, number density, bulk velocity, X-ray and EUV (XUV) volume heating rates, and abundance of the considered species as a function of distance from the planetary center. From these quantities, we estimate the positions of the maximum dissociation and ionisation, the mass-loss rate, and the effective radius of the XUV absorption. We show that our results are in good agreement with previously published studies employing similar codes. We further present an interpolation routine capable to extract the modelling output parameters for any planet lying within the grid boundaries. We use the grid to identify the connection between the system parameters and the resulting atmospheric properties. We finally apply the grid and the interpolation routine to estimate atmospheric evolutionary tracks for the close-in, high-density planets CoRoT-7 b and HD219134 b,c...

Planets that form early enough to be embedded in the circumstellar gas disk accumulate thick atmospheres of nebular gas. Models of these atmospheres need to specify the surface luminosity (i.e. energy loss rate) of the planet. This luminosity is usually associated with a continuous inflow of solid bodies, where the gravitational energy released from these bodies is the source of energy. However, if these bodies release energy in the atmosphere instead of at the surface, this assumption might not be justified. Our aim is to explore the interactions of infalling planetesimals with primordial atmospheres at an embedded phase of evolution. We investigate effects of atmospheric interaction on the planetesimals (mass loss) and the atmosphere (heating/cooling). We used atmospheric parameters from a snapshot of time-dependent evolution simulations for embedded atmospheres and simulated purely radial, infall events of siliceous planetesimals in a 1D, explicit code. We implemented energy transfer between friction, radiation transfer by the atmosphere and the body and thermal ablation; this gives us the possibility to examine the effects on the planetesimals and the atmosphere. We find that a significant amount of gravitational energy is indeed dissipated into the atmosphere, especially for larger planetary cores, which consequently cannot contribute to the atmospheric planetary luminosity. Furthermore, we examine that planetesimal infall events for cores, $M_\mathrm{C} > 2$M$_{\oplus}$, which actually result in a local cooling of the atmosphere; this is totally in contradiction with the classical model.

Context: The thermal and chemical structures of the upper atmospheres of planets crucially influence losses to space and must be understood to constrain the effects of losses on atmospheric evolution. Aims: We develop a 1D first-principles hydrodynamic atmosphere model that calculates atmospheric thermal and chemical structures for arbitrary planetary parameters, chemical compositions, and stellar inputs. We apply the model to study the reaction of the Earth's upper atmosphere to large changes in the CO$_2$ abundance and to changes in the input solar XUV field due to the Sun's activity evolution from 3~Gyr in the past to 2.5~Gyr in the future. Methods: For the thermal atmosphere structure, we consider heating from the absorption of stellar X-ray, UV, and IR radiation, heating from exothermic chemical reactions, electron heating from collisions with non-thermal photoelectrons, Joule heating, cooling from IR emission by several species, thermal conduction, and energy exchanges between the neutral, ion, and electron gases. For the chemical structure, we consider $\sim$500 chemical reactions, including 56 photoreactions, eddy and molecular diffusion, and advection. In addition, we calculate the atmospheric structure by solving the hydrodynamic equations. To solve the equations in our model, we develop the Kompot code and provide detailed descriptions of the numerical methods used in the appendices. Results: We verify our model by calculating the structures of the upper atmospheres of the modern Earth and Venus. By varying the CO$_2$ abundances at the lower boundary (65~km) of our Earth model, we show that the atmospheric thermal structure is significantly altered. [Abstract Truncated]

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.

The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the search for life in the Universe by combining expertise from four NASA science disciplines including heliophysics, astrophysics, planetary science and Earth science. The NASA Nexus for Exoplanetary System Science (NExSS) has provided an efficient environment for such interdisciplinary studies. Solar flares, coronal mass ejections and solar energetic particles produce disturbances in interplanetary space collectively referred to as space weather, which interacts with the Earth upper atmosphere and causes dramatic impact on space and ground-based technological systems. Exoplanets within close in habitable zones around M dwarfs and other active stars are exposed to extreme ionizing radiation fluxes, thus making exoplanetary space weather (ESW) effects a crucial factor of habitability. In this paper, we describe the recent developments and provide recommendations in this interdisciplinary effort with the focus on the impacts of ESW on habitability, and the prospects for future progress in searching for signs of life in the Universe as the outcome of the NExSS workshop held in Nov 29 - Dec 2, 2016, New Orleans, LA. This is one of five Life Beyond the Solar System white papers submitted by NExSS to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe.

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.

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.

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.

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.

The strong, nearly wavelength-independent absorption cross section of aerosols produces featureless exoplanet transmission spectra, limiting our ability to characterize their atmospheres. Here we show that even in the presence of featureless spectra, we can still characterize certain atmospheric properties. Specifically, we constrain the upper and lower pressure boundaries of aerosol layers, and present plausible composition candidates. We study the case of the bloated Saturn-mass planet WASP-49b, where near-infrared observations reveal a flat transmission spectrum between 0.7 and 1.0 {\microns}. First, we use a hydrodynamic upper-atmosphere code to estimate the pressure reached by the ionizing stellar high-energy photons at $10^{-8}$ bar, setting the upper pressure boundary where aerosols could exist. Then, we combine HELIOS and Pyrat Bay radiative-transfer models to constrain the temperature and photospheric pressure of atmospheric aerosols, in a Bayesian framework. For WASP-49b, we constrain the transmission photosphere (hence, the aerosol deck boundaries) to pressures above $10^{-5}$ bar (100$\times$ solar metallicity), $10^{-4}$ bar (solar), and $10^{-3}$ bar (0.1$\times$ solar) as lower boundary, and below $10^{-7}$ bar as upper boundary. Lastly, we compare condensation curves of aerosol compounds with the planet's pressure-temperature profile to identify plausible condensates responsible for the absorption. Under these circumstances, we find as candidates: Na$_{2}$S (at 100$\times$ solar metallicity); Cr and MnS (at solar and 0.1$\times$ solar); and forsterite, enstatite, and alabandite (at 0.1$\times$ solar).

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.

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.

We present a uniform analysis of the atmospheric escape rate of Neptune-like planets with estimated radius and mass (restricted to $M_{\rm p}<30\,M_{\oplus}$). For each planet we compute the restricted Jeans escape parameter, $\Lambda$, for a hydrogen atom evaluated at the planetary mass, radius, and equilibrium temperature. Values of $\Lambda\lesssim20$ suggest extremely high mass-loss rates. We identify 27 planets (out of 167) that are simultaneously consistent with hydrogen-dominated atmospheres and are expected to exhibit extreme mass-loss rates. We further estimate the mass-loss rates ($L_{\rm hy}$) of these planets with tailored atmospheric hydrodynamic models. We compare $L_{\rm hy}$ to the energy-limited (maximum-possible high-energy driven) mass-loss rates. We confirm that 25 planets (15\% of the sample) exhibit extremely high mass-loss rates ($L_{\rm hy}>0.1\,M_{\oplus}{\rm Gyr}^{-1}$), well in excess of the energy-limited mass-loss rates. This constitutes a contradiction, since the hydrogen envelopes cannot be retained given the high mass-loss rates. We hypothesize that these planets are not truly under such high mass-loss rates. Instead, either hydrodynamic models overestimate the mass-loss rates, transit-timing-variation measurements underestimate the planetary masses, optical transit observations overestimate the planetary radii (due to high-altitude clouds), or Neptunes have consistently higher albedos than Jupiter planets. We conclude that at least one of these established estimations/techniques is consistently producing biased values for Neptune planets. Such an important fraction of exoplanets with misinterpreted parameters can significantly bias our view of populations studies, like the observed mass--radius distribution of exoplanets for example.

Aims: We study the acceleration of the stellar winds of rapidly rotating low mass stars and the transition between the slow magnetic rotator and fast magnetic rotator regimes. We aim to understand the properties of stellar winds in the fast magnetic rotator regime and the effects of magneto-centrifugal forces on wind speeds and mass loss rates. Methods: We extend the solar wind model of Johnstone et al. (2015b) to 1D magnetohydrodynamic (MHD) simulations of the winds of rotating stars. We test two assumptions for how to scale the wind temperature to other stars and assume the mass loss rate scales as Mdot ~ Rstar^2 OmegaStar^1.33 Mstar^-3.36, in the unsaturated regime, as estimated by Johnstone et al. (2015a). Results: For 1.0 Msun stars, the winds can be accelerated to several thousand km/s, and the effects of magneto-centrifugal forces are much weaker for lower mass stars. We find that the different assumptions for how to scale the wind temperature to other stars lead to significantly different mass loss rates for the rapid rotators. If we assume a constant temperature, the mass loss rates of solar mass stars do not saturate at rapid rotation, which we show to be inconsistent with observed rotational evolution. If we assume the wind temperatures scale positively with rotation, the mass loss rates are only influenced significantly at rotation rates above 75 OmegaSun. We suggest that models with increasing wind speed for more rapid rotators are preferable to those that assume a constant wind speed. If this conclusion is confirmed by more sophisticated wind modelling. it might provide an interesting observational constraint on the properties of stellar winds.

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.

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.

Terrestrial planets formed within gaseous protoplanetary disks can accumulate significant hydrogen envelopes. The evolution of such an atmosphere due to XUV driven evaporation depends on the activity evolution of the host star, which itself depends sensitively on its rotational evolution, and therefore on its initial rotation rate. In this letter, we derive an easily applicable method for calculating planetary atmosphere evaporation that combines models for a hydrostatic lower atmosphere and a hydrodynamic upper atmosphere. We show that the initial rotation rate of the central star is of critical importance for the evolution of planetary atmospheres and can determine if a planet keeps or loses its primordial hydrogen envelope. Our results highlight the need for a detailed treatment of stellar activity evolution when studying the evolution of planetary atmospheres.

Stellar magnetism, explorable via polarimetry, is a crucial driver of activity, ionization, photodissociation, chemistry and winds in stellar environments. Thus it has an important impact on the atmospheres and magnetospheres of surrounding planets. Modeling of stellar magnetic fields and their winds is extremely challenging, both from the observational and the theoretical points of view, and only recent ground breaking advances in observational instrumentation - as were discussed during this Symposium - and a deeper theoretical understanding of magnetohydrodynamic processes in stars enable us to model stellar magnetic fields and winds and the resulting influence on surrounding planets in more and more detail. We have initiated a national and international research network (NFN): 'Pathways to Habitability - From Disks to Active Stars, Planets to Life', to address questions on the formation and habitability of environments in young, active stellar/planetary systems. In this contribution we discuss the work we are carrying out within this project and focus on how stellar magnetic fields, their winds and the relation to stellar rotation can be assessed observationally with relevant techniques such as Zeeman Doppler Imaging (ZDI), field extrapolation and wind simulations.

We study the short-periodic component of stellar activity with a cycle periods Pcyc up to 1000 days using the Kepler mission photometry of fast-rotating (rotational periods from 1 to 4 days) stars with spectra of M4V to F3V. Applying the originally developed two non-spectral methods, we measured the effective period of stellar cycles in 462 objects. The obtained results are in accordance with previous measurements by Vida et al. (2014), do not seem to result from a beating effect. The performed measurements of Pcyc cluster in a specific branch which covers the previously unstudied region in the Saar-Brandenburg (1999) diagram, and connects the branch of inactive stars with the area populated by super-active objects. It is shown that the formation of the discovered branch is due to the alpha-quenching effect, which saturates the magnetic dynamo and decreases the cycle periods with the increase of inverted Rossby number. This finding is important in the context of the discussion on catastrophic quenching and other heuristic approximations of the non-linear alpha-effect.

Aims: We study the evolution of stellar rotation and wind properties for low-mass main-sequence stars. Our aim is to use rotational evolution models to constrain the mass loss rates in stellar winds and to predict how their properties evolve with time on the main-sequence. Methods: We construct a rotational evolution model that is driven by observed rotational distributions of young stellar clusters. Fitting the free parameters in our model allows us to predict how wind mass loss rate depends on stellar mass, radius, and rotation. We couple the results to the wind model developed in Paper I of this series to predict how wind properties evolve on the main-sequence. Results: We estimate that wind mass loss rate scales with stellar parameters as $\dot{M}_\star \propto R_\star^2 \Omega_\star^{1.33} M_\star^{-3.36}$. We estimate that at young ages, the solar wind likely had a mass loss rate that is an order of magnitude higher than that of the current solar wind. This leads to the wind having a higher density at younger ages; however, the magnitude of this change depends strongly on how we scale wind temperature. Due to the spread in rotation rates, young stars show a large range of wind properties at a given age. This spread in wind properties disappears as the stars age. Conclusions: There is a large uncertainty in our knowledge of the evolution of stellar winds on the main-sequence, due both to our lack of knowledge of stellar winds and the large spread in rotation rates at young ages. Given the sensitivity of planetary atmospheres to stellar wind and radiation conditions, these uncertainties can be significant for our understanding of the evolution of planetary environments.

Aims. We study the X-ray emission of low-mass main-sequence stars to derive a reliable general scaling law between coronal temperature and the level of X-ray activity. Methods. We collect ROSAT measurements of hardness ratios and X-ray luminosities for a large sample of stars to derive which stellar X-ray emission parameter is most closely correlated with coronal temperature. We calculate average coronal temperatures for a sample of 24 low-mass main-sequence stars with measured emission measure distributions (EMDs) collected from the literature. These EMDs are based on high-resolution X-ray spectra measured by XMM-Newton and Chandra. Results. We confirm that there is one universal scaling relation between coronal average temperature and surface X-ray flux, Fx, that applies to all low-mass main-sequence stars. We find that coronal temperature is related to Fx by Tcor=0.11 Fx^0.26, where Tcor is in MK and Fx is in erg/s/cm^2.

Aims. We aim to describe the pre-main sequence and main-sequence evolution of X-ray and extreme-ultaviolet radiation of a solar mass star based on its rotational evolution starting with a realistic range of initial rotation rates. Methods. We derive evolutionary tracks of X-ray radiation based on a rotational evolution model for solar mass stars and the rotation-activity relation. We compare these tracks to X-ray luminosity distributions of stars in clusters with different ages. Results. We find agreement between the evolutionary tracks derived from rotation and the X-ray luminosity distributions from observations. Depending on the initial rotation rate, a star might remain at the X-ray saturation level for very different time periods, approximately from 10 Myr to 300 Myr for slow and fast rotators, respectively. Conclusions. Rotational evolution with a spread of initial conditions leads to a particularly wide distribution of possible X-ray luminosities in the age range of 20 to 500 Myrs, before rotational convergence and therefore X-ray luminosity convergence sets in. This age range is crucial for the evolution of young planetary atmospheres and may thus lead to very different planetary evolution histories.

In this paper, we estimate the X-ray emission from close-in exoplanets. We show that the Solar/Stellar Wind Charge Exchange Mechanism (SWCX) which produces soft X-ray emission is very effective for hot Jupiters. In this mechanism, X-ray photons are emitted as a result of the charge exchange between heavy ions in the solar wind and the atmospheric neutral particles. In the Solar System, comets produce X-rays mostly through the SWCX mechanism, but it has also been shown to operate in the heliosphere, in the terrestrial magnetosheath, and on Mars, Venus and Moon. Since the number of emitted photons is proportional to the solar wind mass flux, this mechanism is not very effective for the Solar system giants. Here we present a simple estimate of the X-ray emission intensity that can be produced by close-in extrasolar giant planets due to charge exchange with the heavy ions of the stellar wind. Using the example of HD~209458b, we show that this mechanism alone can be responsible for an X-ray emission of $\approx 10^{22}$~erg~s$^{-1}$, which is $10^6$ times stronger than the emission from the Jovian aurora. We discuss also the possibility to observe the predicted soft X-ray flux of hot Jupiters and show that despite high emission intensities they are unobservable with current facilities.

Aims: We develop a method for estimating the properties of stellar winds for low-mass main-sequence stars between masses of 0.4 and 1.1 solar masses at a range of distances from the star. Methods: We use 1D thermal pressure driven hydrodynamic wind models run using the Versatile Advection Code. Using in situ measurements of the solar wind, we produce models for the slow and fast components of the solar wind. We consider two radically different methods for scaling the base temperature of the wind to other stars: in Model A, we assume that wind temperatures are fundamentally linked to coronal temperatures, and in Model B, we assume that the sound speed at the base of the wind is a fixed fraction of the escape velocity. In Paper II of this series, we use observationally constrained rotational evolution models to derive wind mass loss rates. Results: Our model for the solar wind provides an excellent description of the real solar wind far from the solar surface, but is unrealistic within the solar corona. We run a grid of 1200 wind models to derive relations for the wind properties as a function of stellar mass, radius, and wind temperature. Using these results, we explore how wind properties depend on stellar mass and rotation. Conclusions: Based on our two assumptions about the scaling of the wind temperature, we argue that there is still significant uncertainty in how these properties should be determined. Resolution of this uncertainty will probably require both the application of solar wind physics to other stars and detailed observational constraints on the properties of stellar winds. In the final section of this paper, we give step by step instructions for how to apply our results to calculate the stellar wind conditions far from the stellar surface.

Context. In binary star systems, the winds from the two components impact each other, leading to strong shocks and regions of enhanced density and temperature. Potentially habitable circumbinary planets must continually be exposed to these interactions regions. Aims. We study, for the first time, the interactions between winds from low-mass stars in a binary system, to show the wind conditions seen by potentially habitable circumbinary planets. Methods. We use the advanced 3D numerical hydrodynamic code Nurgush to model the wind interactions of two identical winds from two solar mass stars with circular orbits and a binary separation of 0.5 AU. As input into this model, we use a 1D hydrodynamic simulation of the solar wind, run using the Versatile Advection Code. We derive the locations of stable and habitable orbits in this system to explore what wind conditions potentially habitable planets will be exposed to during their orbits. Results. Our wind interaction simulations result in the formation of two strong shock waves separated by a region of enhanced density and temperature. The wind-wind interaction region has a spiral shape due to Coriolis forces generated by the orbital motions of the two stars. The stable and habitable zone in this system extends from approximately 1.4 AU to 2.4 AU. (TRUNCATED)

We study the interactions between stellar wind and the extended hydrogen-dominated upper atmospheres of planets and the resulting escape of planetary pick-up ions from the 5 "super-Earths" in the compact Kepler-11 system and compare the escape rates with the efficiency of the thermal escape of neutral hydrogen atoms. Assuming the stellar wind of Kepler-11 is similar to the solar wind, we use a polytropic 1D hydrodynamic wind model to estimate the wind properties at the planetary orbits. We apply a Direct Simulation Monte Carlo Model to model the hydrogen coronae and the stellar wind plasma interaction around Kepler-11b-f within a realistic expected heating efficiency range of 15-40%. The same model is used to estimate the ion pick-up escape from the XUV heated and hydrodynamically extended upper atmospheres of Kepler-11b-f. From the interaction model we study the influence of possible magnetic moments, calculate the charge exchange and photoionization production rates of planetary ions and estimate the loss rates of pick-up H+ ions for all five planets. We compare the results between the five "super-Earths" and in a more general sense also with the thermal escape rates of the neutral planetary hydrogen atoms. Our results show that for all Kepler-11b-f exoplanets, a huge neutral hydrogen corona is formed around the planet. The non-symmetric form of the corona changes from planet to planet and is defined mostly by radiation pressure and gravitational effects. Non-thermal escape rates of pick-up ionized hydrogen atoms for Kepler-11 "super-Earths" vary between approximately 6.4e30 1/s and 4.1e31 1/s depending on the planet's orbital location and assumed heating efficiency. These values correspond to non-thermal mass loss rates of approximately 1.07e7 g/s and 6.8e7 g/s respectively, which is a few percent of the thermal escape rates.

The magnetic fields of young stars set their coronal properties and control their spin evolution via the star-disc interaction and outflows. Using 14 magnetic maps of 10 classical T Tauri stars (CTTSs) we investigate their closed X-ray emitting coronae, their open wind-bearing magnetic fields, and the geometry of magnetospheric accretion flows. The magnetic fields of all the CTTSs are multipolar. Stars with simpler (more dipolar) large-scale magnetic fields have stronger fields, are slower rotators, and have larger X-ray emitting coronae compared to stars with more complex large-scale magnetic fields. The field complexity controls the distribution of open and closed field regions across the stellar surface, and strongly influences the location and shapes of accretion hot spots. However, the higher order field components are of secondary importance in determining the total unsigned open magnetic flux, which depends mainly on the strength of the dipole component and the stellar surface area. Likewise, the dipole component alone provides an adequate approximation of the disc truncation radius. For some stars, the pressure of the hot coronal plasma dominates the stellar magnetic pressure and forces open the closed field inside the disc truncation radius. This is significant as accretion models generally assume that the magnetic field has a closed geometry out to the inner disc edge.

As part of the project Pathways to Habitability (http://path.univie.ac.at/), we study the properties of the stellar winds of low-mass and Sun-like stars, and their influences on the atmospheres of potentially habitable planets. For this purpose, we combine mapping of stellar magnetic fields with magnetohydrodynamic wind models.

This summary reports on papers presented at the Cool Stars-16 meeting in the splinter session "Solar and Stellar flares." Although many topics were discussed, the main themes were the commonality of interests, and of physics, between the solar and stellar flare communities, and the opportunities for important new observations in the near future.

Most stellar flares' soft X-ray lightcurves possess a `typical' morphology, which consists of a rapid rise followed by a slow exponential decay. However, a study of 216 of the brightest flares on 161 pre-main sequence stars, observed during the Chandra Orion-Ultradeep Project (COUP), showed that many flare lightcurves depart from this typical morphology. While this can be attributed to the superposition of multiple typical flares, we explore the possibility that the time-variable eclipsing of flares by their host stars may also be an important factor. We assume each flare is contained within a single, uniform plasma density magnetic loop and specify the intrinsic variation of the flare's emission measure with time. We consider rotational eclipse by the star itself, but also by circumstellar discs and flare-associated prominences. Based on this simple model, we generate a set of flares similar to those observed in the COUP database. Many eclipses simply reduce the flare's maximum emission measure or decay time. We conclude therefore that eclipses often pass undetected, but usually have only a modest influence on the flare emission measure profile and hence the derived loop lengths. We show that eclipsing can easily reproduce the observed atypical flare morphologies. The number of atypical modelled flare morphologies is however much less than that found in the COUP sample. The large number of observed atypical flare morphologies, therefore, must be attributed to other processes such as multiple flaring loops.

Significant progress has been made recently in our understanding of the structure of stellar magnetic fields, thanks to advances in detection methods such as Zeeman-Doppler Imaging. The extrapolation of this surface magnetic field into the corona has provided 3D models of the coronal magnetic field and plasma. This method is sensitive mainly to the magnetic field in the bright regions of the stellar surface. The dark (spotted) regions are censored because the Zeeman signature there is suppressed. By modelling the magnetic field that might have been contained in these spots, we have studied the effect that this loss of information might have on our understanding of the coronal structure. As examples, we have chosen two stars (V374 peg and AB Dor) that have very different magnetograms and patterns of spot coverage. We find that the effect of the spot field depends not only on the relative amount of flux in the spots, but also its distribution across the stellar surface. For a star such as AB Dor with a high spot coverage and a large polar spot, at its greatest effect the spot field may almost double the fraction of the flux that is open (hence decreasing the spindown time) while at the same time increasing the X-ray emission measure by two orders of magnitude and significantly affecting the X-ray rotational modulation.

Recent advances in spectropolarimetry have allowed the reconstruction of stellar coronal magnetic fields. This uses Zeeman-Doppler magnetograms of the surface magnetic field as a lower boundary condition. The ZDI maps, however, suffer from the absence of information about the magnetic field over regions of the surface due to the presence of dark starspots and portions of the surface out of view due to a tilt in the rotation axis. They also suffer from finite resolution which leads to small scale field structures being neglected. This paper explores the effects of this loss of information on the extrapolated coronal fields. For this we use simulated stellar surface magnetic maps for two hypothetical stars. Using the potential field approximation, the coronal fields and emission measures are calculated. This is repeated for the cases of missing information due to, (i) starspots, (ii) a large area of the stellar surface out of view, (iii) a finite resolution. The largest effect on the magnetic field structure arises when a significant portion of the stellar surface remains out of view. This changes the nature of the field lines that connect to this obscured hemisphere. Nonetheless, the field structure in the visible hemisphere is reliably reproduced. Thus the calculation of the locations and surface filling factors of accretion funnels is reasonably well reproduced for the observed hemisphere. The decrease with height of the magnetic pressure, which is important in calculating disc truncation radii for accreting stars, is also largely unaffected in the equatorial plane. The fraction of surface flux that is open and therefore able to supply angular momentum loss in a wind, however, is often overestimated in the presence of missing flux.