Dark Matter Sterile Neutrino from Scalar Decays

: We place constraints on DM sterile neutrino scalar decay production (SDP) assuming that sterile neutrinos representa fraction from the total Cold Dark Matter energy density. For the cosmological analysis we complement the CMB anisotropy measurements with CMB lensing gravitational potential measurements, that are sensitive to the DM distribution to high redshifts and with the cosmic shear data that constrain the gravitational potential at lower redshifts than CMB. We also use the most recent low-redshift BAO measurements that are insensitive to the non-linear effects, providing robust geometrical tests. We show that our datasets have enough sensitivity to constrain the sterile neutrino mass m ν s and the mass fraction f S inside the co-moving free-streaming horizon. We ﬁnd that the best ﬁt value m ν s = 7.88 ± 0.73 keV (68% CL) is in the parameter space of interest for DM sterile neutrino decay interpretation of the 3.5 keV X-ray line and that f S = 0.86 ± 0.07 (68% CL) is in agreement with the upper limit constraint on f S from the X-ray non-detection and Ly- α forest measurements that rejects f S = 1 at 3 σ . However, we expect that the future BAO and weak lensing surveys, such as EUCLID, will provide much more robust constraints.


Introduction
The Cosmic Microwave Background (CMB) measurements from the PLANCK satellite, alone or in combination with other astrophysical datasets, provide no powerful evidence supporting new physics beyond the standard ΛCDM cosmological model [1,2]. With around 5% of the total energy density of the universe representing the baryonic matter, 21% the Dark Matter (DM) and 74% accounting for the Dark Energy (DE), the ΛCDM model is remarkably successful in reproducing the large-scale structure (LSS) of the universe. In addition, the Planck results show that the signature of the neutrino sector is consistent with the ΛCDM model assumptions and that DE is compatible with the Λ cosmological constant.
The nature and composition of DM is still unknown. Attempts involving collision-less DM particles fail to solve the ΛCDM problems in reproducing the cosmological structures at small scales (missing satellite problem [3][4][5], core-cusp problem [6][7][8], too-big-to-fail problem [9]), suggesting that DM particles may also exhibit gravitational properties and requiring the extension of the Standard Model (SM) of particle physics [10][11][12].
One of the theoretically well-motivated DM candidate is the sterile neutrino [13][14][15]. Arising in the minimal extension of SM, the sterile neutrino with mass in keV range can simultaneously explain the active neutrino oscillations, the DM properties and the matter-antimatter asymmetry of the universe [16,17].
Few independent detections of a weak X-ray emission line at an energy of ∼3.5 keV have been found in a stacked XMM-Newton spectrum of 73 galaxy clusters with redshifts in the range 0.01-0.35, in Chandra ACIS-I and ACIS-S spectra of Perseus [18] and in the XMM-Newton spectra of the Andromeda galaxy and the Perseus galaxy cluster [19]. Evidences of this emission have been also searched in the Milky Way with data from XMM-Newton [20,21] and Chandra [22], in deep Suzaku X-ray spectra of the central regions of Perseus, Coma, Virgo and Ophiuchus clusters [23,24], in the stacked XMM-Newton spectra of dwarf spheroidal galaxies in the vicinity of the Milky Way [25] and in the stacked spectra of a sample of galaxies selected from Chandra and XMM-Newton public archives [26]. The origin of this line related either to a signature of decaying DM sterile neutrino with a mass of 7.1 keV [20] or to any known atomic transition in thermal plasma (with special interest in Potassium and Chlorine atomic transitions) [27] raised several controversies [28]. The evidence of this emission line in different astrophysical sites is expected to place strong constraints on both hypotheses.
Several keV sterile neutrino DM production mechanisms have been proposed.
In the Dodelson-Widrow (DW) scenario [29], keV sterile neutrinos are produced by non-resonant oscillations with active neutrinos in presence of negligible leptonic asymmetry. This mechanism is now excluded by the observations of structure formation as it produces too hot sterile neutrino velocity spectra [30,31].
The keV sterile neutrino DM resonant production (RP) via the conversion of active to sterile neutrinos through the Shi-Fuller mechanism in presence of leptonic asymmetry [32] has also been investigated [33]. In this scenario, sterile neutrino parameters required to reproduce the X-ray line of ∼3.5 keV are consistent with the main cosmological parameters inferred from present cosmological measurements, Local Group and high-z galaxy count constraints [33]. However, some tension with Ly-α data still exists [34].
The keV sterile neutrino DM production by particle decays has been also extensively discussed (see [35] and references therein). A particularly interesting case is the DM sterile neutrino production by scalar decay (SDP). This process involves a generic scalar singlet with the vacuum expectation value (VEV) that could be produced via SM Higgs interactions. Depending on the strength of the Higgs coupling λ H , the singlet scalar can be produced like Weakly Interacting Massive Particles (WIMPs) via freeze-out [15,36] or like "Feeble Interacting Massive Particles" (FIMPs) via freeze-in mechanisms and must couple with the right-handed neutrino fields through the Yukawa interaction, leading to sterile neutrino Majorana masses [34,37]. A complete treatment of the SDP mechanism for the whole parameter space, giving the general solution at the level of momentum distributions is presented in [38].
So far, the keV sterile neutrino DM properties have been addressed by evaluating their impact on the co-moving free streaming horizon, that relates to the average of sterile neutrino velocity distribution. However, for such models characterised by a highly nonthermal momentum distribution, the average momentum is subject to uncertainties, leading to a fail of the free-streaming horizon in constraining the sterile neutrino parameters [34]. The existing constraints are, in general, obtained in linear theory under the assumption that sterile neutrinos are all DM [37,39].
The aim of this paper is to place constraints on the properties of DM sterile neutrinos produced by the SDP mechanism through their impact on distance-redshift relations and the growth of structures. We consider models where DM is a mixture of CDM and DM sterile neutrinos and analyse if the existing measurements of the CMB gravitational potential, the baryon acoustic oscillation (BAO) and the weak gravitational lensing of galaxies can constrain DM sterile neutrinos properties.
The paper is organized as follows: Section 2 summarizes the SDP Boltzmann formalism calculations. Section 3 describes the model parameters and the methods involved in the analysis. Section 4 presents the datasets. In Section 5 we present our results and discuss the cosmological implications of the sterile neutrino SDP mechanism. The conclusions are summarised in Section 6.

