We investigate the relic abundance of asymmetric dark matter where the asymmetric dark matter is non–thermally produced from the decay of heavier particles in addition to the usual thermal production. We discuss the relic density of asymmetric dark matter including the decay of heavy particles in low-temperature scenarios. Here, we still assume that the Universe is radiation-dominated and there is asymmetry before the decay of heavy particles. We obtain an increased abundance of asymmetric dark matter when there is an additional contribution from the decay of heavier particles. Finally, we find the constraints on the asymmetry factor and annihilation cross-section using Planck data.
Ma Qingshan, Hoernisa Iminniyaz. Asymmetric dark matter abundance including non-thermal production[J]. Communications in Theoretical Physics, 2023, 75(8): 085401. DOI: 10.1088/1572-9494/acde68
1. Introduction
The dark matter problem has undoubtedly been one of the most challenging puzzles for cosmologists and particle physicists in recent years. Although we have the most striking evidence for dark matter, what constitutes dark matter is still a mystery to us. Asymmetric dark matter is one alternative, which arises from the fact that its abundance is just five times larger than the baryon asymmetry [1]. There may be an indication of a common origin, which is responsible for the baryonic asymmetry and dark matter. In an asymmetric dark matter scheme, dark matter particles have distinct antiparticles.
In a standard cosmological scenario, it is supposed that asymmetric dark matter is produced thermally. Based on the assumption that weakly interacting massive particles constitute dark matter, the relic density of asymmetric dark matter is determined by the freeze-out condition [2, 3]. In this scenario, it is assumed that the reheating temperature TR is much larger than the freeze-out temperature TF and the asymmetric dark matter particles and antiparticles were in thermal equilibrium in the early Universe. The final relic density of asymmetric dark matter mainly depends on the particle–antiparticle asymmetry, and the indirect detection signals from annihilation are suppressed because there are very few antiparticles here in the end. On the other hand, because we have observational evidence from Big Bang Nucleosynthesis for temperatures T ≲ (few) MeV [4–6], it is also interesting to explore the asymmetric dark matter in scenarios with TR ≲ TF [7–9]. Asymmetric dark matter can be produced non-thermally in low reheating temperature scenarios since heavy particles decay (e.g. moduli decay) into asymmetric dark matter. In [7], the authors discuss a cogenesis mechanism in which the baryon asymmetry of the Universe and dark matter abundance are simultaneously produced at low reheating temperature. Reference [8] explored the non-thermal production of dark matter including asymmetric dark matter in the scheme of the Minimal Supersymmetric Standard Model and string-motivated models. Relic abundance of dark matter was calculated in [10] including the decay of unstable heavy particles to dark matter particles for low reheating temperature.
In our work, we investigate the relic density of asymmetric dark matter in a low reheating temperature scenario when the heavy unstable particles φ decay into asymmetric dark matter particles and antiparticles. We suppose that the asymmetry already existed before the decay of heavy particles into asymmetric dark matter and that they decay into particles and antiparticles in the same amount. We assume that the metastable particles never dominate the energy density of the Universe. There are several issues that are not present in the radiation-dominated epoch in low reheating temperature scenarios since the entropy density is not conserved and non-perturbative inflaton decays may be significant [11–13]. We avoid these problems by considering the abundance of asymmetric dark matter at some initial temperature T0 as a free parameter. We find that the decay of heavier particles to asymmetric dark matter changes the total relic density of asymmetric dark matter. The relic abundances of asymmetric dark matter particles and antiparticles are both increased when there is an additional contribution from the decay of heavier particles.
The outline of the paper is as follows. In section 2, we discuss the relic density calculation of asymmetric dark matter including the decay of heavy particles into asymmetric dark matter in low-temperature scenarios. Here, we assume that the Universe is still radiation-dominated and additional entropy production from the decay of heavy particles is negligible. We find the approximate analytical solution for the relic density of asymmetric dark matter when there is thermal and non-thermal production. In section 3, the constraints on parameter spaces are obtained by using Planck data. The last section is devoted to a brief summary and conclusions.
2. Relic abundance of asymmetric dark matter including non-thermal production
We discuss the scenario where the unstable heavy particles φ decay into asymmetric dark matter particles and antiparticles. Here, we assume that the heavy-particle φ decays out of thermal equilibrium. Therefore, φ production is negligible. However, we include the thermal and non-thermal production of asymmetric dark matter particles.
In the standard computation of relic density for the asymmetric dark matter, it was assumed that the reheating temperature of the Universe TR is much higher than the freeze-out temperature. In this case, the reheating era has no effect on the final relic density of asymmetric dark matter. On the other hand, the constraints on the reheating temperature originate from Big Bang Nucleosynthesis and TR ≳ 1 MeV [4–6]. Therefore, we consider the case where TR < TF [7–9]. The asymmetric dark matter particles and antiparticles never reach thermal equilibrium due to the low reheating temperature. Thus, we treat the asymmetric dark matter abundance at some initial temperature T0 as a free parameter.
The evolution of relic abundance including the decay of heavy particles is more complicated than the usual thermal production case. Here, we still assume that the Universe is radiation dominated and the entropy production is negligible. We suppose that asymmetry existed before the decay of heavy particles and they decay into particles and antiparticles in the same amount. The coupled Boltzmann equations for the number densities nχ and ${n}_{\bar{\chi }}$ of asymmetric dark matter particles χ and antiparticles $\bar{\chi }$ are as follows:
where H is the expansion rate of the Universe, and 〈σv〉 is the thermal average of the cross-section multiplied by the relative velocity of annihilating asymmetric dark matter particles and antiparticles. Here, we assume that only $\chi \bar{\chi }$ pairs can be annihilated into Standard Model particles. nχ,eq is the equilibrium number density. N is the average number of asymmetric dark matter particles and antiparticles produced in φ decay, while Γφ is the decay rate. Here, it is assumed that φ does not dominate the total energy density. Therefore, the comoving entropy density stays almost constant throughout. The analytical solution for equation (3) is easily obtained. Then, equations (1) and (2) can be solved approximately in low-temperature scenarios.
In order to simplify equations (1), (2) and (3), we introduce the following dimensionless variables ${Y}_{\chi ,\bar{\chi }}={n}_{\chi ,\bar{\chi }}/s$, Yφ = nφ/s and x = m/T with s being the entropy density, then,
Assuming that entropy per comoving volume is conserved, we obtain,
$\begin{eqnarray}\dot{s}=-3{sH},\end{eqnarray}$
where we used $H=\dot{R}/R$, where R is the scale factor of the Universe. Inserting s = (2π2/45)g*T3 into equation (7), where g* is the effective number of relativistic degrees of freedom, we obtain,
In the radiation-dominated era, $H=\pi {m}_{\chi }^{2}/({M}_{\mathrm{Pl}}\,{x}^{2})\sqrt{{g}_{* }/90}$, where mχ is the mass of asymmetric dark matter particles and MPl = 2.4 × 1018 GeV is the reduced Planck mass. We then further simplify equations (15) and (16), subtracting equation (16) from (15), then,
where ${\lambda }_{1}=\tfrac{4\pi }{\sqrt{90}}{m}_{\chi }{M}_{\mathrm{Pl}}\sqrt{{g}_{* }}\,$, ${\lambda }_{2}={\left(0.145{g}_{\chi }/{g}_{* }\right)}^{2}$ and ${\lambda }_{3}={M}_{\mathrm{Pl}}{{\rm{\Gamma }}}_{\phi }/(\pi {m}_{\chi }^{2})\sqrt{90/{g}_{* }}$. Here, we used the equilibrium abundance ${Y}_{\chi ,\bar{\chi },\mathrm{eq}}=0.145{g}_{\chi }/{g}_{* }{x}^{3/2}{{\rm{e}}}^{-x(1\pm {\mu }_{\chi }/m)}\,$, where μχ is the chemical potential for particles, and gχ is the number of internal degrees of freedom. The thermal average of σv is usually approximated by a nonrelativistic expansion:
$\begin{eqnarray}\langle \sigma v\rangle =a+6\,{{bx}}^{-1}+{ \mathcal O }({x}^{-2}),\end{eqnarray}$
where a is the s-wave contribution for the limit v → 0 and b is the p-wave contribution for the suppressed s-wave annihilation. Figure 1 depicts the evolution of relic abundance ${Y}_{\chi ,\bar{\chi }}$ as a function of x and it is based on the numerical solutions of equations (19) and (20). Here, the solid red curve is the asymmetric dark matter particle abundance Yχ,nt with the decay of heavy particles, and the black dotted curve is the antiparticle abundance ${Y}_{\bar{\chi },\mathrm{nt}}$. The red dot-dashed and black dashed curves are for asymmetric dark matter particle and antiparticle relic abundances Yχ and ${Y}_{\bar{\chi }}$ without including the decay of heavy particles. Figure 1 shows that the abundances of asymmetric dark matter particles and antiparticles are increased due to the contribution from the decay of metastable heavier particles to asymmetric dark matter. The increase in antiparticle abundance is more remarkable. We start from the initial value of ${Y}_{\bar{\chi }}({x}_{0})$ at the initial temperature x0. We assume that asymmetric dark matter particles and antiparticles never reach thermal equilibrium due to the low reheating temperature. The particle and antiparticle abundances are increased quickly, shortly after starting, and then decreased both for thermal and non-thermal production. When x increases, the abundance finally goes to a constant value. The increase in particle and antiparticle abundances is slightly larger in panels (b) and (d) of figure 1 for a smaller decay rate with respect to panels (a) and (c). Comparing frames (a) with (c), and (b) with (d), it can be seen that for the larger initial value of Yφ(x0), there is a larger increase in the particle and antiparticle abundances. In frames (c) and (d), the antiparticle abundance almost reaches the same level as that of the particle abundance. This opens the possibility to detect the asymmetric dark matter indirectly.
Figure 1. Evolution of Yχ and ${Y}_{\bar{\chi }}$ as a function of x for ${Y}_{\bar{\chi }}({x}_{0}=22)=1\,\times \,{10}^{-12}$, N = 1 and b = 0 a = 4 × 10−26 cm3 s−1, η = 3 × 10−12, mχ = 100 GeV, x0 = 22, gχ = 2, g* = 90. Here, nt means non-thermal production.
The changes in relic density Ωh2 as a function of the cross-section are shown in figure 2. Here, s-wave annihilation cross-section with b = 0 when asymmetry factor η = 3 × 10−12. The solid red curve and dotted black curve represent the relic densities of asymmetric dark matter particles and antiparticles with the decay of heavy particles. The double-dotted red curve and black dashed curve represent the particles and antiparticles only with thermal production. Figure 2 shows when the cross-section increases. Here, the depletion of antiparticle abundance is slower than the case without including heavy-particle decay.
Figure 2. Relic density Ωh2 for particle χ and antiparticle $\bar{\chi }$ as a function of the s-wave cross-section. Here, ${Y}_{\bar{\chi }}({x}_{0}=22)=1\,\times {10}^{-12}$, N = 1 and Yφ = 10−11, Γφ = 10−18, η = 3 × 10−12, mχ = 100 GeV, x0 = 22, gχ = 2, g* = 90.
In what follows, we try to find the approximate analytical solutions of equations (19) and (20). Equation (20) can be solved first. Initially, (x < x0), ${Y}_{\bar{\chi },\mathrm{eq}}\gg {Y}_{\bar{\chi }}$, because η is usually around 10−12 [3], and $\eta {Y}_{\bar{\chi }}$ is small as well. Therefore, equation (20) becomes,
where the scaled Hubble constant h in units of 100 km s−1 Mpc−1 is 0.673 ± 0.098. Here, ρ = nm = s0Ym and the critical density is ${\rho }_{c}=3{H}_{0}^{2}{M}_{\mathrm{Pl}}^{2}$, where s0 ≃ 2900 cm−3 is the present entropy density and H0 is the Hubble constant.
We use the measured dark matter relic density to find constraints on the parameter space. For asymmetric dark matter, particle and antiparticle contributions have to be added as ${{\rm{\Omega }}}_{\mathrm{DM}}={{\rm{\Omega }}}_{\chi }+{{\rm{\Omega }}}_{\bar{\chi }}$.
The contour plot of the s-wave annihilation cross-section and asymmetry factor η is shown in figure 3 when ΩDMh2 = 0.1199. This figure is based on the numerical solutions of the Boltzmann equations (19) and (20). The dotted (black) and dashed (black) curves represent the case where the decay of heavier particles to asymmetric dark matter is included. The solid (red) curve represents the standard case. We find that the required cross-section for Yφ(x0) = 10−11 is nearly one order of magnitude larger than the standard one in order to fall into the observational range. The increased relic density of asymmetric dark matter due to the decay of heavier particles is suppressed by the larger cross-section. Furthermore, the figure shows that the required cross-section is slightly larger when the initial value of Yφ is large. We note that the horizontal part of the contours is not affected by the decay of heavier particles. The antiparticle dark matter density of ${Y}_{\bar{\chi }}$ is exponentially suppressed when the asymmetry increases and the total dark matter relic density is determined by the particle density as ΩDM ∼ mYχ ∼ mη. Therefore, the final relic density of dark matter is independent of cross-section a.
Figure 3. Contour plot of s-wave annihilation cross-section a and asymmetry factor η when ΩDMh2 = 0.1199. Other parameters are the same as in figure 2.
4. Summary and conclusion
The relic abundance of asymmetric dark matter, including the decay of heavier particles, is discussed in low reheating temperature scenarios. Here, we assume that there is asymmetry before the unstable heavy particles decay into asymmetric dark matter particles and antiparticles in the same amount. We found that relic abundances of asymmetric dark matter particles and antiparticles are both increased when there is an additional contribution from the decay of heavier particles. The increase in particle and antiparticle abundances due to the non-thermal production depends on the decay rate of heavy particles and slightly on the initial value of heavy-particle abundance. Comparison with the usual thermal production shows that the depletion of the antiparticle abundance is slower for the case of including decay of heavy particles.
Finally, we used the observed dark matter abundance to obtain the constraints on the annihilation cross-section and asymmetry factor when there is a contribution from the decay of heavier particles to asymmetric dark matter. We found that the required annihilation cross-section with the decay of heavy particles is almost one order of magnitude larger than the case without including non-thermal production. These results are important for our understanding of asymmetric dark matter. It is possible to detect asymmetric dark matter through indirect detection due to the increased amount of antiparticle relic density when there is an additional contribution from the decay of heavy particles.
This work is supported by the National Natural Science Foundation of China (Grant Nos. 2020640017, 2022D01C52 and 11765021).
IchikawaKKawasakiMTakahashiF2005 The Oscillation effects on thermalization of the neutrinos in the Universe with low reheating temperature Phys. Rev. D 72 043522
BarmanBBorahDDasS JRoshanR2022 Non-thermal origin of asymmetric dark matter from inflaton and primordial black holes J. Cosmol. Astropart. Phys. JCAP03(2022)031