Strong mean field dynamos require supercritical helicity fluxes

Several one and two dimensional mean field models are analyzed where the effects of current helicity fluxes and boundaries are included within the framework of the dynamical quenching model. In contrast to the case with periodic boundary conditions, the final saturation energy of the mean field decreases inversely proportional to the magnetic Reynolds number. If a nondimensional scaling factor in the current helicity flux exceeds a certain critical value, the dynamo can operate even without kinetic helicity, i.e. it is based only on shear and current helicity fluxes, as first suggested by Vishniac & Cho (2001, ApJ 550, 752). Only above this threshold is the current helicity flux also able to alleviate catastrophic quenching. The fact that certain turbulence simulations have now shown apparently non‐resistively limited mean field saturation amplitudes may be suggestive of the current helicity flux having exceeded this critical value. Even below this critical value the field still reaches appreciable strength at the end of the kinematic phase, which is in qualitative agreement with dynamos in periodic domains. However, for large magnetic Reynolds numbers the field undergoes subsequent variations on a resistive time scale when, for long periods, the field can be extremely weak. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The dual role of shear in large‐scale dynamos

The role of shear in alleviating catastrophic quenching by shedding small‐scale magnetic helicity through fluxes along contours of constant shear is discussed. The level of quenching of the dynamo effect depends on the quenched value of the turbulent magnetic diffusivity. Earlier estimates that might have suffered from the force‐free degeneracy of Beltrami fields are now confirmed for shear flows where this degeneracy is lifted. For a dynamo that is saturated near equipartition field strength those estimates result in a 5‐fold decrease of the magnetic diffusivity as the magnetic Reynolds number based on the wavenumber of the energy‐carrying eddies is increased from 2 to 600. Finally, the role of shear in driving turbulence and large‐scale fields by the magneto‐rotational instability is emphasized. New simulations are presented and the 3π /4 phase shift between poloidal and toroidal fields is confirmed. It is suggested that this phase shift might be a useful diagnostic tool in identifying mean‐field dynamo action in simulations and to distinguish this from other scenarios invoking magnetic buoyancy as a means to explain migration away from the midplane. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The saturation of galactic dynamos

