Miguel Icaza-Lizaola

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 ₀ w _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.

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 velocileptors, PyBird, and 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 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.

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³ 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±0.009, H ₀ = 68.8±1.2 kms^-1Mpc^-1 and A _s × 10⁹ = 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.

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.

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 haloes 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.167log₁₀(M*/M_⊙), 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 haloes of DM-only simulations with galaxies, and we discuss the next steps that are required.