21 cm intensity mapping has emerged as a promising technique to map the large-scale structure of the Universe. However, the presence of foregrounds with amplitudes orders of magnitude larger than the cosmological signal constitutes a critical challenge. Here, we test the sparsity-based algorithm Generalised Morphological Component Analysis (GMCA) as a blind component separation technique for this class of experiments. We test the GMCA performance against realistic full-sky mock temperature maps that include, besides astrophysical foregrounds, also a fraction of the polarized part of the signal leaked into the unpolarized one, a very troublesome foreground to subtract, usually referred to as polarization leakage. To our knowledge, this is the first time the removal of such component is performed with no prior assumption. We assess the success of the cleaning by comparing the true and recovered power spectra, in the angular and radial directions. In the best scenario looked at, GMCA is able to recover the input angular (radial) power spectrum with an average bias of $\sim 5\%$ for $\ell>25$ ($20 - 30 \%$ for $k_{\parallel} \gtrsim 0.02 \,h^{-1}$Mpc), in the presence of polarization leakage. Our results are robust also when up to $40\%$ of channels are missing, mimicking a Radio Frequency Interference (RFI) flagging of the data. In perspective, we endorse the improvement on both cleaning methods and data simulations, the second being more and more realistic and challenging the first ones, to make 21 cm intensity mapping competitive.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

We perform an analysis of the cosmic reionization in the standard cold dark matter (CDM) paradigm and in alternative dark matter scenarios. Building upon the work of Corasaniti et al. (2017), we predict the reionization history for CDM, for warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ADM) models which reproduce state-of-art measurements of the galaxy luminosity function at very high-redshifts $6\le z\le 10$. To this purpose we adopt a reionization model parametrized in terms of the limiting UV-magnitude of galaxies contributing to the reionization $M_{\rm lim}$ and the average effective escape fraction of UV photons reaching the intergalactic medium $\tilde{f}$. For each DM model we compute the redshift evolution of the Thomson scattering optical depth $\tau_e(z)$ and the comoving ionization fraction $Q_{\rm HII}(z)$. We find the different DM models to have similar reionization histories. Differences with respect to the CDM case increase at fainter limiting UV-magnitudes and are degenerate with the effect of varying the reionization model parameters. Using Planck's determination of the integrated optical depth in combination with measurements of the neutral hydrogen fraction at different redshifts, we infer constraints on $\tilde{f}$ and $M_{\rm lim}$. The results are largely independent of the assumed DM scenario, in particular for $M_{\rm lim}\gtrsim -13$ we obtain that the effective escape fraction lies in the range $0.07\lesssim \tilde{f}\lesssim 0.15$ at $2\sigma$.

[Abridged] We study the abundance and clustering properties of HI at redshifts $z\leqslant5$ using TNG100, a large state-of-the-art magneto-hydrodynamic simulation of a 75 Mpc/h box size. We show that most of the HI lies within dark matter halos and quantify the average HI mass hosted by halos of mass M at redshift z. We find that only halos with circular velocities larger than $\simeq$ 30 km/s contain HI. While the density profiles of HI exhibit a large halo-to-halo scatter, the mean profiles are universal across mass and redshift. The HI in low-mass halos is mostly located in the central galaxy, while in massive halos is concentrated in the satellites. We show that the HI and matter density probability distribution functions differ significantly. Our results point out that for small halos the HI bulk velocity goes in the same direction and has the same magnitude as the halo peculiar velocity, while in large halos differences show up. We find that halo HI velocity dispersion follows a power-law with halo mass. We find a complicated HI bias, with HI becoming non-linear already at $k=0.3$ h/Mpc at $z\gtrsim3$. Our simulation reproduces the DLAs bias value from observations. We find that the clustering of HI can be accurately reproduced by perturbative methods. We identify a new secondary bias, by showing that the clustering of halos depends not only on mass but also on HI content. We compute the amplitude of the HI shot-noise and find that it is small at all redshifts. We study the clustering of HI in redshift-space, and show that linear theory can explain the ratio between the monopoles in redshift- and real-space down to small scales at high redshift. We find that the amplitude of the Fingers-of-God effect is larger for HI than for matter. We point out that accurate 21 cm maps can be created from N-body or approximate simulations rather than full hydrodynamic simulations.

The nature of the most abundant components of the Universe, dark energy and dark matter, is still to be uncovered. I tackle this subject considering a novel cosmological probe: the neutral hydrogen emitted 21cm radiation, observed with the intensity mapping technique. I analyse competitive and realistic dark energy and dark matter models and show how they produce distinctive and detectable effects on the 21cm signal. Moreover, I provide radio telescope forecasts showing how these models will be distinguishable in an unprecedented way.

