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.

Predictions of solar cycle 24 obtained from the advection-dominated and diffusion-dominated kinematic dynamo model are different if we consider only the Babcock-Leighton mechanism as a poloidal field generation source term. Yeates et al (2008) argue that the discrepancy between the results is due to the different memory of solar dynamo for advection- and diffusion-dominated solar convection zones. We aim to investigate the discrepancy between the solar cycle memory obtained from advection-dominated and diffusion-dominated kinematic solar dynamo models. Specifically, we want to investigate whether another poloidal field generation mechanism, namely, Parker's mean-field alpha effect 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 mean-field alpha effect. The addition of a mean-field alpha effect reduces the memory of the Babcock Leighton solar dynamo to within one cycle for both advection-and diffusion-dominated dynamo regimes. Our result provides an alternative resolution to the discrepancy in memory of advection- and diffusion-dominated kinematic dynamo models. We also find a moderate decrease in the correlation between the radial flux at the cycle minima and the toroidal flux at the maxima of the next cycle with the increase of mean-field alpha amplitude. We confirm the conclusions of Yeates et al (2008) and Karak and Nandy (2012) regarding the short, one-cycle memory of the solar cycle. We establish that the mean-field alpha effect impacts the memory of the solar cycle. However, correlation coefficients are highly dependent on the assumed amplitude and fluctuation levels of the mean-field alpha

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.

We present Atacama Large Millimetre Array and Atacama Compact Array observations of the Sunyaev-Zel'dovich effect in the z = 2 galaxy cluster Cl J1449+0856, an X-ray-detected progenitor of typical massive clusters in the present day Universe. While in a cleaned but otherwise untouched 92 GHz map of this cluster, little to no negative signal is visible, careful subtraction of known sub-millimetre emitters in the uv plane reveals a decrement at 5$\sigma$ significance. The total signal is -190$\pm$36 $\mu$Jy, with a peak offset by 5"-9" ($\sim$50 kpc) from both the X-ray centroid and the still-forming brightest cluster galaxy. A comparison of the recovered uv-amplitude profile of the decrement with different pressure models allows us to derive total mass constraints consistent with the $\sim$6$\times$10$^{13}$ M$_{\odot}$ estimated from X-ray data. Moreover, we find no strong evidence for a deviation of the pressure profile with respect to local galaxy clusters, although a slight tension at small-to-intermediate spatial scales suggests a flattened central profile, opposite to what seen in a cool core and possibly an AGN-related effect. This analysis of the lowest mass single SZ detection so far illustrates the importance of interferometers when observing the SZ effect in high-redshift clusters, the cores of which cannot be considered quiescent, such that careful subtraction of galaxy emission is necessary.

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.

Feedback processes from baryons are expected to strongly affect weak-lensing observables of current and future cosmological surveys. In this paper we present a new parametrisation of halo profiles based on gas, stellar, and dark matter density components. This parametrisation is used to modify outputs of gravity-only $N$-body simulations (following the prescription of Schneider and Teyssier [1]) in order to mimic baryonic effects on the matter density field. The resulting baryonic correction model relies on a few well motivated physical parameters and is able to reproduce the redshift zero clustering signal of hydrodynamical simulations at two percent accuracy below $k\sim10$ h/Mpc. A detailed study of the baryon suppression effects on the matter power spectrum and the weak lensing shear correlation reveals that the signal is dominated by two parameters describing the slope of the gas profile in haloes and the maximum radius of gas ejection. We show that these parameters can be constrained with the observed gas fraction of galaxy groups and clusters from X-ray data. Based on these observations we predict a beyond percent effect on the power spectrum above $k=0.2-1.0$ h/Mpc with a maximum suppression of 15-25 percent around $k\sim 10$ h/Mpc. As a result, the weak lensing angular shear power spectrum is suppressed by 15-25 percent at scales beyond $\ell\sim 100-600$ and the shear correlations $\xi_{+}$ and $\xi_{-}$ are affected at the 10-25 percent level below 5 and 50 arc-minutes, respectively. The relatively large uncertainties of these predictions are a result of the poorly known hydrostatic mass bias of current X-ray observations as well as the generic difficulty to observe the low density gas outside of haloes.

Rotational scaling relationships are examined for the degree of equipartition between magnetic and kinetic energies in stellar convection zones. These scaling relationships are approached from two paradigms, with first a glance at scaling relationship built upon an energy-balance argument and second a look at a force-based scaling. The latter implies a transition between a nearly-constant inertial scaling when in the asymptotic limit of minimal diffusion and magnetostrophy, whereas the former implies a weaker scaling with convective Rossby number. Both scaling relationships are then compared to a suite of 3D convective dynamo simulations with a wide variety of domain geometries, stratifications, and range of convective Rossby numbers.

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.

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.

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.

With the discovery over the last two decades of a large diversity of exoplanetary systems, it is now of prime importance to characterize star-planet interactions and how such systems evolve. We address this question by studying systems formed by a solar-like star and a close-in planet. We focus on the stellar wind spinning down the star along its main sequence phase and tidal interaction causing orbital evolution of the systems. Despite recent significant advances in these fields, all current models use parametric descriptions to study at least one of these effects. Our objective is to introduce simultaneously ab-initio prescriptions of the tidal and braking torques, so as to improve our understanding of the underlying physics. We develop a 1D numerical model of coplanar circular star-planet systems taking into account stellar structural changes, wind braking and tidal interaction and implement it in a code called ESPEM. We follow the secular evolution of the stellar rotation assuming a bi-layer internal structure, and of the semi-major axis of the orbit. After comparing our predictions to recent observations and models, we perform tests to emphasize the contribution of ab-initio prescriptions. Our secular model of stellar wind braking reproduces well the recent observations of stellar rotation in open clusters. Our results show that a planet can affect the rotation of its host star and that the resulting spin-up or spin-down depends on the orbital semi-major axis and on the joint influence of magnetic and tidal effects. The ab-initio prescription for tidal dissipation that we used predicts fast outward migration of massive planet orbiting fast-rotating young stars. Finally, we provide the reader with a criterion based on the system's characteristics that allows us to assess whether or not the planet will undergo orbital decay due to tidal interaction.