Sterile Neutrino Production by the Scalar Decay (SDP)
In this section, we present the sterile neutrino SDP calculations. The evolution of momentum distributions for scalar, f S , and sterile neutrino, f ν s , are obtained by solving the coupled Boltzmann equations: where C S and C ν s are scalar and sterile neutrino collision terms encoding the effects of different processes that contribute to their production andL is the Liouville operator: In the above equation, p is the particle momentum and H is the Hubble expansion rate. Equation (1) can be written in a more convenient form by performing the following variable transformations [38]: Exploiting the correspondence between temperature T and time t and by using the conservation of the co-moving entropy, the above transformations can be written in the form (for details see Appendix A.2 in [38]): where q = p/T is the co-moving momentum and g s (T) is the effective number of entropy degrees of freedom. We choose m 0 = T 0 = m h where m h = 125 GeV is the Higgs boson mass. In terms of the new variables, the Liouville operator reads as: where denotes the derivative with respect to the temperature. The time-temperature relation is obtained as: Equation (6) assumes the conservation of the co-moving entropy density: where a is the cosmological scale factor (a 0 = 1 today) and s 0 = 2891.2 cm −3 is the present entropy density. We used the fitting formulas from [40] to compute the evolution of the effective number of the entropy degrees of freedom g s (T) and its derivative g s (T).
In this work, we take the leading processes contributing to the collision terms C S and C ν s from Equation (1) of the form: where: C S hh↔SS describes the depletion of scalars due to the annihilation into pairs of SM Higgs particles and the reverse process, C S→ν s ν s describes the decay of scalars into pairs of sterile neutrinos and C ν s S→ν s ν s describes the creation of sterile neutrinos from the decays of scalars. With these assumptions, the SDP scenario is parameterized with respect to the sterile neutrino mass m ν s , the scalar mass M S , the strength of the Higgs coupling λ H and of the Yukawa coupling y k .
Using the explicit forms of the collision terms given in [38], we simultaneously evolve the Equations (1) and (6) to obtain the scalar and sterile neutrino momentum distributions in the expanding universe. Figure 1 presents the dependence of the sterile neutrino final momentum distribution obtained for a given scalar mass M S and different values of λ H and y k . The sterile neutrino number density n ν s (r) and the physical energy density Ω ν s h 2 are then given by: where s(r) is the co-moving entropy density and ρ c /h 2 = 1.054 10 −2 MeV cm −3 is the critical density in terms of reduced Hubble constant h. Figure 2 presents the dependence of the abundances of scalar and sterile neutrino on the time parameter r = m h /T for the same models presented in Figure 1.
We test the production code over a large parameter space and find that our distributions are in agreement with similar distributions obtained in [37,38].
The SDP mechanism relies on the assumption that the Higgs boson, that is sourcing the scalar production, is in thermal equilibrium. To asses the validity of this assumption, we evaluate the total time available for the Higgs boson to decouple from the thermal plasma. Since DM production happens in the early Universe, we consider the era of radiation dominance.  [36]. We evaluate t(T) assuming that the Higgs boson decouples from the thermal plasma at a temperature similar to its mass, namely T h ∼ 125 GeV. Assuming g s (T h ) = 106.5, we find t(T h ) ∼ 1.51 × 10 −11 s. The Higgs vacuum lifetime is τ h ∼ 1.56 × 10 −22 s [41]. During the Higgs lifetime, the scalars are generated through hh → SS annihilation process. Then, the reverse process SS → hh generates Higgs particles in thermal equilibrium if the scalar lifetime τ S = 16π/(y 2 t m S ) is smaller than the time available for the Higgs boson to decouple from the thermal plasma. Through these processes, the Higgs boson population is kept constant and in thermal equilibrium. For example, for the pair (m S , y t ) = (533 GeV, 10 −7.2 ) we obtain τ S = 1.4 × 10 −11 s, while for (m S , y t )=(60 GeV, 10 −6 ) we obtain τ S = 5.5 × 10 −13 s. In both cases, the Higgs boson is in equilibrium sourcing the scalar production.
After the scalar develops a non-zero VEV, the DM sterile neutrinos can be produced through the decay S → ν s ν s , leading to sterile neutrino masses m ν s = y k < S >, where < S > is the scalar VEV. The collision term C ν s in Equation (8) is controlled by the scalar effective decay width function C Γ = M 0 (T)y t /m 2 S [37]. The production of sterile neutrinos occurs mainly at temperatures on the order of the scalar mass, T prod ∼ 10 2 GeV. The reduction in the effective number of degrees of freedom from g * (T prod ) = 106.5 to g * (0.1 MeV) = 3.36 and the entropy production that takes place as the universe cools from T prod causes the DM sterile neutrino population to be diluted and redshifted. Following [36], these processes dilute the population of sterile neutrinos by a factor ξ ≥ 33 in the density and by a factor ξ 1/3 ≥ 3.2 in the average momentum for the SM degrees of freedom, leading to DM sterile neutrinos cold enough to satisfy the bounds on small-scale structures. Such sterile neutrinos can decay into neutrinos and photons with energy E γ = m ν s /2, reaching the sensitivity of X-ray observatories.

