On general features of warm dark matter with reduced relativistic gas

TL;DR

This paper uses the Reduced Relativistic Gas model to analyze the general features of warm dark matter, providing a simplified, model-independent way to understand its cosmological effects and observational signatures.

Contribution

It introduces the RRG approach as a tool to study WDM's features and effects without detailed particle physics assumptions, highlighting its utility in cosmological modeling.

Findings

RRG reproduces the small-scale matter power spectrum with less than 1% error.

High warmness parameter values erase the radiation-dominated epoch.

Low warmness parameter values are consistent with CMB observations.

Abstract

Reduced Relativistic Gas (RRG) is a useful approach to describe the warm dark matter (WDM) or the warmness of baryonic matter in the approximation when the interaction between the particles is irrelevant. The use of Maxwell distribution leads to the complicated equation of state of the J\"{u}ttner model of relativistic ideal gas. The RRG enables one to reproduce the same physical situation but in a much simpler form. For this reason RRG can be a useful tool for the theories with some sort of a "new Physics". On the other hand, even without the qualitatively new physical implementations, the RRG can be useful to describe the general features of WDM in a model-independent way. In this sense one can see, in particular, to which extent the cosmological manifestations of WDM may be dependent on its Particle Physics background. In the present work RRG is used as a complementary approach to derive the main observational exponents for the WDM in a model-independent way. The only assumption concerns a non-negligible velocity $v$ for dark matter particles which is parameterized by the warmness parameter $b$. The relatively high values of $b$ ( $b^2\gtrsim 10^{-6}$) erase the radiation (photons and neutrinos) dominated epoch and cause an early warm matter domination after inflation. Furthermore, RRG approach enables one to quantify the lack of power in linear matter spectrum at small scales and in particular, reproduces the relative transfer function commonly used in context of WDM with accuracy of $\lesssim 1\%$. A warmness with $b^2\lesssim 10^{-6}$ (equivalent to $v\lesssim 300 km/s$) does not alter significantly the CMB power spectrum and is in agreement with the background observational tests.

Tables

Table 1. Table 1: Result of the statistical analysis for the joint analysis using SNIa (JLA), H 0 ​ ( z ) subscript 𝐻 0 𝑧 H_{0}(z) and BAO data sets. The best-fit result was obtained when χ m ​ i ​ n 2 = 685.3 superscript subscript 𝜒 𝑚 𝑖 𝑛 2 685.3 \chi_{min}^{2}=685.3 .

Parameter	Best-fit	95% lower	95% upper
$h$	$0.7386$	$0.6923$	$0.7833$
${\Omega }_{w0}$	$0.2607$	$0.2237$	$0.3025$
$b$	$6.425\times {10}^{-6}$	$0^{\ast }$	$2.163\times {10}^{-5}$
$\alpha$	$0.1417$	$0.1277$	$0.1545$
$\beta$	$3.094$	$2.936$	$3.265$
$M$	$-18.92$	$-19.08$	$-18.78$
${\Delta }_M$	$-0.07238$	$-0.1181$	$-0.02448$
${\Omega }_{m0}$	$0.3107$	$0.2737$	$0.3525$

Table 2. Table 2: Result of the statistical analysis for the joint analysis using SNIa (JLA), H 0 ​ ( z ) subscript 𝐻 0 𝑧 H_{0}(z) and BAO data sets. The best-fit result was obtained when χ m ​ i ​ n 2 = 685.3 superscript subscript 𝜒 𝑚 𝑖 𝑛 2 685.3 \chi_{min}^{2}=685.3 .

$b^2$	$\alpha$	$\nu$
${10}^{-10}$	$2.450$	$2.12$
${10}^{-11}$	$0.510$	$1.64$
${10}^{-12}$	$0.350$	$1.23$
${10}^{-13}$	$0.092$	$1.15$
${10}^{-14}$	$0.028$	$1.10$
${10}^{-15}$	$0.005$	$0.92$

Equations

$$p_w=\frac{{\rho }_w}{3}\left[1-{\left(\frac{m{c}^2}{ϵ}\right)}^2\right],$$

$$ϵ=\frac{m{c}^2}{\sqrt{1-{\beta }^2}}.$$

$$p_w=\frac{{\rho }_w}{3}\left[1-{\left(\frac{{\rho }_c}{{\rho }_w}\right)}^2\right],$$

$${\rho }_w(a)={\rho }_{w0}{a}^{-3}\sqrt{\frac{1+b^2{a}^{-2}}{1+b^2}},b=\frac{\beta }{\sqrt{1-{\beta }^2}}.$$

$$w(a)=\frac{p_w}{{\rho }_w}=\frac{1}{3}-\frac{a^2}{3(a^2+b^2)}.$$

$$H^2=H_0^2\left({\Omega }_{\Lambda 0}+\frac{{\Omega }_{w0}}{a^3}\sqrt{\frac{1+b^2{a}^{-2}}{1+b^2}}+\frac{{\Omega }_{b0}}{a^3}+\frac{{\Omega }_{r0}}{a^4}\right).$$

$$\underset{a\to 0}{\lim }\frac{{\Omega }_r\left(a\right)}{{\Omega }_w\left(a\right)}>1\Rightarrow b^2<\frac{{\Omega }_{r0}^2}{{\Omega }_{w0}^2-{\Omega }_{r0}^2}.$$

$${\Omega }_r(z_{eq})={\Omega }_w(z_{eq})+{\Omega }_b(z_{eq}).$$

