Relic Density of Asymmetric Dark Matter in Modified Cosmological Scenarios
Iminniyaz Hoernisa, Salai Burhan, Lv Guo-Liang
School of Physics Science and Technology, Xinjiang University, Urumqi 830046, China

 

† Corresponding author. E-mail:

Supported by the National Natural Science Foundation of China under Grant Nos. 11365022 and 11765021

Abstract
Abstract

We discuss the relic abundance of asymmetric Dark Matter particles in modified cosmological scenarios where the Hubble rate is changed with respect to the standard cosmological scenario. The modified Hubble rate leaves its imprint on the relic abundance of asymmetric Dark Matter particles if the asymmetric Dark Matter particles freeze–out in this era. For generality we parameterize the modification of the Hubble rate and then calculate the relic abundance of asymmetric Dark Matter particles and anti–particles. We find the abundances for the Dark Matter particles and anti–particles are enhanced in the modified cosmological models. The indirect detection signal is possible for the asymmetric Dark Matter particles due to the increased annihilation rate in the modified cosmological models. Applying Planck data, we find the constraints on the parameters of the modified cosmological models.

1 Introduction

Recent Planck team’s full sky measurements of Cosmic Microwave Background (CMB) temperature anisotropies provide the precise values for the cosmological parameters, the amount of cold Dark Matter (DM) relic density is given as

where h = 0.673 ± 0.098 is the present Hubble expansion rate in units of 100 km · s−1 · Mpc−1.[1]

Although the precise value of the Dark Matter relic density is given by the observations, the nature of Dark Matter is still a mysterious question for scientists. Experimental searches are ongoing to find the Dark Matter using direct detection or indirect detection. Parallel to the experiments, there are also many theoretical works which attempts to disclose the character of Dark Matter. One of the general belief is that the Weakly Interacting Massive Particles (WIMPs) are the good candidates for Dark Matter. Majorana particle neutralino is one of the favorite WIMP which is appeared in supersymmetry and its particle and anti-particle are the same. Because we have no evidence that shows that the Dark Matter should be Majorana particle, we still have other possibilities that the Dark Matter can be Dirac particles for which its particle and anti-particle are different.[23] The relic density of asymmetric Dark Matter is calculated in the standard cosmological scenario.[45] In the standard model of cosmology, it is assumed that Dark Matter particles were in thermal and chemical equilibirum in the radiation-dominated epoch after the period of last entropy production. Dark Matter particles decouple when they become nonrelativistic while the universe cools down. For asymmetric Dark Matter, the asymmetry is created well before the decoupling of the asymmetric Dark Matter particles according to the assumption, therefore there are more particles than the anti-particles in the beginning. Dark Matter anti-particles are annihilated away with the particles and the rest are the Dark Matter particles in the present universe. The direct ditection can be used to find the asymmetric Dark Matter particles in principle.

On the other hand, there is no observational evidence before Big Bang Nucleosynthesis (BBN). The alternative cosmological models like scalar-tensor gravity,[6] quintessence models with a kination phase,[710] brane world cosmology[1114] or the late decaying of inflaton model[1517] and anistropic expansion[1819] predict that the expansion rate may have been different from the one expected by the standard Friedmann-Robertson-Walker model before BBN. If this occurs around the time when the asymmetric Dark Matter particles decouple from the thermal equilibrium, the relic abundance of asymmetric Dark Matter particles are affected by the modified expansion rate of the universe. The papers[2025] already discussed the effects of the modified expansion rate on the relic density of asymmetric Dark Matter particles. In those papers, the authors concluded that the enhanced expansion rate leads to the earlier asymmetric Dark Matter particles freeze-out and enhanced relic density. In the modified cosmological models, the difference of the abundance between the particle and anti-particle is not large for reasonable annihilation cross sections. In order to match with the observed range of the Dark Matter relic density, the annihilation cross section should be large. Then the annihilation rate is increased and it leads to the possible indirect detection signal for asymmetric Dark Matter as well.[21,23,25] In this work, we present the generalized parametrization of the modification of the expansion rate of the early universe and derive the general equation to compute the relic density of asymmetric Dark Matter particles. Then we use the Planck data to find the constraints on the parameter space of the modified cosmological models.

The next section is devoted to the analysis of the relic density of asymmetric Dark Matter particles in modified cosmological scenarios both in numerical and analytical ways. In Sec. 3, using the Planck data, we find the constraints on the parameter space for the modified cosmological models. There is conclusion and summary in the final section.

