Search

Removing contaminants is a delicate yet crucial step in neutral hydrogen (HI) intensity mapping, often considered the technique's greatest challenge. Here, we address this challenge by analysing HI intensity maps of about $100$ deg$^2$ at redshift $z\approx0.4$ collected by the MeerKAT radio telescope, a SKA Observatory (SKAO) precursor, with a combined 10.5-hour observation. Using unsupervised statistical methods, we remove the contaminating foreground emission and systematically test step-by-step common pre-processing choices to facilitate the cleaning process. We also introduce and test a novel multiscale approach, where data is redundantly decomposed into subsets referring to different spatial scales (large and small), and the cleaning procedure is performed independently. We confirm the detection of the HI cosmological signal in cross-correlation with an ancillary galactic data set without the need to correct for signal loss. In the best set-up reached, we constrain the HI distribution through the combination of its cosmic abundance ($\Omega_{\rm HI}$) and linear clustering bias ($b_{\rm HI}$) up to a cross-correlation coefficient ($r$) and measure $\Omega_{\rm HI}b_{\rm HI}r = [0.93 \pm 0.17]\,\times\,10^{-3}$ with $\approx6\sigma$ confidence. The measurement is independent of scale cuts at both edges of the probed scale range ($0.04 \lesssim k \lesssim 0.3 \,h$Mpc$^{-1}$), corroborating its robustness. Our new pipeline has successfully found an optimal compromise in separating contaminants without incurring a catastrophic signal loss, instilling more confidence in the outstanding science we can deliver with MeerKAT on the path towards HI intensity mapping surveys with the full SKAO., Comment: 24 pages (+ appendix), 24 figures. Summary of results in Table 1 and Figure 17 on page 16. Comments welcome

We present results from MeerKAT single-dish HI intensity maps, the final observations to be performed in L-band in the MeerKAT Large Area Synoptic Survey (MeerKLASS) campaign. The observations represent the deepest single-dish HI intensity maps to date, produced from 41 repeated scans over $236\,{\rm deg}^2$, providing 62 hours of observational data for each of the 64 dishes before flagging. By introducing an iterative self-calibration process, the estimated thermal noise of the reconstructed maps is limited to ${\sim}\,1.21\,$mK ($1.2\,\times$ the theoretical noise level). This thermal noise will be sub-dominant relative to the HI fluctuations on large scales ($k\,{\lesssim}\,0.15\,h\,\text{Mpc}^{-1}$), which demands upgrades to power spectrum analysis techniques, particularly for covariance estimation. In this work, we present the improved MeerKLASS analysis pipeline, validating it on both a suite of mock simulations and a small sample of overlapping spectroscopic galaxies from the Galaxy And Mass Assembly (GAMA) survey. Despite only overlapping with ${\sim}\,25\%$ of the MeerKLASS deep field, and a conservative approach to covariance estimation, we still obtain a ${>}\,4\,\sigma$ detection of the cross-power spectrum between the intensity maps and the 2269 galaxies at the narrow redshift range $0.39\,{<}\,z\,{<}\,0.46$. We briefly discuss the HI auto-power spectrum from this data, the detection of which will be the focus of follow-up work. For the first time with MeerKAT single-dish intensity maps, we also present evidence of HI emission from stacking the maps onto the positions of the GAMA galaxies., Comment: 30 pages, 27 figures. See Fig.14 and 15 for auto-HI and HI-galaxy cross power spectra. Fig 19 and 20 for stacking detections. Accepted to MNRAS

Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital sources in the Radio Navigational Satellite System within the 1100-1350 MHz frequency. This simulation can be split into two parts: (I) satellite positioning, emission power, and beam response on the telescope and (II) fitting of the satellite signal to data in order to improve the original model. We use previously observed single dish MeerKAT L-band data which needs to be specially calibrated to include data contaminated by satellite-based RFI. We find that due to non-linearity effects, it becomes non-trivial to fit the satellite power. However, when masking regions where this non-linearity is problematic, we can recreate the satellite contamination with high accuracy around its peak frequencies. The simulation can predict satellite movements and signal for past and future observations, which can help in RFI avoidance and testing novel cleaning methods. The predicted signal from simulations sits below the noise in the target cosmology window for the L-band (970 - 1015 MHz) making it difficult to confirm any out-of-band emission from satellites. However, a power spectrum analysis shows that such signal can still contaminate the 21cm power spectrum at these frequencies. In our simulations, this contamination overwhelms the auto-power spectrum but still allows for a clean detection of the signal in cross-correlations with mild foreground cleaning. Whether such contamination does exist one will require further characterization of the satellite signals far away from their peak frequencies., Comment: 20 pages, 27 figures, 4 tables, Comments welcome

