Dynamo coefficients from local simulations of the turbulent ISM

Observations in polarized emission reveal the existence of large‐scale coherent magnetic fields in a wide range of spiral galaxies. Radio‐polarization data show that these fields are strongly inclined towards the radial direction, with pitch angles up to 35° and thus cannot be explained by differential rotation alone. Global dynamo models describe the generation of the radial magnetic field from the underlying turbulence via the so called α ‐effect. However, these global models still rely on crude assumptions about the small‐scale turbulence. To overcome these restrictions we perform fully dynamical MHD simulations of interstellar turbulence driven by supernova explosions. From our simulations we extract profiles of the contributing diagonal elements of the dynamo α ‐tensor as functions of galactic height. We also measure the coefficients describing vertical pumping and find that the ratio between these two effects has been overestimated in earlier analytical work, where dynamo action seemed impossible. In contradiction to these models based on isolated remnants we always find the pumping to be directed inward. In addition we observe that depends on whether clustering in terms of superbubbles is taken into account. Finally, we apply a test field method to derive a quantitative measure of the turbulent magnetic diffusivity which we determine to be ∼2 kpckms^–1. (© 2008 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)     

Dynamical quenching with non‐local α and downward pumping

In light of new results, the one‐dimensional mean‐field dynamo model of Brandenburg & Käpylä (2007) with dynamical quenching and a nonlocal Babcock‐Leighton α effect is re‐examined for the solar dynamo. We extend the one‐dimensional model to include the effects of turbulent downward pumping (Kitchatinov & Olemskoy 2011), and to combine dynamical quenching with shear. We use both the conventional dynamical quenching model of Kleeorin & Ruzmaikin (1982) and the alternate one of Hubbard & Brandenburg (2011), and confirm that with varying levels of non‐locality in the α effect, and possibly shear as well, the saturation field strength can be independent of the magnetic Reynolds number. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Large‐scale magnetic flux concentrations from turbulent stresses

