Baryon-Dark matter interaction in presence of magnetic fields in light of EDGES signal

TL;DR

This paper explores how cosmic magnetic fields and baryon-dark matter interactions jointly influence the 21 cm signal observed by EDGES, deriving bounds on magnetic field strength, dark matter mass, and interaction cross-section.

Contribution

It investigates the combined effects of magnetic fields and baryon-dark matter interactions on the EDGES signal, providing new bounds and showing the allowed magnetic field range can increase significantly.

Findings

Strong magnetic fields require larger baryon-dark matter cross-sections to explain EDGES.

Upper bound on magnetic field strength is $3.48\times10^{-6}$ G for dark matter mass $\lesssim 10^{-2}$ GeV.

Allowed primordial magnetic field range can increase by three orders of magnitude.

Abstract

We have shown that in presence of a cosmic magnetic field the bounds on baryon dark matter cross-section ($\hat \sigma$), dark-matter mass ($m_d$) and values of the magnetic field ($B_0$) can strongly influence each other. This requires to rework the bounds on $\hat \sigma\,$, $m_d$ and $B_0$ which can explain the observed absorption signal by EDGES collaboration. The upper limit on the magnetic field strength can modify in presence of baryon-dark matter interaction cross-section. In the presence of a strong magnetic field, a large baryon-dark matter interaction cross-section is required to balance magnetic heating of gas to explain the EDGES signal as compared to a weak magnetic field. Subsequently, the strong magnetic-fields can even erase the 21 cm signal--this gives an upper bound on the strength of magnetic-fields, dark-matter mass and baryon-dark matter cross-section. In the special case when $\hat \sigma$=0, one can recover the bound on magnetic field strength calculated in [1]. In this work we find that the allowed range of the primordial magnetic field can increase by three orders of magnitude in comparison with [1]. We get upper bound on the magnetic field strength: $3.48\times10^{-6}$ G for the dark matter mass $\lesssim 10^{-2}$ GeV.

Equations

$$\frac{n_1}{n_0}=\frac{g_1}{g_0}{e}^{-\frac{E_{21}}{T_S}}\simeq 3\left(1-\frac{E_{21}}{T_S}\right),$$

$$T_S^{-1}=\frac{T_{\mathrm{CMB}}^{-1}+x_c{T}_{\mathrm{gas}}^{-1}+x_{\alpha }{T}_{\mathrm{Ly}\alpha }^{-1}}{1+x_c+x_{\alpha }},$$

$$x_c=\frac{E_{21}}{T_{\mathrm{CMB}}}\frac{C_{10}}{A_{10}},x_{\alpha }=\frac{E_{21}}{T_{\mathrm{CMB}}}\frac{P_{01}^w}{A_{10}},$$

$$T_{21}=\frac{1}{1+z}(T_S-T_{\mathrm{CMB}})(1-{\exp }^{-\tau }),$$

$$T_{21}^{\mathrm{obs}}(z=17)=-500_{-500}^{+200}\mathrm{mK},$$

$$⟨{\tilde{\mathbf{B}}}_i(\mathbf{k}){\tilde{\mathbf{B}}}_j^{\ast }(\mathbf{q})⟩=\frac{(2\pi {)}^3}{2}{\delta }_D^3(\mathbf{k}-\mathbf{q})\left({\delta }_{ij}-\frac{k_i{k}_j}{k^2}\right){\mathcal{P}}_B(k),$$

$$\frac{k_{\max }}{2\pi {\mathrm{Mpc}}^{-1}}={\left[1.32\times 10^{-3}{\left(\frac{B_0}{1\mathrm{nG}}\right)}^2{\left(\frac{{\Omega }_b{h}^2}{0.02}\right)}^{-1}\left(\frac{{\Omega }_m{h}^2}{0.15}\right)\right]}^{\frac{-1}{(n_B+5)}},$$

$${\Gamma }_{\mathrm{ambi}}$$

$${\Gamma }_{\mathrm{decay}}$$

$$\gamma =\frac{\frac{1}{2}{n}_H⟨\sigma v{⟩}_{{\mathrm{H}}^{+},\mathrm{H}}+\frac{4}{5}{n}_{\mathrm{He}}⟨\sigma v{⟩}_{{\mathrm{H}}^{+},\mathrm{He}}}{m_H[n_H+4n_{\mathrm{He}}]},$$

$${\Gamma }_{\mathrm{ambi}}\sim \frac{{\rho }_n}{16{\pi }^2\gamma {\rho }_b^2{\rho }_i}\frac{B^4}{L^2},$$

$$\frac{d}{dt}\left(\frac{|\mathbf{B}{|}^2}{8\pi }\right)=-4H(t)\left(\frac{|\mathbf{B}{|}^2}{8\pi }\right)-{\Gamma }_{\mathrm{ambi}}-{\Gamma }_{\mathrm{decay}}.$$

$$\frac{d{T}_{\mathrm{gas}}}{dz}$$