Direct detection of the Cosmic Dawn and Epoch of Reionization via the redshifted 21-cm line of neutral Hydrogen will have unprecedented implications for studying structure formation in the early Universe. This exciting goal is challenged by the difficulty of extracting the faint 21-cm signal buried beneath bright astrophysical foregrounds and contaminated by numerous systematics. Here, we focus on improving the Gaussian Process Regression (GPR) signal separation method originally developed for LOFAR observations. We address a key limitation of the current approach by incorporating covariance prior models learnt from 21-cm signal simulations using Variational Autoencoder (VAE) and Interpolatory Autoencoder (IAE). Extensive tests are conducted to evaluate GPR, VAE-GPR, and IAE-GPR in different scenarios. Our findings reveal that the new method outperforms standard GPR in component separation tasks. Moreover, the improved method demonstrates robustness when applied to signals not represented in the training set. It also presents a certain degree of resilience to data systematics, highlighting its ability to effectively mitigate their impact on the signal recovery process. However, our findings also underscore the importance of accurately characterizing and understanding these systematics to achieve successful detection. Our generative approaches provide good results even with limited training data, offering a valuable advantage when a large training set is not feasible. Comparing the two algorithms, IAE-GPR shows slightly higher fidelity in recovering power spectra compared to VAE-GPR. These advancements highlight the strength of generative approaches and optimise the analysis techniques for future 21-cm signal detection at high redshifts., Comment: 14 pages, 9 figues, submitted to MNRAS, any comments welcome

Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensity mapping experiments. However, assessing whether this approach is viable for future intensity mapping surveys where precision cosmology is the aim, remains unexplored. In this work, using simulations, we validate for the first time the use of a foreground transfer function to reconstruct power spectra of foreground-cleaned low-redshift intensity maps and look to expose any limitations. We reveal that even when aggressive foreground cleaning is required, which causes ${>}\,50\%$ negative bias on the largest scales, the power spectrum can be reconstructed using a transfer function to within sub-percent accuracy. We specifically outline the recipe for constructing an unbiased transfer function, highlighting the pitfalls if one deviates from this recipe, and also correctly identify how a transfer function should be applied in an auto-correlation power spectrum. We validate a method that utilises the transfer function variance for error estimation in foreground-cleaned power spectra. Finally, we demonstrate how incorrect fiducial parameter assumptions (up to ${\pm}100\%$ bias) in the generation of mocks, used in the construction of the transfer function, do not significantly bias signal reconstruction or parameter inference (inducing ${<}\,5\%$ bias in recovered values)., Comment: 25 pages, 20 figures. See Figure 4 for the main demonstration of the transfer function's performance for reconstructing signal loss from foreground cleaning. Accepted for publication in MNRAS

