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)     

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)     

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)     

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)     

Axisymmetric field generation within an ambient axial field

The generation of magnetic field in a homogeneous, electrically conducting fluid – as required for the dynamo generation of the fields of many astrophysical bodies – is normally a threshold process; the dynamo mechanism, applicable to such bodies in unmagnetised environments, requires motions of sufficient strength to overcome the innate magnetic diffusion. In the presence of an ambient field, however, the critical nature of the field generation process is relaxed. Motions can distort and amplify the ambient field for all amplitudes of flow. For motions with appropriate geometries, an internal ‘dynamo‐like’ field of appreciable strength can be generated, even for relatively weak flows. At least a minority of planets, moons and other bodies exist within significant external astrophysical fields. For these bodies, the ambient field problem is more relevant than the classical dynamo problem, yet it remains relatively little studied. In this paper we consider the effect of an axial ambient field on a spherical mean‐field α ²ω dynamo model, through nonlinear calculations with α ‐quenching feedback. Ambient fields of varying strengths, and both stationary and oscillatory in time, are imposed. Particular focus is placed on the effects of these fields on the equatorial symmetry and the time dependence of the preferred solutions. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Can turbulent plane‐shear flows support large‐scale dynamos?

Turbulent plane‐shear flow is found to show same basic effects of mean‐fieldMHD as rotating turbulence. In particular, the mean electromotive force (EMF) includes highly anisotropic turbulent diffusion and alpha‐effect. Only magnetic diffusion remains for spatially‐uniform turbulence. The question is addressed whether in this case a self‐excitation of a magnetic field by so‐called sher‐current dynamo is possible and the quasilinear theory provides a negative answer. The streamaligned component of the EMF has the sign opposite to that required for dynamo. If, however, the turbulence is not uniform across the flow direction then a dynamo‐active α ‐effect emerges. The critical magnetic Reynolds number for the alpha‐shear dynamo is estimated to be slightly above ten. Possibilities for cross‐checking theoretical predictions with MHD experiments are discussed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling spatio‐temporal nonlocality in mean‐field dynamos

When scale separation in space or time is poor, the mean‐field α effect and turbulent diffusivity have to be replaced by integral kernels by which the dependence of the mean electromotive force on the mean magnetic field becomes nonlocal. Earlier work in computing these kernels using the test‐field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel for isotropic stationary turbulence is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean‐field equations are solved for oscillatory α –shear dynamos as well as α² dynamos with α linearly depending on position, which makes this dynamo oscillatory, too. In both cases, the critical values of the dynamo number is lowered due to spatio‐temporal nonlocality.When scale separation in space or time is poor, the mean‐field α effect and turbulent diffusivity have to be replaced by integral kernels by which the dependence of the mean electromotive force on the mean magnetic field becomes nonlocal. Earlier work in computing these kernels using the test‐field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel for isotropic stationary turbulence is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean‐field equations are solved for oscillatory α –shear dynamos as well as α² dynamos     

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)     

Dynamo saturation in direct simulations of the multi‐phase turbulent interstellar medium

The ordered magnetic field observed via polarised synchrotron emission in nearby disc galaxies can be explained by a mean‐field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical‐flux initial conditions are potentially able to influence this dynamo via the occurrence of the magnetorotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean‐field dynamo. This is motivated by the observation that different saturation behaviour was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three‐dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing‐periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean‐field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched α) dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxy's star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub‐grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields. (© 2015 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)     

Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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 2×10¹¹ cm² s⁻¹. The quenching of α, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong α-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.8×10¹² cm² s⁻¹ 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.     

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)     

Non‐holonomic filamentary dynamos as generalized Arnold's maps in Riemannian space

A filamentary non‐holonomic dynamo solution of self‐induction magnetic field equation is found by considering highly conducting filaments. It is shown that planar filaments cannot support dynamo action since the flow along the filament vanishes for torsion‐free filaments. This is a generalization of the Zeldovich theorem for linear magnetic dynamo filaments. The flow of filament is proportionally to the product between Frenet torsion and curvature. This shows that filamentary dynamos must possess Frenet torsion. A well‐known example of this result is the α ‐dynamo in solar physics. Magnetic helicity and magnetic energy for this filamentary dynamo are computed. Magnetic helicity vanishes by construction and the magnetic field decays with torsion energy in helicoidal dynamos. The approach considered here is useful for the investigation of anisotropic turbulent cascades. As a particular simple example it is shown that under certain constraints the solution can be reduced to the Arnold cat dynamo map solution where the non‐holonomic directional mixed derivative, would play the role of the Lyapunov exponent which appears on stretching the magnetic field in Riemannian space. The solution seems to describe marginal slow dynamos when the velocities involved in the dynamo flows are constants. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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)     