$$1+z_{eq}={\left(\frac{{\Omega }_{b0}{\Omega }_{r0}}{{\Omega }_{b0}^2-{\Omega }_{w0}^2}-\frac{\sqrt{(1+b^2){\Omega }_{w0}^2{\Omega }_{r0}^2+b^2{\Omega }_{w0}^2[(1+b^2){\Omega }_{b0}^2-{\Omega }_{w0}^2]}}{(1+b^2){\Omega }_{b0}^2-{\Omega }_{w0}^2}\right)}^{-1}.$$

$${\dot{\delta }}_w+(1+w)\left({\theta }_w+\frac{\dot{h}}{2}\right)+3\mathcal{H}\left(c_s^2-w\right){\delta }_w+9{\mathcal{H}}^2\left(c_s^2-w\right)(1+w)\frac{{\theta }_w}{k^2}+3\mathcal{H}\dot{w}\frac{{\theta }_w}{k^2}=0,$$

$${\dot{\theta }}_w+\mathcal{H}(1-3c_s^2){\theta }_w-\frac{k^2{c}_s^2}{1+w}{\delta }_w=0.$$

$$\delta p_w=c_s^2\delta {\rho }_w,$$

$$\dot{w}=-\frac{\mathcal{H}}{3}\frac{a^2}{(a^2+b^2)}\text{mbox}and{c}_s^2=w-\frac{\dot{w}}{3\mathcal{H}(1+w)}.$$

$${\dot{\delta }}_w+\frac{4}{3}{\theta }_w+\frac{2}{3}\dot{h}=0\text{mbox}and{\dot{\theta }}_w-\frac{k^2}{4}{\delta }_w=0.$$

$${\chi }^2(\theta )=\Delta y(\theta {)}^T{\mathbf{C}}^{-1}\Delta y(\theta ),$$

$${\delta }_m\equiv \frac{\delta {\rho }_m}{{\rho }_m}=\frac{\delta {\rho }_w+\delta {\rho }_b}{{\rho }_w+{\rho }_b}.$$

$${\delta }_m=\frac{{\tilde{\Omega }}_w(a){\delta }_w+{\Omega }_{b0}{\delta }_b}{{\tilde{\Omega }}_w(a)+{\Omega }_{b0}},{\tilde{\Omega }}_{w0}(a)={\Omega }_{w0}\sqrt{\frac{1+b^2{a}^{-2}}{1+b^2}}.$$

$$T(k)\equiv \frac{{\delta }_m(k,z=0){\delta }_m(0,z=0)}{{\delta }_m(k,z\to \infty ){\delta }_m(0,z\to \infty )}.$$

$$\Delta D_l=D_l^{\Lambda CDM}-D_l^{\Lambda WDM},D_l=\frac{l(l+1)C_l}{2\pi }$$

$$\bar{T}(k)\equiv {\left[\frac{P_{\Lambda \mathrm{WDM}}(k)}{P_{\Lambda \mathrm{CDM}}(k)}\right]}^{1/2},$$

$$\bar{T}(k)=[1+(\alpha k{)}^{2\nu }{]}^{-5/\nu },$$

$$\alpha =0.049{\left(\frac{m_w}{1keV}\right)}^{-1.11}{\left(\frac{{\Omega }_{w0}}{0.25}\right)}^{0.11}{\left(\frac{0.01H_0}{0.7}\right)}^{1.22}.$$

$$b^2=\left(1\times 10^{-15},2\times 10^{-15},\dots ,9\times 10^{-12},1\times 10^{-11}\right).$$

$${\chi }^2={\left(\bar{T}{\left(k\right)}^{th.}-\bar{T}{\left(k\right)}^{num.}\right)}^2,$$

$$m_w=4.65\cdot 10^{-6}{\left(b^2\right)}^{-2/5}keV.$$

Full text

I Introduction

In the last decades cosmological observations provided numerous evidence for the two dark components nominated dark matter (DM) and dark energy (DE), which are responsible for $\sim 96\%$ of the content of the universe. In particular, the confirmation of existence of these two dark components come from the measurements of the luminosity redshift of type Ia supernova Riess:1998cb ; Perlmutter:1998np , baryon acoustic oscillations Tegmark:2003uf , anisotropies of the cosmic microwave background (CMB)Aghanim:2015xee ; Ade:2015xua and other observations Moresco:2012jh . The standard interpretation suggests that DE is necessary to accelerate the expansion of the universe. On the other hand the DM has non-baryonic nature and is important, in particular, to describe the formation of cosmic structure. The standard cosmology, $\Lambda$CDM model, assumes that the DE is a cosmological constant, and regards DM as a non-relativistic matter with negligible pressure (cold dark matter). $\Lambda$CDM provides an excellent agreement with most of the data (see, e.g., Bergstrom:2000pn for a general review), however this agreement is not perfect due to the tensions with part of the observational data (see for example Buchert:2015wwr ). In part due to these difficulties, some alternative models have been proposed and studied as possible DE and DM candidates (see for example Caldwell:1997ii ; Liddle:1998xm ; Peebles:2002gy ; Sahni:2002dx ; Lue:2005ya ). Let us note that some of these alternative models aim to describe fluids that replace both DM and DE (see for example Kamenshchik:2001cp ; Fabris:2001tm ; Bilic:2001cg ; HipolitoRicaldi:2009je ; Zimdahl:2011mb ; Marttens:2017njo ) or describe interaction between DE and DM Billyard:2000bh ; Amendola:1999er ; Honorez:2010rr ; Zimdahl:2013yak ; Wang:2016lxa ; Funo:2014poa ; Marttens:2014yja .

