BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Francisco Villaescusa-Navarro

Francisco Villaescusa-Navarro contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.GAastro-ph.IMhep-phgr-qchep-thMachine LearningApplicationsArtificial IntelligenceComputer Vision

Trust snapshot

Quick read

Trust 21 - Emerging
46works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2026arXiv

CAMELS Environments: The Impact of Local Neighbours on Galaxy Evolution across the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE Simulations

Internal feedback from massive stars and active galactic nuclei (AGN) play a key role in galaxy evolution, but external environmental effects can also strongly influence galaxies. We investigate the impact of environment on galaxy evolution, and its dependence on baryonic physics implementation, using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) spanning a wide range of stellar and AGN feedback implementations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. We show that satellite galaxies are significantly affected by the environment in all simulation models, with their gas fraction and star formation rate (SFR) suppressed in overdense regions compared to similar mass satellites in underdense environments at $z=0$. Central galaxies are less sensitive to environment but tend to show lower gas fraction and SFR in overdense regions at low stellar mass, transitioning to higher gas fraction and SFR for massive galaxies in higher-density environments. Halo baryon fraction ($f_{\rm B}$) and circumgalactic medium mass fraction ($f_{\rm CGM}$) at $z=0$ show clear environmental effects. In SIMBA, low-mass haloes in overdense regions have systematically lower $f_{\rm B}$ and $f_{\rm CGM}$ at fixed halo mass, while Swift-EAGLE haloes in overdense regions have systematically higher $f_{\rm B}$ and $f_{\rm CGM}$ across the full halo mass range, and IllustrisTNG and ASTRID show opposite trends at the low and high mass ends. Environmental effects can flip at higher redshift, with SFR and $f_{\rm B}$ increasing with local density in low-mass haloes before quenching at an increasing overdensity threshold. Our results demonstrate that the impact of environment on galaxy evolution depends significantly on galaxy formation model, and higher-density environments can either suppress or enhance star formation depending on galaxy mass and cosmic epoch.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this &#34;back-reaction&#34; of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel&#39;dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmology with one galaxy?

Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine Learning and Cosmology

Methods based on machine learning have recently made substantial inroads in many corners of cosmology. Through this process, new computational tools, new perspectives on data collection, model development, analysis, and discovery, as well as new communities and educational pathways have emerged. Despite rapid progress, substantial potential at the intersection of cosmology and machine learning remains untapped. In this white paper, we summarize current and ongoing developments relating to the application of machine learning within cosmology and provide a set of recommendations aimed at maximizing the scientific impact of these burgeoning tools over the coming decade through both technical development as well as the fostering of emerging communities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Percent-level constraints on baryonic feedback with spectral distortion measurements

High-significance measurements of the monopole thermal Sunyaev-Zel&#39;dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reionization with Simba: How much does astrophysics matter in modeling cosmic reionization?

Traditional large-scale models of reionization usually employ simple deterministic relations between halo mass and luminosity to predict how reionization proceeds. We here examine the impact on modelling reionization of using more detailed models for the ionizing sources as identified within the $100~{\rm Mpc/h}$ cosmological hydrodynamic simulation {\sc Simba}, coupled with post-processed radiative transfer. Comparing with simple (one-to-one) models, the main difference of using {\sc Simba} sources is the scatter in the relation between dark matter halos and star formation, and hence ionizing emissivity. We find that, at the power spectrum level, the ionization morphology remains mostly unchanged, regardless of the variability in the number of sources or escape fraction. In particular, the power spectrum shape remains unaffected and its amplitude changes slightly by less than 5-10\%, throughout reionization, depending on the scale and neutral fraction. Our results show that simplified models of ionizing sources remain viable to efficiently model the structure of reionization on cosmological scales, although the precise progress of reionization requires accounting for the scatter induced by astrophysical effects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robust field-level inference with dark matter halos

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Wavelet Moments for Cosmological Parameter Estimation

Extracting non-Gaussian information from the non-linear regime of structure formation is key to fully exploiting the rich data from upcoming cosmological surveys probing the large-scale structure of the universe. However, due to theoretical and computational complexities, this remains one of the main challenges in analyzing observational data. We present a set of summary statistics for cosmological matter fields based on 3D wavelets to tackle this challenge. These statistics are computed as the spatial average of the complex modulus of the 3D wavelet transform raised to a power $q$ and are therefore known as invariant wavelet moments. The 3D wavelets are constructed to be radially band-limited and separable on a spherical polar grid and come in three types: isotropic, oriented, and harmonic. In the Fisher forecast framework, we evaluate the performance of these summary statistics on matter fields from the Quijote suite, where they are shown to reach state-of-the-art parameter constraints on the base $Λ$CDM parameters, as well as the sum of neutrino masses. We show that we can improve constraints by a factor 5 to 10 in all parameters with respect to the power spectrum baseline.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