Accurate cosmology from upcoming weak lensing surveys relies on knowledge of the total matter power spectrum at percent level at scales $k < 10$ $h$/Mpc, for which modelling the impact of baryonic physics is crucial. We compare measurements of the total matter power spectrum from the Horizon cosmological hydrodynamical simulations: a dark matter-only run, one with full baryonic physics, and another lacking Active Galactic Nuclei (AGN) feedback. Baryons cause a suppression of power at $k\simeq 10$ $h/$Mpc of $<15\%$ at $z=0$, and an enhancement of a factor of a few at smaller scales due to the more efficient cooling and star formation. The results are sensitive to the presence of the highest mass haloes in the simulation and the distribution of dark matter is also impacted up to a few percent. The redshift evolution of the effect is non-monotonic throughout $z=0-5$ due to an interplay between AGN feedback and gas pressure, and the growth of structure. We investigate the effectiveness of an analytic `baryonic correction model' in describing our results. We require a different redshift evolution and propose an alternative fitting function with $4$ free parameters that reproduces our results within $5\%$. Compared to other simulations, we find the impact of baryonic processes on the total matter power spectrum to be smaller at $z=0$. Correspondingly, our results suggest that AGN feedback is not strong enough in the simulation. Total matter power spectra from the Horizon simulations are made publicly available at https://www.horizon-simulation.org/catalogues.html

A well characterised detection pipeline is an important ingredient for X-ray cluster surveys. We present the final development of the XXL Survey pipeline. The pipeline optimally uses X-ray information by combining many overlapping observations of a source when possible, both for its detection and its characterisation. It can robustly detect and characterise several types of X-ray sources: AGNs (point-like), galaxy clusters (extended), galaxy clusters contaminated by a central AGN, and pairs of AGNs close on the sky. We perform a thorough suite of validation tests via realistic simulations of XMM-Newton images and we introduce new selection criteria for various types of sources that will be detected by the survey. We find that the use of overlapping observations allows new clusters to be securely identified that would be missed or less securely identified by using only one observation at the time. We also find that with the new pipeline we can robustly identify clusters with a central AGN that would otherwise have been missed, and we can flag pairs of AGNs close on the sky that might have been mistaken for a cluster.

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 magnetized wind from stars that impact exoplanets should lead to radio emissions. According to the scaling laws derived in the solar system, the radio emission should depend on the stellar wind, interplanetary magnetic field, and topology of the exoplanet magnetosphere. The aim of this study is to calculate the dissipated power and subsequent radio emission from exoplanet magnetospheres with different topologies perturbed by the interplanetary magnetic field and stellar wind, to refine the predictions from scaling laws, and to prepare the interpretation of future radio detections. We use the magnetohydrodynamic (MHD) code PLUTO in spherical coordinates to analyze the total radio emission level resulting from the dissipation of the kinetic and magnetic (Poynting flux) energies inside the exoplanet's magnetospheres. We apply a formalism to infer the detailed contribution in the exoplanet radio emission on the exoplanet's day side and magnetotail. The model is based on Mercury-like conditions, although the study results are extrapolated to exoplanets with stronger magnetic fields, providing the lower bound of the radio emission. The predicted dissipated powers and resulting radio emissions depends critically on the exoplanet magnetosphere topology and interplanetary magnetic field (IMF) orientation. The radio emission on the exoplanet's night and day sides should thus contain information on the exoplanet magnetic field topology. In addition, if the topology of an exoplanet magnetosphere is known, the radio emission measurements can be used as a proxy of the instantaneous dynamic pressure of the stellar wind, IMF orientation, and intensity.

We present a detailed study of the integrated total hydrostatic mass profiles of the five most massive ($M^{\mathrm{SZ}}_{500} > 5 \times 10^{14}$ M$_{\odot}$) galaxy clusters selected at $z\sim1$ via the Sunyaev-Zel'dovich effect. These objects represent an ideal laboratory to test structure formation models where the primary driver is gravity. Optimally exploiting spatially-resolved spectroscopic information from XMM-Newton and Chandra observations, we used both parametric (forward, backward) and non-parametric methods to recover the mass profiles, finding that the results are extremely robust when density and temperature measurements are both available. Our X-ray masses at $R_{500}$ are higher than the weak lensing masses obtained from the Hubble Space Telescope (HST), with a mean ratio of $1.39^{+0.47}_{-0.35}$. This offset goes in the opposite direction to that expected in a scenario where the hydrostatic method yields a biased, underestimated, mass. We investigated halo shape parameters such as sparsity and concentration, and compared to local X-ray selected clusters, finding hints for evolution in the central regions (or for selection effects). The total baryonic content is in agreement with the cosmic value at $R_{500}$. Comparison with numerical simulations shows that the mass distribution and concentration are in line with expectations. These results illustrate the power of X-ray observations to probe the statistical properties of the gas and total mass profiles in this high-mass, high-redshift regime.

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.

Weak gravitational lensing depends on the integrated mass along the line of sight. Baryons contribute to the mass distribution of galaxy clusters and the resulting mass estimates from lensing analysis. We use the cosmo-OWLS suite of hydrodynamic simulations to investigate the impact of baryonic processes on the bias and scatter of weak lensing mass estimates of clusters. These estimates are obtained by fitting NFW profiles to mock data using MCMC techniques. In particular, we examine the difference in estimates between dark matter-only runs and those including various prescriptions for baryonic physics. We find no significant difference in the mass bias when baryonic physics is included, though the overall mass estimates are suppressed when feedback from AGN is included. For lowest-mass systems for which a reliable mass can be obtained ($M_{200} \approx 2 \times 10^{14}$ $M_{\odot}$), we find a bias of $\approx -10$ per cent. The magnitude of the bias tends to decrease for higher mass clusters, consistent with no bias for the most massive clusters which have masses comparable to those found in the CLASH and HFF samples. For the lowest mass clusters, the mass bias is particularly sensitive to the fit radii and the limits placed on the concentration prior, rendering reliable mass estimates difficult. The scatter in mass estimates between the dark matter-only and the various baryonic runs is less than between different projections of individual clusters, highlighting the importance of triaxiality.

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.

This study aims at investigating the ability of strong interactions between magnetic field concentrations during their rise through the convection zone to produce complex active regions at the solar surface. To do so, we perform numerical simulations of buoyant magnetic structures evolving and interacting in a model solar convection zone. We first produce a 3D model of rotating convection and then introduce idealized magnetic structures close to the bottom of the computational domain. These structures possess a certain degree of field line twist and they are made buoyant on a particular extension in longitude. The resulting twisted $\Omega$-loops will thus evolve inside a spherical convective shell possessing large-scale mean flows. We present results on the interaction between two such loops with various initial parameters (mainly buoyancy and twist) and on the complexity of the emerging magnetic field. In agreement with analytical predictions, we find that if the loops are introduced with opposite handedness and same axial field direction or same handedness but opposite axial field, they bounce against each other. The emerging region is then constituted of two separated bipolar structures. On the contrary, if the loops are introduced with the same direction of axial and peripheral magnetic fields and if sufficiently close, they merge while rising. This more interesting case produces complex magnetic structures, with a high degree of non-neutralized currents, especially when the convective motions act significantly on the magnetic field. This indicates that those interactions could be good candidates to produce eruptive events like flares or CMEs.

Observations of Sun-like stars over the last half-century have improved our understanding of how magnetic dynamos, like that responsible for the 11-year solar cycle, change with rotation, mass and age. Here we show for the first time how metallicity can affect a stellar dynamo. Using the most complete set of observations of a stellar cycle ever obtained for a Sun-like star, we show how the solar analog HD 173701 exhibits solar-like differential rotation and a 7.4-year activity cycle. While the duration of the cycle is comparable to that generated by the solar dynamo, the amplitude of the brightness variability is substantially stronger. The only significant difference between HD 173701 and the Sun is its metallicity, which is twice the solar value. Therefore, this provides a unique opportunity to study the effect of the higher metallicity on the dynamo acting in this star and to obtain a comprehensive understanding of the physical mechanisms responsible for the observed photometric variability. The observations can be explained by the higher metallicity of the star, which is predicted to foster a deeper outer convection zone and a higher facular contrast, resulting in stronger variability.

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}$.

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.

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.

[Abridged] Typical disc galaxies forming in a LambdaCDM cosmology encounter a violent environment, where they often experience mergers with massive satellites. The fact that disc galaxies are ubiquitous in the local Universe suggests that a quiescent history is not necessary for their formation. Modern cosmological simulations can now obtain relatively realistic populations of disc galaxies, but it still remains to be clarified how discs manage to survive massive mergers. Here we use a suite of high-resolution hydrodynamical simulations set in a LambdaCDM cosmology to elucidate the fate of discs encountering massive mergers. We extract a sample of approximately 100 disc galaxies and follow the changes in their post-merger morphologies, as tracked by their disc-to-total ratios (D/T). We also examine the relations between their present-day morphology, assembly history and gas fractions. We find that approximately half of present-day disc galaxies underwent at least one merger with a satellite of total mass exceeding the host system's stellar mass, a third had mergers with satellites of mass exceeding 3 times the host's stellar mass, and approximately one-sixth had mergers with satellites of mass exceeding 10 times of the host's stellar mass. These mergers lead to a sharp, but often temporary, decrease in the D/T of the hosts, implying that discs are usually disrupted but then quickly re-grow. To do so, high cold gas fractions are required post-merger, as well as a relatively quiescent recent history (over a few Gyrs before z=0). Our results show that discs can form via diverse merger pathways and that quiescent histories are not the dominant mode of disc formation.

We investigate the evolution of the dark matter density profiles of the most massive galaxy clusters in the Universe. Using a `zoom-in' procedure on a large suite of cosmological simulations of total comoving volume of $3\,(h^{-1}\,\rm Gpc)^3$, we study the 25 most massive clusters in four redshift slices from $z\sim 1$ to the present. The minimum mass is $M_{500} > 5.5 \times 10^{14}$ M$_{\odot}$ at $z=1$. Each system has more than two million particles within $r_{500}$. Once scaled to the critical density at each redshift, the dark matter profiles within $r_{500}$ are strikingly similar from $z\sim1$ to the present day, exhibiting a low dispersion of 0.15 dex, and showing little evolution with redshift in the radial logarithmic slope and scatter. They have the running power law shape typical of the NFW-type profiles, and their inner structure, resolved to $3.8\,h^{-1}$ comoving kpc at $z=1$, shows no signs of converging to an asymptotic slope. Our results suggest that this type of profile is already in place at $z>1$ in the highest-mass haloes in the Universe, and that it remains exceptionally robust to merging activity.

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.

Modern cosmological simulations rely heavily on feedback from active galactic nuclei (AGN) in order to stave off overcooling in massive galaxies and galaxy groups and clusters. An important independent test is whether or not the simulations capture the broad demographics of the observed AGN population. Here, we have used the cosmo-OWLS suite of cosmological hydrodynamical simulations to produce realistic synthetic catalogs of X-ray AGN out to $z$=3, with the aim of comparing the catalogs to the observed X-ray AGN population in the XXL survey and other recent surveys. We focused on the unabsorbed X-ray luminosity function (XLF), the Eddington ratio distribution, the black hole mass function, and the projected clustering of X-ray AGN. To compute the unabsorbed XLF of the simulated AGN, we used recent empirically-determined bolometric corrections. We show that the simulated AGN sample accurately reproduces the observed XLF over 3 orders of magnitude in X-ray luminosity in all redshift bins. To compare to the observed Eddington ratio distribution and the clustering of AGN, we produced detailed 'XMM-Newton-detected' catalogs of the simulated AGN. This requires the production of synthetic X-ray images extracted from light cones of the simulations that fold in the relevant instrumental effects of XMM-Newton. We apply a luminosity- and redshift-dependent obscuration function for the AGN and employ the same AGN detection algorithm as used for the real XXL survey. We demonstrate that the detected population of simulated AGN reproduces the observed Eddington ratio distribution and projected clustering from XXL quite well. We conclude that the simulations have a broadly realistic population of AGN and that our synthetic X-ray AGN catalogs should be useful for interpreting additional trends and as a helpful tool for quantifying AGN contamination in galaxy group and cluster X-ray surveys.

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.

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.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

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.

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.

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

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations to investigate the scatter and evolution of the global hot gas properties of large simulated populations of galaxy groups and clusters. Our aim is to compare the predictions of different physical models and to explore the extent to which commonly-adopted assumptions in observational analyses (e.g. self-similar evolution) are violated. We examine the relations between (true) halo mass and the X-ray temperature, X-ray luminosity, gas mass, Sunyaev-Zel'dovich (SZ) flux, the X-ray analogue of the SZ flux ($Y_X$) and the hydrostatic mass. For the most realistic models, which include AGN feedback, the slopes of the various mass-observable relations deviate substantially from the self-similar ones, particularly at late times and for low-mass clusters. The amplitude of the mass-temperature relation shows negative evolution with respect to the self-similar prediction (i.e. slower than the prediction) for all models, driven by an increase in non-thermal pressure support at higher redshifts. The AGN models predict strong positive evolution of the gas mass fractions at low halo masses. The SZ flux and $Y_X$ show positive evolution with respect to self-similarity at low mass but negative evolution at high mass. The scatter about the relations is well approximated by log-normal distributions, with widths that depend mildly on halo mass. The scatter decreases significantly with increasing redshift. The exception is the hydrostatic mass-halo mass relation, for which the scatter increases with redshift. Finally, we discuss the relative merits of various hot gas-based mass proxies.

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.

The evolution of the large-scale distribution of matter is sensitive to a variety of fundamental parameters that characterise the dark matter, dark energy, and other aspects of our cosmological framework. Since the majority of the mass density is in the form of dark matter that cannot be directly observed, to do cosmology with large-scale structure one must use observable (baryonic) quantities that trace the underlying matter distribution in a (hopefully) predictable way. However, recent numerical studies have demonstrated that the mapping between observable and total mass, as well as the total mass itself, are sensitive to unresolved feedback processes associated with galaxy formation, motivating explicit calibration of the feedback efficiencies. Here we construct a new suite of large-volume cosmological hydrodynamical simulations (called BAHAMAS, for BAryons and HAloes of MAssive Systems) where subgrid models of stellar and Active Galactic Nucleus (AGN) feedback have been calibrated to reproduce the present-day galaxy stellar mass function and the hot gas mass fractions of groups and clusters in order to ensure the effects of feedback on the overall matter distribution are broadly correct. We show that the calibrated simulations reproduce an unprecedentedly wide range of properties of massive systems, including the various observed mappings between galaxies, hot gas, total mass, and black holes, and represent a significant advance in our ability to mitigate the primary systematic uncertainty in most present large-scale structure tests.

We report the discovery of a remarkable concentration of massive galaxies with extended X-ray emission at $z_{spec} = 2.506$, which contains 11 massive ($M_{*} \gtrsim 10^{11} M_{\odot}$) galaxies in the central 80kpc region (11.6$\sigma$ overdensity). We have spectroscopically confirmed 17 member galaxies with 11 from CO and the remaining ones from $H\alpha$. The X-ray luminosity, stellar mass content and velocity dispersion all point to a collapsed, cluster-sized dark matter halo with mass $M_{200c} = 10^{13.9\pm0.2} M_{\odot}$, making it the most distant X-ray-detected cluster known to date. Unlike other clusters discovered so far, this structure is dominated by star-forming galaxies (SFGs) in the core with only 2 out of the 11 massive galaxies classified as quiescent. The star formation rate (SFR) in the 80kpc core reaches $\sim$3400 $M_{\odot}$ yr$^{-1}$ with a gas depletion time of $\sim 200$ Myr, suggesting that we caught this cluster in rapid build-up of a dense core. The high SFR is driven by both a high abundance of SFGs and a higher starburst fraction ($\sim25\%$, compared to 3\%-5\% in the field). The presence of both a collapsed, cluster-sized halo and a predominant population of massive SFGs suggests that this structure could represent an important transition phase between protoclusters and mature clusters. It provides evidence that the main phase of massive galaxy passivization will take place after galaxies accrete onto the cluster, providing new insights into massive cluster formation at early epochs. The large integrated stellar mass at such high redshift challenges our understanding of massive cluster formation.