$$\frac{d{T}_d}{dz}$$

$$\frac{dv}{dz}$$

$${\Gamma }_C=\frac{8{\sigma }_T{a}_r{T}_{\mathrm{CMB}}^4{x}_e}{3(1+x_e+x_{He})m_e}.$$

$$D(v)\equiv \frac{{\rho }_m\hat{\sigma }}{m_H+m_d}\frac{1}{v^2}F(r),$$

$$F(r)\equiv \mathrm{erf}(\frac{r}{\sqrt{2}})-\sqrt{\frac{2}{\pi }}r{e}^{-r^2/2},$$

$$\frac{d{Q_{gas}}}{dt}=\frac{2m_{H}\rho_{d}\hat{\sigma}e^{-r^{2}/2}}{\sqrt{2\pi}(m_{H}+m_{d})^{2}u_{\mathrm{th}}^{3}}\bigg{(}T_{\rm d}-T_{gas}\bigg{)}+\frac{\mu\rho_{d}}{\rho_{M}}vD(v)~{},$$

$$\frac{d{x}_e}{dz}$$

Full text

Magnetic fields, 21 cm signal, Baryon-dark matter interaction, EDGES signal

I Introduction

Within the standard cold dark matter cosmology ($\Lambda$CDM), free electrons and protons cool sufficiently after $3\times 10^{5}$ years of big-bang to form neutral hydrogen atoms. At the end of recombination ($z\approx 1100$), matter decouples from the cosmic microwave background (CMB) photons and the temperature reaches around 3000 K. After that, the dark age begins and the Universe becomes homogeneous. Later during the cosmic dawn era ($15<z<35$), over density grow in the matter perturbations and collapse to form the first star and galaxy in the Universe. During the dark ages, residual electrons from the recombination scatter off the baryons to maintain the thermal equilibrium until $z\approx 200$. Subsequently, gas cools adiabatically due to the thermal expansion of the Universe and the gas temperature reaches to $\approx 6.8$ K at $z=17$.

During the cosmic dawn era, the gas temperature is lower than the CMB temperature hence hyperfine transitions in the neutral hydrogen atoms produce 21 cm absorption spectra. The hyperfine transitions are due to the CMB photons, gas collisions and the Ly-$\alpha$ radiations from the first star. 21 cm absorption line leaves an imprint of spectral distortion in the low-frequency tail of the CMB spectrum. The observation of the 21 cm line can give logical reason behind the density fluctuations Hogan and Rees (1979), cosmic reionization Scott and Rees (1990) and X-ray heating of the Inter-Galactic Medium (IGM) Fialkov et al. (2014). Recently, “Experiment to Detect the Global Epoch of Reionization Signature” (EDGES) collaboration reported the first detection of such absorption signal centered at 78 MHz Bowman et al. (2018). The observed absorption dip at $z\approx 17$ is approximately $2.5$ times larger than the standard $\Lambda$CDM model prediction Bowman et al. (2018). Several attempts have been made in the literature to explain the observed EDGES anomaly.

There are two ways to explain this enhanced 21 cm absorption signal: one is to the heating of CMB photons and another way is to cool the gas in IGM below the standard $\Lambda$CDM prediction. The first possibility has been investigated by authors of the Ref. Moroi et al. (2018); Fraser et al. (2018); Pospelov et al. (2018); Liu et al. (2019). In the second scenario, IGM gas can be cooled by emitting the photons between the Ly-limit to Ly-$\gamma$ wavelengths Chuzhoy and Shapiro (2007, 2006). There are very few mechanisms to cool the gas. Since the dark-matter is colder than the gas, effective cooling of the gas can be obtained by elastic scattering between the dark-matter and baryon particles Barkana et al. (2018); Barkana (2018); Tashiro et al. (2014). Therefore, a non-standard Coulomb interaction between dark-matter and baryon: $\sigma=\hat{\sigma}v^{-4}$ can be considered to explain the EDGES signal, where $v$ is the relative velocity between the dark matter and baryon Bowman et al. (2018); Tashiro et al. (2014); Dvorkin et al. (2014); Barkana et al. (2018). In these cases, dark-matter considered as millicharge-dark-matter Bowman et al. (2018); Barkana (2018). Using this mechanism 21 cm absorption signal can be explained Barkana (2018); Barkana et al. (2018); Bowman et al. (2018); Tashiro et al. (2014). In these scenarios, cooling of the gas, by transferring energy to the dark-matter, is tightly constrained because of constraints on the dark-matter mass and cross-section by cosmological and astrophysical observations Barkana et al. (2018); Barkana (2018); Berlin et al. (2018); Creque-Sarbinowski et al. (2019).