NECOLA: Towards a Universal Field-level Cosmological Emulator

We train convolutional neural networks to correct the output of fast and approximate N-body simulations at the field level. Our model, Neural Enhanced COLA --NECOLA--, takes as input a snapshot generated by the computationally efficient COLA code and corrects the positions of the cold dark matter particles to match the results of full N-body Quijote simulations. We quantify the accuracy of the network using several summary statistics, and find that NECOLA can reproduce the results of the full N-body simulations with sub-percent accuracy down to $k\simeq1~h{\rm Mpc}^{-1}$. Furthermore, the model, that was trained on simulations with a fixed value of the cosmological parameters, is also able to correct the output of COLA simulations with different values of $Ω_{\rm m}$, $Ω_{\rm b}$, $h$, $n_s$, $σ_8$, $w$, and $M_ν$ with very high accuracy: the power spectrum and the cross-correlation coefficients are within $\simeq1\%$ down to $k=1~h{\rm Mpc}^{-1}$. Our results indicate that the correction to the power spectrum from fast/approximate simulations or field-level perturbation theory is rather universal. Our model represents a first step towards the development of a fast field-level emulator to sample not only primordial mode amplitudes and phases, but also the parameter space defined by the values of the cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Removing Astrophysics in 21 cm maps with Neural Networks

Measuring temperature fluctuations in the 21 cm signal from the Epoch of Reionization and the Cosmic Dawn is one of the most promising ways to study the Universe at high redshifts. Unfortunately, the 21 cm signal is affected by both cosmology and astrophysics processes in a non-trivial manner. We run a suite of 1,000 numerical simulations with different values of the main astrophysical parameters. From these simulations we produce tens of thousands of 21 cm maps at redshifts $10\leq z\leq 20$. We train a convolutional neural network to remove the effects of astrophysics from the 21 cm maps, and output maps of the underlying matter field. We show that our model is able to generate 2D matter fields that not only resemble the true ones visually, but whose statistical properties agree with the true ones within a few percent down to pretty small scales. We demonstrate that our neural network retains astrophysical information, that can be used to constrain the value of the astrophysical parameters. Finally, we use saliency maps to try to understand which features of the 21 cm maps is the network using in order to determine the value of the astrophysical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe&#39;s Fundamental Parameters with Artificial Intelligence

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The GIGANTES dataset: precision cosmology from voids in the machine learning era

We present GIGANTES, the most extensive and realistic void catalog suite ever released -- containing over 1 billion cosmic voids covering a volume larger than the observable Universe, more than 20 TB of data, and created by running the void finder VIDE on QUIJOTE&#39;s halo simulations. The expansive and detailed GIGANTES suite, spanning thousands of cosmological models, opens up the study of voids, answering compelling questions: Do voids carry unique cosmological information? How is this information correlated with galaxy information? Leveraging the large number of voids in the GIGANTES suite, our Fisher constraints demonstrate voids contain additional information, critically tightening constraints on cosmological parameters. We use traditional void summary statistics (void size function, void density profile) and the void auto-correlation function, which independently yields an error of $0.13\,\mathrm{eV}$ on $\sum\,m_ν$ for a 1 $h^{-3}\mathrm{Gpc}^3$ simulation, without CMB priors. Combining halos and voids we forecast an error of $0.09\,\mathrm{eV}$ from the same volume. Extrapolating to next generation multi-Gpc$^3$ surveys such as DESI, Euclid, SPHEREx, and the Roman Space Telescope, we expect voids should yield an independent determination of neutrino mass. Crucially, GIGANTES is the first void catalog suite expressly built for intensive machine learning exploration. We illustrate this by training a neural network to perform likelihood-free inference on the void size function. Cosmology problems provide an impetus to develop novel deep learning techniques, leveraging the symmetries embedded throughout the universe from physical laws, interpreting models, and accurately predicting errors. With GIGANTES, machine learning gains an impressive dataset, offering unique problems that will stimulate new techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining $M_ν$ with the Bispectrum I: Breaking Parameter Degeneracies

