Open Access
Issue
A&A
Volume 691, November 2024
Article Number A280
Number of page(s) 12
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/202451139
Published online 19 November 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Local spiral galaxies are strongly magnetised, with a magnetic energy density in equipartition with turbulent kinetic and cosmic ray energy densities (Beck et al. 2019). Galactic dynamo theory (see Brandenburg & Ntormousi 2023, for a review) has been very successful in explaining cosmic growth and the observed large-scale coherence of the field, predicting growth rates of the order of 2 Gyr−1 for the large-scale field (Beck et al. 1996). However, recent observations suggest field strengths > μG at high redshifts as z ∼ 1 (Bernet et al. 2008; Mao et al. 2017), z ∼ 2.6 (Geach et al. 2023), and up to z ∼ 5.6 (Chen et al. 2024). Strongly magnetised galaxies at a cosmic time of about 1 − 2 Gyr pose a critical challenge: either the seed fields for the galactic dynamo are stronger than previously thought or there is an additional magnetic field amplification process at high redshift.

Numerical simulations of cosmic structure formation predict that massive galaxies are fed by filamentary cold accretion, penetrating through the hot circumgalactic medium (CGM) (Fardal et al. 2001; Kereš et al. 2005; Dekel & Birnboim 2006; Dekel et al. 2009). Recent high-resolution simulations of these cold streams have focused on their small-scale structure and hydrodynamical stability (Mandelker et al. 2019, 2020a; Berlok & Pfrommer 2019; Aung et al. 2019, 2024; Ledos et al. 2024; Hong et al. 2024). Using magnetohydrodynamic (MHD) simulations with radiative cooling, Ledos et al. (2024) have revealed that an nG-strength field can be amplified up to μG values in the cold gas, due to the tangling and stretching of the field lines inside the mixing layer between the stream and the CGM.

In light of this result and the observational evidence for high-z galactic magnetisation, for this work we estimated the degree of magnetisation that can be obtained from magnetised gas inflow onto early galaxies. To this end, we constructed a simple analytical model, quantifying the build-up of magnetic energy in the galaxy from magnetised cold inflows.

First, we present the origin of the magnetic field amplification from the simulations of Ledos et al. (2024) and from new 3D simulations in Sect. 2. Then, we present the model in Sect. 3 and discuss the results in Sect. 4.

2. Magnetic field growth in cold streams

Among recent high-resolution simulations of cold streams in the CGM with idealised geometry, Ledos et al. (2024) is to date the only one to cover the evolution of cold streams with magnetic field and radiative cooling over a wide range of parameters. We used their simulation results as input for our analytical model. We briefly present the simulation suite below, along with the new 3D simulations. For an extensive description, we refer to Ledos et al. (2024).

2.1. The simulation set-up

The stream is positioned at the centre of a box defined by Lx × LyLz) = 64 Rs × 32 Rs(×64 Rs) with Rs representing the stream radius. The stream axis is oriented in the y-direction. The stream traverses the simulation box repeatedly due to periodic boundary conditions. Initially, the cold stream and the hot CGM gas are uniformly distributed in isobaric equilibrium, with a density contrast δ = ρs/ρcgm. This set-up is consistent with previous idealised cold stream simulations (Mandelker et al. 2016, 2019, 2020a; Berlok & Pfrommer 2019; Vossberg et al. 2019).

The initial conditions for the density and the magnetic field are shown in the top panels of Fig. 1. The physical properties are assumed to match those of the stream and CGM at the virial radius Rv = 100 kpc of a typical massive halo of 1012M at z ∼ 2, consistent with cosmological simulations (Goerdt et al. 2010), empirical models (Dekel et al. 2013; Mandelker et al. 2020b), and observations (e.g. Martin et al. 2012; Daddi et al. 2021; Emonts et al. 2023; Zhang et al. 2023). This yields ns = 10−3, 10−2, 10−1 cm−3, along with two density ratios δ ≡ ρs/ρcgm = 30, 100, a set of stream and CGM metallicities ( Z s , Z cgm ) = ( 10 1.5 , 10 1 ) Z $ (Z_{\mathrm{s}}, Z_{\mathrm{cgm}}) = \left(10^{-1.5}, 10^{-1}\right)\,\mathrm{Z}_{\odot} $, and three different stream Mach numbers ℳ = ℳcgm = 0.5, 1, 2. We fix the stream radius to Rs = 1 kpc. In all cases, the magnetic field is initially uniform with a ∼nG strength, which yields a plasma beta (ratio of thermal pressure to magnetic pressure) β = p/pmag = 105. Magnetic fields with angle of 0°, 45°, and 90° with respect to the stream are also considered. Initial perturbations of the velocity field are introduced at the interface of the stream and the CGM.

thumbnail Fig. 1.

Time evolution of the stream and the magnetic field lines for a 3D simulation (with ns = 0.01 cm−3). Left: Maps showing volume rendering of the density field and the magnetic field lines. The lines are coloured by magnetic field strength. Right: Same as the left panel, but with only the magnetic field lines.

Since magnetic fields can behave differently in 3D than in 2D, we ran two 3D simulations to validate our 2D results, extending the 2D stream into a 3D cylinder. The 3D simulations are initialised with ns = {0.001,0.01} cm−3, δ = 30, and ℳ = 1, employing a radiative cooling model from Ferland et al. (2017) without self-shielding or thermal conduction to reduce computational time1.

It should be noted that δ = 30 represents a limit for the CGM properties, leading to relatively low temperature (Tcgm ∼ 3 × 105 K) and sound speed (ccgm ∼ 120 km s−1) and high density (ncgm ∼ 3 × 10−5 − 3 × 10−3 cm−3) (Dekel & Birnboim 2006; Mandelker et al. 2020b). However, this value of δ does not significantly affect the magnetic field growth and the stream’s evolution (Ledos et al. 2024), and has the advantage of reducing the computational cost of the simulations.

2.2. Stream evolution

Figure 1 illustrates the evolution of the stream density and of the magnetic field lines in the 3D simulation with ns = 0.01 cm−3. The time is in units of the sound crossing time in the stream t0 = Rs/cs, with cs the sound speed in the stream. For the simulations in Fig. 1, t0 ∼ 45 Myr, with a virial crossing time of tv = Rv/vs, 0 ∼ 18 t0. As the stream progresses, the shear at the interface between the stream and the CGM triggers the Kelvin-Helmholtz instability (KHI). The growth of the KHI creates a mixing layer. The fate of the stream is determined by the ratio of the cooling time in the mixing layer to the dynamical time,

ξ = t cool t dyn , $$ \begin{aligned} \xi =\frac{t_{\mathrm{cool}}}{t_{\mathrm{dyn}}}, \end{aligned} $$(1)

with tcool = (Tmixkb)/[(γ−1)nmixΛnet, mix], and tdyn ∼ Rs/vs, 0, where Tmix, nmix, and Λnet, mix are the mixing layer gas temperature, number density, and net cooling rate, estimated as in Begelman (1990) and Hillier & Arregui (2019). The term vs, 0 is the stream velocity. If cooling dominates (i.e. ξ ≲ 1), then the CGM gas in the mixing layer condenses onto the stream, leading to a growth of the stream mass. This situation is illustrated in Fig. 1, where the simulations shown have ξ ∼ 6 × 10−2. In cases where cooling is not dominant, Ledos et al. (2024) show that magnetic fields and thermal conduction can stabilise the stream against KHI up to ξ ∼ 8. For higher ξ values, the thermal conduction timescale becomes shorter than tcool and tdyn, causing the stream to diffuse into the hot CGM, reducing its cold gas accretion rate.

2.3. Magnetic field amplification mechanisms

Ledos et al. (2024) showed, in the situation described in the previous section, that two mechanisms lead to magnetic field growth: (1) the coherent stretching of the field lines due to the velocity shear at the interface between the stream and the CGM and (2) the tangling of field lines due to the growth of the KHI (including an incoherent stretching of field lines by turbulence or vortexes). We confirm that mechanisms (1) and (2) are dominant and well-resolved in our simulations by analysing the power spectra in Appendix A. Additionally, the condensation of CGM gas onto the stream can lead to a compression mechanism, which we discuss at the end of this section.

Figure 1 shows the stretching and tangling of the magnetic field lines around the stream. The magnetic field is initially amplified by the stretching of the field lines at t ∼ 0.4 t0 ∼ 0.02 tv. Following the growth of the KHI, the field lines tangle, and the magnetic field is mixed within the stream. This tangling has already developed well at the early time of t ∼ 2 t0 ∼ 0.1 tv. Appendix B compares the results between the 2D and 3D simulations. We demonstrate that at later simulation times, the 3D simulations allow more efficient mixing of the magnetic field in the stream compared to 2D simulations with the same flow parameters. Some simulations in Ledos et al. (2024) have a magnetic field parallel to the stream, which results in a magnetic field growth only by tangling as there is no magnetic field component to be coherently stretched by shear at the interface of the CGM and the stream. We refer to this configuration as a no-stretching amplification case, which represents a lower limit for the magnetic field amplification.

