# Magnetic Field Amplification in Young Galaxies

###### Key Words.:

Dynamo – Magnetohydrodynamics (MHD) – Turbulence – Galaxies: magnetic fields – Galaxies: high-redshift^{1}

^{2}

## Abstract

Context:The Universe at present is highly magnetized, with fields of the order of a few G and coherence lengths larger than 10 kpc in typical galaxies like the Milky Way.

Aims:We propose that the magnetic field was amplified to this values already during the formation and the early evolution of galaxies. Turbulence in young galaxies is driven by accretion as well as by supernova (SN) explosions of the first generation of stars. The small-scale dynamo can convert the turbulent kinetic energy into magnetic energy and amplify very weak primordial seed fields on short timescales. Amplification takes place in two phases: in the kinematic phase the magnetic field grows exponentially, with the largest growth rate on the smallest non-resistive scale. In the following non-linear phase the magnetic energy is shifted towards larger scales until the dynamo saturates on the turbulent forcing scale.

Methods:To describe the amplification of the magnetic field quantitatively we model the microphysics in the interstellar medium (ISM) of young galaxies and determine the growth rate of the small-scale dynamo. We estimate the resulting saturation field strengths and dynamo timescales for two turbulent forcing mechanisms: accretion-driven turbulence and SN-driven turbulence. We compare them to the field strength that is reached, when only stellar magnetic fields are distributed by SN explosions.

Results:We find that the small-scale dynamo is much more efficient in magnetizing the ISM of young galaxies. In the case of accretion-driven turbulence a magnetic field strength of the order of is reached after a time of , while in SN-driven turbulence the dynamo saturates at field strengths of typically after only . This is considerably shorter than the Hubble time.

Conclusions:Our work can help to understand why present-day galaxies are highly magnetized.

## 1 Introduction

The present-day Universe is filled with magnetic fields. Observations show that galaxies (Beck et al. 1999; Beck 2011) and stars (Donati & Landstreet 2009; Reiners 2012) are strongly magnetized and there are also hints towards weak magnetic fields in the intergalactic medium (Kim et al. 1989; Kronberg 1994; Neronov et al. 2013). The origin of these strong fields remains an unsolved problem in astrophysics.

Local spiral galaxies have typical turbulent field components of G within the arms and bars, while a field of G is observed in the central starburst regions. Moreover, these fields appear to be coherent on scales larger than 10 kpc, which is the same order of magnitude as the size of the galaxy. The magnetic energy in the galactic interstellar medium is thus approximately in equipartition with the thermal energy and the energy in cosmic rays. The field in the interarm region is usually ordered and has a strength of the order of G (Beck 2011). Also dwarf irregular galaxies have magnetic fields, however they appear not to be ordered on large scales and have a lower strength of G (Chyży et al. 2011).

New observations indicate that even highly redshifted galaxies have magnetic field strengths comparable to present-day galaxies (Bernet et al. 2008). For instance the rotation measure, a quantity depending on the magnetic field along the line of sight, is constant up to redshifts of roughly 5 (Hammond et al. 2012). An important tool is moreover the FIR-radio correlation, which can be interpreted as a relation between the star formation rate and the synchrotron radiation of cosmic ray electrons in magnetic fields (Sargent et al. 2010; Bourne et al. 2011) and appears to be valid until (Murphy 2009). We note however, that one expects a breakdown at higher redshift as a result of inverse Compton scattering with cosmic microwave background photons (Schleicher & Beck 2013). Observations of the intergalactic medium provide further information on primordial seed fields. Detailed analysis of the CMB temperature bispectrum using data from the PLANCK satellite gives an upper limit of the magnetic field strength of the order of a few nG on the Mpc scale (Shiraishi et al. 2012). The increasing evidence for magnetic fields in highly redshifted galaxies and the intergalactic medium indicates an early generation of the magnetic fields.

Theoretically the first seed fields might already have been generated in the very early Universe during inflation leading to a field strength of - G on a scale of Mpc (Turner & Widrow 1988).
Another generation mechanism are first order phase transitions. Sigl et al. (1997) predict a field strength of from the electroweak phase transition and from the QCD phase transition on a scale of Mpc. The correlation length of the primordial magnetic seed fields has been shown to depend crucially on the initial properties of the field, e.g. on the amount of magnetic helicity (Banerjee & Jedamzik 2004). We can determine the typical strength of a statistical seed field by writing the magnetic energy as . Magnetic fields can also be generated as a result of the Biermann term in the generalized Ohm’s law, which takes into account the different behavior of the electron and ion fluid (Biermann 1950; Kulsrud & Zweibel 2008). A typical field strength resulting from this so-called “Biermann battery” is (Xu et al. 2008). Recently, Schlickeiser (2012) has shown that a turbulent magnetic field can be generated in plasma fluctuations within an unmagnetized non-relativistic medium. From this effect we would expect typical seed fields of a few G within the first galaxies. Although this exceeds the resulting field strengths of other generation mechanisms it still cannot explain the typical values in local galaxies. Thus, amplification processes need to take place.

A very efficient mechanism to amplify weak seed fields is the small-scale or turbulent dynamo, which converts kinetic energy from turbulence into magnetic energy by randomly stretching and twisting the field lines. The magnetic energy grows exponentially in the kinematic phase, while the growth rate is largest on the resistive scale (Kazantsev 1968; Subramanian 1997; Brandenburg & Subramanian 2005). There are two dimensionless parameters, which control the efficiency of this process (Schober et al. 2012b; Bovino et al. 2013): the hydrodynamic Reynolds number

(1) |

where is the turbulent velocity on the outer scale of the inertial range and is the viscosity, and the magnetic Reynolds number

(2) |

with being the magnetic resistivity. The ratio of the Reynolds numbers defines the magnetic Prandtl number

(3) |

Furthermore, the dynamo growth rate depends on the type of the turbulence ranging from incompressible Kolmogorov turbulence (Kolmogorov 1941) to highly compressible Burgers turbulence (Burgers 1948). Eventually the magnetic field is strong enough for back reactions to occur on the gas and the non-linear growth sets in (Schekochihin et al. 2002). In this phase the magnetic energy is transported towards larger scales. Now the evolution no longer depends on the Reynolds and Prandtl numbers, but still on the type of turbulence (Schleicher et al. 2013). The non-linear phase comes to an end, when the magnetic field reaches saturation on the turbulent forcing scale .