Massive neutrinos suppress the growth of structure below their free-streaming scale and leave an imprint on large-scale structure. Measuring this imprint allows us to constrain the sum of neutrino masses, $M_ν$, a key parameter in particle physics beyond the Standard Model. However, degeneracies among cosmological parameters, especially between $M_ν$ and $σ_8$, limit the constraining power of standard two-point clustering statistics. In this work, we investigate whether we can break these degeneracies and constrain $\smnu$ with the next higher-order correlation function --- the bispectrum. We first examine the redshift-space halo bispectrum of $800$ $N$-body simulations from the HADES suite and demonstrate that the bispectrum helps break the $M_ν$--$σ_8$ degeneracy. Then using 22,000 $N$-body simulations of the Quijote suite, we quantify for the first time the full information content of the redshift-space halo bispectrum down to nonlinear scales using a Fisher matrix forecast of $\{Ω_m$, $Ω_b$, $h$, $n_s$, $σ_8$, $M_ν\}$. For $k_{\rm max}{=}0.5~h/{\rm Mpc}$, the bispectrum provides $Ω_m$, $Ω_b$, $h$, $n_s$, and $σ_8$ constraints 1.9, 2.6, 3.1, 3.6, and 2.6 times tighter than the power spectrum. For $M_ν$, the bispectrum improves the 1$σ$ constraint from 0.2968 to 0.0572 eV --- over 5 times tighter than the power spectrum. Even with priors from {\em Planck}, the bispectrum improves $M_ν$ constraints by a factor of 1.8. Although we reserve marginalizing over a more complete set of bias parameters to the next paper of the series, these constraints are derived for a $(1~h^{-1}{\rm Gpc})^3$ box, a substantially smaller volume than upcoming surveys. Thus, our results demonstrate that the bispectrum offers significant improvements over the power spectrum, especially for constraining $M_ν$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Super-resolution emulator of cosmological simulations using deep physical models

We present an extension of our recently developed Wasserstein optimized model to emulate accurate high-resolution features from computationally cheaper low-resolution cosmological simulations. Our deep physical modelling technique relies on restricted neural networks to perform a mapping of the distribution of the low-resolution cosmic density field to the space of the high-resolution small-scale structures. We constrain our network using a single triplet of high-resolution initial conditions and the corresponding low- and high-resolution evolved dark matter simulations from the Quijote suite of simulations. We exploit the information content of the high-resolution initial conditions as a well constructed prior distribution from which the network emulates the small-scale structures. Once fitted, our physical model yields emulated high-resolution simulations at low computational cost, while also providing some insights about how the large-scale modes affect the small-scale structure in real space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Effective Halo Model: Creating a Physical and Accurate Model of the Matter Power Spectrum and Cluster Counts

