Showing 1–39 of 39 results for author: Baldauf, T

The two-loop bispectrum in the effective theory of large-scale structure

Authors: Tobias Baldauf, Mathias Garny, Petter Taule, Theo Steele

Abstract: We study the bispectrum of large-scale structure in the EFTofLSS including corrections up to two-loop. We derive an analytic result for the double-hard limit of the two-loop correction, and show that the UV-sensitivity can be absorbed by the same four EFT operators that renormalize the one-loop bispectrum. For the single-hard region, we employ a simplified treatment, introducing one extra EFT para… ▽ More We study the bispectrum of large-scale structure in the EFTofLSS including corrections up to two-loop. We derive an analytic result for the double-hard limit of the two-loop correction, and show that the UV-sensitivity can be absorbed by the same four EFT operators that renormalize the one-loop bispectrum. For the single-hard region, we employ a simplified treatment, introducing one extra EFT parameter. We compare our results to N-body simulations, and show that going from one- to two-loop extends the wavenumber range with percent-level agreement from $k \simeq 0.08$ to $0.15~h/\mathrm{Mpc}$. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 4 pages, 2 figures. Contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond

The two-loop bispectrum of large-scale structure

Authors: Tobias Baldauf, Mathias Garny, Petter Taule, Theo Steele

Abstract: The bispectrum is the leading non-Gaussian statistic in large-scale structure, carrying valuable information on cosmology that is complementary to the power spectrum. To access this information, we need to model the bispectrum in the weakly non-linear regime. In this work we present the first two-loop, i.e., next-to-next-to-leading order perturbative description of the bispectrum within an effecti… ▽ More The bispectrum is the leading non-Gaussian statistic in large-scale structure, carrying valuable information on cosmology that is complementary to the power spectrum. To access this information, we need to model the bispectrum in the weakly non-linear regime. In this work we present the first two-loop, i.e., next-to-next-to-leading order perturbative description of the bispectrum within an effective field theory (EFT) framework. Using an analytic expansion of the perturbative kernels up to $F_6$ we derive a renormalized bispectrum that is demonstrated to be independent of the UV cutoff. We show that the EFT parameters associated with the four independent second-order EFT operators known from the one-loop bispectrum are sufficient to absorb the UV sensitivity of the two-loop contributions in the double-hard region. In addition, we employ a simplified treatment of the single-hard region, introducing one extra EFT parameter at two-loop order. We compare our results to N-body simulations using the realization-based grid-PT method and find good agreement within the expected range, as well as consistent values for the EFT parameters. The two-loop terms start to become relevant at $k\approx 0.07h~\mathrm{Mpc}^{-1}$. The range of wavenumbers with percent-level agreement, independently of the shape, extends from $0.08h~\mathrm{Mpc}^{-1}$ to $0.15h~\mathrm{Mpc}^{-1}$ when going from one to two loops at $z=0$. In addition, we quantify the impact of using exact instead of Einstein-de-Sitter kernels for the one-loop bispectrum, and discuss in how far their impact can be absorbed into a shift of the EFT parameters. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 34 pages, 17 figures

Report number: TUM-HEP-1360/21

Precise Calibration of the One-Loop Trispectrum in the Effective Field Theory of Large Scale Structure

Authors: Theodore Steele, Tobias Baldauf

Abstract: The Large-Scale Structure (LSS) of the universe has the potential to provide decisive answers to the remaining open questions in cosmology. Early attempts at modelling it analytically focused on using perturbation theory. However, small-scale effects introduced by gravitational collapse cannot be described perturbatively and this failure of perturbation theory is reflected even on the largest scal… ▽ More The Large-Scale Structure (LSS) of the universe has the potential to provide decisive answers to the remaining open questions in cosmology. Early attempts at modelling it analytically focused on using perturbation theory. However, small-scale effects introduced by gravitational collapse cannot be described perturbatively and this failure of perturbation theory is reflected even on the largest scales. The Effective Field Theory of Large Scale Structure (EFTofLSS) has emerged as a consistent method for describing LSS on large scales by introducing counterterms that account for the effects of small-scale dynamics. So far studies of the EFT have mostly focused on the two and three point functions with little attention devoted to the four point function or trispectrum. The trispectrum probes cubic interactions arising from non-linear clustering, biasing, and primordial non-Gaussianities, and constitutes a key element of the covariance matrix of the power spectrum. In this paper, we present explicit calibrations of the EFT counterterms for the one-loop trispectrum. Specifically, we find clear evidence for non-zero EFT corrections. We define two one-parameter ansätze for the counterterm of the one-loop propagator and show that they provide a good correction to the residual at scales below k~0.07 h/Mpc. We then take the amplitudes of the linear and quadratic counterkernels calculated in our previous paper on the bispectrum and use them in the remaining counterterms, establishing consistency of the counterterms in the two, three and four point function. We also show that the commonly used EdS approximation for the growth of the density fields leads to errors that are of the same magnitude as loop corrections to the trispectrum on large scales. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: 32 pages, 21 figures, comments welcome

Journal ref: Phys. Rev. D 103, 103518 (2021)

Non-perturbative halo clustering from cosmological density peaks

Authors: Tobias Baldauf, Sandrine Codis, Vincent Desjacques, Christophe Pichon

Abstract: Associating the formation sites of haloes with the maxima of the smoothed linear density field, we present non-perturbative predictions for the Lagrangian and evolved halo correlation functions that are valid at all separations. In Lagrangian space, we find significant deviations from the perturbative bias calculation at small scales, in particular, a pronounced exclusion region where $ξ=-1$ for m… ▽ More Associating the formation sites of haloes with the maxima of the smoothed linear density field, we present non-perturbative predictions for the Lagrangian and evolved halo correlation functions that are valid at all separations. In Lagrangian space, we find significant deviations from the perturbative bias calculation at small scales, in particular, a pronounced exclusion region where $ξ=-1$ for maxima of unequal height. Our predictions are in good agreement with the Lagrangian clustering of dark matter proto-haloes reconstructed from N-body simulations. Our predictions for the mean infall and velocity dispersion of haloes, which differ from the local bias expansion, show a similar level of agreement with simulations. Finally, we displace the initial density peaks according to the Zeldovich approximation in order to predict the late-time clustering of dark matter haloes. While we are able to reproduce the early evolution of this conserved set of tracers, our approximation fails at the collapse epoch (z=0) on non-linear scales r<10Mpc/h, emphasizing the need for a non-perturbative treatment of the halo displacement field. △ Less

Submitted 28 December, 2020; originally announced December 2020.

Comments: 19 pages, 11 figures, comments welcome

Journal ref: Phys. Rev. D 103, 083530 (2021)

Precise Calibration of the One-Loop Bispectrum in the Effective Field Theory of Large Scale Structure

Authors: Theodore Steele, Tobias Baldauf

Abstract: The bispectrum is the leading non-Gaussian statistic in Large-Scale Structure (LSS) clustering and encodes the interactions in the underlying field. It is thus an important diagnostic for primordial non-Gaussianity and higher order galaxy biasing. In this paper we present a detailed test and calibration of the matter bispectrum counterterms in the Effective Field Theory of LSS against a suite of… ▽ More The bispectrum is the leading non-Gaussian statistic in Large-Scale Structure (LSS) clustering and encodes the interactions in the underlying field. It is thus an important diagnostic for primordial non-Gaussianity and higher order galaxy biasing. In this paper we present a detailed test and calibration of the matter bispectrum counterterms in the Effective Field Theory of LSS against a suite of $N$-body simulations. We are going beyond previous studies in employing realisation based perturbation theory that allows for a significant reduction in cosmic variance error bars. This enables the measurement of the low-energy constants on large scales before two-loop corrections become relevant, around $k<0.09 h\mathrm{Mpc}^{-1}$ at $z=0$. We also go beyond previous work in using bispectrum propagator terms, i.e. correlators with linear and second order fields, to quantify the two new counterterms in isolation and to establish consistency with the power spectrum counterterm. By investigating the fully non-linear bispectrum, $B_{\mathrm{nnn}}$, as well as the terms $B_{\mathrm{n}11}$ and $B_{\mathrm{n}21}$, we find evidence for the new counterterms deviating from the shape suggested by the UV-limit of the relevant bispectrum contributions. We also show that the commonly used Einstein-de Sitter approximation for the time dependence of the tree-level bispectrum is insufficient for precise studies of the one-loop bispectrum and that it is necessary to use $Λ$CDM growth factors in order to obtain meaningful one-loop counterterm constraints. Finally, we also find evidence for small deviations in the growth factors that arise from time integration inaccuracies in the $N$-body simulations. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: 29 pages, 16 figures, comments welcome

