The existence of a stably stratified layer underneath the core-mantle boundary (CMB) has been recently revived by corroborating evidences coming from seismic studies, mineral physics and thermal evolution models. Such a layer could find its physical origination either in compositional stratification due to the accumulation of light elements at the top or the core or in thermal stratification due to the heat flux becoming locally sub-adiabatic. The exact properties of this stably-stratified layer, namely its size $\mathcal{H}_S$ and the degree of its stratification characterised by the Brunt-V\"ais\"al\"a frequency $N$, are however uncertain and highly debated. A stable layer underneath the CMB can have crucial dynamical impacts on the geodynamo. Because of the inhibition of the convective motions, a stable layer is expected to primarily act as a low-pass filter on the magnetic field, smoothing out the rapidly-varying and small-scale features by skin effect. To investigate this effect more systematically, we compute 70 global geodynamo models varying the size of the stably-stratified layer from 0 to 300~km and its amplitude from $N/\Omega = 0$ to $N/\Omega \simeq 50$, $\Omega$ being the rotation rate. We show that the penetration of the convective flow in the stably-stratified layer is controlled by the typical size of the convective eddies and by the local variations of the ratio $N/\Omega$. Using quantitative measures of the degree of morphological semblance between the magnetic field obtained in numerical models and the geomagnetic field at the CMB, we establish an upper bound for the stable layer thickness $\mathcal{H}_s < (N/\Omega)^{-1} d_c$, $d_c$ being the horizontal size of the convective flow at the base of the stable layer. This defines a strong geomagnetic constraint on the properties of a stably-stratified layer beneath the CMB.

The release of spin-down energy by a magnetar is a promising scenario to power several classes of extreme explosive transients. However, it lacks a firm basis because magnetar formation still represents a theoretical challenge. Using the first three-dimensional simulations of a convective dynamo based on a protoneutron star interior model, we demonstrate that the required dipolar magnetic field can be consistently generated for sufficiently fast rotation rates. The dynamo instability saturates in the magnetostrophic regime with the magnetic energy exceeding the kinetic energy by a factor of up to 10. Our results are compatible with the observational constraints on galactic magnetar field strength and provide strong theoretical support for millisecond protomagnetar models of gamma-ray burst and superluminous supernova central engines.

The new data delivered by NASA's Juno spacecraft significantly increase our understanding of Jupiter's internal dynamics. The gravity data constrain the depth of the zonal flows observed at cloud level and suggest that they slow down considerably at a depth of about $0.96\,r_J$, where $r_J$ is the mean radius at the one bar level. Juno's magnetometer reveals the planet's internal magnetic field. We combine the new zonal flow and magnetic field models with an updated electrical conductivity profile to assess the zonal wind induced dynamo action, concentrating on the outer part of Jupiter's molecular hydrogen region where the conductivity increases very rapidly with depth. Dynamo action remains quasi-stationary and can thus reasonably be estimated where the magnetic Reynolds number remains smaller than one, which is roughly the region above $0.96\,r_J$. We calculate that the locally induced radial magnetic field reaches rms values of about $10^{-6}\,$T in this region and may just be detectable by the Juno mission. Very localized dynamo action and a distinct pattern that reflects the zonal wind system increases the chance to disentangle this locally induced field from the background field. The estimates of the locally induced currents also allow calculating the zonal flow related Ohmic heating and associated entropy production. The respective quantities remain below new revised predictions for the total dissipative heating and total entropy production in Jupiter for any of the explored model combinations. Thus neither Ohmic heating nor entropy production offer additional constraints on the depth of the zonal winds.

We present a new pseudo-spectral open-source code nicknamed pizza. It is dedicated to the study of rapidly-rotating Boussinesq convection under the 2-D spherical quasi-geostrophic approximation, a physical hypothesis that is appropriate to model the turbulent convection that develops in planetary interiors. The code uses a Fourier decomposition in the azimuthal direction and supports both a Chebyshev collocation method and a sparse Chebyshev integration formulation in the cylindrically-radial direction. It supports several temporal discretisation schemes encompassing multi-step time steppers as well as diagonally-implicit Runge-Kutta schemes. The code has been tested and validated by comparing weakly-nonlinear convection with the eigenmodes from a linear solver. The comparison of the two radial discretisation schemes has revealed the superiority of the Chebyshev integration method over the classical collocation approach both in terms of memory requirements and operation counts. The good parallelisation efficiency enables the computation of large problem sizes with $\mathcal{O}(10^4\times 10^4)$ grid points using several thousands of ranks. This allows the computation of numerical models in the turbulent regime of quasi-geostrophic convection characterised by large Reynolds $Re$ and yet small Rossby numbers $Ro$. A preliminary result obtained for a strongly supercritical numerical model with a small Ekman number of $10^{-9}$ and a Prandtl number of unity yields $Re\simeq 10^5$ and $Ro \simeq 10^{-4}$. pizza is hence an efficient tool to study spherical quasi-geostrophic convection in a parameter regime inaccessible to current global 3-D spherical shell models.

