No Arabic abstract
We derive a simple prescription for including beyond-linear halo bias within the standard, analytical halo-model power spectrum calculation. This results in a corrective term that is added to the usual two-halo term. We measure this correction using data from $N$-body simulations and demonstrate that it can boost power in the two-halo term by a factor of $sim2$ at scales $ksim0.7,h Mpc^{-1}$, with the exact magnitude of the boost determined by the specific pair of fields in the two-point function. How this translates to the full power spectrum depends on the relative strength of the one-halo term, which can mask the importance of this correction to a greater or lesser degree, again depending on the fields. Generally we find that our correction is more important for signals that arise from lower-mass haloes. When comparing our calculation to simulated data we find that the under-prediction of power in the transition region between the two- and one-halo terms, which typically plagues halo-model calculations, is almost completely eliminated when including the full non-linear halo bias. We show improved results for the auto and cross spectra of galaxies, haloes and matter. In the specific case of matter-matter or matter-halo power we note that a large fraction of the improvement comes from the non-linear biasing between low- and high-mass haloes. We envisage our model being useful in the analytical modelling of cross correlation signals. Our non-linear bias halo-model code is available at https://github.com/alexander-mead/BNL
The interpretation of redshift surveys requires modeling the relationship between large-scale fluctuations in the observed number density of tracers, $delta_mathrm{h}$, and the underlying matter density, $delta$. Bias models often express $delta_mathrm{h}$ as a truncated series of integro-differential operators acting on $delta$, each weighted by a bias parameter. Due to the presence of `composite operators (obtained by multiplying fields evaluated at the same spatial location), the linear bias parameter measured from clustering statistics does not coincide with that appearing in the bias expansion. This issue can be cured by re-writing the expansion in terms of `renormalised operators. After providing a pedagogical and comprehensive review of bias renormalisation in perturbation theory, we generalize the concept to non-perturbative dynamics and successfully apply it to dark-matter haloes extracted from a large suite of N-body simulations. When comparing numerical and perturbative results, we highlight the effect of the window function employed to smooth the random fields. We then measure the bias parameters as a function of halo mass by fitting a non-perturbative bias model (both before and after applying renormalisation) to the cross spectrum $P_{delta_mathrm{h}delta}(k)$. Finally, we employ Bayesian model selection to determine the optimal operator set to describe $P_{delta_mathrm{h}delta}(k)$ for $k<0.2,h$ Mpc$^{-1}$ at redshift $z=0$. We find that it includes $delta, abla^2delta, delta^2$ and the square of the traceless tidal tensor, $s^2$. Considering higher-order terms (in $delta$) leads to overfitting as they cannot be precisely constrained by our data. We also notice that next-to-leading-order perturbative solutions are inaccurate for $kgtrsim 0.1,h$ Mpc$^{-1}$.
Secondary halo bias, commonly known as assembly bias, is the dependence of halo clustering on a halo property other than mass. This prediction of the Lambda-Cold Dark Matter cosmology is essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate from halo assembly histories. We show conclusively that this is incorrect for cluster-size haloes. We present an up-to-date summary of secondary halo biases of high-mass haloes due to various halo properties including concentration, spin, several proxies of assembly history, and subhalo properties. While concentration, spin, and the abundance and radial distribution of subhaloes exhibit significant secondary biases, properties that directly quantify halo assembly history do not. In fact, the entire assembly histories of haloes in pairs are nearly identical to those of isolated haloes. In general, a global correlation between two halo properties does not predict whether or not these two properties exhibit similar secondary biases. For example, assembly history and concentration (or subhalo abundance) are correlated for both paired and isolated haloes, but follow slightly different conditional distributions in these two cases. This results in a secondary halo bias due to concentration (or subhalo abundance), despite the lack of assembly bias in the strict sense for cluster-size haloes. Due to this complexity, caution must be exercised in using any one halo property as a proxy to study the secondary bias due to another property.
HI intensity mapping data traces the large-scale structure matter distribution using the integrated emission of neutral hydrogen gas (HI). The cross-correlation of the intensity maps with optical galaxy surveys can mitigate foreground and systematic effects, but has been shown to significantly depend on galaxy evolution parameters of the HI and the optical sample. Previously, we have shown that the shot noise of the cross-correlation scales with the HI content of the optical samples, such that the shot noise estimation infers the average HI masses of these samples. In this article, we present an adaptive framework for the cross-correlation of HI intensity maps with galaxy samples using our implementation of the halo model formalism (Murray et al 2018, in prep) which utilises the halo occupation distribution of galaxies to predict their power spectra. We compare two HI population models, tracing the spatial halo and the galaxy distribution respectively, and present their auto- and cross-power spectra with an associated galaxy sample. We find that the choice of the HI model and the distribution of the HI within the galaxy sample have minor significance for the shape of the auto- and cross-correlations, but highly impact the measured shot noise amplitude of the estimators, a finding we confirm with simulations. We demonstrate parameter estimation of the HI halo occupation models and advocate this framework for the interpretation of future experimental data, with the prospect of determining the HI masses of optical galaxy samples via the cross-correlation shot noise.
The strong dependence of the large-scale dark matter halo bias on the (local) non-Gaussianity parameter, f_NL, offers a promising avenue towards constraining primordial non-Gaussianity with large-scale structure surveys. In this paper, we present the first detection of the dependence of the non-Gaussian halo bias on halo formation history using N-body simulations. We also present an analytic derivation of the expected signal based on the extended Press-Schechter formalism. In excellent agreement with our analytic prediction, we find that the halo formation history-dependent contribution to the non-Gaussian halo bias (which we call non-Gaussian halo assembly bias) can be factorized in a form approximately independent of redshift and halo mass. The correction to the non-Gaussian halo bias due to the halo formation history can be as large as 100%, with a suppression of the signal for recently formed halos and enhancement for old halos. This could in principle be a problem for realistic galaxy surveys if observational selection effects were to pick galaxies occupying only recently formed halos. Current semi-analytic galaxy formation models, for example, imply an enhancement in the expected signal of ~23% and ~48% for galaxies at z=1 selected by stellar mass and star formation rate, respectively.
It has been recently shown that any halo velocity bias present in the initial conditions does not decay to unity, in agreement with predictions from peak theory. However, this is at odds with the standard formalism based on the coupled fluids approximation for the coevolution of dark matter and halos. Starting from conservation laws in phase space, we discuss why the fluid momentum conservation equation for the biased tracers needs to be modified in accordance with the change advocated in Baldauf, Desjacques & Seljak (2014). Our findings indicate that a correct description of the halo properties should properly take into account peak constraints when starting from the Vlasov-Boltzmann equation.