In this work we present updated forecasts on parameterised modifications of gravity that can capture deviations of the behaviour of cosmological density perturbations beyond $\Lambda$CDM. For these forecasts we adopt the SKA Observatory (SKAO) as a benchmark for future cosmological surveys at radio frequencies, combining a continuum survey for weak lensing and angular galaxy clustering with an HI galaxy survey for spectroscopic galaxy clustering that can detect baryon acoustic oscillations and redshift space distortions. Moreover, we also add 21cm HI intensity mapping, which provides invaluable information at higher redshifts, and can complement tomographic resolution, thus allowing us to probe redshift-dependent deviations of modified gravity models. For some of these cases, we combine the probes with other optical surveys, such as the Dark Energy Spectroscopic Instrument (DESI) and the Vera C. Rubin Observatory (VRO). We show that such synergies are powerful tools to remove systematic effects and degeneracies in the non-linear and small-scale modelling of the observables. Overall, we find that the combination of all SKAO radio probes will have the ability to constrain the present value of the functions parameterising deviations from $\Lambda$CDM ($\mu$ and $\Sigma$) with a precision of $2.7\%$ and $1.8\%$ respectively, competitive with the constraints expected from optical surveys and with constraints we have on gravitational interactions in the standard model. Exploring the radio-optical synergies, we find that the combination of VRO with SKAO can yield extremely tight constraints on $\mu$ and $\Sigma$ ($0.9\%$ and $0.7\%$ respectively), which are further improved when the cross-correlation between intensity mapping and DESI galaxies is included., Comment: 24 pages, 10 figures

