Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfv\'en-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $\tau_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( \Omega_\ast / \Omega_\odot \right)^{2.82}$, yielding a spin-down law $\Omega_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfv\'en waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto \Omega_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfv\'en radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfv\'en radius matches previous works which used thermally-driven winds. Our results suggest that the Alfv\'en-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence.

In contrast with the fast solar wind, that originates in coronal holes, the source of the slow solar wind is still debated. Often intermittent and enriched with low FIP elements -- akin to what is observed in closed coronal loops -- the slow wind could form in bursty events nearby helmet streamers. Slow winds also exhibit density perturbations which have been shown to be periodic and could be associated with flux ropes ejected from the tip of helmet streamers, as shown recently by the WISPR white light imager onboard Parker Solar Probe (PSP). In this work, we propose that the main mechanism controlling the release of flux ropes is a flow-modified tearing mode at the heliospheric current sheet (HCS). We use MHD simulations of the solar wind and corona to reproduce realistic configurations and outflows surrounding the HCS. We find that this process is able to explain long ($\sim 10-20$h) and short ($\sim 1-2$h) timescales of density structures observed in the slow solar wind. This study also sheds new light on the structure, topology and composition of the slow solar wind, and could be, in the near future, compared with white light and in situ PSP observations.

The distribution of hot Jupiters, for which star-planet interactions can be significant, questions the evolution of exosystems. We aim to follow the orbital evolution of a planet along the rotational and structural evolution of the host star by taking into account the coupled effects of tidal and magnetic torques from ab initio prescriptions. It allows us to better understand the evolution of star-planet systems and to explain some properties of the distribution of observed close-in planets. To this end we use a numerical model of a coplanar circular star-planet system taking into account stellar structural changes, wind braking and star-planet interactions, called ESPEM (Benbakoura et al. 2019). We find that depending on the initial configuration of the system, magnetic effects can dominate tidal effects during the various phases of the evolution, leading to an important migration of the planet and to significant changes on the rotational evolution of the star. Both kinds of interactions thus have to be taken into account to predict the evolution of compact star-planet systems.

Context. Stellar spin-down is the result of a complex process involving rotation, dynamo, wind and magnetism. Multi-wavelength surveys of solar-like stars have revealed the likely existence of relationships between their rotation, X-ray luminosity, mass-losses and magnetism. Those impose strong constraints on the corona and wind of cool stars. Aims. We aim to provide power-law prescriptions of the mass-loss of stars, of their magnetic field, and of their base coronal density and temperature that are compatible with their observationally-constrained spin-down. Methods. We link the magnetic field and the mass-loss rate from a wind torque formulation in agreement with the distribution of stellar rotation periods in open clusters and the Skumanich law. Given a wind model and an expression of the X-ray luminosity from radiative losses, we constrain the coronal properties by assuming different physical scenarii linking closed loops to coronal holes. Results. We find that the magnetic field and the mass loss are involved in a one-to-one correspondence constrained from spin-down considerations. We show that a magnetic field depending on both the Rossby number and the stellar mass is required to keep a consistent spin-down model. The estimates of the magnetic field and the mass-loss rate obtained from our formalism are consistent with statistical studies as well as individual observations and give new leads to constrain the magnetic field-rotation relation. The set of scaling-laws we derived can be broadly applied to cool stars from the PMS to the end of the MS, and allow for a stellar wind modelling consistent with all the observational constraints available to date.

Though generated deep inside the convection zone, the solar magnetic field has a direct impact on the Earth space environment via the Parker spiral. It strongly modulates the solar wind in the whole heliosphere, especially its latitudinal and longitudinal speed distribution over the years. However the wind also influences the topology of the coronal magnetic field by opening the magnetic field lines in the coronal holes, which can affect the inner magnetic field of the star by altering the dynamo boundary conditions. This coupling is especially difficult to model because it covers a large variety of spatio-temporal scales. Quasi-static studies have begun to help us unveil how the dynamo-generated magnetic field shapes the wind, but the full interplay between the solar dynamo and the solar wind still eludes our understanding. We use the compressible magnetohydrodynamical (MHD) code PLUTO to compute simultaneously in 2.5D the generation and evolution of magnetic field inside the star via an alpha-Omega dynamo process and the corresponding evolution of a polytropic coronal wind over several activity cycles for a young Sun. A multi-layered boundary condition at the surface of the star connects the inner and outer stellar layers, allowing both to adapt dynamically. Our continuously coupled dynamo-wind model allows us to characterize how the solar wind conditions change as a function of the cycle phase, and also to quantify the evolution of integrated quantities such as the Alfv\'en radius. We further assess the impact of the solar wind on the dynamo itself by comparing our results with and without wind feedback.

