Search

We infer the growth of large scale structure over the redshift range $0.4\lesssim z \lesssim 1$ from the cross-correlation of spectroscopically calibrated Luminous Red Galaxies (LRGs) selected from the Dark Energy Spectroscopic Instrument (DESI) legacy imaging survey with CMB lensing maps reconstructed from the latest Planck and ACT data. We adopt a hybrid effective field theory (HEFT) model that robustly regulates the cosmological information obtainable from smaller scales, such that our cosmological constraints are reliably derived from the (predominantly) linear regime. We perform an extensive set of bandpower- and parameter-level systematics checks to ensure the robustness of our results and to characterize the uniformity of the LRG sample. We demonstrate that our results are stable to a wide range of modeling assumptions, finding excellent agreement with a linear theory analysis performed on a restricted range of scales. From a tomographic analysis of the four LRG photometric redshift bins we find that the rate of structure growth is consistent with $\Lambda$CDM with an overall amplitude that is $\simeq5-7\%$ lower than predicted by primary CMB measurements with modest $(\sim2\sigma)$ statistical significance. From the combined analysis of all four bins and their cross-correlations with Planck we obtain $S_8 = 0.765\pm0.023$, which is less discrepant with primary CMB measurements than previous DESI LRG cross Planck CMB lensing results. From the cross-correlation with ACT we obtain $S_8 = 0.790^{+0.024}_{-0.027}$, while when jointly analyzing Planck and ACT we find $S_8 = 0.775^{+0.019}_{-0.022}$ from our data alone and $\sigma_8 = 0.772^{+0.020}_{-0.023}$ with the addition of BAO data. These constraints are consistent with the latest Planck primary CMB analyses at the $\simeq 1.6-2.2\sigma$ level, and are in excellent agreement with galaxy lensing surveys., Comment: 60 pages, 26 figures, comments welcome

We present a high-significance cross-correlation of CMB lensing maps from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) with spectroscopically calibrated luminous red galaxies (LRGs) from the Dark Energy Spectroscopic Instrument (DESI). We detect this cross-correlation at a significance of 38$\sigma$; combining our measurement with the Planck Public Release 4 (PR4) lensing map, we detect the cross-correlation at 50$\sigma$. Fitting this jointly with the galaxy auto-correlation power spectrum to break the galaxy bias degeneracy with $\sigma_8$, we perform a tomographic analysis in four LRG redshift bins spanning $0.4 \le z \le 1.0$ to constrain the amplitude of matter density fluctuations through the parameter combination $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4}$. Prior to unblinding, we confirm with extragalactic simulations that foreground biases are negligible and carry out a comprehensive suite of null and consistency tests. Using a hybrid effective field theory (HEFT) model that allows scales as small as $k_{\rm max}=0.6$ $h/{\rm Mpc}$, we obtain a 3.3% constraint on $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4} = 0.792^{+0.024}_{-0.028}$ from ACT data, as well as constraints on $S_8^\times(z)$ that probe structure formation over cosmic time. Our result is consistent with the early-universe extrapolation from primary CMB anisotropies measured by Planck PR4 within 1.2$\sigma$. Jointly fitting ACT and Planck lensing cross-correlations we obtain a 2.7% constraint of $S_8^\times = 0.776^{+0.019}_{-0.021}$, which is consistent with the Planck early-universe extrapolation within 2.1$\sigma$, with the lowest redshift bin showing the largest difference in mean. The latter may motivate further CMB lensing tomography analyses at $z<0.6$ to assess the impact of potential systematics or the consistency of the $\Lambda$CDM model over cosmic time., Comment: Prepared for submission to JCAP (47 pages, 13 figures)