Journal ref: Phys. Rev. D 103, 023520 (2021)

Projected two- and three-point statistics: Forecasts and mitigation of non-linear RSDs

Authors: Oliver Leicht, Tobias Baldauf, James Fergusson, Paul Shellard

Abstract: The combination of two- and three-point clustering statistics of galaxies and the underlying matter distribution has the potential to break degeneracies between cosmological parameters and nuisance parameters and can lead to significantly tighter constraints on parameters describing the composition of the Universe and the dynamics of inflation. Here we investigate the relation between biases in th… ▽ More The combination of two- and three-point clustering statistics of galaxies and the underlying matter distribution has the potential to break degeneracies between cosmological parameters and nuisance parameters and can lead to significantly tighter constraints on parameters describing the composition of the Universe and the dynamics of inflation. Here we investigate the relation between biases in the estimated parameters and inaccurate modelling of non-linear redshift-space distortions for the power spectrum and bispectrum of projected galaxy density fields and lensing convergence. Non-linear redshift-space distortions are one of the leading systematic uncertainties in galaxy clustering. Projections along the line of sight suppress radial modes and are thus allowing a trade-off between biases due to non-linear redshift-space distortions and statistical uncertainties. We investigate this bias-error trade-off for a CMASS-like survey with a varying number of redshift bins. Improved modelling of the non-linear redshift-space distortions allows the recovery of more radial information when controlling for biases. Not modelling non-linear redshift space distortions inflates error bars for almost all parameters by 20%. The information loss for the amplitude of local non-Gaussianities is smaller, since it is best constrained from large scales. In addition, we show empirically that one can recover more than 99% of the 3D power spectrum information if the depth of the tomographic bins is reduced to 10 $h^{-1}$Mpc. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 20 pages, 15 figures, 6 tables

Density reconstruction from biased tracers and its application to primordial non-Gaussianity

Authors: Omar Darwish, Simon Foreman, Muntazir M. Abidi, Tobias Baldauf, Blake D. Sherwin, P. Daniel Meerburg

Abstract: Large-scale Fourier modes of the cosmic density field are of great value for learning about cosmology because of their well-understood relationship to fluctuations in the early universe. However, cosmic variance generally limits the statistical precision that can be achieved when constraining model parameters using these modes as measured in galaxy surveys, and moreover, these modes are sometimes… ▽ More Large-scale Fourier modes of the cosmic density field are of great value for learning about cosmology because of their well-understood relationship to fluctuations in the early universe. However, cosmic variance generally limits the statistical precision that can be achieved when constraining model parameters using these modes as measured in galaxy surveys, and moreover, these modes are sometimes inaccessible due to observational systematics or foregrounds. For some applications, both limitations can be circumvented by reconstructing large-scale modes using the correlations they induce between smaller-scale modes of an observed tracer (such as galaxy positions). In this paper, we further develop a formalism for this reconstruction, using a quadratic estimator similar to the one used for lensing of the cosmic microwave background. We incorporate nonlinearities from gravity, nonlinear biasing, and local-type primordial non-Gaussianity, and verify that the estimator gives the expected results when applied to N-body simulations. We then carry out forecasts for several upcoming surveys, demonstrating that, when reconstructed modes are included alongside directly-observed tracer density modes, constraints on local primordial non-Gaussianity are generically tightened by tens of percents compared to standard single-tracer analyses. In certain cases, these improvements arise from cosmic variance cancellation, with reconstructed modes taking the place of modes of a separate tracer, thus enabling an effective "multitracer" approach with single-tracer observations. △ Less

Submitted 16 July, 2020; originally announced July 2020.

Comments: 30 pages plus 14 pages appendices, 19 figures

Primordial non-Gaussianity without tails -- how to measure fNL with the bulk of the density PDF

Authors: Oliver Friedrich, Cora Uhlemann, Francisco Villaescusa-Navarro, Tobias Baldauf, Marc Manera, Takahiro Nishimichi

Abstract: We investigate the possibility to detect primordial non-Gaussianity by analysing the bulk of the probability distribution function (PDF) of late-time cosmic density fluctuations. For this purpose we devise a new method to predict the impact of general non-Gaussian initial conditions on the late-time density PDF. At redshift $z=1$ and for a smoothing scale of 30Mpc/$h$ our predictions agree with th… ▽ More We investigate the possibility to detect primordial non-Gaussianity by analysing the bulk of the probability distribution function (PDF) of late-time cosmic density fluctuations. For this purpose we devise a new method to predict the impact of general non-Gaussian initial conditions on the late-time density PDF. At redshift $z=1$ and for a smoothing scale of 30Mpc/$h$ our predictions agree with the high-resolution Quijote N-body simulations to $\sim 0.2\%$ precision. This is within cosmic variance of a $\sim 100(\mathrm{Gpc}/h)^3$ survey volume. When restricting to this 30Mpc/$h$ smoothing scale and to mildly non-linear densities ($δ[30\mathrm{Mpc}/h] \in [-0.3, 0.4]$) and also marginalizing over potential ignorance of the amplitude of the non-linear power spectrum an analysis of the PDF for such a survey volume can still measure the amplitude of different primordial bispectrum shapes to an accuracy of \smash{$Δf_{\mathrm{NL}}^{\mathrm{loc}}=\pm 7.4\ ,\ Δf_{\mathrm{NL}}^{\mathrm{equi}}=\pm 22.0\ ,\ Δf_{\mathrm{NL}}^{\mathrm{ortho}}=\pm 46.0$} . When pushing to smaller scales and assuming a joint analysis of the PDF with smoothing radii of 30Mpc/$h$ and 15Mpc/$h$ ($δ[15\mathrm{Mpc}/h] \in [-0.4, 0.5]$) this improves to \smash{$Δf_{\mathrm{NL}}^{\mathrm{loc}}=\pm 3.3\ ,\ Δf_{\mathrm{NL}}^{\mathrm{equi}}=\pm 11.0\ ,\ Δf_{\mathrm{NL}}^{\mathrm{ortho}}=\pm 17.0\ $} - even when marginalizing over the non-linear variances at both scales as two free parameters. Especially, such an analysis could simultaneously measure $f_{\mathrm{NL}}$ and the amplitude and slope of the non-linear power spectrum. However, at 15Mpc/$h$ our predictions are only accurate to $\lesssim 0.8\%$ for the considered density range. We discuss how this has to be improved in order to push to these small scales and make full use of upcoming surveys with a PDF-based analysis. △ Less

Submitted 30 July, 2021; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: changes to agree with journal version; 22 pages, 12 figures; public code available at https://github.com/OliverFHD/CosMomentum

Report number: YITP-19-122

Journal ref: Monthly Notices of the Royal Astronomical Society, October 2020, Volume 498, Issue 1, pp.464-483

Primordial Non-Gaussianity

Authors: P. Daniel Meerburg, Daniel Green, Muntazir Abidi, Mustafa A. Amin, Peter Adshead, Zeeshan Ahmed, David Alonso, Behzad Ansarinejad, Robert Armstrong, Santiago Avila, Carlo Baccigalupi, Tobias Baldauf, Mario Ballardini, Kevin Bandura, Nicola Bartolo, Nicholas Battaglia, Daniel Baumann, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Colin Bischoff, Jonathan Blazek, J. Richard Bond, Julian Borrill , et al. (153 additional authors not shown)

Abstract: Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with… ▽ More Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale. △ Less

Submitted 14 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

Cubic Halo Bias in Eulerian and Lagrangian Space

Authors: Muntazir Mehdi Abidi, Tobias Baldauf

