Mach number dependence of turbulent magnetic field amplification: solenoidal versus compressive flows

Computational approaches to modeling dynamos in galaxies

Galaxies are observed to host magnetic fields with a typical total strength of around 15  μ G. A coherent large-scale field constitutes up to a few microgauss of the total, while the rest is built from strong magnetic fluctuations over a wide range of spatial scales. This represents sufficient magnetic energy for it to be dynamically significant. Several questions immediately arise: What is the physical mechanism that gives rise to such magnetic fields? How do these magnetic fields affect the formation and evolution of galaxies? In which physical processes do magnetic fields play a role, and how can that role be characterized? Numerical modelling of magnetized flows in galaxies is playing an ever-increasing role in finding those answers. We review major techniques used for these models. Current results strongly support the conclusion that field growth occurs during the formation of the first galaxies on timescales shorter than their accretion timescales due to small-scale turbulent dynamos. The saturated small-scale dynamo maintains field strengths at only a few percent of equipartition with turbulence. This is in contradiction with the observed magnitude of turbulent fields, but may be reconciled by the further contribution to the turbulent field of the large-scale dynamo. The subsequent action of large-scale dynamos in differentially rotating discs produces field strengths observed in low redshift galaxies, where it reaches equipartition with the turbulence and has substantial power at large scales. The field structure resulting appears consistent with observations including Faraday rotation and polarisation from synchrotron and dust thermal emission. Major remaining challenges include scaling numerical models toward realistic scale separations and Prandtl and Reynolds numbers.

Generation of Near-Equipartition Magnetic Fields in Turbulent Collisionless Plasmas

The mechanisms that generate "seed" magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ_{B}≃0.4u_{rms}/L, where L/2π is the driving scale and u_{rms} is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k_{int}≃12π/L. Our results show that turbulence-induced by, e.g., the gravitational buildup of galaxies and galaxy clusters-can magnetize collisionless plasmas with large-scale near-equipartition fields.

Inefficient Magnetic-Field Amplification in Supersonic Laser-Plasma Turbulence

We report a laser-plasma experiment that was carried out at the LMJ-PETAL facility and realized the first magnetized, turbulent, supersonic (Ma_{turb}≈2.5) plasma with a large magnetic Reynolds number (Rm≈45) in the laboratory. Initial seed magnetic fields were amplified, but only moderately so, and did not become dynamically significant. A notable absence of magnetic energy at scales smaller than the outer scale of the turbulent cascade was also observed. Our results support the notion that moderately supersonic, low-magnetic-Prandtl-number plasma turbulence is inefficient at amplifying magnetic fields compared to its subsonic, incompressible counterpart.

3D simulations of oxygen shell burning with and without magnetic fields

We present a first 3D magnetohydrodynamic (MHD) simulation of convective oxygen and neon shell burning in a non-rotating [Formula: see text] star shortly before core collapse to study the generation of magnetic fields in supernova progenitors. We also run a purely hydrodynamic control simulation to gauge the impact of the magnetic fields on the convective flow and on convective boundary mixing. After about 17 convective turnover times, the magnetic field is approaching saturation levels in the oxygen shell with an average field strength of [Formula: see text], and does not reach kinetic equipartition. The field remains dominated by small-to-medium scales, and the dipole field strength at the base of the oxygen shell is only [Formula: see text]. The angle-averaged diagonal components of the Maxwell stress tensor mirror those of the Reynolds stress tensor, but are about one order of magnitude smaller. The shear flow at the oxygen-neon shell interface creates relatively strong fields parallel to the convective boundary, which noticeably inhibit the turbulent entrainment of neon into the oxygen shell. The reduced ingestion of neon lowers the nuclear energy generation rate in the oxygen shell and thereby slightly slows down the convective flow. Aside from this indirect effect, we find that magnetic fields do not appreciably alter the flow inside the oxygen shell. We discuss the implications of our results for the subsequent core-collapse supernova and stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects for a better understanding of magnetic fields in supernova progenitors.

Magnetic field amplification in accretion discs around the first stars: implications for the primordial IMF

Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal magneto-hydrodynamic simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called small-scale turbulent dynamo and the large-scale mean-field dynamo. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.

Efficient Highly Subsonic Turbulent Dynamo and Growth of Primordial Magnetic Fields