It is observed that the magnetic fields (MFs) are present on all length-scales (i.e. at the length-scales of galaxies and the clusters). These MFs have a strength of the order of few $\mu$G in the intergalactic medium Kronberg (1994); Neronov and Vovk (2010). The strength of these MFs is constrained from the structure formation, big bang nucleosynthesis (BBN), temperature anisotropies and polarization of CMB Trivedi et al. (2012, 2014); Sethi and Subramanian (2005). Further, individual limits of the order of $\lesssim$ nG on primordial magnetic fields (PMFs) for various cases has been reported by Planck 2015 results Ade et al. (2016). From the data of Fermi and High Energy Stereoscopic System (HESS) gamma-ray telescopes, the peak strength of the magnetic field at a length scale of 1 Mpc can be the order of few nG Neronov and Vovk (2010). However, the upper bound on the amplitude of the magnetic field, obtained from BBN is $\sim 10^{-6}$ G at a comoving length scale of $\sim 100$ pc Cheng et al. (1996); Grasso and Rubinstein (2001). If strong MFs were present at the time of nucleosynthesis, the abundance of relic 4He and other light elements could have different values than we observe today Matese and O’Connell (1969); GREENSTEIN (1969). To satisfy the current observational limits on the light element abundances, the strength of these MFs can have a value of the order of $\approx 10^{-7}$ G, at present time Tashiro and Sugiyama (2006a); Sethi and Subramanian (2005).

In the presence of PMFs, decaying magneto-hydrodynamics effects can heat the IGM gas Sethi and Subramanian (2005); Schleicher et al. (2008); Chluba et al. (2015). Hence, the thermal evolution of the gas can modify and it can erase the 21 cm absorption signal Minoda et al. (2019). Therefore, magnetic heating of the IGM gas can be constrained by the EDGES signal. Energy deposition into the IGM can be done by ambipolar diffusion. To dissipate the energy into IGM gas, for small-scales, can be done by producing decaying magneto-hydrodynamics turbulence Sethi and Subramanian (2005); Sethi et al. (2008). Heating due to the turbulent-decay starts after recombination when radiative viscosity decreases and velocity perturbations are no longer damped Sethi and Subramanian (2005); Schleicher et al. (2008). Energy dissipation due to the ambipolar diffusion is important in neutral IGM gas Sethi and Subramanian (2005), for details see Shu (1992). At the late time, heating of the gas due to the turbulent-decay decreases, but due to the ambipolar diffusion it continues Chluba et al. (2015); Sethi and Subramanian (2005). In the context of EDGES anomaly, this magnetic heating of the IGM gas has been studied by authors of the Ref. Minoda et al. (2019). By the constraint on gas temperature during redshift $15\lesssim z\lesssim 20$ (EDGES result), they put a constraint on the upper limit of the PMFs strength: $B\lesssim 10^{-10}$ Gauss at the length-scale of 1 Mpc.

If one invokes cooling of gas beyond the standard scenario in the presence of baryon dark matter interaction, discussed above, this bound on the strength of PMFs (obtained in the Ref. Minoda et al. (2019)) may change by transferring energy of the gas to the dark matter (DM) using drag between gas and DM. This transfer of energy depends on the dark-matter mass and cross-section. Therefore, we also get the constraint on DM mass and cross-section to get the 21 cm absorption signal.

There are also several other ways to heat the gas: Heating of the neutral hydrogen in the IGM due to Ly-$\alpha$ photons has been investigated in some Refs. Oklopčić and Hirata (2013); Chuzhoy et al. (2007); Chuzhoy and Shapiro (2007); Ghara and Mellema (2019), annihilation of DM can inject the energy into IGM Slatyer (2016a, b); D’Amico et al. (2018); Liu and Slatyer (2018), DM decay can also heat the IGM gas via energy deposition Mitridate and Podo (2018); Slatyer and Wu (2017); Liu and Slatyer (2018). In Ref. Bhatt et al. (2019), authors have considered the dark matter viscosity to heat the gas and dark matter in the context of observed 21 cm signal Sethi and Subramanian (2005); Schleicher et al. (2008); Chluba et al. (2015). Here it is to be noted that, we are interested only in the magnetic heating of the gas in the presence of DM-baryon interaction due to the effective cooling of IGM by DM. The origin of these PMFs fields could be due to some high energy process in the very early universe Turner and Widrow (1988); Sharma et al. (2018); Bhatt and Pandey (2016); Anand et al. (2017); Subramanian (2016).

The present work is divided into the following sections: in section (II), we have revisited the 21 cm observed signal from the EDGES; in section (III), a brief description of the PMFs and it’s decay is discussed; section (IV) contains the result obtained and a detailed discussion. In the end, we have summarized and concluded our work in section (VI).

II 21 cm signal and EDGES observation

At the end of recombination, the baryon number density of the Universe mostly dominated by the neutral hydrogen, small fraction of helium, residual free electrons and protons. The hyperfine interaction between spin of the electron and proton split the ground state of neutral hydrogen atom into singlet and triplet states with an energy difference of $E_{21}=5.9\times 10^{-6}$ eV $=2\pi/(21~{}\text{cm})$. Relative number densities of neutral hydrogen in the singlet and triplet state, define the spin temperature ($T_{S}$) of the gas, and is given by the following relation:

