Constraining Galaxy-Halo Connection Using Machine Learning

TL;DR

This paper explores machine learning methods to model small-scale galaxy clustering for constraining Halo Occupation Distribution parameters, emphasizing the importance of data processing and algorithm choice for unbiased results.

Contribution

It demonstrates that with proper data handling and algorithm selection, ML can efficiently and accurately constrain HOD parameters, outperforming traditional methods in computational efficiency.

Findings

ANNs outperform RF and ridge regression in predictions

Restricting HOD prior space improves ML robustness

Combining clustering statistics enhances parameter constraints

Abstract

We investigate the potential of machine learning (ML) methods to model small-scale galaxy clustering for constraining Halo Occupation Distribution (HOD) parameters. Our analysis reveals that while many ML algorithms report good statistical fits, they often yield likelihood contours that are significantly biased in both mean values and variances relative to the true model parameters. This highlights the importance of careful data processing and algorithm selection in ML applications for galaxy clustering, as even seemingly robust methods can lead to biased results if not applied correctly. ML tools offer a promising approach to exploring the HOD parameter space with significantly reduced computational costs compared to traditional brute-force methods if their robustness is established. Using our ANN-based pipeline, we successfully recreate some standard results from recent literature. Properly restricting the HOD parameter space, transforming the training data, and carefully selecting ML algorithms are essential for achieving unbiased and robust predictions. Among the methods tested, artificial neural networks (ANNs) outperform random forests (RF) and ridge regression in predicting clustering statistics, when the HOD prior space is appropriately restricted. We demonstrate these findings using the projected two-point correlation function ($w_\mathrm{p}(r_\mathrm{p})$), angular multipoles of the correlation function ($\xi_\ell(r)$), and the void probability function (VPF) of Luminous Red Galaxies from Dark Energy Spectroscopic Instrument mocks. Our results show that while combining $w_\mathrm{p}(r_\mathrm{p})$ and VPF improves parameter constraints, adding the multipoles $\xi_0$, $\xi_2$, and $\xi_4$ to $w_\mathrm{p}(r_\mathrm{p})$ does not significantly improve the constraints.