It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We present a novel UV freeze-in mechanism for enhancing the dark matter abundance in cosmologies where inflation is followed by an epoch dominated by a fluid stiffer than radiation. In such scenarios, even a small radiation abundance, produced for instance by instantaneous preheating effects, will eventually dominate the total energy density of the Universe without the need for a complete inflaton decay. For the sake of concreteness, we focus in non-oscillatory (NO) quintessential inflation models, albeit our treatment is rather general and can be extended to other scenarios. The high temperature of the initial thermal bath together with the absence of subsequent entropy injections into the Standard Model plasma translates into a highly-effective UV DM freeze-in. In particular, we find that, during kination, an enhancement in the DM abundance is generically obtained for any production cross-section not-decreasing with temperature.

In the Starobinsky model of inflation, the observed dark matter abundance can be produced from the direct decay of the inflaton field only in a very narrow spectrum of close-to-conformal scalar fields and spinors of mass $\sim 10^7$ GeV. This spectrum can be, however, significantly broadened in the presence of effective non-renormalizable interactions between the dark and the visible sectors. In particular, we show that UV freeze-in can efficiently generate the right dark matter abundance for a large range of masses spanning from the keV to the PeV scale and arbitrary spin, without significantly altering the heating dynamics. We also consider the contribution of effective interactions to the inflaton decay into dark matter.

We study for the first time the case in which Dark Matter (DM) is made of Feebly Interacting Massive Particles (FIMP) interacting just gravitationally with the standard model particles in an extra-dimensional Randall-Sundrum scenario. We assume that both the dark matter and the standard model are localized in the IR-brane and only interact via gravitational mediators, namely the graviton, the Kaluza-Klein gravitons and the radion. We found that in the early Universe DM could be generated via two main processes: the direct freeze-in and the sequential freeze-in. The regions where the observed DM relic abundance is produced are largely compatible with cosmological and collider bounds.

The existence of dark matter particles that carry phenomenologically relevant self-interaction cross sections mediated by light dark sector states is considered to be severely constrained through a combination of experimental and observational data. The conclusion is based on the assumption of specific dark matter production mechanisms such as thermal freeze-out together with an extrapolation of a standard cosmological history beyond the epoch of primordial nucleosynthesis. In this work, we drop these assumptions and examine the scenario from the perspective of the current firm knowledge we have: results from direct and indirect dark matter searches and cosmological and astrophysical observations, without additional assumptions on dark matter genesis or the thermal state of the very early universe. We show that even in the minimal set-up, where dark matter particles self-interact via a kinetically mixed vector mediator, a significant amount of parameter space remains allowed. Interestingly, however, these parameter regions imply a meta-stable, light mediator, which in turn calls for modified search strategies.

Currently, the standard cosmological model faces some tensions and discrepancies between observations at early and late cosmological time. One of them concerns the well-known $H_0$-tension problem, i.e., a $\sim4.4\sigma$-difference between the early-time estimate and late-time measurements of the Hubble constant, $H_0$. Another puzzling question rests in the cosmological lithium abundance, where again local measurements differ from the one predicted by Big Bang Nucleosynthesis (BBN). In this work, we show that a mechanism of light dark matter production might hold the answer for these questions. If dark matter particles are sufficiently light and a fraction of them was produced non-thermally in association with photons, this mechanism has precisely what is needed to destroy Lithium without spoiling other BBN predictions. Besides, it produces enough radiation that leads to a larger $H_0$ value, reconciling early and late-time measurements of the Hubble expansion rate without leaving sizable spectral distortions in the Cosmic Microwave Background spectrum.

A notable feature of UV freeze-in is that the relic density is strongly dependent on the highest temperatures of the thermal bath, and a common assumption is that the relevant 'highest temperature' should be the reheating temperature after inflation $T_\text{RH}$. However, the temperature of the thermal bath can be significantly higher in certain scenarios, reaching a value denoted T max , a fact which is only apparent away from the instantaneous decay approximation. Interestingly, it has been shown that if the operators are of sufficiently high mass dimension then the dark matter abundance can be enhanced by a 'boost factor' depending on ($T_\text{max}/T_\text{RH}$) relative to naive estimates assuming instantaneous reheating. We highlight here that in non-standard cosmological histories the critical mass dimension of the operator above at which the instantaneous decay approximation breaks down, and the exponent of the boost factor, depend on the equation of state $\omega$ prior to reheating. We highlight four examples in which the dark matter abundance receives a significant enhancement in the context of gravitino dark matter, the moduli portal, the Higgs portal, and the spin-2 portal (as might arise in bimetric gravity models). We comment on the transition from kination domination to radiation domination as a motivated example of non-standard cosmologies.