[TABLE]

here, $n_{0}$ and $n_{1}$ are the number densities in the singlet and triplet states respectively. $g_{0}$ and $g_{1}$ are statistical degeneracy of singlet and triplet states respectively. Spin temperature depends on gas collision, emission/absorption of CMB photon and Ly-$\alpha$ radiation from the first star. Equilibrium balance between the populations of singlet and triplet state describes the spin temperature Field (1958); Pritchard and Loeb (2012)

[TABLE]

where, $T_{\rm CMB}$ is the CMB temperature, $T_{\rm gas}$ is the kinetic temperature of the gas and $T_{\rm Ly\alpha}$ is the colour temperature of Ly$\alpha$ radiation from the first star. Here, $x_{c}$ and $x_{\alpha}$ are the collisional and Ly$\alpha$ coupling respectively Field (1958),

[TABLE]

where, $C_{10}$ is the collision rate between $H-H,~{}H-e,~{}p-H$ and $P_{01}^{w}$ is the excitation rate due to Ly-$\alpha$ radiation and $A_{10}=2.9\times 10^{-15}$sec*-1* is the Einstein coefficient for spontaneous emission. After the first star formation, a large number of Ly$\alpha$ photons scattered with the gas, and brought the radiation and the gas into a local thermal equilibrium Pritchard and Loeb (2012). Hence during the cosmic dawn era $T_{\rm gas}\approx T_{\rm Ly\alpha}$. The 21 cm signal can be described, in terms of the redshifted differential brightness temperature Pritchard and Loeb (2012)

[TABLE]

where, optical depth $\tau=\frac{3\lambda_{21}^{2}A_{10}n_{H}}{16T_{S}H(z)}$, $n_{H}$ is the hydrogen number density, $\lambda_{21}=21$ cm and $H(z)$ is the Hubble rate. Depending on the spin and the CMB temperature, three scenarios arise: (i) when $T_{S}~{}=T_{\rm CMB}$, then no signal is observed; (ii) if $T_{S}<T_{\rm CMB}$, photon get absorbed by the gas and absorption spectra is observed; and (iii) if $T_{S}>T_{\rm CMB}$, then it leaves an imprint of emission spectra.

Evolution of the 21 cm signal is as follows: after recombination at $z\sim 1100$ down to $z\sim 200$, the residual free electrons undergo Compton scattering to maintain thermal equilibrium between the gas and CMB, and collisions among the gas is dominant, i.e. $x_{c}\gg 1,x_{\alpha}$ Barkana et al. (2018); Pritchard and Loeb (2012), which set $T_{S}=T_{\rm CMB}$. Hence, 21 cm signal is not observed during this era. Below $z\sim 200$ until $z\sim 40$, the gas cools adiabatically and it falls below CMB temperature, which implies the early 21 cm absorption signal. The sensitivity of radio antennas below 50 Mhz falls dramatically and collisional absorption signal can not be observed. Below $z\sim 40$ down to the first star formation, gas cools sufficiently due to the expansion and the collisional coupling becomes very small due to the dilution, i.e. $x_{c},~{}x_{\alpha}\rightarrow 0$ Barkana et al. (2018); Pritchard and Loeb (2012). This implies, $T_{S}=T_{\rm CMB}$, hence no $21$ cm signal during this period. After the first star formation, a transition between the singlet and triplet states occurs due to the Ly$\alpha$ photons emitted from the first stars via Wouthuysen-Field (WF) effect Wouthuysen (1952); Field (1959). Ly$\alpha$ photons couple the spin temperature to the gas temperature. In this era, $x_{\alpha}\gg 1$, hence spin temperature and the gas temperature become equal to each other i.e., $T_{S}=T_{\rm gas}<T_{\rm CMB}$. Thus an imprint of the 21 cm absorption signal can be seen at the low-frequency tail of the CMB spectrum. Below $z\sim 15$ until $z\sim 7$, X-ray from the active galactic nuclei heats the gas above the CMB temperature and we observe an emission signal Pritchard and Loeb (2012). Below $z\sim 7$ neutral hydrogen became ionized and the signal disappears.

Recently, the EDGES collaboration reported the global brightness temperature at $z=17$

[TABLE]

and corresponding gas temperature is $3.26^{+1.94}_{-1.58}$ K Bowman et al. (2018). On the contrary, standard $\Lambda$CDM predicts the gas temperature 6.8 K at $z=17$ and corresponding brightness temperature $T_{21}\geq-220$ mK Barkana et al. (2018). In order to explain the EDGES absorption signal, the gas temperature needs to be cooler than the $\Lambda$CDM prediction. During the Cosmic dawn era, the Universe was at its coldest phase, and the relative velocity between the DM and baryon was very small [$\mathcal{O}(10^{-6})$]. Also, the temperature of the dark matter is colder than the baryon temperature during this period, so an interaction of the baryon with dark matter can cool the gas temperature. Since the relative velocity is small, scattering cross section of the type $\sigma=\hat{\sigma}v^{-4}$ can enhance the interaction rate and cool the gas sufficiently to explain EDGES absorption dip Tashiro et al. (2014); Dvorkin et al. (2014). In this work, we consider magnetic heating of the gas and the DM via ambipolar and turbulent decay.

