The growth of large-scale structure, as revealed in the anisotropic of clustering of galaxies in the low redshift Universe at z < 2, provides a stringent test of our cosmological model. The strongest current constraints come from the BOSS and eBOSS surveys, with uncertainties on σ 8, the amplitude of clustering on an 8 h -1Mpc scale, of less than 10 per cent. A number of different approaches have been taken to fitting this signal, leading to discrepancies of up to 1σ in the measurements of the amplitude of fluctuations at late times. We compare in some detail two of the leading approaches, one based on fitting a template cosmology whose amplitude and length scales are allowed to float with one based on varying the underlying parameters of a cosmological model directly, 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.

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 DarkSky simulation [1], 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.

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 μ 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.

Abstract We present a new method for consistent, joint analysis of the pre- and post-reconstruction galaxy 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 Ωm = 0.303 ± 0.0082, H 0 = 69.23 ± 0.77 and σ8 = 0.733 ± 0.047 within the context of ΛCDM. These constraints are consistent both with the Planck primary CMB anisotropy data and recent cosmic shear surveys.

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.

The optical depth τ 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 τ 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 21 cm experiments could place tight constraints on the timing and duration of reionization and hence constraints on τ that are competitive with proposed, space-based CMB missions provided they can measure k ≈ 0.1 h Mpc-1 with a clean foreground wedge across redshifts spanning the most active periods of reionization, corresponding to ionization fractions 0.2 ≲ x ≲ 0.8. Significantly improving upon existing CMB-based measurements with next-generation 21 cm surveys would require substantially longer observations (∼ 5 years) than standard O(1000 hour) integration times. Precise measurements of smaller scales will not improve constraints on τ until a better understanding of the astrophysics of reionization is achieved. In the presence of noise and foregrounds even future 21 cm experiments will struggle to constrain τ if the ionization evolution deviates significantly from simple parametric forms.

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 40963 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.

Abstract: 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.

The Lyα 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α 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.

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 (vr2) 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 σ, fitting for these effects as additional parameters increases the error bar by less than 30% across a wide range of bias values.