It has been recently suggested \cite{Bezrukov:2017ike,Bezrukov:2018wvd} that a cosmic scalar field can completely change the keV-scale sterile neutrino production in the early Universe. Its effect may, for various parameter choices, either suppress sterile neutrino production and make moderate active-sterile mixing cosmologically acceptable, or increase the production and generate considerable dark matter component out of sterile neutrino with otherwise negligible mixing with SM. In this paper we provide analytic estimates complementing and providing details of the numerical calculations performed in \cite{Bezrukov:2018wvd} in the case of resonant amplification of the sterile neutrino production. We also discuss phenomenological and theoretical issues related to the successful implementation of this idea in fully realistic extensions of the Standard Model of particle physics.

We study the implications of the Planck temperature power spectrum at low multipoles, $\ell<1000$, and SPTPol data. We show that this combination predicts consistent lensing-induced smoothing of acoustic peaks within $\Lambda$CDM cosmology and yields the robust predictions of the cosmological parameters. Combining only the Planck large-scale temperature data and the SPTPol polarization and lensing measurements within $\Lambda$CDM model we found substantially lower values of linear matter density perturbation $\sigma_8$ which bring the late-time parameter $S_8=\sigma_8\sqrt{\Omega_m/0.3}=0.763\pm0.022$ into accordance with galaxy clustering and weak lensing measurements. It also raises up the Hubble constant $H_0=69.68\pm1.00{\rm \,\,km\,s^{-1}Mpc^{-1}}$ that reduces the Hubble tension to the $2.5\sigma$ level. We examine the residual tension in the Early Dark Energy (EDE) model which produces the brief energy injection prior to recombination. We implement both the background and perturbation evolutions of the scalar field which potential scales as $V(\phi)\propto \phi^{2n}$. Including cosmic shear measurements (KiDS, VIKING-450, DES) and local distance-ladder data (SH0ES) to the combined fit we found that EDE completely alleviates the Hubble tension while not degradating the fit to large-scale structure data. The EDE scenario significantly improves the goodness-of-fit by $2.9\sigma$ in comparison with the concordance $\Lambda$CDM model. The account for the intermediate-redshift data (the supernova dataset and baryon acoustic oscillation data) fits perfectly to our parameter predictions and indicates the preference of EDE over $\Lambda$CDM at $3\sigma$.

We present a new open-source code that calculates one-loop power spectra and cross spectra for matter fields and biased tracers in real and redshift space. These spectra incorporate all ingredients required for a direct application to data: non-linear bias and redshift-space distortions, infra-red resummation, counterterms, and the Alcock-Paczynski effect. Our code is based on the Boltzmann solver CLASS and inherits its advantageous properties: user friendliness, ease of modification, high speed, and simple interface with other software. We present detailed descriptions of the theoretical model, the code structure, approximations, and accuracy tests. A typical end-to-end run for one cosmology takes 0.3 seconds, which is sufficient for Markov Chain Monte Carlo parameter extraction. As an example, we apply the code to the Baryon Oscillation Spectroscopic Survey (BOSS) data and infer cosmological parameters from the shape of the galaxy power spectrum.

We present a Markov-Chain Monte-Carlo (MCMC) forecast for the precision of neutrino mass and cosmological parameter measurements with a Euclid-like galaxy clustering survey. We use a complete perturbation theory model for the galaxy one-loop power spectrum and tree-level bispectrum, which includes bias, redshift space distortions, IR resummation for baryon acoustic oscillations and UV counterterms. The latter encapsulate various effects of short-scale dynamics which cannot be modeled within perturbation theory. Our MCMC procedure consistently computes the non-linear power spectra and bispectra as we scan over different cosmologies. The second ingredient of our approach is the theoretical error covariance which captures uncertainties due to higher-order non-linear corrections omitted in our model. Having specified characteristics of a Euclid-like spectroscopic survey, we generate and fit mock galaxy power spectrum and bispectrum likelihoods. Our results suggest that even under very agnostic assumptions about non-linearities and short-scale physics a future Euclid-like survey will be able to measure the sum of neutrino masses with a standard deviation of 28 meV. When combined with the Planck cosmic microwave background likelihood, this uncertainty decreases to 13 meV. Over-optimistically reducing the theoretical error on the bispectrum down to the two-loop level marginally tightens this bound to 11 meV. Moreover, we show that the future large-scale structure (LSS) spectroscopic data will greatly improve constraints on the other cosmological parameters, e.g. reaching a percent (per mille) error on the Hubble constant with LSS alone (LSS + Planck).

We describe two new generation mechanisms for Dark Matter composed of sterile neutrinos with ${\cal O}(1)$ keV mass. The model contains a light scalar field which coherently oscillates in the early Universe and modulates the Majorana mass of the sterile neutrino. In a region of model parameter space, the oscillations between active and sterile neutrinos are resonantly enhanced. This mechanism allows us to produce sterile neutrino DM with small mixing angle with active neutrinos, thus evading the X-ray constraints. At the same time the spectrum of produced DM is much cooler, than in the case of ordinary oscillations in plasma, opening a window of lower mass DM, which is otherwise forbidden by structure formation considerations. In other regions of the model parameter space, where the resonance does not appear, another mechanism can operate: large field suppresses the active-sterile oscillations, but instead sterile neutrinos are produced by the oscillating scalar field when the effective fermion mass crosses zero. In this case DM component is cold, and even 1 keV neutrino is consistent with the cosmic structure formation.

