The most common statistic used to analyze large-scale structure surveys is the correlation function, or power spectrum. Here, we show how `slicing' the correlation function on local density brings sensitivity to interesting non-Gaussian features in the large-scale structure, such as the expansion or contraction of baryon acoustic oscillations (BAO) according to the local density. The sliced correlation function measures the large-scale flows that smear out the BAO, instead of just correcting them as reconstruction algorithms do. Thus, we expect the sliced correlation function to be useful in constraining the growth factor, and modified gravity theories that involve the local density. Out of the studied cases, we find that the run of the BAO peak location with density is best revealed when slicing on a $\sim 40$ Mpc/$h$ filtered density. But slicing on a $\sim100$ Mpc/$h$ filtered density may be most useful in distinguishing between underdense and overdense regions, whose BAO peaks are separated by a substantial $\sim 5$ Mpc/$h$ at $z=0$. We also introduce `curtain plots' showing how local densities drive particle motions toward or away from each other over the course of an $N$-body simulation.

Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSS-IV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [OII] emitters, which dominate optical ELG selections at $z\simeq 1$. Model [OII] emitters at $0.5<z<1.5$ are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [OII] emitters are in reasonable agreement with observations. The selected [OII] emitters are hosted by haloes with $M_{\rm halo}\geq 10^{10.3}h^{-1}{\rm M}_{\odot}$, with ~90% of them being central star-forming galaxies. The predicted mean halo occupation distributions of [OII] emitters has a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, $\langle N \rangle_{\left[OII\right]\, cen}$, being far from the canonical step function. The $\langle N \rangle_{\left[OII\right]\, cen}$ can be described as the sum of an asymmetric Gaussian for disks and a step function for spheroids, which plateaus below unity. The model [OII] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At $z\sim 1$, a comparison with observed g-band selected galaxy, which are expected to be dominated by [OII] emitters, indicates that our model produces too few [OII] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.

The standard Halo Occupation Distribution (HOD) models were originally developed based on results from semi-analytic and hydrodynamical galaxy formation models. Those models have since progressed, in particular to include AGN feedback to match the galaxy luminosity function in a universe with the observed baryon fraction. AGN feedback affects the relationship between galaxy stellar mass and luminosity, in particular making the relationship non-monotonic. For matched number density samples, galaxies in luminosity-threshold samples occupy a different range of halo masses from those in stellar-mass-threshold samples. We find that the shapes of the HODs of luminosity-threshold samples are slightly more complicated in semi-analytic galaxy formation models that include AGN feedback than are assumed by standard HOD models. We also find that subhalo abundance matching (SHAM) does not preserve these non-standard shapes. We show that catalogues created using SHAM and the semi-analytic model Galform that have the same large-scale 2-point clustering by construction have different void probability functions (VPFs) in both real and redshift space. We find that these differences arise from the different HOD shapes, as opposed to assembly bias, which indicates that the VPF could be used to test the suitability of an HOD model with real data.

In view of future high-precision large-scale structure surveys, it is important to quantify the percent and subpercent level effects in cosmological $N$-body simulations from which theoretical predictions are drawn. One such effect involves deciding whether to zero all modes above the one-dimensional Nyquist frequency, the so-called "corner" modes, in the initial conditions. We investigate this effect by comparing power spectra, density distribution functions, halo mass functions, and halo profiles in simulations with and without these modes. For a simulation with a mass resolution of $m_p \sim 10^{11}\,h^{-1}\,M_{\odot}$, we find that at $z>6$, the difference in the matter power spectrum is large at wavenumbers above $\sim 80$\% of $k_{\rm{Ny}}$, reducing to below 2\% at all scales by $z\sim 3$. Including corner modes results in a better match between low- and high-resolution simulations at wavenumbers around the Nyquist frequency of the low-resolution simulation, but the effect of the corner modes is smaller than the effect of particle discreteness. The differences in mass functions are 3\% for the smallest halos at $z=6$ for the $m_p \sim 10^{11}\,h^{-1}\,M_{\odot}$ simulation, but we find no significant difference in the stacked profiles of well-resolved halos at $z \leq 6$. Thus removing power at $|\mathbf{k}|>k_{\rm{Ny}}$ in the initial conditions of cosmological simulations has a small effect on small scales and high redshifts, typically below a few percent.

When it comes to searches for extensions to general relativity, large efforts are being dedicated to accurate predictions for the power spectrum of density perturbations. While this observable is known to be sensitive to the gravitational theory, its efficiency as a diagnostic for gravity is significantly reduced when Solar System constraints are strictly adhered to. We show that this problem can be overcome by studying weigthed density fields. We propose a transformation of the density field for which the impact of modified gravity on the power spectrum can be increased by more than a factor of three. The signal is not only amplified, but the modified gravity features are shifted to larger scales which are less affected by baryonic physics. Furthermore, the overall signal-to-noise increases, which in principle makes identifying signatures of modified gravity with future galaxy surveys more feasible. While our analysis is focused on modified gravity, the technique can be applied to other problems in cosmology, such as the detection of neutrinos, the effects of baryons or baryon acoustic oscillations.

