Search

M-type stars are the most common stars in the universe. They are ideal hosts for the search of exoplanets in the habitable zone (HZ), as their small size and low temperature make the HZ much closer in than their solar twins. Harboring very deep convective layers, they also usually exhibit very intense magnetic fields. Understanding their environment, in particular their coronal and wind properties, is thus very important, as they might be very different from what is observed in the solar system. The mass loss rate of M-type stars is poorly known observationally, and recent attempts to estimate it for some of them (TRAPPIST-1, Proxima Cen) can vary by an order of magnitude. In this work, we revisit the stellar wind properties of M-dwarfs in the light of the latest estimates of $\dot{M}$ through Lyman-$\alpha$ absorption at the astropause and slingshot prominences. We outline a modeling strategy to estimate the mass loss rate, radiative loss and wind speed, with uncertainties, based on an Alfv\'en wave driven stellar wind model. We find that it is very likely that several TRAPPIST-1 planets lie within the Alfv\'en surface, which imply that these planets experience star-planet magnetic interactions (SPMI). We also find that SPMI between Proxima Cen b and its host star could be the reason of recently observed radio emissions., Comment: 15 pages, 9 figures, accepted for publication in ApJ

We present a new multi-dimensional, robust, and cell-centered finite-volume scheme for the ideal MHD equations. This scheme relies on relaxation and splitting techniques and can be easily used at high order. A fully conservative version is not entropy satisfying but is observed experimentally to be more robust than standard constrained transport schemes at low plasma beta. At very low plasma beta and high Alfv\'en number, we have designed an entropy-satisfying version that is not conservative for the magnetic field but preserves admissible states and we switch locally a-priori between the two versions depending on the regime of plasma beta and Alfv\'en number. This strategy is robust in a wide range of standard MHD test cases, all performed at second order with a classic MUSCL-Hancock scheme., Comment: accepted in JCP

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

Solar Orbiter provides unique capabilities to understand the heliosphere. In particular, it has made observations of the far-side of the Sun and provides unique information to improve space weather monitoring. We aim to quantify how far-side data will affect simulations of the corona and the interplanetary medium, especially in the context of space weather forecasting. We focused on a time period with a single sunspot emerging on the far-side in February 2021. We used two different input magnetic maps: one with the far-side active region and one without. We used three different coronal models: a semi-empirical model (potential field source surface or PFSS) and two different magnetohydrodynamic models (Wind Predict and Wind Predict-AW). We compared all the models with both remote sensing and in situ observations. We find that the inclusion of the far-side active region in the various models has a small local impact due to the limited amount of flux of the sunspot (at most 8% of the total map flux), which leads to coronal hole changes of around 7% for all models. Interestingly, there is a more global impact on the magnetic structure seen in the current sheet, with clear changes in the coronal hole boundaries visible in extreme ultra-violet (EUV) on the western limb. For the Wind Predict-AW model, we demonstrate that the inclusion of the far-side data improves both the structure of the streamers and the connectivity to the spacecraft. In conclusion, the inclusion of a single far-side active region may have a small local effect with respect to the total magnetic flux, but it has global effects on the magnetic structure, and thus it must be taken into account to accurately describe the Sun-Earth connection. The flattening of the heliospheric current sheet for all models reveals an increase of the source surface height, which affects the open and closed magnetic field line distributions., Comment: 22 pages, 12 figures, accepted in A&A

Aims. In this study, we analyse the magnetic field properties of a set of 15 global magnetohydrodynamics (MHD) simulations of solar-type star dynamos conducted using the ASH code. Our objective is to enhance our understanding of these properties by comparing theoretical results to current observations, and to finally provide fresh insights into the field. Methods. We analysed the rotational and magnetic properties as a function of various stellar parameters (mass, age and rotation rate) in a 'Sun in time' approach in our extended set of 3D MHD simulations. To facilitate direct comparisons with stellar magnetism observations using various Zeeman-effect techniques, we decomposed the numerical data into vectorial spherical harmonics. Results. A comparison of the trends we find in our simulations set reveals a promising overall agreement with the observational context of stellar magnetism, enabling us to suggest a plausible scenario for the magneto-rotational evolution of solar-type stars. In particular, we find that the magnetic field may reach a minimum amplitude at a transition value of the Rossby number near unity. This may have important consequences on the long-term evolution of solar-type stars, by impacting the relation between stellar age, rotation and magnetism. This supports the need for future observational campaigns, especially for stars in the high Rossby number regime., Comment: 10 pages + Appendix, 10 Figures, Accepted for publication in A&A

The solar tachocline is an internal region of the Sun possessing strong radial and latitudinal shears straddling the base of the convective envelope. Based on helioseismic inversions, the tachocline is known to be thin (less than 5\% of the solar radius). Since the first theory of the solar tachocline in 1992, this thinness has not ceased to puzzle solar physicists. In this review, we lay out the grounds of our understanding of this fascinating region of the solar interior. We detail the various physical mechanisms at stake in the solar tachocline, and put a particular focus on the mechanisms that have been proposed to explain its thinness. We also examine the full range of MHD processes including waves and instabilies that are likely to occur in the tachocline, as well as their possible connection with active region patterns observed at the surface. We reflect on the most recent findings for each of them, and highlight the physical understanding that is still missing and that would allow the research community to understand, in a generic sense, how the solar tachocline and stellar tachocline are formed, are sustained, and evolve on secular timescales., Comment: 46 pages, 10 figures, accepted for publication in Space Science Reviews