We study the imprint of non-standard dark energy (DE) and dark matter (DM) models on the 21cm intensity map power spectra from high-redshift neutral hydrogen (HI) gas. To this purpose we use halo catalogs from N-body simulations of dynamical DE models and DM scenarios which are statistically indistinguishable from the standard Cold Dark Matter model with Cosmological Constant (LCDM) using currently available cosmological observations. We limit our analysis to halo catalogs at redshift z = 1 and 2.3 which are common to all simulations. For each catalog we model the HI distribution by using a simple prescription to associate the HI gas mass to N-body halos. We find that the DE models leave a distinct signature on the HI spectra across a wide range of scales, which correlates with differences in the halo mass function and the onset of the non-linear regime of clustering. In the case of the non-standard DM model significant differences of the HI spectra with respect to the LCDM model only arise from the suppressed abundance of low mass halos. These cosmological model dependent features also appear in the 21cm spectra. In particular, we find that future SKA measurements can distinguish the imprints of DE and DM models at high statistical significance.

We investigate the cross-correlation signal between 21cm intensity mapping maps and the Lyman-alpha forest in the fully non-linear regime using state-of-the-art hydrodynamic simulations. The cross-correlation signal between these fields can provide a coherent and comprehensive picture of the neutral hydrogen (HI) content of our Universe in the post-reionization era, probing both its mass content and volume distribution. We compute the auto-power spectra of both fields together with their cross-power spectrum at z = 2.4 and find that on large scales the fields are completely anti-correlated. This anti-correlation arises because regions with high (low) 21cm emission, such as those with a large (low) concentration of damped Lyman-alpha systems, will show up as regions with low (high) transmitted flux. We find that on scales smaller than k = 0.2 h/Mpc the cross-correlation coefficient departs from -1, at a scale where non-linearities show up. We use the anisotropy of the power spectra in redshift-space to determine the values of the bias and of the redshift-space distortion parameters of both fields; we find that the errors on the value of these parameters could decrease by 30% when adding data from the cross-power spectrum in a conservative analysis. Our results point out that linear theory is capable of reproducing the shape and amplitude of the cross-power up to rather non-linear scales. Finally, we find that the 21cm-Lya cross-power spectrum can be detected by combining data from a BOSS-like survey together with 21cm intensity mapping observations by SKA1-MID with a S/N ratio higher than 3 in the range 0.06< k <1 h/Mpc. We emphasize that while the shape and amplitude of the 21cm auto-power spectrum can be severely affected by residual foreground contamination, cross-power spectra will be less sensitive to that and therefore can be used to identify systematics in the 21cm maps.

We investigate the impact that warm dark matter (WDM) has in terms of 21cm intensity mapping in the post-reionization Universe at z = 3 - 5. We perform hydrodynamic simulations for 5 different models: cold dark matter and WDM with 1,2,3,4 keV (thermal relic) mass and assign the neutral hydrogen a-posteriori using two different methods that both reproduce observations in terms of column density distribution function of neutral hydrogen systems. Contrary to naive expectations, the suppression of power present in the linear and non-linear matter power spectra, results in an increase of power in terms of neutral hydrogen and 21cm power spectra. This is due to the fact that there is a lack of small mass halos in WDM models with respect to cold dark matter: in order to distribute a total amount of neutral hydrogen within the two cosmological models, a larger quantity has to be placed in the most massive halos, that are more biased compared to the cold dark matter cosmology. We quantify this effect and address significance for the telescope SKA1-LOW, including a realistic noise modeling. The results indicate that we will be able to rule out a 4 keV WDM model with 5000 hours of observations at z > 3, with a statistical significance of > 3 sigma, while a smaller mass of 3 keV, comparable to present day constraints, can be ruled out at more than 2 sigma confidence level with 1000 hours of observations at z > 5.

Dark matter halos are built from accretion and merging. During merging some of the dark matter particles may be ejected with velocities higher than the escape velocity. We use both N-body simulations and single-particle smooth-field simulations to demonstrate that rapid changes to the mean field potential are responsible for such ejection, and in particular that dynamical friction plays no significant role in it. Studying a range of minor mergers, we find that typically between 5-15% of the particles from the smaller of the two merging structures are ejected. We also find that the ejected particles originate essentially from the small halo, and more specifically are particles in the small halo which pass later through the region in which the merging occurs.

In electrovacuum stationary, asymptotically flat black holes in scalar-tensor theories of gravity are described by the Kerr-Newman family of solutions, just as in general relativity. We show that there exist two mechanisms which can render Kerr black holes unstable when matter is present in the vicinity of the black hole, as this induces an effective mass for the scalar. The first mechanism is a tachyonic instability that appears when the effective mass squared is negative, triggering the development of scalar hair --- a black hole version of "spontaneous scalarization". The second instability is associated with superradiance and is present when the effective mass squared is positive and when the black hole spin exceeds a certain threshold. The second mechanism is also responsible for a resonant effect in the superradiant scattering of scalar waves, with amplification factors as large as 10^5 or more.

We uncover two mechanisms that can render Kerr black holes unstable in scalar-tensor gravity, both associated to the presence of matter in the vicinity of the black hole and the fact that this introduces an effective mass for the scalar. Our results highlight the importance of understanding the structure of spacetime in realistic, astrophysical black holes in scalar-tensor theories.