We present a new perturbative full-shape analysis of BOSS galaxy clustering data, including the full combination of the galaxy power spectrum and bispectrum multipoles, baryon acoustic oscillations, and cross-correlations with the gravitational lensing of cosmic microwave background measured from \textit{Planck}. Assuming the $\Lambda$CDM model, we constrain the matter density fraction $\Omega_m = 0.3154\pm 0.0089$, the Hubble constant $H_0=68.34\pm 0.77\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, and the mass fluctuation amplitude $\sigma_8=0.686\pm 0.027$ (equivalent to $S_8 = 0.704\pm 0.031$). Cosmic structure at low redshifts appears suppressed with respect to the Planck $\Lambda$CDM concordance model at $4.5\sigma$. We explore whether this tension can be explained by the recent DESI preference for dynamical dark energy (DDE): the BOSS data combine with DESI BAO and PantheonPlus supernovae competitively compared to the CMB, yielding no preference for DDE, but the same $\sim 10\%$ suppression of structure, with dark energy being consistent with a cosmological constant at 68\% CL. Our results suggest that either the data contains residual systematics, or more model-building efforts may be required to restore cosmological concordance., Comment: 9 pages, 2 figures, to be submitted to PRL

We study the bispectrum in Lagrangian perturbation theory. Extending past results for the power spectrum, we describe a method to efficiently compute the bispectrum in LPT, focusing on the Zeldovich approximation, in which contributions due to linear displacements are captured to all orders in a manifestly infrared (IR) safe way. We then isolate the effects of these linear displacements on oscillatory components of the power spectrum like baryon acoustic oscillations or inflationary primordial features and show that the Eulerian perturbation theory (EPT) prescription wherein their effects are resummed by a Gaussian damping of the oscillations arise as a saddle-point approximation of our calculation. These two methods of IR resummation are in excellent agreement at 1-loop in the bispectrum. At tree level, resummed EPT does less well to capture the nonlinear damping of the oscillations, and the LPT calculation does not require an artificial split of the power spectrum into smooth and oscillatory components, making the latter particularly useful for modeling exotic features. We finish by extending our analysis of IR resummation in LPT to N-point functions of arbitrary order., Comment: 47 pages, 13 figures, updated to match version accepted by JCAP

We derive new constraints on the $E_G$ statistic as a test of gravity, combining the CMB lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with SDSS BOSS CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the $E_G$ statistic can distinguish $\Lambda$CDM from alternative models of gravity. We find $E_G= 0.31^{+0.06}_{-0.05}$ for ACT and CMASS data at 68.28\% confidence level, and $E_G = 0.49^{+0.14}_{-0.11}$ for ACT and LOWZ. Systematic errors are estimated to be 3\% and 4\% respectively. Including CMB lensing information from Planck PR4 results in $E_G = 0.34^{+0.05}_{-0.05}$ with CMASS and $E_G= 0.43^{+0.11}_{-0.09}$ with LOWZ. These are consistent with predictions for the $\Lambda$CDM model that best fits the Planck CMB anisotropy and SDSS BOSS BAO, where $E_G^{\rm GR} (z_{\rm eff} = 0.555) = 0.401\pm 0.005$ for CMB lensing combined with CMASS and $E_G^{\rm GR} (z_{\rm eff} = 0.316) = 0.452\pm0.005$ combined with LOWZ. We also find $E_G$ to be scale independent, with PTE $>5\%$, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements., Comment: 33 pages, 21 figures, prepared for submission to PRD

Massive particles leave imprints on primordial non-Gaussianity via couplings to the inflaton, even despite their exponential dilution during inflation: practically, the Universe acts as a Cosmological Collider. We present the first dedicated search for spin-zero particles using BOSS redshift-space galaxy power spectrum and bispectrum multipoles, as well as Planck CMB non-Gaussianity data. We demonstrate that some Cosmological Collider models are well approximated by the standard equilateral and orthogonal parametrization; assuming negligible inflaton self-interactions, this facilitates us translating Planck non-Gaussianity constraints into bounds on Collider models. Many models have signatures that are not degenerate with equilateral and orthogonal non-Gaussianity and thus require dedicated searches. Here, we constrain such models using BOSS three-dimensional redshift-space galaxy clustering data, focusing on spin-zero particles in the principal series and constraining their couplings to the inflaton at varying speed and mass, marginalizing over the unknown inflaton self-interactions. This is made possible through an improvement in Cosmological Bootstrap techniques and the combination of perturbation theory and halo occupation distribution models for galaxy clustering. Our work sets the standard for inflationary spectroscopy with cosmological observations, providing the ultimate link between physics on the largest and smallest scales., Comment: 35 pages, 15 figures, 6 tables

