Search

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., Comment: 14 pages, 13 figures, submitted to MNRAS

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., Comment: 5 pages, 1 figure. Submitted to MNRAS Letters

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., Comment: Matches revision accepted in MNRAS

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., Comment: 15 pages, 15 figures, submitted to MNRAS

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., Comment: 14 pages, 11 figures, submitted to MNRAS. V2: fixed an error in model predictions

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., Comment: 8 pages; Accepted by Astrophysical Journal

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., Comment: 5 pages, 3 figures, Accepted by ApjL

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., Comment: Accepted to ApJ. 7 pages, 2 figures

The Primordial Inflation Polarization Explorer (PIPER) is it balloon-borne instrument designed to search for the faint signature of inflation in the polarized component of the cosmic microwave background (C-N-113). Each flight will be configured for a single frequency, but in order to aid in the removal of the polarized foreground signal due to Galactic dust, the filters will be changed between flights. In this way, the CMB polarization at a total of four different frequencies (200, 270, 350, and 600 GHz) will be, measured on large angular scales. PIPER consists of a pair of cryogenic telescopes, one for measuring each of Stokes Q and U in the instrument frame. Each telescope receives both linear orthogonal polarizations in two 32 x 40 element planar arrays that utilize Transition-Edge Sensors (TES). The first element in each telescope is a variable-delay polarization modulator (VPM) that fully modulates the linear Stokes parameter to which the telescope is sensitive. There are several advantages to this architecture. First, by modulating at the front of the optics, instrumental polarization is unmodulated and is therefore cleanly separated from source polarization. Second, by implementing this system with the appropriate symmetry, systematic effects can be further mitigated. In the PIPER design, many of the. systematics are manifest in the unmeasured linear Stokes parameter for each telescope and this can be separated from the desired signal. Finally, the modulation cycle never mixes the Q and U linear Stokes parameters, and thus residuals in the modulation do not twist the observed polarization vector. This is advantageous because measuring the angle of linear polarization is critical for separating the inflationary signal from other polarized components.