2 Relic Density of Asymmetric Dark Matter in Modified Cosmological Scenarios

In this section, we introduce the parameterization of modification of the expansion rate of the universe and then derive the generalized equation to calculate the relic density for asymmetric Dark Matter. To describe the effects of the revised Hubble expansion rate in the early universe, we add the new dark density to the Friedmann equation as

where

is the radiation energy density with g* being the effective number of relativistic degrees of freedom and

here ρD parameterizes the modification of the expansion rate. nD is the constant which parametrizes the behavior of the energy density. We introduce

where T0 is the reference temperature, which approaches to the freeze-out temperature, then

The approximation holds in a range of temperatures where g* does not change sizably with respect to its value at T0. At BBN time, the radiation energy density remains dominant, and the new term in the Hubble rate should vanish in order not to conflict with the BBN prediction. When nD = 6, the new dark density corresponds to the quintessence model with kination phase; and nD = 8 refers to the brane world cosmology or the late decaying of inflaton field.

The number densities of the asymmetric Dark Matter particles are resolved by the following Boltzmann equations,

Here only pairs can annihilate into Standard Model (SM) particles according to the assumption. nχ, are the number densities of particle and anti-particle. ⟨ σ v ⟩ represents the product of thermal average of the annihilation cross section and the relative velocity of the Dark Matter particle and anti-particle. The equilibrium number densities of χ and are (in equilibrium, the chemical potential of the particles )

where gχ is the number of intrinsic degrees of freedom of the particles.

New quantities Yχ = nχ/s, , and x = m/T are introduced in order to solve Eqs. (7), where

is the entropy density with the effective number of entropic degrees of freedom g*s. Assuming that the universe expands adiabatically and using the entropy conservation s R3 where R is the scale factor of the universe, we obtain

Inserting Eq. (9) into Eq. (10) and using x = m/T, we have

assuming g*g*s and dg*s/dx ≃ 0, we have

Then using Eq. (12), with the modified expansion rate, the Boltzmann equations (7) are rewritten as

where , MPl is the reduced Planck mass.

From Eqs. (13), (14), we obtain

This hints that the difference of the abundance of the particles and anti–particles are conserved,

where Λ is constant. Then Eqs. (13), (14) become

where . The annihilation cross section of WIMP is expanded in the relative velocity v between the annihilating WIMPs in most cases. The thermal average of σv is

where a, b represent the s-wave contribution for the limit v → 0 and the p-wave contribution for the suppressed s-wave annihilation.

2.1 Numerical Solutions

In the standard frame of the Dark Matter particle evolution, it is assumed the particles χ and were in thermal equilibrium at high temperature, when the temperature drops below the mass (T ˂ m) for m ˃ |μχ|, the equilibrium number densities for particles and anti-particles are decreased exponentially.[5,26] The decline of the interaction rates Γ for particle and for anti-particle below the Hubble expansion rate H leads to the decoupling of the particles and anti-particles from the equilibrium state where the temperature at the decoupling is called the freeze-out temperature. After the decoupling, the number densities of the particles and anti-particles in a co-moving space are nearly constant.

Figure 1 demonstrates the evolution of the abundances for particle and anti–particle for different nD and η. These plots are based on the numerical solutions of Eqs. (17), (18). As we mentioned earlier nD = 6 corresponds to the quintessence scenario with kination phase while nD = 8 is for the brane world cosmology or the late inflaton decay model. We find the abundances for particle and anti-particle are increased in those models. The increase depends on the enhancement factor η, when nD = 6, we see that the increase is larger for η = 104 in Fig. 1(b) compared to η = 103 in Fig. 1(a). The situation is the same for nD = 8 in Figs. 1(c) and 1(d). On the other hand, it is shown the increase of the relic abundance is larger for smaller nD when η takes the same value. The reason for the increases of the relic abundances for particle and anti–particle are the enhanced Hubble expansion rate in those modified cosmological models. The enhanced Hubble expansion rate of the universe leads to the early decay of the asymmetric Dark Matter particles from the equilibrium state. Therefore, there are more particles and anti-particles in the modified scenarios with respect to the standard cosmological scenario.