This paper provides a comprehensive overview of how fitting of Baryon Acoustic Oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument's (DESI) 2024 results using its DR1 dataset, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential biases in measurements of the isotropic ($\alpha_{\mathrm{iso}}$) and anisotropic ($\alpha_{\mathrm{ap}}$) BAO distance scales, and how these can be effectively removed with an appropriate choice of reconstruction algorithm. We then demonstrate how theory leads to a clear choice for how to model the BAO and develop, implement and validate a new model for the remaining smooth-broadband (i.e., without BAO) component of the galaxy clustering. Finally, we explore the impact of all remaining modelling choices on the BAO constraints from DESI using a suite of high-precision simulations, arriving at a set of best-practices for DESI BAO fits, and an associated theory and modelling systematic error. Overall, our results demonstrate the remarkable robustness of the BAO to all our modelling choices and motivate a combined theory and modelling systematic error contribution to the post-reconstruction DESI BAO measurements of no more than $0.1\%$ ($0.2\%$) for its isotropic (anisotropic) distance measurements. We expect the theory and best-practices laid out to here to be applicable to other BAO experiments in the era of DESI and beyond., Comment: 30 pages, 18 figures, 1 table, updated to match version accepted by MNRAS

The $E_G$ statistic is a discriminating probe of gravity developed to test the prediction of general relativity (GR) for the relation between gravitational potential and clustering on the largest scales in the observable universe. We present a novel high-precision estimator for the $E_G$ statistic using CMB lensing and galaxy clustering correlations that carefully matches the effective redshifts across the different measurement components to minimize corrections. A suite of detailed tests is performed to characterize the estimator's accuracy, its sensitivity to assumptions and analysis choices and the non-Gaussianity of the estimator's uncertainty is characterized. After finalization of the estimator, it is applied to $\textit{Planck}$ CMB lensing and SDSS CMASS and LOWZ galaxy data. We report the first harmonic space measurement of $E_G$ using the LOWZ sample and CMB lensing and also updated constraints using the final CMASS sample and the latest $\textit{Planck}$ CMB lensing map. We find $E_G^{Planck+CMASS} = 0.36^{+0.06}_{-0.05}$ (68.27%) and $E_G^{\rm \textit{Planck}+LOWZ} = 0.40^{+0.11}_{-0.09} $ (68.27%), with additional subdominant systematic error budget estimates of 2% and 3% respectively. Using $\Omega_{\rm m,0}$ constraints from $\textit{Planck}$ and SDSS BAO observations, $\Lambda$CDM-GR predicts $E_G^{\rm GR} (z = 0.555) = 0.401 \pm 0.005$ and $E_G^{\rm GR} (z = 0.316) = 0.452 \pm 0.005$ at the effective redshifts of the CMASS and LOWZ based measurements. We report the measurement to be in good statistical agreement with the $\Lambda$CDM-GR prediction, and report that the measurement is also consistent with the more general GR prediction of scale-independence for $E_G$. This work provides a carefully constructed and calibrated statistic with which $E_G$ measurements can be confidently and accurately obtained with upcoming survey data., Comment: 36 pages, 20 figures, Accepted for publication in PRD

We formulate the Lagrangian perturbation theory of galaxy intrinsic alignments and compute the resulting auto and cross power spectra of galaxy shapes, densities and matter to 1-loop order. Our model represents a consistent effective-theory description of galaxy shape including the resummation of long-wavelength displacements which damp baryon acoustic oscillations, and includes one linear, three quadratic and two cubic dimensionless bias coefficients at this order, along with counterterms and stochastic contributions whose structure we derive. We compare this Lagrangian model against the three-dimensional helicity spectra of halo shapes measured in N-body simulations by Akitsu et al (2023) and find excellent agreement on perturbative scales while testing a number of more restrictive bias parametrizations. The calculations presented are immediately relevant to analyses of both cosmic shear surveys and spectroscopic shape measurements, and we make a fast FFTLog-based code spinosaurus publicly available with this publication., Comment: 41 pages, 8 figures, updated to match version accepted in JCAP. Fast Python code to predict IA correlations computed in the text can be found at https://github.com/sfschen/spinosaurus

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological $N$-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis. We also provide an optimized package for computing the power spectra and other two-point statistics of an arbitrary galaxy catalog as well as a pipeline for obtaining CV-reduced measurements on any of the AbacusSummit cubic box outputs. We make our scripts publicly available and report a speed improvement of $\sim$10 for a grid size of $N_{\rm mesh} = 256^3$ compared with \texttt{nbodykit}., Comment: 17 pages, 9 figures, public package (for power spectrum and control variates estimation)