Self‐similar evolutionary solutions for an accreting magneto‐fluid around a compact object with finite electrical conductivity

In this paper, we investigate the time evolution of an accreting magneto‐fluid with finite conductivity. For the case of a thin disk, the fluid equations along with Maxwell's equations are derived in a simplified, one‐dimensional model that neglects the latitudinal dependence of the flow. The finite electrical conductivity of the plasma is taken into account by Ohm's law; however, the shear viscous stress is neglected, as well as the self‐gravity of the disk. In order to solve the integrated equations that govern the dynamical behaviour of the magneto‐fluid, we have used a self‐similar solution. We introduce two dimensionless variables, S₀ and ε_ϱ, which represent the size of the electrical conductivity and the density behaviour with time, respectively. The effect of each of these on the structure of the disk is studied. While the pressure is obtained simply by solving an ordinary differential equation, the density, the magnetic field, the radial velocity, and the rotational velocity are presented analytically. The solutions show that the S0 and ε_ϱ parameters affect the radial thickness of the disk. Also, radial velocity and gas pressure are more sensitive to the electrical conductivity in the inner regions of disk. Moreover, the parameter ε_ϱ has a more significant effect on the physical quantities for small radii. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Differential rotation and meridional flow of Arcturus

The spectroscopic variability of Arcturus hints at cyclic activity cycle and differential rotation. This could provide a test of current theoretical models of solar and stellar dynamos. To examine the applicability of current models of the flux transport dynamo to Arcturus, we compute a mean‐field model for its internal rotation, meridional flow, and convective heat transport in the convective envelope. We then compare the conditions for dynamo action with those on the Sun. We find solar‐type surface rotation with about 1/10th of the shear found on the solar surface. The rotation rate increases monotonically with depth at all latitudes throughout the whole convection zone. In the lower part of the convection zone the horizontal shear vanishes and there is a strong radial gradient. The surface meridional flow has maximum speed of 170 m/s and is directed towards the equator at high and towards the poles at low latitudes. Turbulent magnetic diffusivity is of the order 10¹⁵–10¹⁶ cm²/s. The conditions on Arcturus are not favorable for a circulation‐dominated dynamo (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Kinematic Numerical Simulations of the Solar Dynamo: Dependence on α and Ω Values

We investigated the generation of dynamo waves in the solar convection zone through a numerical simulation. We integrated the axisymmetric α–Ω kinematic dynamo equations in a spherical geometry, where the α- and Ω-profiles depend on the spatial coordinates. The model results show that the fundamental parameter that determines the behavior of the system is the product between the characteristic intensities of the α and Ω contributions. In particular, we found three different regimes in which the system exhibits different behaviors: a regime without a dynamo effect, one with an exponential amplification of the magnetic field, and one with dynamo waves.     

A two-layer αω-dynamo model with dynamic feedback on the ω-effect

In an attempt to produce a simple representation of an interface dynamo, I examine a dynamo model composed of two one-dimensional (radially averaged) pseudo-spherical layers, one in the convection zone and possessing an α-effect, and the other in the tachocline and possessing an ω-effect. The two layers communicate by means of an analogue of Newton's law of cooling, and a dynamical back-reaction of the magnetic field on ω is provided. Extensive bifurcation diagrams are calculated for three separate values of η, the ratio of magnetic diffusivities of the two layers. I find recognizable similarities to, but also dramatic differences from, the comparable one-layer model examined by Roald &38; Thomas. In particular, the solar-like dynamo mode found previously is no longer stable in the two-layer version; in its place there is a sequence of periodic, quasi-periodic and chaotic modes probably created in a homoclinic bifurcation. These differences are important enough to provide support for the view that the solar dynamo cannot be meaningfully modelled in one dimension.     

The high redshift intergalactic and interstellar media with X‐shooter

I review here a few important questions that X‐shooter can help tackle and answer in the field of quasar absorption lines. This includes (i) determine the ionizing background and the physical state of the inter‐galactic medium (IGM) by analysing the characteristics of the Lyman‐α forest and the proximity effect; (ii) investigate the metal content of the high redshift IGM; (iii) study the small scale transverse correlation in the IGM by observing pairs of quasars with small separation in the sky; (iv) study the galaxy‐IGM relations by detecting the counterpart of damped Lyman‐α systems (DLAs) or determining the correlation between the properties of galaxies and absorption lines; (v) detect and characterize the long‐sought cold diffuse molecular (H₂ and CO) interstellar medium (ISM) of high redshift galaxies and study its dust content (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)