Parameterization and Methods
The baseline model is an extension of the flat ΛCDM model, described by the following cosmological parameters: where: Ω b h 2 is the present baryon energy density, Ω c h 2 is the present CDM energy density, θ s is the ratio of sound horizon to angular diameter distance at decoupling, τ is the optical depth at reionization, A s and n s are amplitude and spectral index of primordial power spectrum of curvature perturbations at pivot scale k = 0.05 Mpc −1 , ∑ m ν is the total mass of three active neutrino flavors and N e f f the number of relativistic degrees of freedom that parameterize the contributions from non-interacting relativistic particles. In the SM with three active neutrino flavours, N e f f = 3.046 due to non-instantaneous decoupling corrections [42].
The SDP model extends the baseline model by including the following parameters: where: m ν s is the sterile neutrino mass, M S is the scalar mass, y k is the Yukawa strengths coupling and λ H the Higgs strengths coupling. The sterile neutrino mass fraction is a derived parameter, f S = Ω ν s /Ω c , where the sterile neutrino energy density Ω ν s is computed using Equation (10). The matter energy density is then given by We modify the Boltzmann code camb 1 [43] to allow the calculation of sterile neutrino SDP formalism presented in the previous section.
The parameter extraction from cosmological datasets is based on the Monte-Carlo Markov Chain (MCMC) technique. We modify the publicly available version of the package CosmoMC 2 [44] to sample from the space of cosmological and sterile neutrino SDP mechanism parameters and generate estimates of their posterior mean and confidence intervals. We first run the modified versions of CosmoMC and camb setting to zero the additional parameters of the SDP model. We find good consistency between our bounds and the existing constraints for the ΛCDM model [1]. We use the default convergence settings implemented in CosmoMC: MPI_Limit_Converge = 0.025 and MPI_Limit_Converge_Err = 0.2. With these choices, the CosmoMC run stops when the confidence interval for each parameter at 95% C.L. is accurate at 0.2 σ. This error can be reduced, but in this case the computing time increases to reach the convergence limit. This is critical for non-standard models as SDP for which the execution time is larger than in the standard case because of numerical evolution of momentum distributions in camb. The final runs are based on 120 independent channels, reaching the convergence criterion (R − 1) = 0.007, defined as the ratio between the variance of the means and the mean of variances for the second half of chains [44].
We assume a flat Universe and uniform priors for all parameters adopted in the analysis in the intervals listed in Tables 1 and 2. The Hubble expansion rate H 0 and sterile neutrino energy density Ω ν s h 2 are derived parameters in our analysis. We constrained H 0 values to reject the extreme models and restrict the values of Ω ν s h 2 to the interval of Ω c h 2 .
The sterile neutrino lower mass limit is restricted by the universal Tremain-Gunn bound [45] valid for the case in which all DM is made by sterile neutrinos.
A stronger bound is coming from the analysis of DM phase-space distribution in dwarf spheroidal galaxies, leading to m ν s > 1.8 keV for the sterile neutrino non-resonant production (NRP) [46]. This bound was revisited in [47] for the Λ Cold plus Warm Dark Matter model (ΛCWDM) where WDM in the form of sterile neutrinos is a fraction f WDM from the total DM. For f WDM = 1 (pure ΛWDM model), the combined analysis of WMAP5 and Ly−α datasets led to a lower bound on sterile neutrino mass m ν s > 1.6 keV (at 95%CL) in the NRP case. The same analysis shows that after a certain value of the lower bound of sterile neutrino mass f WDM approaches a constant value because the contribution of the WDM component becomes too small to be constrained by the data. We made MCMC convergence tests with different priors on the parameters of SDP models leading to m ν s = 2 keV as the smallest mass of sterile neutrino included in the present analysis . The upper bound value of the sterile neutrino mass is restricted by the nondetection of the emission lines from X-ray observations [48]. Table 1. Priors and constraints for the ΛCDM-ext parameters adopted in the analysis. All priors are uniform in the listed intervals. We assume a flat Universe.