We summarize the current state of the long term discussion about the saturation mechanisms associated with rapid growth of small‐scale magnetic field, that operate in large‐scale galactic dynamos, and related problems with magnetic helicity conservation. Our general conclusion is that, taking into account magnetic helicity fluxes, large‐scale magnetic field can be amplified up to about the equipartition level. In contrast, models without helicity fluxes give an initial temporal magnetic field growth, but then decay. In our opinion, it is more appropriate to refer to the situation as a “potentially catastrophic scenario” rather than as “catastrophic α‐quenching” (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A magnetohydrodynamical model for the formation of episodic jets

Episodic ejection of plasma blobs has been observed in many black hole systems. While steady, continuous jets are believed to be associated with large-scale open magnetic fields, what causes the episodic ejection of blobs remains unclear. Here by analogy with the coronal mass ejection on the Sun, we propose a magnetohydrodynamical model for episodic ejections from black holes associated with the closed magnetic fields in an accretion flow. Shear and turbulence of the accretion flow deform the field and result in the formation of a flux rope in the disc corona. Energy and helicity are accumulated and stored until a threshold is reached. The system then loses its equilibrium and the flux rope is thrust outward by the magnetic compression force in a catastrophic way. Our calculations show that for parameters appropriate for the black hole in our Galactic centre, the plasmoid can attain relativistic speeds in about 35 min.     

Magnetic helicity in stellar dynamos: new numerical experiments

The theory of large scale dynamos is reviewed with particular emphasis on the magnetic helicity constraint in the presence of closed and open boundaries. In the presence of closed or periodic boundaries, helical dynamos respond to the helicity constraint by developing small scale separation in the kinematic regime, and by showing long time scales in the nonlinear regime where the scale separation has grown to the maximum possible value. A resistively limited evolution towards saturation is also found at intermediate scales before the largest scale of the system is reached. Larger aspect ratios can give rise to different structures of the mean field which are obtained at early times, but the final saturation field strength is still decreasing with decreasing resistivity. In the presence of shear, cyclic magnetic fields are found whose period is increasing with decreasing resistivity, but the saturation energy of the mean field is in strong super‐equipartition with the turbulent energy. It is shown that artificially induced losses of small scale field of opposite sign of magnetic helicity as the large scale field can, at least in principle, accelerate the production of large scale (poloidal) field. Based on mean field models with an outer potential field boundary condition in spherical geometry, we verify that the sign of the magnetic helicity flux from the large scale field agrees with the sign of α. For solar parameters, typical magnetic helicity fluxes lie around 10⁴⁷ Mx² per cycle.     

Magnetic Energy and Helicity Fluxes at the Photospheric Level

The source of coronal magnetic energy and helicity lies below the surface of the Sun, probably in the convective zone dynamo. Measurements of magnetic and velocity fields can capture the fluxes of both magnetic energy and helicity crossing the photosphere. We point out the ambiguities which can occur when observations are used to compute these fluxes. In particular, we show that these fluxes should be computed only from the horizontal motions deduced by tracking the photospheric cut of magnetic flux tubes. These horizontal motions include the effect of both the emergence and the shearing motions whatever the magnetic configuration complexity is. We finally analyze the observational difficulties involved in deriving such fluxes, in particular the limitations of the correlation tracking methods.     

Alleviation of catastrophic quenching in solar dynamo model with nonlocal alpha‐effect

The nonlocal alpha‐effect of Babcock‐Leighton type is not prone to the catastrophic quenching due to conservation of magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha‐effect, the diamagnetic pumping, and dynamical equation for the magnetic alpha‐effect. The same model shows catastrophic quenching when the alpha‐effect is changed to its local formulation. The nonlocal model shows a preferred excitation of magnetic fields of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal‐to‐polar fields ratio of about a thousand (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Catastrophe of coronal magnetic flux ropes in partially open magnetic fields

Using a 2.5-D, time-dependent ideal MHD model in Cartesian coordinates, a numerical study is carried out to find equilibrium solutions associated with a magnetic flux rope in the corona. The ambient magnetic field is partially open, consisting of a closed arcade in the center and an open field at the flank. The coronal magnetic flux rope is characterized by its magnetic properties, including the axial and annular magnetic fluxes and the magnetic helicity, and its geometrical features, including the height of the rope axis, the halfwidth of the rope and the length of the vertical current sheet below the rope. It is shown that for a given partially open ambient magnetic field, the dependence of the geometrical features on the magnetic properties displays a catastrophic behavior, namely, there exists a certain critical point, across which an infinitesimal enhancement of the magnetic parameters causes a finite jump of the geometrical parameters for the rope. The amplitude of the jump depends on the extent to which the ambient magnetic field in open, and approaches to zero when the ambient magnetic field becomes completely closed. The implication of such a catastrophe in solar active phenomena is briefly discussed.     

Basic Properties of Mutual Magnetic Helicity

We derive the magnetic helicity for configurations formed by flux tubes contained fully or only partially in the spatial domain considered (called closed and open configurations, respectively). In both cases, magnetic helicity is computed as the sum of mutual helicity over all possible pairs of magnetic flux tubes weighted by their magnetic fluxes. We emphasize that these mutual helicities have properties which are not those of mutual inductances in classical circuit theory. For closed configurations, the mutual helicity of two closed flux tubes is their relative winding around each other (known as the Gauss linkage number). For open configurations, the magnetic helicity is derived directly from the geometry of the interlaced flux tubes so it can be computed without reference to a ground state (such as a potential field). We derive the explicit expression in the case of a planar and spherical boundary. The magnetic helicity has two parts. The first one is given only by the relative positions of the flux tubes on the boundary. It is the only part if all flux tubes are arch-shaped. The second part counts the integer number of turns each pair of flux tubes wind about each other. This provides a general method to compute the magnetic helicity with discrete or continuous distributions of magnetic field. The method sets closed and open configurations on an equal level within the same theoretical framework.     

Predictions of Energy and Helicity in Four Major Eruptive Solar Flares

In order to better understand the solar genesis of interplanetary magnetic clouds (MCs), we model the magnetic and topological properties of four large eruptive solar flares and relate them to observations. We use the three-dimensional Minimum Current Corona model (Longcope, 1996, Solar Phys. 169, 91) and observations of pre-flare photospheric magnetic field and flare ribbons to derive values of reconnected magnetic flux, flare energy, flux rope helicity, and orientation of the flux-rope poloidal field. We compare model predictions of those quantities to flare and MC observations, and within the estimated uncertainties of the methods used find the following: The predicted model reconnection fluxes are equal to or lower than the reconnection fluxes inferred from the observed ribbon motions. Both observed and model reconnection fluxes match the MC poloidal fluxes. The predicted flux-rope helicities match the MC helicities. The predicted free energies lie between the observed energies and the estimated total flare luminosities. The direction of the leading edge of the MC’s poloidal field is aligned with the poloidal field of the flux rope in the AR rather than the global dipole field. These findings compel us to believe that magnetic clouds associated with these four solar flares are formed by low-corona magnetic reconnection during the eruption, rather than eruption of pre-existing structures in the corona or formation in the upper corona with participation of the global magnetic field. We also note that since all four flares occurred in active regions without significant pre-flare flux emergence and cancelation, the energy and helicity that we find are stored by shearing and rotating motions, which are sufficient to account for the observed radiative flare energy and MC helicity.     

Solar dynamo model with nonlocal alpha-effect

The first results of the solar dynamo model that allows for the diamagnetic effect of inhomogeneous turbulence and the nonlocal alpha-effect due to the rise ofmagnetic loops are discussed. The nonlocal alpha-effect is not subject to the catastrophic quenching related to the conservation of magnetic helicity. Given the diamagnetic pumping, the magnetic fields are concentrated near the base of the convection zone, although the distributed-typemodel covers the entire thickness of the convection zone. Themagnetic cycle period, the equatorial symmetry of the field, its meridional drift, and the polar-to-toroidal field ratio obtained in the model are in agreement with observations. There is also some disagreement with observations pointing the ways of improving the model.     

Magnetic flux separation in photospheric convection

Three-dimensional non-linear magnetoconvection in a strongly stratified compressible layer exhibits different patterns as the strength of the imposed magnetic field is reduced. There is a transition from a magnetically dominated regime, with small-scale convection in slender hexagonal cells, to a convectively dominated regime, with clusters of broad rising plumes that confine the magnetic flux to narrow lanes where fields are locally intense. Both patterns can coexist for intermediate field strengths, giving rise to flux separation: clumps of vigorously convecting plumes, from which magnetic flux has been excluded, are segregated from regions with strong fields and small-scale convection. A systematic numerical investigation of these different states shows that flux separation can occur over a significant parameter range and that there is also hysteresis. The results are related to the fine structure of magnetic fields in sunspots and in the quiet Sun.     

Galactic dynamo with helicity fluxes

We construct an approximation for the magnetic field of galaxies that takes into account the magnetic helicity conservation law. In our calculations, we use the fact that the galactic disk is fairly thin and, therefore, the magnetic field component perpendicular to the galactic disk can be neglected (the so-called no-z approximation). However, an averaging of the magnetic field over the entire galaxy, as was done in previous works, is not performed. Our results are compared both with the approximation that disregards the helicity flux and with the results obtained in models with helicity fluxes but without averaging. We show that, compared to the classical model, there are a number of new effects (for example, magnetic field oscillations) and, compared to the model with averaging, the behavior of the magnetic field “softens”: its stationary value is reached more slowly and the oscillation amplitude decreases. This is because the dissipative processes changing the magnetic field growth rate are taken into account in our model. In contrast to the model with averaging, here it becomes possible to construct the dependence of the magnetic field and helicity on the distance from the galactic center.     

The Formation of a Small-Scale Filament After Flux Emergence on the Quiet Sun

We present observations of the formation process of a small-scale filament on the quiet Sun during 5?–?6 February 2016 and investigate its formation cause. Initially, a small dipole emerged, and its associated arch filament system was found to reconnect with overlying coronal fields accompanied by numerous extreme ultraviolet bright points. When the bright points faded, many elongated dark threads formed and bridged the positive magnetic element of the dipole and the external negative network fields. Interestingly, an anticlockwise photospheric rotational motion (PRM) set in within the positive endpoint region of the newborn dark threads following the flux emergence and lasted for more than 10 hours. Under the drive of the PRM, these dispersive dark threads gradually aligned along the north-south direction and finally coalesced into an inverse S-shaped filament. Consistent with the dextral chirality of the filament, magnetic helicity calculations show that an amount of negative helicity was persistently injected from the rotational positive magnetic element and accumulated during the formation of the filament. These observations suggest that twisted emerging fields may lead to the formation of the filament via reconnection with pre-existing fields and release of its inner magnetic twist. The persistent PRM might trace a covert twist relaxation from below the photosphere to the low corona.     

The Magnetic Helicity Injected by Shearing Motions

Photospheric shearing motions are one of the possible ways to inject magnetic helicity into the corona. We explore their efficiency as a function of their particular properties and those of the magnetic field configuration. Based on the work of M. A. Berger, we separate the helicity injection into two terms: twist and writhe. For shearing motions concentrated between the centers of two magnetic polarities the helicity injected by twist and writhe add up, while for spatially more extended shearing motions, such as differential rotation, twist and writhe helicity have opposite signs and partially cancel. This implies that the amount of injected helicity can change in sign with time even if the shear velocity is time independent. We confirm the amount of helicity injected by differential rotation in a bipole in the two particular cases studied by DeVore (2000), and further explore the parameter space on which this injection depends. For a given latitude, tilt and magnetic flux, the generation of helicity is slightly more efficient in young active regions than in decayed ones (up to a factor 2). The helicity injection is mostly affected by the tilt of the AR with respect to the solar equator. The total helicity injected by shearing motions, with both spatial and temporal coherence, is at most equivalent to that of a twisted flux tube having the same magnetic flux and a number of turns of 0.3. In the solar case, where the motions have not such global coherence, the injection of helicity is expected to be much smaller, while for differential rotation this maximum value reduces to 0.2 turns. We conclude that shearing motions are a relatively inefficient way to bring magnetic helicity into the corona (compared to the helicity carried by a significantly twisted flux tube).     

Helicity as a Component of Filament Formation

In this paper we seek the origin of the axial component of the magnetic field in filaments by adapting theory to observations. A previous paper (Mackay, Gaizauskas, and van Ballegooijen, 2000) showed that surface flows acting on potential magnetic fields for 27 days – the maximum time between the emergence of magnetic flux and the formation of large filaments between the resulting activity complexes – cannot explain the chirality or inverse polarity nature of the observed filaments. We show that the inclusion of initial helicity, for which there is observational evidence, in the flux transport model results in sufficiently strong dextral fields of inverse polarity to account for the existence and length of an observed filament within the allotted time. The simulations even produce a large length of dextral chirality when just small amounts of helicity are included in the initial configuration. The modeling suggests that the axial field component in filaments can result from a combination of surface (flux transport) and sub-surface (helicity) effects acting together. Here surface effects convert the large-scale helicity emerging in active regions into a smaller-scale magnetic-field component parallel to the polarity inversion line so as to form a magnetic configuration suitable for a filament.     

The Magnetic Field, Horizontal Motion and Helicity in a Fast Emerging Flux Region which Eventually forms a Delta Spot

We studied the behavior of magnetic field, horizontal motion and helicity in a fast emerging flux region NOAA 10488 which eventually forms a δ spot. It is found that the rotation of photospheric footpoints forms in the earlier stage of magnetic flux emergence and the relative shear motion of different magnetic flux systems appears later in this active region (AR). Therefore the emerging process of the AR can be separated into two phases: rotation and shear. We have computed the magnetic helicity injected into the corona using the local correlation tracking (LCT) technique. Furthermore we determined the vertical component of current helicity density and the vertical component of induction electric fields E_z = (V× B)_z in the photosphere. Particularly we have presented the comparison of the injection rate of magnetic helicity and the variation of the current helicity density. The main results are as follows: (1) The strong shear motion (SSM) between the new emerging flux system and the old one brings more magnetic helicity into the corona than the twisting motions. (2) After the maturity of the main bipolar spots, their twist decreases and the SSM becomes dominant and the major contributor of magnetic non-potentiality in the solar atmosphere in this AR. (3) The positions of the maxima of E_z (about 0.1 ∼ 0.2 V cm⁻¹) shift from the twisting areas to the areas showing SSMs as the AR evolved from the rotation phase to the shear one, but no obvious correlation is found between the kernels of Hα flare and E_z for the M1.6 flare in this AR. (4) The coronal helicity inferred from the horizontal motion of this AR amounts to −6 × 10⁴³ Mx². It is comparable with the coronal helicity of ARs producing flares with coronal mass ejections (CMEs) or helicity carried away by magnetic clouds (MCs) reported in previous studies (Nindos, Zhang, and Zhang, 2003; Nindos and Andrews, 2004). In addition, the formation of the δ configuration in this AR belongs to the third formation type indicated by Zirin and Liggett (1987), i.e., collision of opposite polarities from different dipoles, and can be naturally explained by the SSM.     

The Roles of Reconnected Flux and Overlying Fields in CME Speeds

Researchers have reported i) correlations of coronal mass ejection (CME) speeds and the total photospheric magnetic flux swept out by flare ribbons in flare-associated eruptive events, and, separately, ii) correlations of CME speeds and more rapid decay, with height, of magnetic fields in potential-field coronal models above eruption sites. Here, we compare the roles of both ribbon fluxes and the decay rates of overlying fields in a set of 16 eruptive events. We confirm previous results that higher CME speeds are associated with both higher ribbon fluxes and more rapidly decaying overlying fields. We find the association with ribbon fluxes to be weaker than a previous report, but stronger than the dependence on the decay rate of overlying fields. Since the photospheric ribbon flux is thought to approximate the amount of coronal magnetic flux reconnected during the event, the correlation of speeds with ribbon fluxes suggests that reconnection plays some role in accelerating CMEs. One possibility is that reconnected fields that wrap around the rising ejection produce an increased outward hoop force, thereby increasing CME acceleration. The correlation of CME speeds with more rapidly decaying overlying fields might be caused by greater downward magnetic tension in stronger overlying fields, which could act as a source of drag on rising ejections.     

EVOLUTION OF MAGNETIC HELICITY IN MAGNETIC RECONNECTION

This paper presents a definition of magnetic helicity specifically for two-dimensional magnetic fields and derives the associated helicity equation. The newly defined helicity is closely related to its three-dimensional counterpart and serves as a measure of the shear of magnetic field. Based on this, a numerical simulation is carried out on magnetic reconnection occurring in the lower solar atmosphere. It is found that the helicity dissipation due to magnetic reconnection is very small. A large amount of helicity is transferred upward and escapes from the domain of the solution, and the total helicity is approximately conserved during the magnetic reconnection and helicity transfer. This is in support of the applicability of a postulate, which was proposed by Taylor (1974, 1986) concerning the approximate conservation of magnetic helicity in the presence of resistive dissipation and magnetic reconnection in a highly conductive laboratory plasma, to the solar atmosphere.