Abstract: Predictions of the next-to-leading order, i.e. one-loop, halo power spectra depend on local and non-local bias parameters up to cubic order. The linear bias parameter can be estimated from the large scale limit of the halo-matter power spectrum, and the second order bias parameters from the large scale, tree-level, bispectrum. Cubic operators would naturally be quantified using the tree-level tris… ▽ More Predictions of the next-to-leading order, i.e. one-loop, halo power spectra depend on local and non-local bias parameters up to cubic order. The linear bias parameter can be estimated from the large scale limit of the halo-matter power spectrum, and the second order bias parameters from the large scale, tree-level, bispectrum. Cubic operators would naturally be quantified using the tree-level trispectrum. As the latter is computationally expensive, we extent the quadratic field method proposed in Schmittfull et al. 2014 to cubic fields in order to estimate cubic bias parameters. We cross-correlate a basis set of cubic bias operators with the halo field and express the result in terms of the cross-spectra of these operators in order to cancel cosmic variance. We obtain significant detections of local and non-local cubic bias parameters, which are partially in tension with predictions based on local Lagrangian bias schemes. We directly measure the Lagrangian bias parameters of the protohaloes associated with our halo sample and clearly detect a non-local quadratic term in Lagrangian space. We do not find a clear detection of non-local cubic Lagrangian terms for low mass bins, but there is some mild evidence for their presence for the highest mass bin. While the method presented here focuses on cubic bias parameters, the approach could also be applied to quantifications of cubic primordial non-Gaussianity. △ Less

Submitted 20 July, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 37 pages, 16 figures, 4 tables; figures 14 and 15 are slightly changed compare to the previous version; comments welcome

Journal ref: JCAP07(2018)029

Cosmological Perturbation Theory Using the FFTLog: Formalism and Connection to QFT Loop Integrals

Authors: Marko Simonović, Tobias Baldauf, Matias Zaldarriaga, John Joseph Carrasco, Juna A. Kollmeier

Abstract: We present a new method for calculating loops in cosmological perturbation theory. This method is based on approximating a $Λ$CDM-like cosmology as a finite sum of complex power-law universes. The decomposition is naturally achieved using an FFTLog algorithm. For power-law cosmologies, all loop integrals are formally equivalent to loop integrals of massless quantum field theory. These integrals ha… ▽ More We present a new method for calculating loops in cosmological perturbation theory. This method is based on approximating a $Λ$CDM-like cosmology as a finite sum of complex power-law universes. The decomposition is naturally achieved using an FFTLog algorithm. For power-law cosmologies, all loop integrals are formally equivalent to loop integrals of massless quantum field theory. These integrals have analytic solutions in terms of generalized hypergeometric functions. We provide explicit formulae for the one-loop and the two-loop power spectrum and the one-loop bispectrum. A chief advantage of our approach is that the difficult part of the calculation is cosmology independent, need be done only once, and can be recycled for any relevant predictions. Evaluation of standard loop diagrams then boils down to a simple matrix multiplication. We demonstrate the promise of this method for applications to higher multiplicity/loop correlation functions. △ Less

Submitted 27 August, 2017; originally announced August 2017.

Comments: 42 pages, 12 figures

Iterative initial condition reconstruction

Authors: Marcel Schmittfull, Tobias Baldauf, Matias Zaldarriaga

Abstract: Motivated by recent developments in perturbative calculations of the nonlinear evolution of large-scale structure, we present an iterative algorithm to reconstruct the initial conditions in a given volume starting from the dark matter distribution in real space. In our algorithm, objects are first moved back iteratively along estimated potential gradients, with a progressively reduced smoothing sc… ▽ More Motivated by recent developments in perturbative calculations of the nonlinear evolution of large-scale structure, we present an iterative algorithm to reconstruct the initial conditions in a given volume starting from the dark matter distribution in real space. In our algorithm, objects are first moved back iteratively along estimated potential gradients, with a progressively reduced smoothing scale, until a nearly uniform catalog is obtained. The linear initial density is then estimated as the divergence of the cumulative displacement, with an optional second-order correction. This algorithm should undo nonlinear effects up to one-loop order, including the higher-order infrared resummation piece. We test the method using dark matter simulations in real space. At redshift $z=0$, we find that after eight iterations the reconstructed density is more than $95\%$ correlated with the initial density at $k\le 0.35\; h\mathrm{Mpc}^{-1}$. The reconstruction also reduces the power in the difference between reconstructed and initial fields by more than 2 orders of magnitude at $k\le 0.2\; h\mathrm{Mpc}^{-1}$, and it extends the range of scales where the full broadband shape of the power spectrum matches linear theory by a factor of 2-3. As a specific application, we consider measurements of the baryonic acoustic oscillation (BAO) scale that can be improved by reducing the degradation effects of large-scale flows. In our idealized dark matter simulations, the method improves the BAO signal-to-noise ratio by a factor of 2.7 at $z=0$ and by a factor of 2.5 at $z=0.6$, improving standard BAO reconstruction by $70\%$ at $z=0$ and $30\%$ at $z=0.6$, and matching the optimal BAO signal and signal-to-noise ratio of the linear density in the same volume. For BAO, the iterative nature of the reconstruction is the most important aspect. △ Less

Submitted 27 November, 2020; v1 submitted 21 April, 2017; originally announced April 2017.

Comments: 26 pages, 14 figures, published version. Code available at https://github.com/mschmittfull/iterrec

Journal ref: Phys. Rev. D 96, 023505 (2017)

Phenomenology of BAO evolution from Lagrangian to Eulerian Space

Authors: Tobias Baldauf, Vincent Desjacques

Abstract: The baryon acoustic oscillation (BAO) feature provides an important distance scale for the measurement of the expansion history of the Universe. Theoretical models of the BAO in the distribution of biased tracers of the large scale structure usually rely on an initially linear BAO. With aid of N-body simulations, we demonstrate that the BAO in the initial (Lagrangian) halo 2-point function is sign… ▽ More The baryon acoustic oscillation (BAO) feature provides an important distance scale for the measurement of the expansion history of the Universe. Theoretical models of the BAO in the distribution of biased tracers of the large scale structure usually rely on an initially linear BAO. With aid of N-body simulations, we demonstrate that the BAO in the initial (Lagrangian) halo 2-point function is significantly sharper than in the linear matter distribution, in agreement with peak theory. Using this approach, we delineate the scale-dependence induced by the higher-derivative and velocity bias before assessing how much of the initial BAO enhancement survives until the collapse epoch. Finally, we discuss the extent to which the velocity or gravity bias, which is also imprinted in the displacement field of halos, affects the contrast of the BAO obtained with a reconstruction. △ Less

Submitted 14 December, 2016; originally announced December 2016.

Comments: 17 pages, 5 figures

Journal ref: Phys. Rev. D 95, 043535 (2017)

Verifying the consistency relation for the scale-dependent bias from local primordial non-Gaussianity

Authors: Matteo Biagetti, Titouan Lazeyras, Tobias Baldauf, Vincent Desjacques, Fabian Schmidt

Abstract: 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 Fri… ▽ More 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. △ Less

Submitted 15 May, 2017; v1 submitted 15 November, 2016; originally announced November 2016.

Comments: 12 pages, 6 figures. v2. several improvements: new theory section, new fig.2, added references, matches version published in MNRAS

Journal ref: Mon. Not. Roy. Astron. Soc.,468,3, p. 3277-3288 (2017)

LSS constraints with controlled theoretical uncertainties

Authors: Tobias Baldauf, Mehrdad Mirbabayi, Marko Simonović, Matias Zaldarriaga

Abstract: Forecasts and analyses of cosmological observations often rely on the assumption of a perfect theoretical model over a defined range of scales. We explore how model uncertainties and nuisance parameters in perturbative models of the matter and galaxy spectra affect constraints on neutrino mass and primordial non-Gaussianities. We provide a consistent treatment of theoretical errors and argue that… ▽ More Forecasts and analyses of cosmological observations often rely on the assumption of a perfect theoretical model over a defined range of scales. We explore how model uncertainties and nuisance parameters in perturbative models of the matter and galaxy spectra affect constraints on neutrino mass and primordial non-Gaussianities. We provide a consistent treatment of theoretical errors and argue that their inclusion is a necessary step to obtain realistic cosmological constraints. We find that galaxy surveys up to high redshifts will allow a detection of the minimal neutrino mass and local non-Gaussianity of order unity, but improving the constraints on equilateral non-Gaussianity beyond the CMB limits will be challenging. We argue that similar considerations apply to analyses where theoretical models are based on simulations. △ Less

