We investigate the space environment conditions near the Earth-size planet TOI~700~d using a set of numerical models for the stellar corona and wind, the planetary magnetosphere, and the planetary ionosphere. We drive our simulations using a scaled-down stellar input and a scaled-up solar input in order to obtain two independent solutions. We find that for the particular parameters used in our study, the stellar wind conditions near the planet are not very extreme -- slightly stronger than that near the Earth in terms of the stellar wind ram pressure and the intensity of the interplanetary magnetic field. Thus, the space environment near TOI700-d may not be extremely harmful to the planetary atmosphere, assuming the planet resembles the Earth. Nevertheless, we stress that the stellar input parameters and the actual planetary parameters are unconstrained, and different parameters may result in a much greater effect on the atmosphere of TOI700-d. Finally, we compare our results to solar wind measurements in the solar system and stress that modest stellar wind conditions may not guarantee atmospheric retention of exoplanets.

We present a 0.3-5 micron transmission spectrum of the hot Jupiter HAT-P-32Ab observed with the Space Telescope Imaging Spectrograph (STIS) and Wide Field Camera 3 (WFC3) instruments mounted on the Hubble Space Telescope, combined with Spitzer Infrared Array Camera (IRAC) photometry. The spectrum is composed of 51 spectrophotometric bins with widths ranging between 150 and 400 \AA, measured to a median precision of 215 ppm. Comparisons of the observed transmission spectrum to a grid of 1D radiative-convective equilibrium models indicate the presence of clouds/hazes, consistent with previous transit observations and secondary eclipse measurements. To provide more robust constraints on the planet's atmospheric properties, we perform the first full optical to infrared retrieval analysis for this planet. The retrieved spectrum is consistent with a limb temperature of 1248$\pm$92 K, a thick cloud deck, enhanced Rayleigh scattering, and $\sim$10x solar H2O abundance. We find log($Z/Z_{\odot}$) = 2.41$_{-0.07}^{+0.06}$, in agreement with the mass-metallicity relation derived for the Solar System.

Coronal mass ejections (CMEs) on stars other than the Sun have proven very difficult to detect. One promising pathway lies in the detection of type II radio bursts. Their appearance and distinctive properties are associated with the development of an outward propagating CME-driven shock. However, dedicated radio searches have not been able to identify these transient features in other stars. Large Alfv\'en speeds and the magnetic suppression of CMEs in active stars have been proposed to render stellar eruptions "radio-quiet". Employing 3D magnetohydrodynamic simulations, we study here the distribution of the coronal Alfv\'en speed, focusing on two cases representative of a young Sun-like star and a mid-activity M-dwarf (Proxima Centauri). These results are compared with a standard solar simulation and used to characterize the shock-prone regions in the stellar corona and wind. Furthermore, using a flux-rope eruption model, we drive realistic CME events within our M-dwarf simulation. We consider eruptions with different energies to probe the regimes of weak and partial CME magnetic confinement. While these CMEs are able to generate shocks in the corona, those are pushed much farther out compared to their solar counterparts. This drastically reduces the resulting type II radio burst frequencies down to the ionospheric cutoff, which impedes their detection with ground-based instrumentation.

The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox.

We perform atmospheric retrievals on the full optical to infrared ($0.3-5 \, \mu \mathrm{m}$) transmission spectrum of the inflated hot Jupiter WASP-52b by combining HST/STIS, WFC3 IR, and Spitzer/IRAC observations. As WASP-52 is an active star which shows both out-of-transit photometric variability and starspot crossings during transits, we account for the contribution of non-occulted active regions in the retrieval. We recover a $0.1-10\times$ solar atmospheric composition, in agreement with core accretion predictions for giant planets, and a weak contribution of aerosols. We also obtain a $<3000$ K temperature for the starspots, a measure which is likely affected by the models used to fit instrumental effects in the transits, and a 5% starspot fractional coverage, compatible with expectations for the host star's spectral type. Such constraints on the planetary atmosphere and on the activity of its host star will inform future JWST GTO observations of this target.

The stellar magnetic field completely dominates the environment around late-type stars. It is responsible for driving the coronal high-energy radiation (e.g. EUV/X-rays), the development of stellar winds, and the generation transient events such as flares and coronal mass ejections (CMEs). While progress has been made for the first two processes, our understanding of the eruptive behavior in late-type stars is still very limited. One example of this is the fact that despite the frequent and highly energetic flaring observed in active stars, direct evidence for stellar CMEs is almost non-existent. Here we discuss realistic 3D simulations of stellar CMEs, analyzing their resulting properties in contrast with solar eruptions, and use them to provide a common framework to interpret the available stellar observations. Additionally, we present results from the first 3D CME simulations in M-dwarf stars, with emphasis on possible observable signatures imprinted in the stellar corona.

Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can be used locally only, since the limitations in available computational resources do not allow for their use throughout a global computational domain. In the present work, we present a new coupled scheme which allows to switch blocks in the block-adaptive grids from fluid MHD to hybrid simulations, without modifying the self-consistent computation of the electromagnetic fields acting on fluids (in MHD simulation) or charged ion macroparticles (in hybrid simulation). In this way, the hybrid scheme can refine the description in specified regions of interest without compromising the efficiency of the global MHD code.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

We report the results of the first state-of-the-art numerical simulations of Coronal Mass Ejections (CMEs) taking place in realistic magnetic field configurations of moderately active M-dwarf stars. Our analysis indicates that a clear, novel, and observable, coronal response is generated due to the collapse of the eruption and its eventual release into the stellar wind. Escaping CME events, weakly suppressed by the large-scale field, induce a flare-like signature in the emission from coronal material at different temperatures due to compression and associated heating. Such flare-like profiles display a distinctive temporal evolution in their Doppler shift signal (from red to blue), as the eruption first collapses towards the star and then perturbs the ambient magnetized plasma on its way outwards. For stellar fields providing partial confinement, CME fragmentation takes place, leading to rise and fall flow patterns which resemble the solar coronal rain cycle. In strongly suppressed events, the response is better described as a gradual brightening, in which the failed CME is deposited in the form of a coronal rain cloud leading to a much slower rise in the ambient high-energy flux by relatively small factors ($\sim2-3$). In all the considered cases (escaping/confined) a fractional decrease in the emission from mid-range coronal temperature plasma occurs, similar to the coronal dimming events observed on the Sun. Detection of the observational signatures of these CME-induced features requires a sensitive next generation X-ray space telescope.

The sun's chromosphere is a highly dynamic, partially-ionized region where spicules (hot jets of plasma) form. Here we present a two-fluid MHD model to study the chromosphere, which includes ion-neutral interaction and frictional heating. Our simulation recovers a magnetic canopy shape that forms quickly, but is also quickly disrupted by the formation of a jet. Our simulation produces a shock self-consistently, where the jet is driven by the frictional heating, which is much greater than the ohmic heating. Thus, our simulation demonstrates that the jet could be driven purely by thermal effects due to ion-neutral collisions and not by magnetic reconnection. We plan to improve the model to include photo-chemical effects and radiation.

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.