Gravity modes (g modes), mixed gravito-acoustic modes (mixed modes), and gravito-inertial modes (gi modes) possess unmatched properties as probes for stars with radiative interiors. The structural and dynamical constraints that they are able to provide cannot be accessed by other means. While they provide precious insights into the internal dynamics of evolved stars as well as massive and intermediate-mass stars, their non-detection in main sequence (MS) solar-type stars make them a crucial missing piece in our understanding of angular momentum transport in radiative zones and stellar rotational evolution. In this work, we aim to apply certain analysis tools originally developed for helioseismology in order to look for g-mode signatures in MS solar-type stars. We select a sample of the 34 most promising MS solar-type stars with Kepler four-year long photometric time series. All these stars are well-characterised late F-type stars with thin convective envelopes, fast convective flows, and stochastically excited acoustic modes (p modes). For each star, we compute the background noise level of the Fourier power spectrum to identify significant peaks at low frequency. After successfully detecting individual peaks in 12 targets, we further analyse four of them and observe distinct patterns of surrounding peaks with a low probability of being noise artifacts. Comparisons with the predictions from reference models suggest that these patterns are compatible with the presence of non-asymptotic low-order pure g modes, pure p modes, and mixed modes. Given their sensitivity to both the convective core interface stratification and the coupling between p- and g-mode resonant cavities, such modes are able to provide strong constraints on the structure and evolutionary states of the related targets. [abridged], Comment: 19 pages, 19 figures, accepted for publication in A&A

The strength and morphology of the Sun's magnetic field evolves significantly during the solar cycle, with the overall polarity of the Sun's magnetic field reversing during the maximum of solar activity. Long-term changes are also observed in sunspot and geomagnetic records, however systematic magnetic field observations are limited to the last four cycles. We investigate the long-term evolution of the Sun's magnetic field, and the influence this has on the topology and rotation of the solar corona. The Sun's photospheric magnetic field was decomposed into spherical harmonics using synoptic Carrington magnetograms from 1) WSO, 2) MDI onboard the SOHO, and 3) HMI onboard the SDO. The time-evolution of the spherical harmonic coefficients was used to explore the variation of the Sun's magnetic field, focusing on the large-scale modes. PFSS extrapolations of the photospheric field were computed to follow topological changes in the corona. The footpoints of the Sun's open magnetic field vary between the polar coronal holes and activity driven features such as active regions, and equatorial coronal holes. Consequently, the mean rotation rate of the solar wind is modulated during each cycle by the latitudinal variation of open field footpoints, with slower rotation during minima and faster (Carrington-like) rotation during maxima. Thisc variation is sensitive to cycle to cycle differences in the polar field strengths and hemispherical flux emergence rates, with the ratio of quadrupole to dipole energy following a similar variation. Cycle 23 maintained a larger fraction of quadrupolar energy in the declining phase, which kept the sources of open magnetic flux closer to the equator, extending the period of faster equator-ward connectivity. The ratio of quadrupole to dipole energy could be a useful proxy when examining the impact of differential rotation on the coronae of other Sun-like stars., Comment: Accepted to A&A. 13 Pages + Appendix. 11 Figures + 7 Appendix Figures

We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scale dynamo. We continue to studies that model the Sun at different ages and to studies of stars of different masses and evolutionary stages. Both simulations and observations indicate that rotation, measured by the Rossby number which is the ratio of rotation period and convective turnover time, is a key ingredient in setting the overall level and characteristics of magnetic activity. Finally, efforts to understand global 3D simulations in terms of mean-field dynamo theory are discussed., Comment: 70 pages, 9 figures, submitted to Space Science Reviews

The arrival of powerful instruments will provide valuable data for the characterization of rocky exoplanets. It is then crucial to accurately model the dynamical state of exoplanets. Rocky planets with sufficiently large orbits should have non-zero eccentricities and/or obliquities. Realistic models of tides for rocky planets can allow for higher spin states than the synchronization state in the presence of eccentricities or obliquities. This work explores the secular evolution of a star-planet system under tidal interactions, both gravitational and thermal, induced respectively by the quadrupolar component of the gravitational potential and the irradiation of the planet's surface. We use the formalism of Kaula associated with an Andrade rheology to model a relevant response of a rocky planet to gravitational tides and a prescription of thermal tides fitted for Venus to model the response of the atmosphere to the thermal tides. We implemented the general secular evolution equations of tidal interactions in the secular code ESPEM (French acronym for Evolution of Planetary System and Magnetism). We show the possible spin-orbit evolution and resonances for eccentric orbits and explore the possible spin orbit resonances raised by the obliquity of the planet. Our simulations have shown that the secular evolution of the spin and obliquity can lead to the retrograde spin of the Venus-like planet if the system starts from a high spin obliquity, in agreement to previous studies. Taking into account the luminosity evolution of the Sun changes the picture. We find that the planet never reaches the equilibrium: the timescale of rotation evolution is longer than the luminosity variation timescale, which suggests that Venus may never reach a spin equilibrium state but may still evolve., Comment: 19 pages, 16 figues