The Gaussianization transform has been proposed as a method to remove the issues of scale-dependent galaxy bias and nonlinearity from galaxy clustering statistics, but these benefits have yet to be thoroughly tested for realistic galaxy samples. In this paper, we test the effectiveness of the Gaussianization transform for different galaxy types by applying it to realistic simulated blue and red galaxy samples. We show that in real space, the shapes of the Gaussianized power spectra of both red and blue galaxies agree with that of the underlying dark matter, with the initial power spectrum, and with each other to smaller scales than do the statistics of the usual (untransformed) density field. However, we find that the agreement in the Gaussianized statistics breaks down in redshift space. We attribute this to the fact that red and blue galaxies exhibit very different fingers of god in redshift space. After applying a finger-of-god compression, the agreement on small scales between the Gaussianized power spectra is restored. We also compare the Gaussianization transform to the clipped galaxy density field and find that while both methods are effective in real space, they have more complicated behaviour in redshift space. Overall, we find that Gaussianization can be useful in recovering the shape of the underlying dark matter power spectrum to k ~ 0.5 h/Mpc and of the initial power spectrum to k ~ 0.4 h/Mpc in certain cases at z = 0.

Accurate modeling of nonlinearities in the galaxy bispectrum, the Fourier transform of the galaxy three-point correlation function, is essential to fully exploit it as a cosmological probe. In this paper, we present numerical and theoretical challenges in modeling the nonlinear bispectrum. First, we test the robustness of the matter bispectrum measured from N-body simulations using different initial conditions generators. We run a suite of N-body simulations using the Zel'dovich approximation and second-order Lagrangian perturbation theory (2LPT) at different starting redshifts, and find that transients from initial decaying modes systematically reduce the nonlinearities in the matter bispectrum. To achieve 1% accuracy in the matter bispectrum for $z\le3$ on scales $k<1$ $h$/Mpc, 2LPT initial conditions generator with initial redshift of $z\gtrsim100$ is required. We then compare various analytical formulas and empirical fitting functions for modeling the nonlinear matter bispectrum, and discuss the regimes for which each is valid. We find that the next-to-leading order (one-loop) correction from standard perturbation theory matches with N-body results on quasi-linear scales for $z\ge1$. We find that the fitting formula in Gil-Mar\'{\i}n et al. (2012) accurately predicts the matter bispectrum for $z\le1$ on a wide range of scales, but at higher redshifts, the fitting formula given in Scoccimarro & Couchman (2001) gives the best agreement with measurements from N-body simulations.

Baryon acoustic oscillations (BAO) are a powerful probe of the expansion history of the universe, which can tell us about the nature of dark energy. In order to accurately characterize the dark energy equation of state using BAO, we must understand the effects of both nonlinearities and redshift space distortions on the location and shape of the acoustic peak. In a previous paper we introduced a novel approach to 2nd order perturbation theory in configuration space using the Zel'dovich approximation, and presented a simple result for the first nonlinear term of the correlation function. In this paper, we extend this approach to redshift space. We show how perform the computation, and present the analytic result for the first nonlinear term in the correlation function. Finally, we validate our result through comparison to numerical simulations.

Large-scale bulk flows in the Universe distort the initial density field, broadening the baryon-acoustic-oscillation (BAO) feature that was imprinted when baryons were strongly coupled to photons. Additionally, there is a small shift inward in the peak of the conventional overdensity correlation function, a mass-weighted statistic. This shift occurs when high density peaks move toward each other. We explore whether this shift can be removed by applying to the density field a transform (such as a logarithm) that gives fairer statistical weight to fluctuations in underdense regions. Using configuration-space perturbation theory in the Zel'dovich approximation, we find that the log-density correlation function shows a much smaller inward shift in the position of the BAO peak at low redshift than is seen in the overdensity correlation function. We also show that if the initial, Lagrangian density of matter parcels could be estimated at their Eulerian positions, giving a displaced-initial-density field, its peak shift would be even smaller. In fact, a transformed field that accentuates underdensities, such as the reciprocal of the density, pushes the peak the other way, outward. In our model, these shifts in the peak position can be attributed to shift terms, involving the derivative of the linear correlation function, that entirely vanish in this displaced-initial-density field.

Baryon acoustic oscillations are an excellent technique to constrain the properties of dark energy in the Universe. In order to accurately characterize the dark energy equation of state, we must understand the effects of both the nonlinearities and redshift space distortions on the location and shape of the acoustic peak. In this paper, we consider these effects using the Zel'dovich approximation and a novel approach to 2nd order perturbation theory. The second order term of the Zel'dovich power spectrum is built from convolutions of the linear power spectrum with polynomial kernels in Fourier space, suggesting that the corresponding term of the the Zel'dovich correlation function can be written as a sum of quadratic products of a broader class of correlation functions, expressed through simple spherical Bessel transforms of the linear power spectrum. We show how to systematically perform such a computation. We explicitly prove that our result is the Fourier transform of the Zel'dovich power spectrum, and compare our expressions to numerical simulations. Finally, we highlight the advantages of writing the nonlinear expansion in configuration space, as this calculation is easily extended to redshift space, and the higher order terms are mathematically simpler than their Fourier counterparts.