Assuming that inflation is followed by a phase where the energy density of the Universe is dominated by a component with a general equation of state, we evaluate the spectrum of primordial gravitational waves induced in the post-inflationary Universe. We show that if the energy density of the Universe is dominated by a component $\phi$ before Big Bang nucleosynthesis, its equation of state could be constrained by gravitational wave experiments depending on the ratio of energy densities of $\phi$ and radiation, and also the temperature at the end of the $\phi$ dominated era. Also, we discuss the impact of scale dependence of tensor modes on the primordial gravitational wave spectrum during the $\phi$-domination. These models are motivated by beyond Standard Model physics and scenarios for non-thermal production of dark matter in the early Universe. We also constrain the parameter space of the tensor spectral index and the tensor-to-scalar ratio, using the experimental limits from gravitational wave experiments.

If cosmic inflation was driven by an electrically neutral scalar field stable on cosmological time scales, the field necessarily constitutes all or part of dark matter (DM). We study this possibility in a scenario where the inflaton field $s$ resides in a hidden sector, which is coupled to the Standard Model sector through the Higgs portal $\lambda_{hs} s^2\mathcal{H}^\dagger\mathcal{H}$ and non-minimally to gravity via $\xi_s s^2 R$. We study scenarios where the field $s$ first drives inflation, then reheats the Universe, and later constitutes all DM. We consider two benchmark scenarios where the DM abundance is generated either by production during reheating or via non-thermal freeze-in. In both cases, we take into account all production channels relevant for DM in the mass range from keV to PeV scale. On the inflationary side, we compare the dynamics and the relevant observables in two different but well-motivated theories of gravity (metric and Palatini), discuss multifield effects in case both fields ($s$ and $h$) were dynamical during inflation, and take into account the non-perturbative nature of particle production during reheating. We find that, depending on the initial conditions for inflation, couplings and the DM mass, the scenario works well especially for large DM masses, $10^2$ GeV$\lesssim m_{s}\lesssim 10^6$ GeV, although there are also small observationally allowed windows at the keV and MeV scales. We discuss how the model can be tested through astrophysical observations.

We study production of self-interacting dark matter (DM) during an early matter-dominated phase. As a benchmark scenario, we consider a model where the DM consists of singlet scalar particles coupled to the visible Standard Model (SM) sector via the Higgs portal. We consider scenarios where the initial DM abundance is set by either the usual thermal freeze-out or an alternative freeze-in mechanism, where DM was never in thermal equilibrium with the SM sector. For the first time, we take the effect of self-interactions within the hidden sector into account in determining the DM abundance, reminiscent to the Strongly Interacting Massive Particle (SIMP) scenario. In all cases, the number density of DM may change considerably compared to the standard radiation-dominated case, having important observational and experimental ramifications.

We study production of dark matter (DM) in models with a non-standard expansion history. We consider both freeze-out and freeze-in mechanisms for producing the observed DM abundance in a model where the DM consists of scalar singlet particles coupled to the Standard Model sector via the Higgs portal. We show that a non-standard expansion phase can lead to a significant change in the DM abundance and therefore to observational ramifications. For example, for DM freeze-in the required portal coupling can be much larger, whereas for DM freeze-out much smaller values become allowed. We evaluate the relevant constraints and discuss prospects for direct detection of such DM.

The existence of Dark Matter (DM) has been well established from various cosmological and astrophysical evidences. However, the particle properties of DM are largely undetermined and attempts to probe its interactions with the Standard Model (SM) particles have, so far, not met with any success. The stringent constraints on the DM-SM interactions, while does not exclude the standard lore of producing weakly massive interacting particle DM candidates through thermal freeze-out mechanism in its entirety, have certainly cast shadow on the same. In this work, we consider non-thermal production of DM within a simple extension of the SM including an inflaton field and a scalar DM candidate. Assuming negligible interactions between the SM particles and the DM, we study the production of the latter at the end of inflation, during the (p)reheating epoch. In this context, we explore the role of DM self-interactions and its interaction with the inflaton field, and find that DM can be over produced in a significant region of the parameter space. We further demonstrate that large self-interaction of the DM can suppress its abundance during preheating and to a certain extent helps to achieve the observed relic abundance via cannibalization.