The Juno mission is delivering spectacular data of Jupiter's magnetic field, while the gravity measurements finally allow constraining the depth of the winds observed at cloud level. However, to which degree the zonal winds contribute to the planet's dynamo action remains an open question. Here we explore numerical dynamo simulations that include an Jupiter-like electrical conductivity profile and successfully model the planet's large scale field. We concentrate on analyzing the dynamo action in the Steeply Decaying Conductivity Region (SDCR) where the high conductivity in the metallic Hydrogen region drops to the much lower values caused by ionization effects in the very outer envelope of the planet. Our simulations show that the dynamo action in the SDCR is strongly ruled by diffusive effects and therefore quasi stationary. The locally induced magnetic field is dominated by the horizontal toroidal field, while the locally induced currents flow mainly in the latitudinal direction. The simple dynamics can be exploited to yield estimates of surprisingly high quality for both the induced field and the electric currents in the SDCR. These could be potentially be exploited to predict the dynamo action of the zonal winds in Jupiter's SDCR but also in other planets.

Recent interferometric data have been used to constrain the brightness distribution at the surface of nearby stars, in particular the so-called gravity darkening that makes fast rotating stars brighter at their poles than at their equator. However, good models of gravity darkening are missing for stars that posses a convective envelope. In order to better understand how rotation affects the heat transfer in stellar convective envelopes, we focus on the heat flux distribution in latitude at the outer surface of numerical models. We carry out a systematic parameter study of three-dimensional, direct numerical simulations of anelastic convection in rotating spherical shells. Restricting our investigations to hydrodynamical models with a thermal Prandtl number fixed to unity, we consider both thick and thin (solar-like) shells, and vary the stratification over three orders of magnitude. We measure the heat transfer efficiency in terms of the Nusselt number, defined as the output luminosity normalised by the conductive state luminosity. We report diverse Nusselt number profiles in latitude, ranging from brighter (usually at the onset of convection) to darker equator and uniform profiles. We find that the variations of the surface brightness are mainly controlled by the surface value of the local Rossby number: when the Coriolis force dominates the dynamics, the heat flux is weakened in the equatorial region by the zonal wind and enhanced at the poles by convective motions inside the tangent cylinder. In the presence of a strong background density stratification however, as expected in real stars, the increase of the local Rossby number in the outer layers leads to uniformisation of the surface heat flux distribution.

The Juno mission will measure Jupiter's magnetic field with unprecedented precision and provide a wealth of additional data that will allow to constrain the planet's interior structure and dynamics. Here we analyse 66 numerical simulations in order to explore the sensitivity of the dynamo-generated magnetic field to the planets interior properties. The degree l=4 field model VIP4 and up-to-date interior models based on ab initio simulations serve as benchmarks. Our results suggest that VIP4-like magnetic fields can be found for a number of different models. We find that whether we assume an ideal gas or use the more realistic interior model based on ab initio simulations makes no difference. However, two other factors are important. Low Rayleigh number leads to strong axial dipole contribution while the axial dipole dominance is lost when the convective driving is too strong. The required intermediate range that yields Jupiter-like magnetic fields depends on the other system properties. The second factor is the convective magnetic Reynolds number profile Rmc(r), a product of the non-axisymmetric flow velocity and electrical conductivity. We find that the depth where Rmc exceeds about 50 is a good proxy for the top of the dynamo region. When the dynamo region sits too deep, the axial dipole is once more too dominant due to geometric reasons. We extrapolated our results to Jupiter and the result suggests that the Jovian dynamo extends to 95% of the planetary radius. The zonal flows in our simulations are dominated by an equatorial jet largely confined to the molecular layer. Where the jet reaches down to higher electrical conductivities, it is gives rise to a secondary alpha-Omega dynamo that modifies the dipole-dominated field produced deeper in the planet. This secondary dynamo can lead to strong magnetic field patches at low latitudes that seem compatible with the VIP4 model.

Earth sustains its magnetic field by a dynamo process driven by convection in the liquid outer core. Geodynamo simulations have been successful in reproducing many observed properties of the geomagnetic field. However, while theoretical considerations suggest that flow in the core is governed by a balance between Lorentz force, rotational force and buoyancy (called MAC balance for Magnetic, Archimedean, Coriolis) with only minute roles for viscous and inertial forces, dynamo simulations must employ viscosity values that are many orders of magnitude larger than in the core due to computational constraints. In typical geodynamo models viscous and inertial forces are not much smaller than the Coriolis force and the Lorentz force plays a sub-dominant role. This has led to conclusions that these simulations are viscously controlled and do not represent the physics of the geodynamo. Here we show by a direct analysis of the relevant forces that a MAC balance can be achieved when the viscosity is reduced to values close to the current practical limit. Lorentz force, buoyancy and the uncompensated (by pressure) part of the Coriolis force are of very similar strength, whereas viscous and inertia are smaller by a factor of at least 20 in the bulk of the fluid volume. Compared to non-magnetic convection at otherwise identical parameters, the dynamo flow is of larger scale, less invariant parallel to the rotation axis (less geostrophic) and convection transports twice as much heat, all of which is expected when the Lorentz force strongly influences the convection properties.