Some of the mentioned $\Lambda$CDM difficulties are related with the choice of the cold dark matter (CDM) paradigm Warren:2005ey . For instance, at small scales the issues such as the missing satellites problem Klypin:1999uc , core/cusp problem deBlok:2009sp , and the Too big to fail problem BoylanKolchin:2011de , can be alleviated by assuming that the DM is not completely cold. In contrast to the CDM, the Hot Dark Matter (HDM) scenario implies that the free streaming due to a thermal motion of particles is important to suppress structure formation at small scales. Nevertheless this scenario was ruled out and opens the space for the Warm Darm Matter (WDM) scenario. The main feature of WDM models is that thermal velocities of the DM particles are not so high as in the HDM scenario and, on the other hand, not negligible like in the CDM scenario. Typically, the WDM models assume that it is composed by particles of mass about $keV$ instead of $GeV$ which is “typical” for CDM and $eV$ which is the standard case for the HDM. The standard approach to explore the possible warmness of DM and its consequences for structure formation are based on to solution of the Hierarchy Bolztmann equation, taking into account the specific properties of the given WDM candidate Bond:1983hb ; Hannestad:2000gt ; Viel:2005qj ; Bode:2000gq ; Barkana:2001gr ; Ma:1995ey 666One has to remember that equations for DM are always coupled to the Bolztmann equations for other components of the universe.. For example, relation between mass and warmness for each WDM candidate comes from the particle physics arguments. This is in fact very good, because the ultimate knowledge of the DM nature may be achieved only within the particle physics and, more concrete, by means of laboratory experiments.

Until the moment when the DM will be detected in the laboratory experiments, one can always assume that the properties of DM derived within a particle physics models may be violated by some the qualitatively new scenarios for the DM, which can be never ruled out completely Bergstrom:2000pn . From this perspective, it is useful to develop also model-independent approaches to investigate the cosmological features of a WDM. In the present work we will explore the consequences and impacts of warmness in the process of structure formation and CMB anisotropies, but using a model-independent approach which is based on the RRG approximation. The RRG is a model of ideal gas of relativistic particles, which has a very simple equation of state. This nice property is due to the main assumption – that the particles of the ideal gas have non-negligible but equal thermal velocities. Regardless of this simplicity, the model has long history which started in a glorious way. The RRG equation of state was first introduced by A.D. Sakharov in order to explore the acoustic features of Cosmic Microwave Background (CMB) in the early universe sakharov1966initial . Using this model Sakharov predicted the existence of oscillations in CMB temperature spectra long before its observational discovery (see Grishchuk:2011wk for the historical review).

Recently, RRG was reinvented in deBerredoPeixoto:2004ux , where the derivation of its equation of state was first presented explicitly. The simplicity of the equation of state is due to the assumption that all particles of relativistic gas have equal kinetic energies, i.e., equal velocities. Therefore RRG is a reduced version of well-known Jüttner model of relativistic ideal gas Juttnerttner:1911 ; pauli1981theory . A comparison between the equations of state of the relativistic ideal gas and RRG shows that the difference does not exceed $2.5\%$ even in the low-energy region deBerredoPeixoto:2004ux and becomes completely negligible at higher energied. Further considerations have shown that RRG model enables one to achieve a simple and reliable description of the matter warmness in cosmology. In Ref. Fabris:2008qg RRG was used to decribe WDM and its perturbations were compared with the Large Scale Structure data. Furthermore, the general analytic solutions for the several background cosmological models involving RRG were discussed in Ref. Medeiros:2012ud .

The RRG was successfully used in deBerredoPeixoto:2004ux ; Fabris:2008qg as an interpolation between radiation and dark matter eras in cosmology. An upper bound on the warmness coming from RRG Fabris:2008qg ; Fabris:2011am is very close to the one obtained from much more complicated analysis based on a complete WDM treatment, based on the Boltzmann equation. This standard approach requires specifying the nature of the particle physics candidate for the WDM contents Bond:1983hb ; Hannestad:2000gt ; Viel:2005qj ; Bode:2000gq ; Barkana:2001gr ; Ma:1995ey , while the approach based on RRG requires only one parameter, that is the warmness of DM. In this sense RRG represents a really useful tool for exploring WDM cosmology without specifying a particular candidate for the WDM. Such a model may be helpful for better understanding of the model-dependence or independence of the cosmological features of WDM.

The main goal of the present work is to take advantage of the RRG and its analytical solutions for the background cosmology and apply it to WDM, instead of considering full set of WDM hierarchical Boltzmann equations. The RRG enables one to make greater part of considerations analytically and hence provide better qualitative understanding of the results. In this way we consider the formation of large-scale structure and the problem of CMB anisotropy in a model-independent manner. Following Fabris:2008qg , we shall establish the bounds for the thermal velocities of the WDM particles in a general way. With this objective in mind we consider the model of the spatially flat Friedmann-Roberston-Walker universe filled by radiation 777Of course, WDM has a radiation behavior in early epochs but their physical processes are different than photons or neutrinos. For this reason we will diferenciate in all paper WDM in early stages from “standard“ radiation (photons and neutrinos). and RRG, representing WDM. Gravity is described by the general relativity, with the cosmological constant representing DE. We shall refer to this model as to $\Lambda$WDM. All the perturbative treatments will be performed at the linear order only, and using the normalization with the scale factor at present $a_{0}=1$. With these notations, the WDM space of parameters is reduced by using the most recent data from SNIa, $H(z)$ and BAO.

