Search

Ly$\alpha$ nebulae ubiquitously found around z>2 quasars can supply unique constraints on the properties of the Circumgalactic Medium, such as its density distribution, provided the quasar halo mass is known. We present a new method to constrain quasar halo masses based on the line-of-sight velocity dispersion maps of Ly$\alpha$ nebulae. By using MUSE-like mock observations obtained from cosmological hydrodynamic simulations under the assumption of maximal quasar fluorescence, we show that the velocity dispersion radial profiles of Ly$\alpha$-emitting gas are strongly determined by gravity and that they are thus self-similar with respect to halo mass when rescaled by the virial radius. Through simple analytical arguments and by exploiting the kinematics of HeII1640\.A emission for a set of observed nebulae, we show that Ly$\alpha$ radiative transfer effects plausibly do not change the shape of the velocity dispersion profiles but only their normalisation without breaking their self-similarity. Taking advantage of these results, we define the variable $\eta^{140-200}_{40-100}$ as the ratio of the median velocity dispersion in two specifically selected annuli and derive an analytical relation between $\eta^{140-200}_{40-100}$ and the halo mass which can be directly applied to observations. We apply our method to 37 observed quasar Ly$\alpha$ nebulae at 3

We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_\phi$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_\phi$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_\phi$ in a similar manner, roughly preserving the shape of the linear $b_\phi(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_\phi$ in opposite directions. This induces significant changes to the $b_\phi(b_1)$ relation, with higher concentration halos having higher amplitude of $b_\phi(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_\phi \approx 3$, $b_\phi \approx -1$ and $b_\phi \approx 9$, respectively. Varying the halo concentration can make $b_\phi$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_\phi(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $\sigma_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_\phi(b_1)$ relation of halos and galaxies., Comment: 21 pages, 4 figures, comments welcome; V2 minor modifications throughout, published in JCAP

We study the statistics of various large-scale structure tracers in gravity-only cosmological simulations including baryons and cold dark matter (CDM) initialized with two different transfer functions, and simulated as two distinct fluids. This allows us to study the impact of baryon-CDM relative perturbations on these statistics. In particular, we focus on the statistics of cosmic voids, as well as on the matter and halo real-space 2-point correlation function and baryon acoustic oscillations (BAO) peak. We find that the void size function is affected at the 1-2% level at maximum, and that the impact is more important at higher redshift, while the void density profile and void bias are roughly unaffected. We do not detect a sizeable impact of relative baryon-CDM perturbations on the real-space correlation functions of matter and halos or the BAO peak, which is in line with results from previous works. Our results imply that it would be hard to use voids or real-space correlation functions to constrain baryon-CDM relative perturbations, but also that we might not have to include them in models for the analysis of future cosmological surveys data., Comment: 17 pages, 12 figures, v2: version published in MNRAS

We use the forward modeling approach to galaxy clustering combined with the likelihood from the effective-field theory of large-scale structure to measure assembly bias, i.e. the dependence of halo bias on properties beyond the total mass, in the linear ($b_1$) and second order bias parameters ($b_2$ and $b_{K^2}$) of dark matter halos in $N$-body simulations. This is the first time that assembly bias in the tidal bias parameter $b_{K^2}$ is measured. We focus on three standard halo properties: the concentration $c$, spin $\lambda$, and sphericity $s$, for which we find an assembly bias signal in $b_{K^2}$ that is opposite to that in $b_1$. Specifically, at fixed mass, halos that get more (less) positively biased in $b_1$, get less (more) negatively biased in $b_{K^2}$. We also investigate the impact of assembly bias on the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations, and find that while the $b_2(b_1)$ relation stays roughly unchanged, assembly bias strongly impacts the $b_{K^2}(b_1)$ relation. This impact likely extends also to the corresponding relation for galaxies, which motivates future studies to design better priors on $b_{K^2}(b_1)$ for use in cosmological constraints from galaxy clustering data., Comment: 22 pages, 5 figures ; v2: version published in JCAP

We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear $b_1$ and second-order bias parameters $b_2$ and $b_{K^2}$. The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us to infer these bias parameters with high signal-to-noise, even from relatively small volumes ($L_{\rm box} = 205{\rm Mpc}/h$). We consider galaxies from the IllustrisTNG simulations, and our main result is that the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations obtained from gravity-only simulations for total mass selected objects are broadly preserved for simulated galaxies selected by stellar mass, star formation rate, color and black hole accretion rate. We also find good agreement between the bias relations of the simulated galaxies and a number of recent estimates for observed galaxy samples. The consistency under different galaxy selection criteria suggests that theoretical priors on these bias relations may be used to improve cosmological constraints based on observed galaxy samples. We do identify some small differences between the bias relations in the hydrodynamical and gravity-only simulations, which we show can be linked to the environmental dependence of the relation between galaxy properties and mass. We also show that the EFT likelihood recovers the value of $\sigma_8$ to percent-level from various galaxy samples (including splits by color and star formation rate) and after marginalizing over 8 bias parameters. This demonstration using simulated galaxies adds to previous works based on halos as tracers, and strengthens further the potential of forward models to infer cosmology from galaxy data., Comment: 19 pages, 9 figures. Version 2 matches published version and contains a comparison to bias estimates from observations. Comments welcomed!

Baryons and cold dark matter (CDM) did not comove prior to recombination. This leads to differences in the local baryon and CDM densities, the so-called baryon-CDM isocurvature perturbations $\delta_{bc}$. These perturbations are usually neglected in the analysis of Large-Scale Structure data but taking them into account might become important in the era of high precision cosmology. Using gravity-only 2-fluid simulations we assess the impact of such perturbations on the dark matter halos distribution. In particular, we focus on the baryon fraction in halos as a function of mass and large-scale $\delta_{bc}$, which also allows us to study details of the nontrivial numerical setup required for such simulations. We further measure the cross-power spectrum between the halo field and $\delta_{bc}$ over a wide range of mass. This cross-correlation is nonzero and negative which shows that halo formation is impacted by $\delta_{bc}$. We measure the associated bias parameter $b_{\delta_{bc}}$ and compare it to recent results, finding good agreement. Finally we quantify the impact of such perturbations on the halo-halo power spectrum and show that this effect can be degenerate with the one of massive neutrinos for surveys like DESI., Comment: 27 pages, 11 figures, Version published in JCAP

Using the publicly available Quijote simulations, we present the first measurements of the assembly bias of dark matter halos in N-body simulations which include massive neutrinos. We focus on the dependence of the linear bias $b_1$ on three halo properties: 1) concentration $c$, 2) spin $\lambda$, and 3) ellipticity $s$. Although these simulations cover a large volume, superior to any future surveys, we do not detect any effect of neutrinos on the relations $b_1(c)$, $b_1(\lambda)$ and $b_1(s)$ at fixed halo mass. We further study the dependence of halo properties and environment on neutrinos, finding these quantities to be impacted by neutrino masses at the same level as assembly bias. We find that the effect of neutrinos on spin and shape can be largely attributed to the change in the cold dark matter $\sigma_8$ in neutrinos simulations, which is not the case for concentration., Comment: 25 pages, 9 figures, v2: minor modifications through the paper. Version accepted by JCAP