III Primordial magnetic fields after the recombination era

In this section, we present the evolution of PMFs via two processes i.e. ambipolar diffusion and decaying turbulent. We assume that, due to some early Universe process, tangled magnetic fields were present at a sufficiently large length scale after the recombination era Turner and Widrow (1988); Sharma et al. (2018); Bhatt and Pandey (2016); Anand et al. (2017); Subramanian (2016). We also consider a very small velocity induced by the PMFs to avoid any dissipations of the initial magnetic fields due to viscosity and other dissipation Jedamzik et al. (1998); Subramanian and Barrow (1998). This is applicable in linear regime and magnetic fields evolve adiabatically as ${\bf B}(t,{\bf x})=\tilde{{\bf B}}({\bf x})/a(t)^{2}$, where ${\bf x}$ is the comoving coordinate, $\tilde{{\bf B}}$ is the comoving strength of the magnetic fields and $a(t)$ is the scale factor. Since plasma in the early Universe remains highly conductive, the adiabatic evolution of the magnetic fields is true. However, at a sufficiently small scale, when non-linear effects operate, adiabatic decay no longer satisfy. In this case, it is needed to consider Euler Eqs. along with the magnetic induction equation and Maxwell’s Eqs. to understand the dynamics of the fluid. We consider an isotropic and homogeneous Gaussian random magnetic field, whose power spectrum is given by the following equation

[TABLE]

where, $\mathcal{P}_{B}(k)$ is the magnetic power spectrum and $k=|{\bf k}|$ is the comoving wave number. For simplicity, we consider a power law spectrum of the magnetic fields $\mathcal{P}_{B}(k)=A\,k^{n_{B}}$ for $k<k_{\rm max}$ ($k_{\rm max}$ is calculated by the damping at recombination due to the viscosity) Jedamzik et al. (1998); Subramanian and Barrow (1998). Here $n_{B}$ and $A$ are spectral index and the normalization constants respectively. In particular, $n_{B}=2$ for white noise Hogan (1983), $n_{B}=4$ for Batchelor spectrum Durrer and Caprini (2003) and $n_{B}=-2.9$ for nearly scale invariant spectrum Sethi and Subramanian (2005). The amplitude ($A$) can be obtained by demanding the magnetic fields are smooth over the cut off scale, and after that $\mathcal{P}_{B}(k)=0$. Once the recombination period end, baryons and CMB photons decouple and their velocity start to increase. Eventually, these particles achieve a common velocity, determined by the equipartition between the magnetic field and kinetic energy of the baryon gas. This velocity is given by the Alfv́en velocity $v_{A}=c\,B_{0}/\sqrt{4\pi\rho_{b_{0}}a(t)}$, here $\rho_{b_{0}}$ is the present baryon density and $B_{0}$ is the currently observed strength of the magnetic field at Mpc scale. The cut off scale is defined as $k_{{\rm max}}\approx 2\pi\frac{Ha}{v_{A}}$. Therefore, the cut off value from this relation can be written as Sethi and Subramanian (2005); Tashiro and Sugiyama (2006b):

[TABLE]

here, $\Omega_{b}$, $\Omega_{m}$, and $h$ are the cosmological parameters and have the mathematical value $h=0.674$ ($H_{0}=100\,h$ Km/s/Mpc), $h^{2}\Omega_{b}=0.0224\pm 0.0001$, $h^{2}\Omega_{m}=0.143$ Aghanim et al. (2018).

Dissipation of magnetic energy

After recombination, any present magnetic field dissipates its energy and heats the gas through two mechanisms, namely the ambipolar diffusion and the turbulent decay Sethi and Subramanian (2005); Tashiro and Sugiyama (2006b). The heating of the gas gives a considerable change in the thermal evolution of neutral atoms. The velocity difference in the ionized and neutral particles after the recombination leads to the ambipolar diffusion of the magnetic energy. Direct cascade happens due to the non-linear processes, which couples the different modes and the cascading of magnetic energy from large to small scale. This happens through the breaking of the larger eddy into the smaller eddies, when eddy turn over time $t_{\rm eddy}$ is equal to the Hubble time i.e. ($t_{\rm eddy}\sim H^{-1}$). The energy dissipation due to ambipolar and direct cascade can be given by Cowling (1956); Shu (1992),

[TABLE]