Fig. 1 (Color online) Evolution of Yχ and as a function of the inverse-scaled temperature x for nD = 6 and nD = 8. η = 103 for panels (a), (c) and η = 104 for panels (b), (d). Here Λ = 5 × 10−12, m = 100 GeV, gχ = 2, g* = 90, a = 1.31 × 10−25 cm3 · s−1, b = 0, x0 = 25.
2.2 Analytical Solutions

First, we find the analytic solution for Eq. (18), and then obtain the solution for Eq. (17) using the relation . In terms of , Eq. (18) is written as

Repeating the same method which is mentioned in Ref. [5], we obtain

for high temperature. The inverse-scaled freeze-out temperatures xF and are determined by this solution.[5] For sufficiently low temperature, we obtain

where

The relic abundance for χ particle is

For , the solutions for Yχ and match with the constraint . The final abundance is expressed as

where

The predicted present total Dark Matter relic density is then

According to the assumption that the deviation is of the same order of the equilibrium value of , the freeze-out temperature for is fixed using the equality where ξ being a constant, usually .[26]

3 Constraints

The Dark Matter relic density is given by Planck date as in Eq. (1). The contour plot of the s-wave annihilation cross section a and asymmetry factor Λ is displayed in Fig. 2 when the Dark Matter relic density is ΩDM h2 = 0.1199. Here the dashed (black) line is for nD = 6, the double dot–dashed (blue) line is for nD = 8 and the dotted (black) line is for the standard cosmology. We find the cross sections are larger for the modified cosmological models with respect to the standard cosmology. The relic density is larger in modified cosmology due to the enhanced Hubble rate, therefore the cross section should be increased in order to satisfy the observed range of the Dark Matter relic density. Of course the increased value of the cross section depends on the parameter nD. For smaller nD (dashed line), the larger cross section is needed to satisfy the observed range of Dark Matter relic density, for example, the cross section changes from 4.6 × 10−25 cm3 · s−1 to 3.5 × 10−23 cm3 · s−1 for nD = 6 and 1.8 × 10−25 cm3 · s−1 to 3.5 × 10−23 cm3 · s−1 for nD = 8. Here Λ takes the value from Λ = 9.0 × 10−14 to Λ = 4.4 × 10−12. We can analyze this from panels (a) and (c) in Fig. 1. For the same η, the increase of relic abundance is larger for nD = 6 with respect to the case nD = 8, therefore the annihilation cross section should be large to satisfy the observed range of the Dark Matter abundance for smaller nD.

Fig. 2 (Color online) The allowed region in the (a,Λ) plane for different nD and standard cosmology when the Dark Matter relic density Ω h2 = 0.1199. Here m = 100 GeV, gχ = 2 and g* = 90, x0 = 25, η = 103. The three dot-dashed (red) line is the limiting cross section, which corresponds to the Fermi-LAT bounds for channel and the thick (red) line is for channel.[27]

Since the annihilation cross sections of the asymmetric Dark Matter particles in the modified cosmological scenarios are increased, then the annihilation rate is enhanced and the indirect detection is possible for the asymmetric Dark Matter particles, which is originally assumed to be detected only through the direct detection due to the suppressed abundance of the anti-particles. Therefore we can use the Fermi Large Area Telescope (Fermi-LAT)[27] data to find the limiting cross sections for the asymmetric Dark Matter in nonstandard cosmological scenarios. In the case of the symmetric Dark Matter, the annihilation rate is[21]

where ρDM = ρcrit ΩDM and ⟨ σselfv ⟩ is the thermal average of the self annihilation cross section. For asymmetric Dark Matter, the annihilation rate is

here we use and nχ = ρχ/m, then Eq. (29) becomes

According to the analysis of Ref. [21], the ratio of Γasym and ΓFermi must be less than one,

where ⟨ σ vFermi is the derived cross section bound from the Fermi-LAT data.

Then using Eq. (27), we obtain

Fermi-LAT provides the following uppor bounds on the annihilation cross section a = 1.34 × 10−25 cm3 · s−1 for channel and a = 1.38 × 10−24 cm3 · s−1 for channel while m = 100 GeV.[27] Applying Planck data Eq. (1) and Fermi-LAT data to Eq. (32) together, we obtain the limiting bound for the annihilation cross sections in the modified cosmological scenarios. We find when Λ = 9 × 10−14, for channel, the corresponding constrained cross sections are ⟨ σ v ⟩ ≤ 2.7 × 10−25 cm3 · s−1 (three dot–dashed line in Fig. 2); for channel, it is ⟨ σ v ⟩ ≤ 2.8 × 10−24 cm3 · s−1 (thick red line in Fig. 2).