Submitted 1 February, 2016; originally announced February 2016.

Comments: 18 pages, 9 figures

Linear response to long wavelength fluctuations using curvature simulations

Authors: Tobias Baldauf, Uroš Seljak, Leonardo Senatore, Matias Zaldarriaga

Abstract: We study the local response to long wavelength fluctuations in cosmological $N$-body simulations, focusing on the matter and halo power spectra, halo abundance and non-linear transformations of the density field. The long wavelength mode is implemented using an effective curved cosmology and a mapping of time and distances. The method provides an alternative, most probably more precise, way to mea… ▽ More We study the local response to long wavelength fluctuations in cosmological $N$-body simulations, focusing on the matter and halo power spectra, halo abundance and non-linear transformations of the density field. The long wavelength mode is implemented using an effective curved cosmology and a mapping of time and distances. The method provides an alternative, most probably more precise, way to measure the isotropic halo biases. Limiting ourselves to the linear case, we find generally good agreement between the biases obtained from the curvature method and the traditional power spectrum method at the level of a few percent. We also study the response of halo counts to changes in the variance of the field and find that the slope of the relation between the responses to density and variance differs from the naive derivation assuming a universal mass function by 18%. This has implications for measurements of the amplitude of local non-Gaussianity using scale dependent bias. We also analyze the halo power spectrum and halo-dark matter cross-spectrum response to long wavelength fluctuations and derive second order halo bias from it, as well as the super-sample variance contribution to the galaxy power spectrum covariance matrix. △ Less

Submitted 4 November, 2015; originally announced November 2015.

Comments: 29 pages, 14 figures

Precision measurement of the local bias of dark matter halos

Authors: Titouan Lazeyras, Christian Wagner, Tobias Baldauf, Fabian Schmidt

Abstract: 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 spe… ▽ More 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. △ Less

Submitted 14 April, 2016; v1 submitted 3 November, 2015; originally announced November 2015.

Comments: 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

Journal ref: JCAP 02 (2016) 018

Peak exclusion, stochasticity and convergence of perturbative bias expansions in 1+1 gravity

Authors: Tobias Baldauf, Sandrine Codis, Vincent Desjacques, Christophe Pichon

Abstract: The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spe… ▽ More The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spectrum of these tracers is investigated and shows clear deviations from Poisson noise at low frequencies. On large scales, the convergence of a perturbative bias expansion is discussed. Finally, we go beyond Gaussian statistics for the initial conditions and investigate the subsequent evolution of the two-point clustering of peaks through their Zel'dovich ballistic displacement, to clarify how exclusion effects mix up with scale-dependencies induced by nonlinear gravitational evolution. While the expected large-scale separation limit is recovered, significant deviations are found in the exclusion zone that tends in particular to be reduced at later times. Even though these findings apply to the clustering of one-dimensional tracers, they provide useful insights into halo exclusion and its impact on the two-halo term. △ Less

Submitted 18 December, 2015; v1 submitted 30 October, 2015; originally announced October 2015.

Comments: 16 pages, 9 figures, accepted for publication in MNRAS

The Effective Field Theory of Large Scale Structure at Two Loops: the apparent scale dependence of the speed of sound

Authors: Tobias Baldauf, Lorenzo Mercolli, Matias Zaldarriaga

Abstract: We study the Effective Field Theory of Large Scale Structure for cosmic density and momentum fields. We show that the finite part of the two-loop calculation and its counterterms introduce an apparent scale dependence for the leading order parameter $c_\text{s}^2$ of the EFT starting at k=0.1 h/Mpc. These terms limit the range over which one can trust the one-loop EFT calculation at the 1 % level… ▽ More We study the Effective Field Theory of Large Scale Structure for cosmic density and momentum fields. We show that the finite part of the two-loop calculation and its counterterms introduce an apparent scale dependence for the leading order parameter $c_\text{s}^2$ of the EFT starting at k=0.1 h/Mpc. These terms limit the range over which one can trust the one-loop EFT calculation at the 1 % level to k<0.1 h/Mpc at redshift z=0. We construct a well motivated one parameter ansatz to fix the relative size of the one- and two-loop counterterms using their high-k sensitivity. Although this one parameter model is a very restrictive choice for the counterterms, it explains the apparent scale dependence of $c_\text{s}^2$ seen in simulations. It is also able to capture the scale dependence of the density power spectrum up to k$\approx$ 0.3 h/Mpc at the 1 % level at redshift $z=0$. Considering a simple scheme for the resummation of large scale motions, we find that the two loop calculation reduces the need for this IR-resummation at k<0.2 h/Mpc. Finally, we extend our calculation to momentum statistics and show that the same one parameter model can also describe density-momentum and momentum-momentum statistics. △ Less

Submitted 4 March, 2016; v1 submitted 8 July, 2015; originally announced July 2015.

Comments: 22 pages, 13 figures, fixed typos, results unchanged, published version

Journal ref: Phys. Rev. D 92, 123007 (2015)

On the reach of perturbative methods for dark matter density fields

Authors: Tobias Baldauf, Emmanuel Schaan, Matias Zaldarriaga

Abstract: We study the mapping from Lagrangian to Eulerian space in the context of the Effective Field Theory (EFT) of Large Scale Structure. We compute Lagrangian displacements with Lagrangian Perturbation Theory (LPT) and perform the full non-perturbative transformation from displacement to density. When expanded up to a given order, this transformation reproduces the standard Eulerian Perturbation Theory… ▽ More We study the mapping from Lagrangian to Eulerian space in the context of the Effective Field Theory (EFT) of Large Scale Structure. We compute Lagrangian displacements with Lagrangian Perturbation Theory (LPT) and perform the full non-perturbative transformation from displacement to density. When expanded up to a given order, this transformation reproduces the standard Eulerian Perturbation Theory (SPT) at the same order. However, the full transformation from displacement to density also includes higher order terms. These terms explicitly resum long wavelength motions, thus making the resulting density field better correlated with the true non-linear density field. As a result, the regime of validity of this approach is expected to extend that of the Eulerian EFT, and match that of the IR-resummed Eulerian EFT. This approach thus effectively enables a test of the IR-resummed EFT at the field level. We estimate the size of stochastic, non-perturbative contributions to the matter density power spectrum. We find that in our highest order calculation, at redshift z=0 the power spectrum of the density field is reproduced with an accuracy of 1 % (10 %) up to k=0.25 h/Mpc (k=0.46 h/Mpc). We believe that the dominant source of the remaining error is the stochastic contribution. Unfortunately, on these scales the stochastic term does not yet scale as $k^4$ as it does in the very low-k regime. Thus, modeling this contribution might be challenging. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: 22 pages, 10 figures

On the reach of perturbative descriptions for dark matter displacement fields

Authors: Tobias Baldauf, Emmanuel Schaan, Matias Zaldarriaga

Abstract: We study Lagrangian Perturbation Theory (LPT) and its regularization in the Effective Field Theory (EFT) approach. We evaluate the LPT displacement with the same phases as a corresponding $N$-body simulation, which allows us to compare perturbation theory to the non-linear simulation with significantly reduced cosmic variance, and provides a more stringent test than simply comparing power spectra.… ▽ More We study Lagrangian Perturbation Theory (LPT) and its regularization in the Effective Field Theory (EFT) approach. We evaluate the LPT displacement with the same phases as a corresponding $N$-body simulation, which allows us to compare perturbation theory to the non-linear simulation with significantly reduced cosmic variance, and provides a more stringent test than simply comparing power spectra. We reliably detect a non-vanishing leading order EFT coefficient and a stochastic displacement term, uncorrelated with the LPT terms. This stochastic term is expected in the EFT framework, and, to the best of our understanding, is not an artifact of numerical errors or transients in our simulations. This term constitutes a limit to the accuracy of perturbative descriptions of the displacement field and its phases, corresponding to a $1\%$ error on the non-linear power spectrum at $k=0.2 h$/Mpc at $z=0$. Predicting the displacement power spectrum to higher accuracy or larger wavenumbers thus requires a model for the stochastic displacement. △ Less

Submitted 26 May, 2015; originally announced May 2015.

Comments: 48 pages, 29 figures, comments welcome

