A coupled model of magnetic flux generation and transport in stars

We present a combined model for magnetic field generation and transport in cool stars with outer convection zones. The mean toroidal magnetic field, which is generated by a cyclic thin-layer α Ω dynamo at the bottom of the convection zone is taken to determine the emergence probability of magnetic flux tubes in the photosphere. Following the nonlinear rise of the unstable thin flux tubes, emergence latitudes and tilt angles of bipolar magnetic regions are determined. These quantities are put into a surface flux transport model, which simulates the surface evolution of magnetic flux under the effects of large-scale flows and turbulent diffusion. First results are discussed for the case of the Sun and for more rapidly rotating solar-type stars. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stability of magnetic flux tubes in close binary stars

Photometric and Doppler imaging observations of active binaries indicate the existence of starspots at preferred longitudes (position angles with respect to the companion star). We investigate the stability of magnetic flux tubes in the convection zone of close, fast‐rotating binary stars and explore whether the observed preferred longitudes could be caused by tidal forces and the deformation of the active star. We assume a synchronized binary system with spin axes perpendicular to the orbital plane and a rotation period of a few days. The tidal force and the deviation from spherical structure are considered in lowest‐order perturbation theory. The magnetic field is in the form of toroidal magnetic flux rings, which are stored in mechanical equilibrium within the stably stratified overshoot region beneath the convection zone until the field has grown sufficiently strong for the undulatory instability to initiate the formation of rising loops. Frequencies and geometry of stable as well as growth rates of unstable eigenmodes are determined by linear stability analysis. Particular consideration is given to the question whether the effects of tidal forces and perturbations of the stellar structure can force a rising flux loop to enter the convection zone at specific longitudes.     

Comparing Simulations of Rising Flux Tubes Through the Solar Convection Zone with Observations of Solar Active Regions: Constraining the Dynamo Field Strength