In addition to the intrinsic clustering of galaxies themselves, the spatial distribution of galaxies observed in surveys is modulated by the presence of weak lensing due to matter in the foreground. This effect, known as magnification bias, is a significant contaminant to analyses of galaxy-lensing cross-correlations and must be carefully modelled. We present a method to estimate the magnification bias in spectroscopically confirmed galaxy samples based on finite differences of galaxy catalogues while marginalizing over errors due to finite step size. We use our estimator to measure the magnification biases of the CMASS and LOWZ samples in the SDSS BOSS galaxy survey, analytically taking into account the dependence on galaxy shape for fiber and PSF magnitudes, finding $\alpha_{\rm CMASS} = 2.71 \pm 0.02$ and $\alpha_{\rm LOWZ} = 2.45 \pm 0.02$ and quantify modelling uncertainties in these measurements. Finally, we quantify the redshift evolution of the magnification bias within the CMASS and LOWZ samples, finding a difference of up to a factor of three between the lower and upper redshift bounds for the former. We discuss how to account for this evolution in modelling and its interaction with commonly applied redshift-dependent weights. Our method should be readily-applicable to upcoming surveys and we make our code publicly available as part of this work., Comment: 14 pages, 9 figures, 3 tables, Accepted in MNRAS

We present the Aemulus $\nu$ simulations: a suite of 150 $(1.05 h^{-1}\rm Gpc)^3$ $N$-body simulations with a mass resolution of $3.51\times 10^{10} \frac{\Omega_{cb}}{0.3} ~ h^{-1} M_{\odot}$ in a $w\nu$CDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in $\sigma_8$ to facilitate investigations of tension between large scale structure and cosmic microwave background cosmological probes. Neutrinos are treated as a second particle species to ensure accuracy to $0.5\, \rm eV$, the maximum neutrino mass that we have simulated. By employing Zel'dovich control variates, we increase the effective volume of our simulations by factors of $10-10^5$ depending on the statistic in question. As a first application of these simulations, we build new hybrid effective field theory and matter power spectrum surrogate models, demonstrating that they achieve $\le 1\%$ accuracy for $k\le 1\, h\,\rm Mpc^{-1}$ and $0\le z \le 3$, and $\le 2\%$ accuracy for $k\le 4\, h\,\rm Mpc^{-1}$ for the matter power spectrum. We publicly release the trained surrogate models, and estimates of the surrogate model errors in the hope that they will be broadly applicable to a range of cosmological analyses for many years to come., Comment: 37 pages, 15 figures, matching version accepted by JCAP

The growth of large-scale structure, as revealed in the anisotropic of clustering of galaxies in the low redshift Universe, provides a stringent test of our cosmological model. The strongest current constraints come from the BOSS and eBOSS surveys, with uncertainties on the amplitude of clustering of less than 10 per cent. A number of different approaches have been taken to fitting this signal, leading to apparently discrepant conclusions about the amplitude of fluctuations at late times. We compare in some detail two of the leading approaches, one based on a fitting a template cosmology whose amplitude and length scales are allowed to float with one based on a more traditional forward modeling approach, when fitting to the BOSS DR12 data. Holding the input data, scale cuts, window functions and modeling framework fixed we are able to isolate the cause of the differences and discuss the implications for future surveys.

