On general features of warm dark matter with reduced relativistic gas

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.

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ü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 modelindependent 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 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 1%. A warmness with b 2 10 −6 (equivalent to v 300km/s) does not alter significantly the CMB power spectrum and is in agreement with the background observational tests.

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 ∼ 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 [1,2], baryon acoustic oscillations [3], anisotropies of the cosmic microwave background (CMB) [4,5] and other observations [6]. 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, ΛCDM model, assumes that the DE is a cosmological constant, and regards DM as a non-relativistic matter with negligible pressure (cold dark matter). ΛCDM provides an excellent agreement with most of the data (see, e.g., [7] for a general review), however this agreement is not perfect due to the tensions with part of the observational data (see for example [8]). In part due to these difficulties, some alternative models have been proposed and studied as possible DE and DM candidates (see for example [9][10][11][12][13]). Let us note that some of these alternative models aim to describe fluids that replace both DM and DE (see for example [14][15][16][17][18][19]) or describe interaction between DE and DM [20][21][22][23][24][25][26].
Some of the mentioned ΛCDM difficulties are related with the choice of the cold dark matter (CDM) paradigm [27]. For instance, at small scales the issues such as the missing satellites problem [28], core/cusp problem [29], and the Too big to fail problem [30], 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 [31][32][33][34][35][36] 1 . 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 [7]. From this perspective, it is useful to develop also modelindependent 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 [37]. Using this model Sakharov predicted the existence of oscillations in CMB temperature spectra long before its observational discovery (see [38] for the historical review).
Recently, RRG was reinvented in [39], 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 [40,41]. 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 [39] 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. [42] 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. [43].
The RRG was successfully used in [39,42] as an interpolation between radiation and dark matter eras in cosmology. An upper bound on the warmness coming from RRG [42,44] 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 [31][32][33][34][35][36], 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 [42], 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 2 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 Λ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 [39,42] or recent [45] 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 β = v/c [39] (c is the light speed). This leads to the following relation between WDM pressure p w and WDM energy density ρ w , where m is the WDM particle mass, and is the kinetic energy of each particle of the system which is given by, Here ρ c is introduced as a notation for the rest energy density, i.e., the energy density for the v = 0 case. Thus, ρ c = ρ 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 which can be regarded as equation of state (EoS) of the WDM fluid. Using Eq. (3) in the energy conservation relation, the solution for ρ w has the form Thus, b parameter measures velocity and warmness of the WDM particles at present. In the limit v c we have b ≈ 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, 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 In the last equation Ω x0 (with x = Λ, w, b and r) is the value of the DE, WDM, baryons and radiation density parameters at present, while Ω Λ0 = 1 − Ω w0 − Ω b0 − Ω 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 → 0 and radiation at b → ∞ cases. Because of this interpolation feature, RRG can be used to investigate the cosmological consequences of the transition between epochs of radiation and dust [37][38][39]. 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 density 3 . This requirement leads to an upper bound on the warmness b-parameter, 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 Ω r0 ∼ 10 −4 and Ω w0 ∼ 10 −1 , we expect that b 2 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 with the following solution, 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 ∼ 10 −6 there is no radiationdominated era. After inflation the universe is always dominated by WDM. Moreover, for any value of b 2 -parameter smaller than ∼ 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. Ω i (a)/Ω T (a)) for radiation, baryons and WDM (here Ω 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 [42], 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: For the sake of convenience we used synchronous gauge and hence h is the trace of the scalar metric perturbations, δ w ≡ δρ w /ρ w is the WDM density contrast and θ w is the velocity. We have followed conventions for metric signature and Fourier transform of [36], and the dot represents derivative with respect to conformal time η and H =ȧ/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 [46,47]. 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 where c 2 s =ṗ w /ρ 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,ẇ and 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 ∝ η. In the fluid description, for the early radiation era, WDM case can be described by the following equations: By solving equation for h in the super-horizon limit and in the radiation era we arrive to the well-known solution h ∝ (kη) 2 . With the last solution we found, for the relevant limits, that δ w = − 2 3 C(kη) 2 and θ w = − 1 18 Ck(kη) 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 [36] and we will not repeat them here. To integrate the system including baryons, radiation, WDM and cosmological constant, we modify the Boltzmann CAMB code [48]. The initial value Ω b0 = 0.0223h −2 is chosen to provide the agreement with Big Bang nucleosynthesis [49], while Ω r0 is taken to agree with CMB measurements [5]. The free parameters related to WDM are H 0 , Ω w0 and b, and in principle they have as priors 0 < H 0 < 100, 0 < Ω w0 < 1 and 0 < b 2 10 −6 to ensure a radiation dominated era. In this way we consider

A. Background tests
The background tests related to SNIa, BAO and H 0 are based on the likelihood computed using the χ 2 function, where θ = (H 0 , Ω m0 , b) and ∆y(θ) = y i − y(x i ; θ). Here y(x i ; θ) represents the theoretical predictions for a given set of parameters, y i the data and C is the covariance matrix. Note that, for convenience, the total matter density parameter Ω m0 = Ω w0 + Ω b0 was used here as a free parameter instead of Ω w0 . In order to perform the background statistical analysis it was used the numerical code CLASS [50] combined with the statistical code MontePython [51]. For the data set we have used the complete SNeIa data and correlation matrix from the JLA sample [52], H 0 is considered from [53], and for BAO test we have used data from 6dFGS [54], SDSS [55], BOSS CMASS [56] and WiggleZ survey [57]. 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 [57]. The set of free parameters θ can be divided in two parts: the cosmological free parameters Ω m0 , h and b; and the nuissance parameters α, β, M and ∆ M , related to SNe Ia data. The results of the complete statistical analysis is presented in TABLE I and the contour curves are presented in FIG. 2. These results are in agreement with the previous results [42,44] but here we have updated the results and error was reduced due to the improved quality of observational data in recent years.

B. 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, By recalling that for each component δρ x = ρ x δ x , one can write analytic expression for the total matter density contrast in the RRG-based model,  Let us note that in this case, different from Λ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 ∼ 100% for a 1 (δ m ≈ δ w ) to ∼ 87% when a = 1 (i.e δ m ≈ 0.87δ w + 0.13δ b ) while in ΛCDM the contribution is always constant and of the order ∼ 85% (i.e δ m ≈ 0.85δ w + 0.15δ 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 δ m for scale k = 2hM pc −1 and in the bottom panel it is shown δ m for scale k = 5hM pc −1 . In the first case the difference with CDM case is minimal and ∼ 5% at maximum. However, in the second case, this difference goes to ∼ 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) ∝ k ns T (k) 2 , where n s is the scalar spectral index and T (k) is the transfer function. The transfer function is defined as Note that suppression in small scales is proportional to b and is more evident in these quantities.
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 ΛCDM result, while at the small scales there is considerably lack of power proportional to the value of warmness b in relation to the Λ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 10 −10 . At large scales, when l 30, all curves coincide. The differences only appear at the scales smaller than l ∼ 30. The most of the difference is at the intermediate scales 30 l 1300. In order to quantify deviations from ΛCDM we compute difference where C l represents the CMB temperature power spectrum. Bottom panel of Fig. 4 shows ∆D l . Notice that the maximum difference is 0.015% and takes place for b 2 = 10 −10 . However, ∆D l could be slightly higher for the values of b larger than b 2 = 10 −10 . Even though large scales l 30 are not influenced by thermal velocities of dark matter, the rest of the spectrum does. It is possible to show that some velocities 30km/s (b 2 10 −8 ) produce strong distortions in the interval 30 l 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 T (k) which is defined asT where P ΛWDM and P ΛCDM are linear matter power spectra for ΛWDM and ΛCDM cases, respectively. The function T (k) can be approximated by the following fitting expression [34],  (19). Note that for l 30 all curves are indistinguishable and differences appear after l ∼ 30. For velocities of the order of v 3km/s, differences with CDM case is ∼ 0, 015% at maximum.
where α and ν are fitting parameters. For the sake of comparison, let us denote the relative function transfer computed via RRG byT RRG (k). After computingT RRG (k), we perform a fit for Eq. 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 bound 4 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 ν = 1.12 and the parameter α, in units of h −1 Mpc, is related to the mass scale m w via [33,34,59], 0.01H 0 0.7 1.22 .
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 , b 2 = 1 × 10 −15 , 2 × 10 −15 , . . . , 9 × 10 −12 , 1 × 10 −11 .
For each point it was defined a χ 2 function, where theT (k) th. is given by the equation (21) and theT (k) 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 α-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, m w = 4.65 · 10 −6 b 2 −2/5 keV.
By using eq.(25) we found that a mass scale of 3.3 keV in thermal relics is equivalent to b 2 = 2.36 × 10 −15 ∼ 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 10 −6 the expansion dynamics after matter-radiation equality is indistinguishable from CDM case and also, the best fit value for Ω m0 is almost the same as ΛCDM (see FIG. 2). On the other hand, at linear regime, CMB signal via RRG for b 2 10 −10 is completly indistinguishable of Λ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 ∆ 2 (k) = k 3 2π 2 P (k) . The left panel of the FIG. 7 shows ∆ 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 ∆ 2 (k) ∼ 30 where ∆ 2 (k) 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-α forest data, for masses larger than 3.3 keV , the WDM appears no better than CDM in solving the small scale anomalies [60]. 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 σ 2 (R) = 1, being σ 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 ΛCDM for scales R 1 M pc h −1 .
Finally, by using ∆ 2 (k), z nl and a top-hat filter normalized with the results from Planck [5] 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 [61]. 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 Λ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 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 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σ CL, i.e b ∈ [0, 2.1 × 10 −5 ] and Ω 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 ν 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 [31][32][33][34][35]. 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 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 ∼ 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 ∆ 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 betweenT 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.