We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7\sigma$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($\Omega_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constraint $\Omega_{\rm HI} b_{\rm HI} r\,{=}\,[0.86\,{\pm}\,0.10\,({\rm stat})\,{\pm}\,0.12\,({\rm sys})]\,{\times}\,10^{-3}$, at an effective scale of $k_{\rm eff}\,{\sim}\,0.13\,h\,{\rm Mpc}^{-1}$. The intensity maps were obtained from a pilot survey with the MeerKAT telescope, a 64-dish pathfinder array to the SKA Observatory (SKAO). The data were collected from 10.5 hours of observations using MeerKAT's L-band receivers over six nights covering the 11hr field of WiggleZ, in the frequency range $1015-973\,{\rm MHz}$ (0.400$\,{<}\,z\,{<}\,$0.459 in redshift). This detection is the first practical demonstration of the multi-dish auto-correlation intensity mapping technique for cosmology. This marks an important milestone in the roadmap for the cosmology science case with the full SKAO., Comment: 11 pages, 7 figures. See Fig.5 for the main result, showing the cross-power spectrum between MeerKAT and WiggleZ, revealing the 7.7\sigma detection. v2 is the version accepted for publication by MNRAS. The main changes from v1 were in introducing the radio intensity mapping data

The cross-correlation between 21-cm intensity mapping experiments and photometric surveys of galaxies (or any other cosmological tracer with a broad radial kernel) is severely degraded by the loss of long-wavelength radial modes due to Galactic foreground contamination. Higher-order correlators are able to restore some of these modes due to the non-linear coupling between them and the local small-scale clustering induced by gravitational collapse. We explore the possibility of recovering information from the bispectrum between a photometric galaxy sample and an intensity mapping experiment, in the context of the clustering-redshifts technique. We demonstrate that the bispectrum is able to calibrate the redshift distribution of the photometric sample to the required accuracy of future experiments such as the Rubin Observatory, using future single-dish and interferometric 21-cm observations, in situations where the two-point function is not able to do so due to foreground contamination. We also show how this calibration is affected by the photometric redshift width $\sigma_{z,0}$ and maximum scale $k_{\mathrm{max}}$. We find that it is important to reach scales $k \gtrsim 0.3\,h\,\mathrm{Mpc}^{-1}$, with the constraints saturating at around $k\sim 1\,h\,\mathrm{Mpc}^{-1}$ for next-generation experiments., Comment: 20 pages, 8 figures

Neutral Hydrogen Intensity Mapping (HI IM) surveys will be a powerful new probe of cosmology. However, strong astrophysical foregrounds contaminate the signal and their coupling with instrumental systematics further increases the data cleaning complexity. In this work, we simulate a realistic single-dish HI IM survey of a $5000$~deg$^2$ patch in the $950 - 1400$ MHz range, with both the MID telescope of the SKA Observatory (SKAO) and MeerKAT, its precursor. We include a state-of-the-art HI simulations and explore different foreground models and instrumental effects such as non-homogeneous thermal noise and beam side-lobes. We perform the first Blind Foreground Subtraction Challenge for HI IM on these synthetic data-cubes, aiming to characterise the performance of available foreground cleaning methods with no prior knowledge of the sky components and noise level. Nine foreground cleaning pipelines joined the Challenge, based on statistical source separation algorithms, blind polynomial fitting, and an astrophysical-informed parametric fit to foregrounds. We devise metrics to compare the pipeline performances quantitatively. In general, they can recover the input maps' 2-point statistics within 20 per cent in the range of scales least affected by the telescope beam. However, spurious artefacts appear in the cleaned maps due to interactions between the foreground structure and the beam side-lobes. We conclude that it is fundamental to develop accurate beam deconvolution algorithms and test data post-processing steps carefully before cleaning. This study was performed as part of SKAO preparatory work by the HI IM Focus Group of the SKA Cosmology Science Working Group., Comment: Accepted version for publication in MNRAS

The success of HI intensity mapping is largely dependent on how well 21cm foreground contamination can be controlled. In order to progress our understanding further, we present a range of simulated foreground data from four different $\sim3000$\,deg$^2$ sky regions, with and without effects from polarization leakage. Combining these with underlying cosmological HI simulations creates a range of single-dish intensity mapping test cases that require different foreground treatments. This allows us to conduct the most generalized study to date into 21cm foregrounds and their cleaning techniques for the post-reionization era. We first provide a pedagogical review of the most commonly used blind foreground removal techniques (PCA/SVD, FASTICA, GMCA). We also trial a non-blind parametric fitting technique and discuss potential hybridization of methods. We highlight the similarities and differences in these techniques finding that the blind methods produce near equivalent results, and we explain the fundamental reasons for this. The simulations allow an exact decomposition of the resulting cleaned data and we analyse the contribution from foreground residuals. Our results demonstrate that polarized foreground residuals should be generally subdominant to HI on small scales ($k\gtrsim0.1\,h\,\text{Mpc}^{-1}$). However, on larger scales, results are more region dependent. In some cases, aggressive cleans severely damp HI power but still leave dominant foreground residuals. We also demonstrate the gain from cross-correlations with optical galaxy surveys, where extreme levels of residual foregrounds can be circumvented. However, these residuals still contribute to errors and we discuss the optimal balance between over- and under-cleaning., Comment: 21 pages, 16 figures. Changes made in response to reviewer report. Noted changes from V1: analysis now only performed on two foreground regions (excluding Galactic Plane and Barnard's loop since these were very high emission regions). General results and conclusions unchanged though. Version 2 accepted for publication in MNRAS

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. Having quantified the notable effect of polarisation leakage on our results, in perspective we advocate the use of more realistic simulations when testing 21 cm intensity mapping capabilities., Comment: 16 pages and 17 figures. This manuscript matches the published paper. Code and demonstrative notebooks are available at https://github.com/isab3lla/gmca4im and data-set to reproduce our results is available at https://doi.org/10.5281/zenodo.3991818

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., Comment: Science white paper submitted to Decadal 2020 survey

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., Comment: 8 pages, 4 figures

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., Comment: 5 pages, 1 figure, submitted to the Astro2020 Decadal Survey Science White Paper call

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$., Comment: Summary plot is Fig. 8: constraints on escape fraction and UV limiting magnitude for cold and non-cold dark matter. This version matches the published one

[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., Comment: 41 pages. 30 figures. Analysis tools publicly available at https://github.com/franciscovillaescusa/Pylians/tree/master/HI_Illustris

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., Comment: 6 pages, 3 figures. Proceedings of the 8th Young Researcher Meeting, Cagliari, May 29 to June 1, 2017. To be published in Journal of Physics: Conference Series (IOP)

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., Comment: to be submitted, comments welcome

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., Comment: Accepted for publication in JCAP

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., Comment: 25 pages, 11 figures, comments are welcome

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., Comment: 18 pages, 12 figures. Accepted for publication in JCAP