We present a new simulation technique in which any chosen mode $\mathbf{k}$ of the density contrast field can be amplified by an amplitude $\Delta$. These amplified-mode simulations allow us to study the response of the halo density field to a long-wavelength mode other than the DC mode. In this sense they are a generalization of the separate-universe simulations to finite-wavelength modes. In particular, we use these simulations to obtain robust measurements of the first higher-derivative bias of dark matter halos $b_{\nabla^2 \delta}$. We find a negative bias at all mass considered, roughly following the $-R_L^2(M)$ relation, the Lagrangian radius of halos squared, as naively expected. We compare our results with those obtained from a fit to the 1-loop halo-matter power spectrum, as well as with the recent results from Abidi and Baldauf (2018), and to the prediction from the peak theory., Comment: 26 pages, 9 figures; v2: minor modifications through the paper and on fig.7. Version accepted by JCAP

We present direct measurements of cubic bias parameters of dark matter halos from the halo-matter-matter-matter trispectrum. We measure this statistic efficiently by cross-correlating the halo field measured in N-body simulations with specific third-order nonlocal transformations of the initial density field in the same simulation. Together with the recent Abidi & Baldauf (2018), these are the first measurements of halo bias using the four-point function that have been reported to date. We also obtain constraints on the quadratic bias parameters. For all individual cubic parameters involving the tidal field $\mathcal{K}_{ij}$, we find broad consistency with the prediction of the Lagrangian local-in-matter-density ansatz, with some indications of a positive Lagrangian coefficient $b_{\rm td}^L$ multiplying the time derivative of $\mathcal{K}_{ij}$. For the quadratic tidal bias ($b_{K^2}$), we obtain a significant detection of a negative Lagrangian tidal bias., Comment: 29 pages, 12 figures; v2: new renormalization procedure, added figure 12, results at higher redshifts on figs 1-5, added appendix B, C and F, clarifications throughout; v3: clarifications throughout, version accepted by JCAP