The paper is organized as follows. Sec. II describes the dynamics of the WDM in the framework of RRG, both at the background and perturbative levels. It is shown that high level of warmness may erase ”standard” radiation era from the cosmic history. Starting from this point one can establish an upper bound for the velocity of the RRG particles, which preserves the “standard“ primordial scenario for the universe. This bound is used as a physical prior in the consideration of Sec. III, devoted to the statistical analysis using the background data. In this framework we reduce the space of parameters for WDM and use this reduced space in the consequent analysis. At the next stage the CAMB code is modified and used to quantify the relation between the DM warmness and the total matter density contrast, linear matter power spectrum and CMB power spectrum. We show that the RRG is capable to reproduce the main feature of the WDM, i.e., the suppression of matter over-densities at small scales. Furthermore, in Sec. IV we discuss proprierties of thermal relics via RRG. Finally, Sec. V includes discussions and conclusions.

II A description for a warm dark matter fluid

Let us start with the background notions. The reader can consult deBerredoPeixoto:2004ux ; Fabris:2008qg or recent Reis:2017bjf for further details.

In the RRG approach WDM is treated as an approximation of a Maxwell-distributed ideal gas formed by massive particles. All these particles have equal kinetic energies, or equal velocity $\beta=v/c$ deBerredoPeixoto:2004ux ($c$ is the light speed). This leads to the following relation between WDM pressure $p_{w}$ and WDM energy density $\rho_{w}$,

[TABLE]

where $m$ is the WDM particle mass, and $\epsilon$ is the kinetic energy of each particle of the system which is given by,

[TABLE]

Here $\rho_{c}$ is introduced as a notation for the rest energy density, i.e., the energy density for the $v=0$ case. Thus, $\rho_{c}=\rho_{c0}\,a^{-3}=n\,mc^{2}$, where $n$ is the number density $a=a(t)$ is the scale factor of the metric. Using this relation, Eq. (1) can be cast into the form

[TABLE]

which can be regarded as equation of state (EoS) of the WDM fluid.

Using Eq. (3) in the energy conservation relation, the solution for $\rho_{w}$ has the form

[TABLE]

Thus, $b$ parameter measures velocity and warmness of the WDM particles at present. In the limit $v\ll c$ we have $b\approx v/c$. Note also that for $b=0$ the CDM case is recovered. Combining Eqs. (3) and (4), one can find a posteriori state parameter,

[TABLE]

Here we called this term as a state parameter a posteriori because the ”natural” EoS for the RRG description, given by Eq. (3), implicitly depends on the scale factor and on the nowadays WDM energy density. However after the integration of continuity equation it is possible to write the state parameter that depends only on the scale factor. This form will prove useful in the perturbative analysis.

In what follows we consider the model with cosmological constant, which does not agglomerate, WDM described by a RRG, baryons and radiation. All of them are assumed being interacting gravitationally and only photons and baryons also interacting via Thomson scattering before recombination. In this situation Hubble rate takes the form

[TABLE]

In the last equation $\Omega_{x0}$ (with $x=\Lambda,w,b$ and $r$) is the value of the DE, WDM, baryons and radiation density parameters at present, while $\Omega_{\Lambda 0}=1-\Omega_{w0}-\Omega_{b0}-\Omega_{r0}$ since we deal with a spatially flat universe. It is easy to see that the expressions (1), (3), (4) and (5) interpolate between the dust at $b\to 0$ and radiation at $b\to\infty$ cases. Because of this interpolation feature, RRG can be used to investigate the cosmological consequences of the transition between epochs of radiation and dust sakharov1966initial ; Grishchuk:2011wk ; deBerredoPeixoto:2004ux .

One can note that the WDM with EoS (3) has a remarkable consequence at early times, when RRG becomes very close to radiation. This feature could cause an early warm matter domination and erase the ”standard” radiation dominated epoch. In order to ensure the existence of a “standard” radiation dominated era we must impose that in the very early universe the radiation energy density is bigger than WDM energy density888We must emphasize that even though in a primordial universe RRG behaves like radiation at background, this is not true at perturbative level. On the other hand, processes involving WDM radiation limit will be, in general, different that those involving the standard radiation (photons and neutrinos). The case where WDM dominates even in early times deserves a more carefully study of earlier processes like nucleosynthesis, reionization, reheating, etc and is beyond the scope of this paper.. This requirement leads to an upper bound on the warmness $b$-parameter,

[TABLE]

Note that in early times, radiation dominates over baryons which decay as $a^{-3}$. For this reason we do not take them into account in Eq. (7). Since the present-day values are $\Omega_{r0}\sim 10^{-4}$ and $\Omega_{w0}\sim 10^{-1}$, we expect that $\,b^{2}\lesssim 10^{-6}$, which corresponds to a DM particle velocity approximately equal to $300\,km/s$. Mathematically WDM dominating over ”standard”radiation means the absence of a real value for $z_{eq}$, that is the redshift at the point of radiation and matter equilibrium. One can evaluate $z_{eq}$ from the relation

[TABLE]

with the following solution,

[TABLE]