Numerical simulations in cosmology require trade-offs between volume, resolution and run-time that limit the volume of the Universe that can be simulated, leading to sample variance in predictions of ensemble-average quantities such as the power spectrum or correlation function(s). Sample variance is particularly acute at large scales, which is also where analytic techniques can be highly reliable. This provides an opportunity to combine analytic and numerical techniques in a principled way to improve the dynamic range and reliability of predictions for clustering statistics. In this paper we extend the technique of Zel'dovich control variates, previously demonstrated for 2-point functions in real space, to reduce the sample variance in measurements of 2-point statistics of biased tracers in redshift space. We demonstrate that with this technique, we can reduce the sample variance of these statistics down to their shot-noise limit out to $k \sim 0.2\, h\rm Mpc^{-1}$. This allows a better matching with perturbative models and improved predictions for the clustering of e.g.~quasars, galaxies and neutral Hydrogen measured in spectroscopic redshift surveys at very modest computational expense. We discuss the implementation of ZCV, give some examples and provide forecasts for the efficacy of the method under various conditions., Comment: 17 pages main text, 5 pages of appendices, 9 figures

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization., Comment: Updated to reflect published version

The optical depth $\tau$ is the least well determined parameter in the standard model of cosmology, and one whose precise value is important for both understanding reionization and for inferring fundamental physics from cosmological measurements. We forecast how well future epoch of reionization experiments could constraint $\tau$ using a symmetries-based bias expansion that highlights the special role played by anisotropies in the power spectrum on large scales. Given a parametric model for the ionization evolution inspired by the physical behavior of more detailed reionization simulations, we find that future 21cm experiments could place tight constraints on the timing and duration of reionization and hence constraints on $\tau$ that are competitive with proposed, space-based CMB missions provided they can measure $k\approx 0.1\,h\,\text{Mpc}^{-1}$ with a clean foreground wedge across redshifts spanning the most active periods of reionization, corresponding to ionization fractions $0.2 \lesssim x \lesssim 0.8$. Significantly improving upon existing CMB-based measurements with next-generation 21cm surveys would require substantially longer observations ($\sim5$ years) than standard $\mathcal{O}(1000 \,\,\text{hour})$ integration times. Precise measurements of smaller scales will not improve constraints on $\tau$ until a better understanding of the astrophysics of reionization is achieved. In the presence of noise and foregrounds even future 21cm experiments will struggle to constrain $\tau$ if the ionization evolution deviates significantly from simple parametric forms., Comment: published in JCAP

We present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain $\Omega_\mathrm{m} = 0.303\pm 0.008$, $H_0 = 69.21\pm 0.78$ and $\sigma_8 = 0.743\pm 0.043$. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to $\Omega_m = 0.300 \pm 0.008$, $H_0 = 69.21 \pm 0.77$ and $\sigma_8 = 0.707 \pm 0.035$, in tension with CMB and cosmic shear constraints. The combination of $\Omega_m$ and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $\sigma_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($z\simeq 0.38$) and at large angular scales, which show a $20\%$ deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions., Comment: 46 pages, 15 figures, updated to match version accepted by JCAP

We make use of neural networks to accelerate the calculation of power spectra required for the analysis of galaxy clustering and weak gravitational lensing data. For modern perturbation theory codes, evaluation time for a single cosmology and redshift can take on the order of two seconds. In combination with the comparable time required to compute linear predictions using a Boltzmann solver, these calculations are the bottleneck for many contemporary large-scale structure analyses. In this work, we construct neural network-based surrogate models for Lagrangian perturbation theory (LPT) predictions of matter power spectra, real and redshift space galaxy power spectra, and galaxy--matter cross power spectra that attain $\sim 0.1\%$ (at one sigma) accuracy over a broad range of scales in a $w$CDM parameter space. The neural network surrogates can be evaluated in approximately one millisecond, a factor of 1000 times faster than the full Boltzmann code and LPT computations. In a simulated full-shape redshift space galaxy power spectrum analysis, we demonstrate that the posteriors obtained using our surrogates are accurate compared to those obtained using the full LPT model. We make our surrogate models public at https://github.com/sfschen/EmulateLSS, so that others may take advantage of the speed gains they provide to enable rapid iteration on analysis settings, something that is essential in complex contemporary large-scale structure analyses., Comment: 12 pages, 4 figures, models available at https://github.com/sfschen/EmulateLSS

We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use Halo Occupation Distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with datasets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of hybrid EFT methods to characterize models of the galaxy--halo connection at the field level, revealing new connections between once-disparate modelling frameworks., Comment: 16 pages, 10 figures. Revised version accepted to MNRAS. Revisions include a new appendix on the covariance of the field-level bias estimator