Sun-like stars shed angular momentum due to the presence of magnetised stellar winds. Magnetohydrodynamic models have been successful in exploring the dependence of this "wind-braking torque" on various stellar properties, however the influence of surface differential rotation is largely unexplored. As the wind-braking torque depends on the rotation rate of the escaping wind, the inclusion of differential rotation should effectively modulate the angular momentum-loss rate based on the latitudinal variation of wind source regions. In order to quantify the influence of surface differential rotation on the angular momentum-loss rate of the Sun, we exploit the dependence of the wind-braking torque on the effective rotation rate of the coronal magnetic field. This quantity is evaluated by tracing field lines through a Potential Field Source Surface (PFSS) model, driven by ADAPT-GONG magnetograms. The surface rotation rates of the open magnetic field lines are then used to construct an open-flux weighted rotation rate, from which the influence on the wind-braking torque can be estimated. During solar minima, the rotation rate of the corona decreases with respect to the typical solid-body rate (the Carrington rotation period is 25.4 days), as the sources of the solar wind shift towards the slowly-rotating poles. With increasing activity, more solar wind emerges from the Sun's active latitudes which enforces a Carrington-like rotation. The effect of differential rotation on the Sun's current wind-braking torque is found to be small. The wind-braking torque is ~10-15% lower during solar minimum, than assuming solid body rotation, and a few percent larger during solar maximum. For more rapidly-rotating Sun-like stars, differential rotation may play a more significant role, depending on the configuration of the large-scale magnetic field., Comment: Accepted to A&A. 12 Pages + Appendix. 9 Figures + 4 Appendix Figures

Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere., Comment: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

There are no strong constraints placed thus far on the amplitude of internal gravity waves (IGWs) that are stochastically excited in the radiative interiors of solar-type stars. Late F-type stars have relatively thin convective envelopes with fast convective flows and tend to be fast rotators compared to solar-type stars of later spectral types. These two elements are expected to directly impact the IGW excitation rates and properties. We want to estimate the amplitude of stochastically excited gravity modes (g-modes) in F-type stars for different rotational regimes. We used the ASH code to perform 3D simulations of deep-shell models of 1.3 $M_\odot$ F-type solar-type stars, including the radiative interior and the shallow convective envelope. The IGWs are excited by interface interactions between convective plumes and the top of the radiative interior. We were able to characterise the IGWs and g-mode properties in the radiative interior, and we compared these properties using the computation from the 1D oscillation code GYRE. The amplitude of low-frequency modes is significantly higher in fast-rotating models and the evolution of the period spacing of consecutive modes exhibits evidence of a behaviour that is modified by the influence of the Coriolis force. For our fastest rotating model, we were able to detect the intermediate degree g-mode signature near the top of the simulation domain. Nevertheless, the predicted luminosity perturbations from individual modes still remain at small amplitudes. We obtained mode amplitudes that are several orders of magnitude higher than those of prior 3D simulations of solar models. Our simulations suggest that g-mode signatures could be detectable in late F-type stars. [abridged], Comment: 22 pages, 28 figures. Accepted in A&A

Anti-solar differential rotation profiles have been found for decades in numerical simulations of convective envelopes of solar-type stars. These profiles are characterized by a slow equator and fast poles (i.e., reversed with respect to the Sun) and have been found in simulations for high Rossby numbers (slow rotators). Rotation profiles like this have been reported observationally in evolved stars, but have never been unambiguously observed for cool solar-type stars on the main sequence. In this context, detecting this regime in main-sequence solar-type stars would improve our understanding of their magnetorotational evolution. The goal of this study is to identify the most promising cool main-sequence stellar candidates for anti-solar differential rotation in the \textit{Kepler} sample. First, we introduce a new theoretical formula to estimate fluid Rossby numbers, $Ro_{\rm f}$, of main-sequence solar-type stars, from observational quantities, and taking the influences of the internal structure and metallicity into account. We obtain a list of the most promising stars that are likely to show anti-solar differential rotation. We identify two samples: one at solar metallicity, including 14 targets, and another for other metallicities, including 8 targets. We find that the targets with the highest $Ro_{\rm f}$ are likely to be early-G or late-F stars at about log$_{10}g=4.37$~dex. We conclude that cool main-sequence stellar candidates for anti-solar differential rotation exist in the \textit{Kepler} sample. The most promising candidate is KIC~10907436, and two other particularly interesting candidates are the solar analog KIC~7189915 and the seismic target KIC~12117868. Future characterization of these 22 stars is expected to help us understand how dynamics can impact magnetic and rotational evolution of old solar-type stars at high Rossby number., Comment: 12 pages + Appendix ; 9 Figures ; 3 Tables ; Accepted for publication in A&A