Parameter
Prior

Cosmological Datasets
For the cosmological analysis, we use the following datasets: The CMB measurements: We use the CMB angular power spectra from the PLANCK experiment [1] and the PLANCK likelihood codes corresponding to different multipole ranges: Commander for 2 ≤ l ≤ 29, CamSpec 50 ≤ l ≤ 2500, LowLike for 2 ≤ l ≤ 32 for polarization data and Lensing for 40 ≤ l ≤ 400 for lensing data. As sterile neutrino production is expected to affect the redshift-distance relations and the growth of structures, we include in the analysis the PLANCK power spectrum of the reconstructed lensing potential [2]. We will refer to the combination of these measurements as the PLANCK+lens dataset.
The baryonic acoustic oscillations (BAO): BAO measurements are low-redshift probes insensitive to non-linear effects because their characteristic acoustic scale, of around 147 Mpc, is much larger than the scale of the virialized cosmological structures. Moreover, as BAO features in the matter power spectrum can be observed as a function of both angular and redshift separations [49,50], these measurements are robust geometrical tests. We include in the analysis the BAO characteristic parameter measurements from the Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ at z e f f = 0.32, CMASS at z e f f = 0.57 [50] and 6dF Galaxy Survey (6dFGS) at z e f f = 0.1 [51]. We will refer to the combination of these measurements as the BAO dataset.
Cosmic shear: Weak lensing of galaxies, or cosmic shear, probes the gravitational potential at redshifts lower than the CMB lensing. We use the COSMOMC implementation of the one-year Dark Energy Survey (DES) cosmic shear measurements described in [52], referred hereafter as DES dataset.

