Mock Galaxy Catalogs Research Articles

Inferring Cosmological Parameters on SDSS via Domain-generalized Neural Networks and Light-cone Simulations

We present a proof-of-concept simulation-based inference on Ωm and σ 8 from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ Northern Galactic Cap (NGC) catalog using neural networks and domain generalization techniques without the need of summary statistics. Using rapid light-cone simulations L-picola, mock galaxy catalogs are produced that fully incorporate the observational effects. The collection of galaxies is fed as input to a point cloud-based network, Minkowski-PointNet . We also add relatively more accurate Gadget mocks to obtain robust and generalizable neural networks. By explicitly learning the representations that reduce the discrepancies between the two different data sets via the semantic alignment loss term, we show that the latent space configuration aligns into a single plane in which the two cosmological parameters form clear axes. Consequently, during inference, the SDSS BOSS LOWZ NGC catalog maps onto the plane, demonstrating effective generalization and improving prediction accuracy compared to non-generalized models. Results from the ensemble of 25 independently trained machines find Ωm = 0.339 ± 0.056 and σ 8 = 0.801 ± 0.061, inferred only from the distribution of galaxies in the light-cone slices without relying on any indirect summary statistics. A single machine that best adapts to the Gadget mocks yields a tighter prediction of Ωm = 0.282 ± 0.014 and σ 8 = 0.786 ± 0.036. We emphasize that adaptation across multiple domains can enhance the robustness of the neural networks in observational data.

Secondary halo bias through cosmic time II. Reconstructing halo properties using clustering information

When constructing mock galaxy catalogs based on suites of dark matter halo catalogs generated with approximated, calibrated, or machine-learning approaches, assigning intrinsic properties for these tracers is a step of paramount importance, given that they can shape the abundance and spatial distribution of mock galaxies and galaxy clusters. We explore the possibility of assigning properties of dark matter halos within the context of calibrated or learning approaches, explicitly using clustering information. The goal is to retrieve the correct signal of primary and secondary large-scale effective bias as a function of properties reconstructed solely based on phase-space properties of the halo distribution and dark matter density field. The algorithm reconstructs a set of halo properties (such as virial mass, maximum circular velocity, concentration, and spin) constrained to reproduce both primary and secondary (or assembly) bias. The key ingredients of the algorithm are the implementation of individually-assigned large-scale effective bias, a multi-scale approach to account for halo exclusion, and a hierarchical assignment of halo properties. The method facilitates the assignment of halo properties, aiming to replicate the large-scale effective bias, both primary and secondary. This constitutes an improvement over previous methods in the literature, especially for the high-mass end population. We have designed a strategy for reconstructing the main properties of dark matter halos obtained using calibrated or learning algorithms, such that the one- and two-point statistics (on large scales) replicate the signal from detailed $N$-body simulations. We encourage the application of this strategy (or the implementation of our algorithm) for the generation of mock catalogs of dark matter halos based on approximated methods.

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

ABSTRACT We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A bacusS ummitN-body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

Percolation Statistics in the MillenniumTNG Simulations

The statistical analysis of cosmic large-scale structure is most often based on simple two-point summary statistics, like the power spectrum or the two-point correlation function of a sample of galaxies or other types of tracers. In contrast, topological measures of clustering are also sensitive to higher-order correlations and thus offer the prospect to access additional information that may harbor important constraining power. We here revisit one such geometric measure of the cosmic web in the form of the so-called percolation analysis, using the recent MillenniumTNG simulation suite of the ΛCDM paradigm. We analyze continuum percolation statistics both for high-resolution dark matter particle distributions and for galaxy mock catalogs from a semianalytic galaxy formation model within a periodic simulation volume of 3000 Mpc on a side. For comparison, we also investigate the percolation statistics of random particle sets and neutrino distributions with two different summed particle masses. We find that the percolation statistics of the dark matter distribution evolves strongly with redshift and thus clustering strength, yielding a progressively lower percolation threshold toward later times. However, there is a sizable residual dependence on numerical resolution, which we interpret as a residual influence of different levels of shot noise. This is corroborated by our analysis of galaxy mock catalogs, whose results depend on sampling density more strongly than on galaxy selection criteria. While this limits the discriminative power of percolation statistics, our results suggest that it still remains useful as a complementary cosmological test when controlled for sampling density.