The early domination of WDM is shown in Fig. 1. The left panel of the Fig. 1 shows the densities of radiation and WDM for different values of parameter $b$. We can see that for $b^{2}$-values higher than $\sim 10^{-6}$ there is no radiation-dominated era. After inflation the universe is always dominated by WDM. Moreover, for any value of $b^{2}$-parameter smaller than $\sim 10^{-6}$, the equality between WDM and radiation happens earlier compared to the CDM case. In the right panel of the Fig. 1 one can see the plot for the scale factor dependence of fractional abundances (i.e. $\Omega_{i}(a)/\Omega_{T}(a)$) for radiation, baryons and WDM (here $\Omega_{T}(a)$ is the total density parameter). The plot in the top panel corresponds to the case $b^{2}=10^{-5}$ and clearly shows that WDM always dominates, while in the bottom panel, for $b^{2}=10^{-6}$, we still have an epoch dominated by standard radiation. In both cases the baryons contribution is subdominant. Consider now the structure formation process, which is strongly dependent on the behavior of WDM both at the background and perturbative level. The dynamics of WDM perturbations has been described in Fabris:2008qg , so we can just write down the main result for the dynamics of WDM perturbations. The energy and momentum balance equation, in Fourier space for each $k$-mode in flat universe lead to following equations:

[TABLE]

For the sake of convenience we used synchronous gauge and hence $h$ is the trace of the scalar metric perturbations, $\delta_{w}\equiv\delta\rho_{w}/\rho_{w}$ is the WDM density contrast and $\theta_{w}$ is the velocity. We have followed conventions for metric signature and Fourier transform of Ma:1995ey , and the dot represents derivative with respect to conformal time $\eta$ and ${\cal{H}}=\dot{a}/a$. Note that for $w=0$, CDM case is reproduced in the above equations. The equations are written in the frame which is co-moving to the WDM fluid and hence here was considered the rest-frame sound speed $c^{2}_{s}$ Weller:2003hw ; Gordon:2004ez . We shall consider WDM as adiabatic fluid. Actually, as far as we are dealing with thermal systems, it is possible that some intrinsic non-adiabaticity traces could be present. However, as a first approximation, we suppose that they are negligible. Then one can use the relation

[TABLE]

where $c^{2}_{s}=\dot{p}_{w}/\dot{\rho}_{w}$. Eqs. (10) and (11) require analytical expression for the rest-frame sound speed and the derivative of the state parameter with respect to the conformal time. Using the background quantities it is straightforward to obtain,

[TABLE]

In order to solve the system it is necessary to fix initial conditions. The WDM initial conditions can be implemented in the super-horizon regime and deep into the radiation-dominated epoch, i.e., $a\propto\eta$. In the fluid description, for the early radiation era, WDM case can be described by the following equations:

[TABLE]

By solving equation for $h$ in the super-horizon limit and in the radiation era we arrive to the well-known solution $h\propto(k\eta)^{2}$. With the last solution we found, for the relevant limits, that $\delta_{w}=-\frac{2}{3}C(k\eta)^{2}$ and $\theta_{w}=-\frac{1}{18}Ck(k\eta)^{3}$ are appropriate initial conditions. Of course, equations (10) and (11) are coupled with DE via background solutions and with baryons and radiation both at the background and perturbative level. One has to solve the complete system in order to analyse the consequences of DM warmness for the observables such as, e.g., CMB power spectrum, linear matter power spectrum and the transfer function.

III Consequences of DM warmness via RRG

In addition to Eqs. (10) and (11) we need also the equations describing perturbations for baryons and radiation. These equations can be found for example in Ma:1995ey and we will not repeat them here. To integrate the system including baryons, radiation, WDM and cosmological constant, we modify the Boltzmann CAMB code Lewis:1999bs . The initial value $\Omega_{b0}=0.0223h^{-2}$ is chosen to provide the agreement with Big Bang nucleosynthesis Pettini:2012ph , while $\Omega_{r0}$ is taken to agree with CMB measurements Ade:2015xua . The free parameters related to WDM are $H_{0}$, $\Omega_{w0}$ and $b$, and in principle they have as priors $0<H_{0}<100$, $0<\Omega_{w0}<1$ and $0<b^{2}\lesssim 10^{-6}$ to ensure a radiation dominated era. In this way we consider a reduction of this WDM space of parameters by using background observational tests. Thus, in what follows we limit our analysis to the values of WDM parameters such that they are in $1\sigma$ CL (confidence level) region of the joint analysis based on SNIa, BAO and $H(z)$ data. This shall help us to get more realistic and measurable warmness effects that do not contradict observations, at least at the background level.

III.1 Background tests

The background tests related to SNIa, BAO and $H_{0}$ are based on the likelihood computed using the $\chi^{2}$ function,

[TABLE]

where $\theta=(H_{0},\Omega_{m0},b)$ and $\Delta y(\theta)=y_{i}-y(x_{i};\theta)$. Here $y(x_{i};\theta)$ represents the theoretical predictions for a given set of parameters, $y_{i}$ the data and $\mathbf{C}$ is the covariance matrix. Note that, for convenience, the total matter density parameter $\Omega_{m0}=\Omega_{w0}+\Omega_{b0}$ was used here as a free parameter instead of $\Omega_{w0}$.