The dissipation of the kinetic energy of wave-like tidal flows within the convective envelope of low-mass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of short-period exoplanetary systems. Although low-mass stars are magnetically active objects, the question of how the star's magnetic field impacts large-scale tidal flows and the excitation, propagation and dissipation of tidal waves still remains open. Our goal is to investigate the impact of stellar magnetism on the forcing of tidal waves, and their propagation and dissipation in the convective envelope of low-mass stars as they evolve. We have estimated the amplitude of the magnetic contribution to the forcing and dissipation of tidally induced magneto-inertial waves throughout the structural and rotational evolution of low-mass stars (from M to F-type). For this purpose, we have used detailed grids of rotating stellar models computed with the stellar evolution code STAREVOL. The amplitude of dynamo-generated magnetic fields is estimated via physical scaling laws at the base and the top of the convective envelope. We find that the large-scale magnetic field of the star has little influence on the excitation of tidal waves in the case of nearly-circular orbits and coplanar hot-Jupiter planetary systems, but that it has a major impact on the way waves are dissipated. Our results therefore indicate that a full magneto-hydrodynamical treatment of the propagation and dissipation of tidal waves is needed to properly assess the impact of star-planet tidal interactions throughout the evolutionary history of low-mass stars hosting short-period massive planets.

Current burning issues in stellar physics, for both hot and cool stars, concern their magnetism. In hot stars, stable magnetic fields of fossil origin impact their stellar structure and circumstellar environment, with a likely major role in stellar evolution. However, this role is complex and thus poorly understood as of today. It needs to be quantified with high-resolution UV spectropolarimetric measurements. In cool stars, UV spectropolarimetry would provide access to the structure and magnetic field of the very dynamic upper stellar atmosphere, providing key data for new progress to be made on the role of magnetic fields in heating the upper atmospheres, launching stellar winds, and more generally in the interaction of cool stars with their environment (circumstellar disk, planets) along their whole evolution. UV spectropolarimetry is proposed on missions of various sizes and scopes, from POLLUX on the 15-m telescope LUVOIR to the Arago M-size mission dedicated to UV spectropolarimetry.

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

Observational evidence of star-planet magnetic interactions (SPMI) in compact exo-systems have been looked for in the past decades. Indeed planets in close-in orbit can be magnetically connected to their host star, and channel Alfv\'en waves carrying large amounts of energy towards the central star. The strength and temporal modulation of SPMIs are primarily set by the magnetic topology of the host star and the orbital characteristics of the planet. As a result, SPMI signals can be modulated over the rotational period of the star, the orbital period of the planet, or a complex combination of the two. The detection of SPMI thus have to rely on multiple-epochs and multiple-wavelengths observational campaigns. We present a new method to characterize SPMIs and apply it to Kepler-78, a late G star with a super-Earth on an 8.5 hours orbit. We model the corona of Kepler-78 using the large-scale magnetic topology of the star observed with Zeeman-Doppler-Imaging. We show that the closeness of Kepler-78b allows the interaction to channel energy flux densities up to a few kW m$^{-2}$ towards the central star. We show that this flux is large enough to be detectable in classical activity tracers such as H$\alpha$. It is nonetheless too weak to explain the modulation observed by \citet{Moutou2016}. We furthermore demonstrate how to predict the temporal modulation of SPMI signals in observed systems such as Kepler-78. The methodology presented here thus paves the road towards denser, specific observational campaigns that would allow a proper identification of SPMIs in compact star-planet systems.

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

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