Rayleigh-B\'enard convection in rotating spherical shells can be considered as a simplified analogue of many astrophysical and geophysical fluid flows. Here, we use three-dimensional direct numerical simulations to study this physical process. We construct a dataset of more than 200 numerical models that cover a broad parameter range with Ekman numbers spanning $3\times 10^{-7} \leq E \leq 10^{-1}$, Rayleigh numbers within the range $10^3 < Ra < 2\times 10^{10}$ and a Prandtl number unity. We investigate the scaling behaviours of both local (length scales, boundary layers) and global (Nusselt and Reynolds numbers) properties across various physical regimes from onset of rotating convection to weakly-rotating convection. Close to critical, the convective flow is dominated by a triple force balance between viscosity, Coriolis force and buoyancy. For larger supercriticalities, a subset of our numerical data approaches the asymptotic diffusivity-free scaling of rotating convection $Nu\sim Ra^{3/2}E^{2}$ in a narrow fraction of the parameter space delimited by $6\,Ra_c \leq Ra \leq 0.4\,E^{-8/5}$. Using a decomposition of the viscous dissipation rate into bulk and boundary layer contributions, we establish a theoretical scaling of the flow velocity that accurately describes the numerical data. In rapidly-rotating turbulent convection, the fluid bulk is controlled by a triple force balance between Coriolis, inertia and buoyancy, while the remaining fraction of the dissipation can be attributed to the viscous friction in the Ekman layers. Beyond $Ra \simeq E^{-8/5}$, the rotational constraint on the convective flow is gradually lost and the flow properties vary to match the regime changes between rotation-dominated and non-rotating convection. The quantity $Ra E^{12/7}$ provides an accurate transition parameter to separate rotating and non-rotating convection.

Zonal flows in rapidly-rotating celestial objects such as the Sun, gas or ice giants form in a variety of surface patterns and amplitudes. Whereas the differential rotation on the Sun, Jupiter and Saturn features a super-rotating equatorial region, the ice giants, Neptune and Uranus harbour an equatorial jet slower than the planetary rotation. Global numerical models covering the optically thick, deep-reaching and rapidly rotating convective envelopes of gas giants reproduce successfully the prograde jet at the equator. In such models, convective columns shaped by the dominant Coriolis force typically exhibit a consistent prograde tilt. Hence angular momentum is pumped away from the rotation axis via Reynolds stresses. Those models are found to be strongly geostrophic, hence a modulation of the zonal flow structure along the axis of rotation, e.g. introduced by persistent latitudinal temperature gradients, seems of minor importance. Within our study we stimulate these thermal gradients and the resulting ageostrophic flows by applying an axisymmetric and equatorially symmetric outer boundary heat flux anomaly ($Y_{20}$) with variable amplitude and sign. Such a forcing pattern mimics the thermal effect of intense solar or stellar irradiation. Our results suggest that the ageostrophic flows are linearly amplified with the forcing amplitude $q^\star$ leading to a more pronounced dimple of the equatorial jet (alike Jupiter). The geostrophic flow contributions, however, are suppressed for weak $q^\star$, but inverted and re-amplified once $q^\star$ exceeds a critical value. The inverse geostrophic differential rotation is consistently maintained by now also inversely tilted columns and reminiscent of zonal flow profiles observed for the ice giants. Analysis of the main force balance and parameter studies further foster these results.

We discuss here a purely hydrodynamical mechanism to invert the sign of the kinetic helicity, which plays a key role in determining the direction of propagation of cyclical magnetism in most models of dynamo action by rotating convection. Such propagation provides a prominent, and puzzling constraint on dynamo models. In the Sun, active regions emerge first at mid-latitudes, then appear nearer the equator over the course of a cycle, but most previous global-scale dynamo simulations have exhibited poleward propagation (if they were cyclical at all). Here, we highlight some simulations in which the direction of propagation of dynamo waves is altered primarily by an inversion of the kinetic helicity throughout much of the interior, rather than by changes in the differential rotation. This tends to occur in cases with a low Prandtl number and internal heating, in regions where the local density gradient is relatively small. We analyse how this inversion arises, and contrast it to the case of convection that is either highly columnar (i.e., rapidly rotating) or locally very stratified; in both of those situations, the typical profile of kinetic helicity (negative throughout most of the northern hemisphere) instead prevails.

We study rotating thermal convection in spherical shells. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-5}$. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-B\'enard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Furthermore, we explore the influence of the magnetic field on the {\em non-zonal} flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at $E=10^{-5}$ magnetic field significantly affects the convection and a maximum gain of about 30\% (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at $E=10^{-5}$ is probably in a `magnetostrophic' regime.

Despite the lack of a shear-rich tachocline region low-mass fully convective stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three dimensional model of magnetic field generation in low-mass fully convective stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman-Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass fully convective stars.

We simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9$. In this geometry, curvature and radial variations of the gravitationnal acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from $\eta=r_i/r_o=0.2$ to $\eta=0.95$) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is comparable between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame, that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the bulk and boundary layer contributions to the thermal and viscous dissipation rates using numerical models with $\eta=0.6$ and a gravity proportional to $1/r^2$. We show that our spherical models are consistent with the predictions of Grossmann \& Lohse's (2000) theory and that $Nu(Ra)$ and $Re(Ra)$ scalings are in good agreement with plane layer results.

Recent spectropolarimetric surveys of main-sequence intermediate-mass stars have exhibited a dichotomy in the distribution of the observed magnetic field between the kG dipoles of Ap/Bp stars and the sub-Gauss magnetism of Vega and Sirius. We would like to test whether this dichotomy is linked to the stability versus instability of large-scale magnetic configurations in differentially rotating radiative zones. We computed the axisymmetric magnetic field obtained from the evolution of a dipolar field threading a differentially rotating shell. A full parameter study including various density profiles and initial and boundary conditions was performed with a 2D numerical code. We then focused on the ratio between the toroidal and poloidal components of the magnetic field and discuss the stability of the configurations dominated by the toroidal component using local stability criteria and insights from recent 3D numerical simulations. The numerical results and a simple model show that the ratio between the toroidal and the poloidal magnetic fields is highest after an Alfv\'en crossing time of the initial poloidal field. For high density contrasts, this ratio converges towards an asymptotic value that can thus be extrapolated to realistic stellar cases. We then consider the stability of the magnetic configurations to non-axisymmetric perturbations and find that configurations dominated by the toroidal component are likely to be unstable if the shear strength is significantly higher than the poloidal Alfv\'en frequency. An expression for the critical poloidal field below which magnetic fields are likely to be unstable is found and is compared to the lower bound of Ap/Bp magnetic fields.