The physical properties responsible for the formation and evolution of the corona and heliosphere are still not completely understood. 3D MHD global modeling is a powerful tool to investigate all the possible candidate processes. To fully understand the role of each of them, we need a validation process where the output from the simulations is quantitatively compared to the observational data. In this work, we present the results from our validation process applied to the wave turbulence driven 3D MHD corona-wind model WindPredict-AW. At this stage of the model development, we focus the work to the coronal regime in quiescent condition. We analyze three simulations results, which differ by the boundary values. We use the 3D distributions of density and temperature, output from the simulations at the time of around the First Parker Solar Probe perihelion (during minimum of the solar activity), to synthesize both extreme ultraviolet (EUV) and white light polarized (WL pB) images to reproduce the observed solar corona. For these tests, we selected AIA 193 A, 211 A and 171 A EUV emissions, MLSO K-Cor and LASCO C2 pB images obtained the 6 and 7 November 2018. We then make quantitative comparisons of the disk and off limb corona. We show that our model is able to produce synthetic images comparable to those of the observed corona., Comment: in press

We use the ASH code to model the convective dynamo of solar-type stars. Based on a series of 15 3-D MHD simulations spanning 4 bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simulations. We find a weaker trend between differential rotation and stellar rotation rate, ($\Delta \Omega \sim (|\Omega|/\Omega_{\odot})^{0.46}$) in the MHD solutions than in their HD counterpart $(|\Omega|/\Omega_{\odot})^{0.66})$, yielding a better agreement with the observational trends based on power laws. We find that for a fluid Rossby number between $0.15 \lesssim Ro_f \lesssim 0.65$ the solutions possess long magnetic cycle, if $Ro_f \lesssim 0.42$ a short cycle and if $Ro_f \gtrsim 1$ (anti-solar-like differential rotation) a statistically steady state. We show that short-cycle dynamos follow the classical Parker-Yoshimura rule whereas the long-cycle period ones do not. We further demonstrate that the Rossby number dependency of the large-scale surface magnetic field in the simulation ($B_{L,surf} \sim Ro_{f}^{-1.26}$) agrees better with observations ($B_{V} \sim Ro_{s}^{-1.4 \pm 0.1}$) and differs from dynamo scaling based on the global magnetic energy ($B_{bulk} \sim Ro_{f}^{-0.5}$). We also show that up to few percents of the stellar luminosity can be channelled into the star's magnetism, hence providing a large energy reservoir for possible surface eruptive events., Comment: 45 pages, 24 figures, 7 tables, 3 appendices, accepted by ApJ on 12/23/21, published on 02/03/22. (1) DAp/AIM, CEA/IRFU, CNRS/INSU, Univ. Paris-Saclay \& Univ. de Paris, France, (2) CMPA, KU Leuven, Belgium, (3) Universidad Carlos III de Madrid, Spain, (4) Dept. de physique, Univ. de Montr\'eal, Canada, (5) JILA and Dept. of Astrophysical and Planetary Sciences, University of Colorado, USA

Over the course of their lifetimes, the rotation of solar-type stars goes through different phases. Once they reach the zero-age main sequence, their global rotation rate decreases during the main sequence until at least the solar age, approximately following the empirical Skumanich's law and enabling gyrochronology. Older solar-type stars might then reach a point of transition when they stop braking, according to recent results of asteroseismology. Additionally, recent 3D numerical simulations of solar-type stars show that different regimes of differential rotation can be characterized with the Rossby number. In particular, anti-solar differential rotation (fast poles, slow equator) may exist for high Rossby number (slow rotators). If this regime occurs during the main sequence and, in general, for slow rotators, we may consider how magnetic generation through the dynamo process might be impacted. In particular, we consider whether slowly rotating stars are indeed subject to magnetic cycles. We find that kinematic $\alpha$ $\Omega$ dynamos allow for the presence of magnetic cycles and global polarity reversals for both rotation regimes, but only if the $\alpha$-effect is saddled on the tachocline. If it is distributed in the convection zone, solar-type cases still possess a cycle and anti-solar cases do not. Conversely, we have not found any possibility for sustaining a magnetic cycle with the traditional Babcock-Leighton flux-transport dynamos in the anti-solar differential rotation regime due to flux addition. Graphic interpretations are proposed in order to illustrate these cases. However, we find that hybrid models containing both prescriptions can still sustain local polarity reversals at some latitudes. We conclude that stars in the anti-solar differential rotation regime can sustain magnetic cycles only for very specific dynamo processes., Comment: 18 pages and 16 Figures, accepted in A&A

The oscillations of a slowly rotating star have long been classified into spheroidal and toroidal modes. The spheroidal modes include the well-known 5-min acoustic modes used in helioseismology. Here we report observations of the Sun's toroidal modes, for which the restoring force is the Coriolis force and whose periods are on the order of the solar rotation period. By comparing the observations with the normal modes of a differentially rotating spherical shell, we are able to identify many of the observed modes. These are the high-latitude inertial modes, the critical-latitude inertial modes, and the equatorial Rossby modes. In the model, the high-latitude and critical-latitude modes have maximum kinetic energy density at the base of the convection zone, and the high-latitude modes are baroclinically unstable due to the latitudinal entropy gradient. As a first application of inertial-mode helioseismology, we constrain the superadiabaticity and the turbulent viscosity in the deep convection zone., Comment: Forthcoming article, Astronomy and Astrophysics