In order to perform the background statistical analysis it was used the numerical code CLASS Blas:2011rf combined with the statistical code MontePython Audren:2012wb . For the data set we have used the complete SNeIa data and correlation matrix from the JLA sample Betoule:2014frx , $H_{0}$ is considered from Riess:2011yx , and for BAO test we have used data from 6dFGSBeutler:2011hx , SDSS Anderson:2013zyy , BOSS CMASS Padmanabhan:2012hf and WiggleZ survey Kazin:2014qga . The 6dFGS, SDSS and BOSS CMASS data are mutually uncorrelated and also they are not correlated with WiggleZ data, however we must take into account correlation beetween WiggleZ data points given in Kazin:2014qga . The set of free parameters $\theta$ can be divided in two parts: the cosmological free parameters $\Omega_{m0}$, $h$ and $b$; and the nuissance parameters $\alpha$, $\beta$, $M$ and $\Delta_{M}$, related to SNe Ia data. The results of the complete statistical analysis is presented in TABLE 1 and the contour curves are presented in FIG. 2. These results are in agreement with the previous results Fabris:2008qg ; Fabris:2011am but here we have updated the results and error was reduced due to the improved quality of observational data in recent years.

III.2 Perturbative analysis

The reduced space of parameters is in agreement with the SNIa+BAO+$H(z)$ data, which we will use to study consequences of the DM warmness in the two relevant observables, namely the structure formation and CMB anisotopies. Before starting the coprresponding consideration, let us illustrate the consequences of the free-streaming of WDM in the matter perturbations.

Concerning the structure formation, a relevan quantity is the total matter density contrast,

[TABLE]

By recalling that for each component $\delta\rho_{x}=\rho_{x}\,\delta_{x}$, one can write analytic expression for the total matter density contrast in the RRG-based model,

[TABLE]

Let us note that in this case, different from $\Lambda$CDM, the total matter density contrast depends on the scale factor. Furthermore, after decoupling, the contribution of warm matter to total matter density varies from $\sim 100\%$ for $a\ll 1$ ($\delta_{m}\approx\delta_{w}$) to $\sim 87\%$ when $a=1$ (i.e $\delta_{m}\approx 0.87\delta_{w}+0.13\delta_{b}$) while in $\Lambda$CDM the contribution is always constant and of the order $\sim 85\%$ (i.e $\delta_{m}\approx 0.85\delta_{w}+0.15\delta_{b}$).

The left panel of the FIG. 3 shows the total matter density contrast for different scales and for $b^{2}=10^{-14}$. In the top panel it is shown $\delta_{m}$ for scale $k=2hMpc^{-1}$ and in the bottom panel it is shown $\delta_{m}$ for scale $k=5hMpc^{-1}$. In the first case the difference with CDM case is minimal and $\sim 5\%$ at maximum. However, in the second case, this difference goes to $\sim 20\%$. These results indicate a strong suppression of the growth of matter perturbations at the small scales, in contrast with the CDM case. Once again, one can see that the RRG enables one to reproduce known features of WDM in a very economic way.

One should expect that the suppression in the total matter density contrast caused by DM warmness also appears in the linear matter power spectrum and in its transfer function. The linear matter power spectrum is computed as $P(k)\propto k^{n_{s}}T(k)^{2}$, where $n_{s}$ is the scalar spectral index and $T(k)$ is the transfer function. The transfer function is defined as

[TABLE]

At the next stage we use our modified CAMB code to compute the linear matter power spectrum and transfer function. The results for $b^{2}=10^{-13},10^{-14}$ and $10^{-15}$ are shown in the right panel of the FIG. 3. The top right panel shows the linear matter power spectrum while bottom right panel shows the transfer function for these cases. From these plots one can conclude that at large scales there is no much deviation from the $\Lambda$CDM result, while at the small scales there is considerably lack of power proportional to the value of warmness $b$ in relation to the $\Lambda$CDM. This situation by itself is not new at all, it is regarded as one of the main features of WDM models. However, it is remarkable that one can reproduce it by using the simple RRG description, in a model-independent way and without any supposition about particle physics models.

One can wonder how CMB power spectrum is affected by the suppression on matter overdensities in small scales. The FIG. 4 shows the CMB temperature power spectrum for different $b$-values . One can observe that even with the strong suppression in $P(k)$, the CMB temperature power spectrum is not considerably affected for $b^{2}\lesssim 10^{-10}$. At large scales, when $l\lesssim 30$, all curves coincide. The differences only appear at the scales smaller than $l\sim 30$. The most of the difference is at the intermediate scales $30\lesssim l\lesssim 1300$. In order to quantify deviations from $\Lambda$CDM we compute difference

[TABLE]

where $C_{l}$ represents the CMB temperature power spectrum. Bottom panel of Fig. 4 shows $\Delta D_{l}$. Notice that the maximum difference is $\lesssim 0.015\%$ and takes place for $b^{2}=10^{-10}$. However, $\Delta D_{l}$ could be slightly higher for the values of $b$ larger than $b^{2}=10^{-10}$. Even though large scales $l\lesssim 30$ are not influenced by thermal velocities of dark matter, the rest of the spectrum does. It is possible to show that some velocities $\gtrsim 30km/s$ ($b^{2}\gtrsim 10^{-8}$) produce strong distortions in the interval $30\lesssim l\lesssim 1300$, such that it would hardly fit the data.

It is interesting to compare the RRG-based results with the ones which are based on different approaches. In the context of WDM, the effect of the free-streaming on matter distribution is quantified by a relative function transfer $\bar{T}(k)$ which is defined as

[TABLE]

where $\,P_{\Lambda\rm{WDM}}\,$ and $\,P_{\Lambda\rm{CDM}}\,$ are linear matter power spectra for $\Lambda\rm{WDM}$ and $\Lambda\rm{CDM}$ cases, respectively. The function $\bar{T}(k)$ can be approximated by the following fitting expression Bode:2000gq ,

[TABLE]

where $\alpha$ and $\nu$ are fitting parameters.