Figure 3 demonstrates the relation between the s–wave annihilation cross section a and the enhancement factor η when Ω h2 = 0.1199. In this plot, we take the value for asymmetry factor as Λ = 3.5. × 10−12 which is in the range of the values to obtain the observed Dark Matter abundance. It is indeed shown in Fig. 2. Here the thick (red) line is for nD = 8 and the dashed (blue) line is for nD = 6. While nD = 8, for the cross sections from 6.2 × 10−26 cm3 · s−1 to 1.9 × 10−24 cm3 · s−1, η changes from 0.1 to 3.1 × 105 and to 9.0 × 103 when nD = 6. For the same cross section, the enhancement factor η is increased for larger nD. We can understand the reason from Fig. 1. In the frame (b) and (c) of Fig. 1, the cross section is same for nD = 6 and nD = 8. We find for smaller nD, there are more increase of asymmetric Dark Matter relic density for η = 104 in Fig. 1(b) than the case for η = 103 in Fig. 1(c). In order to fall the observed range of Dark Matter relic density, η should be decreased (increased) for small (large) value of nD.

Fig. 3 (Color online) The allowed region in the (a,η) plane for b = 0 when the Dark Matter relic density Ω h2 = 0.1199. Here we take m = 100 GeV, Λ = 3.5 × 10−12, gχ = 2 and g* = 90, x0 = 25. The dashed (blue) line is for nD = 6 and the thick (red) line is for nD = 8.

In this work, we use the observed Dark Matter abundance to derive the constraints on the enhancement factor for the fixed asymmetry factor Λ. In the quintessence model (nD = 6), the constraint on the enhancement factor from BBN is η ˂ 30,[29] which matches partly with the limit η ˂ 9.1 × 103 that we obtained. In paper,[29] the bound on η was obtained at the temperature T = 10 MeV, here we take this bound as an example of a reference. As we mentioned earlier nD = 8 corresponds to the brane Randall-Sundrum II model,[11] which in our notations represents η = ρrad(T0)/(2λ) with λ being the tension of the brane which is related to the 5-dimensional Planck mass M5 by , here MP is the 4-dimensional Planck mass.[23] From Fig. 3, we see that η ˂ 3.1 × 105, this implies λ ˃ 1.9 GeV4 and it leads M5 ˃ 103 TeV. The BBN bound is λ ˃ 1 MeV4 and so M5 ˃ 30 TeV.[30]

4 Summary and Conclusions

We have studied the relic density of asymmetric Dark Matter particles in the modified cosmological scenarios. We have no observational evidence for the early pre-BBN universe. There are different cosmological models, which predicted the Hubble expansion rate before BBN may be faster or slower than the standard expansion law. If the Hubble expansion rate is changed in the modified cosmological scenarios, it affects the relic density of asymmetric Dark Matter particles. We have discussed in which extent the modification of the Hubble expansion rate affects the relic density of asymmetric Dark Matter in the modified cosmological models.

The Hubble expansion rate of the universe is increased in the modified cosmological scenarios like quintessence with kination phase, brane world cosmology or the late inflaton decay model and etc. For generality we found the parameterization of the enhancement of the Hubble expansion rate. Then we calculated the relic density of asymmetric Dark Matter particles in the modified cosmological scenarios where the expansion rate is faster than the standard cosmological scenarios. We found the resulting relic abundances are increased. Since the expansion rate of the universe is enhanced during the era of asymmetric Dark Matter decoupling, the asymmetric Dark Matter particles decouple from the thermal equilibrium earlier than the standard cosmological scenario and there are increased relic densities for particles and anti-particles. The increases depend on the enhancement factor η as well as the cross section and the factor nD. While using Planck date for the relic density of asymmetric Dark Matter, we found the cross sections are boosted for the modified cosmological models in order to satisfy the observed Dark Matter relic density. From the contour plot of (a, η), we found η is decreased (increased) for small (large) value of nD when the cross section is same for different nD.