MAMBO: An empirical galaxy and AGN mock catalogue for the exploitation of future surveys

Current and future large surveys will produce unprecedented amounts of data. Realistic simulations have become essential for the design and development of these surveys, as well as for the interpretation of the results. We present MAMBO, a flexible and efficient workflow to build empirical galaxy and active galactic nucleus (AGN) mock catalogues that reproduce the physical and observational properties of these sources. We started with simulated dark matter (DM) haloes, to preserve the link with the cosmic web, and we populated them with galaxies and AGN using abundance matching techniques. We followed an empirical methodology, using stellar mass functions, host galaxy AGN mass functions, and AGN accretion rate distribution functions studied at different redshifts to assign, among other properties, stellar masses, the fraction of quenched galaxies, or the AGN activity (demography, obscuration, multiwavelength emission, etc.). As a proof test, we applied the method to a Millennium DM lightcone of 3.14 $ deg^2$ up to a redshift of $z=10$ and down to stellar masses of $ M \, M_ We show that the AGN population from the mock lightcone presented here reproduces with good accuracy various observables, such as state-of-the-art luminosity functions in the X-ray up to $z 7$ and in the ultraviolet up to $z 5$, optical/near-infrared colour-colour diagrams, and narrow emission line diagnostic diagrams. Finally, we demonstrate how this catalogue can be used to make useful predictions for large surveys. Using Euclid as a case example, we compute, among other forecasts, the expected surface densities of galaxies and AGN detectable in the Euclid $H_ E $ band. We find that Euclid might observe (on $H_ E $ only) about $10^ $ and $8 $ type 1 and 2 AGN, respectively, and $2 $ galaxies at the end of its $ deg^2$ Wide survey, in good agreement with other published forecasts.

The cosmic microwave background lensing imprint of cosmic voids detected in the WISE-Pan-STARRS luminous red galaxy catalog

The cross-correlation of cosmic voids with the lensing convergence (kappa ) map of cosmic microwave background (CMB) fluctuations offers a powerful tool with which to refine our understanding of the dark sector in the consensus cosmological model. Our principal aim is to compare the lensing signature of our galaxy dataset with simulations based on the concordance model and to characterize the results with an $A_ kappa $ consistency parameter normalized to unity. In particular, our measurements contribute to the understanding of the “lensing-is-low” tension of the Lambda CDM model. In this analysis, we selected luminous red galaxies (LRGs) from the WISE- Pan-STARRS dataset, enabling an extended cross-correlation measurement using a 14,200 deg$^2$ sky area, which offers a more precise measurement than previous studies. We created 2D and 3D void catalogs to cross-correlate their locations with the Planck CMB lensing map and studied their average imprint signal using a stacking methodology. Applying the same procedure, we also generated a mock galaxy catalog from the WebSky simulation to serve as a basis for comparison. The 2D void analysis reveals a good agreement with the standard cosmological model, with $A_ kappa 0.08$ amplitude; that is, $S/N=13.3$, showing a higher signal-to-noise than previous studies using voids detected in the Dark Energy Survey (DES) dataset. The 3D void analysis exhibited a lower signal-to-noise ratio and demonstrated worse agreement with our mock catalog than the 2D voids. These deviations might be attributed to limitations in the mock catalog, such as imperfections in the LRG selection, as well as a potential asymmetry between the northern and southern patches of the Pan-STARRS dataset in terms of data quality. Overall, we present a significant detection of a CMB lensing signal associated with cosmic voids, largely consistent with the concordance model. Future analyses using even larger datasets also hold great promise of further sharpening these results, given their complementary nature to large-scale structure analyses.

Generating mock galaxy catalogues for flux-limited samples like the DESI Bright Galaxy Survey