The question of the origin and evolution of magnetic fields in stars possessing a radiative envelope, like the A-type stars, is still regarded as a challenge for stellar physics. Those zones are likely to be differentially rotating, which suggests that strong interactions between differential rotation and magnetic fields could be at play. We numerically compute the joint evolution of the magnetic and velocity fields in a 3D spherical shell starting from an initial profile for the poloidal magnetic field and differential rotation. The poloidal magnetic field is initially wound-up by the differential rotation to produce a toroidal field which becomes unstable. In the particular setup studied here where the differential rotation is dominant, the magneto-rotational instability is triggered. The growth rate of the instability depends mainly on the initial rotation rate, while the background state typically oscillates over a poloidal Alfv\'en time. We thus find that the axisymmetric magnetic configuration is strongly modified by the instability only if the ratio between the poloidal Alfv\'en frequency and the rotation rate is sufficiently small. An enhanced transport of angular momentum is found in the most unstable cases: the typical time to flatten the rotation profile is then much faster than the diffusion time scale. We conclude that the magneto-rotational instability is always favored (over the Tayler instability) in unstratified spherical shells when an initial poloidal field is sheared by a sufficiently strong cylindrical differential rotation. A possible application to the magnetic desert observed among A stars is given. We argue that the dichotomy between stars exhibiting strong axisymmetric fields (Ap stars) and those harboring a sub-Gauss magnetism could be linked to the threshold for the instability.

Observations of cool stars reveal dark spot-like features on their surfaces. Compared to sunspots, starspots can be bigger or cover a larger fraction of the stellar surface. While sunspots appear only at low latitudes, starspots are also found in polar regions, in particular on rapidly rotating stars. Sunspots are believed to result from the eruption of magnetic flux-tubes rising from the deep interior of the Sun. The strong magnetic field locally reduces convective heat transport to the solar surface. Such flux-tube models have also been invoked to explain starspot properties. However, these models use several simplifications and so far the generation of either sunspots or starspots has not been demonstrated in a self-consistent simulation of stellar magnetic convection. Here we show that direct numerical simulations of a distributed dynamo operating in a density-stratified rotating spherical shell can spontaneously generate cool spots. Convection in the interior of the model produces a large scale magnetic field which interacts with near surface granular convection leading to strong concentrations of magnetic flux and formation of starspots. Prerequisites for the formation of sizeable high-latitude spots in the model are sufficiently strong density stratification and rapid rotation. Our model presents an alternate mechanism for starspot formation by distributed dynamo action.

Spacecraft data reveal a very Earth-like Jovian magnetic field. This is surprising since numerical simulations have shown that the vastly different interiors of terrestrial and gas planets can strongly affect the internal dynamo process. Here we present the first numerical dynamo that manages to match the structure and strength of the observed magnetic field by embracing the newest models for Jupiter's interior. Simulated dynamo action primarily occurs in the deep high electrical conductivity region while zonal flows are dynamically constrained to a strong equatorial jet in the outer envelope of low conductivity. Our model reproduces the structure and strength of the observed global magnetic field and predicts that secondary dynamo action associated to the equatorial jet produces banded magnetic features likely observable by the Juno mission. Secular variation in our model scales to about 2000 nT per year and should also be observable during the one year nominal mission duration.

The surface winds of Jupiter and Saturn are primarily zonal. Each planet exhibits strong prograde equatorial flow flanked by multiple alternating zonal winds at higher latitudes. The depth to which these flows penetrate has long been debated and is still an unsolved problem. Previous rotating convection models that obtained multiple high latitude zonal jets comparable to those on the giant planets assumed an incompressible (Boussinesq) fluid, which is unrealistic for gas giant planets. Later models of compressible rotating convection obtained only few high latitude jets which were not amenable to scaling analysis. Here we present 3-D numerical simulations of compressible convection in rapidly-rotating spherical shells. To explore the formation and scaling of high-latitude zonal jets, we consider models with a strong radial density variation and a range of Ekman numbers, while maintaining a zonal flow Rossby number characteristic of Saturn. All of our simulations show a strong prograde equatorial jet outside the tangent cylinder. At low Ekman numbers several alternating jets form in each hemisphere inside the tangent cylinder. To analyse jet scaling of our numerical models and of Jupiter and Saturn, we extend Rhines scaling based on a topographic $\beta$-parameter, which was previously applied to an incompressible fluid in a spherical shell, to compressible fluids. The jet-widths predicted by this modified Rhines length are found to be in relatively good agreement with our numerical model results and with cloud tracking observations of Jupiter and Saturn.