Over the last two decades, a large population of close-in planets has been detected around a wide variety of host stars. Such exoplanets are likely to undergo planetary migration through magnetic and tidal interactions. We aim to follow the orbital evolution of a planet along the structural and rotational evolution of its host star by taking into account simultaneously tidal and magnetic torques, in order to explain some properties of the distribution of observed close-in planets. We rely on a numerical model of a coplanar circular star-planet system taking into account stellar structural changes, wind braking and star-planet interactions, called ESPEM. We browse the parameter space of star-planet systems' configurations and assess the relative influence of magnetic and tidal torques on secular evolution. We then synthesize star-planet populations and confront their distribution in orbital and stellar rotation periods to Kepler satellite data. First, we find that after the dissipation of the protoplanetary disk, both types of interactions can dominate secular evolution depending on the initial configuration of the system and the evolutionary phase considered. Moreover, different populations of star-planet systems emerge from the combined action of both kinds of interactions, according to the evolutionary phase during which the planet migrates significantly. This may affect significantly the detectability of star-planet systems as well as the validity of gyrochonology. All in all, star-planet interactions significantly impact the global distribution in orbital periods during the main sequence, while the global distribution in stellar rotation periods is marginally affected. More precisely, star-planet magnetic interactions significantly affect the distribution of super-Earths around slowly rotating stars, while tidal effects are found to shape the distribution of giant planets., Comment: 28 pages, 24 figures, accepted for publication in A&A

We focus on the connection between the internal dynamo magnetic field and the stellar wind. If the star has a cyclic dynamo, the modulations of the magnetic field can affect the wind, which in turn can back-react on the boundary conditions of the star, creating a feedback loop. We have developed a 2.5-dimensional numerical set-up to model this essential coupling. We have implemented an alpha-Omega mean-field dynamo in the PLUTO code, and then coupled it to a spherical polytropic wind model via an interface composed of four grid layers with dedicated boundary conditions. We present here a dynamo model close to a young Sun with cyclic magnetic activity. First we show how this model allows to track the influence of the dynamo activity on the corona by displaying the correlation between the activity cycle, the coronal structure and the time evolution of integrated quantities. Then we add the feedback of the wind on the dynamo and discuss the changes observed in the dynamo symmetry and the wind variations. We explain these changes in terms of dynamo modes : in this parameter regime, the feedback loop leads to a coupling between the dynamo families via a preferred growth of the quadrupolar mode. We also study our interface in terms of magnetic helicity, and show that it leads to a small injection in the dynamo. This model confirms the importance of coupling physically internal and external stellar layers, as it has a direct impact on both the dynamo and the wind., Comment: 30 pages, 19 figures, accepted in the Astrophysical Journal

We study the behaviour and properties of the solar wind using a 2.5D Alfv\'en wave driven wind model. We first systematically compare the results of an Alfv\'en wave (AW) driven wind model with a polytropic approach. Polytropic magnetohydrodynamic wind models are thermally driven, while Alfv\'en waves act as additional acceleration and heating mechanisms in the Alfv\'en wave driven model. We confirm that an AW-driven model is required to reproduce the observed bimodality of slow and fast solar winds. We are also able to reproduce the observed anti-correlation between the terminal wind velocity and the coronal source temperature with the AW-driven wind model. We also show that the wind properties along an eleven-year cycle differ significantly from one model to the other. The AW-driven model again shows the best agreement with observational data. Indeed, solar surface magnetic field topology plays an important role in the Alfv\'en wave driven wind model, as it enters directly into the input energy sources via the Poynting flux. On the other hand, the polytropic wind model is driven by an assumed pressure gradient; thus it is relatively less sensitive to the surface magnetic field topology. Finally, we note that the net torque spinning down the Sun exhibits the same trends in the two models, showing that the polytropic approach still captures correctly the essence of stellar winds., Comment: 33 Pages, 16 Figures, Accepted in Astrophysical Journal