The field amplification is approximated as a stretching flux tube (see e.g. Spruit 2013). The details of the model are derived and discussed in Ledos et al. (2024, Appendix D). The magnetic field growth takes the form of

B th ( t ) B 0 = 1 2 [ ( β R s R s + v c t ) 2 + 4 ( β + 1 ) ] 1 / 2 1 2 β R s R s + v c t , $$ \begin{aligned} \frac{B_{\mathrm{th}}(t)}{B_0} = \frac{1}{2}\left[\left(\beta \frac{R_{\mathrm{s}}}{R_{\mathrm{s}} + v_{\mathrm{c} }t}\right)^2 + 4(\beta +1)\right]^{1/2} - \frac{1}{2} \beta \frac{R_{\mathrm{s}}}{R_{\mathrm{s}} + v_{\mathrm{c} }t}, \end{aligned} $$(2)

where Bth is the theoretical magnetic field, and β = 105 the plasma beta. The term Rs/(Rs+vct) indicates the stretching of the field lines over time, with vc the characteristic velocity of the mechanism responsible for the amplification. In the case where field lines are stretched due to shear (mechanism 1), vc = vs, 0 (i.e. the shear velocity between the stream and the CGM). In the case of the tangling of the field lines due to mixing (no-stretching case), vc = vturb (Eq. 34 Ledos et al. 2024), where vturb denotes the turbulent velocity in the mixing layer.

In Fig. 2 we present the time profile of the mean magnetic field in the stream for the 3D and 2D simulations, along with the model derived from Eq. (2). The magnetic field growth follows the same functional dependence on time in 2D and 3D. However, when cooling dominates (ns = 10−2 cm−3), the 3D simulation predicts a higher magnetic field strength compared to its 2D counterpart. This occurs because in 3D the mixing layer surrounds the stream, leading to a larger surface area for the stream to incorporate magnetic fields through the condensation of CGM gas. When condensation is subdominant (ns = 10−3 cm−3), there is no significant difference of magnetic strength between the 2D and 3D simulations. Given that the model (Eq. 2) is independent of ns, it reasonably reproduces the magnetic field growth in terms of shape and order of magnitude. Figure 2 also depicts the magnetic field growth for a 2D simulation in the no-stretching case (only tangling and incoherent stretching, corresponding to a magnetic field initially parallel to the stream). By replacing the shear velocity vs, 0 in Eq. (2) with the turbulent velocity, Bth closely matches the simulation result (compare the grey dotted and dashed lines). This highlights that in the absence of stretching, the primary mechanism for the field amplification is the tangling of the magnetic field lines due to turbulence.

thumbnail Fig. 2.

Time profiles of the magnetic field in the stream for 2D and 3D simulations (dotted and solid lines, respectively). Simulations are shown for ns = 10−2 and 10−3 cm−3 in purple and orange, respectively. The dotted grey line stands for a 2D simulation with ns = 10−2 cm−3 and with a magnetic field parallel to the stream (no amplification due to stretching). The dashed lines show the model results from Eq. (2) with the same Mach number as the simulations (ℳ = 1).

Following these results, we consider that the magnetic field growth due to gas condensation is, in most cases, small compared to the effects of stretching and tangling. If compression was the dominant mechanism, the magnetic field growth would primarily depend on the efficiency of cooling in the mixing layer (i.e. the parameter ξ), and there would be almost no magnetic field amplification when the stream does not condense CGM gas (ξ > 1). However, (1) the magnetic field growth is both significant and independent of ξ when the magnetic field is not parallel to the stream (Ledos et al. 2024, Fig. 12), and (2) when the magnetic field is parallel to the stream the magnetic field growth can be well described by our model that accounts for the tangling of the field lines (see Fig. 2). The dependence on ξ then arises from the turbulent velocity vturb. Based on these two arguments, we argue that compression is subdominant, and we do not include it in our simple model.

3. A simple analytical model

In this section we present a simple analytical model to estimate the potential contribution of magnetised cold streams to the galactic magnetic field.

3.1. Model summary

We illustrate the main idea of the model in Fig. 3. The model can be decomposed into three physical steps:

  • Step 1: Growth and accumulation of the magnetic field in the stream. Once the cold stream enters the hot halo, the velocity shear at the interface of the cold stream and the hot CGM triggers the KHI. From the growth of the KHI, a mixing layer forms at the interface. The magnetic field lines are tangled and stretched by shear and KHI-generated turbulence, leading to an amplification of the magnetic field in the mixing layer. Concurrently, if the cooling time in the mixing layer is shorter than the shearing time, the CGM gas condenses onto the stream. The condensation incorporates the magnetic field from the mixing layer into the cold stream. This process occurs relatively quickly, with t ≲ 8 t0 ∼ 0.4 tv. As a consequence of the interaction between the stream and the CGM, the amplified magnetic field Bs within the stream and the stream’s cold mass flux Ės are expressed as follows:

    B s = B s ( X ) , $$ \begin{aligned} B_{\mathrm{s}}&= B_{\mathrm{s}}(\boldsymbol{X}), \end{aligned} $$(3)

    d M s d t = F ( X ) . $$ \begin{aligned} \dfrac{\mathrm{d} M_{\mathrm{s} }}{\mathrm{d} t}&= F(\boldsymbol{X}). \end{aligned} $$(4)

    To keep generality, one may use any functions for Bs(X) and F(X) with X being a vector of any relevant inputs variables. The specific form of these functions is defined in Sect. 3.3 within the context of our simulations.

  • Step 2: Growth of the magnetic field in the stream from compression. The cold stream falling towards the central galaxy undergoes compression due to the gravitational potential of the dark matter halo. This compression is expected to amplify the magnetic field within the stream further. Additionally, as the stream grows due to the condensation of CGM gas, its velocity increases, resulting in a higher rate of cold gas accretion onto the galaxy. Our suite of simulations, presented in Sect. 2, does not account for this process, but we included it in the analytical model by considering the results from other works. Recent studies from Hong et al. (2024) and Aung et al. (2024) show that the potential can indeed enhance the cold gas accretion rate by a factor of ≲3.5. Because of the halo potential’s influence, we defined the cold mass flux in reaching the galaxy as

    d M in d t = η F ( X ) , $$ \begin{aligned} \dfrac{\mathrm{d} M_{\mathrm{in} }}{\mathrm{d} t} = \eta F(\boldsymbol{X}), \end{aligned} $$(5)

    and its associated mean magnetic field Bin as

    B in = α B s ( X ) , $$ \begin{aligned} B_{\mathrm{in} } = \alpha B_{\mathrm{s}}(\boldsymbol{X}), \end{aligned} $$(6)

    where η ∈ [0, ∞] represents an arbitrary efficiency rate for the cold inflow due to the halo potential, and α is the compression coefficient that accounts for the magnetic field amplification resulting from compression.

  • Step 3: Accumulation of magnetic energy in the galaxy. We assumed that the galaxy is fed by magnetised cold inflows with accretion rate in and magnetic field Bin. Over time, this inflow accumulates magnetic energy within the galaxy, thereby increasing the mean galactic magnetic field Bg. Clearly, galaxies can also lose magnetic flux through outflows, an effect not included in the model. Since outflow properties depend on the galactic properties in a non-trivial way, including them would require extensions to the model that are beyond the scope of this work. However, in Appendix C we considered the simple case of a continuous outflow term. More detailed estimates of the effect of outflows on the galaxy’s magnetisation will be the subject of future work.

thumbnail Fig. 3.

Illustrative sketch of the main processes of the model. The cold stream arrives initially from the intergalactic medium (IGM) with a mass flow rate Ms, 0 and a magnetic field B0 ∼ nG. Step 1: As the cold stream enters the halo, a mixing layer forms at the interface of the CGM and the stream. Magnetic field lines are tangled and stretched in the mixing layer, resulting in magnetic field amplification. Simultaneously, CGM gas in the mixing layer condenses onto the stream, accumulating magnetic fields. From the interaction of the stream and the CGM, the cold mass flux and the magnetic field of the stream become F and Bs. Step 2: Upon falling into the halo, the cold stream narrows and is squeezed by the dark matter potential. Consequently, magnetic fields in the stream are also compressed, further amplifying it by a factor α. The stream is also expected to entrain CGM gas while falling through the halo potential, increasing its mass flux by an efficiency factor η. Step 3: Ultimately, if the stream reaches the galaxy, it continuously accretes cold gas at a rate in with an associated magnetic energy inflow rate Ėm,in to the galaxy, elevating the galactic mean magnetic field. By default, we did not consider outflow; however, we introduced a simple outflow term in Appendix C.

3.2. Derivation

A general definition of the evolution of the magnetic energy in the galaxy can be written as