Sensitivity of Cosmological Data to m ν s and f S
A change in sterile neutrino mass fraction f S leads first to a change in the Hubble expansion rate H. At the CMB photons decoupling (T cmb = 0.26 eV), the changes in H lead to changes in the sound horizon distance r s (that scales as 1/H) and in the photon diffusion damping r d (that scales as √ 1/H).
The CMB anisotropy spectrum is sensitive to both r s and r d changes as projected angles over the co-moving angular distance, D A , to the last scattering: θ s = r s /D A and θ d = r d /D A . While a change in θ s shifts the position of the acoustic Doppler peaks through the CMB anisotropy spectrum, a change in θ d relative to θ s modifies its amplitude. CMB can measure the expansion rate by measuring the ratio θ d /θ s ∼ √ H. For models sharing the same value of f S , the ratio θ d /θ s is a measure of the relativistic energy density at T cmb , usually parameterized by N e f f (T cmb ) [53]. The contribution of sterile neutrinos to N e f f (T cmb ) encodes information on their mass and production mechanism parameters that lies in the momentum distribution function and can be computed by comparing the sterile neutrino kinetic energy density to the energy density of other relativistic particles in equilibrium at T cmb [37].
The BAO measurements lead to joint constraints on the co-moving angular diameter distance, D A (z), and on the Hubble parameter, H(z), at smaller redshifts than CMB. At these redshifts, for models sharing the same sterile neutrino mass fraction, both D A (z) and H(z) are degenerated. The BAO observations typically constrain the quantity r s (z d )/D V (z), r s (z d ) is the sound horizon distance at the drag epoch redshift z d and D V (z) is given by: where c is the speed of light. Figure 3 presents the redshift dependence of the variation ∆(r s (z d )/D V (z)) for models sharing the same f S obtained in the SDP scenario, showing that the BAO measurements break the degeneracy between these models, leading to constraints on sterile neutrino mass and production parameters.  The gravitational lensing of the CMB photons provides direct measurements on scales where the effects of sterile neutrino are expected to manifest. The largest scale affected is the present value of the co-moving free-streaming horizon given by [47]: where < v(T) > is the velocity dispersion of sterile neutrinos, T prod is the sterile neutrino production temperature and T 0 is the present temperature. We compute λ 0 f sh following the analytical approach from [47]: where a eq = Ω R /Ω m is the scale factor at matter-radiation equality, a nr = T 0 /m ν s is the scale factor at the time of sterile neutrino non-relativistic transition, Ω R the is the radiation energy density and Ω m is the matter energy density. The analytical approach assumes that the Universe is completely radiation-dominated until a eq and neglects the small contribution to the integral (12) from dark energy. We account for the entropy dilution from T prod until T 0 rescaling λ 0 f sh by a factor ξ 1/3 s = g s (T prod )/g s (T 0 ) [34], where g s is the effective number of relativistic entropy degrees of freedom (g s (T EWFT ) ≈ 109.5 and g s (T 0 ) ≈ 3.36). Figure 4 presents the likelihood probability distributions of the free-streaming horizon wave-number k min f s , obtained for our models. One should note that these wave-numbers, k ∼ 1 hMpc −1 , correspond to the typical size of dwarf galaxies [37], while the observations of Ly-α absorption in spectra of distant quasars are tracers of linear density fluctuations on scales 0.1 < k < 3 hMpc −1 . On the other hand, the power spectrum of the CMB-projected gravitational potential, C φφ l , is sensitive to both geometry and growth of structures. In the Limber approximation, C φφ l can be written as: where: χ is the co-moving coordinate distance, χ * is the co-moving coordinate distance to the last scattering surface, k is the wave-number, D A (χ) is the co-moving angular diameter distance and P Ψ (z, k) is the power spectrum of the gravitational potential. Figure 5 presents the deflection angle power spectra obtained for models sharing the same values of f S and m ν s . We indicate the contributions of different wave-numbers k (in Mpc −1 ) to the deflection angle power spectra. The figure shows that for the multipole range involved in this analysis, 40 ≤ l ≤ 400, the deflection angle power spectrum of the CMB lensing potential is sensitive to the DM sterile neutrino parameters. We conclude that PLANCK + lens + BAO + DES datasets have enough sensitivity to constrain the sterile neutrino mass and mass fraction inside the co-moving free-streaming horizon. Table 3 presents the mean values and the absolute errors of the main cosmological parameters from the fit of ΛCDM-ext and SDP models with the PLANCK + lens + BAO + DES datasets. We found a good agreement of the cosmological predictions for all cosmological parameters within 1−σ.

