Solar dynamo models with realistic internal rotation

Solar dynamo models based on differential rotation inferred from helioseismology tend to produce rather strong magnetic activity at high solar latitudes, in contrast to the observed fact that sunspots appear at low latitudes. We show that a meridional circulation penetrating below the tachocline can solve this problem.

Quenching of Meridional Circulation in Flux Transport Dynamo Models

Guided by the recent observational result that the meridional circulation of the Sun becomes weaker at the time of the sunspot maximum, we have included a parametric quenching of the meridional circulation in solar dynamo models such that the meridional circulation becomes weaker when the magnetic field at the base of the convection zone is stronger. We find that a flux transport solar dynamo tends to become unstable on including this quenching of meridional circulation if the diffusivity in the convection zone is less than about 2x10(11) cm(2) s(-1). The quenching of alpha, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong alpha-quenching. For dynamo models with high diffusivity, the quenching of meridional circulation does not produce a large effect and the dynamo remains stable. We present a solar-like solution from a dynamo model with diffusivity 2.8x10(12) cm(2) s(-1) in which the quenching of meridional circulation makes the meridional circulation vary periodically with solar cycle as observed and does not have any other significant effect on the dynamo.

Dynamo models of grand minima

Since a universally accepted dynamo model of grand minima does not exist at the present time, we concentrate on the physical processes which may be behind the grand minima. After summarizing the relevant observational data, we make the point that, while the usual sources of irregularities of solar cycles may be sufficient to cause a grand minimum, the solar dynamo has to operate somewhat differently from the normal to bring the Sun out of the grand minimum. We then consider three possible sources of irregularities in the solar dynamo: (i) nonlinear effects; (ii) fluctuations in the poloidal field generation process; (iii) fluctuations in the meridional circulation. We conclude that (i) is unlikely to be the cause behind grand minima, but a combination of (ii) and (iii) may cause them. If fluctuations make the poloidal field fall much below the average or make the meridional circulation significantly weaker, then the Sun may be pushed into a grand minimum.

THE ROLE OF DIFFUSIVITY QUENCHING IN FLUX-TRANSPORT DYNAMO MODELS

In the nonlinear phase of a dynamo process, the back-reaction of the magnetic field upon the turbulent motion results in a decrease of the turbulence level and therefore in a suppression of both the magnetic field amplification (the alpha-quenching effect) and the turbulent magnetic diffusivity (the eta-quenching effect). While the former has been widely explored, the effects of eta-quenching in the magnetic field evolution have rarely been considered. In this work, we investigate the role of the suppression of diffusivity in a flux-transport solar dynamo model that also includes a nonlinear alpha-quenching term. Our results indicate that, although for alpha-quenching the dependence of the magnetic field amplification with the quenching factor is nearly linear, the magnetic field response to eta-quenching is nonlinear and spatially nonuniform. We have found that the magnetic field can be locally amplified in this case, forming long-lived structures whose maximum amplitude can be up to similar to 2.5 times larger at the tachocline and up to similar to 2 times larger at the center of the convection zone than in models without quenching. However, this amplification leads to unobservable effects and to a worse distribution of the magnetic field in the butterfly diagram. Since the dynamo cycle period increases when the efficiency of the quenching increases, we have also explored whether the eta-quenching can cause a diffusion-dominated model to drift into an advection-dominated regime. We have found that models undergoing a large suppression in eta produce a strong segregation of magnetic fields that may lead to unsteady dynamo-oscillations. On the other hand, an initially diffusion-dominated model undergoing a small suppression in eta remains in the diffusion-dominated regime.

An elementary introduction to solar dynamo theory

The cyclically varying magnetic field of the Sun is believed to be produced by the hydromagnetic dynamo process. We first summarize the relevant observational data pertaining to sunspots and solar cycle. Then we review the basic principles of MHD needed to develop the dynamo theory. This is followed by a discussion how bipolar sunspots form due to magnetic buoyancy of flux tubes formed at the base of the solar convection zone. Following this, we come to the heart of dynamo theory. After summarizing the basic ideas of a turbulent dynamo and the basic principles of its mean field formulation, we present the famous dynamo wave solution, which was supposed to provide a model for the solar cycle. Finally we point out how a flux transport dynamo can circumvent some of the difficulties associated with the older dynamo models.

Importance of meridional circulation i flux transport dynamo: the possibility of a maunder-like grand minimum