d E m , g d t = d E m , in d t + d E m , amp d t d E m , diff d t d E m , out d t , $$ \begin{aligned} \dfrac{\mathrm{d} E_{\mathrm{m,g} }}{\mathrm{d} t} = \dfrac{\mathrm{d} E_{\mathrm{m,in} }}{\mathrm{d} t} + \dfrac{\mathrm{d} E_{\mathrm{m,amp} }}{\mathrm{d} t} - \dfrac{\mathrm{d} E_{\mathrm{m,diff} }}{\mathrm{d} t} - \dfrac{\mathrm{d} E_{\mathrm{m,out} }}{\mathrm{d} t}, \end{aligned} $$(7)

where the terms Ėm,in, Ėm,amp, Ėm,diff, and Ėm,out stand for the variation of galactic magnetic energy by inflows, the amplification inside the galaxy, diffusion processes, and removal from outflows, respectively. For simplicity, in this paper we focus only on the inflow component Em, in, leaving all other terms equal to zero. We provide the derivation of the model with a simple outflow component Em, out in Appendix C. We can then rewrite Ėm,g as

d E m , g d t = d E m , in d t = d E m , in d M in d M in d t · $$ \begin{aligned} \dfrac{\mathrm{d} E_{\mathrm{m,g} }}{\mathrm{d} t} = \dfrac{\mathrm{d} E_{\mathrm{m,in} }}{\mathrm{d} t} = \dfrac{\mathrm{d} E_{\mathrm{m,in} }}{\mathrm{d} M_{\mathrm{in} }}\dfrac{\mathrm{d} M_{\mathrm{in} }}{\mathrm{d} t}\cdot \end{aligned} $$(8)

Knowing that,

E m , in = B in 2 2 V in = B in 2 2 M in ρ in , $$ \begin{aligned} E_{\mathrm{m,in} } = \frac{B_{\mathrm{in} }^2}{2} V_{\mathrm{in} } = \frac{B_{\mathrm{in} }^2}{2} \frac{M_{\mathrm{in} }}{\rho _{\mathrm{in} }}, \end{aligned} $$(9)

where Vin, Min, and ρin are the volume, mass, and density of the gas inflow reaching the galaxy, respectively. We end up with a general formulation as

d E m , g d t = B in 2 2 ρ in d M in d t = C in . $$ \begin{aligned} \dfrac{\mathrm{d} E_{\mathrm{m,g}}}{\mathrm{d} t} = \frac{B_{\mathrm{in} }^2}{2 \rho _{\mathrm{in} }} \dfrac{\mathrm{d} M_{\mathrm{in} }}{\mathrm{d} t} = C_{\mathrm{in} }. \end{aligned} $$(10)

Assuming that Cin is invariant in time yields

E m , g ( t ) = C in t + E m , g , 0 . $$ \begin{aligned} E_{\mathrm{m,g} }(t) = C_{\mathrm{in} } t + E_{\mathrm{m,g,0} }. \end{aligned} $$(11)

For Cin, using Eqs. (5) and (6) for Bin and in gives

C in = α B s 2 2 ρ s η F ( X ) , $$ \begin{aligned} C_{\mathrm{in} } = \frac{\alpha B_{\mathrm{s} }^2}{2 \rho _{\mathrm{s} }} \eta F(\boldsymbol{X}), \end{aligned} $$(12)

where we use the compression coefficient such that ρin ≡ αρs. Finally, the mean magnetic field in the galactic volume Vg can be recovered as

B g ( t ) = ( 2 E m , g ( t ) V g ) 1 / 2 . $$ \begin{aligned} B_{\mathrm{g}}(t) = \left(2\frac{E_{\mathrm{m,g} }(t)}{V_{\mathrm{g}}}\right)^{1/2}. \end{aligned} $$(13)

From Eq. (13), the model predicts a continuous growth of the galactic magnetic field.

The above derivation is kept general, such that any estimation of the cold mass inflow in and the inflow magnetic field Bin can be used and plugged in, either in Eq. (10) or Eq. (13).

3.3. Step 1: Input values from simulations

We now outline how we utilised the simulations to determine values for the cold stream mass flux F and its magnetic field Bs resulting from the interaction with the CGM (Step 1). To determine the cold mass flux subsequent to the interaction with the CGM, we employed a fitting function based on the cold mass flux rate obtained from the 2D simulations suite (Ledos et al. 2024). This function is expressed as

F ( X ) M ˙ s , 0 = 1 a ξ b , $$ \begin{aligned} \frac{F(\boldsymbol{X})}{\dot{M}_{\mathrm{s},0}} = 1 - a\,\xi ^b, \end{aligned} $$(14)

where a ∼ {0.03,0.11,0.05} and b ∼ {0.95,0.24,0.12}, for ℳ = {0.5,1,2}, respectively. Here X = (ℳ,ξ). It is important to note that the fits were performed based on the mass flux at t = tend, which therefore represents a lower limit scenario for the cold mass survival rate.

The comparison of the fitted model and the simulations is shown in Fig. 4. As described in Sect. 2.1, the cold stream mass flux is expected to decrease with increasing ξ. Moreover, for ξ ≳ 1, higher Mach numbers correspond to lower KHI growth rates, resulting in higher values of F in the simulations. The fitted model from Eq. (14) effectively predicts the simulation outcomes, demonstrating its suitability as a definition for F.

thumbnail Fig. 4.

Cold mass flux rate normalised by the initial cold mass flux of the simulations. Each point represents one simulation with a fixed ξ, ℳ, and initial magnetic angle ϕ with respect to the stream axis. Red dots, blue squares, and cyan diamonds refer to 2D simulations with ℳ = 2, 1, 0.5, respectively. Using the same colour legend, the fitted model is represented by coloured curves. For a fixed ξ and ℳ, the two points represent simulations with an initial magnetic field angle of 90° and 45° with the highest values associated with the 90° angle in all cases.

For Bs, we assumed the magnetic field to be at the mean saturated value following the initial rapid growth, Bs = ⟨Bth⟩, with the average done over the full simulation time t = 22 t0 and with Bth expressed in Eq. (2) for ξ < 8. For ξ > 8, to account for the special case where the stream is diffused inside the CGM due to thermal conduction, we used a linear function dependent on Mach numbers fitted from values of Ledos et al. (2024). This gives,