The energy dissipation of wave-like tidal flows in the convective envelope of low-mass stars is one of the key physical mechanisms that shape the orbital and rotational dynamics of short-period planetary systems. Tidal flows, and the excitation, propagation, and dissipation of tidally-induced inertial waves can be modified by stellar magnetic fields (e.g., Wei 2016, 2018, Lin and Ogilvie 2018). It is thus important to assess for which stars, at which location of their internal structure, and at which phase of their evolution, one needs to take into account the effects of magnetic fields on tidal waves. Using scaling laws that provide the amplitude of dynamo-generated magnetic fields along the rotational evolution of these stars (e.g., Christensen et al. 2009, Brun et al. 2015), combined with detailed grids of stellar rotation models (e.g., Amard et al. 2016), we examine the influence of a magnetic field on tidal forcing and dissipation near the tachocline of solar-like stars. We show that full consideration of magnetic fields is required to compute tidal dissipation, but not necessarily for tidal forcing.

We want to study the connections between the magnetic field generated inside the Sun and the solar wind impacting Earth, especially the influence of north-south asymmetry on the magnetic and velocity fields. We study a solar-like 11-year cycle in a quasi-static way: an asymmetric dynamo field is generated through a 2.5-dimensional (2.5-D) flux-transport model with the Babcock-Leighton mechanism, and then is used as bottom boundary condition for compressible 2.5-D simulations of the solar wind. We recover solar values for the mass loss rate, the spin-down time scale and the Alfv\'en radius, and are able to reproduce the observed delay in latitudinal variations of the wind and the general wind structure observed for the Sun. We show that the phase lag between the energy of the dipole component and the total surface magnetic energy has a strong influence on the amplitude of the variations of global quantities. We show in particular that the magnetic torque variations can be linked to topological variations during a magnetic cycle, while variations in the mass loss rate appear to be driven by variations of the magnetic energy.

The periods of magnetic activity cycles in the Sun and solar-type stars do not exhibit a simple or even single trend with respect to rotation rate or luminosity. Dynamo models can be used to interpret this diversity, and can ultimately help us understand why some solar-like stars do not exhibit a magnetic cycle, whereas some do, and for the latter what physical mechanisms set their magnetic cycle period. Three-dimensional non-linear magnetohydrodynamical simulations present the advantage of having only a small number of tunable parameters, and produce in a dynamically self-consistent manner the flows and the dynamo magnetic fields pervading stellar interiors. We conducted a series of such simulations within the EULAG-MHD framework, varying the rotation rate and luminosity of the modeled solar-like convective envelopes. We find decadal magnetic cycles when the Rossby number near the base of the convection zone is moderate (typically between 0.25 and 1). Secondary, shorter cycles located at the top of the convective envelope close to the equator are also observed in our numerical experiments, when the local Rossby number is lower than 1. The deep-seated dynamo sustained in these numerical experiments is fundamentally non-linear, in that it is the feedback of the large-scale magnetic field on the large-scale differential rotation that sets the magnetic cycle period. The cycle period is found to decrease with the Rossby number, which offers an alternative theoretical explanation to the variety of activity cycles observed in solar-like stars.

Planets in close-in orbits interact magnetically and tidally with their host stars. These interactions lead to a net torque that makes close-in planets migrate inward or outward depending on their orbital distance. We compare systematically the strength of magnetic and tidal torques for typical observed star-planet systems (T-Tauri & hot Jupiter, M dwarf & Earth-like planet, K star & hot Jupiter) based on state-of-the-art scaling-laws. We find that depending on the characteristics of the system, tidal or magnetic effects can dominate. For very close-in planets, we find that both torques can make a planet migrate on a timescale as small as 10 to 100 thousands of years. Both effects thus have to be taken into account when predicting the evolution of compact systems.

The magnetic fields of solar-type stars are observed to cycle over decadal periods -11 years in the case of the Sun. The fields originate in the turbulent convective layers of stars and have a complex dependency upon stellar rotation rate. We have performed a set of turbulent global simulations that exhibit magnetic cycles varying systematically with stellar rotation and luminosity. We find that the magnetic cycle period is inversely proportional to the Rossby number, which quantifies the influence of rotation on turbulent convection. The trend relies on a fundamentally non-linear dynamo process and is compatible with the Sun's cycle and those of other solar-type stars.