We use luminous red galaxies selected from the imaging surveys that are being used for targeting by the Dark Energy Spectroscopic Instrument (DESI) in combination with CMB lensing maps from the Planck collaboration to probe the amplitude of large-scale structure over $0.4\le z\le 1$. Our galaxy sample, with an angular number density of approximately $500\,\mathrm{deg}^{-2}$ over 18,000 sq.deg., is divided into 4 tomographic bins by photometric redshift and the redshift distributions are calibrated using spectroscopy from DESI. We fit the galaxy autospectra and galaxy-convergence cross-spectra using models based on cosmological perturbation theory, restricting to large scales that are expected to be well described by such models. Within the context of $\Lambda$CDM, combining all 4 samples and using priors on the background cosmology from supernova and baryon acoustic oscillation measurements, we find $S_8=\sigma_8(\Omega_m/0.3)^{0.5}=0.73\pm 0.03$. This result is lower than the prediction of the $\Lambda$CDM model conditioned on the Planck data. Our data prefer a slower growth of structure at low redshift than the model predictions, though at only modest significance., Comment: 44 pages, 16 figures. Matches version accepted by journal: more details on analysis, updated references, link to data added

We present a new method for consistent, joint analysis of the pre- and post-reconstruction two-point functions of the BOSS survey. The post-reconstruction correlation function is used to accurately measure the distance-redshift relation and expansion history, while the pre-reconstruction power spectrum multipoles constrain the broad-band shape and the rate-of-growth of large-scale structure. Our technique uses Lagrangian perturbation theory to self-consistently work at the level of two-point functions, i.e.\ directly with the measured data, without approximating the constraints with summary statistics normalized by the drag scale. Combining galaxies across the full redshift range and both hemispheres we constrain $\Omega_m=0.303 \pm 0.0082$, $H_0=69.23 \pm 0.77$ and $\sigma_8=0.733 \pm 0.047$ within the context of $\Lambda$CDM. These constraints are in good agreement both with the Planck primary CMB anisotropy data and recent cosmic shear surveys., Comment: 39 pages, 10 figures, updated to correct minor typos and match version in JCAP. Note updated results from v1 fixing bug in code