Cosmological Constraints
Constraints on several cosmological parameters are highly dependent on the number of priors and cosmological datasets. As we do not include the high-l CMB polarization data in the analysis, we find a higher value for τ than that obtained by the PLANCK Collaboration [1].
Our datasets prove scales where the effects of sterile neutrino are expected to manifest. Moreover, the damped amplitude of the CMB anisotropy at small scales, quantified by e −2τ , is accurately constrained by the CMB lensing power spectra, while the lensing potential reconstruction provides the determination of the amplitude independent on the optical depth. As the CMB lensing power spectrum constraints the matter density fluctuations along the line of sight, the present-day rms matter density power, σ 8 , is also determined. On the other hand, the CMB small-scale anisotropy directly fixes the combination σ 8 e −τ that is tightly constrained by the CMB data. We find a very good consistency between the values of σ 8 e −τ = 0.741 ± 0.021 (68% CL) obtained from our analysis and σ 8 e −τ = 0.768 ± 0.009 (68% CL) obtained from the analysis of PLANCK CMB data alone, in the case of the ΛCDM model. Table 3 for the ΛCDM model are in good agreement with the PLANCK Collaboration results obtained from the joint analysis of the CMB-lensing, BAO and galaxy-lensing (see Table 2 in [2]). Figure 6 shows the likelihood probability distributions and the joint confidence regions obtained for the SDP mechanism parameters.

The values of other cosmological parameters presented in
The dominant effect on DM sterile neutrino production is given by the strength of the Higgs coupling λ H , that sets the scalar mass at M S 533 ± 47 GeV and the strength of Yukawa coupling y k sets the sterile neutrino mass at m ν s = 7.88 ± 0.73 keV, leading to sterile neutrino DM mass fraction f S = 0.86 ± 0.07 (errors at 68% C.L.). This value clearly indicates that a dominant component of DM can be produced by the SDP mechanism.

Conclusions
We place constraints on the DM sterile neutrino produced by the scalar decay (SDP), assuming that sterile neutrino represents a fraction f S from the total CDM energy density. So far, the keV sterile neutrino properties have been addressed under the assumption that sterile neutrinos are all DM, by evaluating their impact on the co-moving free streaming horizon that relates on the average velocity distribution.
For such models, characterized by highly non-thermal momentum distributions, the average momentum is subject to uncertainties and the free-streaming horizon fails to constrain the DM sterile neutrino parameters.
For our cosmological analysis, we complement the CMB anisotropy data with the CMB lensing gravitational potential measurements, that are sensitive to the DM distribution to high redshifts and with the cosmic shear data, that constrain the gravitational potential at lower redshifts than the CMB anisotropy. We also use the most recent low-redshift BAO measurements that are insensitive to the non-linear effects, providing robust geometrical tests.
We show that for models sharing the same DM sterile neutrino mass fraction f S , the accurate determination of the acoustic scale from CMB anisotropy measurements breaks the degeneracy of Hubble parameter at the photon decoupling, constraining the sterile neutrino mass m ν s , while the BAO measurements constrain m ν s at lower redshifts. We evaluate the co-moving free-streaming horizon showing that the deflection angle power spectrum of the CMB lensing potential, C φφ l , is sensitive to sterile neutrino production mechanism parameters for the multipole range involved in this analysis (40 ≤ l ≤ 400). Depending on both the angular diameter distance and the matter density fluctuations, we show that C φφ l can break the degeneracy between sterile neutrino mass m ν s and mass fraction f S . We also show that our datasets have enough sensitivity to constrain m ν s and f S inside the co-moving free-streaming horizon.
We find that the best fit values of m ν s are in the parameter space of interest for sterile neutrino DM decay interpretation of the 3.5 keV X-ray line [19] and theDM sterile neutrino mass fraction is f S = 0.86 ± 0.07 (at 68% CL), in agreement with the upper limit constraint on f S from the X-ray non-detection and Ly-α forest measurements that rejects f S = 1 at 3 σ level [48]. We expect that the future BAO and weak lensing surveys, such as EUCLID, will provide much more robust constraints on the DM sterile neutrino properties.