Equivalence Principle and the Baryon Acoustic Peak

Authors: Tobias Baldauf, Mehrdad Mirbabayi, Marko Simonović, Matias Zaldarriaga

Abstract: We study the dominant effect of a long wavelength density perturbation $δ(λ_L)$ on short distance physics. In the non-relativistic limit, the result is a uniform acceleration, fixed by the equivalence principle, and typically has no effect on statistical averages due to translational invariance. This same reasoning has been formalized to obtain a "consistency condition" on the cosmological correla… ▽ More We study the dominant effect of a long wavelength density perturbation $δ(λ_L)$ on short distance physics. In the non-relativistic limit, the result is a uniform acceleration, fixed by the equivalence principle, and typically has no effect on statistical averages due to translational invariance. This same reasoning has been formalized to obtain a "consistency condition" on the cosmological correlation functions. In the presence of a feature, such as the acoustic peak at $l_{\rm BAO}$, this naive expectation breaks down for $λ_L<l_{\rm BAO}$. We calculate a universal piece of the three-point correlation function in this regime. The same effect is shown to underlie the spread of the acoustic peak, and is calculable to all orders in the long modes. This can be used to improve the result of perturbative calculations - a technique known as "infra-red resummation" - and is explicitly applied to the one-loop calculation of power spectrum. Finally, the success of BAO reconstruction schemes is argued to be another empirical evidence for the validity of the results. △ Less

Submitted 25 January, 2016; v1 submitted 16 April, 2015; originally announced April 2015.

Comments: 8 pages, 6 figures. Minor changes. Version published in PRD

Journal ref: Phys. Rev. D 92, 043514 (2015)

Testing Inflation with Large Scale Structure: Connecting Hopes with Reality

Authors: Marcelo Alvarez, Tobias Baldauf, J. Richard Bond, Neal Dalal, Roland de Putter, Olivier Doré, Daniel Green, Chris Hirata, Zhiqi Huang, Dragan Huterer, Donghui Jeong, Matthew C. Johnson, Elisabeth Krause, Marilena Loverde, Joel Meyers, P. Daniel Meerburg, Leonardo Senatore, Sarah Shandera, Eva Silverstein, Anže Slosar, Kendrick Smith, Matias Zaldarriaga, Valentin Assassi, Jonathan Braden, Amir Hajian , et al. (3 additional authors not shown)

Abstract: The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing… ▽ More The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $Δf_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014. △ Less

Submitted 15 December, 2014; originally announced December 2014.

Comments: 27 pages + references

Near optimal bispectrum estimators for large-scale structure

Authors: Marcel Schmittfull, Tobias Baldauf, Uroš Seljak

Abstract: Clustering of large-scale structure provides significant cosmological information through the power spectrum of density perturbations. Additional information can be gained from higher-order statistics like the bispectrum, especially to break the degeneracy between the linear halo bias $b_1$ and the amplitude of fluctuations $σ_8$. We propose new simple, computationally inexpensive bispectrum stati… ▽ More Clustering of large-scale structure provides significant cosmological information through the power spectrum of density perturbations. Additional information can be gained from higher-order statistics like the bispectrum, especially to break the degeneracy between the linear halo bias $b_1$ and the amplitude of fluctuations $σ_8$. We propose new simple, computationally inexpensive bispectrum statistics that are near optimal for the specific applications like bias determination. Corresponding to the Legendre decomposition of nonlinear halo bias and gravitational coupling at second order, these statistics are given by the cross-spectra of the density with three quadratic fields: the squared density, a tidal term, and a shift term. For halos and galaxies the first two have associated nonlinear bias terms $b_2$ and $b_{s^2}$, respectively, while the shift term has none in the absence of velocity bias (valid in the $k \rightarrow 0$ limit). Thus the linear bias $b_1$ is best determined by the shift cross-spectrum, while the squared density and tidal cross-spectra mostly tighten constraints on $b_2$ and $b_{s^2}$ once $b_1$ is known. Since the form of the cross-spectra is derived from optimal maximum-likelihood estimation, they contain the full bispectrum information on bias parameters. Perturbative analytical predictions for their expectation values and covariances agree with simulations on large scales, $k\lesssim 0.09h/\mathrm{Mpc}$ at $z=0.55$ with Gaussian $R=20h^{-1}\mathrm{Mpc}$ smoothing, for matter-matter-matter, and matter-matter-halo combinations. For halo-halo-halo cross-spectra the model also needs to include corrections to the Poisson stochasticity. △ Less

Submitted 6 March, 2015; v1 submitted 24 November, 2014; originally announced November 2014.

Comments: 22+6 pages, 11 figures, included minor text improvements to match published version

Journal ref: Phys. Rev. D 91, 043530 (2015)

Lagrangian perturbation theory at one loop order: successes, failures, and improvements

Authors: Zvonimir Vlah, Uroš Seljak, Tobias Baldauf

Abstract: We apply the convolved Lagrangian perturbation theory (CLPT) formalism, in which one can express the matter density power spectrum in terms of integrals over a function of cumulants of the displacement field, allowing for a resummation of the terms, to evaluate the full one loop power spectrum. We keep the cumulants up to third order, extending the Zel'dovich approximation and providing the power… ▽ More We apply the convolved Lagrangian perturbation theory (CLPT) formalism, in which one can express the matter density power spectrum in terms of integrals over a function of cumulants of the displacement field, allowing for a resummation of the terms, to evaluate the full one loop power spectrum. We keep the cumulants up to third order, extending the Zel'dovich approximation and providing the power spectrum analogous to the calculations recently performed for the correlation function. We compare the results to the N-body simulations and to the Lagrangian perturbation simulations up to the second order. We find that the analytic calculations are in a good agreement with the LPT simulations, but when compared to full N-body simulations, we find that while one loop calculations improve upon the Zel'dovich approximation in the power spectrum, they still significantly lack power. As found previously in the correlation function one loop CLPT improves slightly against Zel'dovich above 30Mpc/h, but is actually worse than Zel'dovich below that. We investigate the deficiencies of the CLPT approach and argue that main problem of CLPT is its inability to trap particles inside dark matter halos, which leads to an overestimate of the small scale power of the displacement field and to an underestimate of the small scale power from one halo term effects. We model this using the displacement field damped at a nonlinear scale (CLPTs). To explore this in more detail we decompose the power spectrum and correlation function into three additive components: Zel'dovich, residual BAO wiggle, and residual broad band. One loop CLPT predicts small modifications to BAO wiggles that are enhanced in CLPTs, with up to 5\% corrections to correlation function around BAO scale. △ Less

Submitted 7 October, 2014; originally announced October 2014.

Comments: 19 pages, 8 figures

Journal ref: Phys. Rev. D 91, 023508 (2015)

The Bispectrum in the Effective Field Theory of Large Scale Structure

Authors: Tobias Baldauf, Lorenzo Mercolli, Mehrdad Mirbabayi, Enrico Pajer

Abstract: We study the bispectrum in the Effective Field Theory of Large Scale Structure, consistently accounting for the effects of short-scale dynamics. We begin by proving that, as long as the theory is perturbative, it can be formulated to arbitrary order using only operators that are local in time. We then derive all the new operators required to cancel the UV-divergences and obtain a physically meanin… ▽ More We study the bispectrum in the Effective Field Theory of Large Scale Structure, consistently accounting for the effects of short-scale dynamics. We begin by proving that, as long as the theory is perturbative, it can be formulated to arbitrary order using only operators that are local in time. We then derive all the new operators required to cancel the UV-divergences and obtain a physically meaningful prediction for the one-loop bispectrum. In addition to new, subleading stochastic noises and the viscosity term needed for the one-loop power spectrum, we find three new effective operators. The three new parameters can be constrained by comparing with N-body simulations. The best fit is precisely what is suggested by the structure of UV-divergences, hence justifying a formula for the EFTofLSS bispectrum whose only fitting parameter is already fixed by the power spectrum. This result predicts the bispectrum of N-body simulations up to $k \approx 0.22\, h\, \text{Mpc}^{-1}$ at $z=0$, an improvement by nearly a factor of two as compared to one-loop standard perturbation theory. △ Less

Submitted 14 April, 2015; v1 submitted 16 June, 2014; originally announced June 2014.