Cosmic rays (CRs) are tracers of solar events when they are associated with solar flares, but also galactic events when they come from outside our solar system. SEPs are correlated with the 11-year solar cycle while GCRs are anti-correlated due to their interaction with the heliospheric magnetic field and the solar wind. Our aim is to quantify separately the impact of the amplitude and the geometry of the magnetic field on the propagation of CRs of various energies in the inner heliosphere. We focus especially on the diffusion caused by the magnetic field along and across the field lines. To do so, we use the results of 3D MHD wind simulations running from the lower corona up to 1 AU. The wind is modeled using a polytropic approximation, and fits and power laws are used to account for the turbulence. Using these results, we compute the parallel and perpendicular diffusion coefficients of the Parker CR transport equation, yielding 3D maps of the diffusion of CRs in the inner heliosphere. By varying the amplitude of the magnetic field, we change the amplitude of the diffusion by the same factor, and the radial gradients by changing the spread of the current sheet. By varying the geometry of the magnetic field, we change the latitudinal gradients of diffusion by changing the position of the current sheets. By varying the energy, we show that the distribution of values for SEPs is more peaked than GCRs. For realistic solar configurations, we show that diffusion is highly non-axisymmetric due to the configuration of the current sheets, and that the distribution varies a lot with the distance to the Sun with a drift of the peak value. This study shows that numerical simulations and theory can help quantify better the influence of the various magnetic field parameters on the propagation of CRs. This study is a first step towards generating synthetic CR rates from numerical simulations., Comment: 29 pages, 12 figures, accepted for topical issue "Space climate: The past and future of solar activity" of Journal of Space Weather and Space Climate

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., Comment: accepted for publication in The Astrophysical Journal

Predictions of solar cycle 24 obtained from advection-dominated and diffusion-dominated kinematic dynamo models are different if the Babcock-Leighton mechanism is the only source of the poloidal field. Yeates et al. (2008) argue that the discrepancy arises due to different memories of the solar dynamo for advection- and diffusion-dominated solar convection zones. We aim to investigate the differences in solar cycle memory obtained from advection-dominated and diffusion-dominated kinematic solar dynamo models. Specifically, we explore whether the inclusion of Parker's mean-field $\alpha$ effect, in addition to the Babcock-Leighton mechanism, has any impact on the memory of the solar cycle. We used a kinematic flux transport solar dynamo model where poloidal field generation takes place due to both the Babcock-Leighton mechanism and the mean-field $\alpha$ effect. We additionally considered stochastic fluctuations in this model and explored cycle-to-cycle correlations between the polar field at the minima and toroidal field at cycle maxima. Solar dynamo memory is always limited to only one cycle in diffusion-dominated dynamo regimes while in advection-dominated regimes the memory is distributed over a few solar cycles. However, the addition of a mean-field alpha effect reduces the memory of the solar dynamo to within one cycle in the advection-dominated dynamo regime when there are no fluctuations in the mean-field $\alpha$ effect. When fluctuations are introduced in the mean-field poloidal source a more complex scenario is evident, with very weak but significant correlations emerging across a few cycles. Our results imply that inclusion of a mean-field alpha effect in the framework of a flux transport Babcock-Leighton dynamo model leads to additional complexities that may impact memory and predictability of predictive dynamo models of the solar cycle., Comment: Accepted in Astronomy and Astrophysics, 12 pages, 5 figures

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

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., Comment: 5 pages, 2 figures, to appear in the Proceedings of IAU Symposium No. 354 - Solar and Stellar Magnetic Fields: Origins and Manifestations

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., Comment: 9 pages, 5 figures, IAUS 354

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., Comment: White paper submitted on the use of HST for TESS follow-up observations, 5 pages, 2 figures

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., Comment: To appear in the Proceedings of IAU Symposium No. 345 - Origins: from the Protosun to the First Steps of Life

Differential rotation plays a key role in stellar evolution by triggering hydrodynamical instabilities and large-scale motions that induce transport of chemicals and angular momentum and by modifying the propagation and the frequency spectrum of gravito-inertial waves. It is thus crucial to investigate its effect on the propagation of gravity waves to build reliable seismic diagnostic tools, especially for fast rotating stars, where perturbative treatments of rotation fail. Generalising a previous work done in the case of uniform rotation, we derived a local dispersion relation for gravity waves in a differentially rotating star, taking the full effect of rotation (both Coriolis and centrifugal accelerations) into account. Then we modelled the propagation of axisymmetric waves as the propagation of rays. This allowed us to efficiently probe the properties of the waves in various regimes of differential rotation., Comment: 4 pages, 5 figures, to appear in the proceedings of the PHOST conference in honour of Pr. Shibahashi

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., Comment: 27 pages, 11 figures, Space Weather Special Issue

Context. Differential rotation has a strong influence on stellar internal dynamics and evolution, notably by triggering hydrodynamical instabilities, by interacting with the magnetic field, and more generally by inducing transport of angular momentum and chemical elements. Moreover, it modifies the way waves propagate in stellar interiors and thus the frequency spectrum of these waves, the regions they probe, and the transport they generate. Aims. We investigate the impact of a general differential rotation (both in radius and latitude) on the propagation of axisymmetric gravito-inertial waves. Methods. We use a small-wavelength approximation to obtain a local dispersion relation for these waves. We then describe the propagation of waves thanks to a ray model that follows a Hamiltonian formalism. Finally, we numerically probe the properties of these gravito-inertial rays for different regimes of radial and latitudinal differential rotation. Results. We derive a local dispersion relation that includes the effect of a general differential rotation. Subsequently, considering a polytropic stellar model, we observe that differential rotation allows for a large variety of resonant cavities that can be probed by gravito-inertial waves. We identify that for some regimes of frequency and differential rotation, the properties of gravito-inertial rays are similar to those found in the uniformly rotating case. Furthermore, we also find new regimes specific to differential rotation, where the dynamics of rays is chaotic. Conclusions. As a consequence, we expect modes to follow the same trend. Some parts of oscillation spectra corresponding to regimes similar to those of the uniformly rotating case would exhibit regular patterns, while parts corresponding to the new regimes would be mostly constituted of chaotic modes with a spectrum rather characterised by a generic statistical distribution., Comment: 17 pages, 12 figures, accepted for publication in A&A