In this study we provide the first numerical demonstration of the effects of turbulence on the mean Lorentz force and the resulting formation of large‐scale magnetic structures. Using three‐dimensional direct numerical simulations (DNS) of forced turbulence we show that an imposed mean magnetic field leads to a decrease of the turbulent hydromagnetic pressure and tension. This phenomenon is quantified by determining the relevant functions that relate the sum of the turbulent Reynolds and Maxwell stresses with the Maxwell stress of the mean magnetic field. Using such a parameterization, we show by means of two‐dimensional and three‐dimensional mean‐field numerical modelling that an isentropic density stratified layer becomes unstable in the presence of a uniform imposed magnetic field. This large‐scale instability results in the formation of loop‐like magnetic structures which are concentrated at the top of the stratified layer. In three dimensions these structures resemble the appearance of bipolar magnetic regions in the Sun. The results of DNS and mean‐field numerical modelling are in good agreement with theoretical predictions. We discuss our model in the context of a distributed solar dynamo where active regions and sunspots might be rather shallow phenomena (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling the total and polarized emission in evolving galaxies: “Spotty” magnetic structures

Future radio observations with the Square Kilometre Array (SKA) and its precursors will be sensitive to trace spiral galaxies and their magnetic field configurations up to redshift z ≈ 3. We suggest an evolutionary model for the magnetic configuration in star‐forming disk galaxies and simulate the magnetic field distribution, the total and polarized synchrotron emission, and the Faraday rotation measures for disk galaxies at z ≲ 3. Since details of dynamo action in young galaxies are quite uncertain, we model the dynamo action heuristically relying only on well‐established ideas of the form and evolution of magnetic fields produced by the mean‐field dynamo in a thin disk. We assume a small‐scale seed field which is then amplified by the small‐scale turbulent dynamo up to energy equipartition with kinetic energy of turbulence. The large‐scale galactic dynamo starts from seed fields of 100 pc and an averaged regular field strength of 0.02 μG, which then evolves to a “spotty” magnetic field configuration in about 0.8 Gyr with scales of about one kpc and an averaged regular field strength of 0.6 μG. The evolution of these magnetic spots is simulated under the influence of star formation, dynamo action, stretching by differential rotation of the disk, and turbulent diffusion. The evolution of the regular magnetic field in a disk of a spiral galaxy, as well as the expected total intensity, linear polarization and Faraday rotation are simulated in the rest frame of a galaxy at 5GHz and 150 MHz and in the rest frame of the observer at 150 MHz. We present the corresponding maps for several epochs after disk formation. Dynamo theory predicts the generation of large‐scale coherent field patterns (“modes”). The timescale of this process is comparable to that of the galaxy age. Many galaxies are expected not to host fully coherent fields at the present epoch, especially those which suffered from major mergers or interactions with other galaxies. A comparison of our predictions with existing observations of spiral galaxies is given and discussed (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent dynamo action in a shear flow

We consider the mean electromotive force and a dynamo-generated magnetic field, taking into account the stretching of turbulent magnetic field lines by a shear flow. Calculations are performed by making use of the second-order correlation approximation. In the presence of shear, the mirror symmetry of turbulence can be broken; thus turbulent motions become suitable for the generation of a large-scale magnetic field. Regardless of the shear law, turbulence can lead to a rapid amplification of the mean magnetic field. The growth rate of the mean magnetic field depends on the length-scale: it is faster for the fields with smaller length-scale. The mechanism considered is qualitatively different from the alpha dynamo, and can generate only a magnetic field that is inhomogeneous in the direction of flow. In contrast to the alpha dynamo, this mechanism also allows the generation of two-dimensional fields. The suggested mechanism may play an important role in the generation of magnetic fields in accretion discs, galaxies and jets.     

Can the turbulent galactic dynamo generate large-scale magnetic fields?

Large-scale magnetic fields in galaxies are thought to be generated by a turbulent dynamo. However, the same turbulence also leads to a small-scale dynamo which generates magnetic noise at a more rapid rate. The efficiency of the large-scale dynamo depends on how this noise saturates. We examine this issue, taking into account ambipolar drift, which obtains in a galaxy with significant neutral gas. We argue as follows.
(i) The small-scale dynamo generated field does not fill the volume, but is concentrated into intermittent rope-like structures. The flux ropes are curved on the turbulent eddy scales. Their thickness is set by the diffusive scale determined by the effective ambipolar diffusion.
(ii) For a largely neutral galactic gas, the small-scale dynamo saturates, as a result of inefficient random stretching, when the peak field in a flux rope has grown to a few times the equipartition value.
(iii) The average energy density in the saturated small-scale field is subequipartition, since it does not fill the volume.
(iv) Such fields neither drain significant energy from the turbulence nor convert eddy motion of the turbulence on the outer scale into wave-like motion. The diffusive effects needed for the large-scale dynamo operation are then preserved until the large-scale field itself grows to near equipartition levels.     

Hydrostatic equilibrium in a magnetized, warped Galactic disc

Hydrostatic equilibrium of the multiphase interstellar medium in the solar vicinity is reconsidered, with the regular and turbulent magnetic fields treated separately. The regular magnetic field strength required to support the gas is consistent with independent estimates, provided that energy equipartition is maintained between turbulence and random magnetic fields. Our results indicate that a mid-plane value of B ₀=4 μG for the regular magnetic field near the Sun leads to more attractive models than B ₀=2 μG . The vertical profiles of both the regular and random magnetic fields contain disc and halo components, the parameters of which we have determined. The layer at 1≲| z |≲4 kpc can be overpressured and an outflow at a speed of about 50 km s⁻¹ may occur there, presumably associated with a Galactic fountain flow, if B ₀≃2 μG .
We show that hydrostatic equilibrium in a warped disc must produce asymmetric density distributions in z , in rough agreement with H  i observations in the outer Galaxy. This asymmetry may be a useful diagnostic of the details of the warping mechanism in the Milky Way and other galaxies. We find indications that gas and magnetic field pressures are different above and below the warped midplane in the outer Galaxy, and quantify the difference in terms of turbulent velocity and/or magnetic field strength.     

Dynamo action in an ambient field

The origin of global magnetic fields in celestial bodies is generally ascribed to dynamo action by fluid motions in their electrically conducting interiors. Some objects – e.g. close‐in extra‐solar planets or the moons of some giant planets – are embedded in ambient magnetic fields which modify the generation of the internal field in these bodies. Recently, the feedback of the magnetospheric field by Chapman‐Ferraro currents in the magnetopause onto the interior dynamo has been proposed to explain the observed weakness of the intrinsic magnetic field of planet Mercury. We study a simplified mean‐field dynamo model which allows us to analytically address various issues like positive and negative feedback situations, stationary versus time‐dependent solutions, and the stability of weak and strong field branches. We discuss the influence of the response function on the solutions when the external field depends on the strength of the intrinsic field like in the situation of the feedback dynamo of Mercury. We find that the feedback mechanism works only for a narrow range of dynamo numbers in the case of Mercury which makes him unique in our solar system. We conclude with some implications for extra‐solar planets (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Predicting solar ‘climate’ by assimilating magnetic data into a flux-transport dynamo

Observational and theoretical knowledge about global-scale solar dynamo ingredients have reached the stage that it is possible to calibrate a flux-transport dynamo for the Sun by adjusting only a few tunable parameters. The important ingredients in this class of model are differential rotation (Omega-effect), helical turbulence (alpha-effect), meridional circulation and turbulent diffusion. The meridional circulation works as a conveyor belt and governs the dynamo cycle period. Meridional circulation and magnetic diffusivity together govern the memory of the Sun's past magnetic fields. After describing the physical processes involved in a flux-transport dynamo, we will show that a predictive tool can be built from it to predict mean solar cycle features by assimilating magnetic field data from previous cycles. We will discuss the theoretical and observational connections among various predictors, such as dynamo-generated toroidal flux integral, cross-equatorial flux, polar fields and geomagnetic indices. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Reconstructing the IRAS Point Source Catalog Redshift Survey with a generalized PIZA

We study the influence of different magnetic boundary conditions on the generation of magnetic fields by turbulent convection. It is found that the structure and strength of the generated field in the vicinity of the boundary is strongly dependent on the choice of boundary conditions. In the convective interior, however, the solutions remain largely insensitive to the boundary conditions. In all cases the overall efficiency of the dynamo process remains high with a steady state magnetic energy density between 12 and 25 per cent of the turbulent kinetic energy, and peak field values exceeding the equipartition level. These results support the idea that the solar granulation may constitute a dynamo source for magnetic fields in the quiet photosphere.     

The relation between magnetic and gas arms in spiral galaxies

The intriguing question of the origin of the arm-like magnetic structures seen between the optical spiral arms of certain spiral galaxies is addressed. Using a two-dimensional approximation to the non-linear disc dynamo equation, it is shown that gas streaming along the arms may produce such a field configuration. Another possibility is a spiral modulation of the turbulent diffusivity, associated with an enhancement of turbulence in the interstellar medium within the arms. The effects of a similar modulation of the alpha-effect are also examined. Finally, the consequences of a non-linear feedback of the large-scale field on the turbulent diffusivity ('η-quenching') are briefly discussed.     

Nonhelical mean‐field dynamos in a sheared turbulence

Mechanisms of nonhelical large‐scale dynamos (shear‐current dynamo and effect of homogeneous kinetic helicity fluctuations with zero mean) in a homogeneous turbulence with large‐scale shear are discussed. We have found that the shearcurrent dynamo can act even in random flows with small Reynolds numbers. However, in this case mean‐field dynamo requires small magnetic Prandtl numbers (i.e., when Pm < Pm^cr < 1). The threshold in the magnetic Prandtl number, Pm^cr = 0.24, is determined using second order correlation approximation (or first‐order smoothing approximation) for a background random flow with a scale‐dependent viscous correlation time τ_c = (νk 2)^–1 (where ν is the kinematic viscosity of the fluid and k is the wave number). For turbulent flows with large Reynolds numbers shear‐current dynamo occurs for arbitrary magnetic Prandtl numbers. This dynamo effect represents a very generic mechanism for generating large‐scale magnetic fields in a broad class of astrophysical turbulent systems with large‐scale shear. On the other hand, mean‐field dynamo due to homogeneous kinetic helicity fluctuations alone in a sheared turbulence is not realistic for a broad class of astrophysical systems because it requires a very specific random forcing of kinetic helicity fluctuations that contains, e.g., low‐frequency oscillations. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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)     

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age‐old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucial role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large‐scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean‐field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some current issues of solar simulations are put forward (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Do mean‐field dynamos in nonrotating turbulent shear‐flows exist?

A plane‐shear flow in a fluid with forced turbulence is considered. If the fluid is electrically‐conducting then a mean electromotive force (EMF) results even without basic rotation and the magnetic diffusivity becomes a highly anisotropic tensor. It is checked whether in this case self‐excitation of a large‐scale magnetic field is possible (so‐called W̄ × J̄ ‐dynamo) and the answer is NO. The calculations reveal the cross‐stream components of the EMF perpendicular to the mean current having the wrong signs, at least for small magnetic Prandtl numbers. After our results numerical simulations with magnetic Prandtl number of about unity have only a restricted meaning as the Prandtl number dependence of the diffusivity tensor is rather strong. If, on the other hand, the turbulence field is strati.ed in the vertical direction then a dynamo‐active α ‐effect is produced. The critical magnetic Reynolds number for such a self‐excitation in a simple shear flow is slightly above 10 like for the other – but much more complicated – flow patterns used in existing dynamo experiments with liquid sodium or gallium. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent dynamo near tachocline and reconstruction of azimuthal magnetic field in the solar convection zone

In order to extend the abilities of the αΩ dynamo model to explain the observed regularities and anomalies of the solar magnetic activity, the negative buoyancy phenomenon and the magnetic quenching of the α effect were included in the model, as well as newest helioseismically determined inner rotation of the Sun were used. Magnetic buoyancy constrains the magnitude of toroidal field produced by the Ω effect near the bottom of the solar convection zone (SCZ). Therefore, we examined two “antibuoyancy” effects: i) macroscopic turbulent diamagnetism and ii) magnetic advection caused by vertical inhomogeneity of fluid density in the SCZ, which we call the ∇ρ effect. The Sun's rotation substantially modifies the ∇ρ effect. The reconstruction of the toroidal field was examined assuming the balance between mean‐field magnetic buoyancy, turbulent diamagnetism and the rotationally modified ∇ρ effect. It is shown that at high latitudes antibuoyancy effects block the magnetic fields in the deep layers of the SCZ, and so the most likely these deep‐rooted fields could not become apparent at the surface as sunspots. In the near‐equatorial region, however, the upward ∇ρ effect can facilitate magnetic fields of about 3000 – 4000 G to emerge through the surface at the sunspot belt. Allowance for the radial inhomogeneity of turbulent velocity in derivations of the helicity parameter resulted in a change of sign of the α effect from positive to negative in the northern hemisphere near the bottom of the SCZ. The change of sign is very important for direction of the Parker's dynamo‐waves propagation and for parity of excited magnetic fields. The period of the dynamo‐wave calculated with allowance for the magnetic quenching is about seven years, that agrees by order of magnitude with the observed mean duration of the sunspot cycles. Using the modern helioseismology data to define dynamo‐parameters, we conclude that north‐south asymmetry should exist in the meridional field. At low latitudes in deep layers of the SCZ, the αΩ dynamo excites most efficiency the dipolar mode of the meridional field. Meanwhile, in high‐latitude regions a quadrupolar mode dominates in the meridional field. The obtained configuration of the net meridional field is likely to explain the magnetic anomaly of polar fields (the apparent magnetic “monopole”) observed near the maxima of solar cycles. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo-generated magnetic fields at the surface of a massive star

Spruit has shown that an astrophysical dynamo can operate in the non-convective material of a differentially rotating star as a result of a particular instability in the magnetic field (the Tayler instability). By assuming that the dynamo operates in a state of marginal instability, Spruit has obtained formulae which predict the equilibrium strengths of azimuthal and radial field components in terms of local physical quantities. Here, we apply Spruit's formulae to our previously published models of rotating massive stars in order to estimate Tayler dynamo field strengths. There are no free parameters in Spruit's formulae. In our models of 10- and  50-M_⊙  stars on the zero-age main sequence, we find internal azimuthal fields of up to 1 MG, and internal radial components of a few kG. Evolved models contain weaker fields. In order to obtain estimates of the field strength at the stellar surface, we examine the conditions under which the Tayler dynamo fields are subject to magnetic buoyancy. We find that conditions for Tayler instability overlap with those for buoyancy at intermediate to high magnetic latitudes. This suggests that fields emerge at the surface of a massive star between magnetic latitudes of about 45° and the poles. We attempt to estimate the strength of the field which emerges at the surface of a massive star. Although these estimates are very rough, we find that the surface field strengths overlap with values which have been reported recently for line-of-sight fields in several O and B stars.     

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)