Comments: 63 pages, 18 figures, as published in JCAP

Velocity bias in the distribution of dark matter halos

Authors: Tobias Baldauf, Vincent Desjacques, Uroš Seljak

Abstract: The standard formalism for the co-evolution of halos and dark matter predicts that any initial halo velocity bias rapidly decays to zero. We argue that, when the purpose is to compute statistics like power spectra etc., the coupling in the momentum conservation equation for the biased tracers must be modified. Our new formulation predicts the constancy in time of any statistical halo velocity bias… ▽ More The standard formalism for the co-evolution of halos and dark matter predicts that any initial halo velocity bias rapidly decays to zero. We argue that, when the purpose is to compute statistics like power spectra etc., the coupling in the momentum conservation equation for the biased tracers must be modified. Our new formulation predicts the constancy in time of any statistical halo velocity bias present in the initial conditions, in agreement with peak theory. We test this prediction by studying the evolution of a conserved halo population in N-body simulations. We establish that the initial simulated halo density and velocity statistics show distinct features of the peak model and, thus, deviate from the simple local Lagrangian bias. We demonstrate, for the first time, that the time evolution of their velocity is in tension with the rapid decay expected in the standard approach. △ Less

Submitted 5 January, 2016; v1 submitted 22 May, 2014; originally announced May 2014.

Comments: 6+ pages, extended to match published version, conclusions unchanged

Journal ref: Phys. Rev. D 92, 123507 (2015)

Understanding higher-order nonlocal halo bias at large scales by combining the power spectrum with the bispectrum

Authors: Shun Saito, Tobias Baldauf, Zvonimir Vlah, Uroš Seljak, Teppei Okumura, Patrick McDonald

Abstract: Understanding the relation between underlying matter distribution and biased tracers such as galaxy or dark matter halo is essential to extract cosmological information from ongoing or future galaxy redshift surveys. At sufficiently large scales such as the BAO scale, a standard approach for the bias problem on the basis of the perturbation theory (PT) is to assume the `local bias' model in which… ▽ More Understanding the relation between underlying matter distribution and biased tracers such as galaxy or dark matter halo is essential to extract cosmological information from ongoing or future galaxy redshift surveys. At sufficiently large scales such as the BAO scale, a standard approach for the bias problem on the basis of the perturbation theory (PT) is to assume the `local bias' model in which the density field of biased tracers is deterministically expanded in terms of matter density field at the same position. The higher-order bias parameters are then determined by combining the power spectrum with higher-order statistics such as the bispectrum. As is pointed out by recent studies, however, nonlinear gravitational evolution naturally induces nonlocal bias terms even if initially starting only with purely local bias. As a matter of fact, previous works showed that the second-order nonlocal bias term, which corresponds to the gravitational tidal field, is important to explain the characteristic scale-dependence of the bispectrum. In this paper we extend the nonlocal bias term up to third order, and investigate whether the PT-based model including nonlocal bias terms can simultaneously explain the power spectrum and the bispectrum of simulated halos in $N$-body simulations. We show that the power spectrum, including density and momentum, and the bispectrum between halo and matter in $N$-body simulations can be simultaneously well explained by the model including up to third-order nonlocal bias terms up to k~0.1h/Mpc. Also, the results seem in a good agreement with theoretical predictions of a simple coevolution picture, although the agreement is not perfect. These demonstration clearly shows a failure of the local bias model even at such large scales, and we conclude that nonlocal bias terms should be consistently included in order to model statistics of halos. [abridged] △ Less

Submitted 17 December, 2014; v1 submitted 6 May, 2014; originally announced May 2014.

Comments: 29 pages, matched to the accepted version

Report number: IPMU14-0105

Journal ref: Phys. Rev. D 90, 123522 (2014)

Halo Stochasticity from Exclusion and non-linear Clustering

Authors: Tobias Baldauf, Uroš Seljak, Robert E. Smith, Nico Hamaus, Vincent Desjacques

Abstract: The clustering of galaxies in ongoing and upcoming galaxy surveys contains a wealth of cosmological information, but extracting this information is a non-trivial task since galaxies and their host haloes are stochastic tracers of the matter density field. This stochasticity is usually modeled as Poisson shot noise, which is constant as a function of wavenumber with amplitude given by 1/n, where n… ▽ More The clustering of galaxies in ongoing and upcoming galaxy surveys contains a wealth of cosmological information, but extracting this information is a non-trivial task since galaxies and their host haloes are stochastic tracers of the matter density field. This stochasticity is usually modeled as Poisson shot noise, which is constant as a function of wavenumber with amplitude given by 1/n, where n is the number density of galaxies. Here we use dark matter haloes in N-body simulations to show evidence for deviations from this simple behaviour and develop models that explain the behaviour of the large scale stochasticity. First, haloes are extended, non-overlapping objects, i.e., their correlation function needs to go to -1 on small scales. This leads to a negative correction to the stochasticity relative to the Poisson value at low wavenumber k, decreasing to zero for wavenumbers large compared to the inverse exclusion scale. Second, haloes show a non-linear enhancement of clustering outside the exclusion scale, leading to a positive stochasticity correction. Both of these effects go to zero for high-k, making the stochasticity scale dependent even for k<0.1 h/Mpc. We show that the corrections in the low-k regime are the same in Eulerian and Lagrangian space, but that the transition scale is pushed to smaller scales for haloes observed at present time, relative to the initial conditions. These corrections vary with halo mass and redshift. We also discuss simple applications of these effects to the galaxy samples with non-vanishing satellite fraction, where the stochasticity can again deviate strongly from the fiducial Poisson expectation. Overall these effects affect the clustering of galaxies at a level of a few percent even on very large scales and need to be modelled properly if we want to extract high precision cosmological information from the upcoming galaxy surveys. △ Less

Submitted 13 May, 2013; originally announced May 2013.

Comments: 27 pages, 15 figures

Journal ref: Phys. Rev. D 88, 083507 (2013)

Cosmological parameter constraints from galaxy-galaxy lensing and galaxy clustering with the SDSS DR7

Authors: Rachel Mandelbaum, Anze Slosar, Tobias Baldauf, Uros Seljak, Christopher M. Hirata, Reiko Nakajima, Reinabelle Reyes, Robert E. Smith

Abstract: Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy halos, independent of the details of how galaxies populate dark matter halos. This finding makes it possible to determine the dark matter clustering from measurements of galaxy-galaxy weak lensing and galaxy clustering. We present n… ▽ More Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy halos, independent of the details of how galaxies populate dark matter halos. This finding makes it possible to determine the dark matter clustering from measurements of galaxy-galaxy weak lensing and galaxy clustering. We present new cosmological parameter constraints based on large-scale measurements of spectroscopic galaxy samples from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We generalise the approach of Baldauf et al. (2010) to remove small scale information (below 2 and 4 Mpc/h for lensing and clustering measurements, respectively), where the cross-correlation coefficient differs from unity. We derive constraints for three galaxy samples covering 7131 sq. deg., containing 69150, 62150, and 35088 galaxies with mean redshifts of 0.11, 0.28, and 0.40. We clearly detect scale-dependent galaxy bias for the more luminous galaxy samples, at a level consistent with theoretical expectations. When we vary both σ_8 and Ω_m (and marginalise over non-linear galaxy bias) in a flat LCDM model, the best-constrained quantity is σ_8 (Ω_m/0.25)^{0.57}=0.80 +/- 0.05 (1-sigma, stat. + sys.), where statistical and systematic errors have comparable contributions, and we fixed n_s=0.96 and h=0.7. These strong constraints on the matter clustering suggest that this method is competitive with cosmic shear in current data, while having very complementary and in some ways less serious systematics. We therefore expect that this method will play a prominent role in future weak lensing surveys. When we combine these data with WMAP7 CMB data, constraints on σ_8, Ω_m, H_0, w_{de} and \sum m_ν become 30--80 per cent tighter than with CMB data alone, since our data break several parameter degeneracies. △ Less

Submitted 9 March, 2013; v1 submitted 4 July, 2012; originally announced July 2012.

Comments: 34 pages, 26 figures; submitted to MNRAS; textual revisions in response to referee comments, but no changes in results. The data used for cosmological parameter constraints and the resulting MCMC chains can be downloaded from the first author's website, http://www.andrew.cmu.edu/user/rmandelb/data.html