Context: Active stars possess magnetized wind that has a direct impact on planets that can lead to radio emission. Mercury is a good test case to study the effect of the solar wind and interplanetary magnetic field on radio emission driven in the planet magnetosphere. Such studies could be used as proxies to characterize the magnetic field topology and intensity of exoplanets. Aims: The aim of this study is to quantify the radio emission in the Hermean magnetosphere. Methods: We use the MHD code PLUTO in spherical coordinates with an axisymmetric multipolar expansion for the Hermean magnetic field, to analyze the effect of the interplanetary magnetic field (IMF) orientation and intensity, as well as the hydrodynamic parameters of the solar wind (velocity, density and temperature), on the net power dissipated on the Hermean day and night side. We apply the formalism derived by Zarka [2001, 2007] to infer the radio emission level from the net dissipated power. We perform a set of simulations with different hydrodynamic parameters of the solar wind, IMF orientations and intensities, that allow us to calculate the dissipated power distribution and infer the existence of radio emission hot spots on the planet day side, and to calculate the integrated radio emission of the Hermean magnetosphere. Results: The obtained radio emission distribution of dissipated power is determined by the IMF orientation (associated with the reconnection regions in the magnetosphere), although the radio emission strength is dependent on the IMF intensity and solar wind hydro parameters. The calculated total radio emission level is in agreement with the one estimated in [Zarka, 2001], between $5 \times 10^{5}$ and $2 \times 10^{6}$ W.

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.

We present the discovery of a giant $\gtrsim$100~kpc Ly$\alpha$ nebula detected in the core of the X-ray emitting cluster CL~J1449+0856 at $z=1.99$ through Keck/LRIS narrow-band imaging. This detection extends the known relation between Ly$\alpha$ nebulae and overdense regions of the Universe to the dense core of a $5-7\times10^{13}$ M$_{\odot}$ cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members, while cooling from the X-ray phase and cosmological cold flows are disfavored primarily because of the high Ly$\alpha$ to X-ray luminosity ratio ($L_{\mathrm{Ly\alpha}}/L_{\mathrm{X}} \approx0.3$, $\gtrsim10-1000\times$ higher than in local cool-core clusters) and by current modeling. Given the physical conditions of the Ly$\alpha$-emitting gas and the possible interplay with the X-ray phase, we argue that the Ly$\alpha$ nebula would be short-lived ($\lesssim10$ Myr) if not continuously replenished with cold gas at a rate of $\gtrsim1000$ M$_{\odot}$ yr$^{-1}$. We investigate the possibility that cluster galaxies supply the required gas through outflows and we show that their total mass outflow rate matches the replenishment necessary to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium, as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of $\thickapprox2$ keV per particle in the intracluster medium over a $2$ Gyr interval. In our baseline calculation AGN provide up to $85$% of the injected energy and 2/3 of the mass, while the rest is supplied by supernovae-driven winds.

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.

The solar wind speed at 1 AU shows variations in latitude and in time which reflect the evolution of the global background magnetic field during the activity cycle. It is commonly accepted that the terminal wind speed in a magnetic flux-tube is anti-correlated with its expansion ratio, which motivated the definition of widely-used semi-empirical scaling laws relating one to the other. In practice, such scaling laws require ad-hoc corrections. A predictive law based solely on physical principles is still missing. We test whether the flux-tube expansion is the controlling factor of the wind speed at all phases of the cycle and at all latitudes using a very large sample of wind-carrying open magnetic flux-tubes. We furthermore search for additional physical parameters based on the geometry of the coronal magnetic field which have an influence on the terminal wind flow speed. We use MHD simulations of the corona and wind coupled to a dynamo model to provide a large statistical ensemble of open flux-tubes which we analyse conjointly in order to identify relations of dependence between the wind speed and geometrical parameters of the flux-tubes which are valid globally (for all latitudes and moments of the cycle). Our study confirms that the terminal speed of the solar wind depends very strongly on the geometry of the open magnetic flux-tubes through which it flows. The total flux-tube expansion is more clearly anti-correlated with the wind speed for fast rather than for slow wind flows, and effectively controls the locations of these flows during solar minima. Overall, the actual asymptotic wind speeds attained are also strongly dependent on field-line inclination and magnetic field amplitude at the foot-points. We suggest ways of including these parameters on future predictive scaling-laws for the solar wind speed.

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.

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations, which includes different galactic feedback models, to predict the cross-correlation signal between weak gravitational lensing and the thermal Sunyaev-Zeldovich (tSZ) $y$-parameter. The predictions are compared to the recent detection reported by van Waerbeke and collaborators. The simulations reproduce the weak lensing-tSZ cross-correlation, $\xi_{y\kappa}(\theta)$, well. The uncertainty arising from different possible feedback models appears to be important on small scales only ($\theta \lesssim 10$ arcmin), while the amplitude of the correlation on all scales is sensitive to cosmological parameters that control the growth rate of structure (such as $\sigma_8$, $\Omega_m$ and $\Omega_b$). This study confirms our previous claim (in Ma et al.) that a significant proportion of the signal originates from the diffuse gas component in low-mass ($M_{\rm{halo}} \lesssim 10^{14} M_{\odot}$) clusters as well as from the region beyond the virial radius. We estimate that approximately 20$\%$ of the detected signal comes from low-mass clusters, which corresponds to about 30$\%$ of the baryon density of the Universe. The simulations also suggest that more than half of the baryons in the Universe are in the form of diffuse gas outside halos ($\gtrsim 5$ times the virial radius) which is not hot or dense enough to produce a significant tSZ signal or be observed by X-ray experiments. Finally, we show that future high-resolution tSZ-lensing cross-correlation observations will serve as a powerful tool for discriminating between different galactic feedback models.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

Traditionally, galaxy clusters have been expected to retain all the material accreted since their formation epoch. For this reason, their matter content should be representative of the Universe as a whole, and thus their baryon fraction should be close to the Universal baryon fraction. We make use of the sample of the 100 brightest galaxy clusters discovered in the XXL Survey to investigate the fraction of baryons in the form of hot gas and stars in the cluster population. We measure the gas masses of the detected halos and use a mass--temperature relation directly calibrated using weak-lensing measurements for a subset of XXL clusters to estimate the halo mass. We find that the weak-lensing calibrated gas fraction of XXL-100-GC clusters is substantially lower than was found in previous studies using hydrostatic masses. Our best-fit relation between gas fraction and mass reads $f_{\rm gas,500}=0.055_{-0.006}^{+0.007}\left(M_{\rm 500}/10^{14}M_\odot\right)^{0.21_{-0.10}^{+0.11}}$. The baryon budget of galaxy clusters therefore falls short of the Universal baryon fraction by about a factor of two at $r_{\rm 500}$. Our measurements require a hydrostatic bias $1-b=M_X/M_{\rm WL}=0.72_{-0.07}^{+0.08}$ to match the gas fraction obtained using lensing and hydrostatic equilibrium. Comparing our gas fraction measurements with the expectations from numerical simulations, our results favour an extreme feedback scheme in which a significant fraction of the baryons are expelled from the cores of halos. This model is, however, in contrast with the thermodynamical properties of observed halos, which might suggest that weak-lensing masses are overestimated. We note that a mass bias $1-b=0.58$ as required to reconcile Planck CMB and cluster counts should translate into an even lower baryon fraction, which poses a major challenge to our current understanding of galaxy clusters. [Abridged]

We study the magnetic flux carried by pores located outside active regions with sunspots and investigate their possible contribution to the reversal of the global magnetic field of the Sun. We find that they contain a total flux of comparable amplitude to the total magnetic flux contained in polar caps. The pores located at distances of 40--100~Mm from the closest active region have systematically the correct sign to contribute to the polar cap reversal. These pores can predominantly be found in bipolar magnetic regions. We propose that during grand minima of solar activity, such a systematic polarity trend, akin to a weak magnetic (Babcock-Leighton-like) source term could still be operating but was missed by the contemporary observers due to the limited resolving power of their telescopes.

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.

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.

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.