Turbulence is driven efficiently for the first time in the history of the Universe when dark matter halos become massive enough that gas begins to cool efficiently and flows into the potential wells of dark matter halos. This leads to the formation of the first generation of stars and the subsequent build-up of galaxies. The first stars form at redhifts between 20 and 15 within primordial minihalos, which have typical masses of more than solar masses () (see e.g. Abel et al. (2002); Bromm & Larson (2004); Clark et al. (2011)). In recent publications it was shown, numerically as well as semi-analytically, that during the formation of the first stars, dynamically important magnetic fields can be generated by the small-scale dynamo on short timescales (Sur et al. 2012; Turk et al. 2012; Schober et al. 2012a). According to the theory of hierarchical structure formation the first galaxies, also called protogalaxies, form at redshifts smaller than about 10 in massive dark matter halos with more than (Greif et al. 2008; Bromm et al. 2009). In young galaxies accretion as well as the penetration of supernovae (SN) shocks through the gas generate turbulence, which initiates small-scale dynamo action (Beck et al. 2012; Latif et al. 2013).

In this paper, we follow the evolution of the magnetic field in an initially weakly magnetized young galaxy. Because the true dynamical nature of the first galaxies is not well known, we adopt two simplified complementary models: a spherical galaxy as well as a disk-like system, both with constant density and temperature. We model microphysical processes, such as the diffusion of the kinematic and magnetic energy, in order to find the magnetohydrodynamical (MHD) quantities, which determine the growth rate of the small-scale dynamo. Turbulence can be generated by accretion flows into the center of the halo, for which we estimate the typical Reynolds numbers. Then we follow the evolution of the magnetic field strength in the kinematic and the non-linear phase, until saturation on the driving scale of the turbulence is reached. Also stellar feedback, in particular SN explosions, influences the evolution of the magnetic field. On the one hand supernovae distribute stellar magnetic fields in the interstellar medium (ISM) (Rees 1987), on the other hand they drive turbulence, which again leads to dynamo action (Balsara et al. 2004). We compare the resulting magnetic field strengths from both mechanisms with the field strength gained by an accretion-driven small-scale dynamo.

The outline of the paper is as follows: In section 2 we describe our model. We determine the values of viscosity and magnetic diffusivity in the interstellar medium and estimate the evolution of SN explosions. Driving mechanisms of turbulence are discussed in general. In the last part of this section we summarize the main points of a mathematical description of magnetic field amplification by the small-scale dynamo. The kinematic phase described by the so-called Kazantsev theory and a model for the non-linear growth phase are introduced. In section 3 we present our results for the evolution of the magnetic field in the different types of models. First we discuss the generation of turbulence by accretion and the resulting efficiencies of the dynamo, i.e. the saturation magnetic field strength and the time until saturation occurs. Second, we analyze the effect of stellar feedback. We compare the efficiency of distributing stellar magnetic fields by SN with the one of the SN-driven turbulent dynamo. We draw our conclusions in section 4.

## 2 Modeling Physical Processes in a Protogalaxy

### 2.1 General Aspects

The nature of young galaxies is still an active topic of research (see Bromm & Yoshida (2011) for a review). For our order of magnitude estimate of the magnetic field evolution we use a very simplified model, with the choice of parameters being motivated from numerical simulations (Greif et al. 2008; Bromm et al. 2009; Latif et al. 2013). We are interested in massive protogalactic objects at redshifts of roughly 10.

In our model we assume a mean particle density of

(4) |

and a temperature of

(5) |

The density as well as the temperature are, as first approximation, constant throughout the whole galaxy. For simplicity we take into account a gas that only consists of hydrogen, which is at the given values of and mostly ionized.

The mean shape of the primordial galaxies differs most probably from the one of present-day galaxies. Due to a significant amount of angular momentum the protogalaxies may form in a spherical way and develop a more disk-like structure at later stages. To account for the unknown typical shape, we model two extreme cases, a spherical and a disk-like galaxy, which have the same gas mass.

**Spherical Galaxy** In the case of a spherical protogalaxy we assume the radius to be

(6) |

As within this radius the density as well as the temperature are constant we find a total mass of the baryonic gas of

(7) |

**Disk-Like Galaxy** As our second fiducial model we use a galaxy with disk scale height of ten percent of the radius, i.e.

(8) |

With the condition that the gas mass of the disk needs to be the same as in the spherical case, the disk radius is

(9) |

### 2.2 Microphysics in the ISM

As the temperature in the primordial ISM is very high, we can assume the gas to be (at least partially) ionized. We thus need to deal with the full plasma equations, i.e. the continuity, the momentum and the energy equations for both the ions and the electrons. Closures of these equations were found by Braginskii (1965), who used the Chapman-Enskog scheme (Chapman et al. 1953). The closure is based on the assumption that the macroscopic scale of the plasma is large compared to the mean-free path

(10) |

or compared to the gyro-radii of the electrons and the ions

(11) |

Here is the distance of closest particle approach with being the elementary charge and the Boltzmann constant. The mass of the species is labeled , where the s stands for electrons (e) or protons (p), and is the speed of light. Further, we use here the thermal velocity and assume that the temperatures of the ions and electrons are equal (). In principle, the components of a plasma can have unequal temperatures as during plasma heating the different fluids are heated differently. However, after a certain time , an equilibrium will be reached. The electron-proton equilibrium time can be computed by (Spitzer 1956)

(12) |

where is the charge of species s, its initial temperature and the Coulomb logarithm is defined by

(13) | |||||

If we assume and to be extremely different, e.g. , the typical for our model is of the order of 440 yr. It will be shown later that this is way below the typical dynamo timescales, which can be up to many Myr. Thus, the electron and the proton temperature can be assumed to be equal in our calculation.

A comparison of the length scales (10) and (11) in our model can be found in figure 1. When the gyro-radius becomes smaller than the mean-free path, the magnetic field dominates the dynamics of the plasma, i.e. it becomes “magnetized”. In our model the electron fluid becomes magnetized at a magnetic field strength of roughly G, the ion fluid at G.

#### Viscosity