We study a model of a keV-scale sterile neutrino with a relatively large mixing with the Standard Model sector. Usual considerations predict active generation of such particles in the early Universe, which leads to constraints from the total Dark Matter density and absence of X-ray signal from sterile neutrino decay. These bounds together may deem any attempt of creation of the keV scale sterile neutrino in the laboratory unfeasible. We argue that for models with a hidden sector coupled to the sterile neutrino these bounds can be evaded, opening new perspectives for the direct studies at neutrino experiments such as Troitsk $\nu$-mass and KATRIN. We estimate the generation of sterile neutrinos in scenarios with the hidden sector dynamics keeping the sterile neutrinos either massless or superheavy in the early Universe. In both cases the generation by oscillations from active neutrinos in plasma is suppressed.

It has been recently suggested~\cite{Berezhiani:2015yta} that a subdominant fraction of dark matter decaying after recombination may alleviate tension between high-redshift (CMB anisotropy) and low-redshift (Hubble constant, cluster counts) measurements. In this report, we continue our previous study~\cite{Chudaykin:2016yfk} of the decaying dark matter (DDM) model adding all available recent baryon acoustic oscillation (BAO) and redshift space distortions (RSD) measurements. We find, that the BAO/RSD measurements generically prefer the standard $\Lambda$CDM and combined with other cosmological measurements impose an upper limit on the DDM fraction at the level of $\sim$\,5\,\%, strengthening by a factor of 1.5 limits obtained in \cite{Chudaykin:2016yfk} mostly from CMB data. However, the numbers vary from one analysis to other based on the same Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) galaxy sample. Overall, the model with a few percent DDM fraction provides a better fit to the combined cosmological data as compared to the $\Lambda$CDM: the cluster counting and direct measurements of the Hubble parameter are responsible for that. The improvement can be as large as 1.5\,$\sigma$ and grows to 3.3\,$\sigma$ when the CMB lensing power amplitude ${\rm A_L}$ is introduced as a free fitting parameter.

It has been recently suggested~\cite{Berezhiani:2015yta} that emerging tension between cosmological parameter values derived in high-redshift (CMB anisotropy) and low-redshift (cluster counts, Hubble constant) measurements can be reconciled in a model which contains subdominant fraction of dark matter decaying after recombination. We check the model against the CMB Planck data. We find that lensing of the CMB anisotropies by the large-scale structure gives strong extra constraints on this model, limiting the fraction as $F<8\%$ at 2\,$\sigma$ confidence level. However, investigating the combined data set of CMB and conflicting low-$z$ measurements, we obtain that the model with $F\approx2\!-\!5$\% exhibits better fit (by 1.5-3\,$\sigma$ depending on the lensing priors) compared to that of the concordance $\Lambda$CDM cosmological model.

We address the cosmological role of an additional ${\cal O}(1)$ eV sterile neutrino in modified gravity models. We confront the present cosmological data with predictions of the FLRW cosmological model based on a variant of $f(R)$ modified gravity proposed by one of the authors previously. This viable cosmological model which deviation from general relativity with a cosmological constant $\Lambda$ decreases as $R^{-2n}$ for large, but not too large values of the Ricci scalar $R$ provides an alternative explanation of present dark energy and the accelerated expansion of the Universe. Various up-to-date cosmological data sets exploited include Planck CMB anisotropy, CMB lensing potential, BAO, cluster mass function and Hubble constant measurements. We find that the CMB+BAO constraints strongly the sum of neutrino masses from above. This excludes values $\lambda\sim 1$ for which distinctive cosmological features of the model are mostly pronounced as compared to the $\Lambda$CDM model, since then free streaming damping of perturbations due to neutrino rest masses is not sufficient to compensate their extra growth occurring in $f(R)$ gravity. Thus, we obtain $\lambda>8.2$ ($2\sigma$) with cluster systematics and $\lambda>9.4$ ($2\sigma$) without that. In the latter case we find for the sterile neutrino mass $0.47\,\,\rm{eV}$$\,<\,$$m_{\nu,\,\rm{sterile}}$$\,<\,$$1\,\,\rm{eV}$ ($2\sigma$) assuming the active neutrinos are massless, not significantly larger than in the standard $\Lambda$CDM with the same data set: $0.45\,\,\rm{eV}$$\,<\,$$m_{\nu,\,\rm{sterile}}$$\,<\,$$0.92\,\,\rm{eV}$ ($2\sigma$). However, a possible discovery of a sterile neutrino with the mass $m_{\nu,\,\rm{sterile}} \approx 1.5\,$eV motivated by various anomalies in neutrino oscillation experiments would favor cosmology based on $f(R)$ gravity rather than the $\Lambda$CDM model.