Laboratory plasma dynamos, astrophysical dynamos and magnetic helicity evolution

The term 'dynamo' means different things to the laboratory fusion plasma and astrophysical plasma communities. To alleviate the resulting confusion and to facilitate interdisciplinary progress, we pinpoint conceptual differences and similarities between laboratory plasma dynamos and astrophysical dynamos. We can divide dynamos into three types: 1. magnetically dominated helical dynamos which sustain a large-scale magnetic field against resistive decay and drive the magnetic geometry towards the lowest energy state, 2. flow-driven helical dynamos which amplify or sustain large-scale magnetic fields in an otherwise turbulent flow and 3. flow-driven non-helical dynamos which amplify fields on scales at or below the driving turbulence. We discuss how all three types occur in astrophysics whereas plasma confinement device dynamos are of the first type. Type 3 dynamos require no magnetic or kinetic helicity of any kind. Focusing on Types 1 and 2 dynamos, we show how different limits of a unified set of equations for magnetic helicity evolution reveal both types. We explicitly describe a steady-state example of a Type 1 dynamo, and three examples of Type 2 dynamos: (i) closed volume and time dependent; (ii) steady state with open boundaries; (iii) time dependent with open boundaries.     

The solar dynamo

The solar dynamo continues to pose a challenge to observers and theoreticians. Observations of the solar surface reveal a magnetic field with a complex, hierarchical structure consisting of widely different scales. Systematic features such as the solar cycle, the butterfly diagram, and Hale's polarity laws point to the existence of a deep-rooted large-scale magnetic field. At the other end of the scale are magnetic elements and small-scale mixed-polarity magnetic fields. In order to explain these phenomena, dynamo theory provides all the necessary ingredients including the effect, magnetic field amplification by differential rotation, magnetic pumping, turbulent diffusion, magnetic buoyancy, flux storage, stochastic variations and nonlinear dynamics. Due to advances in helioseismology, observations of stellar magnetic fields and computer capabilities, significant progress has been made in our understanding of these and other aspects such as the role of the tachocline, convective plumes and magnetic helicity conservation. However, remaining uncertainties about the nature of the deep-seated toroidal magnetic field and the effect, and the forbidding range of length scales of the magnetic field and the flow have thus far prevented the formulation of a coherent model for the solar dynamo. A preliminary evaluation of the various dynamo models that have been proposed seems to favor a buoyancy-driven or distributed scenario. The viewpoint proposed here is that progress in understanding the solar dynamo and explaining the observations can be achieved only through a combination of approaches including local numerical experiments and global mean-field modeling.Received: 5 May 2003, Published online: 15 July 2003     

A model of magnetically induced disc–corona for black hole binaries

We propose a model of magnetic connection (MC) of a black hole with its surrounding accretion disc based on large-scale magnetic field. The MC gives rise to transport of energy and angular momentum between the black hole and the disc, and the closed field lines pipe the hot matter evaporated from the disc, and shape it in the corona above the disc to form a magnetically induced disc–corona system, in which the corona has the same configuration as the large-scale magnetic field. We numerically solve the dynamic equations in the context of the Kerr metric, in which the large-scale magnetic field is determined by dynamo process and equipartition between magnetic pressure and gas pressure. Thus we can obtain a global solution rather than assuming the distribution of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) The rotational energy of a fast-spinning black hole can be extracted, enhancing the dissipation in the accretion disc, (2) the closed field lines provide a natural channel for corona matter escaping from disc and finally falling into black hole and (3) the scope of the corona can be bounded by the conservation of magnetic flux. We simulate the high-energy spectra of this system by using Monte Carlo method, and find that the relative hardness of the spectra decreases as accretion rate or black hole spin a _* increases. We fit the typical X-ray spectra of three black hole binaries  (GRO J1655−40, XTE 1118+480 and GX 339−4)  in the low/hard or very high state.     

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.     

Non-ideal evolution of non-axisymmetric, force-free magnetic fields in a magnetar