Solar CMEs and flares have a statistically well defined relation, with more energetic X-ray flares corresponding to faster and more massive CMEs. How this relation extends to more magnetically active stars is a subject of open research. Here, we study the most probable stellar CME candidates associated with flares captured in the literature to date, all of which were observed on magnetically active stars. We use a simple CME model to derive masses and kinetic energies from observed quantities, and transform associated flare data to the GOES 1--8~\AA\ band. Derived CME masses range from $\sim 10^{15}$ to $10^{22}$~g. Associated flare X-ray energies range from $10^{31}$ to $10^{37}$~erg. Stellar CME masses as a function of associated flare energy generally lie along or below the extrapolated mean for solar events. In contrast, CME kinetic energies lie below the analogous solar extrapolation by roughly two orders of magnitude, indicating approximate parity between flare X-ray and CME kinetic energies. These results suggest that the CMEs associated with very energetic flares on active stars are more limited in terms of the ejecta velocity than the ejecta mass, possibly because of the restraining influence of strong overlying magnetic fields and stellar wind drag. Lower CME kinetic energies and velocities present a more optimistic scenario for the effects of CME impacts on exoplanets in close proximity to active stellar hosts.

A physically realistic stellar wind model based on Alfv\'en wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distributions of three different stars that constitute admissible magnetic proxies of this object. Under these considerations, Barnard's Star b experiences less intense wind pressure than the much more close-in planet Proxima~b and the planets of the TRAPPIST-1 system. The milder wind conditions are more a result of its much greater orbital distance rather than in differences in the surface magnetic field strengths of the host stars. The dynamic pressure experienced by the planet is comparable to present-day Earth values, but it can undergo variations by factors of several during current sheet crossings in each orbit. The magnetospause standoff distance would be $\sim$\,$20 - 40$\,\% smaller than that of the Earth for an equivalent planetary magnetic field strength.

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.

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.

Planets in close proximity to their parent star, such as those in the habitable zones around M dwarfs, could be subject to particularly high doses of particle radiation. We have carried out test-particle simulations of ~GeV protons to investigate the propagation of energetic particles accelerated by flares or travelling shock waves within the stellar wind and magnetic field of a TRAPPIST-1-like system. Turbulence was simulated with small-scale magnetostatic perturbations with an isotropic power spectrum. We find that only a few percent of particles injected within half a stellar radius from the stellar surface escape, and that the escaping fraction increases strongly with increasing injection radius. Escaping particles are increasingly deflected and focused by the ambient spiralling magnetic field as the superimposed turbulence amplitude is increased. In our TRAPPIST-1-like simulations, regardless of the angular region of injection, particles are strongly focused onto two caps within the fast wind regions and centered on the equatorial planetary orbital plane. Based on a scaling relation between far-UV emission and energetic protons for solar flares applied to M dwarfs, the innermost putative habitable planet, TRAPPIST-1e, is bombarded by a proton flux up to 6 orders of magnitude larger than experienced by the present-day Earth. We note two mechanisms that could strongly limit EP fluxes from active stars: EPs from flares are contained by the stellar magnetic field; and potential CMEs that might generate EPs at larger distances also fail to escape.

GJ 3470b is a warm Neptune transiting an M-dwarf star at the edge of the evaporation desert. It offers the possibility of investigating how low-mass, close-in exoplanets evolve under the irradiation from their host stars. We observed three transits of GJ 3470b in the Lyman-alpha line with the Hubble Space Telescope (HST) as part of the Panchromatic Comparative Exoplanet Treasury (PanCET) program. Absorption signatures are detected with similar properties in all three independent epochs, with absorption depths of 35+-7% in the blue wing of the line, and 23+-5% in the red wing. The repeatability of these signatures, their phasing with the planet transit, and the radial velocity of the absorbing gas allow us to conclude that there is an extended upper atmosphere of neutral hydrogen around GJ 3470 b. We determine from our observations the stellar radiation pressure and XUV irradiation from GJ 3470 and use them to perform numerical simulations of the upper atmosphere of GJ 3470b with the EVaporating Exoplanets (EVE) code. The unusual redshifted signature can be explained by the damping wings of dense layers of neutral hydrogen that extend beyond the Roche lobe and are elongated in the direction of the planet motion. This structure could correspond to a shocked layer of planetary material formed by the collision of the expanding thermosphere with the wind of the star. The blueshifted signature is well explained by neutral hydrogen atoms escaping at rates of about 1e10 g s-1 that are blown away from the star by its strong radiation pressure and are quickly photoionized, resulting in a smaller exosphere than that of the warm Neptune GJ 436b. The stronger escape from GJ 3470b, however, may have led to the loss of about 4-35% of its current mass over its 2 Gyr lifetime.

One of the most important problems in the context of cataclysmic variables (CVs) is the lack of observations of systems with periods between 2 and 3.12 hours, known as the period gap. The orbital evolution of CVs with periods shorter than those in the gap is dominated by gravitational radiation while for periods exceeding those of the gap it is dominated by magnetic braking of the secondary star. Spruit & Ritter (1983) showed that as periods approach 3 hours and secondary stars become fully convective a sharp decline in magnetic dynamo and braking efficiency would result in such a gap. Recent X-ray observations finding coronal magnetic energy dissipation is similar in fully convective and partly radiative M dwarfs cast this theory into doubt. In this work, we use Zeeman-Doppler imaging observations culled from the literature to show that the complexity of the surface magnetic fields of rapidly rotating M dwarfs increases with decreasing rotation period. Garraffo et al. (2018) have shown that the efficiency of angular momentum loss of cool stars declines strongly with increasing complexity of their surface magnetic field. We explore the idea of Taam & Spruit (1989) that magnetic complexity might then explain the period gap. By generating synthetic CV populations using a schematic CV evolutionary approach, we show that the CV period gap can naturally arise as a consequence of a rise in secondary star magnetic complexity near the long period edge of the gap that renders a sharp decline in their angular momentum loss rate.

Radio observations grant access to a wide range of physical processes through different emission mechanisms. These processes range from thermal and quiescent to eruptive phenomena, such as shock waves and particle beams. We present a new synthetic radio imaging tool that calculates and visualizes the Bremsstrahlung radio emission. This tool works concurrently with state-of-the-art Magnetohydrodynamic (MHD) simulations of the solar corona using the code BATS-R-US. Our model produces results that are in good agreement with both high and low frequency observations of the solar disk. In this study, a ray-tracing algorithm is used and the radio intensity is computed along the actual curved ray trajectories. We illustrate the importance of refraction in locating the radio emitting source by comparison of the radio imaging illustrations when the line-of-sight instead of the refracted paths are considered. We are planning to incorporate non-thermal radio emission mechanisms in a future version of the radio imaging tool.

The search for exoplanets in the radio bands has been focused on detecting radio emissions produced by the interaction between magnetized planets and the stellar wind (auroral emission). Here we introduce a new tool, which is part of our MHD stellar corona model, to predict the ambient coronal radio emission and its modulations induced by a close planet. For simplicity, the present work assumes that the exoplanet is stationary in the frame rotating with the stellar rotation. We explore the radio flux modulations using a limited parameter space of idealized cases by changing the magnitude of the planetary field, its polarity, the planetary orbital separation, and the strength of the stellar field. We find that the modulations induced by the planet could be significant and observable in the case of hot Jupiter planets --- above 100% modulation with respect to the ambient flux in the $10-100~MHz$ range in some cases, and 2-10% in the frequency bands above $250~MHz$ for some cases. Thus, our work indicates that radio signature of exoplanets might not be limited to low-frequency radio range. We find that the intensity modulations are sensitive to the planetary magnetic field polarity for short-orbit planets, and to the stellar magnetic field strength for all cases. The new radio tool, when applied to real systems, could provide predictions for the frequency range at which the modulations can be observed by current facilities.

Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters.

