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.

Dynamo action in fully convective stars is a debated issue that also questions our understanding of magnetic field generation in partly convective Sun-like stars. During the past few years, spectropolari- metric observations have demonstrated that fully convective objects are able to trigger strong large-scale and long-lived magnetic fields. We present here the first spectropolarimetric study of a sample of active late M dwarfs (M5-M8) carried out with ESPaDOnS@CFHT. It reveals the co-existence of two distinct types of magnetism among stars having similar masses and rotation rates. A possible explanation for this unexpected discovery is the existence of two dynamo branches in this parameter regime, we discuss here the possible identification with the weak vs strong field bistability predicted for the geodynamo.

We review some of the recent progress on modeling planetary and stellar dynamos. Particular attention is given to the dynamo mechanisms and the resulting properties of the field. We present direct numerical simulations using a simple Boussinesq model. These simulations are interpreted using the classical mean-field formalism. We investigate the transition from steady dipolar to multipolar dynamo waves solutions varying different control parameters, and discuss the relevance to stellar magnetic fields. We show that owing to the role of the strong zonal flow, this transition is hysteretic. In the presence of stress-free boundary conditions, the bistability extends over a wide range of parameters.

Magnetic fields of low-mass stars and planets are thought to originate from self-excited dynamo action in their convective interiors. Observations reveal a variety of field topologies ranging from large-scale, axial dipole to more structured magnetic fields. In this article, we investigate more than 70 three-dimensional, self-consistent dynamo models obtained by direct numerical simulations. The control parameters, the aspect ratio and the mechanical boundary conditions have been varied to build up this sample of models. Both, strongly dipolar and multipolar models have been obtained. We show that these dynamo regimes can in general be distinguished by the ratio of a typical convective length scale to the Rossby radius. Models with a predominantly dipolar magnetic field were obtained, if the convective length scale is at least an order of magnitude larger than the Rossby radius. Moreover, we highlight the role of the strong shear associated with the geostrophic zonal flow for models with stress-free boundary conditions. In this case, the above transition disappears and is replaced by a region of bistability for which dipolar and multipolar dynamos co-exist. We interpret our results in terms of dynamo eigenmodes using the so-called test-field method. We can thus show that models in the dipolar regime are characterized by an isolated 'single mode'. Competing overtones become significant as the boundary to multipolar dynamos is approached. We discuss how these findings relate to previous models and to observations.

Observations of magnetism in very low mass stars recently made important progress, revealing characteristics that are now to be understood in the framework of dynamo theory. In parallel, there is growing evidence that dynamo processes in these stars share many similarities with planetary dynamos. We investigate the extent to which the weak \emph{vs} strong field bistability predicted for the geodynamo can apply to recent observations of two groups of very low mass fully-convective stars sharing similar stellar parameters but generating radically different types of magnetic fields. Our analysis is based on previously published spectropolarimetric and spectroscopic data. We argue that these can be interpreted in the framework of weak and strong field dynamos.

Context: Large-scale magnetic fields resulting from hydromagnetic dynamo action may differ substantially in their time dependence. Cyclic field variations, characteristic for the solar magnetic field, are often explained by an important omega-effect, i.e. by the stretching of field lines due to strong differential rotation. Aims: The dynamo mechanism of a convective, oscillatory dynamo model is investigated. Methods: We solve the MHD-equations for a conducting Boussinesq fluid in a rotating spherical shell. For a resulting oscillatory model, dynamo coefficients have been computed with the help of the so-called test-field method. Subsequently, these coefficients have been used in a mean-field calculation in order to explore the underlying dynamo mechanism. Results: Although the rather strong differential rotation present in this model influences the magnetic field, the omega-effect alone is not responsible for its cyclic time variation. If the omega-effect is suppressed, the resulting alpha^2-dynamo remains oscillatory. Surprisingly, the corresponding alpha-omega dynamo leads to a non-oscillatory magnetic field. Conclusions: The assumption of an alpha-omega mechanism does not explain the occurrence of magnetic cycles satisfactorily.