The Ly$\alpha$ forest provides one of the best means of mapping large-scale structure at high redshift, including our tightest constraint on the distance-redshift relation before cosmic noon. We describe how the large-scale correlations in the Ly$\alpha$ forest can be understood as an expansion in cumulants of the optical depth field, which itself can be related to the density field by a bias expansion. This provides a direct connection between the observable and the statistics of the matter fluctuations which can be computed in a systematic manner. We discuss the way in which complex, small-scale physics enters the predictions, the origin of the much-discussed velocity bias and the `renormalization' of the large-scale bias coefficients. Our calculations are within the context of perturbation theory, but we also make contact with earlier work using the peak-background split. Using the structure of the equations of motion we demonstrate, to all orders in perturbation theory, that the large-scale flux power spectrum becomes the linear spectrum times the square of a quadratic in the cosine of the angle to the line of sight. Unlike the case of galaxies, both the isotropic and anisotropic pieces receive contributions from small-scale physics., Comment: 28 pages, 4 figures, updated to match version accepted by JCAP

We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from $N$-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of $N$-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these nonlinear spectra from redshifts $z = 0$ to $z=2$. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h {\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift $z\simeq 0.4$ to scales of $k_{\rm max} \approx 0.6\, h\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy--galaxy lensing analyses with data drawn from an independent $N$-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming datasets from galaxy surveys., Comment: 19 pages, 17 figures. Updated to reflect the journal version. Code available at https://github.com/kokron/anzu

We present the one-loop 2-point function of biased tracers in redshift space computed with Lagrangian perturbation theory, including a full resummation of both long-wavelength (infrared) displacements and associated velocities. The resulting model accurately predicts the power spectrum and correlation function of halos and mock galaxies from two different sets of N-body simulations at the percent level for quasi-linear scales, including the damping of the baryon acoustic oscillation signal due to the bulk motions of galaxies. We compare this full resummation with other, approximate, techniques including the moment expansion and Gaussian streaming model. We discuss infrared resummation in detail and compare our Lagrangian formulation with the Eulerian theory augmented by an infrared resummation based on splitting the input power spectrum into "wiggle" and "no-wiggle" components. We show that our model is able to recover unbiased cosmological parameters in mock data encompassing a volume much larger than what will be available to future galaxy surveys. We demonstrate how to efficiently compute the resulting expressions numerically, making available a fast Python code capable of rapidly computing these statistics in both configuration and Fourier space., Comment: 37 pages, 10 figures, updated to match version accepted by JCAP

We study the ability of perturbative models with effective field theory contributions and infra-red resummation to model the redshift space clustering of biased tracers in models where the linear power spectrum has "features" -- either imprinted during inflation or induced by non-standard expansion histories. We show that both Eulerian and Lagrangian perturbation theory are capable of reproducing the Fourier space two-point functions of halos up to the non-linear scale from a suite of $4096^3$ particle N-body simulations. This is the first demonstration that perturbative models can accurately fit the redshift-space clustering of biased tracers in N-body simulations of such theories. By comparing different theoretical models and IR resummation schemes we assess the current theoretical uncertainty in predicting power spectra for models with features. Our results suggest that future surveys will be able to detect or tightly constrain features in the primordial spectrum below the one percent level across a wide range of scales., Comment: 22 pages, 8 figures, updated to match version accepted by JCAP

The peculiar velocities of biased tracers of the cosmic density field contain important information about the growth of large scale structure and generate anisotropy in the observed clustering of galaxies. Using N-body data, we show that velocity expansions for halo redshift-space power spectra are converged at the percent-level at perturbative scales for most line-of-sight angles $\mu$ when the first three pairwise velocity moments are included, and that the third moment is well-approximated by a counterterm-like contribution. We compute these pairwise-velocity statistics in Fourier space using both Eulerian and Lagrangian one-loop perturbation theory using a cubic bias scheme and a complete set of counterterms and stochastic contributions. We compare the models and show that our models fit both real-space velocity statistics and redshift-space power spectra for both halos and a mock sample of galaxies at sub-percent level on perturbative scales using consistent sets of parameters, making them appealing choices for the upcoming era of spectroscopic, peculiar-velocity and kSZ surveys., Comment: 63 pages, 11 figures, updated to match version accepted by JCAP

We explore the deep ultraviolet (that is, short-distance) limit of the power spectrum (PS) and of the correlation function of a cold dark matter dominated Universe. While for large scales the PS can be written as a double series expansion, in powers of the linear PS and of the wavenumber $k$, we show that, in the opposite limit, it can be expressed via an expansion in powers of the form $1/k^{d+2n}$, where $d$ is the number of spatial dimensions, and $n$ is a non negative integer. The coefficients of the terms of the expansion are nonperturbative in the linear PS, and can be interpreted in terms of the probability density function for the displacement field, evaluated around specific configurations of the latter, that we identify. In the case of the Zel'dovich dynamics, these coefficients can be determined analytically, whereas for the exact dynamics they can be treated as fit, or nuisance, parameters. We confirm our findings with numerical simulations and discuss the necessary steps to match our results to those obtained for larger scales and to actual measurements., Comment: 27 pages, 12 figures

We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify halos and subhalos, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modeling a more limited range in scale., Comment: 10 pages, 7 figures, updated to reflect version to be published in MNRAS

Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. In this paper we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the {\tt DarkSky} simulation \cite{skillman14}, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical results for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation., Comment: 46 pages, 15 figures, updated to match version accepted by JCAP. Corrected a typo in Equation 4.11

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space., Comment: Science white paper submitted to the Astro2020 Decadal Survey

We explore Lagrangian perturbation theory (LPT) for biased tracers in the presence of two fluids, focusing on the case of cold dark matter (CDM) and baryons. The presence of two fluids induces corrections to the Lagrangian bias expansion and tracer advection, both of which we formulate as expansions in the three linear modes of the Lagrangian equations of motion. We compute the linear-order two-fluid corrections in the Zeldovich approximation, finding that modifications to the bias expansion and tracer advection both enter as percent-level corrections over a large range of wavenumbers at low redshift and draw parallels with the Eulerian formalism. We then discuss nonlinear corrections in the two-fluid picture, and calculate contributions from the relative velocity effect ($\propto v_r^2$) at one loop order. Finally, we conduct an exploratory Fisher analysis to assess the impact of two-fluid corrections on baryon acoustic oscillations (BAO) measurements, finding that while modest values of the relative bias parameters can introduce systematic biases in the measured BAO scale of up to $0.5\, \sigma$, fitting for these effects as additional parameters increases the error bar by less than $30\%$ across a wide range of bias values., Comment: 31 pages, 11 figures. Revised to match version accepted by JCAP

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics., Comment: This white-paper was prepared under auspices of the Department of Energy Cosmic Visions Dark Energy program; v1: initial release in Oct 2018; v2: Jul 2019 update: strawman design matches PUMA Decadal Submission, updated forecasting, new figures and science section; v3: fixed reference, missing author in metadata

We study synergies between three promising methods to measure $22$ and suffer from interloper problems even for spectroscopic surveys. Intensity mapping of the 21-cm signal can cover large volumes with exquisite fidelity, but is limited both by loss of information to foreground cleaning and by lack of knowledge of the mean signal. Cosmic microwave background (CMB) lensing is theoretically very clean, but ultimately measures just the projected variations in density. We find that cross-correlation between optical and radio can significantly improve the measurement of growth rate. Combining these with the CMB provides a promising avenue to detecting modified gravity at high redshifts, in particular by independently probing the Weyl and Newtonian potentials and by strengthening control of systematics. We find that cross-correlating a Stage {\sc ii} 21-cm survey with DESI quasars with a reasonable brightness temperature prior could enable measurements of the growth rate $f\sigma_8$ at sub 3\% and sub 8\% levels at $z = 3, 4$, representing a factor of 4 and 8 improvement over constraints obtainable from DESI quasars alone. Similarly, cross-correlating 21-cm data with a futuristic LBG survey to $m_{UV}<24.5$ over 1000 square degrees will make possible $f\sigma_8$ measurements at close to 1\% at $z = 3$ and 3\% at $z = 4$, and improve similar constraints at $z = 5$ by close to a factor of 3 to sub-10\% precision. Combining the above with CMB lensing from a Stage 4 CMB survey and LSST data can additionally constrain the gravitational slip $\gamma$ parameter to similar precision at these redshifts, enabling us to test the predictions of general relativity at large scales., Comment: Updated to match version accepted by JCAP. Fixed typo in Equation 3.2, updated 21-cm experiment specifications and added Figure 3

We present final statistics from a survey for intervening MgII absorption towards 100 quasars with emission redshifts between $z=3.55$ and $z=7.08$. Using infrared spectra from Magellan/FIRE, we detect 279 cosmological MgII absorbers, and confirm that the incidence rate of $W_r>0.3 \AA$ MgII absorption per comoving path length does not evolve measurably between $z=0.25$ and $z=7$. This is consistent with our detection of seven new MgII systems at $z>6$, a redshift range that was not covered in prior searches. Restricting to relatively strong MgII systems ($W_r>1$\AA), there is significant evidence for redshift evolution. These systems roughly double in number density between $z=0$ and $z=2$-$3$, but decline by an order of magnitude from this peak by $z\sim 6$. This evolution mirrors that of the global star formation rate density, which could reflect a connection between star formation feedback and strong MgII absorbers. We compared our results to the Illustris cosmological simulation at $z=2$-$4$ by assigning absorption to catalogued dark-matter halos and by direct extraction of spectra from the simulation volume. To reproduce our results using the halo catalogs, we require circumgalactic (CGM) MgII envelopes within halos of progressively smaller mass at earlier times. This occurs naturally if we define the lower integration cutoff using SFR rather than mass. MgII profiles calculated directly from the Illustris volume yield far too few strong absorbers. We argue that this arises from unresolved phase space structure of CGM gas, particularly from turbulent velocities on sub-mesh scales. The presence of CGM MgII at $z>6$-- just $\sim 250$ Myr after the reionization redshift implied by Planck--suggests that enrichment of intra-halo gas may have begun before the presumed host galaxies' stellar populations were mature and dynamically relaxed. [abridged], Comment: 23 Pages, 17 figures, accepted for publication in ApJ

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease)., Comment: 6 pages, 4 figures, 1 table. To appear in ApJL