-
How much information can be extracted from galaxy clustering at the field level?
Authors:
Nhat-Minh Nguyen,
Fabian Schmidt,
Beatriz Tucci,
Martin Reinecke,
Andrija Kostić
Abstract:
We present optimal Bayesian field-level cosmological constraints from nonlinear tracers of the large-scale structure, specifically the amplitude $σ_8$ of linear matter fluctuations inferred from rest-frame simulated dark matter halos in a comoving volume of $8\,(h^{-1}\mathrm{Gpc})^3$. Our constraint on $σ_8$ is entirely due to nonlinear information, and obtained by explicitly sampling the initial…
▽ More
We present optimal Bayesian field-level cosmological constraints from nonlinear tracers of the large-scale structure, specifically the amplitude $σ_8$ of linear matter fluctuations inferred from rest-frame simulated dark matter halos in a comoving volume of $8\,(h^{-1}\mathrm{Gpc})^3$. Our constraint on $σ_8$ is entirely due to nonlinear information, and obtained by explicitly sampling the initial conditions along with bias and noise parameters via a Lagrangian EFT-based forward model, LEFTfield. The comparison with a simulation-based inference analysis employing the power spectrum and bispectrum -- likewise using the LEFTfield forward model -- shows that, when including precisely the same modes of the same data up to $k_{\mathrm{max}}= 0.10\,h\,\mathrm{Mpc}^{-1}$ ($0.12\,h\,\mathrm{Mpc}^{-1}$), the field-level approach yields a factor of 3.5 (5.2) improvement on the $σ_8$ constraint, from 20.0% to 5.7% (17.0% to 3.3%). This study provides direct insights into cosmological information encoded in galaxy clustering beyond low-order $n$-point functions.
△ Less
Submitted 2 April, 2024; v1 submitted 5 March, 2024;
originally announced March 2024.
-
No evidence for p- or d-wave dark matter annihilation from local large-scale structure
Authors:
Andrija Kostić,
Deaglan J. Bartlett,
Harry Desmond
Abstract:
If dark matter annihilates into standard model particles with a cross-section which is velocity dependent, then Local Group dwarf galaxies will not be the best place to search for the resulting gamma ray emission. A greater flux would be produced by more distant and massive halos, with larger velocity dispersions. We construct full-sky predictions for the gamma-ray emission from galaxy- and cluste…
▽ More
If dark matter annihilates into standard model particles with a cross-section which is velocity dependent, then Local Group dwarf galaxies will not be the best place to search for the resulting gamma ray emission. A greater flux would be produced by more distant and massive halos, with larger velocity dispersions. We construct full-sky predictions for the gamma-ray emission from galaxy- and cluster-mass halos within $\sim 200 \, {\mathrm{Mpc}}$ using a suite of constrained $N$-body simulations (CSiBORG) based on the Bayesian Origin Reconstruction from Galaxies algorithm. Comparing to observations from the Fermi Large Area Telescope and marginalising over reconstruction uncertainties and other astrophysical contributions to the flux, we obtain constraints on the cross-section which are two (seven) orders of magnitude tighter than those obtained from dwarf spheroidals for $p$-wave ($d$-wave) annihilation. We find no evidence for either type of annihilation from dark matter particles with masses in the range $m_χ= 2-500 \, {\mathrm{GeV}}/c^2$, for any channel. As an example, for annihilations producing bottom quarks with $m_χ= 10 \, {\mathrm{GeV}}/c^2$, we find $a_{1} < 2.4 \times 10^{-21} \, {\mathrm{cm^3 s^{-1}}}$ and $a_{2} < 3.0 \times 10^{-18} \, {\mathrm{cm^3 s^{-1}}}$ at 95% confidence, where the product of the cross-section, $σ$, and relative particle velocity, $v$, is given by $σv = a_\ell (v/c)^{2\ell}$ and $\ell=1, 2$ for $p$-, $d$-wave annihilation, respectively. Our bounds, although failing to exclude the thermal relic cross-section for velocity-dependent annihilation channels, are among the tightest to date.
△ Less
Submitted 20 April, 2023;
originally announced April 2023.
-
Consistency tests of field level inference with the EFT likelihood
Authors:
Andrija Kostić,
Nhat-Minh Nguyen,
Fabian Schmidt,
Martin Reinecke
Abstract:
Analyzing the clustering of galaxies at the field level in principle promises access to all the cosmological information available. Given this incentive, in this paper we investigate the performance of field-based forward modeling approach to galaxy clustering using the effective field theory (EFT) framework of large-scale structure (LSS). We do so by applying this formalism to a set of consistenc…
▽ More
Analyzing the clustering of galaxies at the field level in principle promises access to all the cosmological information available. Given this incentive, in this paper we investigate the performance of field-based forward modeling approach to galaxy clustering using the effective field theory (EFT) framework of large-scale structure (LSS). We do so by applying this formalism to a set of consistency and convergence tests on synthetic datasets. We explore the high-dimensional joint posterior of LSS initial conditions by combining Hamiltonian Monte Carlo sampling for the field of initial conditions, and slice sampling for cosmology and model parameters. We adopt the Lagrangian perturbation theory forward model from [1], up to second order, for the forward model of biased tracers. We specifically include model mis-specifications in our synthetic datasets within the EFT framework. We achieve this by generating synthetic data at a higher cutoff scale $Λ_0$, which controls which Fourier modes enter the EFT likelihood evaluation, than the cutoff $Λ$ used in the inference. In the presence of model mis-specifications, we find that the EFT framework still allows for robust, unbiased joint inference of a) cosmological parameters - specifically, the scaling amplitude of the initial conditions - b) the initial conditions themselves, and c) the bias and noise parameters. In addition, we show that in the purely linear case, where the posterior is analytically tractable, our samplers fully explore the posterior surface. We also demonstrate convergence in the cases of nonlinear forward models. Our findings serve as a confirmation of the EFT field-based forward model framework developed in [2-7], and as another step towards field-level cosmological analyses of real galaxy surveys.
△ Less
Submitted 11 July, 2023; v1 submitted 15 December, 2022;
originally announced December 2022.
-
Constraints on dark matter annihilation and decay from the large-scale structure of the nearby universe
Authors:
Deaglan J. Bartlett,
Andrija Kostić,
Harry Desmond,
Jens Jasche,
Guilhem Lavaux
Abstract:
Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modelling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs, and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analy…
▽ More
Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modelling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs, and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analyses discard gamma-ray data over most of the sky. Here we construct full-sky templates for gamma-ray flux from the large-scale structure within $\sim$200 Mpc by means of a suite of constrained $N$-body simulations (CSiBORG) produced using the Bayesian Origin Reconstruction from Galaxies algorithm. Marginalising over uncertainties in this reconstruction, small-scale structure, and parameters describing astrophysical contributions to the observed gamma-ray sky, we compare to observations from the Fermi Large Area Telescope to constrain dark matter annihilation cross sections and decay rates through a Markov Chain Monte Carlo analysis. We rule out the thermal relic cross section for $s$-wave annihilation for all $m_χ\lesssim 7 {\rm \, GeV}/c^2$ at 95\% confidence if the annihilation produces gluons or quarks less massive than the bottom quark. We infer a contribution to the gamma-ray sky with the same spatial distribution as dark matter decay at $3.3σ$. Although this could be due to dark matter decay via these channels with a decay rate $Γ\approx 6 \times 10^{-28} {\rm \, s^{-1}}$, we find that a power-law spectrum of index $p=-2.75^{+0.71}_{-0.46}$, likely of baryonic origin, is preferred by the data.
△ Less
Submitted 30 November, 2022; v1 submitted 25 May, 2022;
originally announced May 2022.
-
Optimal machine-driven acquisition of future cosmological data
Authors:
Andrija Kostić,
Jens Jasche,
Doogesh Kodi Ramanah,
Guilhem Lavaux
Abstract:
We present maps classifying regions of the sky according to their information gain potential as quantified by the Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculate the response of observed cosmic structures to perturbative changes in the cosmological model and chart their respective contri…
▽ More
We present maps classifying regions of the sky according to their information gain potential as quantified by the Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculate the response of observed cosmic structures to perturbative changes in the cosmological model and chart their respective contributions to the Fisher information. Our physical forward modeling machinery transcends the limitations of contemporary analyses based on statistical summaries to yield detailed characterizations of individual 3D structures. We demonstrate this using galaxy counts data and showcase the potential of our approach by studying the information gain of the Coma cluster. We find that regions in the vicinity of the filaments and cluster core, where mass accretion ensues from gravitational infall, are the most informative about our physical model of structure formation in the Universe. Hence, collecting data in those regions would be most optimal for testing our model predictions. The results presented in this work are the first of their kind and elucidate the inhomogeneous distribution of cosmological information in the Universe. This study paves a new way forward to perform efficient targeted searches for the fundamental physics of the Universe, where search strategies are progressively refined with new cosmological data sets within an active learning framework.
△ Less
Submitted 5 December, 2021; v1 submitted 1 July, 2021;
originally announced July 2021.