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.

Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either alone or in combination, more work is required to distinguish between possible solutions. Here we investigate nuclear reaction rate and opacity modifications on Cepheid evolution using the MESA code. We discuss the effects of opacity increases at envelope temperatures of 200,000-400,000 K proposed to explain the pulsation properties of hybrid main-sequence beta Cep/Slowly Pulsating B (SPB) variables which will evolve into Cepheids. We make use of the RSP nonlinear radial pulsation modeling capability in MESA to calculate periods and radial velocity amplitudes of Galactic Cepheids V1334 Cyg, Polaris, and delta Cep.

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.

V473 Lyr is a classical Cepheid which is unique in having substantial amplitude variations with a period of approximately 3.3 years, thought to be similar to the Blazhko variations in RR Lyrae stars. We obtained an {\it XMM-Newton} observation of this star to followup a previous detection in X-rays. Rather than the X-ray burst and rapid decline near maximum radius seen in $\delta$ Cephei itself, the X-ray flux in V473 Lyr remained constant for a third of the pulsation cycle covered by the observation. Thus the X-rays are most probably not produced by the changes around the pulsation cycle. The X-ray spectrum is soft (kT = 0.6 keV), with X-ray properties which are consistent with a young low mass companion. Previously there was no evidence of a companion in radial velocities or in {\it Gaia} and {\it Hipparcos} proper motions. While this rules out companions which are very close or very distant, a binary companion at a separation between 30 and 300 AU is possible. This is an example of an X-ray observation revealing evidence of a low mass companion, which is important in completing the mass ratio statistics of binary Cepheids. Furthermore, the detection of a young X-ray bright companion is a further indication that the Cepheid (primary) is a Population I star, even though its pulsation behavior differs from other classical Cepheids.

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.

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.

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.

Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years.

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.

X-ray observatories contribute fundamental advances in Solar System studies by probing Sun-object interactions, developing planet and satellite surface composition maps, probing global magnetospheric dynamics, and tracking astrochemical reactions. Despite these crucial results, the technological limitations of current X-ray instruments hinder the overall scope and impact for broader scientific application of X-ray observations both now and in the coming decade. Implementation of modern advances in X-ray optics will provide improvements in effective area, spatial resolution, and spectral resolution for future instruments. These improvements will usher in a truly transformative era of Solar System science through the study of X-ray emission.

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.

Radio observations of cool stellar systems provide unique information on their magnetic fields, high-energy processes, and chemistry. Buoyed by powerful new instruments (e.g. ALMA, JVLA, LOFAR), advances in related fields (e.g., the Gaia astrometric revolution), and above all a renewed interest in the relevant stellar astrophysics, stellar radio astronomy is experiencing a renaissance. In this splinter session, participants took stock of the present state of stellar radio astronomy to chart a course for the field's future.

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 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).

An exospheric kinetic solar wind model is interfaced with an observation-driven single fluid magnetohydrodynamic (MHD) model. Initially, a photospheric magnetogram serves as observational input in the fluid approach to extrapolate the heliospheric magnetic field. Then semi-empirical coronal models are used for estimating the plasma characteristics up to a heliocentric distance of 0.1AU. From there on a full MHD model which computes the three-dimensional time-dependent evolution of the solar wind macroscopic variables up to the orbit of the Earth is used. After interfacing the density and velocity at the inner MHD boundary, we compare with the results of a kinetic exospheric solar wind model based on the assumption of Maxwell and Kappa velocity distribution functions for protons and electrons respectively, as well as with \textit{in situ} observations at 1AU. This provides insight on more physically detailed processes, such as coronal heating and solar wind acceleration, that naturally arise by inclusion of suprathermal electrons in the model. We are interested in the profile of the solar wind speed and density at 1AU, in characterizing the slow and fast source regions of the wind and in comparing MHD with exospheric models in similar conditions. We calculate the energetics of both models from low to high heliocentric distances.

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 present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.

In this paper we present an update on the open source MPI-AMRVAC simulation toolkit where we focus on solar- and non-relativistic astrophysical magneto-fluid dynamics. We highlight recent developments in terms of physics modules such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated at the example of monotonicity preserving fifth order reconstruction. Strong stability preserving high order Runge-Kutta time steppers are used to obtain stable evolutions in multidimensional applications realizing up to fourth order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically, or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed, with an outlook to follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30 000 processors allowing to exploit modern petascale infrastructure.