We consider a scenario in which the inflaton $\phi$ is a pseudoscalar field non-minimally coupled to gravity through a term of the form ${\cal X} R \phi^2$. The pseudoscalar is also coupled to a $U(1)$ gauge field (or an ensemble of ${\cal N}$ gauge fields) through an axial coupling of the form $\phi F \tilde{F}$. After M. M. Anber and L. Sorbo, Phys. Rev. D 81, 043534 (2010), Ref. [1], it is well known that this axial coupling leads to a production of gauge particles which acts as a friction term in the dynamics of the inflaton, producing a slow-roll regime even in presence of a steep potential. A remarkable result in this scenario, is that the spectrum of the chiral gravitational waves sourced by the scalar-gauge field interplay can be enhanced due to the non-minimal coupling with gravity, leading to measurable signatures, while maintaining agreement with current observational constraints on $n_s$ and $r$. The inclusion of non-minimal coupling could be helpful to alleviate tensions with non-Gaussianity bounds in models including axial couplings.

Collider, space, and Earth based experiments are now able to probe several extensions of the Standard Model of particle physics which provide viable dark matter candidates. Direct and indirect dark matter searches rely on inputs of astrophysical nature, such as the local dark matter density or the shape of the dark matter profile in the target in object. The determination of these quantities is highly affected by astrophysical uncertainties. The latter, especially those for our own Galaxy, are ill-known, and often not fully accounted for when analyzing the phenomenology of particle physics models. In this paper we present a systematic, quantitative estimate of how astrophysical uncertainties on Galactic quantities (such as the local galactocentric distance, circular velocity, or the morphology of the stellar disk and bulge) propagate to the determination of the phenomenology of particle physics models, thus eventually affecting the determination of new physics parameters. We present results in the context of two specific extensions of the Standard Model (the Singlet Scalar and the Inert Doublet) that we adopt as case studies for their simplicity in illustrating the magnitude and impact of such uncertainties on the parameter space of the particle physics model itself. Our findings point toward very relevant effects of current Galactic uncertainties on the determination of particle physics parameters, and urge a systematic estimate of such uncertainties in more complex scenarios, in order to achieve constraints on the determination of new physics that realistically include all known uncertainties.

We generalize models invoking a spin-2 particle as a mediator between the dark sector and the Standard Model. We show that a massive spin-2 messenger can efficiently play the role of a portal between the two sectors. The dark matter is then produced via a freeze-in mechanism during the reheating epoch. In a large part of the parameter space, production through the exchange of a massive spin-2 mediator dominates over processes involving a graviton with Planck suppressed couplings. We perform a systematic analysis of such models for different values of the spin-2 mass relative to the maximum and the final temperature attained at reheating.

We present an overview of scenarios where the observed Dark Matter (DM) abundance consists of Feebly Interacting Massive Particles (FIMPs), produced non-thermally by the so-called freeze-in mechanism. In contrast to the usual freeze-out scenario, frozen-in FIMP DM interacts very weakly with the particles in the visible sector and never attained thermal equilibrium with the baryon-photon fluid in the early Universe. Instead of being determined by its annihilation strength, the DM abundance depends on the decay and annihilation strengths of particles in equilibrium with the baryon-photon fluid, as well as couplings in the DM sector. This makes frozen-in DM very difficult but not impossible to test. In this review, we present the freeze-in mechanism and its variations considered in the literature (dark freeze-out and reannihilation), compare them to the standard DM freeze-out scenario, discuss several aspects of model building, and pay particular attention to observational properties and general testability of such feebly interacting DM.

Dark matter (DM) halos have long been known to be triaxial, but in studies of possible annihilation and decay signals they are often treated as approximately spherical. In this work, we examine the asymmetry of potential indirect detection signals of DM annihilation and decay, exploiting the large statistics of the hydrodynamic simulation Illustris. We carefully investigate the effects of the baryons on the sphericity of annihilation and decay signals for both the case where the observer is at 8.5 kpc from the center of the halo (exemplified in the case of Milky Way-like halos), and for an observer situated well outside the halo. In the case of Galactic signals, we find that both annihilation and decay signals are expected to be quite symmetric, with axis ratios very different from 1 occurring rarely. In the case of extragalactic signals, while decay signals are still preferentially spherical, the axis ratio for annihilation signals has a much flatter distribution, with elongated profiles appearing frequently. Many of these elongated profiles are due to large subhalos and/or recent mergers. Comparing to gamma-ray emission from the Milky Way and X-ray maps of clusters, we find that the gamma-ray background appears less spherical/more elongated than the expected DM signal from the large majority of halos, and the Galactic gamma ray excess appears very spherical, while the X- ray data would be difficult to distinguish from a DM signal by elongation/sphericity measurements alone.