We present results from a set of numerical simulations aimed at exploring the mechanism of coronal mass ejection (CME) suppression in active stars by an overlying large-scale magnetic field. We use a state-of-the-art 3D magnetohydrodynamic (MHD) code which considers a self-consistent coupling between an Alfv\'en wave-driven stellar wind solution, and a first-principles CME model based on the eruption of a flux-rope anchored to a mixed polarity region. By replicating the driving conditions used in simulations of strong solar CMEs, we show that a large-scale dipolar magnetic field of $75$ G is able to fully confine eruptions within the stellar corona. Our simulations also consider CMEs exceeding the magnetic energy used in solar studies, which are able to escape the large-scale magnetic field confinement. The analysis includes a qualitative and quantitative description of the simulated CMEs and their dynamics, which reveals a drastic reduction of the radial speed caused by the overlying magnetic field. With the aid of recent observational studies, we place our numerical results in the context of solar and stellar flaring events. In this way, we find that this particular large-scale magnetic field configuration establishes a suppression threshold around $\sim$$3 \times 10^{32}$~erg in the CME kinetic energy. Extending the solar flare-CME relations to other stars, such CME kinetic energies could be typically achieved during erupting flaring events with total energies larger than $6 \times 10^{32}$ erg (GOES class $\sim$X70).

We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospheres) depends on the relation between the conductances and impedances of the two regions. Since the energy flux of the stellar wind is very high at these planetary orbits, we find that for the case of high transmission efficiency (when the conductances and impedances are close in magnitude), the energy dissipation in the upper planetary atmospheres is also very large. On average, the Ohmic energy can reach $0.5-1~W/m^2$, about 1\% of the stellar irradiance and 5-15 times the EUV irradiance. Here, using constant values for the ionospheric conductance, we demonstrate that the stellar wind energy could potentially drive large atmospheric heating in terrestrial planets, as well as in hot jupiters. More detailed calculations are needed to assess the ionospheric conductance and to determine more accurately the amount of heating the stellar wind can drive in close-orbit planets.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

We use magnetograms of 8 solar analogues of ages 30~Myr to 3.6~Gyr obtained from Zeeman Doppler Imaging (ZDI) and taken from the literature, together with two solar magnetograms, to drive magnetohydrodynamical (MHD) wind simulations and construct an evolutionary scenario of the solar wind environment and its angular momentum loss rate. With observed magnetograms of the radial field strength as the only variant in the wind model, we find that power law model fitted to the derived angular momentum loss rate against time, $t$, results in a spin down relation $\Omega\propto t^{-0.51}$, for angular speed $\Omega$, which is remarkably consistent with the well-established Skumanich law $\Omega\propto t^{-0.5}$. We use the model wind conditions to estimate the magnetospheric standoff distances for an Earth-like test planet situated at 1~AU for each of the stellar cases, and to obtain trends of minimum and maximum wind ram pressure and average ram pressure in the solar system through time. The wind ram pressure declines with time as $\overline{P_{ram}}\propto t^{2/3}$, amounting to a factor of 50 or so over the present lifetime of the solar system.

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 evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares or by travelling shock waves and advected by the stellar wind. We have performed test-particle numerical simulations of energetic protons propagating into a realistic T~Tauri stellar wind, including a superposed small-scale magnetostatic turbulence. The isotropic (Kolmogorov power spectrum) turbulent component is synthesised along the individual particle trajectories. We have investigated the energy range $[0.1 - 10]$ GeV, consistent with expectations from {\it Chandra} X-ray observations of large flares on T~Tauri stars and with recent indications by the {\it Herschel} Space Observatory of a significant contribution of energetic particles to the disk ionization of young stars. In contrast with a previous theoretical study finding dominance of energetic particles over X-ray in the ionization throughout the disk, we find that the disk ionization is likely dominated by X-rays over much of its area except within narrow regions where particles are channeled onto the disk by the strongly-tangled and turbulent magnetic field. The radial thickness of such regions is $\sim 5$ stellar radii close to the star and broadens with increasing radial distance. This likely continues out to large distances from the star ($10$ AU or greater) where particles can be copiously advected and diffused by the turbulent wind.

We present a thermal emission spectrum of the bloated hot Jupiter HAT-P-32Ab from a single eclipse observation made in spatial scan mode with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope (HST). The spectrum covers the wavelength regime from 1.123 to 1.644 microns which is binned into 14 eclipse depths measured to an averaged precision of 104 parts-per million. The spectrum is unaffected by a dilution from the close M-dwarf companion HAT-P-32B, which was fully resolved. We complemented our spectrum with literature results and performed a comparative forward and retrieval analysis with the 1D radiative-convective ATMO model. Assuming solar abundance of the planet atmosphere, we find that the measured spectrum can best be explained by the spectrum of a blackbody isothermal atmosphere with Tp = 1995 +/- 17K, but can equally-well be described by a spectrum with modest thermal inversion. The retrieved spectrum suggests emission from VO at the WFC3 wavelengths and no evidence of the 1.4 micron water feature. The emission models with temperature profiles decreasing with height are rejected at a high confidence. An isothermal or inverted spectrum can imply a clear atmosphere with an absorber, a dusty cloud deck or a combination of both. We find that the planet can have continuum of values for the albedo and recirculation, ranging from high albedo and poor recirculation to low albedo and efficient recirculation. Optical spectroscopy of the planet's day-side or thermal emission phase curves can potentially resolve the current albedo with recirculation degeneracy.

We explore the scenario of a Coronal Mass Ejection (CME) being the cause of the observed continuous X-ray absorption of the August 30 1997 superflare on the eclipsing binary Algol (the Demon Star). The temporal decay of the absorption is consistent with absorption by a CME undergoing self-similar evolution with uniform expansion velocity. We investigate the kinematic and energetic properties of the CME using the ice-cream cone model for its three-dimensional structure in combination with the observed profile of the hydrogen column density decline with time. Different physically justified length scales were used that allowed us to estimate lower and upper limits of the possible CME characteristics. Further consideration of the maximum available magnetic energy in starspots leads us to quantify its mass as likely lying in the range $2\times 10^{21}$-$2\times 10^{22}$ g and kinetic energy in the range $7\times 10^{35}$-$3\times 10^{38}$ erg. The results are in reasonable agreement with extrapolated relations between flare X-ray fluence and CME mass and kinetic energy derived for solar CMEs.