We present the first study on the amplification of magnetic fields by the turbulent dynamo in the highly subsonic regime, with Mach numbers ranging from 10^{-3} to 0.4. We find that for the lower Mach numbers the saturation efficiency of the dynamo (E_{mag}/E_{kin})_{sat} increases as the Mach number decreases. Even in the case when injection of energy is purely through longitudinal forcing modes, (E_{mag}/E_{kin})_{sat}≳10^{-2} at a Mach number of 10^{-3}. We apply our results to magnetic field amplification in the early Universe and predict that a turbulent dynamo can amplify primordial magnetic fields to ≳10^{-16}  G on scales up to 0.1 pc and ≳10^{-13}  G on scales up to 100 pc. This produces fields compatible with lower limits of the intergalactic magnetic field inferred from blazar γ-ray observations.

Synthesizing Observations and Theory to Understand Galactic Magnetic Fields: Progress and Challenges

Constraining dynamo theories of magnetic field origin by observation is indispensable but challenging, in part because the basic quantities measured by observers and predicted by modelers are different. We clarify these differences and sketch out ways to bridge the divide. Based on archival and previously unpublished data, we then compile various important properties of galactic magnetic fields for nearby spiral galaxies. We consistently compute strengths of total, ordered, and regular fields, pitch angles of ordered and regular fields, and we summarize the present knowledge on azimuthal modes, field parities, and the properties of non-axisymmetric spiral features called magnetic arms. We review related aspects of dynamo theory, with a focus on mean-field models and their predictions for large-scale magnetic fields in galactic discs and halos. Furthermore, we measure the velocity dispersion of H i gas in arm and inter-arm regions in three galaxies, M 51, M 74, and NGC 6946, since spiral modulation of the root-mean-square turbulent speed has been proposed as a driver of non-axisymmetry in large-scale dynamos. We find no evidence for such a modulation and place upper limits on its strength, helping to narrow down the list of mechanisms to explain magnetic arms. Successes and remaining challenges of dynamo models with respect to explaining observations are briefly summarized, and possible strategies are suggested. With new instruments like the Square Kilometre Array (SKA), large data sets of magnetic and non-magnetic properties from thousands of galaxies will become available, to be compared with theory.

From Primordial Seed Magnetic Fields to the Galactic Dynamo

The origin and maintenance of coherent magnetic fields in the Universe is reviewed with an emphasis on the possible challenges that arise in their theoretical understanding. We begin with the interesting possibility that magnetic fields originated at some level from the early universe. This could be during inflation, the electroweak, or the quark-hadron phase transitions. These mechanisms can give rise to fields which could be strong, but often with much smaller coherence scales than galactic scales. Their subsequent turbulent decay decreases their strength but increases their coherence. We then turn to astrophysical batteries which can generate seed magnetic fields. Here the coherence scale can be large, but the field strength is generally very small. These seed fields need to be further amplified and maintained by a dynamo to explain observed magnetic fields in galaxies. Basic ideas behind both small and large-scale turbulent dynamos are outlined. The small-scale dynamo may help to understand the first magnetization of young galaxies, while the large-scale dynamo is important for the generation of fields with scales larger than the stirring scale, as observed in nearby disk galaxies. The current theoretical challenges that turbulent dynamos encounter and their possible resolution are discussed.

Kazantsev dynamo in turbulent compressible flows

We consider the kinematic fluctuation dynamo problem in a flow that is random, white-in-time, with both solenoidal and potential components. This model is a generalization of the well-studied Kazantsev model. If both the solenoidal and potential parts have the same scaling exponent, then, as the compressibility of the flow increases, the growth rate decreases but remains positive. If the scaling exponents for the solenoidal and potential parts differ, in particular if they correspond to typical Kolmogorov and Burgers values, we again find that an increase in compressibility slows down the growth rate but does not turn it off. The slow down is, however, weaker and the critical magnetic Reynolds number is lower than when both the solenoidal and potential components display the Kolmogorov scaling. Intriguingly, we find that there exist cases, when the potential part is smoother than the solenoidal part, for which an increase in compressibility increases the growth rate. We also find that the critical value of the scaling exponent above which a dynamo is seen is unity irrespective of the compressibility. Finally, we realize that the dimension d = 3 is special, as for all other values of d the critical exponent is higher and depends on the compressibility.

Energy cascade and its locality in compressible magnetohydrodynamic turbulence

We investigate energy transfer across scales in three-dimensional compressible magnetohydrodynamic (MHD) turbulence, a model often used to study space and astrophysical plasmas. Analysis shows that kinetic and magnetic energies cascade conservatively from large to small scales in cases with varying degrees of compression. With more compression, energy fluxes due to pressure dilation and subscale mass flux are greater, but conversion between kinetic and magnetic energy by magnetic line stretching is less efficient. Energy transfer between the same fields is dominated by local contributions regardless of compressive effects. In contrast, the conversion between kinetic and internal energy by pressure dilation is dominated by the largest scale contributions. Energy conversion between the velocity and magnetic fields is weakly local.

Saturation of the turbulent dynamo

