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.
more »
« less
Galaxies and haloes on graph neural networks: Deep generative modelling scalar and vector quantities for intrinsic alignment
ABSTRACT In order to prepare for the upcoming wide-field cosmological surveys, large simulations of the Universe with realistic galaxy populations are required. In particular, the tendency of galaxies to naturally align towards overdensities, an effect called intrinsic alignments (IA), can be a major source of systematics in the weak lensing analysis. As the details of galaxy formation and evolution relevant to IA cannot be simulated in practice on such volumes, we propose as an alternative a Deep Generative Model. This model is trained on the IllustrisTNG-100 simulation and is capable of sampling the orientations of a population of galaxies so as to recover the correct alignments. In our approach, we model the cosmic web as a set of graphs, where the graphs are constructed for each halo, and galaxy orientations as a signal on those graphs. The generative model is implemented on a Generative Adversarial Network architecture and uses specifically designed Graph-Convolutional Networks sensitive to the relative 3D positions of the vertices. Given (sub)halo masses and tidal fields, the model is able to learn and predict scalar features such as galaxy and dark matter subhalo shapes; and more importantly, vector features such as the 3D orientation of the major axis of the ellipsoid and the complex 2D ellipticities. For correlations of 3D orientations the model is in good quantitative agreement with the measured values from the simulation, except for at very small and transition scales. For correlations of 2D ellipticities, the model is in good quantitative agreement with the measured values from the simulation on all scales. Additionally, the model is able to capture the dependence of IA on mass, morphological type, and central/satellite type.
more »
« less
- Award ID(s):
- 2020295
- PAR ID:
- 10371079
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 516
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 2406-2419
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Halotools, originally published in 2017, is a Python package for cosmology and astrophysics designed to generate mock universes using existing catalogs of dark matter halos (Hearin et al., 2017). A theoretical basis of the library is the so-called halo model, which describes the matter distribution of dark matter as gravitationally self-bound clouds of dark matter particles that we call halos. Halotools is designed to take an underlying catalog of dark matter halos and populate them with galaxies using subhalo abundance, or halo occupation distribution (HOD) models, creating catalogs of simulated galaxies for use in research. This release (v0.9) adds functionality to align galaxies, injecting what are known as intrinsic alignments (IA) into these catalogs. As a result, these simulated galaxy catalogs can now be created with realistically complex correlations between galaxies, mimicking some effects seen in more expensive hydrodynamic simulations.more » « less
-
The intrinsic alignments (IA) of galaxies, a key contaminant in weak lensing analyses, arise from correlations in galaxy shapes driven by tidal interactions and galaxy formation processes. Accurate IA modeling is essential for robust cosmological inference, but current approaches rely on perturbative methods that break down on nonlinear scales or on expensive simulations. We introduce IAEmu, a neural network-based emulator that predicts the galaxy position-position ( ), position-orientation ( ), and orientation-orientation ( ) correlation functions and their uncertainties using mock catalogs based on the halo occupation distribution (HOD) framework. Compared to simulations, IAEmu achieves ~3% average error for and ~5% for , while capturing the stochasticity of without overfitting. The emulator provides both aleatoric and epistemic uncertainties, helping identify regions where predictions may be less reliable. We also demonstrate generalization to non-HOD alignment signals by fitting to IllustrisTNG hydrodynamical simulation data. As a fully differentiable neural network, IAEmu enables $ $ speed-ups in mapping HOD parameters to correlation functions on GPUs, compared to CPU-based simulations. This acceleration facilitates inverse modeling via gradient-based sampling, making IAEmu a powerful surrogate model for galaxy bias and IA studies with direct applications to Stage IV weak lensing surveys.more » « less
-
Abstract Line intensity mapping (LIM) is a rapidly emerging technique for constraining cosmology and galaxy formation using multi-frequency, low angular resolution maps.Many LIM applications crucially rely on cross-correlations of two line intensity maps, or of intensity maps with galaxy surveys or galaxy/CMB lensing.We present a consistent halo model to predict all these cross-correlations and enable joint analyses, in 3D redshift-space and for 2D projected maps.We extend the conditional luminosity function formalism to the multi-line case, to consistently account for correlated scatter between multiple galaxy line luminosities.This allows us to model the scale-dependent decorrelation between two line intensity maps,a key input for foreground rejection and for approaches that estimate auto-spectra from cross-spectra.This also enables LIM cross-correlations to reveal astrophysical properties of the interstellar medium inacessible with LIM auto-spectra.We expose the different sources of luminosity scatter or “line noise” in LIM, and clarify their effects on the 1-halo and galaxy shot noise terms.In particular, we show that the effective number density of halos can in some cases exceed that of galaxies, counterintuitively.Using observational and simulation input, we implement this halo model for the Hα, [Oiii], Lyman-α, CO and [Cii] lines.We encourage observers and simulators to measure galaxy luminosity correlation coefficients for pairs of lines whenever possible.Our code is publicly available at https://github.com/EmmanuelSchaan/HaloGen/tree/LIM .In a companion paper, we use this halo model formalism and codeto highlight the degeneracies between cosmology and astrophysics in LIM, and to compare the LIM observables to galaxy detection for a number of surveys.more » « less
-
null (Ed.)ABSTRACT Galaxy intrinsic alignments (IAs) have long been recognized as a significant contaminant to weak lensing-based cosmological inference. In this paper we seek to quantify the impact of a common modelling assumption in analytic descriptions of IAs: that of spherically symmetric dark matter haloes. Understanding such effects is important as the current generation of IA models are known to be limited, particularly on small scales, and building an accurate theoretical description will be essential for fully exploiting the information in future lensing data. Our analysis is based on a catalogue of 113 560 galaxies between z = 0.06 and 1.00 from massiveblack-ii, a hydrodynamical simulation of box length $$100 \, h^{-1}$$ Mpc. We find satellite anisotropy contributes at the level of $$\ge 30\!-\!40{{\ \rm per\ cent}}$$ to the small-scale alignment correlation functions. At separations larger than $$1 \, h^{-1}$$ Mpc the impact is roughly scale independent, inducing a shift in the amplitude of the IA power spectra of $$\sim 20{{\ \rm per\ cent}}$$. These conclusions are consistent across the redshift range and between the massiveblack-ii and the illustris simulations. The cosmological implications of these results are tested using a simulated likelihood analysis. Synthetic cosmic shear data are constructed with the expected characteristics (depth, area, and number density) of a future LSST-like survey. Our results suggest that modelling alignments using a halo model based upon spherical symmetry could potentially induce cosmological parameter biases at the ∼1.5σ level for S8 and w.more » « less
An official website of the United States government