The thermal Sunyaev-Zel'dovich (tSZ) effect offers a means of probing the hot gas in and around massive galaxies and galaxy groups and clusters, which is thought to constitute a large fraction of the baryon content of the Universe. The Planck collaboration recently performed a stacking analysis of a large sample of `locally brightest galaxies' (LBGs) and, surprisingly, inferred an approximately self-similar relation between the tSZ flux and halo mass. At face value, this implies that the hot gas mass fraction is independent of halo mass, a result which is in apparent conflict with resolved X-ray observations. We test the robustness of the inferred trend using synthetic tSZ maps generated from cosmological hydrodynamical simulations and using the same tools and assumptions applied in the Planck study. We show that, while the detection and the estimate of the `total' flux (within $5 r_{500}$) is reasonably robust, the inferred flux originating from within $r_{500}$ (i.e. the limiting radius to which X-ray observations typically probe) is highly sensitive to the assumed pressure distribution of the gas. Using our most realistic simulations with AGN feedback, that reproduce a wide variety of X-ray and optical properties of groups and clusters, we estimate that the derived tSZ flux within $r_{500}$ is biased high by up to to an order of magnitude for haloes with masses $M_{500} \sim 10^{13}$ M$_{\odot}$. Moreover, we show that the AGN simulations are consistent with the total tSZ flux-mass relation observed with Planck, whereas a self-similar model is ruled out.

We discuss how recent advances in observations, theory and numerical simulations have allowed the stellar community to progress in its understanding of stellar convection, rotation and magnetism and to assess the degree to which the Sun and other stars share similar dynamical properties. Ensemble asteroseismology has become a reality with the advent of large time domain studies, especially from space missions. This new capability has provided improved constraints on stellar rotation and activity, over and above that obtained via traditional techniques such as spectropolarimetry or CaII H&K observations. New data and surveys covering large mass and age ranges have provided a wide parameter space to confront theories of stellar magnetism. These new empirical databases are complemented by theoretical advances and improved multi-D simulations of stellar dynamos. We trace these pathways through which a lucid and more detailed picture of magnetohydrodynamics of solar-like stars is beginning to emerge and discuss future prospects.

To better understand the observed distributions of rotation rate and magnetic activity of sun-like and low-mass stars, we derive a physically motivated scaling for the dependence of the stellar-wind torque on Rossby number. The torque also contains an empirically-derived scaling with stellar mass (and radius), which provides new insight into the mass-dependence of stellar magnetic and wind properties. We demonstrate that this new formulation explains why the lowest mass stars are observed to maintain rapid rotation for much longer than solar-mass stars, and simultaneously, why older populations exhibit a sequence of slowly rotating stars, in which the low-mass stars rotate more slowly than solar-mass stars. The model also reproduces some previously unexplained features in the period-mass diagram for the Kepler field, notably: the particular shape of the "upper envelope" of the distribution, suggesting that ~95% of Kepler field stars with measured rotation periods are younger than ~4 Gyr; and the shape of the "lower envelope," corresponding to the location where stars transition between magnetically saturated and unsaturated regimes.

Solar flares are associated with intense soft X-ray emission generated by the hot flaring plasma. Kink unstable twisted flux-ropes provide a source of magnetic energy which can be released impulsively and account for the flare plasma heating. We compute the temporal evolution of the thermal X-ray emission in kink-unstable coronal loops using MHD simulations and discuss the results of with respect to solar flare observations. The model consists of a highly twisted loop embedded in a region of uniform and untwisted coronal magnetic field. We let the kink instability develop, compute the evolution of the plasma properties in the loop (density, temperature) without accounting for mass exchange with the chromosphere. We then deduce the X-ray emission properties of the plasma during the whole flaring episode. During the initial phase of the instability plasma heating is mostly adiabatic. Ohmic diffusion takes over as the instability saturates, leading to strong and impulsive heating (> 20 MK), to a quick enhancement of X-ray emission and to the hardening of the thermal X-ray spectrum. The temperature distribution of the plasma becomes broad, with the emission measure depending strongly on temperature. Significant emission measures arise for plasma at temperatures T > 9 MK. The magnetic flux-rope then relaxes progressively towards a lower energy state as it reconnects with the background flux. The loop plasma suffers smaller sporadic heating events but cools down conductively. The total thermal X-ray emission slowly fades away during this phase, and the high temperature component of emission measure distribution converges to the power-law distribution $EM\propto T^{-4.2}$. The amount of twist deduced directly from the X-ray emission patterns is considerably lower than the maximum magnetic twist in the simulated flux-ropes.

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.

The Kepler mission has collected light curves for almost 4 years. The excellent quality of these data has allowed us to probe the structure and the dynamics of the stars using asteroseismology. With the length of data available, we can start to look for magnetic activity cycles. The Kepler data obtained for the F star, KIC 12009504, shows a rotation period of 9.5 days and additional variability that could be due to the magnetic activity of the star. Here we present recent and preliminary 3D global-scale dynamo simulations of this star with the ASH and STELEM codes, capturing a substantial portion of the convection and the stable radiation zone below it. These simulations reveal a multi-year activity cycle whose length tentatively depends upon the width of the tachocline present in the simulation. Furthermore, the presence of a magnetic field and the dynamo action taking place in the convection zone appears to help confine the tachocline, but longer simulations will be required to confirm this.

We use cosmological hydrodynamical simulations to investigate how the inclusion of physical processes relevant to galaxy formation (star formation, metal-line cooling, stellar winds, supernovae and feedback from Active Galactic Nuclei, AGN) change the properties of haloes, over four orders of magnitude in mass. We find that gas expulsion and the associated dark matter (DM) expansion induced by supernova-driven winds are important for haloes with masses M200 < 10^13 Msun, lowering their masses by up to 20% relative to a DM-only model. AGN feedback, which is required to prevent overcooling, has a significant impact on halo masses all the way up to cluster scales (M200 ~ 10^15 Msun). Baryonic physics changes the total mass profiles of haloes out to several times the virial radius, a modification that cannot be captured by a change in the halo concentration. The decrease in the total halo mass causes a decrease in the halo mass function of about 20%. This effect can have important consequences for abundance matching technique as well as for most semi-analytic models of galaxy formation. We provide analytic fitting formulae, derived from simulations that reproduce the observed baryon fractions, to correct halo masses and mass functions from DM-only simulations. The effect of baryonic physics (AGN feedback in particular) on cluster number counts is about as large as changing the cosmology from WMAP7 to Planck, even when a moderately high mass limit of M500 ~ 10^14 Msun is adopted. Thus, for precision cosmology the effects of baryons must be accounted for.

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.

We present a new suite of large-volume cosmological hydrodynamical simulations called cosmo-OWLS. They form an extension to the OverWhelmingly Large Simulations (OWLS) project, and have been designed to help improve our understanding of cluster astrophysics and non-linear structure formation, which are now the limiting systematic errors when using clusters as cosmological probes. Starting from identical initial conditions in either the Planck or WMAP7 cosmologies, we systematically vary the most important `sub-grid' physics, including feedback from supernovae and active galactic nuclei (AGN). We compare the properties of the simulated galaxy groups and clusters to a wide range of observational data, such as X-ray luminosity and temperature, gas mass fractions, entropy and density profiles, Sunyaev-Zel'dovich flux, I-band mass-to-light ratio, dominance of the brightest cluster galaxy, and central massive black hole (BH) masses, by producing synthetic observations and mimicking observational analysis techniques. These comparisons demonstrate that some AGN feedback models can produce a realistic population of galaxy groups and clusters, broadly reproducing both the median trend and, for the first time, the scatter in physical properties over approximately two decades in mass ($10^{13} \lesssim M_{500} \lesssim 10^{15}~\textrm{M}_{\odot}$) and 1.5 decades in radius ($0.05 \lesssim r/r_{500} \lesssim 1.5$). However, in other models, the AGN feedback is too violent (even though they reproduce the observed BH scaling relations), implying calibration of the models is required. The production of realistic populations of simulated groups and clusters, as well as models that bracket the observations, opens the door to the creation of synthetic surveys for assisting the astrophysical and cosmological interpretation of cluster surveys, as well as quantifying the impact of selection effects.

(Abridged) The amplitude of the thermal Sunyaev Zel'dovich effect (tSZ) power spectrum is extremely sensitive to the abundance of galaxy clusters and therefore to fundamental cosmological parameters that control their growth, such as sigma_8 and Omega_m. Here we explore the sensitivity of the tSZ power spectrum to important non-gravitational ('sub-grid') physics by employing the cosmo-OWLS suite of large-volume cosmological hydrodynamical simulations, run in both the Planck and WMAP7 best-fit cosmologies. On intermediate and small angular scales (ell > ~1000, or theta < ~10 arcmin), accessible with the South Pole Telescope and the Atacama Cosmology Telescope, the predicted tSZ power spectrum is highly model dependent, with AGN feedback having a particularly large effect. However, at large scales, observable with the Planck telescope, the effects of sub-grid physics are minor. Comparing the simulations with observations, we find a significant amplitude offset on all measured angular scales (including large scales), if the Planck best-fit cosmology is assumed by the simulations. This is shown to be a generic result for all current tSZ models. By contrast, if the WMAP7 cosmology is adopted, there is full consistency with the Planck power spectrum measurements on large scales and agreement at the 2 sigma level with the SPT/ACT measurements at intermediate scales for our fiducial AGN model, which Le Brun et al. (2014) have shown reproduces the 'resolved' properties of the local cluster population remarkably well. These findings strongly suggest that there are significantly fewer massive galaxy clusters than expected for the Planck best-fit cosmology, which is consistent with recent measurements of the tSZ number counts. Our findings therefore pose a significant challenge to the cosmological parameter values preferred (and/or the model adopted) by the Planck primary CMB analyses.

We present preliminary results of a 3D MHD simulation of the convective envelope of the giant star Pollux for which the rotation period and the magnetic field intensity have been measured from spectroscopic and spectropolarimetric observations. This giant is one of the first single giants with a detected magnetic field, and the one with the weakest field so far. Our aim is to understand the development and the action of the dynamo in its extended convective envelope.

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

SPIRou is a near-infrared (nIR) spectropolarimeter / velocimeter proposed as a new-generation instrument for CFHT. SPIRou aims in particular at becoming world-leader on two forefront science topics, (i) the quest for habitable Earth-like planets around very- low-mass stars, and (ii) the study of low-mass star and planet formation in the presence of magnetic fields. In addition to these two main goals, SPIRou will be able to tackle many key programs, from weather patterns on brown dwarf to solar-system planet atmospheres, to dynamo processes in fully-convective bodies and planet habitability. The science programs that SPIRou proposes to tackle are forefront (identified as first priorities by most research agencies worldwide), ambitious (competitive and complementary with science programs carried out on much larger facilities, such as ALMA and JWST) and timely (ideally phased with complementary space missions like TESS and CHEOPS). SPIRou is designed to carry out its science mission with maximum efficiency and optimum precision. More specifically, SPIRou will be able to cover a very wide single-shot nIR spectral domain (0.98-2.35 \mu m) at a resolving power of 73.5K, providing unpolarized and polarized spectra of low-mass stars with a ~15% average throughput and a radial velocity (RV) precision of 1 m/s.

We study the influence of a dynamo magnetic field on the buoyant rise and emergence of twisted magnetic flux-ropes, and their influence on the global external magnetic field. We ran 3D MHD numerical simulations using the ASH code and analysed the dynamical evolution of such buoyant flux-ropes from the bottom of the convection zone until the post-emergence phases. The global nature of this model represents very crudely and inaccurately the local dynamics of the buoyant rise, but allows to study the influence of global effects such as self-consistently generated differential rotation, meridional circulation and Coriolis forces. Although motivated by the solar context, this model cannot be thought of as a realistic model of the rise of magnetic structures and their emergence in the Sun where the local dynamics are completely different. The properties of initial phases of the buoyant rise in good agreement with previous studies. However, the effects of the interaction of the background dynamo field become increase as the flux-ropes evolve. During the buoyant rise across the CZ, the flux-rope's magnetic field strength and scales as $rho^a$, with $a\lesssim 1$. An increase of velocity, density and current precedes flux emergence at all longitudes. The geometry, latitude and relative orientation of the flux-ropes with respect to the background magnetic field influences the rise speeds, zonal flow amplitudes (which develop within the flux-ropes) and the corresponding surface signatures. This influences the morphology, duration and amplitude of the associated surface shearing and Poynting flux. The emerged flux influences the system's global polarity, leading in some cases to a polarity reversal while inhibiting background dynamo from doing so in some others. The emerged magnetic flux is slowly advected poleward, while being diffused and assimilated by the background dynamo field.

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.

Our global 3D simulations of convection and dynamo action in a Sun-like star reveal that persistent wreaths of strong magnetism can be built within the bulk of the convention zone. Here we examine the characteristics of buoyant magnetic structures that are self-consistently created by dynamo action and turbulent convective motions in a simulation with solar stratification but rotating at three times the current solar rate. These buoyant loops originate within sections of the magnetic wreaths in which turbulent flows amplify the fields to much larger values than is possible through laminar processes. These amplified portions can rise through the convective layer by a combination of magnetic buoyancy and advection by convective giant cells, forming buoyant loops. We measure statistical trends in the polarity, twist, and tilt of these loops. Loops are shown to preferentially arise in longitudinal patches somewhat reminiscent of active longitudes in the Sun, although broader in extent. We show that the strength of the axisymmetric toroidal field is not a good predictor of the production rate for buoyant loops or the amount of magnetic flux in the loops that are produced.

We present three-dimensional numerical simulations of a magnetic loop evolving in either a convectively stable or unstable rotating shell. The magnetic loop is introduced in the shell in such a way that it is buoyant only in a certain portion in longitude, thus creating an \Omega-loop. Due to the action of magnetic buoyancy, the loop rises and develops asymmetries between its leading and following legs, creating emerging bipolar regions whose characteristics are similar to the ones of observed spots at the solar surface. In particular, we self-consistently reproduce the creation of tongues around the spot polarities, which can be strongly affected by convection. We moreover emphasize the presence of ring-shaped magnetic structures around our simulated emerging regions, which we call "magnetic necklace" and which were seen in a number of observations without being reported as of today. We show that those necklaces are markers of vorticity generation at the periphery and below the rising magnetic loop. We also find that the asymmetry between the two legs of the loop is crucially dependent on the initial magnetic field strength. The tilt angle of the emerging regions is also studied in the stable and unstable cases and seems to be affected both by the convective motions and the presence of a differential rotation in the convective cases.

