-
Gaussian Lagrangian Galaxy Bias
Authors:
Jens Stücker,
Marcos Pellejero-Ibáñez,
Rodrigo Voivodic,
Raul E. Angulo
Abstract:
Understanding $\textit{galaxy bias}$ -- that is the statistical relation between matter and galaxies -- is of key importance for extracting cosmological information from galaxy surveys. While the bias function $f$ -- that is the probability of forming galaxy in a region with a given density field -- is usually approximated through a parametric expansion, we show here, that it can also be measured…
▽ More
Understanding $\textit{galaxy bias}$ -- that is the statistical relation between matter and galaxies -- is of key importance for extracting cosmological information from galaxy surveys. While the bias function $f$ -- that is the probability of forming galaxy in a region with a given density field -- is usually approximated through a parametric expansion, we show here, that it can also be measured directly from simulations in a non-parameteric way. Our measurements show that the Lagrangian bias function is very close to a Gaussian for halo selections of any mass. Therefore, we newly introduce a Gaussian bias model with several intriguing properties: (1) It predicts only strictly positive probabilities $f > 0$ (unlike expansion models), (2) It has a simple analytic renormalized form and (3) It behaves gracefully in many scenarios where the classical expansion converges poorly. We show that the Gaussian bias model describes the galaxy environment distribution $p(δ| \mathrm{g})$, the scale dependent bias function $f$ and the renormalized bias function $F$ of haloes and galaxies generally equally well or significantly better than a second order expansion with the same number of parameters. We suggest that a Gaussian bias approach may enhance the range of validity of bias schemes where the canonical expansion converges poorly and further, that it may make new applications possible, since it guarantees the positivity of predicted galaxy densities.
△ Less
Submitted 19 June, 2024; v1 submitted 3 May, 2024;
originally announced May 2024.
-
Probabilistic Lagrangian bias estimators and the cumulant bias expansion
Authors:
Jens Stücker,
Marcos Pellejero-Ibáñez,
Raul E. Angulo,
Francisco Maion,
Rodrigo Voivodic
Abstract:
The spatial distribution of galaxies is a highly complex phenomenon currently impossible to predict deterministically. However, by using a statistical $\textit{bias}$ relation, it becomes possible to robustly model the average abundance of galaxies as a function of the underlying matter density field. Understanding the properties and parametric description of the bias relation is key to extract co…
▽ More
The spatial distribution of galaxies is a highly complex phenomenon currently impossible to predict deterministically. However, by using a statistical $\textit{bias}$ relation, it becomes possible to robustly model the average abundance of galaxies as a function of the underlying matter density field. Understanding the properties and parametric description of the bias relation is key to extract cosmological information from future galaxy surveys. Here, we contribute to this topic primarily in two ways: (1) We develop a new set of probabilistic estimators for bias parameters using the moments of the Lagrangian galaxy environment distribution. These estimators include spatial corrections at different orders to measure bias parameters independently of the damping scale. We report robust measurements of a variety of bias parameters for haloes, including the tidal bias and its dependence with spin at a fixed mass. (2) We propose an alternative formulation of the bias expansion in terms of "cumulant bias parameters" that describe the response of the logarithmic galaxy density to large-scale perturbations. We find that cumulant biases of haloes are consistent with zero at orders $n > 2$. This suggests that: (i) previously reported bias relations at order $n > 2$ are an artefact of the entangled basis of the canonical bias expansion; (ii) the convergence of the bias expansion may be improved by phrasing it in terms of cumulants; (iii) the bias function is very well approximated by a Gaussian -- an avenue which we explore in a companion paper.
△ Less
Submitted 19 June, 2024; v1 submitted 3 May, 2024;
originally announced May 2024.
-
Characterizing Structure Formation through Instance Segmentation
Authors:
Daniel López-Cano,
Jens Stücker,
Marcos Pellejero Ibañez,
Raúl E. Angulo,
Daniel Franco-Barranco
Abstract:
Dark matter haloes form from small perturbations to the almost homogeneous density field of the early universe. Although it is known how large these initial perturbations must be to form haloes, it is rather poorly understood how to predict which particles will end up belonging to which halo. However, it is this process that determines the Lagrangian shape of protohaloes and is therefore essential…
▽ More
Dark matter haloes form from small perturbations to the almost homogeneous density field of the early universe. Although it is known how large these initial perturbations must be to form haloes, it is rather poorly understood how to predict which particles will end up belonging to which halo. However, it is this process that determines the Lagrangian shape of protohaloes and is therefore essential to understand their mass, spin and formation history. Here, we present a machine-learning framework to learn how the protohalo regions of different haloes emerge from the initial density field. This involves one neural network to distinguish semantically which particles become part of any halo and a second neural network that groups these particles by halo membership into different instances. This instance segmentation is done through the Weinberger method, in which the network maps particles into a pseudo-space representation where different instances can be distinguished easily through a simple clustering algorithm. Our model reliably predicts the masses and Lagrangian shapes of haloes object-by-object, as well as summary statistics like the halo-mass function. We find that our model extracts information close to optimal by comparing it to the degree of agreement between two N-body simulations with slight differences in their initial conditions. We publish our model open-source and suggest that it can be used to inform analytical methods of structure formation by studying the effect of systematic manipulations of the initial conditions.
△ Less
Submitted 21 February, 2024; v1 submitted 20 November, 2023;
originally announced November 2023.
-
Phase-space simulations of prompt cusps: simulating the formation of the first haloes without artificial fragmentation
Authors:
Lurdes Ondaro-Mallea,
Raul E. Angulo,
Jens Stücker,
Oliver Hahn,
Simon D. M. White
Abstract:
The first generation of haloes forms from the collapse of the smallest peaks in the initial density field. $N$-body simulations of this process suggest a prompt formation of a steep power-law cusp, but these calculations are plagued by numerical artifacts which casts some doubt on this result. Here, we develop new simulation methods based on the dark matter phase-space sheet approach and present r…
▽ More
The first generation of haloes forms from the collapse of the smallest peaks in the initial density field. $N$-body simulations of this process suggest a prompt formation of a steep power-law cusp, but these calculations are plagued by numerical artifacts which casts some doubt on this result. Here, we develop new simulation methods based on the dark matter phase-space sheet approach and present results which are entirely free of artificial clumps. We find that a cusp with density $ρ\propto r^{-1.5}$ is indeed formed promptly, subsequently accreting a more extended halo and participating in the hierarchical growth of later halo generations. However, our simulations also suggest that the presence of artificial clumps just before peak collapse can significantly shallow the inner profiles of the cusps. We use $N$-body simulations with controlled amounts of small-scale power to place a conservative upper limit on the scales affected by artificial clumps. Finally, we used these results to simulate the collapse of the first generation of peaks of various types and in different cosmologies, finding prompt cusps to form in all cases. We conclude that prompt cusps are a generic feature of the collapse of peaks on the free-streaming scale of the initial density field, and their structure can safely be studied using $N$-body simulations provided care is taken to excise the region potentially affected by artificial clumps.
△ Less
Submitted 11 September, 2023;
originally announced September 2023.
-
The effect of stellar encounters on the dark matter annihilation signal from prompt cusps
Authors:
Jens Stücker,
Go Ogiya,
Simon D. M. White,
Raul E. Angulo
Abstract:
Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $ρ\propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but fo…
▽ More
Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $ρ\propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but for $>90\%$ of its annihilation luminosity in all but the densest regions, where they are tidally disrupted. Here we demonstrate that individual stellar encounters, rather than the mean galactic tide, are the dominant disruptors of prompt cusps within galaxies. Their cumulative effect is fully (though stochastically) characterised by an impulsive shock strength $B_* = 2πG\intρ_*({\bf x}(t))\, \mathrm{d}t$ where $ρ_*$, the total mass density in stars, is integrated over a cusp's entire post-formation trajectory. Stellar encounters and mean tides have only a small effect on the halo annihilation luminosity seen by distant observers, but this is not true for the Galactic halo because of the Sun's position. For a 100 GeV WIMP, Earth-mass prompt cusps are predicted, and stellar encounters suppress their mean annihilation luminosity by a factor of two already at 20 kpc, so that their annihilation emission is predicted to appear almost uniform over the sky. The Galactic Center $γ$-ray Excess is thus unaffected by cusps. If it is indeed dark matter annihilation radiation, then prompt cusps in the outer Galactic halo and beyond must account for 20-80% of the observed isotropic $γ$-ray background in the 1 to 10 GeV range.
△ Less
Submitted 25 April, 2023; v1 submitted 11 January, 2023;
originally announced January 2023.
-
Shedding light on low-mass subhalo survival and annihilation luminosity with numerical simulations
Authors:
Alejandra Aguirre-Santaella,
Miguel A. Sánchez-Conde,
Go Ogiya,
Jens Stücker,
Raul E. Angulo
Abstract:
In this work, we carry out a suite of specially-designed numerical simulations to shed further light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. Specifically, we have developed and employed an improved version of DASH, a GPU $N$-body code, to study the evolution of low-mass subhaloes inside a Milky Way-like hal…
▽ More
In this work, we carry out a suite of specially-designed numerical simulations to shed further light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. Specifically, we have developed and employed an improved version of DASH, a GPU $N$-body code, to study the evolution of low-mass subhaloes inside a Milky Way-like halo with unprecedented accuracy, reaching solar-mass and sub-parsec resolution in our simulations. We simulate subhaloes with varying mass, concentration, and orbital properties, and consider the effect of the gravitational potential of the Milky Way galaxy itself. More specifically, we analyze the evolution of both the bound mass fraction and annihilation luminosity of subhaloes, finding that most subhaloes survive until present time, even though in some cases they lose more than 99% of their mass at accretion. Baryons in the host induce a much more severe mass loss, especially when the subhalo orbit is more parallel to the galactic disk. Many of these subhaloes cross the solar galactocentric radius, thus making it easier to detect their annihilation fluxes from Earth. We find subhaloes orbiting a DM-only halo with a pericentre in the solar vicinity to lose 70-90% of their initial annihilation luminosity at redshift zero, which increases up to 99% when baryons are also included in the host. We find a strong relation between subhalo's mass loss and the effective tidal field at pericentre. Indeed, much of the dependence on concentration, orbital parameters, host potential and baryonic components can be explained through this single parameter. In addition to shedding light on the survival of low-mass galactic subhaloes, our results can provide detailed predictions that will aid current and future quests for the nature of DM.
△ Less
Submitted 18 July, 2022;
originally announced July 2022.
-
The Bacco Simulation Project: Bacco Hybrid Lagrangian Bias Expansion Model in Redshift Space
Authors:
Marcos Pellejero-Ibanez,
Raul E. Angulo,
Matteo Zennaro,
Jens Stuecker,
Sergio Contreras,
Giovanni Arico,
Francisco Maion
Abstract:
We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We b…
▽ More
We present an emulator that accurately predicts the power spectrum of galaxies in redshift space as a function of cosmological parameters. Our emulator is based on a 2nd-order Lagrangian bias expansion that is displaced to Eulerian space using cosmological $N$-body simulations. Redshift space distortions are then imprinted using the non-linear velocity field of simulated particles and haloes. We build the emulator using a forward neural network trained with the simulations of the BACCO project, which covers an 8-dimensional parameter space including massive neutrinos and dynamical dark energy. We show that our emulator provides unbiased cosmological constraints from the monopole, quadrupole, and hexadecapole of a mock galaxy catalogue that mimics the BOSS-CMASS sample down to nonlinear scales ($k\sim0.6$[$h/$Mpc]$^{3}$). This work opens up the possibility of robustly extracting cosmological information from small scales using observations of the large-scale structure of the Universe.
△ Less
Submitted 27 July, 2022; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Tidal Stripping in the Adiabatic Limit
Authors:
Jens Stücker,
Go Ogiya,
Raul E. Angulo,
Alejandra Aguirre-Santaella,
Miguel A. Sánchez-Conde
Abstract:
We present a model for the remnants of haloes that have gone through an adiabatic tidal stripping process. We show that this model exactly reproduces the remnant of an NFW halo that is exposed to a slowly increasing isotropic tidal field and approximately for an anisotropic tidal field. The model can be used to predict the asymptotic mass loss limit for orbiting subhaloes, solely as a function of…
▽ More
We present a model for the remnants of haloes that have gone through an adiabatic tidal stripping process. We show that this model exactly reproduces the remnant of an NFW halo that is exposed to a slowly increasing isotropic tidal field and approximately for an anisotropic tidal field. The model can be used to predict the asymptotic mass loss limit for orbiting subhaloes, solely as a function of the initial structure of the subhalo and the value of the tidal field at pericentre. Predictions can easily be made for differently concentrated host-haloes with and without baryonic components, which differ most notably in their relation between pericentre radius and tidal field. The model correctly predicts several empirically measured relations such as the `tidal track' and the `orbital frequency relation' that was reported by Errani & Navarro (2021) for the case of an isothermal sphere. Further, we propose applications of the `structure-tide' degeneracy which implies that increasing the concentration of a subhalo has exactly the same impact on tidal stripping as reducing the amplitude of the tidal field. Beyond this, we find that simple relations hold for the bound mass, truncation radius, WIMP annihilation luminosity and tidal ratio of tidally stripped NFW haloes in relation to quantities measured at the radius of maximum circular velocity. Finally, we note that NFW haloes cannot be completely disrupted when exposed adiabatically to tidal fields of arbitrary magnitudes. We provide an open-source implementation of our model and suggest that it can be used to improve predictions of dark matter annihilation.
△ Less
Submitted 17 March, 2023; v1 submitted 1 July, 2022;
originally announced July 2022.
-
Boundaries of chaos and determinism in the cosmos
Authors:
Mark Neyrinck,
Shy Genel,
Jens Stücker
Abstract:
According to the standard model of cosmology, the arrangement of matter in the cosmos on scales much larger than galaxies is entirely specified by the initial conditions laid down during inflation. But zooming in by dozens of orders of magnitude to microscopic (and human?) scales, quantum randomness reigns, independent of the initial conditions. Where is the boundary of determinism, and how does t…
▽ More
According to the standard model of cosmology, the arrangement of matter in the cosmos on scales much larger than galaxies is entirely specified by the initial conditions laid down during inflation. But zooming in by dozens of orders of magnitude to microscopic (and human?) scales, quantum randomness reigns, independent of the initial conditions. Where is the boundary of determinism, and how does that interplay with chaos? Here, we make a first attempt at answering this question in an astronomical context, including currently understood processes. The boundary is a function, at least, of length scale, position, and matter type (dark matter being more simply predictable). In intergalactic voids, the primordial pattern of density fluctuations is largely preserved. But we argue that within galaxies, the conditions are at minimum chaotic, and may even be influenced by non-primordial information, or randomness independent of the initial conditions. Randomness could be supplied by events such as supernovae and jets from active galactic nuclei (AGN) and other accretion disks, which, with the help of chaotic dynamics, could broadcast any possible microscopic randomness to larger scales, eventually throughout a galaxy. This may be generated or amplified by a recently investigated process called spontaneous stochasticity, or effective randomness in turbulent systems arising from arbitrarily small perturbations.
△ Less
Submitted 21 June, 2022;
originally announced June 2022.
-
Simulating the complexity of the dark matter sheet II: halo and subhalo mass functions for non-cold dark matter models
Authors:
Jens Stücker,
Raul E. Angulo,
Oliver Hahn,
Simon D. M. White
Abstract:
We present "sheet+release" simulations that reliably follow the evolution of dark matter structure at and below the dark matter free-streaming scale, where instabilities in traditional N-body simulations create a large population of spurious artificial haloes. Our simulations sample a large range of power-spectrum cutoff functions, parameterized through the half-mode scale $k_{\rm{hm}}$ and a slop…
▽ More
We present "sheet+release" simulations that reliably follow the evolution of dark matter structure at and below the dark matter free-streaming scale, where instabilities in traditional N-body simulations create a large population of spurious artificial haloes. Our simulations sample a large range of power-spectrum cutoff functions, parameterized through the half-mode scale $k_{\rm{hm}}$ and a slope parameter $β$. This parameter space can represent many non-cold dark matter models, including thermal relic warm dark matter, sterile-neutrinos, fuzzy dark matter, and a significant fraction of ETHOS models. Combining these simulations with additional N-body simulations, we find the following results. (1) Even after eliminating spurious haloes, the halo mass function in the strongly suppressed regime ($n_{\rm{X}}/n_{\rm{CDM}} < 5\%$) remains uncertain because it depends strongly on the definition of a halo. At these mass scales traditional halo finders primarily identify overdensities that are unbound, highly elongated, dominated by tidal fields, or far from virialized. (2) The regime where the suppression is smaller than a factor of 20 is quite robust to these uncertainties, however, and can be inferred reliably from suitable N-body simulations. (3) Parameterizing the suppression in the halo- and subhalo mass functions through the scales where the suppression reaches $20\%$, $50\%$ and $80\%$, we provide simple formulae which enable predictions for many non-cold dark matter models. (4) The halo mass-concentration relations in our sheet+release simulations agree well with previous results based on N-body simulations. (5) In general, we confirm the validity of previous N-body studies of warm dark matter models, largely eliminating concerns about the effects of artificial haloes.
△ Less
Submitted 19 October, 2021; v1 submitted 20 September, 2021;
originally announced September 2021.
-
Modelling galaxy clustering in redshift space with a Lagrangian bias formalism and $N$-body simulations
Authors:
Marcos Pellejero-Ibanez,
Jens Stuecker,
Raul E. Angulo,
Matteo Zennaro,
Sergio Contreras,
Giovanni Arico
Abstract:
Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describ…
▽ More
Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describe redshift-space distortions even down to small scales. This model corresponds to a second-order perturbative Lagrangian bias expansion which is advected to Eulerian space employing a displacement field extracted from $N$-body simulations. Eulerian coordinates are then transformed into redshift space by directly employing simulated velocity fields augmented with nuisance parameters capturing various possible satellite fractions and intra-halo small-scale velocities. We quantify the accuracy of our approach against samples of physically-motivated mock galaxies selected according to either Stellar Mass (SM) or Star Formation Rate (SFR) at multiple abundances and at $z=0$ and $1$. We find our model describes the monopole, quadrupole, and hexadecapole of the galaxy-power spectra down to scales of $k\approx 0.6 [h/$Mpc] within the accuracy of our simulations. This approach could pave the way to significantly increase the amount of cosmological information to be extracted from future galaxy surveys.
△ Less
Submitted 17 September, 2021;
originally announced September 2021.
-
The Boosted Potential
Authors:
Jens Stücker,
Raul E. Angulo,
Philipp Busch
Abstract:
The global gravitational potential, $φ$, is not commonly employed in the analysis of cosmological simulations, since its level sets do not show any clear correspondence to the underlying density field and its persistent structures. Here, we show that the potential becomes a locally meaningful quantity when considered from a boosted frame of reference, defined by subtracting a uniform gradient term…
▽ More
The global gravitational potential, $φ$, is not commonly employed in the analysis of cosmological simulations, since its level sets do not show any clear correspondence to the underlying density field and its persistent structures. Here, we show that the potential becomes a locally meaningful quantity when considered from a boosted frame of reference, defined by subtracting a uniform gradient term $φ_{\rm{boost}}(\boldsymbol{x}) = φ(\boldsymbol{x}) + \boldsymbol{x} \cdot \boldsymbol{a}_0$ with acceleration $\boldsymbol{a}_0$. We study this "boosted potential" in a variety of scenarios and propose several applications: (1) The boosted potential can be used to define a binding criterion that naturally incorporates the effect of tidal fields. This solves several problems of commonly-used self-potential binding checks: i) it defines a tidal boundary for each halo, ii) it is much less likely to misidentify caustics as haloes (specially in the context of warm dark matter cosmologies), and iii) performs better at identifying virialized regions of haloes -- yielding to the expected value of 2 for the virial ratio. (2) This binding check can be generalized to filaments and other cosmic structures. (3) The boosted potential facilitates the understanding of the disruption of satellite subhaloes. We propose a picture where most mass loss is explained through a lowering of the escape energy through the tidal field. (4) We discuss the possibility of understanding the topology of the potential field in a way that is independent of constant offsets in the first derivative $\boldsymbol{a}_0$. We foresee that this novel perspective on the potential can help to develop more accurate models and improve our understanding of structure formation.
△ Less
Submitted 5 October, 2021; v1 submitted 27 July, 2021;
originally announced July 2021.
-
The BACCO simulation project: biased tracers in real space
Authors:
Matteo Zennaro,
Raul E. Angulo,
Marcos Pellejero-Ibáñez,
Jens Stücker,
Sergio Contreras,
Giovanni Aricò
Abstract:
We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advecte…
▽ More
We present an emulator for the two-point clustering of biased tracers in real space. We construct this emulator using neural networks calibrated with more than $400$ cosmological models in a 8-dimensional cosmological parameter space that includes massive neutrinos an dynamical dark energy. The properties of biased tracers are described via a Lagrangian perturbative bias expansion which is advected to Eulerian space using the displacement field of numerical simulations. The cosmology-dependence is captured thanks to a cosmology-rescaling algorithm. We show that our emulator is capable of describing the power spectrum of galaxy formation simulations for a sample mimicking that of a typical Emission-Line survey at $z \sim 1$ with an accuracy of $1-2\%$ up to nonlinear scales $k \sim 0.7 h \mathrm{Mpc}^{-1}$.
△ Less
Submitted 2 November, 2023; v1 submitted 28 January, 2021;
originally announced January 2021.
-
Non-Halo Structures and their Effects on Gravitational Lensing
Authors:
T. R. G. Richardson,
J. Stücker,
R. E. Angulo,
O. Hahn
Abstract:
Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is currently possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_χ> 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be…
▽ More
Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is currently possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_χ> 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be correct as other structures, such as filaments and pancakes, exist and make up a significant fraction of the mass in the universe, ranging between 5$\%$ -- 50$\%$ depending on the dark matter model. Using novel fragmentation-free simulations of 1 and 3keV WDM cosmologies we study these "non-halo" structures and estimate their impact on flux-ratio observations. We find that these structures display sharp density gradients with short correlation lengths, and can contribute more to the lensing signal than all haloes up to the half-mode mass combined, thus reducing the differences expected among WDM models. We estimate that non-halo structures can be the dominant cause of line-of-sight flux-ratio anomalies in very warm, but already excluded, $m_x \sim 1 \rm{keV}$ scenarios. For colder cases $m_x \gtrsim 3 \rm{keV}$, we estimate that non-haloes can contribute about $5 - 10\%$ of the total flux-ratio signal.
△ Less
Submitted 10 February, 2022; v1 submitted 19 January, 2021;
originally announced January 2021.
-
The BACCO Simulation Project: Exploiting the full power of large-scale structure for cosmology
Authors:
Raul E. Angulo,
Matteo Zennaro,
Sergio Contreras,
Giovanni Aricò,
Marcos Pellejero-Ibañez,
Jens Stücker
Abstract:
We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear m…
▽ More
We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range $0<z<1.5$ and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build an emulator and compare its predictions against several widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2\% level for parameters in the minimal $Λ$ CDM model and to 3\% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco
△ Less
Submitted 6 July, 2021; v1 submitted 13 April, 2020;
originally announced April 2020.
-
Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM
Authors:
Jens Stücker,
Andreas. S. Schmidt,
Simon D. M. White,
Fabian Schmidt,
Oliver Hahn
Abstract:
We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with…
▽ More
We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value at $z=0$, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ up to the scale where we consider our simulations reliable, $k \leq 7 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.
△ Less
Submitted 13 November, 2020; v1 submitted 13 March, 2020;
originally announced March 2020.
-
Cosmological parameter estimation via iterative emulation of likelihoods
Authors:
Marcos Pellejero-Ibañez,
Raul E. Angulo,
Giovanni Aricó,
Matteo Zennaro,
Sergio Contreras,
Jens Stücker
Abstract:
The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterativ…
▽ More
The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterative Gaussian emulation of the target likelihood function. This requires a minimal number of likelihood evaluations and naturally accommodates for stochasticity in theoretical models. We apply the algorithm to estimate 9 parameters from the monopole and quadrupole of a mock power spectrum in redshift space. We obtain accurate posterior distribution functions with approximately 100 times fewer likelihood evaluations than an affine invariant MCMC, roughly independently from the dimensionality of the problem. We anticipate that our parameter estimation algorithm will accelerate the adoption of more accurate theoretical models in data analysis, enabling more comprehensive exploitation of cosmological observables.
△ Less
Submitted 23 January, 2020; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
Authors:
Jens Stücker,
Oliver Hahn,
Raul E. Angulo,
Simon D. M. White
Abstract:
At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented eithe…
▽ More
At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.
△ Less
Submitted 21 May, 2020; v1 submitted 30 August, 2019;
originally announced September 2019.
-
Cosmological N-Body Simulations with a Large-Scale Tidal Field
Authors:
Andreas S. Schmidt,
Simon D. M. White,
Fabian Schmidt,
Jens Stücker
Abstract:
In this paper we carry out anisotropic "separate universe" simulations by including a large-scale tidal field in the N-body code \textsc{gadget}4 using an anisotropic expansion factor $A_{ij}$. We use the code in a pure \textit{particle-mesh} (PM) mode to simulate the evolution of 16 realizations of an initial density field with and without a large-scale tidal field, which are then used to measure…
▽ More
In this paper we carry out anisotropic "separate universe" simulations by including a large-scale tidal field in the N-body code \textsc{gadget}4 using an anisotropic expansion factor $A_{ij}$. We use the code in a pure \textit{particle-mesh} (PM) mode to simulate the evolution of 16 realizations of an initial density field with and without a large-scale tidal field, which are then used to measure the \textit{response function} describing how the tidal field influences structure formation in the linear and non-linear regimes. Together with the previously measured response to a large scale overdensity, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that, contrary to the density response, the tidal response never significantly exceeds the large-scale perturbation-theory prediction even on nonlinear scales for the redshift range we discuss. We develop a simple halo model that takes into account the effect of the tidal field and compare it with our direct measurement from the anisotropic N-body simulations.
△ Less
Submitted 30 April, 2018; v1 submitted 8 March, 2018;
originally announced March 2018.
-
The median density of the Universe
Authors:
Jens Stücker,
Philipp Busch,
Simon D. M. White
Abstract:
Despite the fact that the mean matter density of the universe has been measured to an accuracy of a few percent within the standard $Λ$CDM paradigm, its median density is not known even to order of magnitude. Typical points lie in low-density regions and are not part of a collapsed structure of any scale. Locally, the dark matter distribution is then simply a stretched version of that in the early…
▽ More
Despite the fact that the mean matter density of the universe has been measured to an accuracy of a few percent within the standard $Λ$CDM paradigm, its median density is not known even to order of magnitude. Typical points lie in low-density regions and are not part of a collapsed structure of any scale. Locally, the dark matter distribution is then simply a stretched version of that in the early universe. In this single-stream regime, the distribution of unsmoothed density is sensitive to the initial power spectrum on all scales, in particular on very small scales, and hence to the nature of the dark matter. It cannot be estimated reliably using conventional cosmological simulations because of the enormous dynamic range involved, but a suitable excursion set procedure can be used instead. For the Planck cosmological parameters, a 100 GeV WIMP, corresponding to a free-streaming mass $\sim 10^{-6}M_\odot$, results in a median density of $\sim 4\times 10^{-3}$ in units of the mean density, whereas a 10 $μ$eV axion with free-streaming mass $\sim 10^{-12}M_\odot$ gives $\sim 3\times 10^{-3}$, and Warm Dark Matter with a (thermal relic) mass of 1 keV gives $\sim 8\times 10^{-2}$. In CDM (but not in WDM) universes, single-stream regions are predicted to be topologically isolated by the excursion set formalism. A test by direct N-Body simulations seems to confirm this prediction, although it is still subject to finite size and resolution effects. Unfortunately, it is unlikely that any of these properties is observable and so suitable for constraining the properties of dark matter.
△ Less
Submitted 26 October, 2017;
originally announced October 2017.