In the transition from an unmagnetized to a magnetized state, the plasma becomes anisotropic, i.e. certain physical quantities then depend on their relative orientation to the magnetic field direction.

In the unmagnetized case the kinematic viscosities for electrons and ions obtained from the Chapman-Enskog closure scheme are (Braginskii 1965)

(14) |

(15) |

with the collision times for electrons and ions

(16) | |||||

(17) |

In the presence of a strong magnetic field the viscosity becomes anisotropic and one has to distinguish between the viscosity along (parallel to) and the one perpendicular to the magnetic field lines. While the parallel viscosity stays the same as in the unmagnetized case (e.g. equations 14 and 15), the viscosity perpendicular to the field is given by (Simon 1955)

(18) | |||||

(19) |

with the gyro-frequencies

(20) | |||||

(21) |

Diffusion perpendicular to the magnetic field lines is also known as “Bohm diffusion”.

The different viscosities as a function of density are shown in figure 2. Note that the perpendicular viscosity becomes only valid when the plasma is magnetized, i.e. when the gyro-radius becomes smaller than the mean-free path. According to figure 1 this is the case above a magnetic field strength of G for the electrons and G for the ions. Thus, the most important part of the viscosity is the parallel one and we will ignore the perpendicular part, which decreases proportional the , from now on. Furthermore, the viscosity of the ions exceeds the electron viscosity by roughly two orders of magnitude. This is caused by the fact that the ions carry the largest part of the momentum. In total, the parallel viscosity of the ions is the crucial quantity and we will refer from now on to

(22) |

#### Magnetic Diffusivity

For the parallel conductivity the closure scheme yields (Spitzer 1956)

(23) |

and for the conductivity perpendicular to the magnetic field

(24) |

The conductivity perpendicular to the magnetic field lines is, contrary to the case of viscosity, no function of the magnetic field strength. The difference between the parallel and the perpendicular component of the conductivity is just approximately a factor of two. Usually, is used to determine the magnetic diffusivity of a plasma. We thus find

(25) |

which is also known as “Spitzer resistivity”.

#### Magnetic Prandtl Number

With these values of viscosity and resistivity the magnetic Prandtl number (see equation 3) is

(26) |

### 2.3 Turbulence

#### Generation of Turbulent Motions by Accretion

Structure formation is always associated with accretion. In order to build up the first stars and galaxies, gas flows into the potential wells of dark matter halos, where it gets compressed and cools. The potential energy released during that process in parts gets converted into turbulent kinetic energy (Klessen & Hennebelle 2010). Simulations by Greif et al. (2008) of atomic cooling halos show that two types of accretion occur: in the so-called “hot accretion” mode gas is accreted directly from the intergalactic medium, while in the “cold accretion” mode gas is cooled down and flows into the central regions of the halo at high velocities (Dekel et al. 2009; Nelson et al. 2013).

#### Generation of Turbulent Motions by Supernova Explosions

Once stars have formed, their feedback strongly influences the ISM in galaxies in terms of ionizing radiation and at later stages by SN explosions, which are especially important for the generation of turbulence.

In order to calculate the corresponding energy input, we need to estimate the rate of SN explosions. The star formation rate () is proportional to the mass density over the free-fall time (Mac Low & Klessen 2004; McKee & Ostriker 2007):

(27) |

From the we can estimate the supernova rate (). For this we divide the by the typical mass of a star that results in a SN (). As not all the gas goes into stars and not all the stars are massive enough to end in a SN we introduce an efficiency factor :

(28) |

The number of supernovae within the whole galaxy with a volume and a time interval is then given by

(29) |

where we assume that the stays constant over time. In general the SN rate is expected to change with time, however modeling this time dependency goes beyond the scope of this work.

### 2.4 Turbulent Magnetic Field Amplification

#### Kinematic Small-Scale Dynamo

The induction equation,

(30) |

describes the time evolution of a magnetic field B, where v is the velocity and the magnetic diffusivity (25). An arbitrary magnetic field can, in general, be separated into a mean component and a fluctuating component with

(31) |

Substituting (31) into the induction equation leads to two equations: an equation for the large-scale field evolution and the Kazantsev equation (Kazantsev 1968; Brandenburg & Subramanian 2005), which describes the small-scale evolution of the field.

The derivation of the Kazantsev equation is based on the assumption that the fluctuations of the magnetic field as well as the fluctuations of the velocity field are homogeneous and isotropic even if the mean fields are not isotropic. Furthermore, the fluctuations are assumed to be Gaussian with a zero mean and the velocity fluctuations are thought to be -correlated in time. For simplicity, any helicity of the magnetic field is neglected. With these assumptions the Kazantsev equation is (Kazantsev 1968)

(32) |

The eigenfunctions of this equation are related to the longitudinal correlation function of the magnetic fluctuations by . We call the growth rate of the small-scale magnetic field. The function is the magnetic diffusion coefficient, which contains besides the magnetic diffusivity also a scale-dependent turbulent diffusivity. is called the “potential” of the Kazantsev equation. Both and depend only on the correlation function of the turbulent velocity field and the magnetic diffusivity (Subramanian 1997; Schober et al. 2012b).

The correlation function of the turbulent velocity field in turn depends on the different types of turbulence, which can be distinguished by the slope of the velocity spectrum in the inertial range, where

(33) |

Here is the velocity of the fluctuations on the scale . The range of goes from incompressible Kolmogorov turbulence with (Kolmogorov 1941) to highly compressible Burgers turbulence with (Burgers 1948). The gas motions during structure formation have high Mach numbers and thus the gas gets strongly compressed within shocks. Observations within present-day molecular clouds by Larson (1981) show that the slope of the turbulence spectrum is and thus deviates from Kolmogorov turbulence. However, other studies (Solomon et al. 1987; Ossenkopf & Mac Low 2002; Heyer & Brunt 2004) find a slope of roughly 0.5, whereas Roman-Duval et al. (2011) show that the variance of is very large. For our fiducial model we choose a value of , which lies in between the extremes.

With a model for the turbulent correlation function, the Kazantsev equation (32) can be solved with the WKB-approximation for very large and low Pm. This method is named after Wentzel, Kramers and Brillouin and is used to find approximative solutions for Schrödinger-type differential equations. In our model we are in the limit of the very high Pm, where Schober et al. (2012b) find the growth rate