where, $\rho_{n}$ and $\rho_{i}$ are the mass densities of neutral and the ionized atoms respectively, $t$ is the cosmological time at a generic red shift $z$, $t_{\rm eddy}$ is the physical decay time scale for the turbulent and $t_{i}$ is initial time at which decay starts. $m=2(n_{B}+3)/(n_{B}+5)$. For the present scenario $\gamma$ is given by Shu (1992); Shang et al. (2002); Draine (1980); Schleicher et al. (2008),

[TABLE]

where, $m_{H}$ is the mass of the hydrogen and $n_{\rm He}$ is the number density of the helium atom. The absolute values of the Lorentz force and magnetic energy in Eqs. (7) and (8) can be obtained by considering a suitable power law spectrum of the magnetic field. This can be done using following correlation integrals: $|(\nabla\times{\bf B})\times{\bf B}|^{2}=\int(dk/2\pi)^{3}\int(dq/2\pi)^{3}\,k^{2}\,\mathcal{P_{B}}(t,k)\mathcal{P_{B}}(t,q)(1+z)^{10}$ and $|{\bf B}|^{2}=\int(dk/2\pi)^{3}\mathcal{P_{B}}(t,k)(1+z)^{4}$. However, in the present work, we have taken approximate value of the ambipolar diffusion term, which is Schleicher et al. (2009)

[TABLE]

here, $L$ is a typical length scale. To sum up, the time evolution of the magnetic energy can be written as Subramanian and Barrow (1998); Sethi and Subramanian (2005)

[TABLE]

IV Baryon and dark matter interaction in presence of magnetic field

In this section, we discuss the effects of magnetic field on the baryon and DM temperature when they are interacting with each other. Temperature evolutions of the DM and baryon, having relative velocity, in the presence of the magnetic field are given below

[TABLE]

where, $T_{d}$ and $m_{d}$ are temperature and mass of the DM respectively, $H$ is the Hubble expansion rate, $n_{b}$ is the baryon number density and $\Gamma_{C}$ is the Compton scattering rate, defined as

[TABLE]

Here, $x_{e}=n_{e}/n_{H}$ is the electron fraction, $x_{He}$ is the helium fraction, $a_{r}$ is the Stefan-Boltzmann radiation constant and $\sigma_{T}$ is the Thomson scattering cross section. Drag term $D(v)$ is given by

[TABLE]

where, $\rho_{M}$ is matter density and

[TABLE]

with $r=\frac{v}{u_{\mathrm{th}}}$ and $u_{\mathrm{th}}=\sqrt{\frac{T_{\rm gas}}{m_{H}}+\frac{T_{d}}{m_{d}}}~{}.$ Heat transfer of the gas per unit time is given by Muñoz et al. (2015)

[TABLE]

where, $\mu$ is the reduced mass of DM and baryon, $\hat{\sigma}$ is DM-gas scattering cross-section and $\rho_{d}$ is the DM energy density. In Eq.(16), the first term represents the baryon’s cooling due to its interaction with the DM and second term represents the heating due to drag term. Relative velocity between DM and gas generates friction between two fluid which is responsible for the drag term. Heat transfer rate for dark matter $\left(\frac{d{Q_{d}}}{dt}\right)$ can be obtained by interchanging $gas\leftrightarrow d$ in Eq. (16). Temperature evolution of the DM and gas require electron ionization fraction Ali-Haimoud and Hirata (2011):

[TABLE]

where, $\beta_{B}$ and $\alpha_{B}$ are the photo-ionization rate and case-B recombination coefficient respectively. $E_{21}$ is energy of Ly$\alpha$ wavelength photon and $\Lambda_{2s,1s}=8.22~{}{\rm sec}^{-1}$ is the hydrogen two photon decay rate. The escape rate of Ly$\alpha$ is given by: $R_{\text{Ly}\alpha}=\frac{8\pi H}{3n_{H}(1-x_{e})\lambda_{\text{Ly}\alpha}^{3}}$, $\lambda_{\rm Ly\alpha}$ is the rest wavelength of Ly$\alpha$ photon. As it has been confirmed in Ref. Minoda et al. (2019), that cooling due to effects like Lyman-$\alpha$ emission, Bremsstrahlung and recombination does not have that much effects on the dynamics of the gas and DM, therefore, we have not considered these effects in our present work.

V Results and Discussion