A former member of the Hyades cluster, $\iota$ Horologii ($\iota$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $\iota$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the origin of cyclic activity and stellar magnetism in late-type stars. In this series of articles, we present the results of a comprehensive project aimed at studying the evolving magnetic field in this star and how this evolution influences its circumstellar environment. This paper summarizes the first stage of this investigation, with results from a long-term observing campaign of $\iota$ Hor using ground-based high-resolution spectropolarimetry. The analysis includes precise measurements of the magnetic activity and radial velocity of the star, and their multiple time-scales of variability. In combination with values reported in the literature, we show that the long-term chromospheric activity evolution of $\iota$ Hor follows a beating pattern, caused by the superposition of two periodic signals of similar amplitude at $P_{1} \simeq 1.97 \pm 0.02$ yr and $P_{2} \simeq 1.41 \pm 0.01$ yr. Additionally, using the most recent parameters for $\iota$ Hor b in combination with our activity and radial velocity measurements, we find that stellar activity dominates the radial velocity residuals, making the detection of additional planets in this system challenging. Finally, we report here the first measurements of the surface longitudinal magnetic field strength of $\iota$ Hor, which displays varying amplitudes within $\pm4$ G and served to estimate the rotation period of the star ($P_{\rm rot} =~7.70^{+0.18}_{-0.67}$ d).

Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between $10^3$ and $10^5$ times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfv\'enic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets.

We address the important question of whether the newly discovered exoplanet, Proxima Centauri b (PCb), is capable of retaining an atmosphere over long periods of time. This is done by adapting a sophisticated multi-species MHD model originally developed for Venus and Mars, and computing the ion escape losses from PCb. The results suggest that the ion escape rates are about two orders of magnitude higher than the terrestrial planets of our Solar system if PCb is unmagnetized. In contrast, if the planet does have an intrinsic dipole magnetic field, the rates are lowered for certain values of the stellar wind dynamic pressure, but they are still higher than the observed values for our Solar system's terrestrial planets. These results must be interpreted with due caution, since most of the relevant parameters for PCb remain partly or wholly unknown.

The development of the Zeeman-Doppler Imagine (ZDI) technique has provided synoptic observations of surface magnetic fields of low-mass stars. This led the Stellar Astrophysics community to adopt modeling techniques that have been used in solar physics using solar magnetograms. However, many of these techniques have been neglected by the solar community due to their failure to reproduce solar observations. Nevertheless, some of these techniques are still used to simulate the coronae and winds of solar analogs. Here we present a comparative study between two MHD models for the solar corona and solar wind. The first type of model is a polytropic wind model, and the second is the physics-based AWSOM model. We show that while the AWSOM model consistently reproduces many solar observations, the polytropic model fails to reproduce many of these observations and in the cases it does, its solutions are unphysical. Our recommendation is that polytropic models, which are used to estimate mass-loss rates and other parameters of solar analogs, must be first calibrated with solar observations. Alternatively, these models can be calibrated with models that capture more detailed physics of the solar corona (such as the AWSOM model), and that can reproduce solar observations in a consistent manner. Without such a calibration, the results of the polytropic models cannot be validated, but they can be wrongly used by others.

A planet orbiting in the "habitable zone" of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is "habitable". Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here we study the stellar wind conditions that Proxima b experiences over its orbit. We construct 3-D MHD models of the wind and magnetic field around Proxima Centauri using a surface magnetic field map for a star of the same spectral type and scaled to match the observed ~ 600 G surface magnetic field strength of Proxima. We examine the wind conditions and dynamic pressure over different plausible orbits that sample the constrained parameters of the orbit of Proxima b. For all the parameter space explored, the planet is subject to stellar wind pressures of more than 2000 times those experienced by Earth from the solar wind. During an orbit, Proxima b is also subject to pressure changes of 1 to 3 orders of magnitude on timescales of a day. Its magnetopause standoff distance consequently undergoes sudden and periodic changes by a factor of 2 to 5. Proxima b will traverse the interplanetary current sheet twice each orbit, and likely crosses into regions of subsonic wind quite frequently. These effects should be taken into account in any physically realistic assessment or prediction of its atmospheric reservoir, characteristics and loss.

Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several kilo Gauss large-scale magnetic fields at high latitudes, the corona and its X-ray emission are dominated by star-size large hot loops, while the smaller, underlying colder loops are not visible much in the X-ray. Based on this result we propose that, in rapidly rotating stars, emission from such coronal structures dominates the quiescent, cooler but saturated X-ray emission.

Flares we observe on stars in white light, UV or soft X-rays are probably harbingers of coronal mass ejections (CMEs). If we use the Sun as a guide, large stellar flares will dissipate two orders of magnitude less X-ray radiative energy than the kinetic energy in the associated CME. Since coronal emission on active stars appears to be dominated by flare activity, CMEs pose a quandary for understanding the fraction of their energy budget stars can spend on magnetic activity. One answer is magnetic suppression of CMEs, in which the strong large-scale fields of active stars entrap and prevent CMEs unless their free energy exceeds a critical value. The CME-less flaring active region NOAA 2192 presents a possible solar analogue of this. Monster CMEs will still exist, and have the potential to ravage planetary atmospheres.

We present the results of a comprehensive numerical simulation of the environment around three exoplanet-host stars (HD 1237, HD 22049, and HD 147513). Our simulations consider one of the latest models currently used for space weather studies in the Heliosphere. Large-scale magnetic field maps, recovered with two implementations of the tomographic technique of Zeeman-Doppler imaging, serve to drive steady-state solutions in each system. This paper contains the description of the stellar wind and inner astrosphere, while the coronal structure was previously discussed in Alvarado-G\'omez et al. (2016). The analysis includes the magneto-hydrodynamical properties of the stellar wind, the associated mass and angular momentum loss rates, as well as the topology of the astrospheric current sheet in each system. A systematic comparison among the considered cases is performed, including two reference solar simulations covering activity minimum and maximum. For HD 1237, we investigate the interactions between the structure of the developed stellar wind, and a possible magnetosphere around the Jupiter-mass planet in this system. We find that the process of particle injection into the planetary atmosphere is dominated by the density distribution rather than velocity profile of the stellar wind. In this context, we predict a maximum exoplanetary radio emission of 12 mJy at 40 MHz in this system, assuming the crossing of a high-density streamer during periastron passage. Furthermore, in combination with the analysis performed in Alvarado-G\'omez et al. (2016), we obtain for the first time a fully simulated mass loss-activity relation, which is compared and discussed in the context of the relation based on astrospheric detections proposed by Wood et al. (2005a). Finally, we provide a characterisation of the 3D properties of the stellar wind of these systems, at the inner edges of their habitable zones.

This study examines the relationship between magnetic field complexity and mass and angular momentum losses. Observations of open clusters have revealed a bimodal distribution of the rotation periods of solar-like stars that has proven difficult to explain under the existing rubric of magnetic braking. Recent studies suggest that magnetic complexity can play an important role in controlling stellar spin-down rates. However, magnetic morphology is still neglected in most rotation evolution models due to the difficulty of properly accounting for its effects on wind driving and angular momentum loss. Using state-of-the-art magnetohydrodynamical magnetized wind simulations we study the effect that different distributions of the magnetic flux at different levels of geometrical complexity have on mass and angular momentum loss rates. Angular momentum loss rates depend strongly on the level of complexity of the field but are independent of the way this complexity is distributed. We deduce the analytical terms representing the magnetic field morphology dependence of mass and angular momentum loss rates. We also define a parameter that best represents complexity for real stars. As a test, we use these analytical methods to estimate mass and angular momentum loss rates for $8$ stars with observed magnetograms and compare them to the simulated results. Magnetic field complexity provides a natural physical basis for stellar rotation evolution models requiring a rapid transition between weak and strong spin-down modes.

We present the results of a detailed numerical simulation of the circumstellar environment around three exoplanet-hosting stars. A state-of-the-art global magnetohydrodynamic (MHD) model is considered, including Alfv\'en wave dissipation as a self-consistent coronal heating mechanism. This paper contains the description of the numerical set-up, evaluation procedure, and the simulated coronal structure of each system (HD 1237, HD 22049 and HD 147513). The simulations are driven by surface magnetic field maps, recovered with the observational technique of Zeeman Doppler Imaging (ZDI). A detailed comparison of the simulations is performed, where two different implementations of this mapping routine are used to generate the surface field distributions. Quantitative and qualitative descriptions of the coronae of these systems are presented, including synthetic high-energy emission maps in the Extreme Ultra-Violet (EUV) and Soft X-rays (SXR) ranges. Using the simulation results, we are able to recover similar trends as in previous observational studies, including the relation between the magnetic flux and the coronal X-ray emission. Furthermore, for HD 1237 we estimate the rotational modulation of the high-energy emission due to the various coronal features developed in the simulation. We obtain variations, during a single stellar rotation cycle, up to 15\% for the EUV and SXR ranges. The results presented here will be used, in a follow-up paper, to self-consistently simulate the stellar winds and inner astrospheres of these systems.

COCOA-PUFS is an energy-diverse, time-domain study of the ultra-fast spinning, heavily spotted, yellow giant FK Com (HD117555; G4 III). This single star is thought to be a recent binary merger, and is exceptionally active by measure of its intense ultraviolet and X-ray emissions, and proclivity to flare. COCOA-PUFS was carried out with Hubble Space Telescope in the UV (120-300 nm), using mainly its high-performance Cosmic Origins Spectrograph, but also high-precision Space Telescope Imaging Spectrograph; Chandra X-ray Observatory in the soft X-rays (0.5-10 keV), utilizing its High-Energy Transmission Grating Spectrometer; together with supporting photometry and spectropolarimetry in the visible from the ground. This is an introductory report on the project. FK Com displayed variability on a wide range of time scales, over all wavelengths, during the week-long main campaign, including a large X-ray flare; "super-rotational broadening" of the far-ultraviolet "hot-lines" (e.g., Si IV 139 nm (T~80,000 K) together with chromospheric Mg II 280 nm and C II 133 nm (10,000-30,000 K); large Doppler swings suggestive of bright regions alternately on advancing and retreating limbs of the star; and substantial redshifts of the epoch-average emission profiles. These behaviors paint a picture of a highly extended, dynamic, hot (10 MK) coronal magnetosphere around the star, threaded by cooler structures perhaps analogous to solar prominences, and replenished continually by surface activity and flares. Suppression of angular momentum loss by the confining magnetosphere could temporarily postpone the inevitable stellar spindown, thereby lengthening this highly volatile stage of coronal evolution.

The results from a spectro-polarimetric study of the planet-hosting Sun-like star, HD 147513 (G5V), are presented here. Robust detections of Zeeman signatures at all observed epochs indicate a surface magnetic field, with longitudinal magnetic field strengths varying between 1.0-3.2 G. Radial velocity variations from night to night modulate on a similar timescale to the longitudinal magnetic field measurements. These variations are therefore likely due to the rotational modulation of stellar active regions rather than the much longer timescale of the planetary orbit (Porb=528 d). Both the longitudinal magnetic field measurements and radial velocity variations are consistent with a rotation period of 10 +/- 2 days, which are also consistent with the measured chromospheric activity level of the star (log R'(HK)=-4.64). Together, these quantities indicate a low inclination angle, i~18 degrees. We present preliminary magnetic field maps of the star based on the above period and find a simple poloidal large-scale field. Chemical analyses of the star have revealed that it is likely to have undergone a barium-enrichment phase in its evolution because of a higher mass companion. Despite this, our study reveals that the star has a fairly typical activity level for its rotation period and spectral type. Future studies will enable us to explore the long-term evolution of the field, as well as to measure the stellar rotation period, with greater accuracy.

Rotation evolution of late-type stars is dominated by magnetic braking and the underlying factors that control this angular momentum loss are important for the study of stellar spin-down. In this work, we study angular momentum loss as a function of two different aspects of magnetic activity using a calibrated Alfv\'en wave-driven magnetohydrodynamic wind model: the strengths of magnetic spots and their distribution in latitude. By driving the model using solar and modified solar surface magnetograms, we show that the topology of the field arising from the net interaction of both small-scale and large-scale field is important for spin-down rates and that angular momentum loss is not a simple function of large scale magnetic field strength. We find that changing the latitude of magnetic spots can modify mass and angular momentum loss rates by a factor of two. The general effect that causes these differences is the closing down of large-scale open field at mid- and high-latitudes by the addition of the small-scale field. These effects might give rise to modulation of mass and angular momentum loss through stellar cycles, and present a problem for ab initio attempts to predict stellar spin-down based on wind models. For all the magnetogram cases considered here, from dipoles to various spotted distributions, we find that angular momentum loss is dominated by the mass loss at mid-latitudes. The spin-down torque applied by magnetized winds therefore acts at specific latitudes and is not evenly distributed over the stellar surface, though this aspect is unlikely to be important for understanding spin-down and surface flows on stars.

Stars orbiting within 1$\arcsec$ of the supermassive black hole in the Galactic Centre, Sgr A*, are notoriously difficult to detect due to obscuration by gas and dust. We show that some stars orbiting this region may be detectable via synchrotron emission. In such instances, a bow shock forms around the star and accelerates the electrons. We calculate that around the 10 GHz band (radio) and at 10$^{14}$ Hz (infrared) the luminosity of a star orbiting the black hole is comparable to the luminosity of Sgr A*. The strength of the synchrotron emission depends on a number of factors including the star's orbital velocity. Thus, the ideal time to observe the synchrotron flux is when the star is at pericenter. The star S2 will be $\sim 0.015\arcsec$ from Sgr A* in 2018, and is an excellent target to test our predictions.

We analyse the magnetic activity characteristics of the planet hosting Sun-like star, HD 1237, using HARPS spectro-polarimetric time-series data. We find evidence of rotational modulation of the magnetic longitudinal field measurements consistent with our ZDI analysis, with a period of 7 days. We investigate the effect of customising the LSD mask to the line depths of the observed spectrum and find that it has a minimal effect on shape of the extracted Stokes V profile but does result in a small increase in the S/N ($\sim$ 7%). We find that using a Milne-Eddington solution to describe the local line profile provides a better fit to the LSD profiles in this slowly rotating star, which also impacts the recovered ZDI field distribution. We also introduce a fit-stopping criterion based on the information content (entropy) of the ZDI maps solution set. The recovered magnetic field maps show a strong (+90 G) ring-like azimuthal field distribution and a complex radial field dominating at mid latitudes ($\sim$45 degrees). Similar magnetic field maps are recovered from data acquired five months apart. Future work will investigate how this surface magnetic field distribution impacts the coronal magnetic field and extended environment around this planet-hosting star.

The potential field approximation has been providing a fast, and computationally inexpensive estimation for the solar corona's global magnetic field geometry for several decades. In contrast, more physics-based global magnetohydrodynamic (MHD) models have been used for a similar purpose, while being much more computationally expensive. Here, we investigate the difference in the field geometry between a global MHD model and the potential field source surface model (PFSSM) by tracing individual magnetic field lines in the MHD model from the Alfven surface (AS), through the source surface (SS), all the way to the field line footpoint, and then back to the source surface in the PFSSM. We also compare the flux-tube expansion at two points at the SS and the AS along the same radial line. We study the effect of solar cycle variations, the order of the potential field harmonic expansion, and different magnetogram sources. We find that the flux-tube expansion factor is consistently smaller at the AS than at the SS for solar minimum and the fast solar wind, but it is consistently larger for solar maximum and the slow solar wind. We use the Wang--Sheeley--Arge (WSA) model to calculate the associated wind speed for each field line, and propagate these solar-wind speeds to 1AU. We find a more than five hours deviation in the arrival time between the two models for 20% of the field lines in the solar minimum case, and for 40% of the field lines in the solar maximum case.

Observations of young open clusters have revealed a bimodal distribution of fast and slower rotation rates that has proven difficult to explain with predictive models of spin down that depend on rotation rates alone. The Metastable Dynamo Model proposed recently by Brown, employing a stochastic transition probability from slow to more rapid spin down regimes, appears to be more successful but lacks a physical basis for such duality. Using detailed 3D MHD wind models computed for idealized multipole magnetic fields, we show that surface magnetic field complexity can provide this basis. Both mass and angular momentum losses decline sharply with increasing field complexity. Combined with observation evidence for complex field morphologies in magnetically active stars, our results support a picture in which young, rapid rotators lose angular momentum in an inefficient way because of field complexity. During this slow spin-down phase, magnetic complexity is eroded, precipitating a rapid transition from weak to strong wind coupling.

We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary distance from the star decreases. Our simulations reveal different topologies of the planetary space environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy deposition in to the atmosphere during the transition along the planetary orbit. We find that the stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a significant amount of the planetary atmosphere could be eroded over the course of a billion years.

Stellar flares, winds and coronal mass ejections form the space weather. They are signatures of the magnetic activity of cool stars and, since activity varies with age, mass and rotation, the space weather that extra-solar planets experience can be very different from the one encountered by the solar system planets. How do stellar activity and magnetism influence the space weather of exoplanets orbiting main-sequence stars? How do the environments surrounding exoplanets differ from those around the planets in our own solar system? How can the detailed knowledge acquired by the solar system community be applied in exoplanetary systems? How does space weather affect habitability? These were questions that were addressed in the splinter session "Cool stars and Space Weather", that took place on 9 Jun 2014, during the Cool Stars 18 meeting. In this paper, we present a summary of the contributions made to this session.

Angular Momentum Loss is important for understanding astrophysical phenomena such as stellar rotation, magnetic activity, close binaries, and cataclysmic variables. Magnetic breaking is the dominant mechanism in the spin down of young late-type stars. We have studied angular momentum loss as a function of stellar magnetic activity. We argue that the complexity of the field and its latitudinal distribution are crucial for angular momentum loss rates. In this work we discuss how angular momentum is modulated by magnetic cycles, and how stellar spin down is not just a simple function of large scale magnetic field strength.

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic (MHD) models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvenic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvenic sectors, while no bow shock forms in the sub-Alfvenic sectors. The planets reside most of the time in the sub-Alfvenic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the planetary interaction with the stellar wind. For the steady-state solution, the heating is about 0.1-3\% of the total incoming stellar irradiation, and it is enhanced by 50\% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport.

Abridged. Here we report on the X-ray activity of the primary star, HD189733 A, using a new XMM-Newton observation and a comparison with the previous X-ray observations. The spectrum in the quiescent intervals is described by two temperatures at 0.2 keV and 0.7 keV, while during the flares a third component at 0.9 keV is detected. We obtain estimates of the electron density in the range $n_e = 1.6 - 13 \times 10^{10}$ cm$^{-3}$ and thus the corona of HD189733 A appears denser than the solar one. {For the third time, we observe a large flare that occurred just after the eclipse of the planet. Together with the flares observed in 2009 and 2011, the events are restricted to a small planetary phase range of $\phi = 0.55-0.65$. Although we do not find conclusive evidence of a significant excess of flares after the secondary transits, we suggest that the planet might trigger such flares when it passes close to locally high magnetic field of the underlying star at particular combinations of stellar rotational phases and orbital planetary phases. For the most recent flares, a wavelet analysis of the light curve suggests a loop of length of four stellar radii at the location of the bright flare, and a local magnetic field of order of 40-100 G, in agreement with the global field measured in other studies. The loop size suggests an interaction of magnetic nature between planet and star, separated by only $\sim8 R_*$. We also detect the stellar companion (HD 189733 B, $\sim12"$ from the primary star) in this XMM observation. Its very low X-ray luminosity ($L_X = 3.4\times 10^{26}$ erg s$^{-1}$) confirms the old age of this star and of the binary system. The high activity of the primary star is best explained by a transfer of angular momentum from the planet to the star.

Stellar winds are believed to be the dominant factor in spin down of stars over time. However, stellar winds of solar analogs are poorly constrained due to the challenges in observing them. A great improvement has been made in the last decade in our understanding of the mechanisms responsible for the acceleration of the solar wind and in the development of numerical models for solar and stellar winds. In this paper, we present a grid of Magnetohydrodynamic (MHD) models to study and quantify the values of stellar mass-loss and angular momentum loss rates as a function of the stellar rotation period, magnetic dipole component, and coronal base density. We derive simple scaling laws for the loss rates as a function of these parameters, and constrain the possible mass-loss rate of stars with thermally-driven winds. Despite the success of our scaling law in matching the results of the model, we find a deviation between the "solar dipole" case and a real case based on solar observations that overestimates the actual solar mass-loss rate by a factor of 3. This implies that the model for stellar fields might require a further investigation with higher complexity which might include the use of a filling factor for active regions, as well as the distribution of the strength of the small-scale fields. Mass loss rates in general are largely controlled by the magnetic field strength, with the wind density varying in proportion to the confining magnetic pressure $B^2$. We also find that the mass-loss rates obtained using our grid models drop much faster with the increase in rotation period than scaling laws derived using observed stellar activity. For main-sequence solar-like stars, our scaling law for angular momentum loss vs. poloidal magnetic field strength retrieves the well-known Skumanich decline of angular velocity with time.

Analysis of a database of solar coronal mass ejections (CMEs) and associated flares over the period 1996-2007 finds well-behaved power law relationships between the 1-8 AA flare X-ray fluence and CME mass and kinetic energy. We extrapolate these relationships to lower and higher flare energies to estimate the mass and energy loss due to CMEs from stellar coronae, assuming that the observed X-ray emission of the latter is dominated by flares with a frequency as a function of energy dn/dE=kE^-alpha. For solar-like stars at saturated levels of X-ray activity, the implied losses depend fairly weakly on the assumed value of alpha and are very large: M_dot ~ 5x10^-10 M_sun/yr and E_dot ~ 0.1L_sun. In order to avoid such large energy requirements, either the relationships between CME mass and speed and flare energy must flatten for X-ray fluence >~ 10^31 erg, or the flare-CME association must drop significantly below 1 for more energetic events. If active coronae are dominated by flares, then the total coronal energy budget is likely to be up to an order of magnitude larger than the canonical 10^-3 L_bol X-ray saturation threshold. This raises the question of what is the maximum energy a magnetic dynamo can extract from a star? For an energy budget of 1% of L_bol, the CME mass loss rate is about 5x10^-11 M_sun/yr.

We study the influence of the spatial resolution on scales of $5\deg$ and smaller of solar surface magnetic field maps on global magnetohydrodynamic solar wind models, and on a model of coronal heating and X-ray emission. We compare the solutions driven by a low-resolution Wilcox Solar Observatory magnetic map, the same map with spatial resolution artificially increased by a refinement algorithm, and a high-resolution Solar and Heliospheric Observatory Michelson Doppler Imager map. We find that both the wind structure and the X-ray morphology are affected by the fine-scale surface magnetic structure. Moreover, the X-ray morphology is dominated by the closed loop structure between mixed polarities on smaller scales and shows significant changes between high and low resolution maps. We conclude that three-dimensional modeling of coronal X-ray emission has greater surface magnetic field spatial resolution requirements than wind modeling, and can be unreliable unless the dominant mixed polarity magnetic flux is properly resolved.

We employ three-dimensional state of the art magnetohydrodynamic models of the early solar wind and heliosphere and a two-dimensional model for cosmic ray transport to investigate the cosmic ray spectrum and flux near the Archean Earth. We assess how sensitive the cosmic ray spectrum is to changes in the sunspot placement and magnetic field strength, the large scale dipole magnetic field strength, the wind ram pressure, and the Sun's rotation period. Overall, our results confirm earlier work that suggested the Archean Earth would have experienced a greatly reduced cosmic ray flux than is the case today. The cosmic ray reduction for the early Sun is mainly due to the shorter solar rotation period and tighter winding of the Parker spiral, and to the different surface distribution of the more active solar magnetic field. These effects lead to a global reduction of the cosmic ray flux at 1AU by up to two orders of magnitude or more. Variations in the sunspot magnetic field have more effect on the flux than variations in the dipole field component. The wind ram pressure affects the cosmic ray flux through its influence on the size of the heliosphere via the pressure balance with the ambient interstellar medium. Variations in the interstellar medium pressure experienced by the solar system in orbit through Galaxy could lead to order of magnitude changes in the cosmic ray flux at Earth on timescales of a few million years.

Atmospheric mass-loss from Hot-jupiters can be large due to the close proximity of these planets to their host star and the strong radiation the planetary atmosphere receives. On Earth, a major contribution to the acceleration of atmospheric ions comes from the vertical separation of ions and electrons, and the generation of the ambipolar electric field. This process, known as the "polar wind", is responsible for the transport of ionospheric constituents to the Earth's magnetosphere, where they are well observed. The polar wind can also be enhanced by a relatively small fraction of super-thermal electrons (photoelectrons) generated by photoionization. We formulate a simplified calculation of the effect of the ambipolar electric field and the photoelectrons on the ion scale-height in a generalized manner. We find that the ion scale-height can be increased by a factor of 2-15 due to the polar wind effects. We also estimate a lower limit of an order of magnitude increase of the ion density and the atmospheric mass-loss rate when polar wind effects are included.

We present a novel method for Ankylography: three-dimensional structure reconstruction from a single shot diffraction intensity pattern. Our approach allows reconstruction of objects containing many more details than was ever demonstrated, in a faster and more accurate fashion

We study the effect of a white dwarf on the spin-down of a cataclysmic variable system using a three-dimensional magnetohydrodynamic numerical model. The model includes the stellar corona, the stellar wind, and the WD mass and magnetic field. The existence of the WD modifies the system spin-down by physically blocking the stellar wind, restructuring the wind, channeling the wind towards the WD surface, and by modifying the shape and size of the Alfv\'en surface. The combination of these processes differs among a set of simple test cases, and the resulting angular momentum loss rates vary by factors of 2-3, and by factors of two relative to a test model with a single M dwarf. While the model employs some simplifications, the results suggest angular momentum loss schemes currently employed in cataclysmic variable studies do not require drastic revision. Insights are also gained on wind accretion. We find that efficient accretion switches on quite rapidly with decreasing orbital separation. Accretion rates depend on magnetic field alignment and should be modulated by magnetic cycles on the M dwarf. For particular values of white dwarf magnetic field strength, an efficient syphoning of coronal plasma from the inward facing M dwarf hemisphere occurs. Wind accretion rates are expected to vary by factors of 10 or more between fairly similar close binaries, depending on magnetic field strengths and orbital separation.

We report on the follow-up XMM-Newton observation of the planet-hosting star HD 189733 we obtained in April 2011. We observe a flare just after the secondary transit of the hot Jupiter. This event shares the same phase and many of the characteristics of the flare we observed in 2009. We suggest that a systematic interaction between planet and stellar magnetic fields when the planet passes close to active regions on the star can lead to periodic variability phased with planetary motion. By mean of high resolution X-ray spectroscopy with RGS we determine that the corona of this star is unusually dense.

Stellar mass loss rates are an important input ingredient for stellar evolution models since they determine stellar evolution parameters such as stellar spin-down and increase in stellar luminosity through the lifetime of a star. Due to the lack of direct observations of stellar winds from Sun-like stars stellar X-ray luminosity and stellar level of activity are commonly used as a proxy for estimating stellar mass loss rates. However, such an intuitive activity --- mass loss rate relation is not well defined. In this paper, I study the mass loss rate of the Sun as a function of its activity level. I compare in-situ solar wind measurements with the solar activity level represented by the solar X-ray flux. I find no clear dependency of the solar mass flux on solar X-ray flux. Instead, the solar mass loss rate is scattered around an average value of $2\cdot10^{-14}\;M_\odot\;yr^{-1}$. This independency of the mass loss rate on level of activity can be explained by the fact that the activity level is governed by the large modulations in the solar close magnetic flux, while the mass loss rate is governed by the rather constant open magnetic flux. I derive a simple expression for stellar mass loss rates as a function of the stellar ambient weak magnetic field, the stellar radius, the stellar escape velocity, and the average height of the Alfv\'en surface. This expression predicts stellar mass loss rates of $10^{-15}-10^{-12}\;M_\odot\;yr^{-1}$ for Sun-like stars.

We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to the direction of propagation of the CME. The combination of this feature and the fact that the CME does not expand much by the time it reaches the planet leads to a unique CME-magnetosphere interaction, where the CME itself is highly affected by the presence of the planetary magnetosphere. We find that the planet is well-shielded from CME penetration, even for a relatively weak internal magnetic field. The planetary angular momentum loss associated with such an event is negligible compared to the total planetary angular momentum. We also find that the energy which is deposited in the magnetosphere is much higher than in the case of the Earth, and our simulation suggests there is a large-scale change in the orientation of the magnetosphere-ionosphere current system during the CME event.

We present a focused parameter study of solar wind - magnetosphere interaction for the young Sun and Earth, $~3.5$ Ga ago, that relies on magnetohydrodynamic (MHD) simulations for both the solar wind and the magnetosphere. By simulating the quiescent young Sun and its wind we are able to propagate the MHD simulations up to Earth's magnetosphere and obtain a physically realistic solar forcing of it. We assess how sensitive the young solar wind is to changes in the coronal base density, sunspot placement and magnetic field strength, dipole magnetic field strength and the Sun's rotation period. From this analysis we obtain a range of plausible solar wind conditions the paleomagnetosphere may have been subject to. Scaling relationships from the literature suggest that a young Sun would have had a mass flux different from the present Sun. We evaluate how the mass flux changes with the aforementioned factors and determine the importance of this and several other key solar and magnetospheric variables with respect to their impact on the paleomagnetosphere. We vary the solar wind speed, density, interplanetary magnetic field strength and orientation as well as Earth's dipole magnetic field strength and tilt in a number of steady-state scenarios that are representative of young Sun-Earth interaction. This study is done as a first step of a more comprehensive effort towards understanding the implications of Sun-Earth interaction for planetary atmospheric evolution.

We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sectors with different plasma properties. The star-planet interaction varies in magnitude and complexity, depending on the planetary phase, planetary magnetic field strength, and the relative orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient energy to produce hot-spots observed in Ca II lines in giant planet hosting stars. However, the short duration of the reconnection event suggests that such SPI cannot be observed persistently.

We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfv\'en surface. This reduction in angular momentum loss is due to the interaction between the stellar and planetary Alfv\'en surfaces, which modifies the global structure of the stellar corona and stellar wind on the hemisphere facing the planet, as well as the opposite hemisphere. The simulation also shows that the magnitude of change in angular momentum loss rate depends mostly on the strength of the planetary magnetic field and not on its polarity. The interaction however, begins at greater separation if the overall field topology of the star and the planet are of anti-aligned. Our results are consistent with evidence for excess angular momentum in stars harboring close-in giant planets, and show that the reduction in wind-driven angular momentum loss can compete with, and perhaps dominate, spin-up due to tidal interaction.

We report on two XMM-Newton observations of the planetary host star HD189733. The system has a close in planet and it can potentially affect the coronal structure via interactions with the magnetosphere. We have obtained X-ray spectra and light curves from EPIC and RGS on board XMM-Newton which we have analyzed and interpreted. We reduced X-ray data from primary transit and secondary eclipse occurred in April 17th 2007 and May 18th 2009, respectively. In the April 2007 observation only variability due to weak flares is recognized. In 2009 HD189733 exhibited a X-ray flux always larger than in the 2007 observation. The average flux in 2009 was higher than in 2007 observation by a factor of 45%. During the 2009 secondary eclipse we observed a softening of the X-ray spectrum significant at level of ~3 sigma. Further, we observed the most intense flare recorded at either epochs. This flare occurred 3 ks after the end of the eclipse.The flare decay shows several minor ignitions perhaps linked to the main event and hinting for secondary loops that emit triggered by the main loop. Magneto-Hydro-Dynamical (MHD) simulations show that the magnetic interaction between planet and star enhances the density and the magnetic field in a region comprised between the planet and the star because of their relative orbital/rotation motion. X-ray observations and model predictions are globally found in agreement, despite the quite simple MHD model and the lack of precise estimate of parameters including the alignment and the intensity of stellar and planetary magnetic fields. Future observations should confirm or disprove this hypothesis, by determining whether flares are systematically recurring in the light curve at the same planetary phase.

We perform a numerical simulation of the corona of the young, rapidly rotating K0 dwarf AB Doradus using a global MHD model. The model is driven by a surface map of the radial magnetic field constructed using Zeeman-Doppler Imaging. We find that the global structure of the stellar corona is dominated by strong azimuthal tangling of the magnetic field due to the rapid rotation. The MHD solution enables us to calculate realistic Alfv\'en surfaces and we can therefore estimate the stellar mass loss rate and angular momentum loss rate without making undue theoretical simplifications. We consider three cases, parametrized by the base density of the corona, that span the range of possible solutions for the system. We find that overall, the mass and angular-momentum loss rates are higher than in the solar case; the mass loss rates are 10 to 500 times higher, and the angular momentum loss rate can be up to $3\times{10}^4$ higher than present day solar values. Our simulations show that this model can be use to constrain the wide parameter space of stellar systems. It also shows that an MHD approach can provide more information about the physical system over the commonly used potential field extrapolation.

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance ($\approx 1R_\odot$) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

We present a three-dimensional simulation of the corona of an FK Com-type rapidly rotating G giant using a magnetohydrodynamic model that was originally developed for the solar corona in order to capture the more realistic, non-potential coronal structure. We drive the simulation with surface maps for the radial magnetic field obtained from a stellar dynamo model of the FK Com system. This enables us to obtain the coronal structure for different field topologies representing different periods of time. We find that the corona of such an FK Com-like star, including the large scale coronal loops, is dominated by a strong toroidal component of the magnetic field. This is a result of part of the field being dragged by the radial outflow, while the other part remains attached to the rapidly rotating stellar surface. This tangling of the magnetic field,in addition to a reduction in the radial flow component, leads to a flattening of the gas density profile with distance in the inner part of the corona. The three-dimensional simulation provides a global view of the coronal structure. Some aspects of the results, such as the toroidal wrapping of the magnetic field, should also be applicable to coronae on fast rotators in general, which our study shows can be considerably different from the well-studied and well-observed solar corona. Studying the global structure of such coronae should also lead to a better understanding of their related stellar processes, such as flares and coronal mass ejections, and in particular, should lead to an improved understanding of mass and angular momentum loss from such systems.

Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities,suggest that these close-in planets significantly affect the structure of the outer atmosphere of the star through interactions between the stellar magnetic field and the planetary magnetosphere. Here we carry out the first detailed three-dimensional MagnetoHydroHynamics (MHD) simulation containing the two magnetic bodies and explore the consequences of such interactions on the steady-state coronal structure. The simulations reproduce the observable features of 1) increase in the total X-ray luminosity, 2) appearance of coronal hot spots, and 3) phase shift of these spots with respect to the direction of the planet. The proximate cause of these is an increase in the density of coronal plasma in the direction of the planet, which prevents the corona from expanding and leaking away this plasma via a stellar wind. The simulations produce significant low temperature heating. By including dynamical effects, such as the planetary orbital motion, the simulation should better reproduce the observed coronal heating.

On 13 February 2009, a coronal wave -- CME -- dimming event was observed in quadrature by the STEREO spacecraft. We analyze this event using a three-dimensional, global magnetohydrodynamic (MHD) model for the solar corona. The numerical simulation is driven and constrained by the observations, and indicates where magnetic reconnection occurs between the expanding CME core and surrounding environment. We focus primarily on the lower corona, extending out to $3R_{\odot}$; this range allows simultaneous comparison with both EUVI and COR1 data. Our simulation produces a diffuse coronal bright front remarkably similar to that observed by STEREO/EUVI at 195 \AA. It is made up of \emph{two} components, and is the result of a combination of both wave and non-wave mechanisms. The CME becomes large-scale quite low ($<$ 200 Mm) in the corona. It is not, however, an inherently large-scale event; rather, the expansion is facilitated by magnetic reconnection between the expanding CME core and the surrounding magnetic environment. In support of this, we also find numerous secondary dimmings, many far from the initial CME source region. Relating such dimmings to reconnecting field lines within the simulation provides further evidence that CME expansion leads to the "opening" of coronal field lines on a global scale. Throughout the CME expansion, the coronal wave maps directly to the CME footprint. Our results suggest that the ongoing debate over the "true" nature of diffuse coronal waves may be mischaracterized. It appears that \emph{both} wave and non-wave models are required to explain the observations and understand the complex nature of these events.

We simulate the effect of latitudinal variations in the location of star spots, as well as their magnetic field strength, on stellar angular momentum loss to the stellar wind. We use the Michigan solar corona global MagnetoHydroDynamic model, which incorporates realistic relation between the magnetic field topology and the wind distribution. We find that the spots location significantly affects the stellar wind structure, and as a result, the total mass loss rate and angular momentum loss rate. In particular, we find that the angular momentum loss rate is controlled by the mass flux when spots are located at low latitudes but is controlled by an increased plasma density between the stellar surface and the Alfven surface when spots are located at high latitudes. Our results suggest that there might be a feedback mechanism between the magnetic field distribution, wind distribution, angular momentum loss through the wind, and the motions at the convection zone that generate the magnetic field. This feedback might explain the role of coronal magnetic fields in stellar dynamos.

This study examines the tail disconnection event on April 20, 2007 on comet 2P/Encke, caused by a coronal mass ejection (CME) at a heliocentric distance of 0.34 AU. During their interaction, both the CME and the comet are visible with high temporal and spatial resolution by the STEREO-A spacecraft. Previously, only current sheets or shocks have been accepted as possible reasons for comet tail disconnections, so it is puzzling that the CME caused this event. The MHD simulation presented in this work reproduces the interaction process and demonstrates how the CME triggered a tail disconnection in the April 20 event. It is found that the CME disturbs the comet with a combination of a $180^\circ$ sudden rotation of the interplanetary magnetic field (IMF), followed by a $90^\circ$ gradual rotation. Such an interpretation applies our understanding of solar wind-comet interactions to determine the \textit{in situ} IMF orientation of the CME encountering Encke.