(34) |

Here is the typical velocity on the largest scale of the turbulent eddies of size . By solving the Kazantsev equation (32) numerically, Bovino et al. (2013) have recently confirmed that equation (34) describes the growth rate of the dynamo in the limit large Pm. For our fiducial model with the growth rate thus scales with .

#### Non-Linear Small-Scale Dynamo

As soon as the magnetic energy is comparable to the kinetic energy of the turbulence on the viscous scale the exponential growth comes to an end. We label this point in time . The dynamo is then saturated on the viscous scale and the non-linear growth begins. In this phase the magnetic energy on the scale of fastest amplification evolves as (Schekochihin et al. 2002)

(35) |

with the non-linear growth rate

(36) |

and

(37) |

Schleicher et al. (2013) find that further evaluation of equation (35) yields

(38) |

Thus, in the case of Kolmogorov turbulence with the magnetic energy grows linear in time, while it grows quadratically in case of Burgers turbulence with . In our fiducial model, where we assume , we find on .

In the non-linear phase the dynamo process shifts the magnetic energy to larger scales with the peak scale evolving as

(39) |

From the peak scale to larger scales we assume the spectrum to drop off with the Kazantsev slope. By this we can determine the magnetic field on the forcing scale at each point in time as

(40) |

The non-linear growth phase comes to an end, when saturation on the turbulent forcing scale is achieved. Now the spectrum of the magnetic energy density scales as the one of the kinetic energy density.

#### Saturation Magnetic Field Strength from Dynamo Amplification

A turbulent dynamo can amplify magnetic fields at most to equipartition with the turbulent kinetic energy. However, high-resolution simulations by Federrath et al. (2011) show that only a certain fraction of the turbulent kinetic energy can be transformed into magnetic energy. This fraction depends on the type of forcing as well as on the Mach number . We show for solenoidal and compressive forcing of turbulence in figure 3. Note, that the efficiency of the small-scale dynamo in case of compressive forcing peaks at a Mach number of 1, i.e. at the transition from the subsonic to the supersonic regime. At this point shocks appear, which generate solenoidal motions that are more efficient for dynamo amplification. At larger Mach numbers the efficiency decreases again and appears to become constant.

According to Federrath et al. (2010) solenoidal forcing leads to a slope of the turbulence spectrum of , while compressive forcing results in . For our fiducial model we choose the saturation efficiency of solenoidal driven turbulence, as we assume a spectrum with .

The resulting saturation magnetic field strength on the forcing scale is

(41) |

where is again the velocity at the forcing scale. If we scale down the turbulent velocity to the viscous scale by

(42) |

the saturation magnetic field strength on the viscous scale is

(43) |

#### Evolution of a Magnetic Field Amplified by the Small-Scale Dynamo

Summarizing the results of this section gives for the magnetic field evolution on the viscous scale

(44) |

i.e. it grows exponentially with rate (34) until saturation on the viscous scale is reached at the time .

The field on the turbulent forcing scale evolves as

(45) |

Until the time the field grows exponentially in the kinematic phase. For the dynamo is in the non-linear phase, in which the peak of the magnetic spectrum, which is given by equation (39), is shifted towards larger scales. The dynamo is saturated on all scales of the turbulent inertial range including the driving scale for times .

The dynamo amplification of a weak magnetic seed field of G is shown in figure 4. We choose here a forcing scale of pc, which is the radius of the spherical halo considered here, and three different turbulent velocities: , and . The microphysical quantities are taken from the calculations in the previous sections. In the figure the dashed lines represent the magnetic field strength on the viscous scale, the solid lines the one on the forcing scale.

## 3 Magnetic Field Evolution in a Protogalaxy

### 3.1 Magnetic Fields from an Accretion-Driven Small-scale Dynamo

#### Forcing Turbulence by Accretion

**Accretion in a Spherical Galaxy** During the formation of the primordial halo turbulence is generated by accretion (Birnboim & Dekel 2003; Semelin & Combes 2005; Wise et al. 2008; Vogelsberger et al. 2013). Simulations show that accretion flows have high Mach numbers with respect to the cold gas even in the central regions of the halo (Greif et al. 2008). The characteristic forcing scale in case of a spherical halo is the radius, where the accretion flow comes to a halt:

(46) |

Latif et al. (2013) show in their simulation of a nearly isothermal protogalaxy that the Mach number in such environment is roughly 2. Thus, the typical turbulent velocities from accretion are of the order of

(47) |

where is the sound speed and we use an adiabatic index of 5/3. Further, we assume here that only a certain fraction of the kinetic energy of the accretion flows goes into turbulence, with typically depending on the density contrast between the accretion flows and the halo (Klessen & Hennebelle 2010). Simulations (Latif et al. 2013) indicated that about five percent of the kinetic energy are in turbulent motions, i.e. .

The resulting turbulent length scales, velocities and Reynolds numbers for a spherical galaxy are given in table 1.

**Accretion in a Disk-Like Galaxy** In the case of a disk-like galaxy we adopt the typical forcing scale of the turbulence by accretion flows to be the scale height

(48) |

We further estimate the typical velocity for accretion flows to be of the order of the Kepler velocity in a disk

(49) |

If a percentage of the kinetic energy goes into turbulence, the resulting turbulent velocity is given by

(50) |

Typical values of the length scale, the velocity scales and the Reynolds numbers with a value of are summarized in table 1.

Accretion-Driven Dynamo | SN-Driven Dynamo | |||
---|---|---|---|---|

spherical galaxy | disk-shaped galaxy | spherical galaxy | disk-shaped galaxy | |

#### Accretion-Driven Small-Scale Dynamo

Based on the discussion of the strength of magnetic seed fields in the introduction, we assume the initial magnetic field strength on the viscous scale to be

(51) |

This is a rather conservative estimate.

The small-scale dynamo amplifies this seed field as soon as sufficient turbulence has evolved. The typical growth rates in the kinematic phase are summarized in table 1. We find 150 Myr for the case of a spherical galaxy and 1400 Myr for a disk. A fraction of the magnetic energy can be dissipated again by Ohmic diffusion. The dissipation rate on the viscous scale is given by