ABSTRACT Accurate mock galaxy catalogues are crucial to validate analysis pipelines used to constrain dark energy models. We present a fast HOD-fitting method which we apply to the AbacusSummit simulations to create a set of mock catalogues for the DESI Bright Galaxy Survey, which contain r-band magnitudes and $(g-r)$ colours. The halo tabulation method fits HODs for different absolute magnitude threshold samples simultaneously, preventing unphysical HOD crossing between samples. We validate the HOD fitting procedure by fitting to real-space clustering measurements and galaxy number densities from the MXXL BGS mock, which was tuned to the SDSS and GAMA surveys. The best-fitting clustering measurements and number densities are mostly within the assumed errors, but the clustering for the faint samples is low on large scales. The best-fitting HOD parameters are robust when fitting to simulations with different realizations of the initial conditions. When varying the cosmology, trends are seen as a function of each cosmological parameter. We use the best-fitting HOD parameters to create cubic box and cut sky mocks from the AbacusSummit simulations, in a range of cosmologies. As an illustration, we compare the ${}^{0.1}M_r\lt -20$ sample of galaxies in the mock with BGS measurements from the DESI one-percent survey. We find good agreement in the number densities, and the projected correlation function is reasonable, with differences that can be improved in the future by fitting directly to BGS clustering measurements. The cubic box and cut-sky mocks in different cosmologies are made publicly available.

Joint inference of multiplicative and additive systematics in galaxy density fluctuations and clustering measurements

ABSTRACT Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical uncertainties get smaller, correcting for these systematic contaminants becomes increasingly important for unbiased cosmological analysis. We present and validate a new method for understanding and mitigating both additive and multiplicative systematics in galaxy clustering measurements (two-point function) by joint inference of contaminants in the galaxy overdensity field (one-point function) using a maximum-likelihood estimator (MLE). We test this methodology with Kilo-Degree Survey-like mock galaxy catalogues and synthetic systematic template maps. We estimate the cosmological impact of such mitigation by quantifying uncertainties and possible biases in the inferred relationship between the observed and the true galaxy clustering signal. Our method robustly corrects the clustering signal to the sub-per cent level and reduces numerous additive and multiplicative systematics from $1.5 \sigma$ to less than $0.1\sigma$ for the scenarios we tested. In addition, we provide an empirical approach to identifying the functional form (additive, multiplicative, or other) by which specific systematics contaminate the galaxy number density. Even though this approach is tested and geared towards systematics contaminating the galaxy number density, the methods can be extended to systematics mitigation for other two-point correlation measurements.

SUNBIRD: a simulation-based model for full-shape density-split clustering

ABSTRACT Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-ΛCDM framework, incorporating the effects of redshift-space, Alcock–Paczynski distortions, and models of the halo–galaxy connection. Our models reach sub-per cent level accuracy down to $1 \, h^{-1}\text{Mpc}$ and are robust against different choices of galaxy–halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on ωcdm, σ8, and ns by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters.

An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations

We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the Illustris-TNG suite. We create an HOD mock with Illustris-like correlations using our method, constraining the associated IA model parameters, with the between our model’s correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation–position and orientation–orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to . We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.

Emission line predictions for mock galaxy catalogues: a new differentiable and empirical mapping from DESI

ABSTRACT We present a simple, differentiable method for predicting emission line strengths from rest-frame optical continua using an empirically determined mapping. Extensive work has been done to develop mock galaxy catalogues that include robust predictions for galaxy photometry, but reliably predicting the strengths of emission lines has remained challenging. Our new mapping is a simple neural network implemented using the JAX Python automatic differentiation library. It is trained on Dark Energy Spectroscopic Instrument Early Release data to predict the equivalent widths (EWs) of the eight brightest optical emission lines (including H α, H β, [O ii], and [O iii]) from a galaxy’s rest-frame optical continuum. The predicted EW distributions are consistent with the observed ones when noise is accounted for, and we find Spearman’s rank correlation coefficient ρs > 0.87 between predictions and observations for most lines. Using a non-linear dimensionality reduction technique, we show that this is true for galaxies across the full range of observed spectral energy distributions. In addition, we find that adding measurement uncertainties to the predicted line strengths is essential for reproducing the distribution of observed line-ratios in the BPT diagram. Our trained network can easily be incorporated into a differentiable stellar population synthesis pipeline without hindering differentiability or scalability with GPUs. A synthetic catalogue generated with such a pipeline can be used to characterize and account for biases in the spectroscopic training sets used for training and calibration of photo-z’s, improving the modelling of systematic incompleteness for the Rubin Observatory LSST and other surveys.

