The experiment to detect the global epoch of reionization signature (EDGES) collaboration reported the detection of a line at 78MHz in the sky-averaged spectrum due to neutral hydrogen (HI) 21-cm hyperfine absorption of cosmic microwave background photons at $z\sim 17$. This requires that the spin temperature of HI be coupled to the kinetic temperature of the gas at this $z$ through the scattering of Ly-$\alpha$ photons emitted by massive stars. To explain the experimental result, star formation needs to be sufficiently efficient at $z\sim 17$ and this can be used to constrain models in which small-scale structure formation is suppressed (DMF models), either due to dark matter free-streaming or non-standard inflationary dynamics. We combine simulations of structure formation with a simple recipe for star formation to investigate whether these models emit enough Ly-$\alpha$ photons to reproduce the experimental signal for reasonable values of the star formation efficiency, $f_\star$. We find that a thermal warm dark matter (WDM) model with mass $m_\mathrm{WDM}\sim 4.3\,\mathrm{keV}$ is consistent with the timing of the signal for $f_\star\lesssim 2\%$. The exponential growth of structure around $z\sim 17$ in such a model naturally generates a sharp onset of the absorption. A warmer model with $m_\mathrm{WDM}\sim 3\,\mathrm{keV}$ requires a higher star formation efficiency, $f_\star\sim 6\%$, which is a factor of few above predictions of current star formation models and observations of satellites in the Milky Way. However, uncertainties in the process of star formation at these $z$ do not allow to derive strong constrains on such models using 21-cm absorption line. The onset of the 21-cm absorption is generally faster in DMF compared to cold dark matter (CDM), unless some process significantly suppresses star formation in halos with masses below $\sim 10^8\,h^{-1}\,\mathrm{M}_\odot$.

Future 21 cm intensity mapping surveys such as SKA can provide precise information on the spatial distribution of the neutral hydrogen (HI) in the post-reionization epoch. This information will allow to test the standard LCDM paradigm and with that the nature of gravity. In this work, we employ the SHYBONE simulations, which model galaxy formation in f(R) modified gravity using the IllustrisTNG model, to study the effects of modified gravity on HI abundance and power spectra. We find that the enhanced growth low-mass dark matter halos experience in f(R) gravity at high redshifts alters the HI power spectrum and can be observable through 21 cm intensity mapping. Our results suggest that the HI power spectrum is suppressed by 13\% on scales $k\lesssim2\,h\,\mathrm{Mpc}^{-1}$ at z=2 for F6, a f(R) model which passes most observational constraints. We show that this suppression can be detectable by SKA1-MID with 1000 hours of exposure time, making HI clustering a novel test of gravity at high redshift.

We have carried out a set of cosmological hydrodynamical simulations that follow galaxy formation in $f(R)$ modified gravity models. Our simulations employ the Illustris-TNG full physics model and a new modified gravity solver in the AREPO code. For the first time we are able to investigate the degeneracy in the matter power spectrum between the effects of $f(R)$-gravity and feedback from active galactic nuclei (AGN), and the imprint of modified gravity on the properties of galaxies and on the distribution of dark matter, gas and stars in the universe. $f(R)$-gravity has an observable effect on the neutral hydrogen power spectrum at high redshift at a level of 20%. For both the F6 and F5 models, this is significantly larger than the predicted errors for the SKA1-MID survey, making this probe a powerful test of gravity on large scales. A similar effect is present in the power spectrum of the stars at high redshift. We also show that rotationally supported disc galaxies can form in $f(R)$-gravity, even in the partially screened regime. Our simulations indicate that there might be more disc galaxies in F6 compared to GR, and fewer in F5. Finally, we show that the back reaction between AGN feedback and modified gravity in the matter power spectrum is not important in the F6 model but has a sizeable effect in F5.