Solar-type stars exhibit a rich variety of magnetic activity. Seeking to explore the convective origins of this activity, we have carried out a series of global 3D magnetohydrodynamic (MHD) simulations with the anelastic spherical harmonic (ASH) code. Here we report on the dynamo mechanisms achieved as the effects of artificial diffusion are systematically decreased. The simulations are carried out at a nominal rotation rate of three times the solar value (3$\Omega_\odot$), but similar dynamics may also apply to the Sun. Our previous simulations demonstrated that convective dynamos can build persistent toroidal flux structures (magnetic wreaths) in the midst of a turbulent convection zone and that high rotation rates promote the cyclic reversal of these wreaths. Here we demonstrate that magnetic cycles can also be achieved by reducing the diffusion, thus increasing the Reynolds and magnetic Reynolds numbers. In these more turbulent models, diffusive processes no longer play a significant role in the key dynamical balances that establish and maintain the differential rotation and magnetic wreaths. Magnetic reversals are attributed to an imbalance in the poloidal magnetic induction by convective motions that is stabilized at higher diffusion levels. Additionally, the enhanced levels of turbulence lead to greater intermittency in the toroidal magnetic wreaths, promoting the generation of buoyant magnetic loops that rise from the deep interior to the upper regions of our simulated domain. The implications of such turbulence-induced magnetic buoyancy for solar and stellar flux emergence are also discussed.

The variable magnetic field of the solar photosphere exhibits periodic reversals as a result of dynamo activity occurring within the solar interior. We decompose the surface field as observed by both the Wilcox Solar Observatory and the Michelson Doppler Imager into its harmonic constituents, and present the time evolution of the mode coefficients for the past three sunspot cycles. The interplay between the various modes is then interpreted from the perspective of general dynamo theory, where the coupling between the primary and secondary families of modes is found to correlate with large-scale polarity reversals for many examples of cyclic dynamos. Mean-field dynamos based on the solar parameter regime are then used to explore how such couplings may result in the various long-term trends in the surface magnetic field observed to occur in the solar case.

Stellar magnetic activity of solar like stars is thought to be due to an internal dynamo. While the Sun has been the subject of intense research for refining dynamo models, observations of magnetic cyclic activity in solar type stars have become more and more available, opening a new path to understand the underlying physics behind stellar cycles. For instance, it is key to understand how stellar rotation rate influences magnetic cycle period $\Pcyc$. Recent numerical simulations of advection-dominated Babcock Leighton models have demonstrated that it is difficult to explain this observed trend given a) the strong influence of the cycle period to the meridional circulation amplitude and b) the fact that 3D models indicate that meridional flows become weaker as the rotation rate increases. In this paper, we introduce the turbulent pumping mechanism as another advective process capable also of transporting the magnetic fields. We found that this model is now able to reproduce the observations under the assumption that this effect increases as $\Omega^2$. The turbulent pumping becomes indeed another major player able to circumvent the meridional circulation. However, for high rotation rates ($\Omega \simeq 5 \Omega_\odot$), its effects dominate those of the meridional circulation, entering a new class of regime dominated by the advection of turbulent pumping and thus leading to a cyclic activity qualitatively different from that of the sun.

The stellar luminosity and depth of the convective envelope vary rapidly with mass for G- and K-type main sequence stars. In order to understand how these properties influence the convective turbulence, differential rotation, and meridional circulation, we have carried out 3D dynamical simulations of the interiors of rotating main sequence stars, using the anelastic spherical harmonic (ASH) code. The stars in our simulations have masses of 0.5, 0.7, 0.9, and 1.1 M_sun, corresponding to spectral types K7 through G0, and rotate at the same angular speed as the sun. We identify several trends of convection zone properties with stellar mass, exhibited by the simulations. The convective velocities, temperature contrast between up- and down-flows, and meridional circulation velocities all increase with stellar luminosity. As a consequence of the trend in convective velocity, the Rossby number (at a fixed rotation rate) increases and the convective turnover timescales decrease significantly with increasing stellar mass. The 3 lowest mass cases exhibit solar-like differential rotation, in a sense that they show a maximum rotation at the equator and minimum at higher latitudes, but the 1.1 M_sun case exhibits anti-solar rotation. At low mass, the meridional circulation is multi-cellular and aligned with the rotation axis; as the mass increases, the circulation pattern tends toward a unicellular structure covering each hemisphere in the convection zone.

The NASA Kepler mission is providing an unprecedented set of asteroseismic data. In particular, short-cadence lightcurves (~60s samplings), allow us to study solar-like stars covering a wide range of masses, spectral types and evolutionary stages. Oscillations have been observed in around 600 out of 2000 stars observed for one month during the survey phase of the Kepler mission. The measured light curves can present features related to the surface magnetic activity (starspots) and, thus we are able to obtain a good estimation of the surface (differential) rotation. In this work we establish the basis of such research and we show a potential method to find stars with fast surface rotations.

The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are believed to arise when sections of these magnetic structures become buoyantly unstable and rise from the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic (MHD) simulation in which turbulent convection, stratification, and rotation combine to yield a dynamo that self-consistently generates buoyant magnetic loops. We simulate stellar convection and dynamo action in a spherical shell with solar stratification, but rotating three times faster than the current solar rate. Strong wreaths of toroidal magnetic field are realized by dynamo action in the convection zone. By turning to a dynamic Smagorinsky model for subgrid-scale turbulence, we here attain considerably reduced diffusion in our simulation. This permits the regions of strongest magnetic field in these wreaths to rise toward the top of the convection zone via a combination of magnetic buoyancy instabilities and advection by convective giant cells. Such a global simulation yielding buoyant loops represents a significant step forward in combining numerical models of dynamo action and flux emergence.

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.

We study the connections between the sun's convection zone and the evolution of the solar wind and corona. We let the magnetic fields generated by a 2.5D axisymmetric kinematic dynamo code (STELEM) evolve in a 2.5D axisymmetric coronal isothermal MHD code (DIP). The computations cover an 11 year activity cycle. The solar wind's asymptotic velocity varies in latitude and in time in good agreement with the available observations. The magnetic polarity reversal happens at different paces at different coronal heights. Overall sun's mass loss rate, momentum flux and magnetic braking torque vary considerably throughout the cycle. This cyclic modulation is determined by the latitudinal distribution of the sources of open flux and solar wind and the geometry of the Alfv\'en surface. Wind sources and braking torque application zones also vary accordingly.

We have developed a variational data assimilation technique for the Sun using a toy {\alpha}{\Omega} dynamo model. The purpose of this work is to apply modern data assimilation techniques to solar data using a physically based model. This work represents the first step toward a complete variational model of solar magnetism. We derive the adjoint {\alpha}{\Omega} dynamo code and use a minimization procedure to invert the spatial dependence of key physical ingredients of the model. We find that the variational technique is very powerful and leads to encouraging results that will be applied to a more realistic model of the solar dynamo. We also note that the continuous adjoint mean-field induction equation is derived in the appendix.

Young solar-type stars rotate rapidly and many are magnetically active; some undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the convective envelope of a solar-type star rotating at 5 times the current solar rotation rate. Striking global-scale magnetic wreaths appear in the midst of the turbulent convection zone and show rich time-dependence. The dynamo exhibits cyclic activity and undergoes quasi-periodic polarity reversals where both the global-scale poloidal and toroidal fields change in sense on a roughly 1500 day time scale. These magnetic activity patterns emerge spontaneously from the turbulent flow and are more organized temporally and spatially than those realized in our previous simulations of the solar dynamo. We assess in detail the competing processes of magnetic field creation and destruction within our simulations that contribute to the global-scale reversals. We find that the mean toroidal fields are built primarily through an $\Omega$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $\alpha$-effect. During a reversal the magnetic wreaths propagate towards the polar regions, and this appears to arise from a poleward propagating dynamo wave. The primary response in the convective flows involves the axisymmetric differential rotation which shows variations associated with the poleward propagating magnetic wreaths. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent poleward propagating magnetic fields deep below the solar surface. [abridged]

We study the influence of stratification on stellar turbulent convection near the stellar surface and in depth by carrying out 3D high resolution hydrodynamic simulations with the ASH code. Four simulations with different radial density contrast corresponding to different aspect ratio for the same underlying 4 Myr $0.7 M_\odot$ pre-main sequence star model are thus performed. We highlight the existence of giant cells which are embedded in the complex surface convective patterns using a wavelet and time correlation analysis. Next, we study their properties such as their lifetime, aspect ratio and spatial extension in the different models both in latitude and depth according to the density contrast. We find that these giants cells have a lifetime larger than the stellar period with a typical longitudinal width of 490 Mm and a latitudinal extension increasing with the radial density contrast overpassing 50\degree of latitude for the thickest convective zone. Their rotation rate is much larger than the local differential rotation rate increasing also with the radial density contrast. However, their spatial coherence as a function of depth decrease with the density contrast due to the stronger shear present in these more stratified cases.

Young solar-type stars rotate rapidly and are very magnetically active. The magnetic fields at their surfaces likely originate in their convective envelopes where convection and rotation can drive strong dynamo action. Here we explore simulations of global-scale stellar convection in rapidly rotating suns using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. We explore one case rotates five times faster than the Sun in detail. This dynamo simulation, called case D5, has repeated quasi-cyclic reversals of global-scale polarity. We compare this case D5 to the broader family of simulations we have been able to explore and discuss how future simulations and observations can advance our understanding of stellar dynamos and magnetism.

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for a system rotating at three times the current solar rotation rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in this system. Striking wreaths of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a surprise, for it has been widely believed that such magnetic structures should be disrupted by magnetic buoyancy or turbulent pumping. Thus, many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. We examine how these persistent magnetic wreaths are maintained by dynamo processes and explore whether a classical mean-field $\alpha$-effect explains the regeneration of poloidal field.

Main-sequence massive stars possess convective cores that likely harbor strong dynamo action. To assess the role of core convection in building magnetic fields within these stars, we employ the 3-D anelastic spherical harmonic (ASH) code to model turbulent dynamics within a 10 solar mass main-sequence (MS) B-type star rotating at 4 times the solar rate. We find that strong (900 kG) magnetic fields arise within the turbulence of the core and penetrate into the stably stratified radiative zone. These fields exhibit complex, time-dependent behavior including reversals in magnetic polarity and shifts between which hemisphere dominates the total magnetic energy.

When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial or temporal organization. Wreath-building dynamos are present to some degree at all rotation rates, but are most evident in the more rapidly rotating simulations.

