-
Gaussian Lagrangian Galaxy Bias
Authors:
Jens Stücker,
Marcos Pellejero-Ibáñez,
Rodrigo Voivodic,
Raul E. Angulo
Abstract:
Understanding $\textit{galaxy bias}$ -- that is the statistical relation between matter and galaxies -- is of key importance for extracting cosmological information from galaxy surveys. While the bias function $f$ -- that is the probability of forming galaxy in a region with a given density field -- is usually approximated through a parametric expansion, we show here, that it can also be measured…
▽ More
Understanding $\textit{galaxy bias}$ -- that is the statistical relation between matter and galaxies -- is of key importance for extracting cosmological information from galaxy surveys. While the bias function $f$ -- that is the probability of forming galaxy in a region with a given density field -- is usually approximated through a parametric expansion, we show here, that it can also be measured directly from simulations in a non-parameteric way. Our measurements show that the Lagrangian bias function is very close to a Gaussian for halo selections of any mass. Therefore, we newly introduce a Gaussian bias model with several intriguing properties: (1) It predicts only strictly positive probabilities $f > 0$ (unlike expansion models), (2) It has a simple analytic renormalized form and (3) It behaves gracefully in many scenarios where the classical expansion converges poorly. We show that the Gaussian bias model describes the galaxy environment distribution $p(δ| \mathrm{g})$, the scale dependent bias function $f$ and the renormalized bias function $F$ of haloes and galaxies generally equally well or significantly better than a second order expansion with the same number of parameters. We suggest that a Gaussian bias approach may enhance the range of validity of bias schemes where the canonical expansion converges poorly and further, that it may make new applications possible, since it guarantees the positivity of predicted galaxy densities.
△ Less
Submitted 19 June, 2024; v1 submitted 3 May, 2024;
originally announced May 2024.
-
Probabilistic Lagrangian bias estimators and the cumulant bias expansion
Authors:
Jens Stücker,
Marcos Pellejero-Ibáñez,
Raul E. Angulo,
Francisco Maion,
Rodrigo Voivodic
Abstract:
The spatial distribution of galaxies is a highly complex phenomenon currently impossible to predict deterministically. However, by using a statistical $\textit{bias}$ relation, it becomes possible to robustly model the average abundance of galaxies as a function of the underlying matter density field. Understanding the properties and parametric description of the bias relation is key to extract co…
▽ More
The spatial distribution of galaxies is a highly complex phenomenon currently impossible to predict deterministically. However, by using a statistical $\textit{bias}$ relation, it becomes possible to robustly model the average abundance of galaxies as a function of the underlying matter density field. Understanding the properties and parametric description of the bias relation is key to extract cosmological information from future galaxy surveys. Here, we contribute to this topic primarily in two ways: (1) We develop a new set of probabilistic estimators for bias parameters using the moments of the Lagrangian galaxy environment distribution. These estimators include spatial corrections at different orders to measure bias parameters independently of the damping scale. We report robust measurements of a variety of bias parameters for haloes, including the tidal bias and its dependence with spin at a fixed mass. (2) We propose an alternative formulation of the bias expansion in terms of "cumulant bias parameters" that describe the response of the logarithmic galaxy density to large-scale perturbations. We find that cumulant biases of haloes are consistent with zero at orders $n > 2$. This suggests that: (i) previously reported bias relations at order $n > 2$ are an artefact of the entangled basis of the canonical bias expansion; (ii) the convergence of the bias expansion may be improved by phrasing it in terms of cumulants; (iii) the bias function is very well approximated by a Gaussian -- an avenue which we explore in a companion paper.
△ Less
Submitted 19 June, 2024; v1 submitted 3 May, 2024;
originally announced May 2024.
-
The Effective Field Theory of Large-Scale Structure and Multi-tracer II: redshift space and realistic tracers
Authors:
Thiago Mergulhão,
Henrique Rubira,
Rodrigo Voivodic
Abstract:
We extend the multi-tracer (MT) formalism of the effective field theory of large-scale structure to redshift space, comparing the results of MT to a single-tracer analysis when extracting cosmological parameters from simulations. We used a sub-halo abundance matching method to obtain more realistic multi-tracer galaxy catalogs constructed from N-body simulations. Considering different values for t…
▽ More
We extend the multi-tracer (MT) formalism of the effective field theory of large-scale structure to redshift space, comparing the results of MT to a single-tracer analysis when extracting cosmological parameters from simulations. We used a sub-halo abundance matching method to obtain more realistic multi-tracer galaxy catalogs constructed from N-body simulations. Considering different values for the sample shot noise and volume, we show that the MT error bars on $A_s$, $ω_{\rm cdm}$, and $h$ in a full-shape analysis are approximately $50\%$ smaller relative to ST. We find that cosmological and bias coefficients from MT are less degenerate, indicating that the MT parameter basis is more orthogonal. We conclude that using MT combined with perturbation theory is a robust and competitive way to accommodate the information present in the mildly non-linear scales.
△ Less
Submitted 8 June, 2023;
originally announced June 2023.
-
The Effective Field Theory of Large-Scale Structure and Multi-tracer
Authors:
Thiago Mergulhão,
Henrique Rubira,
Rodrigo Voivodic,
L. Raul Abramo
Abstract:
We study the performance of the perturbative bias expansion when combined with the multi-tracer technique, and their impact on the extraction of cosmological parameters. We consider two populations of tracers of large-scale structure and perform a series of Markov chain Monte Carlo analysis for those two tracers separately. The constraints in $ω_{\rm cdm}$ and $h$ using multi-tracer are less biase…
▽ More
We study the performance of the perturbative bias expansion when combined with the multi-tracer technique, and their impact on the extraction of cosmological parameters. We consider two populations of tracers of large-scale structure and perform a series of Markov chain Monte Carlo analysis for those two tracers separately. The constraints in $ω_{\rm cdm}$ and $h$ using multi-tracer are less biased and approximately $60\%$ better than those obtained for a single tracer. The multi-tracer approach also provides stronger constraints on the bias expansion parameters, breaking degeneracies between them and with their error being typically half of the single-tracer case. Finally, we studied the impacts caused in parameter extraction when including a correlation between the stochastic field of distinct tracers. We also include a study with galaxies showing that multi-tracer still lead to substantial gains in the cosmological parameters.
△ Less
Submitted 13 July, 2022; v1 submitted 25 August, 2021;
originally announced August 2021.
-
Responses of Halo Occupation Distributions: a new ingredient in the halo model & the impact on galaxy bias
Authors:
Rodrigo Voivodic,
Alexandre Barreira
Abstract:
Halo occupation distribution (HOD) models describe the number of galaxies that reside in different haloes, and are widely used in galaxy-halo connection studies using the halo model (HM). Here, we introduce and study HOD response functions $R_\mathcal{O}^g$ that describe the response of the HODs to long-wavelength perturbations $\mathcal{O}$. The linear galaxy bias parameters $b_\mathcal{O}^g$ are…
▽ More
Halo occupation distribution (HOD) models describe the number of galaxies that reside in different haloes, and are widely used in galaxy-halo connection studies using the halo model (HM). Here, we introduce and study HOD response functions $R_\mathcal{O}^g$ that describe the response of the HODs to long-wavelength perturbations $\mathcal{O}$. The linear galaxy bias parameters $b_\mathcal{O}^g$ are a weighted version of $b_\mathcal{O}^h + R_\mathcal{O}^g$, where $b_\mathcal{O}^h$ is the halo bias, but the contribution from $R_\mathcal{O}^g$ is routinely ignored in the literature. We investigate the impact of this by measuring the $R_\mathcal{O}^g$ in separate universe simulations of the IllustrisTNG model for three types of perturbations: total matter perturbations, $\mathcal{O}=δ_m$; baryon-CDM compensated isocurvature perturbations, $\mathcal{O}=σ$; and potential perturbations with local primordial non-Gaussianity, $\mathcal{O}\propto f_{\rm NL}φ$. Our main takeaway message is that the $R_\mathcal{O}^g$ are not negligible in general and their size should be estimated on a case-by-case basis. For stellar-mass selected galaxies, the responses $R_φ^g$ and $R_σ^g$ are sizeable and cannot be neglected in HM calculations of the bias parameters $b_φ^g$ and $b_σ^g$; this is relevant to constrain inflation using galaxies. On the other hand, we do not detect a strong impact of the HOD response $R_1^g$ on the linear galaxy bias $b_1^g$. These results can be explained by the impact that the perturbations have on stellar-to-total-mass relations. We also look into the impact on the bias of the gas distribution and find similar conclusions. We show that a single extra parameter describing the overall amplitude of $R_\mathcal{O}^g$ recovers the measured $b_\mathcal{O}^g$ well, which indicates that $R_\mathcal{O}^g$ can be easily added to HM/HOD studies as a new ingredient.
△ Less
Submitted 8 May, 2021; v1 submitted 8 December, 2020;
originally announced December 2020.
-
The Effective Field Theory and Perturbative Analysis for Log-Density Fields
Authors:
Henrique Rubira,
Rodrigo Voivodic
Abstract:
A logarithm transformation over the matter overdensity field $δ$ brings information from the bispectrum and higher-order n-point functions to the power spectrum. We calculate the power spectrum for the log-transformed field $A$ at one, two and three loops using perturbation theory (PT). We compare the results to simulated data and give evidence that the PT series is asymptotic already on large sca…
▽ More
A logarithm transformation over the matter overdensity field $δ$ brings information from the bispectrum and higher-order n-point functions to the power spectrum. We calculate the power spectrum for the log-transformed field $A$ at one, two and three loops using perturbation theory (PT). We compare the results to simulated data and give evidence that the PT series is asymptotic already on large scales, where the $k$ modes no longer decouple. This motivates us to build an alternative perturbative series for the log-transformed field that is not constructed on top of perturbations of $δ$ but directly over the equations of motion for $A$ itself. This new approach converges faster and better reproduces the large scales at low $z$. We then show that the large-scale behaviour for the log-transformed field power spectrum can be captured by a small number of free parameters. Finally, we add the counter-terms expected within the effective field theory framework and show that the theoretical model, together with the IR-resummation procedure, agrees with the measured spectrum with percent precision until $k \simeq 0.38 $ Mpc$^{-1}$h at $z=0$. It indicates that the non-linear transformation indeed linearizes the density field and, in principle, allows us to access information contained on smaller scales.
△ Less
Submitted 4 February, 2021; v1 submitted 24 November, 2020;
originally announced November 2020.
-
The Halo Void (Dust) Model of Large Scale Structure
Authors:
Rodrigo Voivodic,
Henrique Rubira,
Marcos Lima
Abstract:
Within the Halo Model of large scale structure, all matter is contained in dark matter halos. This simple yet powerful framework has been broadly applied to multiple data sets and enriched our comprehension of how matter is distributed in the Universe. In this work we extend this assumption by allowing for matter to rest not only inside halos but also within cosmic voids and in between halos and v…
▽ More
Within the Halo Model of large scale structure, all matter is contained in dark matter halos. This simple yet powerful framework has been broadly applied to multiple data sets and enriched our comprehension of how matter is distributed in the Universe. In this work we extend this assumption by allowing for matter to rest not only inside halos but also within cosmic voids and in between halos and voids (which we call 'dust'). This assumption leads to additional contributions (1Void, 2Void, Halo-Void, etc.) to the predictions of correlation functions, spectra and profiles for both halos and voids. Whereas the Halo Model can only make predictions for halo quantities, the Halo Void Model extends those for void statistics and halo-void cross-correlations. We provide recipes for all new ingredients of the Halo Void (Dust) Model, such as the void abundance, linear bias and density profile and test their validity in a N-body simulation. Including voids and dust into the calculations improves the transition between the 1Halo and the 2Halo terms by up to $\sim 6\%$. It also eliminates the need to include low-mass structures on the normalization of large-scale terms, suggesting that halos and voids are complementary cosmic structures to effectively describe matter distribution on large scales of the Universe.
△ Less
Submitted 22 January, 2021; v1 submitted 13 March, 2020;
originally announced March 2020.
-
Excursion Set Halos -- ExSHalos: A New Parameter Free Method for Fast Generation of Halo Catalogues
Authors:
Rodrigo Voivodic,
Marcos Lima,
Luis Raul Abramo
Abstract:
We develop a new, simple, fast and parameter-free method to construct dark matter halo catalogues. This method requires as inputs only the linear matter power spectrum and the threshold density for halo formation in linear theory. It directly uses excursion set ideas and Lagrangian perturbation theory to produce halo catalogues with the correct abundance, large scale power spectrum, bispectrum and…
▽ More
We develop a new, simple, fast and parameter-free method to construct dark matter halo catalogues. This method requires as inputs only the linear matter power spectrum and the threshold density for halo formation in linear theory. It directly uses excursion set ideas and Lagrangian perturbation theory to produce halo catalogues with the correct abundance, large scale power spectrum, bispectrum and velocity field. These halo catalogues can be used for the fast construction of mock galaxy catalogues, allowing for the evaluation of covariance matrices for multiple observables. Because of its robustness and predictive nature, this method can be easily adapted to produce catalogues with e.g. primordial non-Gaussianities, modified theories of gravity and non-standard dark energy models, enabling detailed studies of these models in the context of next-generation surveys. We implement this method in a C code, and present numerical comparisons with theoretical predictions as well as full N-body cosmological simulations.
△ Less
Submitted 15 June, 2019;
originally announced June 2019.
-
Cosmic voids in modified gravity scenarios
Authors:
E L D Perico,
R Voivodic,
M Lima,
D F Mota
Abstract:
Modified gravity (MG) theories aim to reproduce the observed acceleration of the Universe by reducing the dark sector while simultaneously recovering General Relativity (GR) within dense environments. Void studies appear to be a suitable scenario to search for imprints of alternative gravity models on cosmological scales. Voids cover an interesting range of density scales where screening mechanism…
▽ More
Modified gravity (MG) theories aim to reproduce the observed acceleration of the Universe by reducing the dark sector while simultaneously recovering General Relativity (GR) within dense environments. Void studies appear to be a suitable scenario to search for imprints of alternative gravity models on cosmological scales. Voids cover an interesting range of density scales where screening mechanisms fade out, which reaches from a density contrast $δ\approx -1$ close to their centers to $δ\approx 0$ close to their boundaries. We present an analysis of the level of distinction between GR and two modified gravity theories, the Hu-Sawicki $f(R)$ and the symmetron theory. This study relies on the abundance, linear bias, and density profile of voids detected in n-body cosmological simulations. We define voids as connected regions made up of the union of spheres with a {\it \textup{mean}} density given by $\overlineρ_v=0.2\,\overlineρ_m$, but disconnected from any other voids. We find that the height of void walls is considerably affected by the gravitational theory, such that it increases for stronger gravity modifications. Finally, we show that at the level of dark matter n-body simulations, our constraints allow us to distinguish between GR and MG models with $|f_{R0}| > 10^{-6}$ and $z_{SSB} > 1$. Differences of best-fit values for MG parameters that are derived independently from multiple void probes may indicate an incorrect MG model. This serves as an important consistency check.
△ Less
Submitted 14 October, 2019; v1 submitted 29 May, 2019;
originally announced May 2019.
-
Relation between the Turnaround radius and virial mass in $f(R)$ model
Authors:
Rafael C. C. Lopes,
Rodrigo Voivodic,
L. Raul Abramo,
Laerte Sodré Jr
Abstract:
We investigate the relationship between the turnaround radius ($R_t$) and the virial mass ($M_v$) of cosmic structures in the context of $Λ$CDM model and in an $f(R)$ model of modified gravity -- namely, the Hu-Sawicki model. The $R_t$ is the distance from the center of the cosmic structure to the shell that is detaching from the Hubble flow at a given time, while the $M_v$ is defined, for this wo…
▽ More
We investigate the relationship between the turnaround radius ($R_t$) and the virial mass ($M_v$) of cosmic structures in the context of $Λ$CDM model and in an $f(R)$ model of modified gravity -- namely, the Hu-Sawicki model. The $R_t$ is the distance from the center of the cosmic structure to the shell that is detaching from the Hubble flow at a given time, while the $M_v$ is defined, for this work, as the mass enclosed within the volume where the density is $200$ times the background density. We consider that gravitationally bound astrophysical systems follow a Navarro-Frenk-White density profile, while beyond the virial radius ($R_v$) the profile is approximated by the 2-halo term of the matter correlation function. By combining them together with the information drawn from solving the spherical collapse for the structures, we are able to connect two observables: the $R_t$ and the $M_v$. We show that, in $Λ$CDM, the turnaround mass ($M_t$) at $z=0$ is related to the $M_v$ of that same structure by $M_t \simeq 3.07 \, M_v$, while in terms of the radii we have that $R_t \simeq 3.7 \, R_v$ (for $M_v$ of $10^{13} \, h^{-1} \, M_\odot$). In the $f(R)$ model, we have $M_t \simeq 3.43 \, M_v$ and $R_t \simeq 4.1 \, R_v$, for $|f_{R0}|=10^{-6}$ and the same mass scale. Therefore, the difference between $Λ$CDM and $f(R)$ in terms of these observable relations is of order $\sim 10-20\%$ even for a relatively mild strength of the modification of gravity ($|f_{R0}|=10^{-6}$). For the $R_t$ itself we find a difference of $\sim 9\%$ between the weakly modification in gravity considered in this work ($|f_{R0}|=10^{-6}$) and $Λ$CDM for a mass of $10^{13} \, h^{-1} \, M_\odot$. Once observations allow precisions of this order or better in measurements $R_t$, as well as the $M_v$, these quantities will become powerful tests of modified gravity.
△ Less
Submitted 26 September, 2018;
originally announced September 2018.
-
Turnaround radius in $f(R)$ model
Authors:
Rafael C. C. Lopes,
Rodrigo Voivodic,
Luis Raul Abramo,
Laerte Sodré Jr
Abstract:
We investigate the turnaround radius in the spherical collapse model, both in General Relativity and in modified gravity, in particular $f(R)$ scenarios. The phases of spherical collapse are marked by the density contrast in the instant of turnaround $δ_t$, and by the linear density contrast in the moment of collapse, $δ_c$. We find that the effective mass of the extra scalar degree of freedom whi…
▽ More
We investigate the turnaround radius in the spherical collapse model, both in General Relativity and in modified gravity, in particular $f(R)$ scenarios. The phases of spherical collapse are marked by the density contrast in the instant of turnaround $δ_t$, and by the linear density contrast in the moment of collapse, $δ_c$. We find that the effective mass of the extra scalar degree of freedom which arises in modified gravity models has an impact on $δ_t$ of up to $\sim10\%$, and that $δ_c$ can increase by $\sim1.0\%$. We also compute the turnaround radius, $R_t$, which in modified gravity models can increase by up to $\sim 6\%$ at $z \simeq 0$.
△ Less
Submitted 24 May, 2018;
originally announced May 2018.
-
Modelling Void Abundance in Modified Gravity
Authors:
Rodrigo Voivodic,
Marcos Lima,
Claudio Llinares,
David F. Mota
Abstract:
We use a spherical model and an extended excursion set formalism with drifting diffusive barriers to predict the abundance of cosmic voids in the context of general relativity as well as f(R) and symmetron models of modified gravity. We detect spherical voids from a suite of N-body simulations of these gravity theories and compare the measured void abundance to theory predictions. We find that our…
▽ More
We use a spherical model and an extended excursion set formalism with drifting diffusive barriers to predict the abundance of cosmic voids in the context of general relativity as well as f(R) and symmetron models of modified gravity. We detect spherical voids from a suite of N-body simulations of these gravity theories and compare the measured void abundance to theory predictions. We find that our model correctly describes the abundance of both dark matter and galaxy voids, providing a better fit than previous proposals in the literature based on static barriers. We use the simulation abundance results to fit for the abundance model free parameters as a function of modified gravity parameters, and show that counts of dark matter voids can provide interesting constraints on modified gravity. For galaxy voids, more closely related to optical observations, we find that constraining modified gravity from void abundance alone may be significantly more challenging. In the context of current and upcoming galaxy surveys, the combination of void and halo statistics including their abundances, profiles and correlations should be effective in distinguishing modified gravity models that display different screening mechanisms.
△ Less
Submitted 8 September, 2016;
originally announced September 2016.