Abstract
We provide a short review of some aspects of dark matter (DM) production in non-standard cosmology. Considering the simplest Higgs portal model as a definite particle physics setup, we consider the impact on the parameter space corresponding to the correct relic density and the complementary experimental constraints of the presence, during thermal production, of an exotic component dominating the energy density of the universe. In the second part of the work, we will focus on the case that such an exotic component satisfies the equation of state of matter and can produce DM non-thermally.
1 Introduction
The solution to the dark matter (DM) puzzle is one of the biggest challenges of modern particle physics. The determination of the mechanism for DM production is a key ingredient in solving this puzzle. Thermal freeze-out is one of the most popular proposals as it relates the DM relic density to a single-particle physics input, the so-called thermally averaged DM pair annihilation cross section. Furthermore, the value of the latter quantity, favored by cosmological observations of the DM abundance (see [1]), is characteristic of weak interactions, leading to the so-called weakly interacting massive particle (WIMP) miracle. The thermal freeze-out paradigm is, however, in increasing tension with null results from DM searches, especially the ones based on the principle of direct detection (DD) (see, e.g., [2, 3] for some reviews). In light of this, an alternative production mechanism, dubbed freeze-in [4], is gaining increasing attention as it can accommodate the correct relic density for very small values of the coupling between the DM and standard model (SM) states, encompassing the aforementioned experimental tensions. Both conventional freeze-in and freeze-out mechanisms rely on the hypothesis of a standard cosmological history of the Universe implying, in particular, that the DM is produced in a radiation-dominated epoch. There are no reasons, in addition to minimality, to enforce a priori such an assumption as we have no confirmed experimental evidence about the cosmological history prior to the Big Bang Nucleosynthesis (BBN). It is then interesting to consider the impact on DM production of a non-standard cosmological evolution of the Universe. By this, we intend the possibility that at some epoch, comprised between the primordial inflation and the BBN, the energy budget of the Universe was dominated by an exotic component, i.e., different by ordinary (and dark) matter and radiation. Such an exotic component impacts DM production in a two-fold manner: it affects the Hubble expansion parameter and the evolution of the temperature of the Universe with time during DM thermal production; it might be itself a source of (non-thermal) production of DM. In this work, we will provide a brief review of some aspects of thermal and non-thermal production of DM in a non-standard cosmological history (please refer to [5] for more extensive discussions).
The remainder of this paper is organized as follows: in Section 2, we will describe the general set of equations for DM production in the presence of a single exotic component, to the energy budget of the Universe, characterized by an arbitrary equation of state parameter. Some analytical approximations for the solution of such equations, in the case that the new component is not a direct source of DM, will also be provided. In Section 3, a reference particle physics framework, i.e., the Higgs portal with scalar DM, will be introduced. The findings of Section 2 will be applied to it. Finally, in Section 4, the case of non-thermal production from an exotic matter component will be reviewed. Again, some examples of the solutions of Boltzmann’s equations will be applied to the scalar Higgs portal. The final section will be devoted to the conclusions.
2 Boltzmann equations
Following [
6], the most general set of Boltzmann equations describing the scenario under concern can be written as follows
1:
where
.
ρϕrepresents the energy density of an exotic component
, with an equation of state relating its energy density and pressure,
. This component can dominate the energy density of the Universe during certain stages. As further detailed in the next section,
ω= 0 is the most popularly considered option, corresponding to a so-called “early matter domination period,” due, for example, to heavy metastable particles. On a similar footing, one might consider additional radiation components, i.e.,
ω= 1/3. Other popular examples are kination (
ω= 1)-dominated [
7–
9] and quintessence (
ω= −1)-dominated Universe [
10,
11]. Different values of
ωcan be expected by considering scalar fields with suitable choices for their potentials (refer to [
12,
13] for example). Finally, scenarios of braneworld cosmologies [
14] and modified gravity theories [
15,
16] can be described via the formalism illustrated below. For a more extensive review of the possible origins of non-standard cosmologies, we refer to [
17]. The time evolution
is governed by a decay rate Γ
ϕ, which must ensure that it disappears before the onset of BBN. At the moment of decay, the energy stored in
is passed to the primordial plasma, increasing its entropy, and possibly to the DM. As just pointed out, in light of the decay of
, the entropy density
sof the primordial plasma is not a conserved quantity any longer; consequently, we need to explicitly consider a Boltzmann equation (alternatively one can consider, instead, an equation for the energy density of radiation, refer to [
18–
21] for example). The last equation in (
1) tracks the time evolution of the DM number density
nχ. In addition to the Hubble expansion, it is governed by pair annihilation processes into SM pairs, described by the thermally averaged cross section ⟨
σv⟩, which can be computed as shown in detail in [
22] or via public available packages as micrOMEGAs [
23,
24] or DarkSUSY [
25,
26], and, possibly, by a non-thermal production term depending on Γ
ϕand on a parameter
, measuring the fraction of the energy density of
, which gets converted into DM. The system written in Eq.
1should be combined witht the following equation for the Hubble expansion parameter
.
with
Let us consider the case where
so that
influences DM production only indirectly by altering the expansion rate of the Universe. A recent extensive semi-analytical study of the solution of Eq.
1has been presented in
[6],and we summarize below the results. The non-standard cosmology is parameterized mostly via three quantities:
where
ωis the already mentioned equation of the state parameter;
κindicates, at a reference temperature, the amount of energy density of
, compared to the one of radiation, and can be used to set the initial conditions for Eq.
1.
is finally the temperature at which the standard radiation domination era starts again after the
-dominated epoch. In the so-called instantaneous decay approximation, it can be determined by the following condition:
Such conditions can be used to set Γ
ϕas a function of
without relying on a specific model. To avoid tension with BBN,
≥ 4 MeV is required [
27–
30]. The system Eq.
1can be solved as function of the parameters listed in Eq.
2without referring to a specific model. There are three additional temperature (and hence time) scales.
• The standard freeze-out temperature , i.e., the temperature at which DM annihilation becomes inefficient and DM starts to decouple from the primordial plasma, in the absence of exotic components in the energy budget of the Universe. The freeze-out temperature can be determined by solving the following equation:
with
being the internal degrees of freedom of the DM candidate. For values of ⟨
σv⟩ around the thermally favored value,
∼ 20 ÷ 30.
• The temperature from which, give an initial value for κ, the components becomes the dominant contribution to .
• The temperature at which the presence of the exotic component starts altering the evolution of the plasma temperature with the scale factor.
Having in mind these relevant scales, one can achieve a semi-analytical determination of the DM abundance in some limiting regimes:
• Teq ≪ Tf.o.: DM freeze-out occurs similar to that in the standard radiation domination scenario. The main effect from is represented by the entropy injection during its decay, causing a dilution of thermal abundance of the DM. In the instantaneous decay approximation, the DM relic density can be written as
• ≪ ≪ : In such a regime, freeze-out occurs when the Hubble expansion parameter is dominated by the component.
We are, however, far enough from its decay time so that the relation between the temperature and scale factor is the same as in standard cosmology. In such a case, the DM abundance is again given by the ratio of a thermal abundance and the same dilution factor defined in the previous case. However, the thermal abundance differs from the standard computation as a consequence of a different freeze-out time.
where
xend=
mχ/
Tendand
refers to a modified freeze-out time, with respect to a standard cosmological history, given by
• ≪ ≪ : In such a case, DM freeze-out is affected by the different relation between the temperature and scale parameter. The relic density can be approximated, this time, as
where the freeze-out temperature is obtained by solving the following equation:
• ≪ : the presence of an epoch dominated by the exotic component has no impact on DM production. The relic density is determined as in the standard cosmological model.
3 Results in a specific case of study
As evident from the previous discussion, the framework under consideration can be analyzed in terms of a limited set of parameters without relying on a specific particle physics framework: the initial ratio κ between and radiation energy densities, the equation of state parameter ω, , the DM mass , and the annihilation cross section ⟨σv⟩. However, for a better understanding of the impact of non-standard cosmologies on DM production, it is useful to consider a definite particle model. Our choice falls on the Higgs portal (see [31] for a review) with scalar DM as it allows for maintaining a low number of free parameters. Indeed, the latter model is fully characterized by the following Lagrangian equation:Here, χ is a real scalar DM candidate2 and is the Higgs doublet. After the Higgs obtains a vacuum expectation value (vev), the Lagrangian function generates a trilinear interaction between the Higgs boson and a pair of DM particles, which allows for DM annihilations into SM fermion, gauge boson, and Higgs boson pairs. In addition to the DM mass,the coupling is the only free parameter of the theory. By comparing the DM annihilation rate with the Hubble expansion rate during a radiation-dominated era, one finds that for ≳ 10–5, the DM was capable of being in thermal equilibrium in the Early Universe. One could then apply the standard freeze-out paradigm and find that the correct relic density, Ωχh2 ≈ 0.12, is matched for O (0.1–1) values of the coupling, with an exception of , where the s-channel resonant enhancement of the DM annihilation cross section allows for very small values of the couplings (see Figure 1). For λχ ≤ 10–6, the DM was not capable of thermalizing with the primordial plasma. Nevertheless, the correct relic density can be achieved for λχ ∼ O (10–11) via the freeze-in mechanism3. The parameter space corresponding to thermal freeze-out can be effectively probed experimentally. The most relevant constraints come from DM DD. For our analysis, we have considered the combination of the limits given by XENON1T [33] for ≤ 10 GeV and by LZ [34] for ≥ 10, GeV. For , DD constraints are well-complemented by the bounds from searches of invisible decays of the SM Higgs (refer to [35] for most recent results). In this work, we have adopted the limit ≤ 0.11 and also considered projected increased sensitivities to 0.05 and 0.01 [36, 37].
FIGURE 1
Figure 1 illustrates, via some examples, the impact of a non-standard cosmological history, as illustrated in the previous section, on the parameter space of the scalar Higgs portal, with a focus on the ≤ 100 GeV region, which is mostly subjected to experimental constraints. The different panels of the figure consider some assignations of the (ω, κ) pair and show, in the (, ) bi-dimensional plane, isocontours of the correct relic density for different values of ranging from 5 MeV (approximately the lower bound from BBN) to 1 GeV. To be viable, such isocontours should be (at least partially) outside the brown- and gray-colored regions corresponding to, respectively, the exclusion bounds from DD and invisible Higgs decays. For reference, the isocontour corresponding to the standard freeze-out scenario has also been shown. From the outcome of the plots, one notices that non-standard cosmologies sensitively affect the parameter space corresponding to the correct DM relic density, allowing lower values of the DM couplings and, more interestingly, lower values of the mass. In the case of an additional matter component ω = 0, it remains very difficult to overcome experimental constraints for . To achieve a viable parameter space at low DM masses, one needs to rely on more exotic components with ω = −1/3 and ω = −2/3.
4 Thermal and non-thermal DM in the universe with early matter domination
The most commonly considered scenario with bχ ≠ 0 is the one in which ϕ is an additional matter component, i.e., ω = 0. In this setup, ϕ can be interpreted as a particle field that is always thermally decoupled from the primordial plasma. The existence of these fields is motivated by several particle physics frameworks (refer to [38–44] for some examples). In more recent times, primordial black holes have been proposed as this exotic matter component [45]. In this kind of setups, is customarily referred to as reheating temperature . Although we will adopt the phenomenological determination given by Eq. 3 and use, as well, the parameter κ to fix the initial conditions for the Boltzmann equations, one can determine Γϕ from the model parameters as follows:where depends on the specific underlying model and the initial energy density is assigned as . In this setup, can be extrapolated from the numerical solution of the Boltzmann equations (refer to [21] for a discussion). Equation 1 can be solved in the case of a new matter field Φ via the following change of variables [18, 21]:Such a change in variables allows us to gauge out the terms linear with the Hubble expansion rate so that the system of equations can be rewritten aswhere is defined as
In the abovementioned equations, represents the entropy effective degrees of freedom.
In the following, we will present some examples of numerical solutions to the equations above, again adopting the Higgs portal as a model for DM interactions with the SM. A similar scenario has also been considered in [46]. Before doing this, we briefly illustrate some approximate solutions, following the discussion of [20, 21] (detailed studies of Boltzmann’s equations for non-thermal DM production have also been conducted in [47–49]). Assuming that is sufficiently lower than so that the thermally produced component of the DM gets sufficiently diluted to account for Ωχ to a negligible extent, we distinguish two leading regimes for the solution of Boltzmann’s equations. In the case that interactions between the DM and SM are substantial, non-thermal production of DM can lead to a significant DM annihilation rate, i.e., Γann = ⟨σv⟩, becoming efficient again, compared to the Hubble expansion rate at low temperatures so that there is a compensation between non-thermal production and annihilations. This situation occurs when the number density of non-thermally produced DM exceeds the critical value given byIn the instantaneous decay approximation, the condition can be re-expressed as [43]evidencing that for a fixed , as decreases, one would need stronger DM interactions (encoded in ⟨σv⟩) to match this condition.
In this regime, also dubbed as the re-annihilation regime in the literature [15, 50], the DM relic density can be approximated by an analogous expression as in the standard freeze-out case but replacing with :where is the relic density computed according to the conventional freeze-out paradigm assuming standard cosmology (we remember ). If the interactions between the DM and SM states are not efficient enough, a fraction of the energy initially stored in the Φ field, determined by the parameter , is directly transferred into DM particles. In such a regime, the DM relic density is directly proportional to and reheating temperature .
Figures 2, 3 show some examples of the solution to Boltzmann’s equations for non-thermal production of DM.
FIGURE 2
FIGURE 3
Figure 2 considers the variation in DM abundance as a function of and . In all cases, κ = 10 is considered. As far as the particle physics input is concerned, the values of 100 GeV and 10–2 have been considered for, respectively, DM mass and coupling. These parameter assignations comply with the constraints from DD and invisible Higgs decay. The left panel of Figure 2 considers a fixed value of the reheating temperature, namely, 1 GeV, and different values of bχ. The DM abundance is very weakly dependent on the latter parameter for values above 0.01. This is because the latter is generated in the re-annihilation regime. Indeed, when the abundance of the non-thermally produced DM exceeds (or, equivalently, ), it is described as the quasi-statical equilibrium (QSE) [50] number density:until the latter drops below and gets frozen. By further decreasing the value of , the amount of non-thermally produced DM is not sufficient to reactivate annihilation processes; hence, we have ∝ . Moving to the right panel of Figure 2, we see that < 1 GeV increases as decreases, while such a trend is reversed for < 100 GeV. In agreement with the discussion of [19], this is due to the transition from the re-annihilation regime, corresponding to , occurring at a high reheating temperature, to the regime in which annihilations are not active, corresponding to Ωχ ∝ TR.
In Figure 3, we have considered the impact of a variation in the DM coupling λχ. Again, Figure 3 show the evolution of Yχ with DM mass. The values TR = 1 GeV and = 0.01 have been considered. For values of > 10–5, the DM is in thermal equilibrium in the first stages of the evolution of the Universe. For a value of , starts deviating from the equilibrium value and the two competing effects become relevant: the dilution by the entropy injection due to the late decay of Φ and the non-thermal production process. The final DM abundance is inversely proportional to coupling as the solution tracks the re-annihilation regime. For < 10–5, the DM is initially produced via freeze-in. Such abundance is diluted away by the decay of the Φ field and replaced by a non-thermal population of DM without re-annihilation as the DM interaction rate is too suppressed. This explains the fact that the DM abundance is independent of the value of .
We conclude our analysis with a few remarks about the detection prospects of the scenarios discussed in this paper. On general grounds, distinguishing non-standard cosmological scenarios only via earth-scale experiments, such as the ones based on direct/indirect detection and collider searches, is complicated as they can reconstruct DM particle properties, like the size of the interactions with SM states, while being affected to a negligible extent by the cosmological ones. It is nevertheless evident that, in the presence of a non-standard cosmological history, the parameter space, corresponding to the correct relic density, can vary substantially with respect to the case of thermal freeze-out. Consequently, a hypothetical future signal, for example, at a current or next-generation direct detection facility, possibly incompatible with the expectations of the conventional freeze-out paradigm, would represent a very useful indication (refer to [6] for similar ideas). A conclusive statement would nevertheless require a complementary signal from a probe of pre-BBN cosmology. In this context, gravitational wave (GW) detectors capable of probing the primordial GW background can make a difference [51, 52]. As a final remark, we mention that scenarios of non-thermal production of DM, such as the one discussed in Section 4, in the re-annihilation regime, have already been effectively probed. Indeed, indirect detection experiments (refer to [53, 54]) and CMB probes ([1, 55]) are already sensitive to DM annihilation cross sections of the order of the thermally favored one. Consequently, scenarios of non-thermal production in the reannihilation regime, which reproduce the correct relic density for the annihilation cross section greater than the thermal one, might be strongly constrained or ruled out. Notice anyway that this statement is strictly valid for models with the s-wave-dominated annihilation cross section, i.e., the case in which the value of ⟨σv⟩ at freeze-out and CMB/present times substantially coincides. Finally, structure formation could also provide insights about non-thermal production scenarios as DM properties at that time could deviate from the conventional cold dark matter paradigm, leaving an imprint that could be traced by Lyman-α [56].
5 Conclusion
Thermal freeze-out is a very popular framework, leading to predictive models that are testable via a broad variety of complementary experimental search strategies. It relies, however, on the assumption of standard cosmological history during DM production. There are no a priori reasons to enforce such an assumption. We have provided a brief review of the scenario of thermal and non-thermal production of DM in a non-standard cosmological history, represented by an exotic component, possibly dominating the energy budget of the Universe during DM production and prior to BBN. Although the result might be illustrated in a very general perspective, we have found it convenient to identify a reference model corresponding to the Higgs portal with scalar DM. Assuming only the thermal production of DM, the non-standard cosmological evolution enlarges the parameter spaces complying with constraints from DD and invisible Higgs decays. In the second part of this review, we have focused on non-thermal production, focusing on the case in which the Universe encounters an early matter domination epoch. We have illustrated the relevant Boltzmann’s equations and discussed both numerical and analytical approximations of the solutions.
Statements
Author contributions
GA: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
1.^The system of Boltzmann’s equations can actually be written in the considered form if the following assumptions hold: i) DM self-interactions influence to a negligible extent its number density (they should ensure anyway thermalization of the DM particle with themselves). Such an assumption can be easily satisfied by a suitable assignation of λs. ii) The DM is, during its whole production process, at least in kinetic equilibrium with the primordial plasma.
2.^The phenomenology would be totally analogous in the case of complex scalar DM.
3.^Assuming that an initial DM population was produced during inflation, the correct relic density would also be achieved for λχ = 0 [32].
References
1.
AghanimNAkramiYAshdownMAumontJBaccigalupiCBallardiniMet alPlanck 2018 results. VI. Cosmological parameters. Astron Astrophys (2020) 641:A6. [Erratum: Astron.Astrophys. 652, C4 (2021)]. 10.1051/0004-6361/201833910
2.
ArcadiGDutraMGhoshPLindnerMMambriniYPierreMet alThe waning of the WIMP? A review of models, searches, and constraints. Eur Phys J C (2018) 78:203. 10.1140/epjc/s10052-018-5662-y
3.
ArcadiGCabo-AlmeidaDDutraMGhoshPLindnerMMambriniYet alThe waning of the WIMP: endgame? (2024). arXiv:2403.15860.
4.
HallLJJedamzikKMarch-RussellJWestSM. Freeze-in production of FIMP dark matter. JHEP (2010) 03:080. 10.1007/JHEP03(2010)080
5.
KaneGLKumarPNelsonBDZhengB. Dark matter production mechanisms with a nonthermal cosmological history: a classification. Phys Rev D (2016) 93:063527. 10.1103/PhysRevD.93.063527
6.
AriasPBernalNHerreraAMaldonadoC. Reconstructing non-standard cosmologies with dark matter. JCAP (2019) 10:047. 10.1088/1475-7516/2019/10/047
7.
BarrowJD. MASSIVE PARTICLES AS A PROBE OF THE EARLY UNIVERSE. Nucl Phys B (1982) 208:501–8. 10.1016/0550-3213(82)90233-4
8.
FordLH. Gravitational particle creation and inflation. Phys Rev D (1987) 35:2955–60. 10.1103/PhysRevD.35.2955
9.
VisinelliL(Non-)thermal production of WIMPs during kination. Symmetry (2018) 10:546. 10.3390/sym10110546
10.
SalatiP. Quintessence and the relic density of neutralinos. Phys Lett B (2003) 571:121–31. 10.1016/j.physletb.2003.07.073
11.
ProfumoSUllioP. SUSY dark matter and quintessence. JCAP (2003) 11:006. 10.1088/1475-7516/2003/11/006
12.
ChoiKString or M theory axion as a quintessence. Phys Rev D (2000) 62:043509. 10.1103/PhysRevD.62.043509
13.
Di MarcoAPradisiGCabellaP. Inflationary scale, reheating scale, and pre-BBN cosmology with scalar fields. Phys Rev D (2018) 98:123511. 10.1103/PhysRevD.98.123511
14.
OkadaNSetoO. Relic density of dark matter in brane world cosmology. Phys Rev D (2004) 70:083531. 10.1103/PhysRevD.70.083531
15.
CatenaRFornengoNMasieroAPietroniMRosatiF. Dark matter relic abundance and scalar - tensor dark energy. Phys Rev D (2004) 70:063519. 10.1103/PhysRevD.70.063519
16.
MeehanMTWhittinghamIB. Dark matter relic density in Gauss-Bonnet braneworld cosmology. JCAP (2014) 12:034. 10.1088/1475-7516/2014/12/034
17.
AllahverdiRet alThe first three seconds: a review of possible expansion histories of the early Universe (2020). 10.21105/astro.2006.16182
18.
GiudiceGFKolbEWRiottoA. Largest temperature of the radiation era and its cosmological implications. Phys Rev D (2001) 64:023508. 10.1103/PhysRevD.64.023508
19.
GelminiGGondoloPSoldatenkoAYagunaCE. Effect of a late decaying scalar on the neutralino relic density. Phys Rev D (2006) 74:083514. 10.1103/PhysRevD.74.083514
20.
GelminiGBGondoloP. Neutralino with the right cold dark matter abundance in (almost) any supersymmetric model. Phys Rev D (2006) 74:023510. 10.1103/PhysRevD.74.023510
21.
ArcadiGUllioP. Accurate estimate of the relic density and the kinetic decoupling in non-thermal dark matter models. Phys Rev D (2011) 84:043520. 10.1103/PhysRevD.84.043520
22.
GondoloPGelminiG. Cosmic abundances of stable particles: improved analysis. Nucl Phys B (1991) 360:145–79. 10.1016/0550-3213(91)90438-4
23.
BelangerGBoudjemaFPukhovASemenovA. MicrOMEGAs 2.0: a Program to calculate the relic density of dark matter in a generic model. Comput Phys Commun (2007) 176:367–82. 10.1016/j.cpc.2006.11.008
24.
BelangerGBoudjemaFPukhovASemenovA. Dark matter direct detection rate in a generic model with micrOMEGAs 2.2. Comput Phys Commun (2009) 180:747–67. 10.1016/j.cpc.2008.11.019
25.
GondoloPEdsjoJUllioPBergstromLSchelkeMBaltzEA. DarkSUSY: computing supersymmetric dark matter properties numerically. JCAP (2004) 0407:008. 10.2172/827305
26.
BringmannTEdsjöJGondoloPUllioPBergströmL. DarkSUSY 6: an advanced tool to compute dark matter properties numerically. JCAP (2018) 07:033. 10.1088/1475-7516/2018/07/033
27.
KawasakiMKohriKSugiyamaN. MeV-scale reheating temperature and thermalization of the neutrino background. Phys Rev D (2000) 62:023506. 10.1103/PhysRevD.62.023506
28.
HannestadS. What is the lowest possible reheating temperature?Phys Rev D (2004) 70:043506. 10.1103/PhysRevD.70.043506
29.
IchikawaKKawasakiMTakahashiF. Oscillation effects on thermalization of the neutrinos in the universe with low reheating temperature. Phys Rev D (2005) 72:043522. 10.1103/PhysRevD.72.043522
30.
De BernardisFPaganoLMelchiorriA. New constraints on the reheating temperature of the universe after WMAP-5. Astropart Phys (2008) 30:192–5. 10.1016/j.astropartphys.2008.09.005
31.
ArcadiGDjouadiARaidalM. Dark matter through the Higgs portal. Phys Rept (2020) 842:1–180. 10.1016/j.physrep.2019.11.003
32.
ArcadiGLebedevOPokorskiSTomaT. Real scalar dark matter: relativistic treatment. JHEP (2019) 08:050. 10.1007/JHEP08(2019)050
33.
AprileEAalbersJAgostiniFAlfonsiMAlthueserLAmaroFet alLight dark matter search with ionization signals in XENON1T. Phys Rev Lett (2019) 123:251801. 10.1103/PhysRevLett.123.251801
34.
AalbersJAkeribDAkerlofCAl MusalhiAAlderFAlqahtaniAet alFirst dark matter search results from the LUX-ZEPLIN (LZ) experiment. Phys Rev Lett (2023) 131:041002. 10.1103/PhysRevLett.131.041002
35.
AadGAbbottBAbelingKAbidiSAboulhormaAAbramowiczHet alCombination of searches for invisible decays of the Higgs boson using 139 fb−1 of proton-proton collision data at s=13 TeV collected with the ATLAS experiment. Phys Lett B (2023) 842:137963. 10.1016/j.physletb.2023.137963
36.
CepedaMet alReport from working group 2: Higgs physics at the HL-LHC and HE-LHC. CERN Yellow Rep Monogr (2019) 7:221–584. 10.23731/CYRM-2019-007.221
37.
de BlasJCepedaMD’HondtJEllisRGrojeanCHeinemannBet alHiggs boson studies at future particle colliders. JHEP (2020) 01:139. 10.1007/JHEP01(2020)139
38.
MoroiTRandallL. Wino cold dark matter from anomaly mediated SUSY breaking. Nucl Phys B (2000) 570:455–72. 10.1016/S0550-3213(99)00748-8
39.
AllahverdiRDuttaBMohapatraRNSinhaK. Supersymmetric model for dark matter and baryogenesis motivated by the recent CDMS result. Phys Rev Lett (2013) 111:051302. 10.1103/PhysRevLett.111.051302
40.
AcharyaBSKumarPBobkovKKaneGShaoJWatsonS. Non-thermal dark matter and the moduli problem in string frameworks. JHEP (2008) 06:064. 10.1088/1126-6708/2008/06/064
41.
AcharyaBSKaneGWatsonSKumarP. Nonthermal “WIMP miracle”. Phys Rev D (2009) 80:083529. 10.1103/PhysRevD.80.083529
42.
MoroiTNagaiMTakimotoM. Non-thermal production of wino dark matter via the decay of long-lived particles. JHEP (2013) 07:066. 10.1007/JHEP07(2013)066
43.
AparicioLCicoliMDuttaBMuiaFQuevedoF. Light higgsino dark matter from non-thermal cosmology. JHEP (2016) 11:038. 10.1007/JHEP11(2016)038
44.
ChowdhuryDDudasEDutraMMambriniY. Moduli portal dark matter. Phys Rev D (2019) 99:095028. 10.1103/PhysRevD.99.095028
45.
Riajul HaqueMKpatchaEMaityDMambriniY. Primordial black hole reheating. Phys Rev D (2023) 108:063523. 10.1103/PhysRevD.108.063523
46.
HardyE. Higgs portal dark matter in non-standard cosmological histories. JHEP (2018) 06:043. 10.1007/JHEP06(2018)043
47.
PallisC. Massive particle decay and cold dark matter abundance. Astropart Phys (2004) 21:689–702. 10.1016/j.astropartphys.2004.05.006
48.
DreesMHajkarimF. Dark matter production in an early matter dominated era. JCAP (2018) 02:057. 10.1088/1475-7516/2018/02/057
49.
DreesMHajkarimF. Neutralino dark matter in scenarios with early matter domination. JHEP (2018) 12:042. 10.1007/JHEP12(2018)042
50.
CheungCElorGHallLJKumarP. Origins of hidden sector dark matter I: cosmology. JHEP (2011) 03:042. 10.1007/JHEP03(2011)042
51.
FigueroaDGTaninEH. Ability of LIGO and LISA to probe the equation of state of the early Universe. JCAP (2019) 08:011. 10.1088/1475-7516/2019/08/011
52.
LiuLChenZCHuangQG. Probing the equation of state of the early Universe with pulsar timing arrays. JCAP (2023) 11:071. 10.1088/1475-7516/2023/11/071
53.
CharlesESánchez-CondeMAndersonBCaputoRCuocoADi MauroMet alSensitivity projections for dark matter searches with the fermi large area telescope. Phys Rept (2016) 636:1–46. 10.1016/j.physrep.2016.05.001
54.
McDanielAAjelloMKarwinCMDi MauroMDrlica-WagnerASánchez-CondeMA. Legacy analysis of dark matter annihilation from the Milky Way dwarf spheroidal galaxies with 14 years of Fermi-LAT data. Phys Rev D (2024) 109:063024. 10.1103/PhysRevD.109.063024
55.
GalliSIoccoFBertoneGMelchiorriA. CMB constraints on Dark Matter models with large annihilation cross-section. Phys Rev D (2009) 80:023505. 10.1103/PhysRevD.80.023505
56.
BallesterosGGarciaMAGPierreM. How warm are non-thermal relics? Lyman-α bounds on out-of-equilibrium dark matter. JCAP (2021) 03:101. 10.1088/1475-7516/2021/03/101
Summary
Keywords
dark matter–cosmology, dark matter theory, beyond the standard model interactions, dark matter phenomenology, early universe
Citation
Arcadi G (2024) Thermal and non-thermal DM production in non-standard cosmologies: a mini review. Front. Phys. 12:1425838. doi: 10.3389/fphy.2024.1425838
Received
30 April 2024
Accepted
29 May 2024
Published
26 June 2024
Volume
12 - 2024
Edited by
Behzad Eslam Panah, University of Mazandaran, Iran
Reviewed by
Fazlollah Hajkarim, University of Oklahoma, United States
Updates
Copyright
© 2024 Arcadi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giorgio Arcadi, giorgio.arcadi@unime.it
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.