Stellar magnetic field measurements obtained from spectropolarimetry offer key data for activity and dynamo studies, and we present the results of a major high-resolution spectropolarimetric Bcool project magnetic snapshot survey of 170 solar-type stars from observations with the Telescope Bernard Lyot and the Canada-France-Hawaii Telescope. For each target star a high signal-to-noise circularly polarised Stokes V profile has been obtained using Least-Squares Deconvolution, and used to detect surface magnetic fields and measure the corresponding mean surface longitudinal magnetic field ($B_{l}$). Chromospheric activity indicators were also measured. Surface magnetic fields were detected for 67 stars, with 21 of these stars classified as mature solar-type stars, a result that increases by a factor of four the number of mature solar-type stars on which magnetic fields have been observed. In addition, a magnetic field was detected for 3 out of 18 of the subgiant stars surveyed. For the population of K-dwarfs the mean value of $B_{l}$ ($|B_{l}|_{mean}$) was also found to be higher (5.7 G) than $|B_{l}|_{mean}$ measured for the G-dwarfs (3.2 G) and the F-dwarfs (3.3 G). For the sample as a whole $|B_{l}|_{mean}$ increases with rotation rate and decreases with age, and the upper envelope for $|B_{l}|$ correlates well with the observed chromospheric emission. Stars with a chromospheric S-index greater than about 0.2 show a high magnetic field detection rate and so offer optimal targets for future studies. This survey constitutes the most extensive spectropolarimetric survey of cool stars undertaken to date, and suggests that it is feasible to pursue magnetic mapping of a wide range of moderately active solar-type stars to improve understanding of their surface fields and dynamos.

Stellar differential rotation can be separated into two main regimes: solar-like when the equator rotates faster than the poles and anti-solar when the polar regions rotate faster than the equator. We investigate the transition between these two regimes with 3-D numerical simulations of rotating spherical shells. We conduct a systematic parameter study which also includes models from different research groups. We find that the direction of the differential rotation is governed by the contribution of the Coriolis force in the force balance, independently of the model setup (presence of a magnetic field, thickness of the convective layer, density stratification). Rapidly-rotating cases with a small Rossby number yield solar-like differential rotation, while weakly-rotating models sustain anti-solar differential rotation. Close to the transition, the two kinds of differential rotation are two possible bistable states. This study provides theoretical support for the existence of anti-solar differential rotation in cool stars with large Rossby numbers.

Dynamos operating in the interiors of rapidly rotating planets and low-mass stars might belong to a similar category where rotation plays a vital role. We quantify this similarity using scaling laws. We analyse direct numerical simulations of Boussinesq and anelastic spherical shell dynamos. These dynamos represent simplified models which span from Earth-like planets to rapidly rotating low-mass stars. We find that magnetic field and velocity in these dynamos are related to the available buoyancy power via a simple power law which holds over wide variety of control parameters.

Observations of active M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. We detail the analogy between some anelastic dynamo simulations and spectropolarimetric observations of 23 M stars. In numerical models, the relative contribution of inertia and Coriolis force in the global force balance -estimated by the so-called local Rossby number- is known to have a strong impact on the magnetic field geometry. We discuss the relevance of this parameter in setting the large-scale magnetic field of M dwarfs.

The observed surface dynamics of Jupiter and Saturn is dominated by a banded system of zonal winds. Their depth remains unclear but they are thought to be confined to the very outer envelopes where hydrogen remains molecular and the electrical conductivity is small. The dynamo maintaining the dipole-dominated magnetic fields of both gas giants likely operates in the deeper interior where hydrogen assumes a metallic state. Here, we present numerical simulations that attempt to model both the zonal winds and the interior dynamo action in an integrated approach. Using the anelastic version of the MHD code MagIC, we explore the effects of density stratification and radial electrical conductivity variation. The electrical conductivity is mostly assumed to remain constant in the thicker inner metallic region and it decays exponentially towards the outer boundary throughout the molecular envelope. Our results show that the combination of stronger density stratification and weaker conducting outer layer is essential for reconciling dipole dominated dynamo action and a fierce equatorial zonal jet. Previous simulations with homogeneous electrical conductivity show that both are merely exclusive, with solutions either having strong zonal winds and multipolar magnetic fields or weak zonal winds and dipole-dominated magnetic fields. All jets tend to be geostrophic and therefore reach right through the convective shell in our simulations. The particular setup explored here allows a strong equatorial jet to remain confined to the weaker conducting outer region where it does not interfere with the deeper seated dynamo action. The flanking mid to high latitude jets, on the other hand, have to remain faint to yield a strongly dipolar magnetic field. The fiercer jets on Jupiter and Saturn only seem compatible with the observed dipolar fields when they remain confined to a weaker conducting outer layer.

Numerical dynamo models always employ parameter values that differ by orders of magnitude from the values expected in natural objects. However, such models have been successful in qualitatively reproducing properties of planetary and stellar dynamos. This qualitative agreement fuels the idea that both numerical models and astrophysical objects may operate in the same asymptotic regime of dynamics. This can be tested by exploring the scaling behavior of the models. For convection-driven incompressible spherical shell dynamos with constant material properties, scaling laws had been established previously that relate flow velocity and magnetic field strength to the available power. Here we analyze 273 direct numerical simulations using the anelastic approximation, involving also cases with radius-dependent magnetic, thermal and viscous diffusivities. These better represent conditions in gas giant planets and low-mass stars compared to Boussinesq models. Our study provides strong support for the hypothesis that both mean velocity and mean magnetic field strength scale as a function of power generated by buoyancy forces in the same way for a wide range of conditions.