Cosmological constraints from density-split clustering in the BOSS CMASS galaxy sample

ABSTRACT We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1 \, h^{-1}\, \text{Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-ΛCDM framework, including redshift-space and Alcock–Paczynski distortions. Our base-ΛCDM analysis finds ωcdm = 0.1201 ± 0.0022, σ8 = 0.792 ± 0.034, and ns = 0.970 ± 0.018, corresponding to fσ8 = 0.462 ± 0.020 at z ≈ 0.525, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-ΛCDM, varying the running of the spectral index, the dark energy equation of state, and the density of mass-less relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy–halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study.

Improving and extending non-Poissonian distributions for satellite galaxies sampling in HOD: applications to eBOSS ELGs

ABSTRACT Halo occupation distribution (HOD) models help us to connect observations and theory, by assigning galaxies to dark matter haloes. In this work, we study one of the components of HOD models: the probability distribution function (PDF), which is used to assign a discrete number of galaxies to a halo, given a mean number of galaxies. For satellite galaxies, the most commonly used PDF is a Poisson Distribution. PDFs with super-Poisson variances have also been studied, allowing for continuous values of variances. This has not been the case for sub-Poisson variances, for which only the nearest integer distribution, with a single variance, has been used in the past. In this work, we propose a distribution based on the binomial one, which provides continuous sub-Poisson variances. We have generated mock galaxy catalogues from two dark-matter only simulations, unit and outerrim, with HOD models assuming different PDFs. We show that the variance of the PDF for satellite galaxies affects the one-halo term of the projected correlation function, and the count-in-cells (CIC) one point statistics. We fit the clustering of eBOSS emission line galaxies, finding results consistent with the literature, when we only vary the parameter controlling the PDF variance and the fraction of satellites. Using a mock catalogue as a reference, we have also included both the clustering and CIC to constrain the parameters of the HOD model. CIC can provide strong constraints to the PDF variance of satellite galaxies.

Leaving No Branches Behind: Predicting Baryonic Properties of Galaxies from Merger Trees

Galaxies play a key role in our endeavor to understand how structure formation proceeds in the Universe. For any precision study of cosmology or galaxy formation, there is a strong demand for huge sets of realistic mock galaxy catalogs, spanning cosmologically significant volumes. For such a daunting task, methods that can produce a direct mapping between dark matter halos from dark matter-only simulations and galaxies are strongly preferred, as producing mocks from full-fledged hydrodynamical simulations or semi-analytical models is too expensive. Here, we present a graph-neural-network-based model that is able to accurately predict key properties of galaxies such as stellar mass, g − r color, star formation rate, gas mass, stellar metallicity, and gas metallicity, purely from dark matter properties extracted from halos along the full assembly history of the galaxies. Tests based on the TNG300 simulation of the IllustrisTNG project show that our model can recover the baryonic properties of galaxies to high accuracy, over a wide redshift range (z = 0–5), for all galaxies with stellar masses more massive than 109 M ⊙ and their progenitors, with strong improvements over the state-of-the-art methods. We further show that our method makes substantial strides toward providing an understanding of the implications of the IllustrisTNG galaxy formation model.

Fast forward modelling of galaxy spatial and statistical distributions

A forward modelling approach provides simple, fast and realistic simulations of galaxy surveys, without a complex underlying model. For this purpose, galaxy clustering needs to be simulated accurately, both for the usage of clustering as its own probe and to control systematics. We present a forward model to simulate galaxy surveys, where we extend the Ultra-Fast Image Generator to include galaxy clustering. We use the distribution functions of the galaxy properties, derived from a forward model adjusted to observations. This population model jointly describes the luminosity functions, sizes, ellipticities, SEDs and apparent magnitudes. To simulate the positions of galaxies, we then use a two-parameter relation between galaxies and halos with Subhalo Abundance Matching (SHAM). We simulate the halos and subhalos using the fast PINOCCHIO code, and a method to extract the surviving subhalos from the merger history. Our simulations contain a red and a blue galaxy population, for which we build a SHAM model based on star formation quenching. For central galaxies, mass quenching is controlled with the parameter Mlimit, with blue galaxies residing in smaller halos. For satellite galaxies, environmental quenching is implemented with the parameter tquench, where blue galaxies occupy only recently merged subhalos. We build and test our model by comparing to imaging data from the Dark Energy Survey Year 1. To ensure completeness in our simulations, we consider the brightest galaxies with i < 20. We find statistical agreement between our simulations and the data for two-point correlation functions on medium to large scales. Our model provides constraints on the two SHAM parameters Mlimit and tquench and offers great prospects for the quick generation of galaxy mock catalogues, optimized to agree with observations.

A new test of gravity – I. Introduction to the method

ABSTRACT We introduce a new scheme based on the marked correlation function to probe gravity using the large-scale structure of the Universe. We illustrate our approach by applying it to simulations of the metric-variation f(R) modified gravity theory and general relativity (GR). The modifications to the equations in f(R) gravity lead to changes in the environment of large-scale structures that could, in principle, be used to distinguish this model from GR. Applying the Monte Carlo Markov Chain algorithm, we use the observed number density and two-point clustering to fix the halo occupation distribution (HOD) model parameters and build mock galaxy catalogues from both simulations. To generate a mark for galaxies when computing the marked correlation function we estimate the local density using a Voronoi tessellation. Our approach allows us to isolate the contribution to the uncertainty in the predicted marked correlation function that arises from the range of viable HOD model parameters, in addition to the sample variance error for a single set of HOD parameters. This is critical for assessing the discriminatory power of the method. In a companion paper, we apply our new scheme to a current large-scale structure survey.

The PAU Survey: a new constraint on galaxy formation models using the observed colour redshift relation

Abstract We use the GALFORM semi-analytical galaxy formation model implemented in the Planck Millennium N-body simulation to build a mock galaxy catalogue on an observer’s past lightcone. The mass resolution of this N-body simulation is almost an order of magnitude better than in previous simulations used for this purpose, allowing us to probe fainter galaxies and hence build a more complete mock catalogue at low redshifts. The high time cadence of the simulation outputs allows us to make improved calculations of galaxy properties and positions in the mock. We test the predictions of the mock against the Physics of the Accelerating Universe Survey, a narrow band imaging survey with highly accurate and precise photometric redshifts, which probes the galaxy population over a lookback time of 8 billion years. We compare the model against the observed number counts, redshift distribution and evolution of the observed colours and find good agreement; these statistics avoid the need for model-dependent processing of the observations. The model produces red and blue populations that have similar median colours to the observations. However, the bimodality of galaxy colours in the model is stronger than in the observations. This bimodality is reduced on including a simple model for errors in the GALFORM photometry. We examine how the model predictions for the observed galaxy colours change when perturbing key model parameters. This exercise shows that the median colours and relative abundance of red and blue galaxies provide constraints on the strength of the feedback driven by supernovae used in the model.

Mock observatory: Two thousand light-cone mock catalogues of luminous red galaxies from the Hyper Suprime-Cam Survey for the cosmological large-scale analysis

ABSTRACT Estimating a reliable covariance matrix for correlation functions of galaxies is a crucial task to obtain accurate cosmological constraints from galaxy surveys. We generate 2000 independent light-cone mock luminous red galaxy (LRGs) catalogues at 0.3 ≤ z ≤ 1.25, designed to cover CAMIRA LRGs observed by the Subaru Hyper Suprime-Cam Subaru Strategic Programme (HSC SSP). We first produce full-sky light-cone halo catalogues using a COmoving Lagrangian Acceleration (COLA) technique, and then trim them to match the footprints of the HSC SSP S20A Wide layers. The mock LRGs are subsequently populated on to the trimmed halo catalogues according to the halo occupation distribution model constrained by the observed CAMIRA LRGs. The stellar mass (M⋆) is assigned to each LRG by the subhalo abundance-matching technique using the observed stellar-mass functions of CAMIRA LRGs. We evaluate photometric redshifts (photo-z) of mock LRGs by incorporating the photo-z scatter, which is derived from the observed M⋆–photo-z-scatter relations of the CAMIRA LRGs. We validate the constructed full-sky halo and light-cone LRG mock catalogues by comparing their angular clustering statistics (i.e. power spectra and correlation functions) with those measured from the halo catalogues of full N-body simulations and the CAMIRA LRG catalogues from the HSC SSP, respectively. We detect clear signatures of baryon acoustic oscillations (BAOs) from our mock LRGs, whose angular scales are well consistent with theoretical predictions. These results demonstrate that our mock LRGs can be used to evaluate covariance matrices at large scales and provide predictions for the BAO detectability and cosmological constraints.

Differentiable stochastic halo occupation distribution

ABSTRACT In this work, we demonstrate how differentiable stochastic sampling techniques developed in the context of deep reinforcement learning can be used to perform efficient parameter inference over stochastic, simulation-based, forward models. As a particular example, we focus on the problem of estimating parameters of halo occupation distribution (HOD) models that are used to connect galaxies with their dark matter haloes. Using a combination of continuous relaxation and gradient re-parametrization techniques, we can obtain well-defined gradients with respect to HOD parameters through discrete galaxy catalogue realizations. Having access to these gradients allows us to leverage efficient sampling schemes, such as Hamiltonian Monte Carlo, and greatly speed up parameter inference. We demonstrate our technique on a mock galaxy catalogue generated from the Bolshoi simulation using a standard HOD model and find near-identical posteriors as standard Markov chain Monte Carlo techniques with an increase of ∼8× in convergence efficiency. Our differentiable HOD model also has broad applications in full forward model approaches to cosmic structure and cosmological analysis.

Fast Generation of Mock Galaxy Catalogs with COLA

We investigate the feasibility of using the comoving Lagrangian acceleration (COLA) technique to efficiently generate galaxy mock catalogs that can accurately reproduce the statistical properties of observed galaxies. Our proposed scheme combines the subhalo abundance-matching (SHAM) procedure with COLA simulations, using only three free parameters: the scatter magnitude (σ scat) in SHAM, the initial redshift (z init) of the COLA simulation, and the time stride (da) used by COLA. In this proof-of-concept study, we focus on a subset of BOSS CMASS NGC galaxies within the redshift range z ∈ [0.45, 0.55]. We perform GADGET simulation and low-resolution COLA simulations with various combinations of (z init, da), each using 10243 particles in an 800 h −1 Mpc box. By minimizing the difference between COLA mock and CMASS NGC galaxies for the monopole of the two-point correlation function (2PCF), we obtain the optimal σ scat. We have found that by setting z init = 29 and da = 1/30, we achieve a good agreement between COLA mock and CMASS NGC galaxies within the range of 4–20 h −1 Mpc, with a computational cost lower by 2 orders of magnitude than that of the GADGET N-body code. Moreover, a detailed verification is performed by comparing various statistical properties, such as anisotropic 2PCF, three-point clustering, and power spectrum multipoles, which shows a similar performance of the GADGET mock and COLA mock catalogs with the CMASS NGC galaxies. Furthermore, we assess the robustness of the COLA mock catalogs for different cosmological models, demonstrating consistent results in the resulting 2PCFs. Our findings suggest that COLA simulations are a promising tool for efficiently generating mock catalogs for emulators and machine-learning analyses to explore the large-scale structure of the Universe.