Magnetic interactions between a planet and its environment are known to lead to phenomena such as aurorae and shocks in the solar system. The large number of close-in exoplanets that were discovered triggered a renewed interest in magnetic interactions in star-planet systems. Multiple other magnetic effects were then unveiled, such as planet inflation or heating, planet migration, planetary material escape, and even modification of the host star properties. We review here the recent efforts in modelling and understanding magnetic interactions between stars and planets in the context of compact systems. We first provide simple estimates of the effects of magnetic interactions and then detail analytical and numerical models for different representative scenarii. We finally lay out a series of future developments that are needed today to better understand and constrain these fascinating interactions.

The surprising thinness of the solar tachocline is still not understood with certainty today. Among the numerous possible scenarios suggested to explain its radial confinement, one hypothesis is based on Maxwell stresses that are exerted by the cyclic dynamo magnetic field of the Sun penetrating over a skin depth below the turbulent convection zone. Our goal is to assess under which conditions (turbulence level in the tachocline, strength of the dynamo-generated field, spreading mechanism) this scenario can be realized in the solar tachocline. We develop a simplified 1D model of the upper tachocline under the influence of an oscillating magnetic field imposed from above. The turbulent transport is parametrized with enhanced turbulent diffusion (or anti-diffusion) coefficients. Two main processes that thicken the tachocline are considered; either turbulent viscous spreading or radiative spreading. An extensive parameter study is carried out to establish the physical parameter regimes under which magnetic confinement of the tachocline that is due to a surface dynamo field can be realized. We have explored a large range of magnetic field amplitudes, viscosities, ohmic diffusivities and thermal diffusivities. We find that, for large but still realistic magnetic field strengths, the differential rotation can be suppressed in the upper radiative zone (and hence the tachocline confined) if weak turbulence is present (with an enhanced ohmic diffusivity of $\eta > 10^{7-8} \, cm^2/s$), even in the presence of radiative spreading. Our results show that a dynamo magnetic field can, in the presence of weak turbulence, prevent the inward burrowing of a tachocline subject to viscous diffusion or radiative spreading.

We seek to characterize how the change of global rotation rate influences the overall dynamics and large scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the latter profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number $R_{of} = \tilde{\omega}/2\Omega_* \sim \tilde{v}/2\Omega_* R_*$, which we calibrate based on our 3-D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map ($R_{of} \gtrsim 1.5-2$) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

Young and rapidly rotating stars are known for intense, dynamo generated magnetic fields. Spectropolarimetric observations of those stars in precisely aged clusters are key input for gyrochronology and magnetochronology. We use ZDI maps of several young K-type stars of similar mass and radius but with various ages and rotational periods, to perform 3D numerical MHD simulations of their coronae and follow the evolution of their magnetic properties with age. Those simulations yield the coronal structure as well as the instant torque exerted by the magnetized, rotating wind on the star. As stars get older, we find that the angular momentum loss decreases with $\Omega^3$, which is the reason for the convergence on the Skumanich law. For the youngest stars of our sample, the angular momentum loss show signs of saturation around $8\Omega_{\odot}$, which is a common value used in spin evolution models for K-type stars. We compare these results to semi-analytical models and existing braking laws. We observe a complex wind speed distribution for the youngest stars with slow, intermediate and fast wind components, which are the result of the interaction with intense and non axisymmetric magnetic fields. Consequently, in our simulations, the stellar wind structure in the equatorial plane of young stars varies significantly from a solar configuration, delivering insight about the past of the solar system interplanetary medium.

This paper provides a brief overview of the formation of stellar fossil magnetic fields and what potential instabilities may occur given certain configurations of the magnetic field. In particular, a purely magnetic instability can occur for poloidal, toroidal, and mixed poloidal-toroidal axisymmetric magnetic field configurations as originally studied in Tayler (1973), Markey & Tayler (1973), and Tayler (1980). However, most of the magnetic field configurations observed at the surface of massive stars are non- axisymmetric. Thus, extending earlier studies of the axisymmetric Tayler instability in spherical geometry (Goossens, 1980), we introduce a formulation for the global change in the potential energy contained in a convectively-stable region given an arbitrary Lagrangian perturbation, which permits the inclusion of both axisymmetric and non-axisymmetric magnetic fields. With this tool in hand, a path is shown by which more general stability criterion can be established.