Observations of sun-like stars rotating faster than our current sun tend to exhibit increased magnetic activity as well as magnetic cycles spanning multiple years. Using global simulations in spherical shells to study the coupling of large-scale convection, rotation, and magnetism in a younger sun, we have probed effects of rotation on stellar dynamos and the nature of magnetic cycles. Major 3-D MHD simulations carried out at three times the current solar rotation rate reveal hydromagnetic dynamo action that yields wreaths of strong toroidal magnetic field at low latitudes, often with opposite polarity in the two hemispheres. Our recent simulations have explored behavior in systems with considerably lower diffusivities, achieved with sub-grid scale models including a dynamic Smagorinsky treatment of unresolved turbulence. The lower diffusion promotes the generation of magnetic wreaths that undergo prominent temporal variations in field strength, exhibiting global magnetic cycles that involve polarity reversals. In our least diffusive simulation, we find that magnetic buoyancy coupled with advection by convective giant cells can lead to the rise of coherent loops of magnetic field toward the top of the simulated domain.

We review recent insights into the dynamics of the solar convection zone obtained from global numerical simulations, focusing on two recent developments in particular. The first is quasi-cyclic magnetic activity in a long-duration dynamo simulation. Although mean fields comprise only a few percent of the total magnetic energy they exhibit remarkable order, with multiple polarity reversals and systematic variability on time scales of 6-15 years. The second development concerns the maintenance of the meridional circulation. Recent high-resolution simulations have captured the subtle nonlinear dynamical balances with more fidelity than previous, more laminar models, yielding more coherent circulation patterns. These patterns are dominated by a single cell in each hemisphere, with poleward and equatorward flow in the upper and lower convection zone respectively. We briefly address the implications of and future of these modeling efforts.

Aims: To understand stellar magnetism and to test the validity of the Babcock-Leighton flux transport mean field dynamo models with stellar activity observations Methods: 2-D mean field dynamo models at various rotation rates are computed with the STELEM code to study the sensitivity of the activity cycle period and butterfly diagram to parameter changes and are compared to observational data. The novelty is that these 2-D mean field dynamo models incorporate scaling laws deduced from 3-D hydrodynamical simulations for the influence of rotation rate on the amplitude and profile of the meridional circulation. These models make also use of observational scaling laws for the variation of differential rotation with rotation rate. Results: We find that Babcock-Leighton flux transport dynamo models are able to reproduce the change in topology of the magnetic field (i.e. toward being more toroidal with increasing rotation rate) but seem to have difficulty reproducing the cycle period vs activity period correlation observed in solar-like stars if a monolithic single cell meridional flow is assumed. It may however be possible to recover the $P_{cyc}$ vs $P_{rot}$ relation with more complex meridional flows, if the profile changes in a particular assumed manner with rotation rate. Conclusions: The Babcock-Leighton flux transport dynamo model based on single cell meridional circulation does not reproduce the $P_{cyc}$ vs $P_{rot}$ relation unless the amplitude of the meridional circulation is assumed to increase with rotation rate which seems to be in contradiction with recent results obtained with 3-D global simulations.

Context: The solar rotation profile is conical rather than cylindrical as one could expect from classical rotating fluid dynamics (e.g. Taylor-Proudman theorem). Thermal coupling to the tachocline, baroclinic effects and latitudinal transport of heat have been advocated to explain this peculiar state of rotation. Aims: To test the validity of thermal wind balance in the solar convection zone using helioseismic inversions for both the angular velocity and fluctuations in entropy and temperature. Methods: Entropy and temperature fluctuations obtained from 3-D hydrodynamical numerical simulations of the solar convection zone are compared with solar profiles obtained from helioseismic inversions. Results: The temperature and entropy fluctuations in 3-D numerical simulations have smaller amplitude in the bulk of the solar convection zone than those found from seismic inversions. Seismic inversion find variations of temperature from about 1 K at the surface up to 100 K at the base of the convection zone while in 3-D simulations they are of order 10 K throughout the convection zone up to 0.96 $R_{\odot}$. In 3-D simulations, baroclinic effects are found to be important to tilt the isocontours of $\Omega$ away from a cylindrical profile in most of the convection zone helped by Reynolds and viscous stresses at some locations. By contrast the baroclinic effect inverted by helioseismology are much larger than what is required to yield the observed angular velocity profile. Conclusion: The solar convection does not appear to be in strict thermal wind balance, Reynolds stresses must play a dominant role in setting not only the equatorial acceleration but also the observed conical angular velocity profile.

We present the first 3D MHD study in spherical geometry of the non-linear dynamical evolution of magnetic flux tubes in a turbulent rotating convection zone. We study numerically the rise of magnetic toroidal flux ropes from the base of a modelled convection zone up to the top of our computational domain where bipolar patches are formed. We compare the dynamical behaviour of flux tubes in a fully convective shell possessing self-consistently generated mean flows such as meridional circulation and differential rotation, with reference calculations done in a quiet isentropic zone. We find that two parameters influence the tubes during their rise through the convection zone: the initial field strength and amount of twist, thus confirming previous findings in Cartesian geometry. Further, when the tube is sufficiently strong with respect to the equipartition field, it rises almost radially independently of the initial latitude (either low or high). By contrast, weaker field cases indicate that downflows and upflows control the rising velocity of particular regions of the rope and could in principle favour the emergence of flux through Omega-loop structures. For these latter cases, we focus on the orientation of bipolar patches and find that sufficiently arched structures are able to create bipolar regions with a predominantly East-West orientation. Meridional flow seems to determine the trajectory of the magnetic rope when the field strength has been significantly reduced near the top of the domain. Finally differential rotation makes it more difficult for tubes introduced at low latitudes to reach the top of the domain.

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for two systems, rotating at three and five times the current solar rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in these systems. Striking wreathes of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a great surprise, for many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. Some dynamos achieved in these rapidly rotating states build persistent global-scale fields which maintain amplitude and polarity for thousands of days. In the case at five times the solar rate, the dynamo can undergo cycles of activity, with fields varying in strength and even changing polarity. As the magnetic fields wax and wane in strength, the primary response in the convective flows involves the axisymmetric differential rotation, which begins to vary on similar time scales. Bands of relatively fast and slow fluid propagate toward the poles on time scales of roughly 500 days. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent a response to poleward propagating magnetic field deep below the solar surface.

Detection of solar gravity modes remains a major challenge to our understanding of the innerparts of the Sun. Their frequencies would enable the derivation of constraints on the core physical properties while their amplitudes can put severe constraints on the properties of the inner convective region. Our purpose is to determine accurate theoretical amplitudes of solar g modes and estimate the SOHO observation duration for an unambiguous detection. We investigate the stochastic excitation of modes by turbulent convection as well as their damping. Input from a 3D global simulation of the solar convective zone is used for the kinetic turbulent energy spectrum. Damping is computed using a parametric description of the nonlocal time-dependent convection-pulsation interaction. We then provide a theoretical estimation of the intrinsic, as well as apparent, surface velocity. Asymptotic g-mode velocity amplitudes are found to be orders of magnitude higher than previous works. Using a 3D numerical simulation, from the ASH code, we attribute this to the temporal-correlation between the modes and the turbulent eddies which is found to follow a Lorentzian law rather than a Gaussian one as previously used. We also find that damping rates of asymptotic gravity modes are dominated by radiative losses, with a typical life-time of $3 \times 10^5$ years for the $\ell=1$ mode at $\nu=60 \mu$Hz. The maximum velocity in the considered frequency range (10-100 $\mu$Hz) is obtained for the $\ell=1$ mode at $\nu=60 \mu$Hz and for the $\ell=2$ at $\nu=100 \mu$Hz. Due to uncertainties in the modeling, amplitudes at maximum i.e. for $\ell=1$ at 60 $\mu$Hz can range from 3 to 6 mm s$^{-1}$.

We review recent advances in modeling global-scale convection and dynamo processes with the Anelastic Spherical Harmonic (ASH) code. In particular, we have recently achieved the first global-scale solar convection simulations that exhibit turbulent pumping of magnetic flux into a simulated tachocline and the subsequent organization and amplification of toroidal field structures by rotational shear. The presence of a tachocline not only promotes the generation of mean toroidal flux, but it also enhances and stabilizes the mean poloidal field throughout the convection zone, promoting dipolar structure with less frequent polarity reversals. The magnetic field generated by a convective dynamo with a tachocline and overshoot region is also more helical overall, with a sign reversal in the northern and southern hemispheres. Toroidal tachocline fields exhibit little indication of magnetic buoyancy instabilities but may be undergoing magneto-shear instabilities.

In the solar convection zone, rotation couples with intensely turbulent convection to drive a strong differential rotation and achieve complex magnetic dynamo action. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the global-scale patterns of convection in such stars and the flows of differential rotation and meridional circulation which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of localized convection emerges. Convection near the equator in these systems is dominated by one or two nests in longitude of locally enhanced convection, with quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation Delta Omega grows with more rapid rotation while the relative shear Delta Omega/Omega_0 decreases. In contrast, at more rapid rotation the meridional circulations decrease in energy and peak velocities and break into multiple cells of circulation in both radius and latitude.

In the solar convection zone, rotation couples with intensely turbulent convection to build global-scale flows of differential rotation and meridional circulation. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the patterns of convection in such stars and the global-scale flows which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of spatially localized convection emerges. Convection near the equator in these systems is dominated by one or two patches of locally enhanced convection, with nearly quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation, as measured by latitudinal angular velocity contrast, Delta_Omega, increases with more rapid rotation while the relative shear, Delta_Omega/Omega, decreases. In contrast, at more rapid rotation the meridional circulations decrease in both energy and peak velocities and break into multiple cells of circulation in both radius and latitude.

Stellar dynamos are driven by complex couplings between rotation and turbulent convection, which drive global-scale flows and build and rebuild stellar magnetic fields. When stars like our sun are young, they rotate much more rapidly than the current solar rate. Observations generally indicate that more rapid rotation is correlated with stronger magnetic activity and perhaps more effective dynamo action. Here we examine the effects of more rapid rotation on dynamo action in a star like our sun. We find that vigorous dynamo action is realized, with magnetic field generated throughout the bulk of the convection zone. These simulations do not possess a penetrative tachocline of shear where global-scale fields are thought to be organized in our sun, but despite this we find strikingly ordered fields, much like sea-snakes of toroidal field, which are organized on global scales. We believe this to be a novel finding.

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code. We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes.

The Sun is a magnetic star whose magnetism and cyclic activity is linked to the existence of an internal dynamo. We aim to understand the establishment of the solar magnetic 22-yr cycle, its associated butterfly diagram and field parity selection through numerical simulations of the solar global dynamo. Inspired by recent observations and 3D simulations that both exhibit multicellular flows in the solar convection zone, we seek to characterise the influence of various profiles of circulation on the behaviour of solar mean-field dynamo models. We are using 2-D mean field flux transport Babcock-Leighton numerical models in which we test several types of meridional flows: 1 large single cell, 2 cells in radius and 4 cells per hemisphere. We confirm that adding cells in latitude tends to speed up the dynamo cycle whereas adding cells in radius more than triples the period. We find that the cycle period in the four cells model is less sensitive to the flow speed than in the other simpler meridional circulation profiles studied. Moreover, our studies show that adding cells in radius or in latitude seems to favour the parity switching to a quadrupolar solution. According to our numerical models, the observed 22-yr cycle and dipolar parity is easily reproduced by models including multicellular meridional flows. On the contrary, the resulting butterfly diagram and phase relationship between the toroidal and poloidal fields are affected to a point where it is unlikely that such multicellular meridional flows persist for a long period of time inside the Sun, without having to reconsider the model itself.

