Search

We present the integrated 3-point correlation functions (3PCF) involving both the cosmic shear and the galaxy density fields. These are a set of higher-order statistics that describe the modulation of local 2-point correlation functions (2PCF) by large-scale features in the fields, and which are easy to measure from galaxy imaging surveys. Based on previous works on the shear-only integrated 3PCF, we develop the theoretical framework for modelling 5 new statistics involving the galaxy field and its cross-correlations with cosmic shear. Using realistic galaxy and cosmic shear mocks from simulations, we determine the regime of validity of our models based on leading-order standard perturbation theory with an MCMC analysis that recovers unbiased constraints of the amplitude of fluctuations parameter $A_s$ and the linear and quadratic galaxy bias parameters $b_1$ and $b_2$. Using Fisher matrix forecasts for a DES-Y3-like survey, relative to baseline analyses with conventional 3$\times$2PCFs, we find that the addition of the shear-only integrated 3PCF can improve cosmological parameter constraints by $20-40\%$. The subsequent addition of the new statistics introduced in this paper can lead to further improvements of $10-20\%$, even when utilizing only conservatively large scales where the tree-level models are valid. Our results motivate future work on the galaxy and shear integrated 3PCFs, which offer a practical way to extend standard analyses based on 3$\times$2PCFs to systematically probe the non-Gaussian information content of cosmic density fields., Comment: 19 pages, 8 figures + appendix. Comments are welcome!

We build a deep learning framework that connects the local formation process of dark matter halos to the halo bias. We train a convolutional neural network (CNN) to predict the final mass and concentration of dark matter halos from the initial conditions. The CNN is then used as a surrogate model to derive the response of the halos' mass and concentration to long-wavelength perturbations in the initial conditions, and consequently the halo bias parameters following the "response bias" definition. The CNN correctly predicts how the local properties of dark matter halos respond to changes in the large-scale environment, despite no explicit knowledge of halo bias being provided during training. We show that the CNN recovers the known trends for the linear and second-order density bias parameters $b_1$ and $b_2$, as well as for the local primordial non-Gaussianity linear bias parameter $b_\phi$. The expected secondary assembly bias dependence on halo concentration is also recovered by the CNN: at fixed mass, halo concentration has only a mild impact on $b_1$, but a strong impact on $b_\phi$. Our framework opens a new window for discovering which physical aspects of the halo's Lagrangian patch determine assembly bias, which in turn can inform physical models of halo formation and bias., Comment: 11 pages, 5 figures, to be submitted to MNRAS, comments welcome

The integrated shear 3-point correlation function $\zeta_{\pm}$ measures the correlation between the local shear 2-point function $\xi_{\pm}$ and the 1-point shear aperture mass in patches of the sky. Unlike other higher-order statistics, $\zeta_{\pm}$ can be efficiently measured from cosmic shear data, and it admits accurate theory predictions on a wide range of scales as a function of cosmological and baryonic feedback parameters. Here, we develop and test a likelihood analysis pipeline for cosmological constraints using $\zeta_{\pm}$. We incorporate treatment of systematic effects from photometric redshift uncertainties, shear calibration bias and galaxy intrinsic alignments. We also develop an accurate neural-network emulator for fast theory predictions in MCMC parameter inference analyses. We test our pipeline using realistic cosmic shear maps based on $N$-body simulations with a DES Y3-like footprint, mask and source tomographic bins, finding unbiased parameter constraints. Relative to $\xi_{\pm}$-only, adding $\zeta_{\pm}$ can lead to $\approx 10-25\%$ improvements on the constraints of parameters like $A_s$ (or $\sigma_8$) and $w_0$. We find no evidence in $\xi_{\pm} + \zeta_{\pm}$ constraints of a significant mitigation of the impact of systematics. We also investigate the impact of the size of the apertures where $\zeta_{\pm}$ is measured, and of the strategy to estimate the covariance matrix ($N$-body vs. lognormal). Our analysis solidifies the strong potential of the $\zeta_{\pm}$ statistic and puts forward a pipeline that can be readily used to improve cosmological constraints using real cosmic shear data., Comment: 21 pages, 11 figures, 3 tables. Comments welcome