(52) |

In our model is of the order of and thus can be neglected compared the growth rate of the magnetic field.

With these relatively large growth rates, the small-scale dynamo amplification works on very short timescales. We find that in a spherical galaxy a magnetic field of G and be reached on a scale of pc after 270 Myr. In a disk the saturation field strength is larger by a factor of more than 2. However, the field is only on a scale of 240 pc, but it is saturated after already 24 Myr.

The efficiency of the small-scale dynamo, i.e. the saturation magnetic field strength ( or ) that can be achieved and the time on which saturation occurs ( or ), depends strongly on the amount of turbulent kinetic energy, controlled by the parameter (see equations 47 and 50). In our fiducial model we use , however this is a rough assumption. We test how the dynamo efficiency changes when varying in figure 5.

In the upper panel of figure 5 we show the dependency of the viscous scale and the forcing scale on . Of course is not effected by , while , which is a function of the Reynolds number and thus of the turbulent velocity, decreases with increasing . The time until saturation of the dynamo, which is shown in the middle panel of figure 5, also decreases with increasing . This is a natural consequence of the larger amount of turbulent kinetic energy. In the same way the plot in the lower panel can be understood: the more turbulent energy, i.e. the higher , the higher is the saturation field strength. The magnetic field strength on the forcing scale increases as

(53) |

### 3.2 Magnetic Fields from Stellar Feedback

#### Distributing Stellar Magnetic Fields by Supernovae

A natural source for magnetic fields in the ISM of galaxies are stellar magnetic fields that get distributed over large volumes by SN explosions. Schober et al. (2012a) have shown that the small-scale dynamo can produce strong magnetic fields during primordial star formation. Hints to dynamical important magnetic fields during the formation of the first stars come also from high-resolution numerical simulations (Federrath et al. 2011; Turk et al. 2012; Sur et al. 2012) and further semi-analytical calculations (Schleicher et al. 2010). Thus, we expect the first and second generations of stars to be magnetized.

**Properties of Supernova Candidates** We assume that a typical star that ends in a supernova has a mass of

(54) |

and a radius of

(55) |

with the solar mass and radius .

It is very difficult to estimate the magnetic energy in a typical population III star, as there is not much theoretical work on that topic so far. In principle, one could assume that a certain percentage of the total energy of the SN energy is within the magnetic field. If the magnetic energy equals e.g. 0.001 , the stellar magnetic field would have a very high value of .

Here, however, we use as a crude estimate for the magnetic field of population III stars based on observations of present-day massive stars. In most high-mass stars no magnetic fields are detected, there are few percent of stars with an enhanced magnetic field (see Donati & Landstreet (2009)). These so-called “peculiar A or B” stars have a typical dipole field strength of

(56) |

We take this value as an upper limit of magnetic fields in primordial stars, but test also lower stellar field strengths in the following.

##### Evolution of a Supernova Remnant

Stars with masses above 8 solar are expected to explode as a core-collapse supernova, introducing additional turbulent energy into the ISM (Choudhuri 1998; Padmanabhan 2001). Initially the shock front of a SN expands freely, i.e. the pressure of the surrounding ISM is negligible. The shock velocity can then be determined by

(57) |

where is the energy of a SN (neglecting the energy loss by neutrinos) and the ejected mass. The shock radius as a function of time is thus

(58) |

The free expansion phase ends, when the accumulated mass of the ISM in front of the shock is of order of . This happens at the so-called sweep-up radius defined by

(59) |

with being the mean density of the ISM. The shock front reaches at a time

(60) |

which is in our model of the order of 100 yr. For the expansion of the supernova remnant is driven adiabatically by thermal pressure, which is known as the Sedov-Taylor phase (Sedov 1946; Taylor 1950; Sedov 1959). We can estimate the radius of the shock in this case with

(61) |

where the indicates a time derivative and the pressure is given by

(62) |

with for adiabatic expansion. We can solve equation (61) with a simple power-law ansatz and find

(63) |

Thus, the evolution of the SN remnant can be described by (Choudhuri 1998)

(64) |

We assume now that the energy released in a SN explosion is and that about 10 percent of the mass of the progenitor star is ejected, i.e. . In our model the SN remnants evolve as described in equation (64) and shown in figure 6 until they collide. At later stages of shock evolution other energy loss mechanisms become dominant. The electrons lose their energy by ionization, bremsstrahlung, synchrotron emission and inverse Compton scattering. The latter is the most important energy loss channel at high redshifts as here the density of the CMB photons is considerably larger (Schleicher & Beck 2013).

If the SNe are distributed homogeneously in the protogalaxy, each SN shell has a mean maximum radius at the first collision of

(65) |

where the radius of the galaxy is given in equations (6) and (9) and the exponent depends on the geometry of the galaxy. In case of a spherical halo , in case of a thin disk . The maximum expansion radius of the SN shock is shown in figure 6 for the spherical case.

By comparing (64) to (65) we find the typical time scale for SN collisions . At that point the SN bubbles fill approximately the whole galaxy. In the spherical case we find

(66) |

which we take as the typical timescale for SN collisions. Further, we use as the typical length scale of SN shocks.

**Magnetic Field Evolution** If now all the stellar magnetic energy is distributed into the volume available by the SN explosion and no significant magnetic energy is left in the stellar remnant, the resulting magnetic field strength in the ISM after the first SN generation is

(67) |

Here, we assumed a spherical shape of the galaxy and flux freezing. All the following stars will produce roughly the same amount of magnetic energy, that is then distributed in the ISM by SN. Thus, the time evolution of the stellar magnetic fields in galaxies can be approximated by

(68) |

The values of and depend obviously on the SN rate, which is determined by the parameter as defined in (28). We obtain for the spherical case:

(69) | |||

(70) |

leading to a dependency of the magnetic field distributed by SN on the efficiency of the SN rate of

(71) |

In figure 7 we show the evolution of the distributed magnetic fields for different mean magnetic fields of the stars ( to G) and for our fiducial case of . Note that the case of G is an upper limit of magnetic fields in massive stars. We assume the magnetic fields of the first stars to be considerably lower.

The distribution of stellar magnetic fields by SNe explosions thus seems to be not important compared to the dynamo amplification in the ISM. However after a sufficient time the SN could contribute to the magnetic energy in the ISM. If the equipartition field strength is roughly G, the time after which SN distribution becomes important is