Meridional circulation is an important ingredient in flux transport dynamo models. We have studied its importance on the period, the amplitude of the solar cycle, and also in producing Maunder-like grand minima in these models. First, we model the periods of the last 23 sunspot cycles by varying the meridional circulation speed. If the dynamo is in a diffusion-dominated regime, then we find that most of the cycle amplitudes also get modeled up to some extent when we model the periods. Next, we propose that at the beginning of the Maunder minimum the amplitude of meridional circulation dropped to a low value and then after a few years it increased again. Several independent studies also favor this assumption. With this assumption, a diffusion-dominated dynamo is able to reproduce many important features of the Maunder minimum remarkably well. If the dynamo is in a diffusion-dominated regime, then a slower meridional circulation means that the poloidal field gets more time to diffuse during its transport through the convection zone, making the dynamo weaker. This consequence helps to model both the cycle amplitudes and the Maunder-like minima. We, however, fail to reproduce these results if the dynamo is in an advection-dominated regime.

The Waldmeier effect and the flux transport solar dynamo

We confirm that the evidence for the Waldmeier effect WE1 (the anticorrelation between rise times of sunspot cycles and their strengths) and the related effect WE2 (the correlation between rise rates of cycles and their strengths) is found in different kinds of sunspot data. We explore whether these effects can be explained theoretically on the basis of the flux transport dynamo models of sunspot cycles. Two sources of irregularities of sunspot cycles are included in our model: fluctuations in the poloidal field generation process and fluctuations in the meridional circulation. We find WE2 to be a robust result which is produced in different kinds of theoretical models for different sources of irregularities. The Waldmeier effect WE1, on the other hand, arises from fluctuations in the meridional circulation and is found only in the theoretical models with reasonably high turbulent diffusivity which ensures that the diffusion time is not more than a few years.

The solar dynamo with meridional circulation

We show that meridional circulation can have a profound influence on dynamo models for the solar cycle. Motivated by the observed tilt angles of sunspot groups we assume that the generation of the poloidal field takes place near the surface, while a shear layer of radial differential rotation produces the toroidal field at the bottom of the convection zone. Both layers are coupled by a circulation with a poleward directed flow in the upper part and an equatorward flow in the deep layers of the convection zone. The circulation forces the toroidal field belts (which are responsible for the surface activity) to move equatorward. This leads to butterfly diagrams in qualitative agreement with the observations, even if the dynamo wave would propagate poleward in the absence of circulation. This result opens the possibility to construct models for the solar cycle which are based on observational data (tilt angles, differential rotation, and meridional circulation).

A DYNAMO MODEL OF MAGNETIC ACTIVITY IN SOLAR-LIKE STARS WITH DIFFERENT ROTATIONAL VELOCITIES

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of 1M(circle dot) stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude f(m) of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of f(m). Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of the magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes rise radially from the bottom of the convection zone. Taking into account the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

Modélisation dynamo des cycles d’activité stellaire

Des décennies d’observation ont permis d’obtenir différentes relations liées à l’activité stellaire. Cependant, il est difficile de reproduire numériquement celles-ci à partir de modèles dynamo, puisqu’il n’y a pas de consensus sur le processus réellement présent dans les étoiles. Nous tentons de reproduire certaines de ces relations avec un modèle global 3D hydrodynamique qui nous fournit le profil de rotation différentielle et le tenseur-α utilisés en entrée dans un modèle de dynamo αΩ. Nous reproduisons ainsi efficacement la corrélation positive entre le rapport P_cyc⁄P_rot et P_rot^(-1). Par contre, nous échouons à reproduire les relations liant ω_cyc⁄Ω et l’énergie magnétique au nombre de Rossby. Cela laisse croire que la variation de P_cyc⁄P_rot avec la période de rotation est une caractéristique robuste du modèle αΩ, mais que l’effet-α ne serait pas le processus principal limitant l’amplitude du cycle. Cette saturation découlerait plutôt de la réaction magnétique sur l’écoulement à grande échelle.

A Coupled 2 × 2D Babcock–Leighton Solar Dynamo Model. I. Surface Magnetic Flux Evolution

The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.

Can Superflares Occur on Our Sun?