The surface zonal winds observed in the giant planets form a complex jet pattern with alternating prograde and retrograde direction. While the main equatorial band is prograde on the gas giants, both ice giants have a pronounced retrograde equatorial jet. We use three-dimensional numerical models of compressible convection in rotating spherical shells to explore the properties of zonal flows in different regimes where either rotation or buoyancy dominates the force balance. We conduct a systematic parameter study to quantify the dependence of zonal flows on the background density stratification and the driving of convection. We find that the direction of the equatorial zonal wind is controlled by the ratio of buoyancy and Coriolis force. The prograde equatorial band maintained by Reynolds stresses is found in the rotation-dominated regime. In cases where buoyancy dominates Coriolis force, the angular momentum per unit mass is homogenised and the equatorial band is retrograde, reminiscent to those observed in the ice giants. In this regime, the amplitude of the zonal jets depends on the background density contrast with strongly stratified models producing stronger jets than comparable weakly stratified cases. Furthermore, our results can help to explain the transition between solar-like and "anti-solar" differential rotations found in anelastic models of stellar convection zones. In the strongly stratified cases, we find that the leading order force balance can significantly vary with depth (rotation-dominated inside and buoyancy-dominated in a thin surface layer). This so-called "transitional regime" has a visible signature in the main equatorial jet which shows a pronounced dimple where flow amplitudes notably decay towards the equator. A similar dimple is observed on Jupiter, which suggests that convection in the planet interior could possibly operate in this regime.

Numerical simulations of convection driven rotating spherical shell dynamos have often been performed with rigid boundary conditions, as is appropriate for the metallic cores of terrestrial planets. Free-slip boundaries are more appropriate for dynamos in other astrophysical objects, such as gas-giants or stars. Using a set of 57 direct numerical simulations, we investigate the effect of free-slip boundary conditions on the scaling properties of heat flow, flow velocity and magnetic field strength and compare it with earlier results for rigid boundaries. We find that the nature of the mechanical boundary condition has only a minor influence on the scaling laws. We also find that although dipolar and multipolar dynamos exhibit approximately the same scaling exponents, there is an offset in the scaling pre-factors for velocity and magnetic field strength. We argue that the offset can be attributed to the differences in the zonal flow contribution between dipolar and multipolar dynamos.

Context: observations of rapidly rotating M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. In dynamo models, the relative importance of inertia in the force balance -- quantified by the local Rossby number -- is known to have a strong impact on the magnetic field geometry. Aims: we aim to assess the relevance of the local Rossby number in controlling the large-scale magnetic field geometry of M dwarfs. Methods: we explore the similarities between anelastic dynamo models in spherical shells and observations of active M-dwarfs, focusing on field geometries derived from spectropolarimetric studies. To do so, we construct observation-based quantities aimed to reflect the diagnostic parameters employed in numerical models. Results: the transition between dipole-dominated and multipolar large-scale fields in early to mid M dwarfs is tentatively attributed to a Rossby number threshold. We interpret late M dwarfs magnetism to result from a dynamo bistability occurring at low Rossby number. By analogy with numerical models, we expect different amplitudes of differential rotation on the two dynamo branches.

Context: dynamo action in giant planets and rapidly rotating stars leads to a broad variety of magnetic field geometries including small scale multipolar and large scale dipole-dominated topologies. Previous dynamo models suggest that solutions become multipolar once inertia becomes influential. Being tailored for terrestrial planets, most of these models neglected the background density stratification. Aims: we investigate the influence of the density stratification on convection-driven dynamo models. Methods: three-dimensional nonlinear simulations of rapidly rotating spherical shells are employed using the anelastic approximation to incorporate density stratification. A systematic parametric study for various density stratifications and Rayleigh numbers allows to explore the dependence of the magnetic field topology on these parameters. Results: anelastic dynamo models tend to produce a broad range of magnetic field geometries that fall on two distinct branches with either strong dipole-dominated or weak multipolar fields. As long as inertia is weak, both branches can coexist but the dipolar branch vanishes once inertia becomes influential. The dipolar branch also vanishes for stronger density stratifications. The reason is the concentration of the convective columns in a narrow region close to the outer boundary equator, a configuration that favors non-axisymmetric solutions. In multipolar solutions, zonal flows can become significant and participate in the toroidal field generation. Parker dynamo waves may then play an important role close to onset of dynamo action leading to a cyclic magnetic field behavior. Conclusion: Our simulations also suggest that the fact that late M dwarfs have dipolar or multipolar magnetic fields can be explained in two ways. They may differ either by the relative influence of inertia or fall into the regime where both types of solutions coexist.

The banded structures observed on the surfaces of the gas giants are associated with strong zonal winds alternating in direction with latitude. We use three-dimensional numerical simulations of compressible convection in the anelastic approximation to explore the properties of zonal winds in rapidly rotating spherical shells. Since the model is restricted to the electrically insulating outer envelope, we therefore neglect magnetic effects. A systematic parametric study for various density scaleheights and Rayleigh numbers allows to explore the dependence of convection and zonal jets on these parameters and to derive scaling laws. While the density stratification affects the local flow amplitude and the convective scales, global quantities and zonal jets properties remain fairly independent of the density stratification. The zonal jets are maintained by Reynolds stresses, which rely on the correlation between zonal and cylindrically radial flow components. The gradual loss of this correlation with increasing supercriticality hampers all our simulations and explains why the additional compressional source of vorticity hardly affects zonal flows. All these common features may explain why previous Boussinesq models were already successful in reproducing the morphology of zonal jets in gas giants.