(72) |

For our fiducial model we find that this time is about Myr in case of typical stellar field strengths of G. Observations of present-day massive stars indicate that only a few percent have these high field strengths. We thus also consider the more likely case of lower mean stellar fields. For a mean strength of G we find that a micro-Gauss ISM field is only reached after Myr and for a mean strength of G after Myr. Thus, the typical timescales of distribution of stellar magnetic fields by supernovae exceed the age of the Universe by many orders of magnitude and this process cannot be an important contribution for the fields in the ISM, unless the first stars were much stronger magnetized than the present-day stars.

In the case of a flat disk-shaped galaxy, where we assume the parameter in equation (65) to be 1/2, the distribution of stellar magnetic fields proceeds marginally faster. Here the typical time until a field strength of G in the ISM is reached is roughly a factor of 10 shorter.

In reality the evolution of magnetic fields in SN shock fronts is of course more complicated. In addition to simple flux freezing further amplification processes can take place. Miranda et al. (1998) argue that in a multiple explosion scenario of structure formation (Ostriker & Cowie 1981; Miranda & Opher 1997) magnetic seed fields of the order of G can be produced on galactic scales. In their model a Biermann battery is operating in the shock of SN explosions of the first stars as here non-parallel gradients of temperature and density can be established. Recently also Beck et al. (2013) have analyzed the magnetic field evolution in protogalaxies based on SN explosions with the cosmological N-body code GADGET. They find that a combination of SNe and subsequent magnetic field amplification leads to magnetic field strengths of the order of a few G, which is comparable to our results, and that the strength of seed field is coupled to the star formation process.

#### Dynamo Amplification Driven by SN Turbulence

**SN-Driven Dynamo in a Spherical Galaxy** In section 2.3 we discussed the generation of SN turbulence based on numerical simulations. Now we estimate the typical forcing scale and the fluctuation velocity on that scale in order to determine the Reynolds number (1) and the resulting growth rate of the kinematic small-scale dynamo (34). For that we assume that the turbulence driving in the galaxy is in equilibrium.

Then the turbulent pressure, which is roughly

(73) |

balances the hydrostatic pressure determined by

(74) |

The gravitational acceleration in the spherical case is roughly . Solving equation (74) and setting it equal to equation (73) yields the turbulent velocity . We find in the spherical case:

(75) |

The forcing length scale can be estimated by comparing the energy input rate with the dissipation rate:

(76) |

where the dissipation timescale is

(77) |

and is the supernova rate (28). Thus, we find the typical forcing scale of SN-driven turbulence

(78) |

As in the case of the accretion-driven small-scale dynamo we start with an initial magnetic field strength on the viscous scale of

(79) |

The turbulence driven by SN shocks makes dynamo action possible, which leads to rapid amplification of the seed field according to equations (44) and (45). For the case of a spherical galaxy we find that the growth rate in the kinematic amplification phase is . After a time of the saturation field strength of is reached on the forcing scale. The characteristic quantities of our fiducial models for the SN-driven dynamo are summarized in the right part of table 1.

As in case of accretion turbulence, the efficiency of the small-scale dynamo is sensible to the amount of kinetic energy that goes into turbulence . Moreover, when modeling the scale of turbulence forcing we add another uncertainty namely the , which includes the efficiency parameter (see equation 28). In our fiducial model we choose , but there could easily be a variation of a factor 10. The dependency of the quantities most important for the dynamo amplification on and in case of a spherical halo is the following:

(80) | |||||

(81) | |||||

(82) | |||||

(83) | |||||

(84) |

We show the dependency of the length scales, the time until saturation and the saturation magnetic field strength on and for different values of in figure 8.

**SN-Driven Dynamo in a Disk-Like Galaxy** We perform the same analysis for the disk case. Here the gravitational acceleration becomes independent of the radius for a thin disk, i.e. . In that approximation we find , which leads to a turbulent velocity of

(85) |

The forcing scale can be determined by equation (78). In case of a disk-shaped galaxy we find that the typical forcing scale

(86) |

We find that the kinematic growth rate in our fiducial model is . The time until saturation on the forcing scale is then only 3.8 Myr and the saturation field strength is .

In case of a disk-shaped galaxy all the quantities (79) to (84) are independent of and . For comparison with the spherical galaxy we show them however also in figure 8.

## 4 Conclusions

In this paper we model the evolution of the (turbulent) magnetic field in a young galaxy. We find that weak magnetic seed fields get amplified very efficiently by the small-scale dynamo (see table 1), which is driven by turbulence from accretion and from supernova (SN) explosions. Dynamo theory predicts that the magnetic field is amplified in two phases: in the kinematic phase the field grows exponentially until the dynamo is saturated on the viscous scale. Then the non-linear phase begins, where the magnetic energy is shifted towards larger scales until saturation on the turbulent forcing scale occurs.

For our fiducial models of a young galaxy we use a fixed particle density of and a temperature of . We concentrate on two different geometries: a spherical and a disk-shaped galaxy (see section 2.1). We determine the viscosity of the plasma, which becomes anisotropic when the plasma becomes magnetized, and the magnetic diffusivity. Turbulence is generated by accretion flows onto the galactic core and also by SN shocks. By estimating typical driving scales and velocities we can determine the hydrodynamic and the magnetic Reynolds number. The magnetic field evolution depends strongly on the type of turbulence, which we assume to be a mixture of solenoidal and compressive modes.

For our fiducial model we find that the dynamo saturates on the largest scale in accretion-driven turbulence after a time of roughly in case of a spherical galaxy and after in case of a disk. Turbulence generated by SN shocks can amplify the magnetic field on shorter timescales, with saturation reached after in a spherical galaxy and in a disk. The dynamo timescale is thus comparable to the free-fall time . The age of the Universe at the onset of galaxy formation, i.e. at a redshift of 10, is roughly 470 Myr, which is larger than the dynamo timescales by factor of 2 to 120 for our four fiducial models. In the models with the longest amplification times our assumption of constant accretion and supernova rates may thus not be very precise. Nevertheless, these models provide an order of magnitude estimate of the resulting magnetic strength. In case of a disk-like galaxy we can compare the dynamo timescales further to the typical time of one rotation, which turns out to be when using the Kepler velocity (49). Thus, we can expect that the small-scale dynamo is saturated within less then one orbital time, the turbulent magnetic field gets ordered and an dynamo, i.e. a galactic large-scale dynamo, sets in. The role of rotation has been analyzed for example by Kotarba et al. (2009) and Kotarba et al. (2011) in numerical simulations and is extremely important for understanding the present-day large-scale structure of galactic magnetic fields.