Planets in close-in orbit interact with the magnetized wind of their hosting star. This magnetic interaction was proposed to be a source for enhanced emissions in the chromosphere of the star, and to participate in setting the migration time-scale of the close-in planet. The efficiency of the magnetic interaction is know to depend on the magnetic properties of the host star, of the planet, and on the magnetic topology of the interaction. We use a global, three-dimensional numerical model of close-in star planet systems, based on the magnetohydrodynamics approximation, to compute a grid of simulations for varying properties of the orbiting planet. We propose a simple parametrization of the magnetic torque that applies to the planet, and of the energy flux generated by the interaction. The dependancy upon the planet properties and the wind properties are clearly identified in the derived scaling laws, which can be used in secular evolution codes to take into account the effect of magnetic interactions in planet migration. They can also be used to estimate a potential magnetic source of enhanced emissions in observed close-in star-planet systems, in order to constrain observationally possible exoplanetary magnetic fields.

The aim of this article is to study how the differential rotation of solar-like stars is influenced by rotation rate and mass in presence of magnetic fields generated by a convective dynamo. We use the ASH code to model the convective dynamo of solar-like stars at various rotation rates and masses, hence different effective Rossby numbers. We obtained models with either prograde (solar-like) or retrograde (anti-solar-like) differential rotation. The trends of differential rotation versus stellar rotation rate obtained for simulations including the effect of the magnetic field are weaker compared with hydro simulations ($\Delta \Omega \propto (\Omega/\Omega_{\odot})^{0.44}$ in the MHD case and $\Delta \Omega \propto (\Omega/\Omega_{\odot})^{0.89}$ in the hydro case), hence showing a better agreement with the observations. Analysis of angular momentum transport revealed that the simulations with retrograde and prograde differential rotation have opposite distribution of the viscous, turbulent Reynolds stresses and meridional circulation contributions. The thermal wind balance is achieved in the prograde cases. However, in retrograde cases Reynolds stresses are dominant for high latitudes and near the top of the convective layer. Baroclinic effects are stronger for faster rotating models.

The impressive development of global numerical simulations of turbulent stellar interiors unveiled a variety of possible differential rotation (solar or anti-solar), meridional circulation (single or multi-cellular), and dynamo states (stable large scale toroidal field or periodically reversing magnetic fields). Various numerical schemes, based on the so-called anelastic set of equations, were used to obtain these results. It appears today mandatory to assess their robustness with respect to the details of the numerics, and in particular to the treatment of turbulent sub-grid scales. We report on an ongoing comparison between two global models, the ASH and EULAG codes. In EULAG the sub-grid scales are treated implicitly by the numerical scheme, while in ASH their effect is generally modelled by using enhanced dissipation coefficients. We characterize the sub-grid scales effect in a turbulent convection simulation with EULAG. We assess their effect at each resolved scale with a detailed energy budget. We derive equivalent eddy-diffusion coefficients and use the derived diffusivities in twin ASH numerical simulations. We find a good agreement between the large-scale flows developing in the two codes in the hydrodynamic regime, which encourages further investigation in the magnetohydrodynamic regime for various dynamo solutions.

Magnetic interactions between a star and a close-in planet are postulated to be a source of enhanced emissions and to play a role in the secular evolution of the orbital system. Close-in planets generally orbit in the sub-alfv\'enic region of the stellar wind, which leads to efficient transfers of energy and angular momentum between the star and the planet. We model the magnetic interactions occurring in close-in star-planet systems with three-dimensional, global, compressible magneto-hydrodynamic numerical simulations of a planet orbiting in a self-consistent stellar wind. We focus on the cases of magnetized planets and explore three representative magnetic configurations. The Poynting flux originating from the magnetic interactions is an energy source for enhanced emissions in star-planet systems. Our results suggest a simple geometrical explanation for ubiquitous on/off enhanced emissions associated with close-in planets, and confirm that the Poynting fluxes can reach powers of the order of $10^{19}$ W. Close-in planets are also showed to migrate due to magnetic torques for sufficiently strong stellar wind magnetic fields. The topology of the interaction significantly modifies the shape of the magnetic obstacle that leads to magnetic torques. As a consequence, the torques can vary by at least an order of magnitude as the magnetic topology of the interaction varies.