We discuss the potential of the multi-tracer technique to improve observational constraints of the local primordial non-Gaussianity (PNG) parameter $f_{\rm NL}$ from the galaxy power spectrum. For two galaxy samples $A$ and $B$, the constraining power is $\propto |b_1^B b_\phi^A - b_1^Ab_\phi^B|$, where $b_1$ and $b_\phi$ are the linear and PNG galaxy bias parameters. We show this allows for significantly improved constraints compared to the traditional expectation $\propto |b_1^A - b_1^B|$ based on naive universality-like relations where $b_\phi \propto b_1$. Using IllustrisTNG galaxy simulation data, we find that different equal galaxy number splits of the full sample lead to different $|b_1^B b_\phi^A - b_1^Ab_\phi^B|$, and thus have different constraining power. Of all of the strategies explored, splitting by $g-r$ color is the most promising, more than doubling the significance of detecting $f_{\rm NL}b_\phi \neq 0$. Importantly, since these are constraints on $f_{\rm NL}b_\phi$ and not $f_{\rm NL}$, they do not require priors on the $b_\phi(b_1)$ relation. For direct constraints on $f_{\rm NL}$, we show that multi-tracer constraints can be significantly more robust than single-tracer to $b_\phi$ misspecifications and uncertainties; this relaxes the precision and accuracy requirements for $b_\phi$ priors. Our results present new opportunities to improve our chances to detect and robustly constrain $f_{\rm NL}$, and strongly motivate galaxy formation simulation campaigns to calibrate the $b_\phi(b_1)$ relation., Comment: 9 pages, 4 figures, 1 table. Comments welcomed! v2 matches version published in JCAP

We use the BOSS DR12 galaxy power spectrum to constrain compensated isocurvature perturbations (CIP), which are opposite-sign primordial baryon and dark matter perturbations that leave the total matter density unchanged. Long-wavelength CIP $\sigma(\vec{x})$ enter the galaxy density contrast as $\delta_g(\vec{x}) \supset b_\sigma\sigma(\vec{x})$, with $b_\sigma$ the linear CIP galaxy bias parameter. We parameterize the CIP spectra as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ and $P_{\sigma\mathcal{R}} = \xi\sqrt{P_{\sigma\sigma}P_{\mathcal{R}\mathcal{R}}}$, where $A$ is the CIP amplitude and $\xi$ is the correlation with the curvature perturbations $\mathcal{R}$. We find a significance of detection of $Ab_\sigma \neq 0$ of $1.8\sigma$ for correlated ($\xi = 1$) and $3.7\sigma$ for uncorrelated ($\xi = 0$) CIP. Large-scale data systematics have a bigger impact for uncorrelated CIP, which may explain the large significance of detection. The constraints on $A$ depend on the assumed priors for the $b_\sigma$ parameter, which we estimate using separate universe simulations. Assuming $b_\sigma$ values representative of all halos we find $\sigma_A = 145$ for correlated CIP and $\sigma_{|A|} = 475$ for uncorrelated CIP. Our strongest uncorrelated CIP constraint is for $b_\sigma$ representative of the $33\%$ most concentrated halos, $\sigma_{|A|} = 197$, which is better than the current CMB bounds $|A| \lesssim 360$. We also discuss the impact of the local primordial non-Gaussianity parameter $f_{\rm NL}$ in CIP constraints. Our results demonstrate the power of galaxy data to place tight constraints on CIP, and motivate works to understand better the impact of data systematics, as well as to determine theory priors for $b_\sigma$., Comment: 11 pages, 4 figures, 3 tables. Comments welcomed!