Distribution function approach to redshift space distortions. Part IV: perturbation theory applied to dark matter

Authors: Zvonimir Vlah, Uroš Seljak, Patrick McDonald, Teppei Okumura, Tobias Baldauf

Abstract: We develop a perturbative approach to redshift space distortions (RSD) using the phase space distribution function approach and apply it to the dark matter redshift space power spectrum and its moments. RSD can be written as a sum over density weighted velocity moments correlators, with the lowest order being density, momentum density and stress energy density. We use standard and extended perturb… ▽ More We develop a perturbative approach to redshift space distortions (RSD) using the phase space distribution function approach and apply it to the dark matter redshift space power spectrum and its moments. RSD can be written as a sum over density weighted velocity moments correlators, with the lowest order being density, momentum density and stress energy density. We use standard and extended perturbation theory (PT) to determine their auto and cross correlators, comparing them to N-body simulations. We show which of the terms can be modeled well with the standard PT and which need additional terms that include higher order corrections which cannot be modeled in PT. Most of these additional terms are related to the small scale velocity dispersion effects, the so called finger of god (FoG) effects, which affect some, but not all, of the terms in this expansion, and which can be approximately modeled using a simple physically motivated ansatz such as the halo model. We point out that there are several velocity dispersions that enter into the detailed RSD analysis with very different amplitudes, which can be approximately predicted by the halo model. In contrast to previous models our approach systematically includes all of the terms at a given order in PT and provides a physical interpretation for the small scale dispersion values. We investigate RSD power spectrum as a function of μ, the cosine of the angle between the Fourier mode and line of sight, focusing on the lowest order powers of μand multipole moments which dominate the observable RSD power spectrum. Overall we find considerable success in modeling many, but not all, of the terms in this expansion. △ Less

Submitted 5 December, 2012; v1 submitted 3 July, 2012; originally announced July 2012.

Comments: 37 pages, 13 figures, published in JCAP

Journal ref: JCAP 11 (2012) 009

Evidence for Quadratic Tidal Tensor Bias from the Halo Bispectrum

Authors: Tobias Baldauf, Uros Seljak, Vincent Desjacques, Patrick McDonald

Abstract: The relation between the clustering properties of luminous matter in the form of galaxies and the underlying dark matter distribution is of fundamental importance for the interpretation of ongoing and upcoming galaxy surveys. The so called local bias model, where galaxy density is a function of local matter density, is frequently discussed as a means to infer the matter power spectrum or correlati… ▽ More The relation between the clustering properties of luminous matter in the form of galaxies and the underlying dark matter distribution is of fundamental importance for the interpretation of ongoing and upcoming galaxy surveys. The so called local bias model, where galaxy density is a function of local matter density, is frequently discussed as a means to infer the matter power spectrum or correlation function from the measured galaxy correlation. However, gravitational evolution generates a term quadratic in the tidal tensor and thus non-local in the density field, even if this term is absent in the initial conditions (Lagrangian space). Because the term is quadratic, it contributes as a loop correction to the power spectrum, so the standard linear bias picture still applies on large scales, however, it contributes at leading order to the bispectrum for which it is significant on all scales. Such a term could also be present in Lagrangian space if halo formation were influenced by the tidal field. We measure the corresponding coupling strengths from the matter-matter-halo bispectrum in numerical simulations and find a non-vanishing coefficient for the tidal tensor term. We find no scale dependence of the bias parameters up to k=0.1 h/Mpc and that the tidal effect is increasing with halo mass. While the Lagrangian bias picture is a better description of our results than the Eulerian bias picture, our results suggest that there might be a tidal tensor bias already in the initial conditions. We also find that the coefficients of the quadratic density term deviate quite strongly from the theoretical predictions based on the spherical collapse model and a universal mass function. Both quadratic density and tidal tensor bias terms must be included in the modeling of galaxy clustering of current and future surveys if one wants to achieve the high precision cosmology promise of these datasets. △ Less

Submitted 23 January, 2012; originally announced January 2012.

Comments: 14 pages, 4 figures, 1 table

Cluster Density Profiles as a Test of Modified Gravity

Authors: Lucas Lombriser, Fabian Schmidt, Tobias Baldauf, Rachel Mandelbaum, Uros Seljak, Robert E. Smith

Abstract: We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 9… ▽ More We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 95% confidence level. As a model-independent assessment of the constraining power of cluster profiles measured through weak gravitational lensing, we also constrain the amplitude F_0 of a phenomenological modification based on the profile enhancement induced by f(R) gravity when not including effects from the increased cluster abundance in f(R). In both scenarios, dark-matter-only simulations of the concordance model corresponding to f_R0=0 and F_0=0 are consistent with the lensing measurements, i.e., at the 68% and 95% confidence level, respectively. △ Less

Submitted 28 May, 2012; v1 submitted 8 November, 2011; originally announced November 2011.

Comments: 19 pages, 10 figures, 3 tables; new figure added to new version, removed F_0>0 prior

Journal ref: Phys.Rev.D85:102001,2012

Galaxy Bias and non-Linear Structure Formation in General Relativity

Authors: Tobias Baldauf, Uros Seljak, Leonardo Senatore, Matias Zaldarriaga

Abstract: Length scales probed by large scale structure surveys are becoming closer to the horizon scale. Further, it has been recently understood that non-Gaussianity in the initial conditions could show up in a scale dependence of the bias of galaxies at the largest distances. It is therefore important to include General Relativistic effects. Here we provide a General Relativistic generalization of the bi… ▽ More Length scales probed by large scale structure surveys are becoming closer to the horizon scale. Further, it has been recently understood that non-Gaussianity in the initial conditions could show up in a scale dependence of the bias of galaxies at the largest distances. It is therefore important to include General Relativistic effects. Here we provide a General Relativistic generalization of the bias, valid both for Gaussian and non-Gaussian initial conditions. The collapse of objects happens on very small scales, while long-wavelength modes are always in the quasi linear regime. Around every collapsing region, it is therefore possible to find a reference frame that is valid for all times and where the space time is almost flat: the Fermi frame. Here the Newtonian approximation is applicable and the equations of motion are the ones of the N-body codes. The effects of long-wavelength modes are encoded in the mapping from the cosmological frame to the local frame. For the linear bias, the effect of the long-wavelength modes on the dynamics is encoded in the local curvature of the Universe, which allows us to define a General Relativistic generalization of the bias in the standard Newtonian setting. We show that the bias due to this effect goes to zero as the squared ratio of the physical wavenumber with the Hubble scale for modes longer than the horizon, as modes longer than the horizon have no dynamical effects. However, the bias due to non-Gaussianities does not need to vanish for modes longer than the Hubble scale, and for non-Gaussianities of the local kind it goes to a constant. As a further application, we show that it is not necessary to perform large N-body simulations to extract information on long-wavelength modes: N-body simulations can be done on small scales and long-wavelength modes are encoded simply by adding curvature to the simulation and rescaling the coordinates. △ Less

Submitted 30 October, 2011; v1 submitted 27 June, 2011; originally announced June 2011.

Comments: 48 pages, 4 figures; v2: added references, JCAP published version

Primordial non-Gaussianity in the Bispectrum of the Halo Density Field

Authors: Tobias Baldauf, Uros Seljak, Leonardo Senatore

Abstract: The bispectrum vanishes for linear Gaussian fields and is thus a sensitive probe of non-linearities and non-Gaussianities in the cosmic density field. Hence, a detection of the bispectrum in the halo density field would enable tight constraints on non-Gaussian processes in the early Universe and allow inference of the dynamics driving inflation. We present a tree level derivation of the halo bispe… ▽ More The bispectrum vanishes for linear Gaussian fields and is thus a sensitive probe of non-linearities and non-Gaussianities in the cosmic density field. Hence, a detection of the bispectrum in the halo density field would enable tight constraints on non-Gaussian processes in the early Universe and allow inference of the dynamics driving inflation. We present a tree level derivation of the halo bispectrum arising from non-linear clustering, non-linear biasing and primordial non-Gaussianity. A diagrammatic description is developed to provide an intuitive understanding of the contributing terms and their dependence on scale, shape and the non-Gaussianity parameter fNL. We compute the terms based on a multivariate bias expansion and the peak-background split method and show that non-Gaussian modifications to the bias parameters lead to amplifications of the tree level bispectrum that were ignored in previous studies. Our results are in a good agreement with published simulation measurements of the halo bispectrum. Finally, we estimate the expected signal to noise on fNL and show that the constraint obtainable from the bispectrum analysis significantly exceeds the one obtainable from the power spectrum analysis. △ Less