For the sake of comparison, let us denote the relative function transfer computed via RRG by $\bar{T}_{RRG}(k)$. After computing $\bar{T}_{RRG}(k)$, we perform a fit for Eq. (21) and find parameters $\alpha$ and $\nu$ for different values of $b$. The results are summarized in three first entries of Table II. The results show that RRG reproduces the relative transfer function which is considered standard in the WDM framework with accuracy of $\lesssim 1\%$, which can be seen in FIG. 5.

In the left top panel of FIG. 5 we show $\bar{T}_{RRG}(k)$ with $b^{2}=10^{-15}$ and $\bar{T}(k)$ with $\alpha=0.0147$ and $\nu=1.12$, and in its bottom panel it is shown the relative error between $\bar{T}_{RRG}(k)$ and $\bar{T}(k)$. The same plots are shown in right panel of FIG. 5 for the case where the $\bar{T}_{RRG}(k)$ was computed with $b^{2}=10^{-14}$ and $\bar{T}(k)$ was computed with $\alpha=0.0242$ and $\nu=1.12$. Note that, for both cases, the relative error is $\lesssim 1\%$.

In what follows we shall consider a more detailed comparison between RRG approach and the well established particle physics candidate for WDM associated to thermal relics. Our comparison shall include some non-linear features.

IV Thermal relics via RRG

In the context of thermal relics it is well-known that there is a lower bound999This limit comes from high redshift Lyman-alpha forest data Viel:2013apy . of $3.3\,keV$ for the WDM particle mass. Hence it would be interesting to perform a comparison between thermal relics with such a bound and RRG approach. For this reason, it is necessary first to find a equivalence between thermal relics mass scales and RRG $b$-parameter. Thus we recall that for relics we have $\nu=1.12$ and the parameter $\alpha$, in units of $h^{-1}\rm{Mpc}$, is related to the mass scale $m_{w}$ via Viel:2005qj ; Bode:2000gq ; Hansen:2001zv ,

[TABLE]

In order to obtain a complete association of the RRG parameter $b^{2}$ and the mass of WDM particles in the thermal relics context, it was used the following particular set of values for the $b^{2}$,

[TABLE]

For each point it was defined a $\chi^{2}$ function,

[TABLE]

where the $\bar{T}\left(k\right)^{th.}$ is given by the equation (21) and the $\bar{T}\left(k\right)^{num.}$ is obtained with the CAMB code for each value of $b^{2}$ in the set (23). Then we minimize the equation (24) in order to find the best-fit $\alpha$-value for each $b^{2}$.

Using the equation (22) we found, for each value of $b^{2}$, the corresponding mass $m_{w}$. This correspondence is shown in FIG. 6, where the dots indicate the best-fit value for $m_{w}$ found through the equation (24) and the solid line is the linear regression in the loglog frame. This linear regression results in the following fit-formula,

[TABLE]

By using eq.(25) we found that a mass scale of $3.3\,keV$ in thermal relics is equivalent to $b^{2}=2.36\times 10^{-15}\sim 10^{-15}$ in RRG approach. This value for $b^{2}$ brings difficulties in distinguishing between WDM and CDM scenarios both at background and linear level. At background this can be seen in the left panel of the FIG. 1, where for $b^{2}\lesssim 10^{-6}$ the expansion dynamics after matter-radiation equality is indistinguishable from CDM case and also, the best fit value for $\Omega_{m0}$ is almost the same as $\Lambda$CDM (see FIG. 2). On the other hand, at linear regime, CMB signal via RRG for $b^{2}\lesssim 10^{-10}$ is completly indistinguishable of $\Lambda$CDM (see FIG. 4) and linear matter power spectrum has the expected little suppreession in small scales (see FIG. 3). In order to better observe such tiny differences, we can for example, recompute the matter power spectrum in the Fourier space in its dimensionless form, denoted by $\Delta^{2}\left(k\right)=\frac{k^{3}}{2\pi^{2}}P(k)$ .

The left panel of the FIG. 7 shows $\Delta^{2}(k)$ at $z=0$ for both cases: standard treatment for thermal relics with mass $3.3\,keV$ and thermal relics via RRG approach with $b^{2}=10^{-15}$. We can see that CDM and WDM (thermal relics) are indistinguishable until $\Delta^{2}\left(k\right)\sim 30$ where $\Delta^{2}\left(k\right)$ falls off too rapidly. In fact, in the context of thermal-like candidates, this feature has been used as justification to claim that albeit such a cut-off is still consistent with constraints based on Lyman-$\alpha$ forest data, for masses larger than 3.3 $keV$, the WDM appears no better than CDM in solving the small scale anomalies Schneider:2013wwa . Of course, this affirmation needs to be better investigated in the context of others WDM candidates and, if possible, in a model-independent way. Thus, we believe that RRG can be very helpful in such direction.

On the other hand, we can also compute the time scale where the perturbations reach the non-linear regime. In the right panel of the FIG. 7 it is shown the time (redshift) scale $z_{nl}$ in which the matter perturbation scale $R$ reachs the non-linear regime. $z_{nl}$ is the redshift where $\sigma^{2}(R)=1$, being $\sigma^{2}(R)$ the mass variance at the scale $R$. Again both cases are considered: standard treatment for thermal relics with mass $3.3\,keV$ and thermal relics via RRG approach with $b^{2}=10^{-15}$. We can see that, in WDM thermal-like candidates context, the non-linearity is reached more recently than in $\Lambda$CDM for scales $R\lesssim 1\,Mpc\,\,h^{-1}$.