We study the effect of super-sample covariance (SSC) on the power spectrum and higher-order statistics: bispectrum, halo mass function, and void size function. We also investigate the effect of SSC on the cross covariance between the statistics. We consider both the matter and halo fields. Higher-order statistics of the large-scale structure contain additional cosmological information beyond the power spectrum and are a powerful tool to constrain cosmology. They are a promising probe for ongoing and upcoming high precision cosmological surveys such as DESI, PFS, Rubin Observatory LSST, Euclid, SPHEREx, SKA, and Roman Space Telescope. Cosmological simulations used in modeling and validating these statistics often have sizes that are much smaller than the observed Universe. Density fluctuations on scales larger than the simulation box, known as super-sample modes, are not captured by the simulations and in turn can lead to inaccuracies in the covariance matrix. We compare the covariance measured using simulation boxes containing super-sample modes to those without. We also compare with the Separate Universe approach. We find that while the power spectrum, bispectrum and halo mass function show significant scale- or mass-dependent SSC, the void size function shows relatively small SSC. We also find significant SSC contributions to the cross covariances between the different statistics, implying that future joint-analyses will need to carefully take into consideration the effect of SSC. To enable further study of SSC, our simulations have been made publicly available at https://github.com/HalfDomeSims/ssc., Comment: 9+2 pages, 4 figures

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

The scale-dependent bias effect on the galaxy power spectrum is a very promising probe of the local primordial non-Gaussianity (PNG) parameter $f_{\rm NL}$, but the amplitude of the effect is proportional to $f_{\rm NL}b_{\phi}$, where $b_{\phi}$ is the linear PNG galaxy bias parameter. Our knowledge of $b_{\phi}$ is currently very limited, yet nearly all existing $f_{\rm NL}$ constraints and forecasts assume precise knowledge for it. Here, we use the BOSS DR12 galaxy power spectrum to illustrate how our uncertain knowledge of $b_{\phi}$ currently prevents us from constraining $f_{\rm NL}$ with a given statistical precision $\sigma_{f_{\rm NL}}$. Assuming different fixed choices for the relation between $b_{\phi}$ and the linear density bias $b_1$, we find that $\sigma_{f_{\rm NL}}$ can vary by as much as an order of magnitude. Our strongest bound is $f_{\rm NL} = 16 \pm 16\ (1\sigma)$, while the loosest is $f_{\rm NL} = 230 \pm 226\ (1\sigma)$ for the same BOSS data. The impact of $b_{\phi}$ can be especially pronounced because it can be close to zero. We also show how marginalizing over $b_{\phi}$ with wide priors is not conservative, and leads in fact to biased constraints through parameter space projection effects. Independently of galaxy bias assumptions, the scale-dependent bias effect can only be used to detect $f_{\rm NL} \neq 0$ by constraining the product $f_{\rm NL}b_{\phi}$, but the error bar $\sigma_{f_{\rm NL}}$ remains undetermined and the results cannot be compared with the CMB; we find $f_{\rm NL}b_{\phi} \neq 0$ with $1.6\sigma$ significance. We also comment on why these issues are important for analyses with the galaxy bispectrum. Our results strongly motivate simulation-based research programs aimed at robust theoretical priors for the $b_{\phi}$ parameter, without which we may never be able to competitively constrain $f_{\rm NL}$ using galaxy data., Comment: 15 pages, 6 figures, 3 tables. Comments welcomed! v2 Matches version published in JCAP

The integrated shear 3-point correlation function $\zeta_{\pm}$ is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function $\xi_{\pm}$ by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate $\zeta_{\pm}$ that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale $\zeta_{\pm}$ is dominated by the nonlinear matter bispectrum in the squeezed limit, which can be evaluated accurately using the nonlinear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale $\zeta_{\pm}$ measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the nonlinear matter power spectrum, thereby permitting to account for these astrophysical effects on $\zeta_{\pm}$ similarly to how they are currently accounted for on $\xi_{\pm}$. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to $\xi_{\pm}$, a combined $\xi_{\pm}\ \&\ \zeta_{\pm}$ analysis can lead to improvements of order $20-40\%$ on the constraints of cosmological parameters such as $\sigma_8$ or the dark energy equation of state parameter $w_0$. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These are encouraging results that motivate future works on the integrated shear 3-point correlation function towards applications to real survey data., Comment: 15 pages, 9 figures; v2 matches the published version in MNRAS