Submitted 17 August, 2011; v1 submitted 5 November, 2010; originally announced November 2010.

Comments: 34 pages, 15 figures, (v3): matches JCAP published version

Journal ref: JCAP 04 (2011) 006

Minimizing the stochasticity of halos in large-scale structure surveys

Authors: Nico Hamaus, Uros Seljak, Vincent Desjacques, Robert E. Smith, Tobias Baldauf

Abstract: In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of hal… ▽ More In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of halo stochasticity and its dependence on halo mass. We investigate the stochasticity matrix, defined as $C_{ij}\equiv<(δ_i -b_iδ_m)(δ_j-b_jδ_m)>$, where $δ_m$ is the dark matter overdensity in Fourier space, $δ_i$ the halo overdensity of the $i$-th halo mass bin and $b_i$ the halo bias. In contrast to the Poisson model predictions we detect nonvanishing correlations between different mass bins. We also find the diagonal terms to be sub-Poissonian for the highest-mass halos. The diagonalization of this matrix results in one large and one low eigenvalue, with the remaining eigenvalues close to the Poisson prediction $1/\bar{n}$, where $\bar{n}$ is the mean halo number density. The eigenmode with the lowest eigenvalue contains most of the information and the corresponding eigenvector provides an optimal weighting function to minimize the stochasticity between halos and dark matter. We find this optimal weighting function to match linear mass weighting at high masses, while at the low-mass end the weights approach a constant whose value depends on the low-mass cut in the halo mass function. Finally, we employ the halo model to derive the stochasticity matrix and the scale-dependent bias from an analytical perspective. It is remarkably successful in reproducing our numerical results and predicts that the stochasticity between halos and the dark matter can be reduced further when going to halo masses lower than we can resolve in current simulations. △ Less

Submitted 16 August, 2010; v1 submitted 29 April, 2010; originally announced April 2010.

Comments: 17 pages, 14 figures, matched the published version in Phys. Rev. D including one new figure

Journal ref: Phys.Rev.D82:043515,2010

Confirmation of general relativity on large scales from weak lensing and galaxy velocities

Authors: Reinabelle Reyes, Rachel Mandelbaum, Uros Seljak, Tobias Baldauf, James E. Gunn, Lucas Lombriser, Robert E. Smith

Abstract: Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to 'galaxy bias' (the difference between the clustering of visi… ▽ More Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to 'galaxy bias' (the difference between the clustering of visible galaxies and invisible dark matter) and is thus robust to the uncertainty in this parameter. Modified theories of gravity generally predict values of EG different from the general relativistic prediction because, in these theories, the 'gravitational slip' (the difference between the two potentials that describe perturbations in the gravitational metric) is non-zero, which leads to changes in the growth of structure and the strength of the gravitational lensing effect3. Here we report that EG = 0.39 +/- 0.06 on length scales of tens of megaparsecs, in agreement with the general relativistic prediction of EG $\approx$ 0.4. The measured value excludes a model within the tensor-vector-scalar gravity theory, which modifies both Newtonian and Einstein gravity. However, the relatively large uncertainty still permits models within f(R) theory, which is an extension of general relativity. A fivefold decrease in uncertainty is needed to rule out these models. △ Less

Submitted 10 March, 2010; originally announced March 2010.

Comments: Submitted version; 13 pages, 2 figures. Accepted version and supplementary material are available at: http://www.nature.com/nature/journal/v464/n7286/full/nature08857.html.

Journal ref: Reyes, R. et al. 2010, Nature, 464, 256-258.

An algorithm for the direct reconstruction of the dark matter correlation function from weak lensing and galaxy clustering

Authors: Tobias Baldauf, Robert E. Smith, Uros Seljak, Rachel Mandelbaum

Abstract: The clustering of matter on cosmological scales is an essential probe for studying the physical origin and composition of our Universe. To date, most of the direct studies have focused on shear-shear weak lensing correlations, but it is also possible to extract the dark matter clustering by combining galaxy-clustering and galaxy-galaxy-lensing measurements. In this study we develop a method that c… ▽ More The clustering of matter on cosmological scales is an essential probe for studying the physical origin and composition of our Universe. To date, most of the direct studies have focused on shear-shear weak lensing correlations, but it is also possible to extract the dark matter clustering by combining galaxy-clustering and galaxy-galaxy-lensing measurements. In this study we develop a method that can constrain the dark matter correlation function from galaxy clustering and galaxy-galaxy-lensing measurements, by focusing on the correlation coefficient between the galaxy and matter overdensity fields. To generate a mock galaxy catalogue for testing purposes, we use the Halo Occupation Distribution approach applied to a large ensemble of N-body simulations to model pre-existing SDSS Luminous Red Galaxy sample observations. Using this mock catalogue, we show that a direct comparison between the excess surface mass density measured by lensing and its corresponding galaxy clustering quantity is not optimal. We develop a new statistic that suppresses the small-scale contributions to these observations and show that this new statistic leads to a cross-correlation coefficient that is within a few percent of unity down to 5 Mpc/h. Furthermore, the residual incoherence between the galaxy and matter fields can be explained using a theoretical model for scale-dependent bias, giving us a final estimator that is unbiased to within 1%. We also perform a comprehensive study of other physical effects that can affect the analysis, such as redshift space distortions and differences in radial windows between galaxy clustering and weak lensing observations. We apply the method to a range of cosmological models and show the viability of our new statistic to distinguish between cosmological models. △ Less

Submitted 23 July, 2010; v1 submitted 25 November, 2009; originally announced November 2009.

Comments: 23 pages, 14 figures, accepted by PRD; minor changes to V1, 1 new figure, more detailed discussion of the covariance of the new ADSD statistic

Journal ref: Phys.Rev.D81:063531,2010

Precision cluster mass determination from weak lensing

Authors: Rachel Mandelbaum, Uros Seljak, Tobias Baldauf, Robert E. Smith

Abstract: Weak gravitational lensing has been used extensively in the past decade to constrain the masses of galaxy clusters, and is the most promising observational technique for providing the mass calibration necessary for precision cosmology with clusters. There are several challenges in estimating cluster masses, particularly (a) the sensitivity to astrophysical effects and observational systematics t… ▽ More Weak gravitational lensing has been used extensively in the past decade to constrain the masses of galaxy clusters, and is the most promising observational technique for providing the mass calibration necessary for precision cosmology with clusters. There are several challenges in estimating cluster masses, particularly (a) the sensitivity to astrophysical effects and observational systematics that modify the signal relative to the theoretical expectations, and (b) biases that can arise due to assumptions in the mass estimation method, such as the assumed radial profile of the cluster. All of these challenges are more problematic in the inner regions of the cluster, suggesting that their influence would ideally be suppressed for the purpose of mass estimation. However, at any given radius the differential surface density measured by lensing is sensitive to all mass within that radius, and the corrupted signal from the inner parts is spread out to all scales. We develop a new statistic that is ideal for estimation of cluster masses because it completely eliminates mass contributions below a chosen scale (which we suggest should be about 20 per cent of the virial radius), and thus reduces sensitivity to systematic and astrophysical effects. We use simulated and analytical profiles to quantify systematic biases on the estimated masses for several standard methods of mass estimation, finding that these can lead to significant mass biases that range from ten to over fifty per cent. The mass uncertainties when using our new statistic are reduced by up to a factor of ten relative to the standard methods, while only moderately increasing the statistical errors. This new method of mass estimation will enable a higher level of precision in future science work with weak lensing mass estimates for galaxy clusters. △ Less

Submitted 19 February, 2010; v1 submitted 25 November, 2009; originally announced November 2009.

Comments: 27 pages, 7 figures, submitted to MNRAS; v2 has expanded explanation for clarity, no change in results or conclusions

Journal ref: Mon.Not.Roy.Astron.Soc.405:2078,2010