The dynamics of the solar wind depends intrinsically on the structure of the global solar magnetic field, which undergoes fundamental changes over the 11-yr solar cycle. For instance, the wind terminal velocity is thought to be anti-correlated with the expansion factor, a measure of how the magnetic field varies with height in the solar corona, usually computed at a fixed height ($\approx 2.5 R_{\odot}$, the source surface radius which approximates the distance at which all magnetic field lines become open). However, the magnetic field expansion affects the solar wind in a more detailed way, its influence on the solar wind properties remaining significant well beyond the source surface. We demonstrate this using 3D global MHD simulations of the solar corona, constrained by surface magnetograms over half a solar cycle (1989-2001). A self-consistent expansion beyond the solar wind critical point (even up to $10R_{\odot}$) makes our model comply with observed characteristics of the solar wind, namely, that the radial magnetic field intensity becomes latitude independent at some distance from the Sun, and that the mass flux is mostly independent of the terminal wind speed. We also show that near activity minimum, the expansion in the higher corona has more influence on the wind speed than the expansion below $2.5 R_{\odot}$., Comment: 11 pages, 10 figures, accepted for publication in the Astrophysical Journal

We present in this work the development of a solar data assimilation method based on an axisymmetric mean field dynamo model and magnetic surface data, our mid-term goal is to predict the solar quasi cyclic activity. Here we focus on the ability of our algorithm to constrain the deep meridional circulation of the Sun based on solar magnetic observations. To that end, we develop a variational data assimilation technique. Within a given assimilation window, the assimilation procedure minimizes the differences between data and the forecast from the model, by finding an optimal meridional circulation in the convection zone, and an optimal initial magnetic field, via a quasi-Newton algorithm. We demonstrate the capability of the technique to estimate the meridional flow by a closed-loop experiment involving 40 years of synthetic, solar-like data. By assimilating the synthetic magnetic proxies annually, we are able to reconstruct a (stochastic) time-varying meridional circulation which is also slightly equatorially asymmetric. We show that the method is robust in estimating a flow whose level of fluctuation can reach 30% about the average, and that the horizon of predictive capability of the method is of the order of 1 cycle length., Comment: Accepted for publication in ApJ. 23 figures, 20 pages. 1 Institut de Physique du Globe de Paris, 75005 Paris, France, 2 AIM/SAp CEA-Saclay, 91191 Gif-Sur-Yvette, France, 3 Institut de Recherche en Astrophysique et Plan\'{e}tologie, 31028 Toulouse, France, 4 Laboratoire de m\'{e}t\'{e}orologie dynamique, Ecole Normale Sup\'{e}rieure, Paris, France

We show how magnetic observations of the Sun can be used in conjunction with an axisymmetric flux-transport solar dynamo model in order to estimate the large-scale meridional circulation throughout the convection zone. Our innovative approach rests on variational data assimilation, whereby the distance between predictions and observations (measured by an objective function) is iteratively minimized by means of an optimization algorithm seeking the meridional flow which best accounts for the data. The minimization is performed using a quasi-Newton technique, which requires the knowledge of the sensitivity of the objective function to the meridional flow. That sensitivity is efficiently computed via the integration of the adjoint flux-transport dynamo model. Closed-loop (also known as twin) experiments using synthetic data demonstrate the validity and accuracy of this technique, for a variety of meridional flow configurations, ranging from unicellular and equatorially symmetric to multicellular and equatorially asymmetric. In this well-controlled synthetic context, we perform a systematic study of the behavior of our variational approach under different observational configurations, by varying their spatial density, temporal density, noise level, as well as the width of the assimilation window. We find that the method is remarkably robust, leading in most cases to a recovery of the true meridional flow to within better than 1%. These encouraging results are a first step towards using this technique to i) better constrain the physical processes occurring inside the Sun and ii) better predict solar activity on decadal time scales., Comment: Published in ApJ 814, 151. 23 pages, 16 figures. 1 AIM/SAp Paris-Saclay, F-91191 Gif-Sur-Yvette, France, 2 Institut de Physique du Globe de Paris, F-75005 Paris, France, 3 Universit\'{e} de Toulouse, Institut de Recherche en Astrophysique et Plan\'{e}tologie, F-31028 Toulouse, France, 4 Laboratoire de m\'{e}t\'{e}orologie dynamique, Ecole Normale Sup\'{e}rieure, Paris, France