Thermal inflation models (which feature two inflationary stages) can display damped primordial curvature power spectra on small scales which generate damped matter fluctuations. For a reasonable choice of parameters, thermal inflation models naturally predict a suppression of the matter power spectrum on galactic and sub-galactic scales, mimicking the effect of warm or interacting dark matter. Matter power spectra in these models are also characterised by an excess of power (w.r.t. the standard $\Lambda$CDM power spectrum) just below the suppression scale. By running a suite of N-body simulations we investigate the non-linear growth of structure in models of thermal inflation. We measure the non-linear matter power spectrum and extract halo statistics, such as the halo mass function, and compare these quantities with those predicted in the standard $\Lambda$CDM model and in other models with damped matter fluctuations. We find that the thermal inflation models considered here produce measurable differences in the matter power spectrum from $\Lambda$CDM at redshifts $z>5$, while the halo mass functions are appreciably different at all redshifts. The halo mass function at $z=0$ for thermal inflation displays an enhancement of around $\sim 20\%$ w.r.t. $\Lambda$CDM and a damping at lower halo masses, with the position of the enhancement depending on the value of the free parameter in the model. The enhancement in the halo mass function (w.r.t. $\Lambda$CDM ) increases with redshift, reaching $\sim 40\%$ at $z=5$. We also study the accuracy of the analytical Press-Schechter approach, using different filters to smooth the density field, to predict halo statistics for thermal inflation. We find that the predictions with the smooth-$k$ filter agree with the simulation results over a wider range of halo masses than is the case with other filters commonly used in the literature.

We investigate the nonlinear evolution of structure in variants of the standard cosmological model which display damped density fluctuations relative to cold dark matter (e.g. in which cold dark matter is replaced by warm or interacting DM). Using N-body simulations, we address the question of how much information is retained from different scales in the initial linear power spectrum following the nonlinear growth of structure. We run a suite of N-body simulations with different initial linear matter power spectra to show that, once the system undergoes nonlinear evolution, the shape of the linear power spectrum at high wavenumbers does not affect the non-linear power spectrum, while it still matters for the halo mass function. Indeed, we find that linear power spectra which differ from one another only at wavenumbers larger than their half-mode wavenumber give rise to (almost) identical nonlinear power spectra at late times, regardless of the fact that they originate from different models with damped fluctuations. On the other hand, the halo mass function is more sensitive to the form of the linear power spectrum. Exploiting this result, we propose a two parameter model of the transfer function in generic damped scenarios, and show that this parametrisation works as well as the standard three parameter models for the scales on which the linear spectrum is relevant.

We propose a new filter, a smooth-$k$ space filter, to use in the Press-Schechter approach to model the dark matter halo mass function which overcomes shortcomings of other filters. We test this against the mass function measured in N-body simulations. We find that the commonly used sharp-$k$ filter fails to reproduce the behaviour of the halo mass function at low masses measured from simulations of models with a sharp truncation in the linear power spectrum. We show that the predictions with our new filter agree with the simulation results over a wider range of halo masses for both damped and undamped power spectra than is the case with the sharp-$k$ and real-space top-hat filters.

We investigate the impact of thermal velocities in N-body simulations of structure formation in warm dark matter models. Adopting the commonly used approach of adding thermal velocities, randomly selected from a Fermi-Dirac distribution, to the gravitationally-induced velocities of the simulation particles, we compare the matter and velocity power spectra measured from CDM and WDM simulations, in the latter case with and without thermal velocities. This prescription for adding thermal velocities introduces numerical noise into the initial conditions, which influences structure formation. At early times, the noise affects dramatically the power spectra measured from simulations with thermal velocities, with deviations of the order of $\sim \mathcal{O}(10)$ (in the matter power spectra) and of the order of $\sim \mathcal{O}(10^2)$ (in the velocity power spectra) compared to those extracted from simulations without thermal velocities. At late times, these effects are less pronounced with deviations of less than a few percent. Increasing the resolution of the N-body simulation shifts these discrepancies to higher wavenumbers. We also find that spurious haloes start to appear in simulations which include thermal velocities at a mass that is $\sim$3 times larger than in simulations without thermal velocities.

Recently it has been pointed out that a cosmic background of relativistic axion-like particles (ALPs) would be produced by the primordial decays of heavy fields in the post-inflation epoch, contributing to the extra-radiation content in the Universe today. Primordial magnetic fields would trigger conversions of these ALPs into sub-MeV photons during the dark ages. This photon flux would produce an early reionization of the Universe, leaving a significant imprint on the total optical depth to recombination $\tau$. Using the current measurement of $\tau$ and the limit on the extra-radiation content $\Delta N_{\rm eff} $ by the Planck experiment we put a strong bound on the ALP-photon conversions. Namely we obtain upper limits on the product of the photon-ALP coupling constant $g_{a\gamma}$ times the magnetic field strength $B$ down to $g_{a\gamma} B \gtrsim 6 \times 10^{-18} \textrm{GeV}^{-1} \textrm{nG} $ for ultralight ALPs.