The origin of strong magnetic fields in the Universe can be explained by amplifying weak seed fields via turbulent motions on small spatial scales and subsequently transporting the magnetic energy to larger scales. This process is known as the turbulent dynamo and depends on the properties of turbulence, i.e., on the hydrodynamical Reynolds number and the compressibility of the gas, and on the magnetic diffusivity. While we know the growth rate of the magnetic energy in the linear regime, the saturation level, i.e., the ratio of magnetic energy to turbulent kinetic energy that can be reached, is not known from analytical calculations. In this paper we present a scale-dependent saturation model based on an effective turbulent resistivity which is determined by the turnover time scale of turbulent eddies and the magnetic energy density. The magnetic resistivity increases compared to the Spitzer value and the effective scale on which the magnetic energy spectrum is at its maximum moves to larger spatial scales. This process ends when the peak reaches a characteristic wave number k☆ which is determined by the critical magnetic Reynolds number. The saturation level of the dynamo also depends on the type of turbulence and differs for the limits of large and small magnetic Prandtl numbers Pm. With our model we find saturation levels between 43.8% and 1.3% for Pm≫1 and between 2.43% and 0.135% for Pm≪1, where the higher values refer to incompressible turbulence and the lower ones to highly compressible turbulence.

Small-scale dynamo at low magnetic Prandtl numbers

The present-day Universe is highly magnetized, even though the first magnetic seed fields were most probably extremely weak. To explain the growth of the magnetic field strength over many orders of magnitude, fast amplification processes need to operate. The most efficient mechanism known today is the small-scale dynamo, which converts turbulent kinetic energy into magnetic energy leading to an exponential growth of the magnetic field. The efficiency of the dynamo depends on the type of turbulence indicated by the slope of the turbulence spectrum v(ℓ)∝ℓ^{ϑ}, where v(ℓ) is the eddy velocity at a scale ℓ. We explore turbulent spectra ranging from incompressible Kolmogorov turbulence with ϑ=1/3 to highly compressible Burgers turbulence with ϑ=1/2. In this work, we analyze the properties of the small-scale dynamo for low magnetic Prandtl numbers Pm, which denotes the ratio of the magnetic Reynolds number, Rm, to the hydrodynamical one, Re. We solve the Kazantsev equation, which describes the evolution of the small-scale magnetic field, using the WKB approximation. In the limit of low magnetic Prandtl numbers, the growth rate is proportional to Rm^{(1-ϑ)/(1+ϑ)}. We furthermore discuss the critical magnetic Reynolds number Rm_{crit}, which is required for small-scale dynamo action. The value of Rm_{crit} is roughly 100 for Kolmogorov turbulence and 2700 for Burgers. Furthermore, we discuss that Rm_{crit} provides a stronger constraint in the limit of low Pm than it does for large Pm. We conclude that the small-scale dynamo can operate in the regime of low magnetic Prandtl numbers if the magnetic Reynolds number is large enough. Thus, the magnetic field amplification on small scales can take place in a broad range of physical environments and amplify week magnetic seed fields on short time scales.

Magnetic field amplification by small-scale dynamo action: dependence on turbulence models and Reynolds and Prandtl numbers

The small-scale dynamo is a process by which turbulent kinetic energy is converted into magnetic energy, and thus it is expected to depend crucially on the nature of the turbulence. In this paper, we present a model for the small-scale dynamo that takes into account the slope of the turbulent velocity spectrum v(ℓ)proportional ℓ([symbol see text])V}, where ℓ and v(ℓ) are the size of a turbulent fluctuation and the typical velocity on that scale. The time evolution of the fluctuation component of the magnetic field, i.e., the small-scale field, is described by the Kazantsev equation. We solve this linear differential equation for its eigenvalues with the quantum-mechanical WKB approximation. The validity of this method is estimated as a function of the magnetic Prandtl number Pm. We calculate the minimal magnetic Reynolds number for dynamo action, Rm_{crit}, using our model of the turbulent velocity correlation function. For Kolmogorov turbulence ([symbol see text] = 1/3), we find that the critical magnetic Reynolds number is Rm(crit) (K) ≈ 110 and for Burgers turbulence ([symbol see text] = 1/2) Rm(crit)(B) ≈ 2700. Furthermore, we derive that the growth rate of the small-scale magnetic field for a general type of turbulence is Γ proportional Re((1-[symbol see text])/(1+[symbol see text])) in the limit of infinite magnetic Prandtl number. For decreasing magnetic Prandtl number (down to Pm >/~ 10), the growth rate of the small-scale dynamo decreases. The details of this drop depend on the WKB approximation, which becomes invalid for a magnetic Prandtl number of about unity.