Subcritical dynamos in shear flows

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non‐rotating wall‐bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non‐axisymmetric instabilities of shearinduced axisymmetric toroidal velocity or magnetic fields, such as Kelvin‐Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high‐Reynolds number shear flows could act as a large‐scale driving mechanism of turbulent flows that would in turn generate an independent small‐scale dynamo. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of Ohmic diffusion on the excitation and dynamics of MRI

In this paper we make an effort to understand the interaction of turbulence generated by the magnetorotational instability (MRI) with turbulence from other sources, such as supernova explosions (SNe) in galactic disks. First we perform a linear stability analysis (LSA) of non‐ideal MRI to derive the limiting value of Ohmic diffusion that is needed to inhibit the growth of the instability for different types of rotation laws. With the help of a simple analytical expression derived under first‐order smoothing approximation (FOSA), an estimate of the limiting turbulence level and hence the turbulent diffusion needed to damp the MRI is derived. Secondly, we perform numerical simulations in local cubes of isothermal nonstratified gas with external forcing of varying strength to see whether the linear result holds for more complex systems. Purely hydrodynamic calculations with forcing, rotation and shear are made for reference purposes, and as expected, non‐zero Reynolds stresses are found. In the magnetohydrodynamic calculations, therefore, the total stresses generated are a sum of the forcing and MRI contributions. To separate these contributions, we perform reference runs with MRI‐stable shear profiles (angular velocity increasing outwards), which suggest that the MRI‐generated stresses indeed become strongly suppressed as function of the forcing. The Maxwell to Reynolds stress ratio is observed to decrease by an order of magnitude as the turbulence level due to external forcing exceeds the predicted limiting value, which we interpret as a sign of MRI suppression. Finally, we apply these results to estimate the limiting radius inside of which the SN activity can suppress the MRI, arriving at a value of 14 kpc (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Current‐driven instabilities in weakly ionized disks

Cool weakly ionized gaseous rotating disk form the basis for many models in astrophysics objects. Instabilities against perturbations in such disks play an important role in the theory of the formation of stars and planets. Traditionally, axisymmetric magnetohydrodynamic (MHD) and recently Hall‐MHD instabilities have been thoroughly studied as providers of an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, the Hall instability against axisymmetric perturbations in incompressible rotating fluid in external poloidal and toroidal magnetic field is considered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hopf bifurcations with fluctuating gain

A Hopf bifurcation describes a transition from a steady state to an oscillatory one. This report describes how the presence of white Gaussian noise in the term driving the instability (the gain) affects the character of such a transition. We find that for small noise intensities the effective critical value of the bifurcation parameter increases linearly with the noise intensity. This result is used to quantitatively describe the onset of the magneto‐rotational instability in a laboratory experiment using liquid metal with turbulent background flow. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Saturation of the magnetorotational instability at large Elsasser number

The magnetorotational instability is investigated within the shearing box approximation in the large Elsasser number regime. In this regime, which is of fundamental importance to astrophysical accretion disk theory, shear is the dominant source of energy, but the instability itself requires the presence of a weaker vertical magnetic field. Dissipative effects are weaker still but not negligible. The regime explored retains the condition that (viscous and ohmic) dissipative forces do not play a role in the leading order linear instability mechanism. However, they are sufficiently large to permit a nonlinear feedback mechanism whereby the turbulent stresses generated by the MRI act on and modify the local background shear in the angular velocity profile. To date this response has been omitted in shearing box simulations and is captured by a reduced pde model derived here from the global MHD fluid equations using multiscale asymptotic perturbation theory. Results from numerical simulations of the reduced pde model indicate a linear phase of exponential growth followed by a nonlinear adjustment to algebraic growth and decay in the fluctuating quantities. Remarkably, the velocity and magnetic field correlations associated with these algebraic growth and decay laws conspire to achieve saturation of the angular momentum transport. The inclusion of subdominant ohmic dissipation arrests the algebraic growth of the fluctuations on a longer, dissipative time scale. (© 2008 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)     

Secondary tearingmode in the nonlinear evolution of magnetorotational instability

Numerical investigation of the two‐dimensional magnetic reconnection is given in the context of the nonlinear evolution of the Magneto‐Rotational Instability (MRI). With a careful comparison to various theories using both one‐ and twodimensional analysis, it is found that a new stabilizing effect of the centrifugal force on tearing instability must be present in the specific geometry of the MRI.Magnetic reconnection might play a key role to the formation of the nonaxisymmetric structures observed in MRI experiments. The results may also be useful for the estimate of the accretion rate in various astrophysical objects. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulence and its parameterization in accretion discs

Accretion disc turbulence is investigated in the framework of the shearing box approximation. The turbulence is either driven by the magneto‐rotational instability or, in the non‐magnetic case, by an explicit and artificial forcing term in the momentum equation. Unlike the magnetic case, where most of the dissipation occurs in the disc corona, in the forced hydrodynamic case most of the dissipation occurs near the midplane. In the hydrodynamic case evidence is presented for the stochastic excitation of epicycles. When the vertical and radial epicyclic frequencies are different (modeling the properties around rotating black holes), the beat frequency between these two frequencies appear to show up as a peak in the temporal power spectrum in some cases. Finally, the full turbulent resistivity tensor is determined and it is found that, if the turbulence is driven by a forcing term, the signs of its off‐diagonal components are such that this effect would not be capable of dynamo action by the shear–current effect. (© 2005 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)     

Properties of the negative effective magnetic pressure instability

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large‐scale instability. In this study, hydromagnetic mean‐field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two‐dimensional. The growth rate is found to depend on a parameter β_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter β_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The origin of our galactic magnetic field

A serious difficulty with the standard alpha‐omega theory of the origin of galactic magnetic fields involves the question of flux expulsion. This is intimately related to flux freezing. The alpha‐omega theory is shown in the context of the giant superbubble explosions that have a large impact on the physics of the interstellar medium. It is shown that superbubbles alone can duplicate the processes of the alpha‐omega dynamo and produce exponential growth of the galactic magnetic field. The possibility of the blow‐out of pieces of the magnetic field is discussed and it is shown that they have the potential to solve the flux‐expulsion problem. However, such an explanation must lead to apparent ‘gaps’ in the field in the galactic disc. These gaps are probably unavoidable in any dynamo theory and should have important observable consequences, one of which is an explanation for the escape of cosmic rays from the disc (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Shear‐driven magnetic buoyancy oscillations

The effects of uniform horizontal shear on a stably stratified layer of gas is studied. The system is initially destabilized by a magnetically buoyant flux tube pointing in the cross‐stream direction. The shear amplifies the initial field to Lundquist numbers of about 200–400, but then its value drops to about 100–300, depending on the value of the sub‐adiabatic gradient. The larger values correspond to cases where the stratification is strongly stable and nearly isothermal. At the end of the runs the magnetic field is nearly axisymmetric, i.e. uniform in the streamwise direction. In view of Cowling's theorem the sustainment of the field remains a puzzle and may be due to subtle numerical effects that have not yet been identified in detail. In the final state the strength of the magnetic field decreases with height in such a way that the field is expected to be unstable. Low amplitude oscillations are seen in the vertical velocity even at late times, suggesting that they might be persistent (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermomagnetic and magnetorotational instabilities of thin Keplerian discs in poloidal magnetic fields

Stability of thin hot Keplerian discs is investigated asymptotically in small disc's aspect ratio, ε. The study is carried out in the local approximation for short vertical waves in the disc‐thickness scale. Besides the radial rotation shear and the vertical magnetic field, the background configuration is characterized by a vertically near‐constant temperature profile with a small vertical gradient. The temperature‐gradient term in Ohm's law, which characterizes the thermomagnetic transport is found to be of the order of ε. The effect of the thermomagnetic transport slightly modifies the conventional magnetorotational instability (MRI), while a new thermomagnetic instability (TMI) emerges in regions of the wavenumber space where MRI is absent. Explicit solutions are obtained for a wide range of values of plasma beta, β, and thermomagnetic transport coefficient, λ. In particular, it is shown for λ ≪ 1 that the MRI dominates in weak magnetic fields, β ≫ 1, while the TMI is exhibited in strong magnetic fields, β ∼ 1, also with the growth rate of the order of inverse rotation period (© 2011 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)     

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, the amplitude of which depends on the vertical position, z , in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays, we can obtain information about the interaction between the warp and the disc turbulence. A high-amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, α _v, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, α _h, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, α _v is approximately equal to α _h and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness (∼0.01) of α _v and α _h we conclude that for β p _gas p _mag∼10 the time-scale for diffusion or damping of a warp is much shorter than the usual viscous time-scale. Finally, we review the astrophysical implications.     

Linear stability analysis of the Hall magnetorotational instability in a spherical domain

We investigate the stability of the Hall‐MHD system and determine its importance for neutron stars at their birth, when they still consist of differentially rotating plasma permeated by extremely strong magnetic fields. We solve the linearised HallMHD equations in a spherical shell threaded by a homogeneous magnetic field. With the fluid/flow coupling and the Hall effect included, the magnetorotational instability and the Hall effect are both acting together. Results differ for magnetic fields aligned with the rotation axis and anti‐parallel magnetic fields. For a positive alignment of the magnetic field the instability grows on a rotational time‐scale for any sufficiently large magnetic Reynolds number. Even the magnetic fields which are stable against the MRI due to the magnetic diffusion are now susceptible to the shear‐Hall instability. In contrast, the negative alignment places strong restrictions on the growth and the magnitude of the fields, hindering the effectiveness of the Hall‐MRI. While non‐axisymmetric modes of the MRI can be suppressed by strong enough rotation, there is no such restriction when the Hall effect is present. The implications for the magnitude and the topology of the magnetic field of a young neutron star may be significant (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mean electromotive force proportional to mean flow in MILD turbulence

In mean‐field magnetohydrodynamics the mean electromotive force due to velocity and magnetic‐field fluctuations plays a crucial role. In general it consists of two parts, one independent of and another one proportional to the mean magnetic field. The first part may be nonzero only in the presence of mhd turbulence, maintained, e.g., by small‐scale dynamo action. It corresponds to a battery, which lets a mean magnetic field grow from zero to a finite value. The second part, which covers, e.g., the α effect, is important for large‐scale dynamos. Only a few examples of the aforementioned first part of the mean electromotive force have been discussed so far. It is shown that a mean electromotive force proportional to the mean fluid velocity, but independent of the mean magnetic field, may occur in an originally homogeneous isotropic mhd turbulence if there are nonzero correlations of velocity and electric current fluctuations or, what is equivalent, of vorticity and magnetic field fluctuations. This goes beyond the Yoshizawa effect, which consists in the occurrence of mean electromotive forces proportional to the mean vorticity or to the angular velocity defining the Coriolis force in a rotating frame and depends on the cross‐helicity defined by the velocity and magnetic field fluctuations. Contributions to the mean electromotive force due to inhomogeneity of the turbulence are also considered. Possible consequences of the above findings for the generation of magnetic fields in cosmic bodies are discussed (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)