-
A comparison of effective field theory models of redshift space galaxy power spectra for DESI 2024 and future surveys
Authors:
M. Maus,
Y. Lai,
H. E. Noriega,
S. Ramirez-Solano,
A. Aviles,
S. Chen,
S. Fromenteau,
H. Gil-Marín,
C. Howlett,
M. Vargas-Magaña,
M. White,
P. Zarrouk,
J. Aguilar,
S. Ahlen,
O. Alves,
S. Brieden,
D. Brooks,
E. Burtin,
T. Claybaugh,
S. Cole,
K. Dawson,
M. Icaza-Lizaola,
A. de la Macorra,
A. de Mattia,
P. Doel
, et al. (32 additional authors not shown)
Abstract:
In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in th…
▽ More
In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in the pipelines \texttt{velocileptors}, \texttt{PyBird}, and \texttt{Folps$ν$}. While these models have been validated independently, a detailed comparison with consistent choices has not been attempted. After briefly discussing the theoretical differences between the models we describe how to provide a more apples-to-apples comparison between them. We present the results of fitting mock spectra from the \texttt{AbacusSummit} suite of N-body simulations provided in three redshift bins to mimic the types of dark time tracers targeted by the DESI survey. We show that the theories behave similarly and give consistent constraints in both the forward-modeling and ShapeFit compressed fitting approaches. We additionally generate (noiseless) synthetic data from each pipeline to be fit by the others, varying the scale cuts in order to show that the models agree within the range of scales for which we expect 1-loop perturbation theory to be applicable. This work lays the foundation of Full-Shape analysis with DESI Y1 galaxy samples where in the tests we performed, we found no systematic error associated with the modeling of the galaxy redshift space power spectrum for this volume.
△ Less
Submitted 6 June, 2024; v1 submitted 10 April, 2024;
originally announced April 2024.
-
Full Modeling and Parameter Compression Methods in configuration space for DESI 2024 and beyond
Authors:
S. Ramirez-Solano,
M. Icaza-Lizaola,
H. E. Noriega,
M. Vargas-Magaña,
S. Fromenteau,
A. Aviles,
F. Rodriguez-Martinez,
J. Aguilar,
S. Ahlen,
O. Alves,
S. Brieden,
D. Brooks,
T. Claybaugh,
S. Cole,
A. de la Macorra,
Arjun Dey,
B. Dey,
P. Doel,
K. Fanning,
J. E. Forero-Romero,
E. Gaztañaga,
H. Gil-Marín,
S. Gontcho A Gontcho,
K. Honscheid,
C. Howlett
, et al. (27 additional authors not shown)
Abstract:
In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (…
▽ More
In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full-Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying $n_s$ and allowing for $w_0w_a$CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within $1σ$ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models.
△ Less
Submitted 16 April, 2024; v1 submitted 10 April, 2024;
originally announced April 2024.
-
High precision accelerator for our hybrid model of the redshift space power spectrum
Authors:
M. Icaza-Lizaola,
Yong-Seon Song,
Minji Oh,
Yi Zheng
Abstract:
Upcoming Large Scale Structure surveys aim to achieve an unprecedented level of precision in measuring galaxy clustering. However, accurately modeling these statistics may require theoretical templates that go beyond second-order perturbation theory, especially for achieving precision at smaller scales. In our previous work, we introduced a hybrid model for the redshift space power spectrum of gal…
▽ More
Upcoming Large Scale Structure surveys aim to achieve an unprecedented level of precision in measuring galaxy clustering. However, accurately modeling these statistics may require theoretical templates that go beyond second-order perturbation theory, especially for achieving precision at smaller scales. In our previous work, we introduced a hybrid model for the redshift space power spectrum of galaxies. This model combines second-order templates with N-body simulations to capture the influence of scale-independent parameters on the galaxy power spectrum. However, the impact of scale-dependent parameters was addressed by precomputing a set of input statistics derived from computationally expensive N-body simulations. As a result, exploring the scale-dependent parameter space was not feasible in this approach. To address this challenge, we present an accelerated methodology that utilizes Gaussian processes, a machine learning technique, to emulate these input statistics. Our emulators exhibit remarkable accuracy, achieving reliable results with just 13 N-body simulations for training. We reproduce all necessary input statistics for a set of test simulations with an error of approximately 0.1 per cent in the parameter space within $5σ$ of the Planck predictions, specifically for scales around $k > 0.1$ $h$Mpc$^{-1}$. Following the training of our emulators, we can predict all inputs for our hybrid model in approximately 0.2,seconds at a specified redshift. Given that performing 13 N-body simulations is a manageable task, our present methodology enables us to construct efficient and highly accurate models of the galaxy power spectra within a manageable time frame.
△ Less
Submitted 17 November, 2023;
originally announced November 2023.
-
Full Shape Cosmology Analysis from BOSS in configuration space using Neural Network Acceleration
Authors:
Sadi Ramirez,
Miguel Icaza-Lizaola,
Sebastien Fromenteau,
Mariana Vargas-Magaña,
Alejandro Aviles
Abstract:
Recently, a new wave of full modeling analyses have emerged within the Large-Scale Structure community, leading mostly to tighter constraints on the estimation of cosmological parameters, when compared with standard approaches used over the last decade by collaboration analyses of stage III experiments. However, the majority of these full-shape analyses have primarily been conducted in Fourier spa…
▽ More
Recently, a new wave of full modeling analyses have emerged within the Large-Scale Structure community, leading mostly to tighter constraints on the estimation of cosmological parameters, when compared with standard approaches used over the last decade by collaboration analyses of stage III experiments. However, the majority of these full-shape analyses have primarily been conducted in Fourier space, with limited emphasis on exploring the configuration space. Investigating n-point correlations in configuration space demands a higher computational cost compared to Fourier space because it typically requires an additional integration step. This can pose a limitation when using these approaches, especially when considering higher-order statistics. One avenue to mitigate the high computation time is to take advantage of neural network acceleration techniques. In this work, we present a full shape analysis of Sloan Digital Sky Survey III/BOSS in configuration space using a neural network accelerator. We show that the efficacy of the pipeline is enhanced by a time factor $10^{3}$ without sacrificing precision, making it possible to reduce the error associated with the surrogate modeling to below $10^{-2}$ percent which is compatible with the precision required for current stage IV experiments such as DESI. We find $Ω_m=0.286\pm 0.009$, $H_0=68.8\pm 1.2$ $\mathrm{km} \mathrm{s^{-1}}\mathrm{Mpc^{-1}}$ and $A_s \times 10^9 =2.09 ^{+0.25}_{-0.29}$. Our results on public BOSS data are in good agreement with BOSS official results and compatible with other independent full modeling analyses. We explore relaxing the prior on $ω_b$ and varying $n_s$, without significant changes in the mean values of the cosmological parameters posterior distributions, but enlarging their widths. Finally, we explore the information content of the multipoles when constraining cosmological parameters.
△ Less
Submitted 26 October, 2023;
originally announced October 2023.
-
A sparse regression approach for populating dark matter halos and subhalos with galaxies
Authors:
M. Icaza-Lizaola,
Richard G. Bower,
Peder Norberg,
Shaun Cole,
Matthieu Schaller
Abstract:
We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM…
▽ More
We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM methods are simple and depend on a relatively small set of adjustable parameters. We collect data from 35,459 galaxies from the EAGLE simulation using 19 redshift slices between $z=0$ and $z=4$ to parameterize the mass evolution of the host halos. Using an appropriate formulation of input parameters, our methodology can model satellite and central halos using a single predictive model that achieves the same accuracy as when predicted separately. This allows us to remove the somewhat arbitrary distinction between those two galaxy types and model them based only on their halo growth history. Our models can accurately reproduce the total galaxy stellar mass function and the stellar mass-dependent galaxy correlation functions ($ξ(r)$) of EAGLE. We show that our SRM model predictions of $ξ(r)$ is competitive with those from sub-halo abundance matching and might be comparable to results from extremely randomized trees. We suggest SRM as an encouraging approach for populating the halos of dark matter only simulations with galaxies and for generating mock catalogues that can be used to explore galaxy evolution or analyse forthcoming large-scale structure surveys.
△ Less
Submitted 8 November, 2022; v1 submitted 18 May, 2022;
originally announced May 2022.
-
A sparse regression approach to modeling the relation between galaxy stellar masses and their host halos
Authors:
M. Icaza-Lizaola,
Richard G. Bower,
Peder Norberg,
Shaun Cole,
Matthieu Schaller,
Stefan Egan
Abstract:
Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 cen…
▽ More
Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 central galaxies from the EAGLE hydrodynamic simulation. By matching the host halos to a DM-only simulation, we collect the halo mass and specific angular momentum at present time and for their main progenitors in 10 redshift bins from $z=0$ to $z=4$. The principal component of our governing equation is a third-order polynomial of the host halo mass, which models the stellar-mass halo-mass relation. The scatter about this relation is driven by the halo mass evolution and is captured by second and third-order correlations of the halo mass evolution with the present halo mass. An advantage of sparse regression approaches is that unnecessary terms are removed. Although we include information on halo specific angular momentum, these parameters are discarded by our methodology. This suggests that halo angular momentum has little connection to galaxy formation efficiency. Our model has a root mean square error (RMSE) of $0.167 \log_{10}(M^*/M_\odot)$, and accurately reproduces both the stellar mass function and central galaxy correlation function of EAGLE. The methodology appears to be an encouraging approach for populating the halos of DM-only simulations with galaxies, and we discuss the next steps that are required.
△ Less
Submitted 13 August, 2021; v1 submitted 8 January, 2021;
originally announced January 2021.
-
The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 LRG sample: structure growth rate measurement from the anisotropic LRG correlation function in the redshift range 0.6 < z < 1.0
Authors:
M. Icaza-Lizaola,
M. Vargas-Magaña,
S. Fromenteau,
S. Alam,
B. Camacho,
H. Gil-Marin,
R. Paviot,
Ashley Ross,
Donald P. Schneider,
Jeremy Tinker,
Yuting Wang,
Cheng Zhao,
Abhishek Prakash,
G. Rossi,
Gong-Bo Zhao,
Irene Cruz-Gonzalez,
Axel de la Macorra
Abstract:
We analyze the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6$<z<$1.0, which consists of 80,118 galaxies from eBOSS and 46,439 galaxies from the BOSS-CMASS sample. The eB…
▽ More
We analyze the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6$<z<$1.0, which consists of 80,118 galaxies from eBOSS and 46,439 galaxies from the BOSS-CMASS sample. The eBOSS-CMASS Luminous Red Galaxy sample has a sky coverage of 1,844 deg$^2$, with an effective volume of 0.9 Gpc$^3$. The analysis was made in configuration space using a Legendre multipole expansion. The Redshift Space Distortion signal is modeled as a combination of the Convolution Lagrangian Perturbation Model and the Gaussian Streaming Model.
We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f (z_{\rm eff})σ_8(z_{\rm eff})=0.454 \pm0.139 $, and the Alcock-Paczynski dilation scales which constraints the angular diameter distance $D_A(z_{eff})=1466.5 \pm 136.6 (r_s/r_s^{\rm fid})$ and $H(z_{\rm eff})=105.8 \pm 16 (r_s^{\rm fid}/r_s) \mathrm{km\,s^{-1}\,Mpc^{-1}}$, where $r_s$ is the sound horizon at the end of the baryon drag epoch and $r_s^{\rm fid}$ is its value in the fiducial cosmology at an effective redshift $z_{\rm eff}=0.72$. These results are in full agreement with the current $Λ$-Cold Dark Matter ($Λ$-CDM) cosmological model inferred from Planck measurements. This study is the first eBOSS LRG full-shape analysis i.e. including Redshift-Space Distortions (RSD) simultaneously with the Alcock-Paczynski (AP) effect and the Baryon Acoustic Oscillation (BAO) scale.
△ Less
Submitted 22 April, 2020; v1 submitted 17 September, 2019;
originally announced September 2019.