We present the results of three-dimensional simulations of the deep convective envelope of a young (10 Myr) one-solar-mass star, obtained with the Anelastic Spherical Harmonic code. Since young stars are known to be faster rotators than their main sequence counterparts, we have systematically studied the impact of the stellar rotation speed, by considering stars spinning up to five times as fast as the Sun. The aim of these nonlinear models is to understand the complex interactions between convection and rotation. We discuss the influence of the turbulence level and of the rotation rate on the intensity and the topology of the mean flows. For all of the computed models, we find a solar-type superficial differential rotation, with an equatorial acceleration, and meridional circulation that exhibits a multicellular structure. Even if the differential rotation contrast decreases only marginally for high rotation rates, the meridional circulation intensity clearly weakens according to our simulations. We have also shown that, for Taylor numbers above a certain threshold (Ta>10^9), the convection can develop a vacillating behavior. Since simulations with high turbulence levels and rotation rates exhibit strongly cylindrical internal rotation profiles, we have considered the influence of baroclinic effects at the base of the convective envelope of these young Suns, to see whether such effect can modify the otherwise near cylindrical profiles to produce more conical, solar-like profiles.

We examine the MHD instabilities arising in the radiation zone of a differentially rotating star, in which a poloidal field of fossil origin is sheared into a toroidal field. We focus on the non-axisymmetric instability that affects the toroidal magnetic field in a rotating star, which was first studied by Pitts and Tayler in the non-dissipative limit. According to Spruit, it could also drive a dynamo. The Pitts & Tayler instability is manifestly present in our simulations, with its conspicuous m=1 dependence in azimuth. But its analytic treatment used so far is too simplified to be applied to the real stellar situation. Although the instability generated field reaches an energy comparable to that of the mean poloidal field, that field seems unaffected by the instability: it undergoes Ohmic decline, and is neither eroded nor regenerated by the instability. The toroidal field is produced by shearing the poloidal field and it draws its energy from the differential rotation. The small scale motions behave as Alfven waves; they cause negligible eddy-diffusivity and contribute little to the net transport of angular momentum. In our simulations we observe no sign of dynamo action, of either mean field or fluctuation type, up to a magnetic Reynolds number of 10^5. However the Pitts & Tayler instability is sustained as long as the differential rotation acting on the poloidal field is able to generate a toroidal field of sufficient strength.

The global scales of solar convection are studied through three-dimensional simulations of compressible convection carried out in spherical shells of rotating fluid which extend from the base of the convection zone to within 15 Mm of the photosphere. Such modelling at the highest spatial resolution to date allows study of distinctly turbulent convection, revealing that coherent downflow structures associated with giant cells continue to play a significant role in maintaining the strong differential rotation that is achieved. These giant cells at lower latitudes exhibit prograde propagation relative to the mean zonal flow, or differential rotation, that they establish, and retrograde propagation of more isotropic structures with vortical character at mid and high latitudes. The interstices of the downflow networks often possess strong and compact cyclonic flows. The evolving giant-cell downflow systems can be partly masked by the intense smaller scales of convection driven closer to the surface, yet they are likely to be detectable with the helioseismic probing that is now becoming available. Indeed, the meandering streams and varying cellular subsurface flows revealed by helioseismology must be sampling contributions from the giant cells, yet it is difficult to separate out these signals from those attributed to the faster horizontal flows of supergranulation. To aid in such detection, we use our simulations to describe how the properties of giant cells may be expected to vary with depth, how their patterns evolve in time, and analyze the statistical features of correlations within these complex flow fields.

The operation of the solar global dynamo appears to involve many dynamical elements. Self-consistent MHD simulations which realistically incorporate all of these processes are not yet computationally feasible, though some elements can now be studied with reasonable fidelity. Here we consider the manner in which turbulent compressible convection within the bulk of the solar convection zone can generate large-scale magnetic fields through dynamo action. We accomplish this through a series of three-dimensional numerical simulations of MHD convection within rotating spherical shells using our ASH code on massively parallel supercomputers. Since differential rotation is a key ingredient in all dynamo models, we also examine here the nature of the rotation profiles that can be sustained within the deep convection zone as strong magnetic fields are built and maintained. We find that the convection is able to maintain a solar-like angular velocity profile despite the influence of Maxwell stresses which tend to oppose Reynolds stresses and thus reduce the latitudinal angular velocity contrast throughout the convection zone. The dynamo-generated magnetic fields exhibit a complex structure and evolution, with radial fields concentrated in downflow lanes and toroidal fields organized into twisted ribbons which are extended in longitude and which achieve field strengths of up to 5000 G. The flows and fields exhibit substantial kinetic and magnetic helicity although systematic hemispherical patterns are only apparent in the former. Fluctuating fields dominate the magnetic energy and account for most of the back-reaction on the flow via Lorentz forces. Mean fields are relatively weak and do not exhibit systematic latitudinal propagation or periodic polarity reversals as in the sun. This may be attributed to the absence of a tachocline.

Core convection and dynamo activity deep within rotating A-type stars of 2 solar masses are studied with 3--D nonlinear simulations. Our modeling considers the inner 30% by radius of such stars, thus capturing within a spherical domain the convective core and a modest portion of the surrounding radiative envelope. The MHD equations are solved using the ASH code to examine turbulent flows and magnetic fields, both of which exhibit intricate time dependence. By introducing small seed magnetic fields into our progenitor hydrodynamic models rotating at one and four times the solar rate, we assess here how the vigorous convection can amplify those fields and sustain them against ohmic decay. Dynamo action is indeed realized, ultimately yielding magnetic fields that are in energy equipartion with the flow. Such magnetism reduces the differential rotation obtained in the progenitors, partly by Maxwell stresses that transport angular momentum poleward and oppose the Reynolds stresses in the latitudinal balance. In contrast, in the radial direction we find that the Maxwell and Reynolds stresses may act together to transport angular momentum. The central columns of slow rotation established in the progenitors are weakened, with the differential rotation waxing and waning in strength as the simulations evolve. We assess the morphology of the flows and magnetic fields, their complex temporal variations, and the manner in which dynamo action is sustained. Differential rotation and helical convection are both found to play roles in giving rise to the magnetic fields. The magnetism is dominated by strong fluctuating fields throughout the core, with the axisymmetric (mean) fields there relatively weak.

We study the physics of the solar tachocline and related MHD instabilities. We have performed 3-D MHD simulations of the solar radiative interior to check whether a fossil magnetic field is able to prevent the spread of the tachocline. Starting with a purely poloidal magnetic field and a latitudinal shear meant to be imposed by the convection zone at the top of the radiation zone, we have investigated the interactions between magnetic fields, rotation and shear, using the spectral code ASH on massive parallel supercomputers. In all cases we have explored, the fossil field diffuses outward and ends up connecting with the convection zone, whose differential rotation is then imprinted at latitudes above 40 deg throughout the radiative interior, according to Ferraro's law of isorotation. Rotation remains uniform in the lower latitude region which is contained within closed field lines. We find that the meridional flow cannot stop the inward progression of the differential rotation. Further, we observe the development of non-axisymmetric magnetohydrodynamic instabilities, first due to the initial poloidal configuration of the fossil field, and later to the toroidal field produced by shearing the poloidal field through the differential rotation. We do not find dynamo action as such in the radiative interior, since the mean poloidal field is not regenerated. But the instability persists during the whole evolution, while slowly decaying with the mean poloidal field. According to our numerical simulations, a fossil magnetic field cannot prevent the radiative spread of the tachocline, and thus it is unable to enforce uniform rotation in the radiation zone. Neither can the observed thinness of that layer be invoked as a proof for such an internal fossil magnetic field.

We have performed the first three-dimensional non-linear simulation of the turbulent convective envelope of a rotating 0.8 Msun RGB star using the ASH code. Adopting a global typical rotation rate of a tenth of the solar rate, we have analyzed the dynamical properties of the convection and the transport of angular momentum within the inner 50% in radius of the convective envelope. The convective patterns consist of a small number of large cell, associated with fast flows (about 3000 m/s) and large temperature fluctuations (about 300 K) in order to carry outward the large luminosity (L* = 400 Lsun) of the star. The interactions between convection and rotation give rise to a large radial differential rotation and a meridional circulation possessing one cell per hemisphere, the flow being poleward in both hemisphere. By analysing the redistribution of angular momentum, we find that the meridional circulation transports the angular momentum outward in the radial direction, and poleward in the latitudinal direction, and that the transport by Reynolds stresses acts in the opposite direction. From this 3-D simulation, we have derived an average radial rotation profile, that we will ultimately introduce back into 1-D stellar evolution code.

We present the first results from three-dimensional spherical shell simulations of magnetic dynamo action realized by turbulent convection penetrating downward into a tachocline of rotational shear. This permits us to assess several dynamical elements believed to be crucial to the operation of the solar global dynamo, variously involving differential rotation resulting from convection, magnetic pumping, and amplification of fields by stretching within the tachocline. The simulations reveal that strong axisymmetric toroidal magnetic fields (about 3000 G in strength) are realized within the lower stable layer, unlike in the convection zone where fluctuating fields are predominant. The toroidal fields in the stable layer possess a striking persistent antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. The associated mean poloidal magnetic fields there have a clear dipolar geometry, but we have not yet observed any distinctive reversals or latitudinal propagation. The presence of these deep magnetic fields appears to stabilize the sense of mean fields produced by vigorous dynamo action in the bulk of the convection zone.

The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles.

We present the results of 3--D nonlinear simulations of magnetic dynamo action by core convection within A-type stars of 2 solar masses, at a range of rotation rates. We consider the inner 30% by radius of such stars, with the spherical domain thereby encompassing the convective core and a portion of the surrounding radiative envelope. The compressible Navier-Stokes equations, subject to the anelastic approximation, are solved to examine highly nonlinear flows that span multiple scale heights, exhibit intricate time dependence, and admit magnetic dynamo action. Small initial seed magnetic fields are found to be amplified greatly by the convective and zonal flows. The central columns of strikingly slow rotation found in some of our progenitor hydrodynamic simulations continue to be realized in some simulations to a lesser degree, with such differential rotation arising from the redistribution of angular momentum by the nonlinear convection and magnetic fields. We assess the properties of the magnetic fields thus generated, the extent of convective penetration, the magnitude of the differential rotation, and the excitation of gravity waves within the radiative envelope.