Recent numerical magnetohydrodynamic calculations by Braithwaite and collaborators support the 'fossil field' hypothesis regarding the origin of magnetic fields in compact stars and suggest that the resistive evolution of the fossil field can explain the reorganization and decay of magnetar magnetic fields. Here, these findings are modelled analytically by allowing the stellar magnetic field to relax through a quasi-static sequence of non-axisymmetric, force-free states, by analogy with spheromak relaxation experiments, starting from a random field. Under the hypothesis that the force-free modes approach energy equipartition in the absence of resistivity, the output of the numerical calculations is semiquantitatively recovered: the field settles down to a linked poloidal–toroidal configuration, which inflates and becomes more toroidal as time passes. A qualitatively similar (but not identical) end state is reached if the magnetic field evolves by exchanging helicity between small and large scales according to an α-dynamo-like, mean-field mechanism, arising from the fluctuating electromotive force produced by the initial random field. The impossibility of matching a force-free internal field to a potential exterior field is discussed in the magnetar context.     

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.     

Mean-field dynamo in partially ionized plasmas – I

There are several astrophysical situations where one needs to study the dynamics of magnetic flux in partially ionized turbulent plasmas. In a partially ionized plasma, the magnetic induction is subjected to the ambipolar diffusion and the Hall effect in addition to the usual resistive dissipation. In this paper, we initiate the study of the kinematic dynamo in a partially ionized turbulent plasma. The Hall effect arises from the treatment of the electrons and the ions as two separate fluids and the ambipolar diffusion due to the inclusion of neutrals as the third fluid. It is shown that these non-ideal effects modify the so-called α effect and the turbulent diffusion coefficient β in a rather substantial way. The Hall effect may enhance or quench the dynamo action altogether. The ambipolar diffusion brings in an α which depends on the mean magnetic field. The new correlations embodying the coupling of the charged fluids and the neutral fluid appear in a decisive manner. The turbulence is necessarily magnetohydrodynamic with new spatial and time-scales. The nature of the new correlations is demonstrated by taking the Alfvénic turbulence as an example.     

Accretion disc coronae as magnetic reservoirs

Most astrophysical sources powered by accretion on to a black hole, either of stellar mass or supermassive, when observed with hard X-rays show signs of a hot Comptonizing component in the flow, the so-called corona , with observed temperatures and optical depths lying in a narrow range (0.1≲ τ ≲1 and 1×10⁹ K≲ T ≲3×10⁹ K). Here we argue that these facts constitute strong supporting evidence for a magnetically dominated corona. We show that the inferred thermal energy content of the corona, in all black hole systems, is far too low to explain their observed hard X-ray luminosities, unless either the size of the corona is at least of the order of 10³ Schwarzschild radii, or the corona itself is in fact a reservoir , where the energy is mainly stored in the form of a magnetic field generated by a sheared rotator (probably the accretion disc). We briefly outline the main reasons why the former possibility is to be discarded, and the latter preferred.     

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)     

Comparisons and connections between mean field dynamo theory and accretion disc theory