During stellar evolution, especially in the PMS, stellar structure and rotation evolve significantly causing major changes in the dynamics and global flows of the star. We wish to assess the consequences of these changes on stellar dynamo, internal magnetic field topology and activity level. To do so, we have performed a series of 3D HD and MHD simulations with the ASH code. We choose five different models characterized by the radius of their radiative zone following an evolutionary track computed by a 1D stellar evolution code. These models characterized stellar evolution from 1 Myr to 50 Myr. By introducing a seed magnetic field in the fully convective model and spreading its evolved state through all four remaining cases, we observe systematic variations in the dynamical properties and magnetic field amplitude and topology of the models. The five MHD simulations develop strong dynamo field that can reach equipartition state between the kinetic and magnetic energy and even super- equipartition levels in the faster rotating cases. We find that the magnetic field amplitude increases as it evolves toward the ZAMS. Moreover the magnetic field topology becomes more complex, with a decreasing axisymmetric component and a non-axisymmetric one becoming predominant. The dipolar components decrease as the rotation rate and the size of the radiative core increase. The magnetic fields possess a mixed poloidal-toroidal topology with no obvious dominant component. Moreover the relaxation of the vestige dynamo magnetic field within the radiative core is found to satisfy MHD stability criteria. Hence it does not experience a global reconfiguration but slowly relaxes by retaining its mixed stable poloidal-toroidal topology., Comment: 26 pages, 23 figures, 10 tables

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., Comment: Published in ApJ, 836, 192. 22 Colour figures, 28 pages. $^1$ AIM/SAp, CEA-Saclay, 91191 Gif-sur-Yvette, France, $^2$ Physics and Astronomy, Univ. of Exeter, Exeter, UK, $^3$ Astronomy Dept., Univ. of Montreal, Montreal, Canada, $^4$ LASP and Dept. of Astrophysical & Planetary Sciences, Univ. of Colorado, Boulder, USA, $^5$ JILA, Univ. of Colorado, Boulder, USA

This paper provides a brief look at dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Two simple models are examined, where one that assumes magnetostrophy and another that includes the effects of inertia. These models are then compared to a suite of convective dynamo simulations of the convective core of a main-sequence B-type star and applied to its later evolutionary stages., Comment: 4 pages, 2 figures, IAU Symposium 329: The lives and death-throes of massive stars

This paper provides a brief overview of dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Three basic approaches are adopted to explore these scaling relationships, with a first look at two simple models: one assuming magnetostrophy and another that includes the effects of inertia. Next, a third scaling relationship is derived that utilizes the assumptions that the dynamo possesses two integral spatial scales and that it is driven by the balance of buoyancy work and ohmic dissipation as studied in Davidson 2013. The results of which are then compared to a suite of convective dynamo simulations that possess a fully convective domain with a weak density stratification and that captured the behavior of the resulting dynamo for a range of convective Rossby numbers (Augustson et al. 2016)., Comment: 7 pages, 3 figures, in the Proceedings of the 19th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun

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., Comment: 15 pages, 9 figures, published in the Astrophysical Journal (update)

The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through 3D global simulations of convective core dynamos operating within a young $\Msun{10}$ B-type star, using the anelastic spherical harmonic (ASH) code. These simulations capture the inner 65\% of this star by radius, encompassing the convective nuclear-burning core (about 23\% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05\% to 16\% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces., Comment: 16 pages, 13 figures

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., Comment: Accepted for publication in the Astrophysical Journal

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., Comment: 6 pages, 3 figures, to appear in the proceedings of the French Astronomy and Astrophysics Society (SF2A) conference

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., Comment: 5 pages, 4 figures, sf2a conference proceedings. appears in Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics (2014)

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)., Comment: 17 pages, 13 figures, accepted for publication in ApJ (10/29/2014)

The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulations relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the $\alpha$-$\Omega$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified., Comment: 24 pages, 13 figures, ApJ, Accepted

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., Comment: 18 pages, 11 figures, 4 tables, accepted for publication in The Astrophysical Journal

Internal gravity waves (hereafter IGWs) are studied for their impact on the angular momentum transport in stellar radiation zones and the information they provide about the structure and dynamics of deep stellar interiors. We here present the first 3D nonlinear numerical simulations of IGWs excitation and propagation in a solar-like star. The aim is to study the behavior of waves in a realistic 3D nonlinear time dependent model of the Sun and to characterize their properties. We compare our results with theoretical and 1D predictions. It allows us to point out the complementarity between theory and simulation and to highlight the convenience but also the limits of the asymptotic and linear theories. We show that a rich spectrum of IGWs is excited by the convection, representing about 0.4% of the total solar luminosity. We study the spatial and temporal properties of this spectrum, the effect of thermal damping and nonlinear interactions between waves. We give quantitative results about the modes frequencies, evolution with time and rotational splitting and we discuss the amplitude of IGWs considering different regimes of parameters. This work points out the importance of high performance simulation for its complementarity with observation and theory. It opens a large field of investigation concerning IGWs propagating nonlinearly in 3D spherical structures. The extension of this work to other types of stars, with different masses, structures and rotation rates will lead to a deeper and more accurate comprehension of IGWs in stars., Comment: 27 pages, 29 figures, accepted for publication in A&A (13/03/2014)

We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles., Comment: Submitted to ApJL 6 pages, 4 figures