We investigate the possible development of magnetohydrodynamical instabilities in the EULAG-MHD "millenium simulation" of Passos & Charbonneau (2014). This simulation sustains a large-scale magnetic cycle characterized by solar-like polarity reversals taking place on a regular multidecadal cadence, and in which zonally-oriented bands of strong magnetic field accumulate below the convective layers, in response to turbulent pumping from above in successive magnetic half-cycles. Key aspects of this simulation include low numerical dissipation and a strongly subadiabatic fluid layer underlying the convectively unstable layers corresponding to the modeled solar convection zone. These properties are conducive to the growth and development of two-dimensional instabilities otherwise suppressed by stronger dissipation. We find evidence for the action of a non-axisymmetric magnetoshear instability operating in the upper portions of the stably stratified fluid layers. We also investigate the possibility that the Tayler instability may be contributing to the destabilization of the large-scale axisymmetric magnetic component at high latitudes. On the basis of our analyses, we propose a global dynamo scenario whereby the magnetic cycle is driven primarily by turbulent dynamo action in the convecting layers, but MHD instabilities accelerate the dissipation of the magnetic field pumped down into the overshoot and stable layers, thus perhaps significantly influencing the magnetic cycle period. Support for this scenario is found in the distinct global dynamo behaviors observed in an otherwise identical EULAG-MHD simulations, using a different degree of subadiabaticity in the stable fluid layers underlying the convection zone.

Observations of surface magnetic fields are now within reach for many stellar types thanks to the development of Zeeman-Doppler Imaging. These observations are extremely useful for constraining rotational evolution models of stars, as well as for characterizing the generation of magnetic field. We recently demonstrated that the impact of coronal magnetic field topology on the rotational braking of a star can be parametrized with a scalar parameter: the open magnetic flux. However, without running costly numerical simulations of the stellar wind, reconstructing the coronal structure of the large scale magnetic field is not trivial. An alternative -broadly used in solar physics- is to extrapolate the surface magnetic field assuming a potential field in the corona, to describe the opening of the field lines by the magnetized wind. This technique relies on the definition of a so-called source surface radius, which is often fixed to the canonical value of 2.5Rsun. However this value likely varies from star to star. To resolve this issue, we use our extended set of 2.5D wind simulations published in 2015, to provide a criteria for the opening of field lines as well as a simple tool to assess the source surface radius and the open magnetic flux. This allows us to derive the magnetic torque applied to the star by the wind from any spectropolarimetric observation. We conclude by discussing some estimations of spin-down time scales made using our technique, and compare them to observational requirements.

The revolution of helio- and asteroseismology provides access to the detailed properties of stellar interiors by studying the star's oscillation modes. Among them, gravity (g) modes are formed by constructive interferences between progressive internal gravity waves (IGWs), propagating in stellar radiative zones. Our new 3D nonlinear simulations of the interior of a solar-like star allows us to study the excitation, propagation, and dissipation of these waves. The aim of this article is to clarify our understanding of the behavior of IGWs in a 3D radiative zone and to provide a clear overview of their properties. We use a method of frequency filtering that reveals the path of {individual} gravity waves of different frequencies in the radiative zone. We are able to identify the region of propagation of different waves in 2D and 3D, to compare them to the linear raytracing theory and to distinguish between propagative and standing waves (g modes). We also show that the energy carried by waves is distributed in different planes in the sphere, depending on their azimuthal wave number. We are able to isolate individual IGWs from a complex spectrum and to study their propagation in space and time. In particular, we highlight in this paper the necessity of studying the propagation of waves in 3D spherical geometry, since the distribution of their energy is not equipartitioned in the sphere.