The origin of large scale magnetic fields in astrophysical rotators, and the conversion of gravitational energy into radiation near stars and compact objects via accretion have been subjects of active research for a half century. Magnetohydrodynamic turbulence makes both problems highly nonlinear, so both subjects have benefitted from numerical simulations.However, understanding the key principles and practical modeling of observations warrants testable semi‐analytic mean field theories that distill the essential physics. Mean field dynamo (MFD) theory and alpha‐viscosity accretion disc theory exemplify this pursuit. That the latter is a mean field theory is not always made explicit but the combination of turbulence and global symmetry imply such. The more commonly explicit presentation of assumptions in 20th century textbook MFDT has exposed it to arguably more widespread criticism than incurred by 20th century alpha‐accretion theory despite complementary weaknesses. In the 21st century however, MFDT has experienced a breakthrough with a dynamical saturation theory that consistently agrees with simulations. Such has not yet occurred in accretion disc theory, though progress is emerging. Ironically however, for accretion engines, MFDT and accretion theory are presently two artificially uncoupled pieces of what should be a single coupled theory. Large scale fields and accretion flows are dynamically intertwined because large scale fields likely play a key role in angular momentum transport. I discuss and synthesize aspects of recent progress in MFDT and accretion disc theory to suggest why the two likely conspire in a unified theory (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic field bending in accretion discs with dynamos

We consider the problem of poloidal magnetic field advection and bending of an initially vertical field owing to radial inflow in thin accretion discs. For a ratio of kinematic viscosity to magnetic diffusivity of order unity, significant bending of an externally applied vertical field cannot occur in a disc with no internal dynamo. However, we show that if poloidal field is generated by a dynamo operating near its critical state, then significant field bending may be possible. Our results are of particular relevance to wind launching from accretion discs.     

Small-scale magnetic helicity losses from a mean-field dynamo

Using mean-field models with a dynamical quenching formalism, we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity density or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong for both large- and small-scale fields. In the presence of shear, there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region, there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behaviour emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain.     

Relating stellar cycle periods to dynamo calculations

Stellar magnetic activity in slowly rotating stars is often cyclic, with the period of the magnetic cycle depending critically on the rotation rate and the convective turnover time of the star. Here we show that the interpretation of this law from dynamo models is not a simple task. It is demonstrated that the period is (unsurprisingly) sensitive to the precise type of non-linearity employed. Moreover the calculation of the wave-speed of plane-wave solutions does not (as was previously supposed) give an indication of the magnetic period in a more realistic dynamo model, as the changes in length-scale of solutions are not easily captured by this approach. Progress can be made, however, by considering a realistic two-dimensional model, in which the radial length-scale of waves is included. We show that it is possible in this case to derive a more robust relation between cycle period and dynamo number. For all the non-linearities considered in the most realistic model, the magnetic cycle period is a decreasing function of | D | (the amplitude of the dynamo number). However, discriminating between different non-linearities is difficult in this case and care must therefore be taken before advancing explanations for the magnetic periods of stars.     

Intermittent Behaviour in Axisymmetric Mean-Field Dynamo Models in Spherical Shells

Axisymmetric mean-field dynamo models in spherical shells are shown to be capable of producing temporally intermittent behaviour. This is of potential importance since (i) it is, as far as we are aware, the first time such behaviour has been produced internally by a mean-field dynamo model in a spherical shell, without requiring any additional assumptions or truncations, and (ii) it may be characteristic of the type of behaviour observed in the long-term record of solar activity, such as Maunder minima. We also show that these types of behaviour persist when the functional form of the alpha quenching is altered and also occur over intervals of the shell thickness and the dynamo number.     

Dimensionless measures of turbulent magnetohydrodynamic dissipation rates

The magnetic Reynolds number, R _M, is defined as the product of a characteristic scale and associated flow speed divided by the microphysical magnetic diffusivity. For laminar flows, R _M also approximates the ratio of advective to dissipative terms in the total magnetic energy equation, but for turbulent flows this latter ratio depends on the energy spectra and approaches unity in a steady state. To generalize for flows of arbitrary spectra we define an effective magnetic dissipation number,   R _M,e  , as the ratio of the advection to microphysical dissipation terms in the total magnetic energy equation, incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, and distinct pairs of inner and outer scales for magnetic and kinetic spectra. As expected, for a substantial parameter range   R _M,e∼ O (1) ≪ R _M  . We also distinguish   R _M,e  from     where the latter is an effective magnetic Reynolds number for the mean magnetic field equation when a turbulent diffusivity is explicitly imposed as a closure. That   R _M,e  and     approach unity even if   R _M≫ 1  highlights that, just as in hydrodynamic turbulence, energy dissipation of large-scale structures in turbulent flows via a cascade can be much faster than the dissipation of large-scale structures in laminar flows. This illustrates that the rate of energy dissipation by magnetic reconnection is much faster in turbulent flows, and much less sensitive to microphysical reconnection rates compared to laminar flows.     

On accretion flow penetration of magnetospheres

We address the problem of plasma penetration of astrophysical magnetospheres, an important issue in a wide variety of contexts, ranging from accretion in cataclysmic variables to flows in protostellar systems. We point out that under well-defined conditions, penetration can occur without any turbulent mixing (driven, for example, by Rayleigh–Taylor or Kelvin–Helmholtz instabilities) caused by charge polarization effects, if the inflowing plasma is bounded in the direction transverse to both the flow velocity and the magnetic field. Depolarization effects limit the penetration depth, which nevertheless can, under specific circumstances, be comparable to the size of the magnetosphere. We discuss the effect of ambient medium on plasma propagation across the stellar magnetic field and determine the criteria for deep magnetosphere penetration. We show that, under conditions appropriate to magnetized white dwarfs in AM Her type cataclysmic variables, charge polarization effects can lead to deep penetration of the magnetosphere.     

On the meaning and inapplicability of the Zeldovich relations of magnetohydrodynamics

Considering a plasma with an initially weak large scale field subject to nonhelical turbulent stirring, Zeldovich (1957), for two‐dimensions, followed by others for three dimensions, have presented formulae of the form 〈b²〉 = f(R_M) . Such “Zeldovich relations” have sometimes been interpreted to provide steady‐state relations between the energy associated with the fluctuating magnetic field and that associated with a large scale or mean field multiplied by a function f that depends on spatial dimension and a magnetic Reynolds number R_M. Here we dissect the origin of these relations and pinpoint pitfalls that show why they are inapplicable to realistic, dynamical MHD turbulence and that they disagree with many numerical simulations. For 2D, we show that when the total magnetic field is determined by a vector potential, the standard Zeldovich relation applies only transiently, characterizing a maximum possible value that the field energy can reach before necessarily decaying. In 3D, we show that the standard Zeldovich relations are derived by balancing subdominant terms. In contrast, balancing the dominant terms shows that the fluctuating field can grow to a value independent of R_M and the initially imposed , as seen in numerical simulations. We also emphasize that these Zeldovich relations of nonhelical turbulence imply nothing about the amount mean field growth in a helical dynamo. In short, by re‐analyzing the origin of the Zeldovich relations, we highlight that they are inapplicable to realistic steady‐states of large R_M MHD turbulence. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo action in the presence of an imposed magnetic field

Dynamo action within the cores of Ap stars may offer intriguing possibilities for understanding the persistent magnetic fields observed on the surfaces of these stars. Deep within the cores of Ap stars, the coupling of convection with rotation likely yields magnetic dynamo action, generating strong magnetic fields. However, the surface fields of the magnetic Ap stars are generally thought to be of primordial origin. Recent numerical models suggest that a primordial field in the radiative envelope may possess a highly twisted toroidal shape. We have used detailed 3-D simulations to study the interaction of such a twisted magnetic field in the radiative envelope with the core-dynamo operating in the interior of a 2 solar mass A-type star. The resulting dynamo action is much more vigorous than in the absence of such a fossil field, yielding magnetic field strengths (of order 100 kG) much higher than their equipartition values relative to the convective velocities. We examine the generation of these fields, as well as the growth of large-scale magnetic structure that results from imposing a fossil magnetic field. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A comparison of DA white dwarf temperatures and gravities from Lyman and Balmer line studies

The stability of turbulent accretion discs is considered, in which a magnetically influenced wind plays a major role in driving the inflow. The magnetic field is generated by a dynamo operating in the disc, involving radial shear and turbulence. The steady angular momentum balance is found to be linearly stable for a range of radial boundary conditions, and an expression is derived for the adjustment time-scale as a function of the equilibrium ratio of the magnetic and viscous disc torques.     

An unconventional accretion disc dynamo

In the light of recent results from numerical simulations of accretion disc MHD turbulence, we revisit the problem of the configuration of large-scale magnetic fields resulting from an α Ω dynamo operating in a thin accretion disc. In particular, we analyse the consequences of the peculiar sign of the α -effect suggested by numerical simulations . We determine the symmetry of the fastest-growing modes in the kinematic dynamo approximation and, in the framework of an ' α -quenched' dynamo model, study the evolution of the magnetic field. We find that the resulting field for this negative polarity of the α -effect generally has dipole symmetry with respect to the disc midplane, although the existence of an equilibrium configuration depends on the properties of the turbulence. The role of magnetic field dragging is discussed and, finally, the presence of an external uniform magnetic field is included to address the issue of magneto centrifugal wind launching from accretion discs.