We present the results of 3--D simulations of core convection within A-type stars of 2 solar masses, at a range of rotation rates. We consider the inner 30% by radius of such stars, thereby encompassing the convective core and some of the surrounding radiative envelope. We utilize our anelastic spherical harmonic (ASH) code, which solves the compressible Navier-Stokes equations in the anelastic approximation, to examine highly nonlinear flows that can span multiple scale heights. The cores of these stars are found to rotate differentially, with central cylindrical regions of strikingly slow rotation achieved in our simulations of stars whose convective Rossby number (R_{oc}) is less than unity. Such differential rotation results from the redistribution of angular momentum by the nonlinear convection that strongly senses the overall rotation of the star. Penetrative convective motions extend into the overlying radiative zone, yielding a prolate shape (aligned with the rotation axis) to the central region in which nearly adiabatic stratification is achieved. This is further surrounded by a region of overshooting motions, the extent of which is greater at the equator than at the poles, yielding an overall spherical shape to the domain experiencing at least some convective mixing. We assess the overshooting achieved as the stability of the radiative exterior is varied, and the weak circulations that result in that exterior. The convective plumes serve to excite gravity waves in the radiative envelope, ranging from localized ripples of many scales to some remarkable global resonances.

We discuss recent advances made in modelling the complex magnetohydrodynamics of the Sun using our anelastic spherical harmonics (ASH) code. We have conducted extensive 3--D simulations of compressible convection in rotating spherical shells with and without magnetic fields, to study the coupling between global-scale convection and rotation in seeking to understand how the solar differential rotation is established and maintained. Such simulations capable of studying fairly turbulent convection have been enabled by massively parallel supercomputers. The resulting convection within domains that capture a good fraction of the bulk of the convection zone is highly time dependent and intricate, and is dominated by intermittent upflows and networks of strong downflows. A high degree of coherent structures involving downflowing plumes can be embedded in otherwise chaotic flow fields. These vortical structures play a significant role in yielding Reynolds stresses that serve to redistribute angular momentum, leading to differential rotation profiles with pole to equator contrasts of about 30% in angular velocity $\Omega$ and some constancy along radial line at mid latitudes, thereby making good contact with deductions from helioseismology. When a magnetic field is introduced, a dynamo regime can be found that does not destroy the strong differential rotation achieved in pure hydrodynamics cases. The magnetic fields are found to concentrate around the downflowing networks and to have significant north-south asymmetry and helicity.

The advent of massively parallel supercomputing has begun to permit explicit 3--D simulations of turbulent convection occurring within the cores of early-type main sequence stars. Such studies should complement the stellar structure and evolution efforts that have so far largely employed 1--D nonlocal mixing length descriptions for the transport, mixing and overshooting achieved by core convection. We have turned to A-type stars as representative of many of the dynamical challenges raised by core convection within rotating stars. The differential rotation and meridional circulations achieved deep within the star by the convection, the likelihood of sustained magnetic dynamo action there, and the bringing of fresh fuel into the core by overshooting motions, thereby influencing main sequence lifetimes, all constitute interesting dynamical questions that require detailed modelling of global-scale convection. Using our anelastic spherical harmonic (ASH) code tested on the solar differential rotation problem, we have conducted a series of 3--D spherical domain simulations that deal with a simplified description of the central regions of rotating A-type stars, i.e a convectively unstable core is surrounded by a stable radiative envelope. A sequence of 3--D simulations are used to assess the properties of the convection (its global patterns, differential rotation, meridional circulations, extent and latitudinal variation of the overshooting) as transitions are made between laminar and turbulent states by changing the effective diffusivities, rotation rates, and subadiabaticity of the radiative exterior. We report on the properties deduced from these models for both the extent of penetration and the profile of rotation sustained by the convection.

Many stars exhibit strong magnetic fields, some of which are thought to be of primordial origin and others a sign of magnetic dynamo processes. We briefly review the results of observational studies of solar-type stars seeking to evaluate the linkage between rotation rate and possible magnetic cycles of activity. Clearly turbulent convection and rotation within spherical shell geometries provide ingredients essential for dynamo action. However, intensive efforts over several decades in solar research have demonstrated that it is no easy matter to achieve cyclic magnetic activity that is in accord with observations. Helioseismology has revealed that an essential element for the global solar dynamo is the presence of a tachocline of shear at the base of the solar convection zone, leading to the likely operation of an interface dynamo. We review the crucial elements for achieving a cyclic magnetic activity. We then discuss some of our current 3--D MHD simulations of solar turbulent convection in spherical shells that yield differential rotation profiles which make good contact with some of the helioseismic findings. We show that such turbulent motions can amplify and sustain magnetic field in the bulk of the convective zone whose strength are sufficient to feed back both upon the convection and its global circulations.

The intense turbulence present in the solar convection zone is a major challenge to both theory and simulation as one tries to understand the origins of the striking differential rotation profile with radius and latitude that has been revealed by helioseismology. The differential rotation must be an essential element in the operation of the solar magnetic dynamo and its cycles of activity, yet there are many aspects of the interplay between convection, rotation and magnetic fields that are still unclear. We have here carried out a series of 3--D numerical simulations of turbulent convection within deep spherical shells using our anelastic spherical harmonic (ASH) code on massively parallel supercomputers. These studies of the global dynamics of the solar convection zone concentrate on how the differential rotation and meridional circulation are established. We have analyzed the transport of angular momentum in establishing such differential rotation, and clarified the roles played by Reynolds stresses and the meridional circulation in this process. We have found that the Reynolds stresses are crucial in transporting angular momentum toward the equator. The effects of baroclinicity (thermal wind) have been found to have a modest role in the resulting mean zonal flows. The simulations have produced differential rotation profiles within the bulk of the convection zone that make reasonable contact with ones inferred from helioseismic inversions, namely possessing a fast equator, an angular velocity difference of about 30% from equator to pole, and some constancy along radial lines at mid-latitudes.

We have computed accurate 1-D solar models including both a macroscopic mixing process in the solar tachocline as well as up-to-date microscopic physical ingredients. Using sound speed and density profiles inferred through primary inversion of the solar oscillation frequencies coupled with the equation of thermal equilibrium, we have extracted the temperature and hydrogen abundance profiles. These inferred quantities place strong constraints on our theoretical models in terms of the extent and strength of our macroscopic mixing, on the photospheric heavy elements abundance, on the nuclear reaction rates such as $S_{11}$ and $S_{34}$ and on the efficiency of the microscopic diffusion. We find a good overall agreement between the seismic Sun and our models if we introduce a macroscopic mixing in the tachocline and allow for variation within their uncertainties of the main physical ingredients. From our study we deduce that the solar hydrogen abundance at the solar age is $X_{\rm inv}=0.732\pm 0.001$ and that based on the $^9$Be photospheric depletion, the maximum extent of mixing in the tachocline is 5% of the solar radius. The nuclear reaction rate for the fundamental $pp$ reaction is found to be $S_{11}(0)=4.06\pm 0.07$ $10^{-25}$ MeV barns, i.e., 1.5% higher than the present theoretical determination. The predicted solar neutrino fluxes are discussed in the light of the new SNO/SuperKamiokande results.

Recently helioseismic observations have revealed the presence of a shear layer at the base of the convective zone related to the transition from differential rotation in the convection zone to almost uniform rotation in the radiative interior, the tachocline. At present, this layer extends only over a few percent of the solar radius and no definitive explanations have been given for this thiness. Following Spiegel and Zahn (1992, Astron. Astrophys.), who invoke anisotropic turbulence to stop the spread of the tachocline deeper in the radiative zone as the Sun evolves, we give some justifications for their hypothesis by taking into account recent results on rotating shear instability (Richard and Zahn 1999, Astron. Astrophys.). We study the impact of the macroscopic motions present in this layer on the Sun's structure and evolution by introducing a macroscopic diffusivity $D_T$ in updated solar models. We find that a time dependent treatment of the tachocline significantly improves the agreement between computed and observed surface chemical species, such as the $^7$Li and modify the internal structure of the Sun (Brun, Turck-Chi\`eze and Zahn, 1999, in Astrophys. J.).

Recent studies have established that peaks in solar oscillation power spectra are not Lorentzian in shape, but have a distinct asymmetry. Fitting a symmetric Lorentzian profile to the peaks therefore produces a shift in frequency of the modes. Accurate determination of low-frequency modes is essential to infer the structure of the solar core by inversion of the mode frequencies. In this paper we investigate how the changes in frequencies of low-degree modes obtained by fitting symmetric and asymmetric peak profiles change the inferred properties of the solar core. We use data obtained by the Global Oscillations at Low Frequencies (GOLF) project on board the SoHO spacecraft. Two different solar models and inversion procedures are used to invert the data to determine the sound speed in the solar core. We find that for a given set of modes no significant difference in the inferred sound-speed results from taking asymmetry into account when fitting the low-degree modes.

In previous work, we have shown that recent updated standard solar models cannot reproduce the radial profile of the sound speed at the base of the convective zone (CZ) and fail to predict the Li7 depletion. In parallel, helioseismology has shown that the transition from differential rotation in the CZ to almost uniform rotation in the radiative solar interior occurs in a shallow layer called the tachocline. This layer is presumably the seat of large scale circulation and of turbulent motions. Here, we introduce a macroscopic transport term in the structure equations, which is based on a hydrodynamical description of the tachocline proposed by Spiegel and Zahn, and we calculate the mixing induced within this layer. We discuss the influence of different parameters that represent the tachocline thickness, the Brunt-Vaissala frequency at the base of the CZ, and the time dependence of this mixing process along the Sun's evolution. We show that the introduction of such a process inhibits the microscopic diffusion by about 25%. Starting from models including a pre-main sequence evolution, we obtain: a) a good agreement with the observed photospheric chemical abundance of light elements such as He3, He4, Li7 and Be9, b) a smooth composition gradient at the base of the CZ, and c) a significant improvement of the sound speed square difference between the seismic sun and the models in this transition region, when we allow the phostospheric heavy element abundance to adjust, within the observational incertitude, due to the action of this mixing process. The impact on neutrino predictions is also discussed.

With the recent results of helioseismology aboard SOHO, solar models are more and more constrained (Brun, Turck-Chieze and Morel 1998 (astro-ph/9806272)). New physical processes, mainly connected to macroscopic motions, must be introduced to understand these new observations. In this poster, we present solar models including a turbulent pressure in the outer layers and mixing due to the tachocline (Spiegel and Zahn 1992). Our results lead us to conclude that: - Turbulent pressure improves the absolute value of the acoustic mode frequencies (~ 10 microHz at 4mHz), - Mixing in a tachocline thickness of 0.05 +/- 0.03 Ro (Corbard et al. 1997) looks promising.

In this paper, we examine a new updated solar model taking advantage of recent reexamination of the nuclear reactions rates and of the microscopic diffusion of helium and heavy elements. Our best model fits the helioseismic data reasonably well: the base of the convective zone is at r=0.715 Ro, the photospheric helium in mass fraction is 0.243 and the sound speed square difference between the Sun and the model is < 1%. This model leads to a reestimate of the neutrinos fluxes : 7.18 SNU for the chlorine experiment, 127.2 SNU for the gallium detector and 4.82 10^6/cm2/s for the B8 neutrino flux. Acoustic and gravity mode predictions are also estimated. We present 3 models which respect the helioseismic constraints and fit the seismic observations better by taking advantage of the known physical uncertainties: nuclear reaction rates, CNO abundances and microscopic diffusion. We conclude that the standard model in side its inherent uncertainties, is robust in the light of the present acoustic mode detection and that the neutrino fluxes prediction is not well constrained.