No Arabic abstract
Small- and intermediate-scale galaxy clustering can be used to establish the galaxy-halo connection to study galaxy formation and evolution and to tighten constraints on cosmological parameters. With the increasing precision of galaxy clustering measurements from ongoing and forthcoming large galaxy surveys, accurate models are required to interpret the data and extract relevant information. We introduce a method based on high-resolution N-body simulations to accurately and efficiently model the galaxy two-point correlation functions (2PCFs) in projected and redshift spaces. The basic idea is to tabulate all information of haloes in the simulations necessary for computing the galaxy 2PCFs within the framework of halo occupation distribution or conditional luminosity function. It is equivalent to populating galaxies to dark matter haloes and using the mock 2PCF measurements as the model predictions. Besides the accurate 2PCF calculations, the method is also fast and therefore enables an efficient exploration of the parameter space. As an example of the method, we decompose the redshift-space galaxy 2PCF into different components based on the type of galaxy pairs and show the redshift-space distortion effect in each component. The generalizations and limitations of the method are discussed.
We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.
We investigate how a property of a galaxy correlates most tightly with a property of its host dark matter halo, using state-of-the-art hydrodynamical simulations of galaxy formation EAGLE, Illustris, and IllustrisTNG. Unlike most of the previous work, our analyses focus on all types of galaxies, including both central and satellite galaxies. We find that the stellar mass of a galaxy at the epoch of the peak circular velocity with an evolution correction gives the tightest such correlation to the peak circular velocity $V_{rm peak}$ of the galaxys underling dark matter halo. The evolution of galaxy stellar mass reduces rather than increases scatter in such a relation. We also find that one major source of scatter comes from star stripping due to the strong interactions between galaxies. Even though, we show that the size of scatter predicted by hydrodynamical simulations has a negligible impact on the clustering of dense $V_{rm peak}$-selected subhalo from simulations, which suggests that even the simplest subhalo abundance matching (SHAM), without scatter and any additional free parameter, can provide a robust prediction of galaxy clustering that can agree impressively well with the observations from the SDSS main galaxy survey.
We study how well we can reconstruct the 2-point clustering of galaxies on linear scales, as a function of mass and luminosity, using the halo occupation distribution (HOD) in several semi-analytical models (SAMs) of galaxy formation from the Millennium Simulation. We find that HOD with Friends of Friends groups can reproduce galaxy clustering better than gravitationally bound haloes. This indicates that Friends of Friends groups are more directly related to the clustering of these regions than the bound particles of the overdensities. In general we find that the reconstruction works at best to 5% accuracy: it underestimates the bias for bright galaxies. This translates to an overestimation of 50% in the halo mass when we use clustering to calibrate mass. We also found a degeneracy on the mass prediction from the clustering amplitude that affects all the masses. This effect is due to the clustering dependence on the host halo substructure, an indication of assembly bias. We show that the clustering of haloes of a given mass increases with the number of subhaloes, a result that only depends on the underlying matter distribution. As the number of galaxies increases with the number of subhaloes in SAMs, this results in a low bias for the HOD reconstruction. We expect this effect to apply to other models of galaxy formation, including the real universe, as long as the number of galaxies incresases with the number of subhaloes. We have also found that the reconstructions of galaxy bias from the HOD model fails for low mass haloes with M = 3-5x10^11 Msun/h. We find that this is because galaxy clustering is more strongly affected by assembly bias for these low masses.
The shapes of galaxies are not randomly oriented on the sky. During the galaxy formation and evolution process, environment has a strong influence, as tidal gravitational fields in the large-scale structure tend to align nearby galaxies. Additionally, events such as galaxy mergers affect the relative alignments of both the shapes and angular momenta of galaxies throughout their history. These intrinsic galaxy alignments are known to exist, but are still poorly understood. This review will offer a pedagogical introduction to the current theories that describe intrinsic galaxy alignments, including the apparent difference in intrinsic alignment between early- and late-type galaxies and the latest efforts to model them analytically. It will then describe the ongoing efforts to simulate intrinsic alignments using both N-body and hydrodynamic simulations. Due to the relative youth of this field, there is still much to be done to understand intrinsic galaxy alignments and this review summarises the current state of the field, providing a solid basis for future work.
Reliable extraction of cosmological information from clustering measurements of galaxy surveys requires estimation of the error covariance matrices of observables. The accuracy of covariance matrices is limited by our ability to generate sufficiently large number of independent mock catalogs that can describe the physics of galaxy clustering across a wide range of scales. Furthermore, galaxy mock catalogs are required to study systematics in galaxy surveys and to test analysis tools. In this investigation, we present a fast and accurate approach for generation of mock catalogs for the upcoming galaxy surveys. Our method relies on low-resolution approximate gravity solvers to simulate the large scale dark matter field, which we then populate with halos according to a flexible nonlinear and stochastic bias model. In particular, we extend the textsc{patchy} code with an efficient particle mesh algorithm to simulate the dark matter field (the textsc{FastPM} code), and with a robust MCMC method relying on the textsc{emcee} code for constraining the parameters of the bias model. Using the halos in the BigMultiDark high-resolution $N$-body simulation as a reference catalog, we demonstrate that our technique can model the bivariate probability distribution function (counts-in-cells), power spectrum, and bispectrum of halos in the reference catalog. Specifically, we show that the new ingredients permit us to reach percentage accuracy in the power spectrum up to $ksim 0.4; ,h,{rm Mpc}^{-1}$ (within 5% up to $ksim 0.6; ,h,{rm Mpc}^{-1}$) with accurate bispectra improving previous results based on Lagrangian perturbation theory.