We study how active-region-scale flux tubes rise buoyantly from the base of the convection zone to near the solar surface by embedding a thin flux tube model in a rotating spherical shell of solar-like turbulent convection. These toroidal flux tubes that we simulate range in magnetic field strength from 15 kG to 100 kG at initial latitudes of 1^° to 40^° in both hemispheres. This article expands upon Weber, Fan, and Miesch (Astrophys. J. 741, 11, 2011) (Article 1) with the inclusion of tubes with magnetic flux of 10²⁰ Mx and 10²¹ Mx, and more simulations of the previously investigated case of 10²² Mx, sampling more convective flows than the previous article, greatly improving statistics. Observed properties of active regions are compared to properties of the simulated emerging flux tubes, including: the tilt of active regions in accordance with Joy’s Law as in Article 1, and in addition the scatter of tilt angles about the Joy’s Law trend, the most commonly occurring tilt angle, the rotation rate of the emerging loops with respect to the surrounding plasma, and the nature of the magnetic field at the flux tube apex. We discuss how these diagnostic properties constrain the initial field strength of the active-region flux tubes at the bottom of the solar convection zone, and suggest that flux tubes of initial magnetic field strengths of ≥?40 kG are good candidates for the progenitors of large (10²¹ Mx to 10²² Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 10²² Mx. With the addition of more magnetic flux values and more simulations, we find that for all magnetic field strengths, the emerging tubes show a positive Joy’s Law trend, and that this trend does not show a statistically significant dependence on the magnetic flux.     

Buoyant Magnetic Loops Generated by Global Convective Dynamo Action

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

The evolution of loop structures in flux rings within the solar convection zone

Choudhuri and Gilman (1987) considered certain implications of the hypothesis that the magnetic flux within the Sun is generated at the bottom of the convection zone and then rises through it. Taking flux rings symmetric around the rotation axis and using reasonable values of different parameters, they found that the Coriolis force deflects these flux rings into trajectories parallel to the rotation axis so that they emerge at rather high latitudes. This paper looks into the question of whether the action of the Coriolis force is subdued when the initial configuration of the flux ring has non-axisymmetries in the form of loop structures. The results depend dramatically on whether the flux ring with the loops lies completely within the convection zone or whether the lower parts of it are embedded in the stable layers underneath the convection zone. In the first case, the Coriolis force supresses the non-axisymmetric perturbations so that the flux ring tends to remain symmetric and the trajectories are very similar to those of Choudhuri and Gilman (1987). In the second case, however, the lower parts of the flux ring may remain anchored underneath the bottom of the convection zone, but the upper parts of the loops still tend to move parallel to the rotation axis and emerge at high latitudes. Thus the problem of the magnetic flux not being able to come out at the sunspot latitudes still persists after the non-axisymmetries in the flux rings are taken into account.National Center for Atmospheric Research is sponsored by the National Science Foundation.     

On the Connection Between Mean Field Dynamo Theory and Flux Tubes

Mean field dynamo theory deals with various mean quantities and does not directly throw any light on the question of existence of flux tubes. We can, however, draw important conclusions about flux tubes in the interior of the Sun by combining additional arguments with the insights gained from solar dynamo solutions. The polar magnetic field of the Sun is of order 10 G, whereas the toroidal magnetic field at the bottom of the convection zone has been estimated to be 100000 G. Simple order-of-magnitude estimates show that the shear in the tachocline is not sufficient to stretch a 10 G mean radial field into a 100000 G mean toroidal field. We argue that the polar field of the Sun must get concentrated into intermittent flux tubes before it is advected to the tachocline. We estimate the strengths and filling factors of these flux tubes. Stretching by shear in the tachocline is then expected to produce a highly intermittent magnetic configuration at the bottom of the convection zone. The meridional flow at the bottom of the convection zone should be able to carry this intermittent magnetic field equatorward, as suggested recently by Nandy and Choudhuri (2002). When a flux tube from the bottom of the convection zone rises to a region of pre-existing poloidal field at the surface, we point out that it picks up a twist in accordance with the observations of current helicities at the solar surface.     

On the Dynamic Disconnection of Rising Omega Loops

The tilt-angle variation of solar bipolar magnetic regions (BMRs) in space and time is a probable indicator of large-scale properties of magnetic field distribution, the dynamo and other different processes affecting the rising of magnetic flux tubes. We analyze the tilt-angle distribution and evolution of clearly oriented BMRs using sunspot position and area measurements from the Greenwich Photo-Heliographic Results (GPHR) supplemented by high-quality full-disk sunspot drawings made in the Haynald Observatory in Hungary. Our work is made on the basis of the 14th sunspot cycle and we compare our results with present theories of magnetic flux evolution in the convective zone. Our results support the assumption that the scatter of tilt angles of BMRs around Joy`s law is determined by the convective turbulence. Furthermore, we find that regular BMRs grouped by age generally do not show the phenomenon of toroidal relaxation towards the east–west direction. We interpret this by the disconnection of -loops from the bottom of the convection zone.     

Limits on the magnetic field strength at the base of the solar convection zone

Howard (1993) finds a relationship between the tilt angles of BMRs (Bipolar Magnetic Regions) and the separation between their leading and following polarities; the tilt angle increases with polarity separation. Here we present a more detailed analysis of this relationship and show that this effect constrains the strength of the magnetic field at the bottom of the convection zone to a value between 40 and 150 kG, which confirms the constraints put by D'Silva and Choudhuri (1993) based on Joy's law (the tilt-latitude relationship), through an entirely different approach.Operated by the Association of Universities for Research in Astronomy, Inc., under CooperativeOperated by the Association of Universities for Research in Astronomy, Inc., under Cooperative     

On the appearance of magnetic flux in the solar photosphere

Ideas and models for the appearance of photospheric magnetic structure are confronted with observational data. Some findings are: The magnetic flux emerging in an active region consists of a bundle of flux tubes which were already concentrated before penetrating into the photosphere. A model of a rising bunch of flux tubes joining into a few strands at larger depths describes the coalescence of spots near the leading and following edges of the active region while more flux may surface near the center of the region. There is no observational evidence for appreciable helical twists in the flux bundles.Throughout the region's lifetime the magnetic elements move coherently, the whole magnetic structure rotates faster than the quiet photosphere. In active regions the convective flow at scales larger than the granulation is arrested by the magnetic structure. The long-lived supergranular cells around spots and in the enhanced network in turn determine the decay properties of spots and facular clusters. The modulation of the convective flow by the magnetic structure explains the slow dispersal of faculae.The hierarchy of magnetic elements (sunspots-pores-knots-facular clusters-facular points) may be explained by a set of magnetostatic flux tube models in the top of the convection zone. The underlying assumptions are that the heat flow along the magnetic field is reduced and that there is no heat exchange across the field except by radiation.A tentative model is proposed to account for the amplification, ascent and emergence of intense flux bundles. The assumptions are: (i) the field is concentrated in toroidal bundles by differential rotation, (ii) in the deep convection zone flux bundles are contained by the external turbulent pressure, and (iii) for field strengths up to the equipartition value efficient lateral heat exchange is possible. After a loop has surfaced radiative cooling and subsequent convective downflow reduce the temperature in the top of the flux tubes which then contract to field strengths well above the local equipartition value. There the heat flow is channelled along the field, which creates the conditions for the magnetostatic flux tube models without requiring a blocking of the heat flow somewhere within the tubes.The paper contains a brief review on the evolution of the magnetic field from the emergence in active regions up to the enigmatic disappearance, and a list of topics for further observational investigation.     

Dynamo saturation in rapidly rotating solar-type stars

The magnetic activity of solar-type stars generally increases with stellar rotation rate. The increase, however, saturates for fast rotation. The Babcock-Leighton mechanism of stellar dynamos saturates as well when the mean tilt angle of active regions approaches ninety degrees. Saturation of magnetic activity may be a consequence of this property of the Babcock-Leighton mechanism. Stellar dynamo models with a tilt angle proportional to the rotation rate are constructed to probe this idea.Two versions of the model- treating the tilt angles globally and using Joy's law for its latitude dependence- are considered. Both models show a saturation of dynamogenerated magnetic flux at high rotation rates. The model with latitude-dependent tilt angles also shows a change in dynamo regime in the saturation region. The new regime combines a cyclic dynamo at low latitudes with an(almost) steady polar dynamo.     

A dynamo scenario

In this paper, the term dynamo refers to the complex of physical mechanisms that cause solar magnetic activity in all its manifestations. Properties of that dynamo are inferred from observational indications to fit them into a scenario. Properties and models of the manifestations of strong magnetic field are briefly summarized, together with their formation during the emergence of -shaped loops from the toroidal flux system in the interface below the convection zone. The evolution of magnetic concentrations and the flux removal from the atmosphere, with indications for flux retraction, are considered. Then the weak (INF) fields are discussed, together with the role of upward floating LI- shaped loops in the removal of toroidal flux. Finally features of strong and weak fields are fitted into a scenario for a cyclic dynamo, in which the regeneration of the poloidal field of proper sign relies on the cancellation of magnetic flux in the surface transport interpreted as reconnection, followed by retraction of reconnected loops.Dedicated to Cornelis de JagerBased on an invited talk during Solar Cycle Workshop, March 28–30, 1996, Tucson.     

Dynamical behavior of strong magnetic fields in the solar convection zone

Magnetic buoyancy is thought to play an important role in the dynamical behavior of the Sun's magnetic field in the convection zone. Magnetic buoyancy is commonly thought to cause inescapable rapid loss of toroidal flux from much of the convection zone, thereby suppressing effective operation of a solar dynamo. This paper re-examines the detailed character of magnetic buoyancy, especially as it is influenced by the magnetic field's effect on heat transport and temperature gradients in the convection zone. It is suggested that suppression of convective heat transport across strong magnetic flux tubes can alter the temperature within the tubes and can subdue, or even reverse, the effect of magnetic buoyancy.     

3-D non-linear evolution of a magnetic flux tube in a spherical shell: The isentropic case

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

Helicity of solar active regions from a dynamo model

We calculate helicities of solar active regions based on the idea that poloidal flux lines get wrapped around a toroidal flux tube rising through the convection zone, thereby giving rise to the helicity. We use our solar dynamo model based on the Babcock-Leighton α-effect to study how helicity varies with latitude and time.     

Double maxima of 11-year solar cycles

We propose a scenario to explain the observed phenomenon of double maxima of sunspot cycles, including the generation of a magnetic field near the bottom of the solar convection zone (SCZ) and the subsequent rise of the field from the deep layers to the surface in the royal zone. Five processes are involved in the restructuring of the magnetic field: the Ω-effect, magnetic buoyancy, macroscopic turbulent diamagnetism, rotary ?ρ-effect, and meridional circulation. It is found that the restructuring of magnetism develops differently in high-latitude and equatorial domains of the SCZ. A key role in the proposed mechanism of the double maxima is played by two waves of toroidal fields from the lower base of the SCZ to the solar surface in the equatorial domain. The deep toroidal fields are excited by the Ω-effect near the tachocline at the beginning of the cycle. Then these fields are transported to the surface due to the combined effect of magnetic buoyancy, macroscopic turbulent diamagnetism, and the rotary magnetic ?ρ-flux in the equatorial domain. After a while, these magnetic fragments can be observed as bipolar sunspot groups at the middle latitudes in the royal zone. This first, upward-directed wave of toroidal fields produces the main maximum of sunspot activity. However, the underlying toroidal fields in the high-latitude polar domains are blocked at the beginning of the cycle near the SCZ bottom by two antibuoyancy effects — the downward turbulent diamagnetic transfer and the magnetic ?ρ-pumping. In approximately 1 or 2 years, a deep equatorward meridional flow transfers these fields to low-latitude parts of the equatorial domain (where there are favorable conditions for magnetic buoyancy), and the belated magnetic fields (the second wave of toroidal fields) rise to the surface. When this second batch of toroidal fields comes to the solar surface at low latitudes, it leads to the second sunspot maximum.     

How to Reach Superequipartition Field Strengths in Solar Magnetic Flux Tubes

A number of independent arguments indicate that the toroidal flux system responsible for the sunspot cycle is stored at the base of the convection zone in the form of flux tubes with field strength close to 10⁵ G. Although the evidence for such strong fields is quite compelling, how such field strength can be reached is still a topic of debate. Flux expulsion by convection should lead to about the equipartition field strength, but the magnetic energy density of a 10⁵-G field is two orders of magnitude larger than the mean kinetic energy density of convective motions. Line stretching by differential rotation (i.e., the “Ω effect” in the classical mean-field dynamo approach) probably plays an important role, but arguments based on energy considerations show that it does not seem feasible that a 10⁵-G field can be produced in this way. An alternative scenario for the intensification of the toroidal flux system in the overshoot layer is related to the explosion of rising, buoyantly unstable magnetic flux tubes, which opens a complementary mechanism for magnetic-field intensification. A parallelism is pointed out with the mechanism of “convective collapse” for the intensification of photospheric magnetic flux tubes up to field strengths well above equipartition; both mechanisms, which are fundamentally thermal processes, are reviewed.     

Differential rotation and meridional circulation in global models of solar convection

In the outer envelope of the Sun and in other stars, differential rotation and meridional circulation are maintained via the redistribution of momentum and energy by convective motions. In order to properly capture such processes in a numerical model, the correct spherical geometry is essential. In this paper I review recent insights into the maintenance of mean flows in the solar interior obtained from high-resolution simulations of solar convection in rotating spherical shells. The Coriolis force induces a Reynolds stress which transports angular momentum equatorward and also yields latitudinal variations in the convective heat flux. Meridional circulations induced by baroclinicity and rotational shear further redistribute angular momentum and alter the mean stratification. This gives rise to a complex nonlinear interplay between turbulent convection, differential rotation, meridional circulation, and the mean specific entropy profile. I will describe how this drama plays out in our simulations as well as in solar and stellar convection zones. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mechanism of cyclically polarity reversing solar magnetic cycle as a cosmic dynamo

We briefly describe historical development of the concept of solar dynamo mechanism that generates electric current and magnetic field by plasma flows inside the solar convection zone. The dynamo is the driver of the cyclically polarity reversing solar magnetic cycle. The reversal process can easily and visually be understood in terms of magnetic field line stretching and twisting and folding in three-dimensional space by plasma flows of differential rotation and global convection under influence of Coriolis force. This process gives rise to formation of a series of huge magnetic flux tubes that propagate along iso-rotation surfaces inside the convection zone. Each of these flux tubes produces one solar cycle. We discuss general characteristics of any plasma flows that can generate magnetic field and reverse the polarity of the magnetic field in a rotating body in the Universe. We also mention a list of problems which are currently being disputed concerning the solar dynamo mechanism together with observational evidences that are to be constraints as well as verifications of any solar cycle dynamo theories of short and long term behaviors of the Sun, particularly time variations of its magnetic field, plasma flows, and luminosity.