Stellar winds are thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. 60 simulations made with a 2.5D, cylindrical and axisymmetric set-up and computed with the PLUTO code were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phase as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486- 4943-1) whose topology has been obtained by Zeeman-Doppler Imaging (ZDI).

The braking of main sequence stars originates mainly from their stellar wind. The efficiency of this angular momentum extraction depends on the rotation rate of the star, the acceleration profile of the wind and the coronal magnetic field. The derivation of scaling laws parametrizing the stellar wind torque is important for our understanding of gyro-chronology and the evolution of the rotation rates of stars. In order to understand the impact of complex magnetic topologies on the stellar wind torque, we present three- dimensional, dynamical simulations of the corona of HD 189733. Using the observed complex topology of the magnetic field, we estimate how the torque associated with the wind scales with model parameters and compare those trends to previously published scaling laws.

Stellar winds are thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. 60 simulations made with a 2.5D, cylindrical and axisymmetric set-up and computed with the PLUTO code were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phase as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486-4943-1) whose topology has been obtained by Zeeman-Doppler Imaging (ZDI). Ideas about how to compute the open flux from ZDI Maps are discussed.

This paper explores and compares the pitfalls of modelling the three-dimensional wind of a spherical star with a cartesian grid. Several numerical methods are compared, using either uniform and stretched grid or adaptative mesh refinement (AMR). An additional numerical complication is added, when an orbiting planet is considered. In this case a rotating frame is added to the model such that the orbiting planet is at rest in the frame of work. The three-dimensional simulations are systematically compared to an equivalent two-dimensional, axisymmetric simulation. The comparative study presented here suggests to limit the rotation rate of the rotating frame below the rotating frame of the star and provides guidelines for further three-dimensional modelling of stellar winds in the context of close-in star-planet interactions.

Magnetic interactions between close-in planets and their host star can play an important role in the secular orbital evolution of the planets, as well as the rotational evolution of their host. As long as the planet orbits inside the Alfv\'en surface of the stellar wind, the magnetic interaction between the star and the planet can modify the wind properties and also lead to direct angular momentum transfers between the two. We model these star-planet interactions using compressible magneto-hydrodynamic (MHD) simulations, and quantify the angular momentum transfers between the star, the planet, and the stellar wind. We study the cases of magnetized and non-magnetized planets and vary the orbital radius inside the Alfv\'en surface of the stellar wind. Based on a grid of numerical simulations, we propose general scaling laws for the modification of the stellar wind torque, for the torque between the star and the planet, and for the planet migration associated with the star-planet magnetic interactions. We show that when the coronal magnetic field is large enough and the star is rotating sufficiently slowly, the effect of the magnetic star-planet interaction is comparable to tidal effects and can lead to a rapid orbital decay.

Coronal Mass Ejections (CMEs) and Corotating Interaction Regions (CIRs) are major sources of magnetic storms on Earth and are therefore considered to be the most dangerous space weather events. The Observatories of Solar Corona and Active Regions (OSCAR) mission is designed to identify the 3D structure of coronal loops and to study the trigger mechanisms of CMEs in solar Active Regions (ARs) as well as their evolution and propagation processes in the inner heliosphere. It also aims to provide monitoring and forecasting of geo-effective CMEs and CIRs. OSCAR would contribute to significant advancements in the field of solar physics, improvements of the current CME prediction models, and provide data for reliable space weather forecasting. These objectives are achieved by utilising two spacecraft with identical instrumentation, located at a heliocentric orbital distance of 1~AU from the Sun. The spacecraft will be separated by an angle of 68$^{\circ}$ to provide optimum stereoscopic view of the solar corona. We study the feasibility of such a mission and propose a preliminary design for OSCAR.

We assess the predictive capabilities of various classes of avalanche models for solar flares. We demonstrate that avalanche models cannot generally be used to predict specific events due to their high sensitivity to their embedded stochastic process. We show that deterministically driven models can nevertheless alleviate this caveat and be efficiently used for large events predictions. Our results promote a new approach for large (typically X-class) solar flares predictions based on simple and computationally inexpensive avalanche models.

Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.