B s = { B th if ξ 8 , B 0 ( 34 M 16 ) if ξ > 8 . $$ \begin{aligned} B_{\mathrm{s}} = \left\{ \begin{array}{ll} \left\langle B_{\mathrm{th} }\right\rangle&\mathrm{if} \,\, \xi \le 8,\\ B_0 \left(34\mathcal{M} -16\right)&\mathrm{if} \,\, \xi >8. \end{array}\right. \end{aligned} $$(15)

For ξ > 8, anisotropic thermal conduction is expected to diffuse both the stream and the perturbations, thereby suppressing the KHI. Thermal conduction occurs along the magnetic field lines, which are stretched at the stream’s interface due to the velocity difference. Faster streams result in more stretching and less efficient thermal conduction (see Ledos et al. 2024, Sect. 4.1.2). Consequently, the stream’s diffusion continues over time only for relatively slow streams (ℳ ≤ 0.5). For ℳ > 0.5, the stream’s diffusion becomes rapidly inefficient, so the magnetic field, initially amplified by stretching, remains in the cold phase, leading to a small amplification of the mean magnetic field within the stream.

As discussed in Sect. 2, we also considered a no-stretching case for the magnetic field amplification. In that case, instead of using Eq. (2), we wanted to account for the impact of thermal conduction and used the power-law fit from Ledos et al. (2024), giving

B s 0.13 ξ 0.1 . $$ \begin{aligned} B_{\mathrm{s}} \sim 0.13\,\xi ^{-0.1}. \end{aligned} $$(16)

We note that Bs is based on the 2D simulations results (see Ledos et al. 2024, Fig. 12). As discussed in Sect. 2, this means that our Bs may be an underestimation compared to our 3D results when cooling is dominant (ξ < 1).

4. Results and discussion

We first apply our model to a general case of an accreting galaxy and then consider the particular case of the galaxy studied by Geach et al. (2023).

4.1. Growth of the galactic magnetic field

We consider a scenario without outflows (Eq. 11) where the impact of the halo potential (Step 2) is neglected (i.e. η = α = 1). For the galaxy volume Vg, we first consider its halo mass to be of the order of 1012M to ensure that its CGM is hot and that it is fed by cold streams between z ∼ 3 and z ∼ 2 (Dekel & Birnboim 2006; Aung et al. 2024). This redshift range corresponds to ∼1.15 Gyr, for which we assume for simplicity that the cold accretion lasts for about 1 Gyr. We use the empirical relation of Girelli et al. (2020) for the halo to stellar mass ratio, and the empirical relation of Yang et al. (2021) for the stellar mass-to-galactic radius conversion2. For our considered halo and redshift, it gives us a stellar mass of M* ∼ 2 × 1011M, and Rg ∼ 3 kpc. The galaxy’s disc height is fixed at dg = 0.1 Rg, as suggested by observations of high-redshift galaxies from the Hubble Space Telescope (HST) catalogue (Elmegreen et al. 2017, Fig. 19). Notably, their values of dg/Rg range from ∼2% to ∼25% at z = 2.5 − 2. The effects of an extended height and radius are discussed in Sect. 4.3.

By fixing δ = 30 and ℳ = 1, F and Bs become solely functions of Rs and ns, the cold stream radius and number density at the virial radius of the halo. Fixing Rs, the galactic magnetic field can then be plotted in function of s,0 by varying only ns.

Figure 5 illustrates the galactic magnetic field strength after 1 Gyr of accretion by one cold stream of radius Rs as a function of the cold stream mass flux. The results are plotted for three values of the stream radius and the two values of Bs growth.

thumbnail Fig. 5.

Mean galactic magnetic field from our model as a function of the stream cold mass flux entering the halo s,0. The cold mass flux is obtained by varying its number density ns ∈ [0.001,0.1] cm−3 for different stream radii Rs = (0.3, 3, 30) kpc (plain, dashed and dotted lines) and for different Mach number ℳ = 0.5, 1, 2 (light to dark colours). The red and blue colour maps show results for the default case Bs and the no-stretching case (see Sect. 3.3), respectively. The vertical grey line shows the maximum cold mass flux simulated in the simulations suite.

From the fiducial model (default Bs), even a small cold stream of radius Rs = 0.3 kpc yields μG galactic magnetic field values for M ˙ s , 0 0.1 M $ \dot{M}_{\mathrm{s},0}\gtrsim 0.1\,M_{\odot} $ yr−1. A lower s,0 corresponds to a lower stream number density ns (higher ξ), resulting in a shorter thermal conduction timescale and a longer cooling time3. Consequently, for narrow streams with Rs = 0.3 kpc, thermal conduction only diffuses the stream at relatively low accretion rates per stream of M ˙ s , 0 < 10 2 M $ \dot{M}_{\mathrm{s},0} < 10^{-2}\,M_{\odot} $ yr−1. This stream diffusion inhibits the magnetisation of the galaxy only below this low accretion rate threshold. In the no-stretching case, the magnetic field amplification occurs for Rs = 3 kpc with M ˙ s , 0 1 M $ \dot{M}_{\mathrm{s},0}\gtrsim 1\,M_{\odot} $ yr−1.

Constraining both the minimum size and the number of cold streams that feed massive galaxies remains a significant challenge, primarily due to the current resolution limits for cold gas in the CGM of cosmological simulations (e.g. Bennett & Sijacki 2020). Given the connectivity between galaxies in cosmological simulations, it is expected that massive star-forming galaxies are fed by approximately two or more streams (Galárraga-Espinosa et al. 2023), with relatively large radii when entering the halo (Rs ≳ 10 kpc). The streams become narrower as they approach the galaxy due to the pressure gradient in the CGM and gravitational focusing caused by the halo potential. This cold stream accretion picture is commonly predicted in simulations and empirical models (e.g. Dekel et al. 2009; Rosdahl & Blaizot 2012; Mandelker et al. 2018, 2020b; Lu et al. 2024). With at least two large cold streams accreting onto a galaxy for about 1 Gyr, our model naturally predicts Bg ≳ 10 μG for both cases of Bs.

4.2. Applying the model to the galaxy of Geach et al. (2018)

We now apply our model to the observed lensed star-forming galaxy (Geach et al. 2018) at redshift z ∼ 2.6 whose magnetic field is estimated at ∼500 μG or less (Geach et al. 2023). The galaxy’s dynamical mass is Mdyn ∼ 8.1 × 1010M, with a radius of about 2.5 kpc and a star formation rate (SFR) of the order of ∼1000 M yr−1. We assume that the galaxy’s host dark matter halo follows a Navarro–Frenk–White (NFW) profile. Then, assuming a polytropic equation of state P cgm ρ cgm γ $ P_{\mathrm{cgm}}\propto \rho_{\mathrm{cgm}}^{\gamma\prime} $, the radial density profile follows (Komatsu & Seljak 2001, Eq. 19)

y ( r ) = ρ cgm ( r ) ρ cgm , 0 = [ 1 + 3 σ 0 γ γ 1 c m ( c ) ( log ( x + 1 ) x 1 ) ] 1 γ 1 , $$ \begin{aligned} y(r) = \frac{\rho _{\mathrm{cgm}}(r)}{\rho _{\mathrm{cgm},0}} = \left[1 +\frac{3}{\sigma _0}\frac{\gamma \prime }{\gamma \prime -1} \frac{c}{m(c)}\left(\frac{\log {(x+1)}}{x}-1\right) \right]^{\frac{1}{\gamma \prime -1}}, \end{aligned} $$(17)

where r is the radius, c is the concentration parameter of the halo, x ≡ rc/Rv, m(c) = log(c + 1) + c/(c + 1), and ρcgm, 0 = ρcgm(r = 0). The constant σ0 and the polytropic index γ′ are empirical functions4 of c. Further assuming pressure equilibrium between the stream and the hot CGM gas Ps = Pcgm (Aung et al. 2024), we can define the compression parameter α as the ratio of the stream density between 0.1 Rv and Rv

α = ρ s ( 0.1 R v ) ρ s ( R v ) = [ y ( 0.1 R v ) y ( R v ) ] γ T s ( R v ) T s ( 0.1 R v ) , $$ \begin{aligned} \alpha = \frac{\rho _{\mathrm{s}}\left(0.1\,R_{\mathrm{v} }\right)}{\rho _{\mathrm{s}}(R_{\mathrm{v} })} = \left[\frac{y\left(0.1\,R_{\mathrm{v} }\right)}{y(R_{\mathrm{v} })}\right]^{\gamma \prime } \frac{T_{\mathrm{s}}(R_{\mathrm{v} })}{T_{\mathrm{s}}\left(0.1\,R_{\mathrm{v} }\right)}, \end{aligned} $$(18)

where we assume Ts(Rv)/Ts(0.1 Rv)∼15.

We define two scenarios for our model: a lower limit scenario and a fiducial scenario. For the lower limit scenario, the mass accretion is assumed at its minimum, such that M ˙ s , 0 = SFR = 1000 M $ \dot{M}_{\mathrm{s},0} = {\mathrm{SFR}} = 1000\,M_{\odot} $ yr−1. For the model parameters, the accretion efficiency η and the concentration parameter c are both set to unity leading to α ∼ 22. The no-stretching case is assumed for Bs. We note that c = 1 only represents a numerical lower limit of the model and is a factor of 2 − 4 lower than estimation with Λ-Cold Dark Matter (Λ-CDM) (e.g. Prada et al. 2012; Dutton & Macciò 2014; Ludlow et al. 2016).

For the fiducial scenario, we adopt a ratio of mass accretion rate in to star formation rate ϵ = M ˙ in / SFR 0.2 $ \epsilon = \dot{M}_{\mathrm{in}}/{\mathrm{SFR}} \sim 0.2 $. This ratio accounts for the fact that the cold gas accreted by the galaxy is not instantaneously transformed into stars, but instead contributes to the gas reservoir of the galaxy, sustaining star formation over longer timescales. The value of 0.2 is motivated by those derived from observations by Daddi et al. (2021, ϵ ∼ 0.12), Emonts et al. (2023, ϵ ∼ 0.56), and Zhang et al. (2023, ϵ ∼ 0.12). It should be noted that this ratio is not well constrained and can vary significantly, for example, during a galaxy’s sudden starburst phase, where ϵ can reach ∼12 for SFR ∼ 1200 M yr−1 (Fu et al. 2021). Further quantification of the ratio ϵ is left for future work as part of the development of a less simplistic model. The resulting total accretion rate onto the galaxy is M ˙ in SFR ϵ 1 = 5000 M $ \dot{M}_{\mathrm{in}} \sim \mathrm{SFR} \, \epsilon^{-1} = 5000\,M_{\odot} $ yr−1. Following Aung et al. (2024), we consider an efficiency parameter η = 3, yielding M ˙ s , 0 = M ˙ in η 1 = 5000 η 1 M $ \dot{M}_{\mathrm{s},0}= \dot{M}_{\mathrm{in}} \eta^{-1} = 5000\, \eta^{-1}\,M_{\odot} $ yr−1 at the virial radius. We assume a concentration parameter c = 10, which leads to a compression coefficient of α ∼ 450.

In both scenarios we assume that the galaxy accretes through three cold streams with a number density at the virial radius of ns = 0.1 cm−36. This results in stream radii of 19 kpc and 35 kpc for the lower limit and fiducial scenarios, respectively. We also apply our simple model extension with outflow presented in Appendix C using a SFR = 1000 M yr−1.

Figure 6 illustrates the time evolution of the galactic magnetic field Bg assuming that the accretions start at redshift z = 3. In the lower limit scenario, owing to the extreme assumptions, the galactic magnetic field is raised to ∼250 μG values by z = 2.6. With the reasonable assumptions of the fiducial model, the galaxy reaches 500 μG rapidly at around t ∼ 10 Myr. If we consider the expulsion of galactic magnetic energy by a constant outflow (see Appendix C), the initial growth is the same as before; however, the galactic magnetic field reaches a plateau. This behaviour illustrates the necessity to embed this inflow model into a more sophisticated framework that includes the complex physics of galactic outflows. Within the limit of the model assumptions, cold streams emerge as plausible candidates for rapidly magnetising galaxies, potentially elucidating the puzzlingly large magnetic field estimates by Geach et al. (2023).

thumbnail Fig. 6.

Time evolution of the mean galactic magnetic field of our model. The corresponding redshift is also given at the top axis assuming that the accretion starts at z = 3. The blue and red curves represent the fiducial and lower-limit scenarios. The solid and dotted lines are models without and with outflows. The black line with the arrow shows the target at 500 μG or less from Geach et al. (2023).

4.3. Discussion

Numerous simplifications are assumed in our model, and we discuss here the potential impact of the missing mechanisms. Firstly, the model assumes that the cold streams directly reach the galaxy. However, cosmological simulations suggest that cold streams have an impact parameter and circularise into an extended ring around the galaxy (Danovich et al. 2015; Stewart et al. 2017). Such an extended ring would likely delay the buildup of magnetic energy in the galaxy. Additionally, upon reaching the galaxy or its outer ring, the cold streams are more likely to fragment, as suggested by observations of Hα clumps in galaxies at redshift z ∼ 1.5 − 2 (Fisher et al. 2017) and simulations (Ceverino et al. 2016; Mandelker et al. 2018; Aung et al. 2019). These studies associate these clumps with the creation of turbulent galactic discs (Ginzburg et al. 2022; Forbes et al. 2023). Such fragmentation and turbulence could lead to further amplification of the magnetic field through mechanisms similar to our Step 1 and Step 2, in addition to galactic dynamo processes. Additionally, we assume that cold accretion remains constant across redshifts. Theoretical and empirical models suggest a (1 + z)2.5 dependence of the accretion rate for z ≳ 2 (Neistein & Dekel 2008; Fakhouri et al. 2010; Dekel et al. 2013). While this would introduce only a factor of ∼1.2 to the final Bg within the redshift range considered in Sect. 4, such a redshift dependence might be necessary to accurately trace magnetic field evolution over a longer timescale. Extending the model to include this aspect is left for future work.

Our model only accounts for inflows. We estimated the possible effect of an outflow from supernovae feedback only in a simplified and limited case. Considering a wider framework of galaxy evolution where supernova outflows arise self-consistently is an obvious next step for our model. Active galactic nucleus (AGN) feedback can also remove magnetic energy from the galaxy, dispersing it into the CGM. Incorporating an AGN outflow term into Eq. (10) is a focus for future work. However, in the presence of outflows, it would be necessary to also consider re-accretion through the galactic fountain (e.g. Putman et al. 2012) within the inflow term. Therefore, while we do not account for AGN feedback, not accounting for re-accretion might lead to an overall underestimation of the galactic magnetic field growth by our model.

Our model also does not explicitly solve any diffusion processes. The simulation suite from Ledos et al. (2024) explicitly solves thermal conduction, and its effects are implicitly considered in our definition of the cold gas flux F after interaction with the CGM (Step 1, Eq. 14). As the stream becomes denser throughout the halo, the thermal conduction timescale increases towards the centre of the halo. Consequently, its impact should be negligible for Steps 2 and 3 of the model. Regarding magnetic field diffusion, Ohmic and ambipolar diffusion are also negligible in Step 1. The stream’s compression throughout the halo could lead to a decrease in its ionisation parameter, resulting in shorter ambipolar diffusion times. However, simulations by Aung et al. (2024) indicate that, with or without self-shielding, the streams maintain a temperature of ∼104 K, suggesting that the ionisation parameter should not significantly decrease. Ohmic diffusion is also expected to remain insignificant throughout the halo. Eventually, magnetic diffusion is likely to occur inside the galaxy and should be included in Eq. (10) as a future work.

Finally, our model does not account for the magnetic field structure in the galaxy. Previous theoretical models (e.g. Henriksen & Irwin 2016; Nixon et al. 2018) estimated the resulting structure and polarisation of the magnetic field for a galaxy embedded in a magnetised small halo (∼20 kpc halo radius). Interestingly, while they do not address the origin of the magnetic field in either the halo or the galaxy, a by-product of their model is an inflow or outflow of magnetic fields. An estimation of the galactic magnetic field structure is left for future work.

5. Conclusions

We demonstrate that the stretching and tangling of magnetic field lines are effective mechanisms for amplifying the magnetic field within cold streams through their interaction with the CGM. Within the framework of the simple assumptions of our model, cold streams emerge as plausible candidates for raising the galactic magnetic field from zero to > μG for massive high-redshift star-forming galaxies. With reasonable assumptions regarding the mass accretion rate derived from SFR, our models, both with and without outflows, predict a very rapid magnetisation of the galaxy (within 10 Myr) consistent with the ∼500 μG galactic magnetic field observed by Geach et al. (2023). This magnetisation timescale could also account for the observed galactic magnetic field at z = 5.6 (Chen et al. 2024).


1

As demonstrated in Ledos et al. (2024), thermal conduction has only a significant effect on magnetic field growth when ℳ ≤ 0.5 and ns ≤ 0.001 cm−3.

2

The Girelli et al. (2020) empirical estimation of the stellar-to-halo mass relation is for star-forming galaxies up to redshift z ∼ 4 from the COSMOS catalogue (Scoville et al. 2007) and the DUSTGRAIN–pathfinder simulations (Giocoli et al. 2018). The empirical relation from Yang et al. (2021) is for massive star-forming galaxies up to z ∼ 2.75 and uses the ASTRODEEP catalogue.

3

Lowering the radius, the number density, and the stream velocity can also lower s,0 and increase the parameter ξ. For a detailed quantitative description of the timescale comparison depending on the stream parameters, we refer to Ledos et al. (2024, Sect. 2.5).

4

For c = 10 and 1, this leads to σ0 ∼ 3.28, γ′∼1.19 and σ0 ∼ 1.55, γ′∼1.10, respectively (Komatsu & Seljak 2001, Eqs. 25 and 26).

5

In equilibrium state, the cold stream’s temperature is defined by the balance between radiative cooling and heating, leading to Ts ∼ 104 K. However, as the stream becomes denser inside the halo, Ts can slightly decrease, leading to Ts(Rv)/Ts(0.1 Rv)∼1 − 2 (Aung et al. 2024). This additional factor would slightly increase the parameter α.

6

We fixed the number density for the stream at a relatively high value compared to predictions from empirical models (Mandelker et al. 2020b; Aung et al. 2024) in order to avoid artificially large stream radii or numbers.

7

To save computational time, the power spectra for the 3D simulation are computed at a uniform resolution (2lmax − 1)3, where lmax the maximum level of refinement in the simulation.

8

The simulation uses static mesh refinement. Hence, if the stream grows with time, the mixing layer can be shifted to lower resolution regions, leading to the appearance of artificial peaks at 2 and 4 times lower wavelengths than the resolution limit.

9

Equipartition is not strictly observed in Fig. A.1. In the simulation, the magnetic field is amplified within the mixing layer, resulting in equipartition only in that region. Since the power spectra are computed over the entire computational domain rather than solely within the mixing layer, this introduces excess kinetic energy into the spectra.

Acknowledgments

N.L. acknowledges support from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program grant agreement No. 864361. E.N. acknowledges funding from the Italian Ministry for Universities and Research (MUR) through the “Young Researchers” funding call (Project MSCA 000074). Numerical computations were carried out on the Cray XC50 at the Center for Computational Astrophysics, National Astronomical Observatory of Japan, and the SQUID at the Cybermedia Center, Osaka University as part of the HPCI system Research Project (hp230089, hp240141). This work is supported in part by the MEXT/JSPS KAKENHI grant numbers 20H00180, 22K21349, 24H00002, and 24H00241 (K.N.). K.N. acknowledges the support from the Kavli IPMU, World Premier Research Center Initiative (WPI), UTIAS, the University of Tokyo.

References

  1. Aung, H., Mandelker, N., Nagai, D., Dekel, A., & Birnboim, Y. 2019, MNRAS, 490, 181 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aung, H., Mandelker, N., Dekel, A., et al. 2024, MNRAS, 532, 2965 [NASA ADS] [CrossRef] [Google Scholar]
  3. Beck, R., Brandenburg, A., Moss, D., Shukurov, A., & Sokoloff, D. 1996, ARA&A, 34, 155 [NASA ADS] [CrossRef] [Google Scholar]
  4. Beck, R., Chamandy, L., Elson, E., & Blackman, E. G. 2019, Galaxies, 8, 4 [NASA ADS] [CrossRef] [Google Scholar]
  5. Begelman, M. C. 1990, MNRAS, 244, 26 [NASA ADS] [Google Scholar]
  6. Bennett, J. S., & Sijacki, D. 2020, MNRAS, 499, 597 [NASA ADS] [CrossRef] [Google Scholar]
  7. Berlok, T., & Pfrommer, C. 2019, MNRAS, 489, 3368 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bernet, M. L., Miniati, F., Lilly, S. J., Kronberg, P. P., & Dessauges-Zavadsky, M. 2008, Nature, 454, 302 [NASA ADS] [CrossRef] [Google Scholar]
  9. Brandenburg, A., & Ntormousi, E. 2023, ARA&A, 61, 561 [CrossRef] [Google Scholar]
  10. Ceverino, D., Sánchez Almeida, J., Muñoz Tuñón, C., et al. 2016, MNRAS, 457, 2605 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chen, J., Lopez-Rodriguez, E., Ivison, R. J., et al. 2024, A&A, in press, https://doi.org/10.1051/0004-6361/202450969 [Google Scholar]
  12. Cho, J., Vishniac, E. T., Beresnyak, A., Lazarian, A., & Ryu, D. 2009, ApJ, 693, 1449 [NASA ADS] [CrossRef] [Google Scholar]
  13. Daddi, E., Valentino, F., Rich, R. M., et al. 2021, A&A, 649, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Danovich, M., Dekel, A., Hahn, O., Ceverino, D., & Primack, J. 2015, MNRAS, 449, 2087 [Google Scholar]
  15. Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451 [Google Scholar]
  17. Dekel, A., Zolotov, A., Tweed, D., et al. 2013, MNRAS, 435, 999 [Google Scholar]
  18. Dutton, A. A., & Macciò, A. V. 2014, MNRAS, 441, 3359 [Google Scholar]
  19. Elmegreen, B. G., Elmegreen, D. M., Tompkins, B., & Jenks, L. G. 2017, ApJ, 847, 14 [NASA ADS] [CrossRef] [Google Scholar]
  20. Emonts, B. H., Lehnert, M. D., Yoon, I., et al. 2023, Science, 379, 1323 [NASA ADS] [CrossRef] [Google Scholar]
  21. Fakhouri, O., Ma, C. P., & Boylan-Kolchin, M. 2010, MNRAS, 406, 2267 [CrossRef] [Google Scholar]
  22. Fardal, M. A., Katz, N., Gardner, J. P., et al. 2001, ApJ, 562, 605 [Google Scholar]
  23. Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. Astron. Astrofis., 53, 385 [NASA ADS] [Google Scholar]
  24. Fisher, D. B., Glazebrook, K., Abraham, R. G., et al. 2017, ApJ, 839, L5 [NASA ADS] [CrossRef] [Google Scholar]
  25. Forbes, J. C., Emami, R., Somerville, R. S., et al. 2023, ApJ, 948, 107 [NASA ADS] [CrossRef] [Google Scholar]
  26. Fu, H., Xue, R., Prochaska, J. X., et al. 2021, ApJ, 908, 188 [Google Scholar]
  27. Galárraga-Espinosa, D., Garaldi, E., & Kauffmann, G. 2023, A&A, 671, A160 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Geach, J. E., Ivison, R. J., Dye, S., & Oteo, I. 2018, ApJ, 866, L12 [NASA ADS] [CrossRef] [Google Scholar]
  29. Geach, J. E., Lopez-Rodriguez, E., Doherty, M. J., et al. 2023, Nature, 621, 483 [NASA ADS] [CrossRef] [Google Scholar]
  30. Ginzburg, O., Dekel, A., Mandelker, N., & Krumholz, M. R. 2022, MNRAS, 513, 6177 [NASA ADS] [CrossRef] [Google Scholar]
  31. Giocoli, C., Baldi, M., & Moscardini, L. 2018, MNRAS, 481, 2813 [Google Scholar]
  32. Girelli, G., Pozzetti, L., Bolzonella, M., et al. 2020, A&A, 634, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Goerdt, T., Dekel, A., Sternberg, A., et al. 2010, MNRAS, 407, 613 [Google Scholar]
  34. Henriksen, R. N., & Irwin, J. A. 2016, MNRAS, 458, 4210 [NASA ADS] [CrossRef] [Google Scholar]
  35. Hillier, A., & Arregui, I. 2019, ApJ, 885, 101 [NASA ADS] [CrossRef] [Google Scholar]
  36. Hong, W.-S., Zhu, W., Wang, T.-R., Yang, X., & Feng, L.-L. 2024, MNRAS, 529, 4262 [NASA ADS] [CrossRef] [Google Scholar]
  37. Kazantsev, A. P., Ruzmaikin, A. A., & Sokolov, D. D. 1985, Zhurnal Eksperimentalnoi i Teoreticheskoi Fiziki, 88, 487 [NASA ADS] [Google Scholar]
  38. Kereš, D., Katz, N., Weinberg, D. H., & Davé, R. 2005, MNRAS, 363, 2 [Google Scholar]
  39. Komatsu, E., & Seljak, U. 2001, MNRAS, 327, 1353 [Google Scholar]
  40. Ledos, N., Takasao, S., & Nagamine, K. 2024, MNRAS, 527, 11304 [Google Scholar]
  41. Lu, Y. S., Mandelker, N., Oh, S. P., et al. 2024, MNRAS, 527, 11256 [Google Scholar]
  42. Ludlow, A. D., Bose, S., Angulo, R. E., et al. 2016, MNRAS, 460, 1214 [Google Scholar]
  43. Mandelker, N., Padnos, D., Dekel, A., et al. 2016, MNRAS, 463, 3921 [NASA ADS] [CrossRef] [Google Scholar]
  44. Mandelker, N., van Dokkum, P. G., Brodie, J. P., van den Bosch, F. C., & Ceverino, D. 2018, ApJ, 861, 148 [NASA ADS] [CrossRef] [Google Scholar]
  45. Mandelker, N., Nagai, D., Aung, H., et al. 2019, MNRAS, 484, 1100 [Google Scholar]
  46. Mandelker, N., Nagai, D., Aung, H., et al. 2020a, MNRAS, 494, 2641 [Google Scholar]
  47. Mandelker, N., Van Den Bosch, F. C., Nagai, D., et al. 2020b, MNRAS, 498, 2415 [NASA ADS] [CrossRef] [Google Scholar]
  48. Mao, S. A., Carilli, C., Gaensler, B. M., et al. 2017, Nat. Astron., 1, 621 [NASA ADS] [CrossRef] [Google Scholar]
  49. Martin, C. L., Shapley, A. E., Coil, A. L., et al. 2012, ApJ, 760, 127 [Google Scholar]
  50. Neistein, E., & Dekel, A. 2008, MNRAS, 383, 615 [Google Scholar]
  51. Nixon, C. J., Hands, T. O., King, A. R., & Pringle, J. E. 2018, MNRAS, 477, 3539 [CrossRef] [Google Scholar]
  52. Oku, Y., Tomida, K., Nagamine, K., Shimizu, I., & Cen, R. 2022, ApJS, 262, 9 [NASA ADS] [CrossRef] [Google Scholar]
  53. Pakmor, R., Bieri, R., van de Voort, F., et al. 2024, MNRAS, 528, 2308 [NASA ADS] [CrossRef] [Google Scholar]
  54. Prada, F., Klypin, A. A., Cuesta, A. J., Betancort-Rijo, J. E., & Primack, J. 2012, MNRAS, 423, 3018 [NASA ADS] [CrossRef] [Google Scholar]
  55. Putman, M. E., Peek, J. E. G., & Joung, M. R. 2012, ARA&A, 50, 491 [Google Scholar]
  56. Rosdahl, J., & Blaizot, J. 2012, MNRAS, 423, 344 [Google Scholar]
  57. Ryu, D., Jones, T. W., & Frank, A. 2000, ApJ, 545, 475 [Google Scholar]
  58. Salvesen, G., Beckwith, K., Simon, J. B., O’Neill, S. M., & Begelman, M. C. 2014, MNRAS, 438, 1355 [NASA ADS] [CrossRef] [Google Scholar]
  59. Schekochihin, A. A., Cowley, S. C., Taylor, S. F., Maron, J. L., & McWilliams, J. C. 2004, ApJ, 612, 276 [NASA ADS] [CrossRef] [Google Scholar]
  60. Scoville, N., Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1 [Google Scholar]
  61. Spruit, H. C. 2013, ArXiv e-prints [arXiv:1301.5572] [Google Scholar]
  62. Stewart, K. R., Maller, A. H., Oñorbe, J., et al. 2017, ApJ, 843, 47 [NASA ADS] [CrossRef] [Google Scholar]
  63. Vazza, F., Brunetti, G., Brüggen, M., & Bonafede, A. 2018, MNRAS, 474, 1672 [NASA ADS] [CrossRef] [Google Scholar]
  64. Vossberg, A. C. E., Cantalupo, S., & Pezzulli, G. 2019, MNRAS, 489, 2130 [NASA ADS] [CrossRef] [Google Scholar]
  65. Xu, H., Li, H., Collins, D. C., Li, S., & Norman, M. L. 2009, ApJ, 698, L14 [NASA ADS] [CrossRef] [Google Scholar]
  66. Yang, L., Roberts-Borsani, G., Treu, T., et al. 2021, MNRAS, 501, 1028 [Google Scholar]
  67. Zhang, S., Cai, Z., Xu, D., et al. 2023, Science, 380, 494 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Power spectra

We first present the general evolution of the power spectra of the kinetic energy, the magnetic energy, and the enstrophy, defined as ω2 with ω the magnitude of the vorticity. We then quantify the similarities between the functional form of the magnetic energy and the enstrophy power spectra.

A.1. General analysis

To confirm that the magnetic field is amplified by stretching and tangling, we compare the energy power spectra7 of the magnetic energy e ^ m $ \widehat{e}_{\mathrm{m}} $, kinetic energy e ^ k $ \widehat{e}_{\mathrm{k}} $, and enstrophy ω 2 ^ $ \widehat{\omega^2} $. The time evolution of the power spectra for the kinetic energy, magnetic energy, and enstrophy is shown in Fig. A.1 for a 2D and a 3D simulation. All spectra are normalised to the initial total kinetic energy E ^ k , 0 = e ^ k ( t = 0 ) d k $ \widehat{E}_{\mathrm{k,0}}= \int \widehat{e}_{\mathrm{k}}(t=0) \,\mathrm{d}k $. The evolution of the spectra is relatively similar between 2D and 3D, although the energy spectra of the 3D simulation have a steeper slope than in 2D. The 2D spectra display artificial peaks at high wavenumber k (green dashed lines for t ≥ 12t0), corresponding exactly to the resolution limits of the simulations.8 The 3D simulation does not exhibit any artificial peaks. This can be attributed to the following two points: (1) the higher number of cells per wavenumber bin when computing the power spectra in 3D compared to 2D; and (2) the lower resolution of the snapshots used to compute the 3D spectra.

thumbnail Fig. A.1.

Power spectra of the kinetic energy, the enstrophy, and the magnetic energy normalised by the total initial kinetic energy. Results are shown for both 2D and 3D MHD simulations with ℳ = 1, δ = 30, ns = 0.01 cm−3. The spectra are plotted at different times in each panel. From left to right, the vertical black lines refer to length-scales of Rs, Rs/16, Rs/32, and Rs/64, the minimum cell size. For illustration, the power laws k3/2, k−5/3 and k−3 are shown by grey dotted, plain and dashed lines, respectively.

All spectra exhibit an energy injection scale, an inertial range, and a dissipation range. The magnetic energy spectrum shows a typical inverse cascade for scales ≳Rs, analogous to the k ∝ 3/2 slope predicted by various theoretical models (e.g. Kazantsev et al. 1985; Schekochihin et al. 2004), and as observed in galaxy and cluster simulations (e.g. Xu et al. 2009; Vazza et al. 2018; Pakmor et al. 2024). A near equipartition9 is first achieved at high wavenumbers, as the magnetic field amplification occurs initially on small scales.In agreement with previous numerical studies on MHD turbulence (e.g. Ryu et al. 2000; Cho et al. 2009; Salvesen et al. 2014), the equipartition extends to larger scales over time, resulting in similar slopes for the magnetic and kinetic energy where equipartition is established. The resulting slopes of the kinetic and magnetic energy at high wavenumbers lie between k−5/3 and k−3.

A detailed characterisation of the slope and its origin through a spectral energy transfer function analysis is beyond the scope of this current work. However, to support that stretching and tangling drive the magnetic field amplification, we present a phenomenological argument: at early times, the magnetic energy spectrum is proportional to the enstrophy spectrum. Initially, the magnetic field is uniform and lacks a defined spectrum, but it immediately adopts the shape of the enstrophy spectrum, corresponding to the characteristic timescale of magnetic field amplification. The magnetic energy spectrum closely follows the shape of the enstrophy spectrum for length scales ≳0.1Rs, which also corresponds to the length scales of the mixing layer, where tangling and stretching of the field lines occur. We verify this similarity by quantifying the difference between the magnetic energy and enstrophy spectra in the next section.

A.2. Quantification of the power spectrum difference

To quantify the difference between the enstrophy and magnetic energy spectrum, both spectra are normalised by their maximum value. From the normalised spectrum Em, norm(k) and ωm, norm2(k), we compute the difference of their logarithmic value,

ε = 1 N N ( log 10 ( E m , norm ( k ) ω m , norm 2 ( k ) ) ) , $$ \begin{aligned} \varepsilon = \frac{1}{N}\sum _{\mathrm{N}}\left( \log _{\mathrm{10}}\left(\frac{E_{\mathrm{m,norm}}(k)}{\omega ^2_{\mathrm{m,norm}}(k)} \right) \right), \end{aligned} $$(A.1)

where N is the number of wavenumber on which the average is done. As we are mainly interested in the short wavenumber, the average is performed for k < 100. Figure A.2 shows the normalised spectrum at the wavenumbers of interest from one simulation at different times and the time profile of the deviation ε for different simulations. From the normalised spectrum in Fig. A.2, it appears that the magnetic energy spectrum is nearly identical to the enstrophy one at the early stage of the simulation. The shapes of the spectra slightly deviates from one another towards the end of the simulation. The deviation ε is shown for different times where it appears that ε ∼ 0.3 can still reflect a relatively good agreement between the normalised spectra. In the right panel of Fig. A.2, all simulations exhibit a deviation ε ≲ 0.3 for t ≲ 7t0, which correlates with the growth of the magnetic field being the strongest also for the early stage of the simulation. This small value of ε suggests that the growth of the magnetic field is driven by tangling and stretching of the field lines at length-scales of size l ≳ 0.1Rs.

thumbnail Fig. A.2.

Difference between the enstrophy and magnetic energy spectrum. Left: Illustrative plot of the enstrophy and magnetic energy power spectra normalised by their maximum values at different time (t0, t1, t2)∼(0.02, 0.1, 1)tend. For clarity, the spectra are shifted on the y-axis, and their corresponding deviation ε is shown. The results shown here are for a 2D MHD+TC simulation with ℳ = 1, δ = 100, ns = 0.01 m−3. Right: Time profile of the difference ε for different Mach number ℳ, and density ratio δ = 100. The results shown here are for 2D MHD+TC simulations with ns = 0.01 m−3, and a 3D MHD simulation (black dotted line) with ℳ = 1, δ = 30, ns = 0.01 m−3.

Appendix B: Evolution of the stream in 2D and 3D

The mechanisms leading to the magnetic field growth in 2D and 3D simulations are the same. However, slight differences appear in the efficiency of the mixing. Figure B.1 shows the time evolution of the stream density and magnetic field strength slices through the stream midplane for both 2D and 3D simulations. Initially, the magnetic field is amplified at the CGM and stream interface where the velocity shear and the KHI are most pronounced. Consequently, at t = 4t0 ∼ 0.2tv, in both 2D and 3D simulations, the magnetic field near the centre of the stream remains at its initial value, while it undergoes significant amplification at the interface, reaching up to ≳200 times its initial strength. By t = 11 t0 ∼ 0.6tv, the magnetisation of the stream differs between 2D and 3D. In 3D, the mixing is more efficient due to the KHI (see Mandelker et al. 2019), leading to a uniform magnetic field distribution in the stream, while in 2D, the magnetic field near the stream’s centre remains unchanged.

thumbnail Fig. B.1.

Comparison of magnetic field (top) and density (bottom) in 2D and 3D MHD simulations with ℳ = 1, δ = 30, ns = 0.01 cm−3 at the initial time, at t ∼ 4 t0, and at half simulation time t = 11 t0. Each panel shows the 32 kpc full stream axis length in the horizontal direction and a zoomed region of 10 kpc in the vertical direction. The 3D maps represent slices of the density and magnetic field and not projected quantities.

Appendix C: Derivation of the model with a simple outflow term

We hereby derive the model assuming that the galaxy experiences a constant outflow from supernovae feedback.

Starting from Eq. 7, we focus on the inflow and outflow components Em, in and Em, out, leaving all other terms equal to zero. We can then rewrite Ėm,g as

d E m , g d t = d E m , in d t d E m , out d t = d E m , in d M in d M in d t d E m , out d M out d M out d t . $$ \begin{aligned} \dfrac{\mathrm{d} E_{\mathrm{m,g} }}{\mathrm{d} t} = \dfrac{\mathrm{d} E_{\mathrm{m,in} }}{\mathrm{d} t} - \dfrac{\mathrm{d} E_{\mathrm{m,out} }}{\mathrm{d} t} = \dfrac{\mathrm{d} E_{\mathrm{m,in} }}{\mathrm{d} M_{\mathrm{in} }}\dfrac{\mathrm{d} M_{\mathrm{in} }}{\mathrm{d} t} - \dfrac{\mathrm{d} E_{\mathrm{m,out} }}{\mathrm{d} M_{\mathrm{out} }}\dfrac{\mathrm{d} M_{\mathrm{out} }}{\mathrm{d} t}. \end{aligned} $$(C.1)

Em, in is still given by Eq. 9, and Em, out yields,

E m , out = B g 2 2 V out = E m , g V g M out ρ out , $$ \begin{aligned} E_{\mathrm{m,out} }= \frac{B_{\mathrm{g} }^2}{2} V_{\mathrm{out} } = \frac{E_{\mathrm{m,g} }}{V_{\mathrm{g} }} \frac{M_{\mathrm{out} }}{\rho _{\mathrm{out} }}, \end{aligned} $$(C.2)

where Vout, Mout and ρout are the volume, mass, and density of the gas outflow leaving the galaxy. Vg is the galactic volume such that Em, g = 0.5VgBg2. We end up with a general formulation as

d E m , g d t = B in 2 2 ρ in d M in d t E m , g ρ out V g d M out d t = C in E m , g τ out , $$ \begin{aligned} \dfrac{\mathrm{d} E_{\mathrm{m,g}}}{\mathrm{d} t} = \frac{B_{\mathrm{in} }^2}{2 \rho _{\mathrm{in} }} \dfrac{\mathrm{d} M_{\mathrm{in} }}{\mathrm{d} t} - \frac{E_{\mathrm{m,g} }}{\rho _{\mathrm{out} }V_{\mathrm{g} }}\dfrac{\mathrm{d} M_{\mathrm{out} }}{\mathrm{d} t} = C_{\mathrm{in} } - \frac{E_{\mathrm{m,g} }}{\tau _{\mathrm{out} }}, \end{aligned} $$(C.3)

with, τout representing the magnetic field expulsion timescale. Assuming that τout and Cin are invariant in time, the solution of the first-order differential Eq. C.3 is,

E m , g ( t ) = E m , g , 0 exp ( t τ out ) + τ out C in [ 1 exp ( t τ out ) ] . $$ \begin{aligned} E_{\mathrm{m,g} }(t) = E_{\mathrm{m,g,0} }\exp {\left(-\frac{t}{\tau _{\mathrm{out} }}\right)} + \tau _{\mathrm{out} }C_{\mathrm{in} } \left[1-\exp {\left(-\frac{t}{\tau _{\mathrm{out} }}\right)}\right]. \end{aligned} $$(C.4)

For Cin, Eq. 12 still holds. For τout, we define the outflow density as the average density of the gas flowing out of the galaxy disc surface πRg2 at a velocity vout,

ρ out = 1 π R g 2 v out d M out d t , $$ \begin{aligned} \rho _{\mathrm{out} } = \frac{1}{\pi R_{\mathrm{g} }^2 v_{\mathrm{out} }}\dfrac{\mathrm{d} M_{\mathrm{out} }}{\mathrm{d} t}, \end{aligned} $$(C.5)

which leads to

τ out = 2 d g v out , $$ \begin{aligned} \tau _{\mathrm{out} }= \frac{2 d_{\mathrm{g} }}{v_{\mathrm{out} }}, \end{aligned} $$(C.6)

where dg is the galaxy’s disc height, and we fix vout = 300 km s−1 based on the isolated galaxy simulations of Oku et al. (2022, Fig. 18).

Finally, the mean magnetic field in the galactic volume Vg can be recovered as

B g ( t ) = ( 2 E m , g ( t ) V g ) 1 / 2 . $$ \begin{aligned} B_{\mathrm{g}}(t) = \left(2\frac{E_{\mathrm{m,g} }(t)}{V_{\mathrm{g}}} \right)^{1/2}. \end{aligned} $$(C.7)

From Eq. C.4, the asymptotic value for Em, g is τoutCin which represents the amount of energy brought by inflow during a timescale t = τout. The time for Em, g to reach a value of 0.9τoutCin is given by t ∼ 2.3τout.

All Figures

thumbnail Fig. 1.

Time evolution of the stream and the magnetic field lines for a 3D simulation (with ns = 0.01 cm−3). Left: Maps showing volume rendering of the density field and the magnetic field lines. The lines are coloured by magnetic field strength. Right: Same as the left panel, but with only the magnetic field lines.

In the text
thumbnail Fig. 2.

Time profiles of the magnetic field in the stream for 2D and 3D simulations (dotted and solid lines, respectively). Simulations are shown for ns = 10−2 and 10−3 cm−3 in purple and orange, respectively. The dotted grey line stands for a 2D simulation with ns = 10−2 cm−3 and with a magnetic field parallel to the stream (no amplification due to stretching). The dashed lines show the model results from Eq. (2) with the same Mach number as the simulations (ℳ = 1).

In the text
thumbnail Fig. 3.

Illustrative sketch of the main processes of the model. The cold stream arrives initially from the intergalactic medium (IGM) with a mass flow rate Ms, 0 and a magnetic field B0 ∼ nG. Step 1: As the cold stream enters the halo, a mixing layer forms at the interface of the CGM and the stream. Magnetic field lines are tangled and stretched in the mixing layer, resulting in magnetic field amplification. Simultaneously, CGM gas in the mixing layer condenses onto the stream, accumulating magnetic fields. From the interaction of the stream and the CGM, the cold mass flux and the magnetic field of the stream become F and Bs. Step 2: Upon falling into the halo, the cold stream narrows and is squeezed by the dark matter potential. Consequently, magnetic fields in the stream are also compressed, further amplifying it by a factor α. The stream is also expected to entrain CGM gas while falling through the halo potential, increasing its mass flux by an efficiency factor η. Step 3: Ultimately, if the stream reaches the galaxy, it continuously accretes cold gas at a rate in with an associated magnetic energy inflow rate Ėm,in to the galaxy, elevating the galactic mean magnetic field. By default, we did not consider outflow; however, we introduced a simple outflow term in Appendix C.

In the text
thumbnail Fig. 4.

Cold mass flux rate normalised by the initial cold mass flux of the simulations. Each point represents one simulation with a fixed ξ, ℳ, and initial magnetic angle ϕ with respect to the stream axis. Red dots, blue squares, and cyan diamonds refer to 2D simulations with ℳ = 2, 1, 0.5, respectively. Using the same colour legend, the fitted model is represented by coloured curves. For a fixed ξ and ℳ, the two points represent simulations with an initial magnetic field angle of 90° and 45° with the highest values associated with the 90° angle in all cases.

In the text
thumbnail Fig. 5.

Mean galactic magnetic field from our model as a function of the stream cold mass flux entering the halo s,0. The cold mass flux is obtained by varying its number density ns ∈ [0.001,0.1] cm−3 for different stream radii Rs = (0.3, 3, 30) kpc (plain, dashed and dotted lines) and for different Mach number ℳ = 0.5, 1, 2 (light to dark colours). The red and blue colour maps show results for the default case Bs and the no-stretching case (see Sect. 3.3), respectively. The vertical grey line shows the maximum cold mass flux simulated in the simulations suite.

In the text
thumbnail Fig. 6.

Time evolution of the mean galactic magnetic field of our model. The corresponding redshift is also given at the top axis assuming that the accretion starts at z = 3. The blue and red curves represent the fiducial and lower-limit scenarios. The solid and dotted lines are models without and with outflows. The black line with the arrow shows the target at 500 μG or less from Geach et al. (2023).

In the text
thumbnail Fig. A.1.

Power spectra of the kinetic energy, the enstrophy, and the magnetic energy normalised by the total initial kinetic energy. Results are shown for both 2D and 3D MHD simulations with ℳ = 1, δ = 30, ns = 0.01 cm−3. The spectra are plotted at different times in each panel. From left to right, the vertical black lines refer to length-scales of Rs, Rs/16, Rs/32, and Rs/64, the minimum cell size. For illustration, the power laws k3/2, k−5/3 and k−3 are shown by grey dotted, plain and dashed lines, respectively.

In the text
thumbnail Fig. A.2.

Difference between the enstrophy and magnetic energy spectrum. Left: Illustrative plot of the enstrophy and magnetic energy power spectra normalised by their maximum values at different time (t0, t1, t2)∼(0.02, 0.1, 1)tend. For clarity, the spectra are shifted on the y-axis, and their corresponding deviation ε is shown. The results shown here are for a 2D MHD+TC simulation with ℳ = 1, δ = 100, ns = 0.01 m−3. Right: Time profile of the difference ε for different Mach number ℳ, and density ratio δ = 100. The results shown here are for 2D MHD+TC simulations with ns = 0.01 m−3, and a 3D MHD simulation (black dotted line) with ℳ = 1, δ = 30, ns = 0.01 m−3.

In the text
thumbnail Fig. B.1.

Comparison of magnetic field (top) and density (bottom) in 2D and 3D MHD simulations with ℳ = 1, δ = 30, ns = 0.01 cm−3 at the initial time, at t ∼ 4 t0, and at half simulation time t = 11 t0. Each panel shows the 32 kpc full stream axis length in the horizontal direction and a zoomed region of 10 kpc in the vertical direction. The 3D maps represent slices of the density and magnetic field and not projected quantities.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.