We investigate the nature of the dynamo bifurcation in a configuration applicable to the Earth's liquid outer core, i.e. in a rotating spherical shell with thermally driven motions. We show that the nature of the bifurcation, which can be either supercritical or subcritical or even take the form of isola (or detached lobes) strongly depends on the parameters. This dependence is described in a range of parameters numerically accessible (which unfortunately remains remote from geophysical application), and we show how the magnetic Prandtl number and the Ekman number control these transitions.

We show that a simple, modified version of the Magnetorotational Instability (MRI) can develop in the outer liquid core of the Earth, in the presence of a background shear. It requires either thermal wind, or a primary instability, such as convection, to drive a weak differential rotation within the core. The force balance in the Earth's core is very unlike classical astrophysical applications of the MRI (such as gaseous disks around stars). Here, the weak differential rotation in the Earth core yields an instability by its constructive interaction with the planet's much larger rotation rate. The resulting destabilising mechanism is just strong enough to counteract stabilizing resistive effects, and produce growth on geophysically interesting timescales. We give a simple physical explanation of the instability, and show that it relies on a force balance appropriate to the Earth's core, known as magnetostrophic balance.

The geo and solar magnetic fields have long been thought to be very different objects both in terms of spatial structure and temporal behavior. The recently discovered field structure of a fully convective star is more reminiscent of planetary magnetic fields than the Sun's magnetic field (Donati J.-F. et al., Science, 311 (2006) 633), despite the fact that the physical and chemical properties of these objects clearly differ. This observation suggests that a simple controlling parameter could be responsible for these different behaviors. We report here the results of three-dimensional simulations which show that varying the aspect ratio of the active dynamo region can yield sharp transition from Earth-like steady dynamos to Sun-like dynamo waves.

We present kinematic simulations of a galactic dynamo model based on the large scale differential rotation and the small scale helical fluctuations due to supernova explosions. We report for the first time direct numerical simulations of the full galactic dynamo using an unparameterized global approach. We argue that the scale of helicity injection is large enough to be directly resolved rather than parameterized. While the actual superbubble characteristics can only be approached, we show that numerical simulations yield magnetic structures which are close both to the observations and to the previous parameterized mean field models. In particular, the quadrupolar symmetry and the spiraling properties of the field are reproduced. Moreover, our simulations show that the presence of a vertical inflow plays an essential role to increase the magnetic growth rate. This observation could indicate an important role of the downward flow (possibly linked with galactic fountains) in sustaining galactic magnetic fields.

The study of dynamo action in astrophysical objects classically involves two timescales: the slow diffusive one and the fast advective one. We investigate the possibility of field amplification on an intermediate timescale associated with time dependent modulations of the flow. We consider a simple steady configuration for which dynamo action is not realised. We study the effect of time dependent perturbations of the flow. We show that some vanishing low frequency perturbations can yield exponential growth of the magnetic field on the typical time scale of oscillation. The dynamo mechanism relies here on a parametric instability associated with transient amplification by shear flows. Consequences on natural dynamos are discussed.

We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying conservation law which, in our formulation, results in a ``finite-surface'' scheme for the induction equation. This naturally leads to the well-known ``constrained transport'' method, with additional continuity requirement on the magnetic field representation. The second ingredient in the MUSCL scheme is the predictor step that ensures second order accuracy both in space and time. We explore specific constraints that the mathematical properties of the induction equations place on this predictor step, showing that three possible variants can be considered. We show that the most aggressive formulations (referred to as C-MUSCL and U-MUSCL) reach the same level of accuracy as the other one (referred to as Runge-Kutta), at a lower computational cost. More interestingly, these two schemes are compatible with the Adaptive Mesh Refinement (AMR) framework. It has been implemented in the AMR code RAMSES. It offers a novel and efficient implementation of a second order scheme for the induction equation. We have tested it by solving two kinematic dynamo problems in the low diffusion limit. The construction of this scheme for the induction equation constitutes a step towards solving the full MHD set of equations using an extension of our current methodology.