We introduce a physically-motivated model of the matter power spectrum, based on the halo model and perturbation theory. This model achieves 1\% accuracy on all $k-$scales between $k=0.02h\,\mathrm{Mpc}^{-1}$ to $k=1h\,\mathrm{Mpc}^{-1}$. Our key ansatz is that the number density of halos depends on the non-linear density contrast filtered on some unknown scale $R$. Using the Effective Field Theory of Large Scale Structure to evaluate the two-halo term, we obtain a model for the power spectrum with only two fitting parameters: $R$ and the effective `sound speed&#39;, which encapsulates small-scale physics. This is tested with two suites of cosmological simulations across a broad range of cosmologies and found to be highly accurate. Due to its physical motivation, the statistics can be easily extended beyond the power spectrum; we additionally derive the one-loop covariance matrices of cluster counts and their combination with the matter power spectrum. This yields a significantly better fit to simulations than previous models, and includes a new model for super-sample effects, which is rigorously tested with separate universe simulations. At low redshift, we find a significant ($\sim 10\%$) exclusion covariance from accounting for the finite size of halos which has not previously been modeled. Such power spectrum and covariance models will enable joint analysis of upcoming large-scale structure surveys, gravitational lensing surveys and cosmic microwave background maps on scales down to the non-linear scale. We provide a publicly released Python code.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Using the Marked Power Spectrum to Detect the Signature of Neutrinos in Large-Scale Structure

Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3σ constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6σ constraint within the same conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

BE-HaPPY: Bias Emulator for Halo Power Spectrum including massive neutrinos

We study the clustering properties of dark matter halos in real- and redshift-space in cosmologies with massless and massive neutrinos through a large set of state-of-the-art N-body simulations. We provide quick and easy-to-use prescriptions for the halo bias on linear and mildly non-linear scales, both in real and redshift space, which are valid also for massive neutrinos cosmologies. Finally we present a halo bias emulator,$\textbf{BE-HaPPY}$, calibrated on the N-body simulations, which is fast enough to be used in the standard Markov Chain Monte Carlo approach to cosmological inference. For a fiducial standard $Λ$CDM cosmology $\textbf{BE-HaPPY}$ provides percent or sub-percent accuracy on the scales of interest (linear and well into the mildly non-linear regime), meeting therefore for the halo-bias the accuracy requirements for the analysis of next-generation large--scale structure surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fisher for complements: Extracting cosmology and neutrino mass from the counts-in-cells PDF

We comprehensively analyse the cosmology dependence of counts-in-cell statistics. We focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly nonlinear scales. Based on large-deviation statistics, we parametrise the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the nonlinear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass $M_ν$ and can disentangle it from the clustering amplitude $σ_8$. Using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the $νΛ$CDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first $νΛ$CDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly nonlinear scales. The matter PDF is particularly powerful for constraining the matter density $Ω_m$, clustering amplitude $σ_8$ and the total neutrino mass $M_ν$. Adding the mildly nonlinear matter PDF to the mildly nonlinear matter power spectrum improves constraints on $Ω_m$ by a factor of 5 and $σ_8$ by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc/h)$^3$. We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Weighing neutrinos with the halo environment

Nonlinear objects like halos and voids exhibit a scale-dependent bias on linear scales in massive neutrino cosmologies. The shape of this scale-dependent bias is a unique signature of the neutrino masses, but the amplitude of the signal is generally small, of the order of $f_ν$, the contribution of neutrinos to the total matter content ($\lesssim 1\%$). In this paper, we demonstrate for the first time how the strength of this signal can be substantially enhanced by using information about the halo environment at a range of scales. This enhancement is achieved by using certain combinations of the large scale Cold Dark Matter and total matter environments of halos, both of which are measurable from galaxy clustering and weak lensing surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Beyond $Λ$CDM: Problems, solutions, and the road ahead

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $Λ$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $Λ$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Degeneracies II: Structure formation in joint simulations of Warm Dark Matter and $f(R)$ gravity

We present for the first time the outcomes of a cosmological N-body simulation that simultaneously implements a Warm Dark Matter (WDM) particle candidate and a modified gravitational interaction in the form of $f(R)$ gravity, and compare its results with the individual effects of these two independent extensions of the standard $Λ$CDM scenario, and with the reference cosmology itself. We consider a rather extreme value of the WDM particle mass ($m_{\rm WDM}=0.4$ keV) and a single realisation of $f(R)$ gravity with $|\bar{f}_{R0}|=10^{-5}$, and we investigate the impact of these models and of their combination on a wide range of cosmological observables with the aim to identify possible observational degeneracies. In particular, we focus on the large-scale matter distribution, as well as on the statistical and structural properties of collapsed halos and cosmic voids. Differently from the case of combining $f(R)$ gravity with massive neutrinos -- previously investigated in Baldi et al. (2014) -- we find that most of the considered observables do not show any significant degeneracy due to the fact that WDM and $f(R)$ gravity are characterised by individual observational footprints with a very different functional dependence on cosmic scales and halo masses. In particular, this is the case for the nonlinear matter power spectrum in real space, for the halo and sub-halo mass functions, for the halo density profiles and for the concentration-mass relation. However, other observables -- like e.g. the halo bias -- do show some level of degeneracy between the two models, while a very strong degeneracy is observed for the nonlinear matter power spectrum in redshift space, for the density profiles of small cosmic voids -- with radius below $\approx 5$ Mpc$/h$ -- and for the voids abundance as a function of the void core density.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmology with a SKA HI intensity mapping survey

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cross-correlating 21cm intensity maps with Lyman Break Galaxies in the post-reionization era

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutral hydrogen in galaxy clusters: impact of AGN feedback and implications for intensity mapping

By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $Ω_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $Ω_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Small scales structures and neutrino masses

We review the impact of massive neutrinos on cosmological observables at the linear order. By means of N-body simulations we investigate the signatures left by neutrinos on the fully non-linear regime. We present the effects induced by massive neutrinos on the matter power spectrum, the halo mass function and on the halo-matter bias in massive neutrino cosmologies. We also investigate the clustering of cosmic neutrinos within galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of massive neutrinos on the BAO peak

We study the impact of neutrino masses on the shape and height of the BAO peak of the matter correlation function, both in real and redshift space. In order to describe the nonlinear evolution of the BAO peak we run N-body simulations and compare them with simple analytic formulae. We show that the evolution with redshift of the correlation function and its dependence on the neutrino masses is well reproduced in a simplified version of the Zel&#39;dovich approximation, in which the mode-coupling contribution to the power spectrum is neglected. While in linear theory the BAO peak decreases for increasing neutrino masses, the effect of nonlinear structure formation goes in the opposite direction, since the peak broadening by large scale flows is less effective. As a result of this combined effect, the peak decreases by $\sim 0.6 \%$ for $ \sum m_ν= 0.15$ eV and increases by $\sim1.2 \%$ for $ \sum m_ν= 0.3$ eV, with respect to a massless neutrino cosmology with equal value of the other cosmological parameters. We extend our analysis to redshift space and to halos, and confirm the agreement between simulations and the analytic formulae. We argue that all analytical approaches having the Zel&#39;dovich propagator in their lowest order approximation should give comparable performances, irrespectively to their formulation in Lagrangian or in Eulerian space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The halo model in a massive neutrino cosmology

We provide a quantitative analysis of the halo model in the context of massive neutrino cosmologies. We discuss all the ingredients necessary to model the non-linear matter and cold dark matter power spectra and compare with the results of N-body simulations that incorporate massive neutrinos. Our neutrino halo model is able to capture the non-linear behavior of matter clustering with a $\sim 20\%$ accuracy up to very non-linear scales of $k=10~h/$Mpc (which would be affected by baryon physics). The largest discrepancies arise in the range $k=0.5-1~h/$Mpc where the 1-halo and 2-halo terms are comparable and are present also in a massless neutrino cosmology. However, at scales $k<0.2~h/$Mpc our neutrino halo model agrees with the results of N-body simulations at the level of 8\% for total neutrino masses of $<0.3$ eV. We also model the neutrino non-linear density field as a sum of a linear and clustered component and predict the neutrino power spectrum and the cold dark matter-neutrino cross-power spectrum up to $k=1~h/$Mpc with $\sim$ 30\% accuracy. For masses below 0.15 eV the neutrino halo model captures the neutrino induced suppression, casted in terms of matter power ratios between massive and massless scenarios, with a 2\% agreement with the results of N-body/neutrino simulations. Finally, we provide a simple application of the halo model: the computation of the clustering of galaxies, in massless and massive neutrinos cosmologies, using a simple Halo Occupation Distribution scheme and our halo model extension.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Voids in massive neutrino cosmologies

Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent $N$-body simulations with cold dark matter and neutrino particles, while the galaxy catalogs are modeled by populating the dark matter halos in simulations via a halo occupation distribution (HOD) model to reproduce the clustering properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on the impact of massive neutrinos on the following void statistical properties: number density, ellipticities, two-point statistics, density and velocity profiles. Considering the matter density field, we find that voids in massive neutrino cosmologies are less evolved than those in the corresponding massless neutrinos case: there is a larger number of small voids and a smaller number of large ones, their profiles are less evacuated, and they present a lower wall at the edge. Moreover, the degeneracy between $σ_8$ and $Ω_ν$ is broken when looking at void properties. In terms of the galaxy density field, we find that differences among cosmologies are difficult to detect because of the small number of galaxy voids in the simulations. Differences are instead present when looking at the matter density and velocity profiles around these voids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Warm dark matter signatures on the 21cm power spectrum: Intensity mapping forecasts for SKA

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weighing neutrinos with cosmic neutral hydrogen

We investigate the signatures left by massive neutrinos on the spatial distribution of neutral hydrogen (HI) in the post-reionization era by running hydrodynamic simulations that include massive neutrinos as additional collisionless particles. We find that halos in massive/massless neutrino cosmologies host a similar amount of neutral hydrogen, although for a fixed halo mass, on average, the HI mass increases with the sum of the neutrino masses. Our results show that HI is more strongly clustered in cosmologies with massive neutrinos, while its abundance, $Ω_{\rm HI}(z)$, is lower. These effects arise mainly from the impact of massive neutrinos on cosmology: they suppress both the amplitude of the matter power spectrum on small scales and the abundance of dark matter halos. Modelling the HI distribution with hydrodynamic simulations at $z > 3$, and a simple analytic model at $z<3$, we use the Fisher matrix formalism to conservatively forecast the constraints that Phase 1 of the Square Kilometre Array (SKA) will place on the sum of neutrino masses, $M_ν\equiv Σm_ν$. We find that with 10,000 hours of interferometric observations at $3 \lesssim z \lesssim 6$ from a deep and narrow survey with SKA1-LOW, the sum of the neutrino masses can be measured with an error $σ(M_ν)\lesssim0.3$ eV (95% CL). Similar constraints can be obtained with a wide and deep SKA1-MID survey at $z \lesssim 3$, using the single-dish mode. By combining data from MID, LOW, and Planck, plus priors on cosmological parameters from a Stage IV spectroscopic galaxy survey, the sum of the neutrino masses can be determined with an error $σ(M_ν)\simeq0.06$ eV (95% CL).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A coarse grained perturbation theory for the Large Scale Structure, with cosmology and time independence in the UV

Standard cosmological perturbation theory (SPT) for the Large Scale Structure (LSS) of the Universe fails at small scales (UV) due to strong nonlinearities and to multistreaming effects. In Pietroni et al. 2011 a new framework was proposed in which the large scales (IR) are treated perturbatively while the information on the UV, mainly small scale velocity dispersion, is obtained by nonlinear methods like N-body simulations. Here we develop this approach, showing that it is possible to reproduce the fully nonlinear power spectrum (PS) by combining a simple (and fast) 1-loop computation for the IR scales and the measurement of a single, dominant, correlator from N-body simulations for the UV ones. We measure this correlator for a suite of seven different cosmologies, and we show that its inclusion in our perturbation scheme reproduces the fully non-linear PS with percent level accuracy, for wave numbers up to $k\sim 0.4\, h~{\rm Mpc^{-1}}$ down to $z=0$. We then show that, once this correlator has been measured in a given cosmology, there is no need to run a new simulation for a different cosmology in the suite. Indeed, by rescaling this correlator by a proper function computable in SPT, the reconstruction procedure works also for the other cosmologies and for all redshifts, with comparable accuracy. Finally, we clarify the relation of this approach to the Effective Field Theory methods recently proposed in the LSS context.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining Warm Dark Matter with high-$z$ supernova lensing

We propose a new method to constrain the warm dark matter (WDM) particle mass, $m_χ$, based on the counts of multiply imaged, distant supernovae (SN) produced by strong lensing by intervening cosmological matter fluctuations. The counts are very sensitive to the WDM particle mass, assumed here to be $m_χ=1, 1.5, 2$ keV. We use the analytic approach developed by Das & Ostriker to compute the probability density function of the cold dark matter (CDM) convergence ($κ$) on the lens plane; such method has been extensively tested against numerical simulations. We have extended this method generalizing it to the WDM case, after testing it against WDM $N$-body simulations. Using the observed cosmic star formation history we compute the probability for a distant SN to undergo a strong lensing event in different cosmologies. A minimum observing time of 2 yr (5 yr) is required for a future 100 square degrees survey reaching $z \approx 4$ ($z \approx 3$) to disentangle at 2$σ$ a WDM ($m_χ=1$ keV) model from the standard CDM scenario. Our method is not affected by any astrophysical uncertainty (such as baryonic physics effects), and, in principle, it does not require any particular dedicated survey strategy, as it may come as a byproduct of a future SN survey.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic Degeneracies I: Joint N-body Simulations of Modified Gravity and Massive Neutrinos

We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $Λ$CDM cosmological scenario - namely an $f(R)$ theory of Modified Gravity (MG) and a cosmological background of massive neutrinos - with the aim to investigate their possible observational degeneracies. We focus on three basic statistics of the large-scale matter distribution, more specifically the nonlinear matter power spectrum, the halo mass function, and the halo bias, for which we determine the deviation with respect to the fiducial $Λ$CDM cosmology in the context of both separate and combined simulations of $f(R)$ MG and massive neutrinos scenarios. Our results show that while these two extended models separately determine very prominent and potentially detectable features in all the three statistics, when we allow them to be simultaneously at work these features are strongly suppressed, resulting in much weaker deviations from the standard model&#39;s predictions. In particular, when an $f(R)$ gravity model with $f_{R0}=-1\times 10^{-4}$ is combined with a total neutrino mass of $Σ_{i}m_{ν_{i}}=0.4$ eV, the resulting matter power spectrum, halo mass function, and bias at z=0 are found to be consistent with the standard model&#39;s predictions at the 10%, 20%, and 5% accuracy levels, respectively. Therefore, our results imply an intrinsic theoretical limit to the effective discriminating power of present and future observational data sets with respect to these widely considered extensions of the standard cosmological scenario in the absence of independent measurements of the neutrino masses from laboratory experiments, even though the high-redshift evolution might still allow to partially break the degeneracy [Abridged].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling the neutral hydrogen distribution in the post-reionization Universe: intensity mapping

We model the distribution of neutral hydrogen (HI) in the post-reionization era and investigate its detectability in 21 cm intensity mapping with the future SKA radio telescope. We rely on high resolution hydrodynamical N-body simulations. The HI is assigned a-posteriori to the gas particles following two different approaches: a halo-based method in which HI is assigned only to gas particles residing within dark matter halos; a particle-based method that assigns HI to all gas particles using a prescription based on the physical properties of the particles. The HI statistical properties are then compared to the observational properties of Damped Lyman-$α$ Absorbers (DLAs) and of lower column density systems and reasonable good agreement is found for all the cases. Among the halo-based method, we further consider two different schemes that aim at reproducing the observed properties of DLAs by distributing HI inside halos: one of this results in a much higher bias for DLAs, in agreement with recent observations, which boosts the 21 cm power spectrum by a factor $\sim 4$ with respect to the other recipe. We compute the 21 cm power spectrum from the simulated HI distribution and calculate the expected signal for both SKA1-mid and SKA1-low configurations at $2.4 \leq z \leq 4$. We find that SKA will be able to detect the 21 cm power spectrum, in the non-linear regime, up to $k\sim 1\,h$/Mpc for SKA1-mid and $k\sim 5\,h$/Mpc for SKA1-low with 100 hours of observations. We also investigate the perspective of imaging the HI distribution. Our findings indicate that SKA1-low could detect the most massive HI peaks with a signal to noise ratio (SNR) higher than 5 for an observation time of about 1000 hours at $z=4$, for a synthesized beam width of $2&#39;$. Detection at redshifts $z\geqslant2.4$ with SKA1-mid would instead require a much longer observation time to achieve a comparable SNR level.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Semi-Analytic Galaxy Formation in Massive Neutrino Cosmologies

The constraints on neutrino masses led to the revision of their cosmological role, since the existence of a cosmological neutrino background is a clear prediction of the standard cosmological model. In this paper, we study the impact of such background on the spatial distribution of both Dark Matter (DM) and galaxies, by coupling $N$-body numerical simulations with semi-analytic models (SAMs) of galaxy formation. Cosmological simulations including massive neutrinos predict a slower evolution of DM perturbations with respect to the $Λ$CDM runs with the same initial conditions and a suppression on the matter power spectrum on small and intermediate scales, thus impacting on the predicted properties of galaxy populations. We explicitly show that most of these deviations are driven by the different $σ_8$ predicted for cosmologies including a massive neutrino background. We conclude that independent estimates of $σ_8$ are needed, in order to unambiguously characterise the effect of this background on the growth of structures. Galaxy properties alone are a weak tracer of deviations with respect to the $Λ$CDM run, but their combination with the overall matter distribution at all scales allows to disentangle between different cosmological models. Moreover, these deviations go on opposite direction with respect to competing models like modified gravity, thus weakening any detectable cosmological signal. Given the ubiquitous presence of a neutrino background, these effects have to be taken into account in future missions aimed at constraining the properties of the &#34;Dark&#34; components of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

VIDE: The Void IDentification and Examination toolkit

We present VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations, characterizing their properties, and providing a platform for more detailed analysis. At its core, VIDE uses a substantially enhanced version of ZOBOV (Neyinck 2008) to calculate a Voronoi tessellation for estimating the density field and a performing a watershed transform to construct voids. Additionally, VIDE provides significant functionality for both pre- and post-processing: for example, vide can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE uses the watershed levels to place voids in a hierarchical tree, outputs a summary of void properties in plain ASCII, and provides a Python API to perform many analysis tasks, such as loading and manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. While centered around ZOBOV, the toolkit is designed to be as modular as possible and accommodate other void finders. VIDE has been in development for several years and has already been used to produce a wealth of results, which we summarize in this work to highlight the capabilities of the toolkit. VIDE is publicly available at http://bitbucket.org/cosmicvoids/vide public and http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmology with massive neutrinos I: towards a realistic modeling of the relation between matter, haloes and galaxies

By using a suite of large box-size N-body simulations that incorporate massive neutrinos as an extra set of particles, we investigate the impact of neutrino masses on the spatial distribution of dark matter haloes and galaxies. We compute the bias between the spatial distribution of dark matter haloes and the overall matter and cold dark matter distributions using statistical tools such as the power spectrum and the two-point correlation function. Overall we find a scale-dependent bias on large scales for the cosmologies with massive neutrinos. However, our results indicate that the scale-dependence in the bias is reduced if the latter is computed with respect to the cold dark matter distribution only. We find that the value of the bias on large scales is reasonably well reproduced by the Tinker fitting formula once the linear cold dark matter power spectrum is used, instead of the total matter power spectrum. We investigate whether scale-dependent bias really comes from purely neutrino&#39;s effect or from nonlinear gravitational collapse of haloes. For this purpose, we address the $Ω_ν$-$σ_8$ degeneracy and find that such degeneracy is not perfect, implying that neutrinos imprint a slight scale dependence on the large-scale bias. Finally, by using a simple halo occupation distribution (HOD) model, we investigate the impact of massive neutrinos on the distribution of galaxies within dark matter haloes. We use the main galaxy sample in the Sloan Digital Sky Survey II Data Release 7 to investigate if the small-scale galaxy clustering alone can be used to discriminate among different cosmological models with different neutrino masses. Our results suggest that different choices of the HOD parameters can reproduce the observational measurements relatively well, and we quantify the difference between the values of the HOD parameters between massless and massive neutrino cosmologies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmology with massive neutrinos III: the halo mass function and an application to galaxy clusters

We use a suite of N-body simulations that incorporate massive neutrinos as an extra-set of particles to investigate their effect on the halo mass function. We show that for cosmologies with massive neutrinos the mass function of dark matter haloes selected using the spherical overdensity (SO) criterion is well reproduced by the fitting formula of Tinker et al. (2008) once the cold dark matter power spectrum is considered instead of the total matter power, as it is usually done. The differences between the two implementations, i.e. using $P_{\rm cdm}(k)$ instead of $P_{\rm m}(k)$, are more pronounced for large values of the neutrino masses and in the high end of the halo mass function: in particular, the number of massive haloes is higher when $P_{\rm cdm}(k)$ is considered rather than $P_{\rm m}(k)$. As a quantitative application of our findings we consider a Planck-like SZ-clusters survey and show that the differences in predicted number counts can be as large as $30\%$ for $\sum m_ν= 0.4$ eV. Finally, we use the Planck-SZ clusters sample, with an approximate likelihood calculation, to derive Planck-like constraints on cosmological parameters. We find that, in a massive neutrino cosmology, our correction to the halo mass function produces a shift in the $σ_8(Ω_{\rm m}/0.27)^γ$ relation which can be quantified as $Δγ\sim 0.05$ and $Δγ\sim 0.14$ assuming one ($N_ν=1$) or three ($N_ν=3$) degenerate massive neutrino, respectively. The shift results in a lower mean value of $σ_8$ with $Δσ_8 = 0.01$ for $N_ν=1$ and $Δσ_8 = 0.02$ for $N_ν=3$, respectively. Such difference, in a cosmology with massive neutrinos, would increase the tension between cluster abundance and Planck CMB measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutrino Signatures on the High Transmission Regions of the Lyman-alpha Forest

We quantify the impact of massive neutrinos on the statistics of low density regions in the intergalactic medium (IGM) as probed by the Lyman-alpha forest at redshifts z=2.2--4. Based on mock but realistic quasar (QSO) spectra extracted from hydrodynamic simulations with cold dark matter, baryons and neutrinos, we find that the probability distribution of weak Lyman-alpha absorption features, as sampled by Lyman-alpha flux regions at high transmissivity, is strongly affected by the presence of massive neutrinos. We show that systematic errors affecting the Lyman-alpha forest reduce but do not erase the neutrino signal. Using the Fisher matrix formalism, we conclude that the sum of the neutrino masses can be measured, using the method proposed in this paper, with a precision smaller than 0.4 eV using a catalog of 200 high resolution (S/N~100) QSO spectra. This number reduces to 0.27 eV by making use of reasonable priors in the other parameters that also affect the statistics of the high transitivity regions of the Lyman-alpha forest. The constraints obtained with this method can be combined with independent bounds from the CMB, large scale structures and measurements of the matter power spectrum from the Lyman-alpha forest to produce tighter upper limits on the sum of the masses of the neutrinos.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-linear evolution of the cosmic neutrino background

We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations. Our set of simulations explore the properties of neutrinos in a reference $Λ$CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass $10^{11}-10^{15}$ $h^{-1}$M$_{\odot}$, over a redshift range $z=0-2$. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than $\sim 10^{13.5}h^{-1}$M$_{\odot}$. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Neutrino Halos in Clusters of Galaxies and their Weak Lensing Signature

We study whether non-linear gravitational effects of relic neutrinos on the development of clustering and large-scale structure may be observable by weak gravitational lensing. We compute the density profile of relic massive neutrinos in a spherical model of a cluster of galaxies, for several neutrino mass schemes and cluster masses. Relic neutrinos add a small perturbation to the mass profile, making it more extended in the outer parts. In principle, this non-linear neutrino perturbation is detectable in an all-sky weak lensing survey such as EUCLID by averaging the shear profile of a large fraction of the visible massive clusters in the universe, or from its signature in the general weak lensing power spectrum or its cross-spectrum with galaxies. However, correctly modeling the distribution of mass in baryons and cold dark matter and suppressing any systematic errors to the accuracy required for detecting this neutrino perturbation is severely challenging.

0 citations0 reviewsNo code linkedOpen access