The magnetic field strengths predicted by our fiducial models are very high with values between and in the accretion-driven case and between and in the SN-driven case for a spherical galaxy and a disk, respectively. These field strengths are comparable with the ones observed in the local Universe, where the typical turbulent field component in present-day disk galaxies is in spiral arms and bars and up to in the central starburst regions (Beck 2011). New radio observations detect also magnetic fields in dwarf galaxies. Their field strengths, which seem to be correlated with the SFR, are typically a factor of roughly three lower compared to the one in spiral galaxies (Chyży et al. 2011).

Our calculations suggest that the turbulent magnetic field of a galaxy was very high already at high redshifts. An observational confirmation of this result is very complicated. A hint towards an early generation of the turbulent magnetic field in galaxies comes from Hammond et al. (2012). They analyze the rotation measure of a huge catalog of extragalactic radio sources as a function of redshift and find that it is constant up to redshifts of 5.3, which is the maximum redshift in their dataset. A very powerful tool provides moreover the far-infrared (FIR) - radio correlation, which relates the star formation rate to the synchrotron loss of cosmic ray electrons. It is observed to be constant up to redshifts of roughly 2 (Sargent et al. 2010; Bourne et al. 2011), but is expected to break down at a higher redshift, which depends on the star formation rate and the evolution of typical ISM densities (Schleicher & Beck 2013). With new instruments like SKA and LOFAR our knowledge about the evolution of the cosmic magnetic fields will increase.

Besides our fiducial models we analyze the effect of changing the amount of kinetic energy that goes into turbulence and find that the dynamo is more efficient with increasing turbulent energy, which is intuitively clear. Furthermore we determine the small-scale dynamo evolution for a varying supernova rate (), which is important for estimating the driving scale of SN turbulence in the case of a spherical core. As expected the time until saturation increases with increasing . However, the typical largest scale of the magnetic field decreases with the .

We further estimate the effect of magnetic field enrichment in galaxies by distributing stellar fields by SN explosions. As an estimate of the magnetic energy in the first stars is very hard, we determine the expected magnetic field evolution in the ISM for three different cases. An upper limit of magnetic field strengths of the primordial stars is G, which is a value observed in a the few percent of present-day massive stars that are magnetized. Distributing these mean stellar fields by SNe in the ISM, we find that a ISM field strength of G is reached after Myr, which is already longer than the Hubble time. Thus, the dynamo increases the magnetic field strength much faster.

With our model we have shown that the small-scale dynamo can amplify weak magnetic seed fields in the ISM of early galaxies on relatively short time scales compared to other evolutionary timescales. This leads to the build-up of strong magnetic fields already at very early phases of (proto)galactic evolution. Theoretical models of galaxy evolution describe a collapse of a spherical object to a disk. Comparison of the gravitational energy, , with the magnetic energy at dynamo saturation, , shows that the field is not strong enough to prevent to collapse. Still the magnetic field produced by the small-scale dynamo has potentially strong impact on ISM dynamics and subsequent star formation.

###### Acknowledgements.

We thank the anonymous referee for useful comments on our manuscript. We acknowledge funding through the Deutsche Forschungsgemeinschaft (DFG) in the Schwerpunktprogramm SPP 1573 “Physics of the Interstellar Medium” under grant KL 1358/14-1 and SCHL 1964/1-1. Moreover, we thank for financial support by the Baden-Württemberg-Stiftung via contract research (grant P-LS-SPII/18) in their program “Internationale Spitzenforschung II” as well as the DFG via the SFB 881 “The Milky Way System” in the sub-projects B1 and B2. J. S. acknowledges the support by IMPRS HD. D. R. G. S. thanks for funding via the SFB 963/1 (project A12) on “Astrophysical flow instabilities and turbulence”.### Footnotes

- institutetext: Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, Albert-Ueberle-Strasse 2, D-69120 Heidelberg, Germany
- institutetext: Georg-August-Universität Göttingen, Institut für Astrophysik, Friedrich-Hund-Platz 1, D-37077 Göttingen, Germany

### References

- Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93
- Balsara, D. S., Kim, J., Mac Low, M.-M., & Mathews, G. J. 2004, ApJ, 617, 339
- Banerjee, R. & Jedamzik, K. 2004, Phys. Rev. D, 70, 123003
- Beck, A. M., Dolag, K., Lesch, H., & Kronberg, P. P. 2013, MNRAS
- Beck, A. M., Lesch, H., Dolag, K., et al. 2012, MNRAS, 422, 2152
- Beck, R. 2011, Space Sci. Rev., 135
- Beck, R., Ehle, M., Shoutenkov, V., Shukurov, A., & Sokoloff, D. 1999, Nature, 397, 324
- Bernet, M. L., Miniati, F., Lilly, S. J., Kronberg, P. P., & Dessauges-Zavadsky, M. 2008, Nature, 454, 302
- Biermann, L. 1950, Zeitschrift Naturforschung Teil A, 5, 65
- Birnboim, Y. & Dekel, A. 2003, MNRAS, 345, 349
- Bourne, N., Dunne, L., Ivison, R. J., et al. 2011, MNRAS, 410, 1155
- Bovino, S., Schleicher, D. R. G., & Schober, J. 2013, New Journal of Physics, 15, 013055
- Braginskii, S. I. 1965, Reviews of Plasma Physics, 1, 205
- Brandenburg, A. & Subramanian, K. 2005, Phys. Rep, 417, 1
- Bromm, V. & Larson, R. B. 2004, ARA&A, 42, 79
- Bromm, V. & Yoshida, N. 2011, ARA&A, 49, 373
- Bromm, V., Yoshida, N., Hernquist, L., & McKee, C. F. 2009, Nature, 459, 49
- Burgers, J. 1948, Advances in Applied Mechanics, Vol. 1, A Mathematical Model Illustrating the Theory of Turbulence (Elsevier), 171
- Chapman, S., Cowling, T., & Társulat, M. T. 1953, The Mathematical Theory of Non-uniform Gases: An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion in Gases. Prepared in Co-operation with D. Burnet (University Press)
- Choudhuri, A. R. 1998, The Physics of Fluids and Plasmas: An Introduction for Astrophysicists (Cambrigde University Press)
- Chyży, K. T., Weżgowiec, M., Beck, R., & Bomans, D. J. 2011, A&A, 529, A94
- Clark, P. C., Glover, S. C. O., Smith, R. J., et al. 2011, Science, 331, 1040
- Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451
- Donati, J.-F. & Landstreet, J. D. 2009, ARA&A, 47, 333
- Federrath, C., Chabrier, G., Schober, J., et al. 2011, Phys. Rev. Lett., 107, 114504
- Federrath, C., Roman-Duval, J., Klessen, R. S., Schmidt, W., & Mac Low, M.-M. 2010, A&A, 512, A81
- Federrath, C., Sur, S., Schleicher, D. R. G., Banerjee, R., & Klessen, R. S. 2011, ApJ, 731, 62
- Greif, T. H., Johnson, J. L., Klessen, R. S., & Bromm, V. 2008, MNRAS, 387, 1021
- Hammond, A. M., Robishaw, T., & Gaensler, B. M. 2012, ArXiv e-prints
- Heyer, M. H. & Brunt, C. M. 2004, ApJ, 615, L45
- Kazantsev, A. P. 1968, Soviet Journal of Experimental and Theoretical Physics, 26, 1031
- Kim, K.-T., Kronberg, P. P., Giovannini, G., & Venturi, T. 1989, Nature, 341, 720
- Klessen, R. S. & Hennebelle, P. 2010, A&A, 520, A17+
- Kolmogorov, A. 1941, Akademiia Nauk SSSR Doklady, 30, 301
- Kotarba, H., Lesch, H., Dolag, K., et al. 2011, MNRAS, 415, 3189
- Kotarba, H., Lesch, H., Dolag, K., et al. 2009, MNRAS, 397, 733
- Kronberg, P. P. 1994, Reports on Progress in Physics, 57, 325
- Kulsrud, R. M. & Zweibel, E. G. 2008, Reports on Progress in Physics, 71, 046901
- Larson, R. B. 1981, MNRAS, 194, 809
- Latif, M. A., Schleicher, D. R. G., Schmidt, W., & Niemeyer, J. 2013, MNRAS
- Mac Low, M.-M. & Klessen, R. S. 2004, Rev. Mod. Phys., 76, 125
- McKee, C. F. & Ostriker, E. C. 2007, ARA&A, 45, 565
- Miranda, O. D., Opher, M., & Opher, R. 1998, MNRAS, 301, 547
- Miranda, O. D. & Opher, R. 1997, ApJ, 482, 573
- Murphy, E. J. 2009, ApJ, 706, 482
- Nelson, D., Vogelsberger, M., Genel, S., et al. 2013, MNRAS, 429, 3353
- Neronov, A., Semikoz, D., & Banafsheh, M. 2013, ArXiv e-prints
- Ossenkopf, V. & Mac Low, M.-M. 2002, A&A, 390, 307
- Ostriker, J. P. & Cowie, L. L. 1981, ApJ, 243, L127
- Padmanabhan, T. 2001, Theoretical Astrophysics - Volume 2, Stars and Stellar Systems
- Rees, M. J. 1987, Q. J. R. Astron. Soc., 28, 197
- Reiners, A. 2012, Living Reviews in Solar Physics, 9, 1
- Roman-Duval, J., Federrath, C., Brunt, C., et al. 2011, ApJ, 740, 120
- Sargent, M. T., Schinnerer, E., Murphy, E., et al. 2010, ApJ, 714, L190
- Schekochihin, A. A., Cowley, S. C., Hammett, G. W., Maron, J. L., & McWilliams, J. C. 2002, New Journal of Physics, 4, 84
- Schleicher, D. R. G., Banerjee, R., Sur, S., et al. 2010, A&A, 522, A115
- Schleicher, D. R. G. & Beck, R. 2013, submitted
- Schleicher, D. R. G., Schober, J., Federrath, C., Bovino, S., & Schmidt, W. 2013, New Journal of Physics, 15, 023017
- Schlickeiser, R. 2012, Physical Review Letters, 109, 261101
- Schober, J., Schleicher, D., Federrath, C., et al. 2012a, ApJ, 754, 99
- Schober, J., Schleicher, D., Federrath, C., Klessen, R., & Banerjee, R. 2012b, Phys. Rev. E, 85, 026303
- Sedov, L. I. 1946, Dokl. Akad. Nauk. SSSR,, 52, 17
- Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics
- Semelin, B. & Combes, F. 2005, A&A, 441, 55
- Shiraishi, M., Nitta, D., Yokoyama, S., & Ichiki, K. 2012, J. Cosmology Astropart. Phys., 3, 41
- Sigl, G., Olinto, A. V., & Jedamzik, K. 1997, Phys. Rev. D, 55, 4582
- Simon, A. 1955, Physical Review, 100, 1557
- Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730
- Spitzer, L. 1956, Physics of Fully Ionized Gases (Interscience)
- Subramanian, K. 1997, ArXiv e-prints, arXiv:astro-ph/9708216
- Sur, S., Federrath, C., Schleicher, D. R. G., Banerjee, R., & Klessen, R. S. 2012, MNRAS, 423, 3148
- Taylor, G. 1950, Royal Society of London Proceedings Series A, 201, 159
- Turk, M. J., Oishi, J. S., Abel, T., & Bryan, G. L. 2012, ApJ, 745, 154
- Turner, M. S. & Widrow, L. M. 1988, Phys. Rev. D, 37, 2743
- Vogelsberger, M., Genel, S., Sijacki, D., et al. 2013, ArXiv e-prints
- Wise, J. H., Turk, M. J., & Abel, T. 2008, ApJ, 682, 745
- Xu, H., O’Shea, B. W., Collins, D. C., et al. 2008, ApJ, 688, L57