Solving Eqs. (11-14) and Eq. (17) with initial conditions $T_{\rm gas}(1010)=T_{\rm CMB}(1010)=2749.92$ K, $T_{d}(1010)\sim 0$ K, $x_{e}(1010)=0.057$ and $B(z)=B_{0}\,(1+z)^{2}|_{z=1010}$ initial magnetic field strength, we get the temperature evolution of the DM and gas for different DM masses, DM-baryon interaction cross-sections and MF’s strengths. Figure (1) shows the evolution of the gas and DM temperature with redshift ($z$). The solid blue line in Fig. (1) corresponds to gas temperature when both the magnetic field and DM-baryon interaction are zero. In this case, gas temperature falls as $T_{\rm gas}\propto(1+z)^{2}$ and reaches 6.8 K at $z=17$. In figure (1a), temperature evolution of the gas and DM is given for different MFs at constant $\hat{\sigma}=10^{-41}~{}{\rm cm}^{2}$ and $m_{d}=0.1$ GeV. For $B_{0}=10^{-5}~{}(10^{-6})$ G, gas temperature falls down due to Hubble expansion and DM-baryon interaction till $z\approx 30~{}(z\approx 20)$, then temperature rises due to magnetic heating. For $B_{0}=10^{-5}$ G, temperature of the DM also increases due to the coupling between DM and baryons at lower redshift. Larger the strength of MFs, earlier the heating begins. Although dark matter temperature at $z\sim 1010$ is taken to be zero, it heats up due to the heat transfer from baryons to DM. By increasing $B_{0}$, magnetic-heating of the gas rises. Thus, DM temperature grows due to the drag term in Eq. (15) and it can be seen in Fig. (1), temperature of DM for $B_{0}=10^{-5}$ G is larger compared to $B_{0}=10^{-6}$ G. Figure (1b) shows the temperature evolution of gas and DM for different DM-baryon interaction cross-section when magnetic field $B_{0}=10^{-6}$ G and DM mass $m_{d}=0.1$ GeV are fixed. Larger the $\hat{\sigma}$ between gas and DM, more heat transfers from gas to the DM and cools the gas efficiently. For $B_{0}=10^{-6}$ G and $\hat{\sigma}=10^{-41}~{}{\rm cm}^{2}$, temperature evolution for different dark matter mass is shown in Fig. (1c). As we increase the DM mass from $0.1$ GeV to $1$ GeV, temperature of both the DM and gas increases and becomes more efficient for large dark matter mass Muñoz et al. (2015). This drag heating is important when mass of DM is comparable or greater than 1 GeV Muñoz et al. (2015). Below $z\sim 50$, in addition to heating due to the drag term, magnetic heating also contribute to the gas temperature, hence the gas temperature for $m_{d}=1$ GeV is higher than $m_{d}=0.1$ GeV.

Correlation between mass of DM and baryon-DM cross section from EDGES observation

In this subsection, we analyze the effect of $B_{0}$, $m_{d}$ and $\hat{\sigma}$ on gas and dark-matter temperature. In Fig. (2), we study constraints on $m_{d}$ and $\hat{\sigma}$ for $T_{21}\simeq-500$ mK ($T_{\rm gas}\simeq 3.26$ K) and $-300$ mK ($T_{\rm gas}\simeq 5.2$ K). Here, to calculate $T_{21}$ we have taken $x_{\alpha}\gg 1$. Thus, from equation (2), $T_{S}\approx T_{\rm gas}$ and one can calculate $T_{21}$ from equation (3). Here note that, we do not include heating due to the Ly$\alpha$ background and the effect of this additional heating become significant for redshift $z\lesssim$17 for the fiducial model considered in Refs. Kovetz et al. (2018); Harker et al. (2015); Mirocha et al. (2015). We discuss this point in the next subsection. In Fig. (2), we consider cases $B_{0}=3.48\times 10^{-6}$ G, $2.67\times 10^{-6}$ G, $10^{-6}$ G and $10^{-9}$ G and solve equations (11-17) for $T_{\rm gas}\simeq 3.26$ and $5.2$ K at $z=17$ to get $m_{d}$ vs $\hat{\sigma}$ plots. The solid and dashed lines represent the cases when $T_{21}\simeq-500$ mK and $-300$ mK respectively. The gold and green regions respectively show the CMB-S4 (forecast) and Planck 2015 upper constraint on $\hat{\sigma}-m_{d}$ with 95% C.L. Kovetz et al. (2018); Boddy et al. (2018). The magenta, black, blue and red lines corresponds to $B_{0}=3.48\times 10^{-6}$ G, $2.67\times 10^{-6}$ G, $10^{-6}$ G and $10^{-9}$ G. As we increase the magnetic field strength from $10^{-9}$ G to $\sim 10^{-6}$ G, larger value of $\hat{\sigma}$ is required for $m_{d}\in\{10^{-6},~{}1\}$ GeV to maintain $T_{21}\simeq-500$ or $-300$ mK at z=17. To get EDGES upper limit on $T_{21}$ (i.e. $-300$ mK), required $\hat{\sigma}$ is smaller compared to the case when $T_{21}=-500$ mK. This is because we need to transfer less energy from gas to the DM to obtain EDGES upper limit on $T_{21}$. We get the upper limit on magnetic field strength $2.67\times 10^{-6}$ G by CMB-S4 (forecast) constraint on $\hat{\sigma}-m_{d}$ and maintaining $T_{21}\simeq-300$ mK at z=17. For $B_{0}=2.67\times 10^{-6}$ G, $m_{d}\gtrsim 10^{-2}$ GeV is excluded. By Planck 2015 constraint on $\hat{\sigma}-m_{d}$, the allowed maximum strength of magnetic field is $3.48\times 10^{-6}$ G by requiring EDGES upper constraint on $T_{21}$ at z=17. For $B_{0}=3.48\times 10^{-6}$ G, mass of dark-matter $\gtrsim 1\times 10^{-2}$ GeV is excluded. Similarly, for the $B_{0}=10^{-6}$ G, $m_{d}\gtrsim 8\times 10^{-1}$ GeV is excluded by CMB-S4 forecast. When the DM mass approaches mass of hydrogen, the drag term in equation (15) also starts to contribute in heating of the gas in addition to the magnetic heating. Therefore, higher mass of dark matter is excluded for higher magnetic field. As discussed in Muñoz et al. (2015), when $m_{d}\gg 1$ GeV, the drag term heat up both the gas and DM in such a way that we can not obtain $T_{\rm gas}=3.26$ K at $z=17$ as required for the EDGES signal. There is a independent bound on the primordial magnetic field from CMB of the order of $\lesssim$ nG Trivedi et al. (2012, 2014). This constraint, in our analysis, restricts value of $\hat{\sigma}$. Here, we note that in our analysis further decreasing value of $B_{0}$ below $10^{-9}$ G, does not change our result in significant way.

