Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MRI at generating a dipole field is still unknown. To answer this question, we study the MRI dynamo in an idealized global model of a fast rotating protoneutron star with differential rotation. We perform 3D incompressible MHD simulations in spherical geometry with explicit diffusivities where the differential rotation is forced at the outer boundary. We vary the initial magnetic field and investigated different magnetic boundary conditions. These simulations were compared to local shearing box simulations. We obtain a self-sustained turbulent MRI-driven dynamo, whose saturated state is independent of the initial magnetic field. The MRI generates a strong turbulent magnetic field of $B \geq 2\times 10^{15}$ G and a non-dominant magnetic dipole, which represents systematically about $5\%$ of the averaged magnetic field strength. Interestingly, this dipole is tilted towards the equatorial plane. We find that local shearing box models can reproduce fairly well several characteristics of global MRI turbulence such as the kinetic and magnetic spectra. The turbulence is nonetheless more vigorous in the local models than in the global ones. Overall, our results support the ability of the MRI to form magnetar-like large-scale magnetic fields. They furthermore predict the presence of a stronger small-scale magnetic field. The resulting magnetic field could be important to power outstanding stellar explosions, such as superluminous supernovae and GRBs.

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.

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.

We consider dynamo action driven by three-dimensional rotating anelastic convection in a spherical shell. Motivated by the behaviour of the solar dynamo, we examine the interaction of hydromagnetic modes with different symmetries and demonstrate how complicated interactions between convection, differential rotation and magnetic fields may lead to modulation of the basic cycle. For some parameters, Type 1 modulation occurs by the transfer of energy between modes of different symmetries with little change in the overall amplitude, for other parameters, the modulation is of Type 2, where the amplitude is significantly affected (leading to grand minima in activity) without significant changes in symmetry. Most importantly, we identify the presence of "supermodulation" in the solutions, where the activity switches chaotically between Type 1 and Type 2 modulation, this is believed to be an important process in solar activity.

Observations of low-mass stars reveal a variety of magnetic field topologies ranging from large-scale, axial dipoles to more complex magnetic fields. At the same time, three-dimensional spherical simulations of convectively driven dynamos reproduce a similar diversity, which is commonly obtained either with Boussinesq models or with more realistic models based on the anelastic approximation, which take into account the variation of the density with depth throughout the convection zone. Nevertheless, a conclusion from different anelastic studies is that dipolar solutions seem more difficult to obtain as soon as substantial stratifications are considered. In this paper, we aim at clarifying this point by investigating in more detail the influence of the density stratification on dipolar dynamos. To that end, we rely on a systematic parameter study that allows us to clearly follow the evolution of the stability domain of the dipolar branch as the density stratification is increased. The impact of the density stratification both on the dynamo onset and the dipole collapse is discussed and compared to previous Boussinesq results. Furthermore, our study indicates that the loss of the dipolar branch does not ensue from a specific modification of the dynamo mechanisms related to the background stratification, but could instead result from a bias as our observations naturally favour a certain domain in the parameter space characterized by moderate values of the Ekman number, owing to current computational limitations. Moreover, we also show that the critical magnetic Reynolds number of the dipolar branch is scarcely modified by the increase of the density stratification, which provides an important insight into the global understanding of the impact of the density stratification on the stability domain of the dipolar dynamo branch.

Numerical modelling of convection driven dynamos in the Boussinesq approximation revealed fundamental characteristics of the dynamo-generated magnetic fields and the fluid flow. Because these results were obtained for an incompressible fluid, their validity for gas planets and stars remains to be assessed. A common approach is to take some density stratification into account with the so-called anelastic approximation. The validity of previous results obtained in the Boussinesq approximation is tested for anelastic models. We point out and explain specific differences between both types of models, in particular with respect to the field geometry and the field strength, but we also compare scaling laws for the velocity amplitude, the magnetic dissipation time, and the convective heat flux. Our investigation is based on a systematic parameter study of spherical dynamo models in the anelastic approximation. We make use of a recently developed numerical solver and provide results for the test cases of the anelastic dynamo benchmark. The dichotomy of dipolar and multipolar dynamos identified in Boussinesq simulations is also present in our sample of anelastic models. Dipolar models require that the typical length scale of convection is an order of magnitude larger than the Rossby radius. However, the distinction between both classes of models is somewhat less explicit than in previous studies. This is mainly due to two reasons: we found a number of models with a considerable equatorial dipole contribution and an intermediate overall dipole field strength. Furthermore, a large density stratification may hamper the generation of dipole dominated magnetic fields. Previously proposed scaling laws, such as those for the field strength, are similarly applicable to anelastic models. It is not clear, however, if this consistency necessarily implies similar dynamo processes in both settings.

Context. Three-dimensional spherical dynamo simulations carried out within the framework of the anelastic approximation have revealed that the established distinction between dipolar and multipolar dynamos tends to be less clear than it was in Boussinesq studies. This result was first interpreted as a direct consequence of the existence of a larger number of models with a high equatorial dipole contribution, together with an intermediate dipole field strength. However, this finding has not been clearly related to specific changes that would have been introduced by the use of the anelastic approximation. Aims. In this paper, we primarily focus on the effects of choosing of a different mass distribution. Indeed, it is likely to have as large consequences as would taking a stratified reference state into account, especially when comparing our results to previous Boussinesq studies. Methods. Our investigation is based on a systematic parameter study of weakly stratified anelastic dynamo models. Results. We show that the tendencies highlighted in previous anelastic dynamo simulations are already present in the Boussinesq limit. Thus they cannot be systematically related to anelastic effects. Actually, a central mass distribution can result in changes in the magnetic field topology that are mainly due to the concentration of convective cells close to the inner sphere.