Finally, by using $\Delta^{2}\left(k\right)$, $z_{nl}$ and a top-hat filter normalized with the results from Planck Ade:2015xua it is possible to obtain the density of halos. This quantity gives us the number of collapsed objects above a given mass $M$. Here, the mass function was computed using the method of Sheth:2001dp . In FIG. 8, the mass function obtained following the RRG approach is compared with the one computed for thermal relics in the standard way and with CDM case. Clearly in WDM context, there is less collapsed objetcs than in $\Lambda$CDM. Furthermore, from FIGS. 7 and 8 it is possible to see the huge similarity between the RRG approach and the WDM thermal-like standard description also in the non-linear regime. In fact, the relative error in both cases is $\lesssim 1\%$.

Our results in this and in the previous section indicate that RRG approach is good enough to capture important features of WDM in linear regime. Specially the suppression on small scales structures and lack of power in matter spectrum at such scales. Also, in the particular case of thermal relics, RRG reproduces with high precision, the potentialities and weakness of the candidate in both linear and non linear regime. Thus, RRG approach could be considered as a complementary alternative approach to investigate warm matter and specially for understanding general behavior of the WDM scenario in a model-independent way.

V Discussion and conclusions

We have shown that the RRG approach enable one to model WDM, which is treated as a gas of particles with non-negligible thermal velocities. The simplifying aspect of RRG is that all such velocities are taken to be equal. The presence of warmness produce consequences on the dynamics of the universe both for the background and perturbations.

For the background the most important is that radiation era may be smearing out for greater warmness. In this case the universe is dominated by WDM up to the DE dominated era. This scenario could have deep impact on the primordial nucleosynthesis, reionization, recombination and other effects in the early universe. This new non-standard scenario may be deserving more careful and specific investigation, which we leave for the future work. Our analysis here was limited by the relatively small warmness, with $b^{2}\lesssim 10^{-6}$. In this case the radiation dominated era is still maintained, but the radiation-matter equality takes place before than in the CDM models. After the equality point, there are no serious differences with expansion is dominated by cold matter, as it is shown at the left panel of FIG. 1.

Instead of dealing with the full space of the WDM parameters, we restrict consideration to a more reduced set. At the background level this is achieved by using recent data from SNIa, $H_{0}$ and BAO. As a result we arrive at reduced space which does not contradict more recent background observations at $1\sigma$ CL, i.e $b\in[0,2.1\times 10^{-5}]$ and $\Omega_{m0}=[0.27,0.35]$. Taking this into account one can expect that the quantification of imprint of warmness on observables should become more significant. Our analysis in this reduced space of parameters shows that velocities which satisfy $\nu\lesssim 3\,km/s$ would agree with the CMB observations. This limit is essentially smaller than the bound for HDM, which may have velocities which are just two order of magnitude smaller than the speed of light. The velocities bound which were found here agree with the ones found earlier by other approaches Bond:1983hb ; Hannestad:2000gt ; Viel:2005qj ; Bode:2000gq ; Barkana:2001gr . This fact shows that, regardless of that the RRG is technically simple, it is a sufficiently reliable approach to probe new physics within the WDM approach.

The consequences of a warmness of DM are is more evident in the dynamics of cosmic perturbations. Since DM is supposed to be the main source of forming gravity potentials and overdensities, the impact of the DM warmness on structure formation and CMB anisotropies is evident. In the case of WDM thermal velocities cause free-streaming out from overdense regions, delaying and inhibiting the growth of fluctuations at certain scales. Another way to interpret this effect is by relating velocity to pressure. The non-negligible pressure of WDM, together with the radiation pressure, are resisting the gravitational compression and therefore suppress the power. This effect is stronger in small scales, as can be seen in Fig. 3. Furthermore, from Fig. 4 one can conclude that thermal velocities do not affect considerably the CMB temperature power spectrum for $b\lesssim 10^{-10}$. Indeed, the situation can be opposite for sufficiently large values of $b$.

As a next step we reproduced features of thermal relics by using RRG prescription. First it was necessary to find a relation between the mass scale in relics context and the $b^{2}$ parameter of RRG. The $b^{2}\sim 10^{-15}$ equivalent to the lower bound known for thermal relics ($3.3keV$) brings the necessity to look more carefully matter perturbations. Thus, we computed the matter perturbations in the Fourier space $\Delta^{2}(k)$, the time scale where the perturbations reach the non-linear regime ($z_{nl}$) and the mass function. In all those cases, our results indicate that RRG approach is good enough to capture important features of WDM even in non-linear regime. Therefore, we have proved that RRG is a reliable model, and can be considered as a complementary, greatly simplified alternative approach to investigate warm matter, in particular for understanding the behavior of the WDM in a totally model-independent way.

One can foresee the possibility of detailed investigations of the new Particle Physics candidates to WDM, which includes the comparison with model-independent RRG. In our opinion some aspects of this possibility would be quite interesting. For example, let us mention a relation between $\bar{T}_{RRG}(k)$ and WDM candidate models, a more complete and comprehensive exploration of the space of parameters via Markov Chain Monte Carlo (MCMC) or verifying how RRG would work in the nonlinear regime of structure formation through numerical simulations and the possibility if WDM is really capable of solve small scale problems. The work on these aspects of the model is currently in progress.

Acknowledgments

Authors are grateful to Winfried Zimdahl for useful discussions. J.F. and R.M. wish to thank CAPES, CNPq and FAPES for partial financial support. I.Sh. acknowledges the partial support from CNPq, FAPEMIG and ICTP.