Recent observations of Sun-like stars, similar to our Sun in their surface temperature (5600-6000 K) and slow rotation (rotational period > 10 d), using the Kepler satellite by Maehara et al. (2012, Nature, 485, 478) have revealed the existence of superflares (with energy of 10(33)-10(35) erg). From statistical analyses of these superflares, it was found that superflares with energy of 10(34) erg occur once in 800 yr, and superflares with 10(35) erg occur once in 5000 yr. In this paper, we examine whether superflares with energy of 10(33)-10(35) erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas related to the mechanisms of the solar dynamo. If magnetic flux is generated by differential rotation at the base of the convection zone, as assumed in typical dynamo models, it is possible that the present Sun would generate a large sunspot with a total magnetic flux of similar to 2 x 10(23) Mx (= G cm(2)) within one solar cycle period, and lead to superflares with an energy of 10(34) erg. To store a total magnetic flux of similar to 10(24) Mx, necessary for generating 10(35) erg superflares, it would take similar to 40 yr. Hot Jupiters have often been argued to be a necessary ingredient for the generation of superflares, but we found that they do not play any essential role in the generation of magnetic flux in the star itself, if we consider only the magnetic interaction between the star and the hot Jupiter. This seems to be consistent with Maehara et al.'s finding of 148 superflare-generating solar-type stars that do not have a hot Jupiter-like companion. Altogether, our simple calculations, combined with Maehara et al.'s analysis of superflares on Sun-like stars, show that there is a possibility that superflares of 10(34) erg would occur once in 800 yr on our present Sun.

On the control of rapidly rotating convection by an axially varying magnetic field

The magnetic field in rapidly rotating dynamos is spatially inhomogeneous. The axial variation of the magnetic field is of particular importance because tall columnar vortices aligned with the rotation axis form at the onset of convection. The classical picture of magnetoconvection with constant or axially varying magnetic fields is that the Rayleigh number and wavenumber at onset decrease appreciably from their non-magnetic values. Nonlinear dynamo simulations show that the axial lengthscale of the self-generated azimuthal magnetic field becomes progressively smaller as we move towards a rapidly rotating regime. With a small-scale field, however, the magnetic control of convection is different from that in previous studies with a uniform or large-scale field. This study looks at the competing viscous and magnetic mode instabilities when the Ekman number E (ratio of viscous to Coriolis forces) is small. As the applied magnetic field strength (measured by the Elsasser number Lambda) increases, the critical Rayleigh number for onset of convection initially increases in a viscous branch, reaches an apex where both viscous and magnetic instabilities co-exist, and then falls in the magnetic branch. The magnetic mode of onset is notable for its dramatic suppression of convection in the bulk of the fluid layer where the field is weak. The viscous-magnetic mode transition occurs at Lambda similar to 1, which implies that small-scale convection can exist at field strengths higher than previously thought. In spherical shell dynamos with basal heating, convection near the tangent cylinder is likely to be in the magnetic mode. The wavenumber of convection is only slightly reduced by the self-generated magnetic field at Lambda similar to 1, in agreement with previous planetary dynamo models. The back reaction of the magnetic field on the flow is, however, visible in the difference in kinetic helicity between cyclonic and anticyclonic vortices.

Dynamos driven by weak thermal convection and heterogeneous outer boundary heat flux

We use numerical dynamo models with heterogeneous core-mantle boundary (CMB) heat flux to show that lower mantle lateral thermal variability may help support a dynamo under weak thermal convection. In our reference models with homogeneous CMB heat flux, convection is either marginally supercritical or absent, always below the threshold for dynamo onset. We find that lateral CMB heat flux variations organize the flow in the core into patterns that favour the growth of an early magnetic field. Heat flux patterns symmetric about the equator produce non-reversing magnetic fields, whereas anti-symmetric patterns produce polarity reversals. Our results may explain the existence of the geodynamo prior to inner core nucleation under a tight energy budget. Furthermore, in order to sustain a strong geomagnetic field, the lower mantle thermal distribution was likely dominantly symmetric about the equator. (C) 2015 Elsevier B.V. All rights reserved.

Imaging the cool stars in the interacting binaries AE Aqr, BV Cen and V426 Oph

It is well known that magnetic activity in late-type stars increases with increasing rotation rate. Using inversion techniques akin to medical imaging, the rotationally broadened profiles from such stars can be used to reconstruct `Doppler images' of the distribution of cool, dark starspots on their stellar surfaces. Interacting binaries, however, contain some of the most rapidly rotating late-type stars known and thus provide important tests of stellar dynamo models. Furthermore, magnetic activity is thought to play a key role in their evolution, behaviour and accretion dynamics. Despite this, we know comparatively little about the magnetic activity and its influence on such binaries. In this review we summarise the concepts behind indirect imaging of these systems, and present movies of the starspot distributions on the cool stars in some interacting binaries. We conclude with a look at the future opportunities that such studies may provide.