Short baseline neutrino experiments, like LSND and MiniBooNE experiments, pointed towards the existence of eV mass scale sterile neutrinos. To reconcile sterile neutrinos with cosmology self interaction between sterile neutrinos has been studied. We analysed Planck cosmic microwave background (CMB) data with self-interacting sterile neutrino (SI$\nu$) and study their impact on inflation models. The fit to the CMB data in SI$\nu$ model is as good as the fit to $\Lambda$CDM model. We find that the spectral index ($n_s$) values shift to $0.9361\pm 0.0055$ in SI$\nu$ model. This has significant impact on the validity of different inflation models. For example the Starobinsky and quartic hilltop model, which were allowed within $\Lambda$CDM cosmology, are ruled out. On the other hand some models like natural and Coleman-Weinberg inflation are now favoured. Therefore, the existence of self interacting sterile neutrinos with eV order of mass will play an important role in the selection of correct inflation model.

For primordial black holes (PBHs) to form a considerable fraction of cold dark matter, the required amplitude of primordial scalar perturbations is quite large ($P_{\zeta}(k) \sim 10^{-2}$) if PBH is formed in radiation epoch. In alternate cosmological histories, where additional epochs of arbitrary equation of state precede radiation epoch, the dynamics of PBH formation and relevant mass ranges can be different leading to requirement of lower primordial power at smaller scales of inflation. Moreover, this alternate history can modify the predictions for the gravitational wave (GW) spectrum, which can be probed by upcoming GW observations. In this paper we show that an early kination epoch can lead to percent level abundance of PBH for a lower amplitude of $P_{\zeta}(k)$ as compared to PBH formation in a standard radiation epoch. Moreover, we calculate the effect of early kination epoch on the GW spectrum for first and second orders in perturbation which show enhancement in the amplitude of the GW spectrum in a kination epoch with respect to that in a standard radiation epoch.

It has been pointed out that there exists a tension in $\sigma_8-\Omega_m$ measurement between CMB and LSS observation. In this paper we show that $\sigma_8-\Omega_m$ observations can be used to test the dark energy theories. We study two models, (1) Hu-Sawicki(HS) Model of $f(R)$ gravity and (2) Chavallier-Polarski-Linder(CPL) parametrization of dynamical dark energy (DDE), both of which satisfy the constraints from supernovae. We compute $\sigma_8$ consistent with the parameters of these models. We find that the well known tension in $\sigma_8$ between Planck CMB and large scale structure (LSS) observations is (1) exacerbated in the HS model and (2) somewhat alleviated in the DDE model. We illustrate the importance of the $\sigma_8$ measurements for testing modified gravity models. Modified gravity models change the matter power spectrum at cluster scale which also depends upon the neutrino mass. We present the bound on neutrino mass in the HS and DDE model.

Effective field theory of dark matter fluid on large scales predicts the presence of viscosity of the order of $10^{-6} H_0 M_P^2$. It has been shown that this magnitude of viscosities can resolve the discordance between large scale structure observations and Planck CMB data in the $\sigma_8$-$\Omega_m^0$ and $H_0$-$\Omega_m^0$ parameters space. Massive neutrinos suppresses the matter power spectrum on the small length scales similar to the viscosities. We show that by including the effective viscosity, which arises from summing over non linear perturbations at small length scales, severely constrains the cosmological bound on neutrino masses. Under a joint analysis of Planck CMB and different large scale observation data, we find that upper bound on the sum of the neutrino masses at 2-$\sigma$ level, decreases from $\sum m_\nu \le 0.396\,$eV (normal hierarchy) and $\sum m_\nu \le 0.378 \,$eV (inverted hierarchy) to $\sum m_\nu \le 0.267\,$eV (normal hierarchy) and $\sum m_\nu \le 0.146\,$eV (inverted hierarchy) when the effective viscosities are included.