We develop and discuss the properties of a new class of lattice-based avalanche models of solar flares. These models are readily amenable to a relatively unambiguous physical interpretation in terms of slow twisting of a coronal loop. They share similarities with other avalanche models, such as the classical stick--slip self-organized critical model of earthquakes, in that they are driven globally by a fully deterministic energy loading process. The model design leads to a systematic deficit of small scale avalanches. In some portions of model space, mid-size and large avalanching behavior is scale-free, being characterized by event size distributions that have the form of power-laws with index values, which, in some parameter regimes, compare favorably to those inferred from solar EUV and X-ray flare data. For models using conservative or near-conservative redistribution rules, a population of large, quasiperiodic avalanches can also appear. Although without direct counterparts in the observational global statistics of flare energy release, this latter behavior may be relevant to recurrent flaring in individual coronal loops. This class of models could provide a basis for the prediction of large solar flares.

We model the magnetized interaction between a star and a close-in planet (SPMIs), using global, magnetohydrodynamic numerical simulations. In this proceedings, we study the effects of the numerical boundary conditions at the stellar surface, where the stellar wind is driven, and in the planetary interior. We show that is it possible to design boundary conditions that are adequate to obtain physically realistic, steady-state solutions for cases with both magnetized and unmagnetized planets. This encourages further development of numerical studies, in order to better constrain and understand SPMIs, as well as their effects on the star-planet rotational evolution.

We present 2D magnetohydrodynamic simulations performed with the PLUTO code to model magnetized star-planet interactions. We study two simple scenarios of magnetized star-planet interactions: the unipolar and dipolar interactions. Despite the simplified hypotheses we consider in the model, the qualitative behavior of the interactions is very well recovered. These encouraging results promote further developments of the model to obtain predictions on the effect and the physical manifestation of magnetized star--close-in planet interactions.

We present a method to characterize the spectral transfers of magnetic energy between scales in simulations of stellar convective dynamos. The full triadic transfer functions are computed thanks to analytical coupling relations of spherical harmonics based on the Clebsch-Gordan coefficients. The method is applied to mean field $\alpha\Omega$ dynamo models as benchmark tests. From the physical standpoint, the decomposition of the dynamo field into primary and secondary dynamo families proves very instructive in the $\alpha\Omega$ case. The same method is then applied to a fully turbulent dynamo in a solar convection zone, modeled with the 3D MHD ASH code. The initial growth of the magnetic energy spectrum is shown to be non-local. It mainly reproduces the kinetic energy spectrum of convection at intermediate scales. During the saturation phase, two kinds of direct magnetic energy cascades are observed in regions encompassing the smallest scales involved in the simulation. The first cascade is obtained through the shearing of magnetic field by the large scale differential rotation that effectively cascades magnetic energy. The second is a generalized cascade that involves a range of local magnetic and velocity scales. Non-local transfers appear to be significant, such that the net transfers cannot be reduced to the dynamics of a small set of modes. The saturation of the large scale axisymmetric dipole and quadrupole are detailed. In particular, the dipole is saturated by a non-local interaction involving the most energetic scale of the magnetic energy spectrum, which points out the importance of the magnetic Prandtl number for large-scale dynamos.

The reason for the observed thinness of the solar tachocline is still not well understood. One of the explanations that have been proposed is that a primordial magnetic field renders the rotation uniform in the radiation zone. We test here the validity of this magnetic scenario through 3D numerical MHD simulations that encompass both the radiation zone and the convection zone. The numerical simulations are performed with the anelastic spherical harmonics (ASH) code. The computational domain extends from $0.07\;R_\odot$ to $0.97\;R_\odot$. In the parameter regime we explored, a dipolar fossil field aligned with the rotation axis can not remain confined in the radiation zone. When the field lines are allowed to interact with turbulent unstationary convective motions at the base of the convection zone, 3D effects prevent the field confinement. In agreement with previous work, we find that a dipolar fossil field, even when it is initially buried deep inside the radiation zone, will spread into the convective zone. According to Ferraro's law of iso-rotation, it then imprints on the radiation zone the latitudinal differential rotation of the convection zone, which is not observed.