Aims. We aim to investigate the long-term temporal evolution of the magnetic field of the solar-type star \xi Bootis A, both from direct magnetic field measurements and from the simultaneous estimate of indirect activity indicators. Methods. We obtained seven epochs of high-resolution, circularly-polarized spectra from the NARVAL spectropolarimeter between 2007 and 2011, for a total of 76 spectra. Using approximately 6,100 photospheric spectral lines covering the visible domain, we employed a cross-correlation procedure to compute a mean polarized line profile from each spectrum. The large-scale photospheric magnetic field of the star was then modelled by means of Zeeman-Doppler Imaging, allowing us to follow the year-to-year evolution of the reconstructed magnetic topology. Simultaneously, we monitored the width of several magnetically sensitive spectral lines, the radial velocity, the line asymmetry of intensity line profiles, and the chromospheric emission in the cores of the Ca II H and Halpha lines. Results. During the highest observed activity states, in 2007 and 2011, the large-scale field of \xi Boo A is almost completely axisymmetric and is dominated by its toroidal component. The magnetic topologies reconstructed for these activity maxima are very similar, suggesting a form of short cyclicity in the large-scale field distribution. Correlated temporal evolution, due to both rotational modulation and seasonal variability, is observed between the Ca II emission, the Halpha emission and the width of magnetically sensitive lines. When measurable, the differential rotation reveals a strong latitudinal shear in excess of 0.2 rad/d.

In cold Cepheids close to the red edge of the classical instability strip, a strong coupling between the stellar pulsations and the surface convective motions occurs. This coupling is by now poorly described by 1-D models of convection, the so-called "time-dependent convection models" (TDC). The intrinsic weakness of such models comes from the large number of unconstrained free parameters entering in the description of turbulent convection. A way to overcome these limits is to compute two-dimensional direct simulations (DNS), in which all the nonlinearities are correctly solved. Two-dimensional DNS of the convection-pulsation coupling are presented here. In an appropriate parameter regime, convective motions can actually quench the radial pulsations of the star, as suspected in Cepheids close to the red edge of the instability strip. These nonlinear simulations can also be used to determine the limits and the relevance of the TDC models.

Context: In Cepheids close to the red edge of the classical instability strip, a coupling occurs between the acoustic oscillations and the convective motions close to the surface.The best topical models that account for this coupling rely on 1-D time-dependent convection (TDC) formulations. However, their intrinsic weakness comes from the large number of unconstrained free parameters entering in the description of turbulent convection. Aims: We compare two widely used TDC models with the first two-dimensional nonlinear direct numerical simulations (DNS) of the convection-pulsation coupling in which the acoustic oscillations are self-sustained by the kappa-mechanism. Methods: The free parameters appearing in the Stellingwerf and Kuhfuss TDC recipes are constrained using a chi2-test with the time-dependent convective flux that evolves in nonlinear simulations of highly-compressible convection with kappa-mechanism. Results: This work emphasises some inherent limits of TDC models, that is, the temporal variability and non-universality of their free parameters. More importantly, within these limits, Stellingwerf's formalism is found to give better spatial and temporal agreements with the nonlinear simulation than Kuhfuss's one. It may therefore be preferred in 1-D TDC hydrocodes or stellar evolution codes.

Context: we study the convection-pulsation coupling that occurs in cold Cepheids close to the red edge of the classical instability strip. In these stars, the surface convective zone is supposed to stabilise the radial oscillations excited by the kappa-mechanism. Aims: we study the influence of the convective motions onto the amplitude and the nonlinear saturation of acoustic modes excited by kappa-mechanism. We are interested in determining the physical conditions needed to lead to a quenching of oscillations by convection. Methods: we compute two-dimensional nonlinear simulations (DNS) of the convection-pulsation coupling, in which the oscillations are sustained by a continuous physical process: the kappa-mechanism. Thanks to both a frequential analysis and a projection of the physical fields onto an acoustic subspace, we study how the convective motions affect the unstable radial oscillations. Results: depending on the initial physical conditions, two main behaviours are obtained: (i) either the unstable fundamental acoustic mode has a large amplitude, carries the bulk of the kinetic energy and shows a nonlinear saturation similar to the purely radiative case; (ii) or the convective motions affect significantly the mode amplitude that remains very weak. In this second case, convection is quenching the acoustic oscillations. We interpret these discrepancies in terms of the difference in density contrast: larger stratification leads to smaller convective plumes that do not affect much the purely radial modes, while large-scale vortices may quench the oscillations.

Phase-resolved observations of the solar-type star Ksi Bootis A were obtained using the NARVAL spectropolarimeter at the Telescope Bernard Lyot (Pic du Midi, France) during years 2007, 2008, 2009 and 2010. The data sets enable us to study both the rotational modulation and the long-term evolution of various magnetic and activity tracers. Here, we focus on the large-scale photospheric magnetic field (reconstructed by Zeeman-Doppler Imaging), the Zeeman broadening of the FeI 846.84 nm magnetic line, and the chromospheric CaII H and H alpha emission.

Magnetic fields of cool stars can be directly investigated through the study of the Zeeman effect on photospheric spectral lines using several approaches. With spectroscopic measurement in unpolarised light, the total magnetic flux averaged over the stellar disc can be derived but very little information on the field geometry is available. Spectropolarimetry provides a complementary information on the large-scale component of the magnetic topology. With Zeeman-Doppler Imaging (ZDI), this information can be retrieved to produce a map of the vector magnetic field at the surface of the star, and in particular to assess the relative importance of the poloidal and toroidal components as well as the degree of axisymmetry of the field distribution. The development of high-performance spectropolarimeters associated with multi-lines techniques and ZDI allows us to explore magnetic topologies throughout the Hertzsprung-Russel diagram, on stars spanning a wide range of mass, age and rotation period. These observations bring novel constraints on magnetic field generation by dynamo effect in cool stars. In particular, the study of solar twins brings new insight on the impact of rotation on the solar dynamo, whereas the detection of strong and stable dipolar magnetic fields on fully convective stars questions the precise role of the tachocline in this process.