Effect of strong magnetic field on brightness temperature

We have discussed above that, as we increase the strength of the magnetic field, for a fix DM mass and interaction cross section, temperature of the gas increases. In figure (3), we plot 21 cm differential brightness temperature with redshift for different magnetic field strengths. This figure is obtained by keeping $m_{d}=0.1$ GeV and $\hat{\sigma}=6.22\times 10^{-15}$ GeV*-2* constant. In figure (3a), to plot $T_{21}$ we assume infinite Ly$\alpha$ coupling ($x_{\alpha}\rightarrow\infty\Rightarrow T_{S}\simeq T_{\rm gas}$) and do not include the X-ray heating. For $B_{0}=1\times 10^{-6}$ G, the 21 cm line absorption signal reported by EDGES (i.e. $-500$ mK) can be explained. In figure (3b), we include the X-ray heating and consider finite Ly$\alpha$ coupling ($x_{\alpha}$) Kovetz et al. (2018); Harker et al. (2015); Mirocha et al. (2015); Zygelman (2005). As we decrease $B_{0}$ from $2\times 10^{-6}$ G, the minimum value of $T_{21}$ profile decreases. For the case when $B_{0}=1\times 10^{-6}$ G (blue solid line), minimum of $T_{21}$ profile is well below the EDGES upper limit on $T_{21}$ (i.e. $-300$ mK—magenta dashed line). In figure (3a), when there is no effect of X-ray heating, $T_{21}=-300$ mK corresponds to $B_{0}=1.35\times 10^{-6}$ G. Thus, we need to lower $B_{0}$ values when the X-ray heating is included to get desired value of $T_{21}$. As shown in Fig.(3), brightness temperature is suppressed by the increase of the strength of the magnetic field and it can even erase the standard 21 cm signal when the magnetic field strength increases above $\sim 1\times 10^{-6}$ G. This sets the upper limit on the strength of the magnetic field for $m_{d}=0.1$ GeV and $\hat{\sigma}=6.22\times 10^{-15}~{}{\rm GeV}^{-2}$.

VI Conclusion

Magnetic field in Minoda et al. (2019) have shown to heat the cosmic gas during the cosmic dawn era by the ambipolar diffusion and the turbulence decay. Since, it could erase the observed 21 cm absorption signal, one can calculate the upper bound on the magnetic field. One of the promising mechanism to explain the absorption signal of the 21 cm line is to have interaction between the dark-matter and baryons Barkana (2018); Barkana et al. (2018). In this work, we have shown that in the presence of such an interaction the upper bound on the strength of magnetic fields can significantly be altered. The magnetic-energy converted to the thermal energy, heat both the gas and dark matter. This is an extra heating effect in addition to the drag heating. The drag term heats the DM and baryons, but in the lower range of dark-matter mass ($\ll 1$ GeV) it favors cooling of gas compared to heating due to the relative motion between DM and gas. To explain the observed anomaly in the 21 cm signal by the EDGES, a large baryon-DM scattering cross-section is required to balance the magnetic heating effect. An earlier saturation occurs in baryon-DM cross-section with respect to the DM mass in the presence of the strong magnetic fields. Considering upper bound on $\hat{\sigma}-m_{d}$ by Planck 2015 Boddy et al. (2018) and EDGES upper constraint on $T_{21}$ (-300 mK) at $z=17$ Bowman et al. (2018), we found upper bound on the magnetic field strength: $B_{0}=3.48\times 10^{-6}$ G, while considering CMB-S4 forecast constraint Kovetz et al. (2018) we get $B_{0}=2.67\times 10^{-6}$ G for the dark matter mass $\lesssim 10^{-2}$ GeV.

Acknowledgments: We would like to thank the anonymous reader whose comments has helped us in improving presentation of our results.