We use separate universe simulations with the IllustrisTNG galaxy formation model to predict the local PNG bias parameters $b_\phi$ and $b_{\phi\delta}$ of atomic neutral hydrogen, ${\rm H_I}$. These parameters and their relation to the linear density bias parameter $b_1$ play a key role in observational constraints of the local PNG parameter $f_{\rm NL}$ using the ${\rm H_I}$ power spectrum and bispectrum. Our results show that the popular calculation based on the universality of the halo mass function overpredicts the $b_\phi(b_1)$ and $b_{\phi\delta}(b_1)$ relations measured in the simulations. In particular, our results show that at $z \lesssim 1$ the ${\rm H_I}$ power spectrum is more sensitive to $f_{\rm NL}$ compared to previously thought ($b_\phi$ is more negative), but is less sensitive at other epochs ($b_\phi$ is less positive). We discuss how this can be explained by the competition of physical effects such as that large-scale gravitational potentials with local PNG (i) accelerate the conversion of hydrogen to heavy elements by star formation, (ii) enhance the effects of baryonic feedback that eject the gas to regions more exposed to ionizing radiation, and (iii) promote the formation of denser structures that shield the ${\rm H_I}$ more efficiently. Our numerical results can be used to revise existing forecast studies on $f_{\rm NL}$ using 21cm line-intensity mapping data. Despite this first step towards predictions for the local PNG bias parameters of ${\rm H_I}$, we emphasize that more work is needed to assess their sensitivity on the assumed galaxy formation physics and ${\rm H_I}$ modeling strategy., Comment: 13 pages, 5 figures, 1 table. Comments welcomed! v2 matches version published in JCAP

We use hydrodynamical separate universe simulations with the IllustrisTNG model to predict the local primordial non-Gaussianity (PNG) bias parameters $b_{\phi}$ and $b_{\phi\delta}$, which enter at leading order in the galaxy power spectrum and bispectrum. This is the first time that $b_{\phi\delta}$ is measured from either gravity-only or galaxy formation simulations. For dark matter halos, the popular assumption of universality overpredicts the $b_{\phi\delta}(b_1)$ relation in the range $1 \lesssim b_1 \lesssim 3$ by up to $\Delta b_{\phi\delta} \sim 3$ ($b_1$ is the linear density bias). The adequacy of the universality relation is worse for the simulated galaxies, with the relations $b_{\phi}(b_1)$ and $b_{\phi\delta}(b_1)$ being generically redshift-dependent and very sensitive to how galaxies are selected (we test total, stellar and black hole mass, black hole mass accretion rate and color). The uncertainties on $b_{\phi}$ and $b_{\phi\delta}$ have a direct, often overlooked impact on the constraints of the local PNG parameter $f_{\rm NL}$, which we study and discuss. For a survey with $V = 100{\rm Gpc}^3/h^3$ at $z=1$, uncertainties $\Delta b_{\phi} \lesssim 1$ and $\Delta b_{\phi\delta} \lesssim 5$ around values close to the fiducial can yield relatively unbiased constraints on $f_{\rm NL}$ using power spectrum and bispectrum data. We also show why priors on galaxy bias are useful even in analyses that fit for products $f_{\rm NL} b_{\phi}$ and $f_{\rm NL} b_{\phi\delta}$. The strategies we discuss to deal with galaxy bias uncertainties can be straightforwardly implemented in existing $f_{\rm NL}$ constraint analyses (we provide fits for some of the bias relations). Our results motivate more works with galaxy formation simulations to refine our understanding of $b_{\phi}$ and $b_{\phi\delta}$ towards improved constraints on $f_{\rm NL}$., Comment: 18 pages + 2 appendices; 8 figures. Comments welcomed! v2: added references and the new Fig. 6 with an additional forecast (same conclusions). v3: minor comments added following the referee report

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!