We present precise measurements of the assembly bias of dark matter halos, i.e. the dependence of halo bias on other properties than the mass, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength matter overdensity into the background density. This method measures the LIMD (local-in-matter-density) bias parameters $b_n$ in the large-scale limit. We focus on the dependence of the first two Eulerian biases $b^E_1$ and $b^E_2$ on four halo properties: the concentration, spin, mass accretion rate, and ellipticity. We quantitatively compare our results with previous works in which assembly bias was measured on fairly small scales. Despite this difference, our findings are in good agreement with previous results. We also look at the joint dependence of bias on two halo properties in addition to the mass. Finally, using the excursion set peaks model, we attempt to shed new insights on how assembly bias arises in this analytical model., Comment: 30 pages, 21 figures ; v2 : added references (sec. 1, 5), clarifications throughout ; v3 : clarifications throughout, version accepted by JCAP

We measure the large-scale bias of dark matter halos in simulations with non-Gaussian initial conditions of the local type, and compare this bias to the response of the mass function to a change in the primordial amplitude of fluctuations. The two are found to be consistent, as expected from physical arguments, for three halo-finder algorithms which use different Spherical Overdensity (SO) and Friends-of-Friends (FoF) methods. On the other hand, we find that the commonly used prediction for universal mass functions, that the scale-dependent bias is proportional to the first-order Gaussian Lagrangian bias, does not yield a good agreement with the measurements. For all halo finders, high-mass halos show a non-Gaussian bias suppressed by 10-15% relative to the universal mass function prediction. For SO halos, this deviation changes sign at low masses, where the non-Gaussian bias becomes larger than the universal prediction., Comment: 12 pages, 6 figures. v2. several improvements: new theory section, new fig.2, added references, matches version published in MNRAS

The dark matter halos that host galaxies and clusters form out of initial high-density patches, providing a biased tracer of the linear matter density field. In the simplest local bias approximation, the halo field is treated as a perturbative series in the average overdensity of the Lagrangian patch. In more realistic models, however, additional quantities will affect the clustering of halo-patches, and this expansion becomes a function of several stochastic variables. In this paper, we present a general multivariate expansion scheme that can parametrize the clustering of any biased Lagrangian tracer, given only the variables involved and their symmetry (in our case rotational invariance). This approach is based on an expansion in the orthonormal polynomials associated with the relevant variables, so that no renormalization of the coefficients ever occurs. We provide explicit expression for the series coefficients, or Lagrangian bias parameters, in the case of peaks of the linear density field. As an application of our formalism, we present a simple derivation of the original BBKS formula, and compute the non-Gaussian bias in the presence of a primordial trispectrum of the local shape., Comment: 16 pages, v2: minor modifications, version accepted by PRD

We present accurate measurements of the linear, quadratic, and cubic local bias of dark matter halos, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength overdensity. This can be seen as an exact implementation of the peak-background split argument. We compare the results with the linear and quadratic bias measured from the halo-matter power spectrum and bispectrum, and find good agreement. On the other hand, the standard peak-background split applied to the Sheth & Tormen (1999) and Tinker et al. (2008) halo mass functions matches the measured linear bias parameter only at the level of 10%. The prediction from the excursion set-peaks approach performs much better, which can be attributed to the stochastic moving barrier employed in the excursion set-peaks prediction. We also provide convenient fitting formulas for the nonlinear bias parameters $b_2(b_1)$ and $b_3(b_1)$, which work well over a range of redshifts., Comment: 23 pages, 8 figures; v2 : added references (sec. 1, 4, 5), results at higher redshifts on fig. 4 and updated fitting formulas (eqs 5.2-5.3), v3 : clarifications throughout, version accepted by JCAP