Although commonly assumed to be spherical, dark matter halos are predicted to be non-spherical by N-body simulations and their asphericity has a potential impact on the systematic uncertainties in dark matter searches. The evaluation of these uncertainties is the main aim of this work, where we study the impact of aspherical dark matter density distributions in Milky-Way-like halos on direct and indirect searches. Using data from the large N-body cosmological simulation Bolshoi, we perform a statistical analysis and quantify the systematic uncertainties on the determination of local dark matter density and the so-called $J$ factors for dark matter annihilations and decays from the galactic center. We find that, due to our ignorance about the extent of the non-sphericity of the Milky Way dark matter halo, systematic uncertainties can be as large as 35%, within the 95% most probable region, for a spherically averaged value for the local density of 0.3-0.4 GeV/cm$^3$. Similarly, systematic uncertainties on the $J$ factors evaluated around the galactic center can be as large as 10% and 15%, within the 95% most probable region, for dark matter annihilations and decays, respectively.

In this Letter, we study the gamma-ray signatures subsequent to the production of a Higgs boson in space by dark matter annihilations. We investigate the cases where the Higgs boson is produced at rest or slightly boosted and show that such configurations can produce characteristic bumps in the gamma-ray data. These results are relevant in the case of the Standard Model-like Higgs boson provided that the dark matter mass is about 63 GeV, 109 GeV or 126 GeV, but can be generalised to any other Higgs boson masses. Here, we point out that it may be worth looking for a 63 GeV line since it could be the signature of the decay of a Standard Model-like Higgs boson produced in space, as in the case of a di-Higgs final state if m_DM ~ 126 GeV. We show that one can set generic constraints on the Higgs boson production rates using its decay properties. In particular, using the Fermi-LAT data from the galactic center, we find that the dark matter annihilation cross section into gamma + a Standard Model-like Higgs boson produced at rest or near rest cannot exceed <sigma v> ~ a few 10^-25 cm^3/s or <sigma v> ~ a few 10^-27 cm^3/s respectively, providing us with information on the Higgs coupling to the dark matter particle. We conclude that Higgs bosons can indeed be used as messengers to explore the dark matter mass range.

A possible connection between the abundances of baryonic and dark matter (DM) has been explored so far mostly in the context of the so-called asymmetric DM. Recently, a very different mechanism, dubbed "WIMPy baryogenesis", has been proposed to relate the baryon asymmetry to DM annihilation. The DM candidate is a weakly interacting massive particle (WIMP), and the usual WIMP scenario is slightly extended to accommodate baryogenesis, which is accomplished around the time of DM freeze-out. We construct an effective field theory that encompasses a quite general class of models which implement the WIMPy baryogenesis. Under some reasonable, simplifying assumptions, we show that a good portion of the parameter space is allowed for these models, after experimental constraints are taken into account. Bounds from the LHC require that the WIMP be heavier than 400 GeV.

We study the abilities of the Fermi-LAT instrument on board of the Fermi mission to simultaneously constrain the Milky Way dark matter density profile and some dark matter particle properties, as annihilation cross section, mass and branching ratio into dominant annihilation channels. A single dark matter density profile is commonly assumed to determine the capabilities of gamma-ray experiments to extract dark matter properties or to set limits on them. However, our knowledge of the Milky Way halo is far from perfect, and thus in general, the obtained results are too optimistic. Here, we study the effect these astrophysical uncertainties would have on the determination of dark matter particle properties and conversely, we show how gamma-ray searches could also be used to learn about the structure of the Milky Way halo, as a complementary tool to other type of observational data that study the gravitational effect caused by the presence of dark matter. In addition, we also show how these results would improve if external information on the annihilation cross section and on the local dark matter density were included and compare our results with the predictions from numerical simulations.