The generalised formula obtained provides the relic density calculation of the asymmetric Dark Matter particles for the different modified cosmological scenarios where the expansion rate is enhanced. In the standard cosmological model, the abundance of minority component of asymmetric Dark Matter is suppressed. For the modified cosmological models, we obtained the enhanced relic densities both for particles and anti-particles, of course the increased amount of the relic abundance depends on the size of the annihilation cross section. The annihilation cross section should be boosted in order to provide the observed range of Dark Matter relic density (1). The corresponding annihilation rate is increased too. This effect is contrary to the general belief that the asymmetric Dark Matter is only detected by direct detection. In modified cosmological scenarios, we have the possibility to detect the asymmetric Dark Matter particle in indirect way. The asymmetric Dark Matter particle properties may be constrained by the indirect detection signals.[27] In addition to the observational Dark Matter relic density,[1] there are collider, direct and indirect detection signals, using those data the parameter space is examined before BBN.

Reference
[1] Ade P. A. R. et al. [Planck Collaboration) Astron. Astro-phys. 594 2016 A13 arXiv:1502.01589 [astro-ph.CO]
[2] Nussinov S. Phys. Lett. B 165 1985 55
[3] Belyaev A. Frandsen M. T. Sannino F. Sarkar S. Phys. Rev. D 83 2011 015007 arXiv:1007.4839
[4] Graesser M. L. Shoemaker I. M. Vecchi L. J. High Energy Phys. 1110 2011 110 arXiv:1103.2771 [hep-ph]
[5] Iminniyaz H. Drees M. Chen X. JCAP 1107 2011 003 arXiv:1104.5548 [hep-ph]
[6] Catena R. Fornengo N. Masiero A. et al. Phys. Rev. D 70 2004 063519 arXiv:astro-ph/0403614
[7] Salati P. Phys. Lett. B 571 2003 121 astro- ph/0207396
[8] Rosati F. hep-ph/0309124
[9] Profumo S. Ullio P. JCAP 0311 2003 006 hep-ph/0309220
[10] Pallis C. JCAP 0510 2005 015 hep-ph/0503080
[11] Randall L. Sundrum R. 1999 Phys. Rev. Lett. 83 4690 hep- th/9906064
[12] Stoica H. Tye S. H. H. Wasserman I. Phys. Lett. B 482 2000 205 hep-th/0004126
[13] Abou El Dahab E. Khalil S. J. High Energy Phys. 0609 2006 042 hep-ph/0607180
[14] Okada N. Seto O. Phys. Rev. D 70 2004 083531 hep-ph/0407092
[15] Arbey A. Deandrea A. Tarhini A. J. High Energy Phys. 1105 2011 078 arXiv:1103.3244 [hep-ph]
[16] Arbey A. Mahmoudi F. Phys. Lett. B 669 2008 46 arXiv:0803.0741 [hep-ph]
[17] D’Eramo F. Fernandez N. Profumo S. 2017 JCAP 1705 012 arXiv:
[18] Kamionkowski M. Turner M. S. Phys. Rev. D 42 1990 3310
[19] Barrow J. D. Nucl. Phys. B 208 1982 501
[20] Iminniyaz H. Chen X. Astropart. Phys. 54 2014 125 arXiv:1308.0353 [hep-ph]
[21] Gelmini G. B. Huh J. H. Rehagen T. JCAP 1308 2013 003 arXiv:1304.3679 [hep-ph]
[22] Wang S. Z. Iminniyaz H. Mamat M. Int. J. Mod. Phys. A 31 2016 1650021 arXiv:1503.06519 [hep-ph]
[23] Meehan M. T. Whittingham I. B. JCAP 1406 2014 018 arXiv:1403.6934 [astro-ph.CO]
[24] Abdusattar H. Iminniyaz H. Commun. Theor. Phys. 66 2016 363 arXiv:1505.03716 [hep-ph]
[25] Iminniyaz H. Phys. Lett. B 765 2017 6 arXiv:1604.04251 [hep-ph]
[26] Scherrer R. J. Turner M. S. Phys. Rev. D 33 1986 1585
[27] Ackermann M. et al. [Fermi-LAT Collaboration) Phys. Rev. D 89 2014 042001 arXiv:1310.0828 [astro-ph.HE]
[28] Schelke M. Catena R. Fornengo N. et al. 2006 Phys. Rev. 74 083505 hep- ph/0605287
[29] Dutta S. Scherrer R. J. 2010 Phys. Rev. 82 083501 arXiv:
[30] Durrer R. 2005 AIP Conf. Proc. 782 202 hep-th/0507006