A strong coupling between convection and pulsations is known to play a major role in the disappearance of unstable modes close to the red edge of the classical Cepheid instability strip. As mean-field models of time-dependent convection rely on weakly-constrained parameters, we tackle this problem by the means of 2-D Direct Numerical Simulations (DNS) of kappa-mechanism with convection. Using a linear stability analysis, we first determine the physical conditions favourable to the kappa-mechanism to occur inside a purely-radiative layer. Both the instability strips and the nonlinear saturation of unstable modes are then confirmed by the corresponding DNS. We next present the new simulations with convection, where a convective zone and the driving region overlap. The coupling between the convective motions and acoustic modes is then addressed by using projections onto an acoustic subspace.

We present a purely-radiative hydrodynamic model of the kappa-mechanism that sustains radial oscillations in Cepheid variables. We determine the physical conditions favourable for the kappa-mechanism to occur by the means of a configurable hollow in the radiative conductivity profile. By starting from these most favourable conditions, we complete nonlinear direct numerical simulations (DNS) and compare them with the results given by a linear-stability analysis of radial modes. We find that well-defined instability strips are generated by changing the location and shape of the conductivity hollow. For a given position in the layer, the hollow amplitude and width stand out as the key parameters governing the appearance of unstable modes driven by the kappa-mechanism. The DNS confirm both the growth rates and structures of the linearly-unstable modes. The nonlinear saturation that arises is produced by intricate couplings between the excited fundamental mode and higher damped overtones. These couplings are measured by projecting the DNS fields onto an acoustic subspace built from regular and adjoint eigenvectors and a 2:1 resonance is found to be responsible for the saturation of the kappa-mechanism instability.

Context: We study the kappa-mechanism that excites radial oscillations in Cepheid variables. Aims: We address the mode couplings that manages the nonlinear saturation of the instability in direct numerical simulations (DNS). Methods: We project the DNS fields onto an acoustic subspace built from the regular and adjoint eigenvectors that are solutions to the linear-oscillations equations. Results: We determine the time evolution of both the amplitude and kinetic energy of each mode that propagates in the DNS. More than 98% of the total kinetic energy is contained in two modes that correspond to the linearly-unstable fundamental mode and the linearly-stable second overtone. Because the eigenperiods ratio is close to 1/2, we discover that the nonlinear saturation is due to a 2:1 resonance between these two modes. An interesting application of this result concerns the reproduction of Hertzsprung's progression observed in Bump Cepheids.

We present here additional results of a spectropolarimetric survey of a small sample of stars ranging from spectral type M0 to M8 aimed at investigating observationally how dynamo processes operate in stars on both sides of the full convection threshold (spectral type M4). The present paper focuses on early M stars (M0--M3), i.e. above the full convection threshold. Applying tomographic imaging techniques to time series of rotationally modulated circularly polarised profiles collected with the NARVAL spectropolarimeter, we determine the rotation period and reconstruct the large-scale magnetic topologies of 6 early M dwarfs. We find that early-M stars preferentially host large-scale fields with dominantly toroidal and non-axisymmetric poloidal configurations, along with significant differential rotation (and long-term variability); only the lowest-mass star of our subsample is found to host an almost fully poloidal, mainly axisymmetric large-scale field ressembling those found in mid-M dwarfs. This abrupt change in the large-scale magnetic topologies of M dwarfs (occuring at spectral type M3) has no related signature on X-ray luminosities (measuring the total amount of magnetic flux); it thus suggests that underlying dynamo processes become more efficient at producing large-scale fields (despite producing the same flux) at spectral types later than M3. We suspect that this change relates to the rapid decrease in the radiative cores of low-mass stars and to the simultaneous sharp increase of the convective turnover times (with decreasing stellar mass) that models predict to occur at M3; it may also be (at least partly) responsible for the reduced magnetic braking reported for fully-convective stars.

We present in this paper the first results of a spectropolarimetric analysis of a small sample (~ 20) of active stars ranging from spectral type M0 to M8, which are either fully-convective or possess a very small radiative core. This study aims at providing new constraints on dynamo processes in fully-convective stars. The present paper focuses on 5 stars of spectral type ~M4, i.e. with masses close to the full convection threshold (~ 0.35 Msun), and with short rotational periods. Tomographic imaging techniques allow us to reconstruct the surface magnetic topologies from the rotationally modulated time-series of circularly polarised profiles. We fnd that all stars host mainly axisymmetric large-scale poloidal fields. Three stars were observed at two different epochs separated by ~1 yr; we find the magnetic topologies to be globally stable on this timescale. We also provide an accurate estimation of the rotational period of all stars, thus allowing us to start studying how rotation impacts the large-scale magnetic field.

Context: Hydrodynamical model of the kappa-mechanism in a purely radiative case. Aims: First, to determine the physical conditions propitious to kappa-mechanism in a layer with a configurable conductivity hollow and second, to perform the (nonlinear) direct numerical simulations (DNS) from the most favourable setups. Methods: A linear stability analysis applied to radial modes using a spectral solver and DNS thanks to a high-order finite difference code are compared. Results: Changing the hollow properties (location and shape) lead to well-defined instability strips. For a given position in the layer, the amplitude and width of the hollow appear to be the key parameters to get unstable modes driven by kappa-mechanism. The DNS achieved from these more auspicious configurations confirm the growth rates as well as structures of linearly unstable modes. The nonlinear saturation follows through intricate couplings between the excited fundamental mode and higher damped overtones.