We study the capabilities of the Fermi-LAT instrument on board of the Fermi mission to constrain particle dark matter properties, as annihilation cross section, mass and branching ratio into dominant annihilation channels, with gamma-ray observations from the galactic center. Besides the prompt gamma-ray flux, we also take into account the contribution from the electrons/positrons produced in dark matter annihilations to the gamma-ray signal via inverse Compton scattering off the interstellar photon background, which turns out to be crucial in the case of dark matter annihilations into mu+mu- and e+e- pairs. We study the signal dependence on different parameters like the region of observation, the density profile, the assumptions for the dark matter model and the uncertainties in the propagation model. We also show the effect of the inclusion of a 20% systematic uncertainty in the gamma-ray background. If Fermi-LAT is able to distinguish a possible dark matter signal from the large gamma-ray background, we show that for dark matter masses below ~200 GeV, Fermi-LAT will likely be able to determine dark matter properties with good accuracy.

We study the capabilities of the Fermi-LAT instrument for identifying particle Dark Matter properties as mass, annihilation cross section and annihilation channels, with gamma-ray observations from the Galactic Center. For the potential Dark Matter signal, besides the prompt gamma-ray flux produced in Dark Matter annihilations, we also take into account the flux produced by inverse Compton scattering of the electrons and positrons generated in Dark Matter annihilations off the interstellar photon background. We show that the addition of this contribution is crucial in the case of annihilations into e^+ e^- and mu^+ mu^- pairs. In addition to the diffuse galactic and extragalactic background, we also consider the full catalog of high-energy gamma-ray point sources detected by Fermi. The impact of the degeneracies between the different Dark Matter annihilation channels has been studied. We find that for Dark Matter masses below ~200 GeV and for typical thermal annihilation cross sections, it will be possible to obtain stringent bounds on the Dark Matter properties.

The addition of non-renormalizable terms involving the Higgs fields to the MSSM (BMSSM) ameliorates the little hierarchy problem of the MSSM. For neutralino dark matter, new regions for which the relic abundance of the LSP is consistent with WMAP (as the bulk region and the stop coannihilation region) are now permitted. In this framework, we analyze in detail the direct dark matter detection prospects in a Xenon-like experiment. On the other hand, we study the capability of detecting gamma-rays, antiprotons and positrons produced in the annihilation of neutralino LSPs in the Fermi and oncoming AMS-02 experiments.

We explore the dark matter detection prospects in the Minimal Supersymmetric Standard Model in the scenario where the scalar partners of the fermions and the Higgs particles (except for the Standard-Model-like one) are assumed to be very heavy and are removed from the low-energy spectrum. We analyse the neutralino LSP ($\chi_1^0$) in scenarios where the gaugino mass parameters are universal at the GUT scale and also the case where they are non-universal. This analysis is carried out in the framework of a Xenon-like 100 kg experiment. In general, an important fraction of the parameter space giving rise to the dark matter relic density measured by WMAP can be probed and excluded in the case of not detecting any WIMP. In the opposite case, once a WIMP signal has been found, we show that for a light $\chi_1^0$ which is a higgsino-gaugino mixture it is possible to reconstruct efficiently the mass and the scattering cross-section of the neutralino LSP. Moreover, we show that it is also feasible to put strong constraints over some of the parameters of the Lagrangian, e.g. the higgsino and the gaugino mass parameters.

We study the possibility of identifying dark matter properties from XENON-like 100 kg experiments and the GLAST satellite mission. We show that whereas direct detection experiments will probe efficiently light WIMPs, given a positive detection (at the 10% level for $m_{\chi} \lesssim 50$ GeV), GLAST will be able to confirm and even increase the precision in the case of a NFW profile, for a WIMP-nucleon cross-section $\sigma_{\chi-p} \lesssim 10^{-8}$ pb. We also predict the rate of production of a WIMP in the next generation of colliders (ILC), and compare their sensitivity to the WIMP mass with the XENON and GLAST projects.

We study the possibility of identifying dark matter properties from direct (XENON100) and indirect (GLAST) detection experiments. In the same way, we examine the perspectives given by the next generation of colliders (ILC). All this analysis is done following a model-independent approach. We have shown that the three detection techniques can act in a highly complementary way. whereas direct detection experiments will probe efficiently light WIMPs, given a positive detection (at the 10% level for $m_{\chi} \sim 50$ GeV), GLAST will be able to confirm and even increase the precision in the case of a NFW profile, for a WIMP-nucleon cross-section $\sigma_{\chi-p}<10^{-8}$ pb. However, for heavier WIMP ($\sim 175$ GeV), the ILC will lead the reconstruction of the mass.