The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $\xi_{\text{gg}}(r)$ and galaxy-matter cross-correlation $\xi_{\text{gm}}(r)$ as a function of $\sigma_8$, $\Omega_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $\sigma_8\Omega_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $\gamma_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $\sigma_8\Omega_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey.

In cosmological $N$-body simulations, the representation of dark matter as discrete "macroparticles" suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near $k_{\rm Nyquist}$. Marcos et al. demonstrate that this is due to sparse sampling of modes near $k_{\rm Nyquist}$ and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate $N$-body code Abacus to implement 2LPT in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1% accuracy in the power spectrum down to $k_{\rm Nyquist}$, versus $k_{\rm Nyquist}/4$ without rescaling or $k_{\rm Nyquist}/13$ without 2LPT, relative to an oversampled reference simulation. We anticipate that our 2LPT will be useful for large simulations where FFTs are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogs. Code to generate initial conditions is available at https://github.com/lgarrison/zeldovich-PLT

The baryon acoustic oscillation (BAO) feature in the clustering of matter in the universe serves as a robust standard ruler and hence can be used to map the expansion history of the universe. We use high force resolution simulations to analyze the effects of galaxy bias on the measurements of the BAO signal. We apply a variety of Halo Occupation Distributions (HODs) and produce biased mass tracers to mimic different galaxy populations. We investigate whether galaxy bias changes the non-linear shifts on the acoustic scale relative to the underlying dark matter distribution presented by Seo et al (2009). For the less biased HOD models (b < 3), we do not detect any shift in the acoustic scale relative to the no-bias case, typically 0.10% \pm 0.10%. However, the most biased HOD models (b > 3) show a shift at moderate significance (0.79% \pm 0.31% for the most extreme case). We test the one-step reconstruction technique introduced by Eisenstein et al. (2007) in the case of realistic galaxy bias and shot noise. The reconstruction scheme increases the correlation between the initial and final (z = 1) density fields achieving an equivalent level of correlation at nearly twice the wavenumber after reconstruction. Reconstruction reduces the shifts and errors on the shifts. We find that after reconstruction the shifts from the galaxy cases and the dark matter case are consistent with each other and with no shift. The 1-sigma systematic errors on the distance measurements inferred from our BAO measurements with various HODs after reconstruction are about 0.07% - 0.15%.

We measure shifts of the acoustic scale due to nonlinear growth and redshift distortions to a high precision using a very large volume of high-force-resolution simulations. We compare results from various sets of simulations that differ in their force, volume, and mass resolution. We find a consistency within 1.5-sigma for shift values from different simulations and derive shift alpha(z) -1 = (0.300\pm 0.015)% [D(z)/D(0)]^{2} using our fiducial set. We find a strong correlation with a non-unity slope between shifts in real space and in redshift space and a weak correlation between the initial redshift and low redshift. Density-field reconstruction not only removes the mean shifts and reduces errors on the mean, but also tightens the correlations: after reconstruction, we recover a slope of near unity for the correlation between the real and redshift space and restore a strong correlation between the low and the initial redshifts. We derive propagators and mode-coupling terms from our N-body simulations and compared with Zeldovich approximation and the shifts measured from the chi^2 fitting, respectively. We interpret the propagator and the mode-coupling term of a nonlinear density field in the context of an average and a dispersion of its complex Fourier coefficients relative to those of the linear density field; from these two terms, we derive a signal-to-noise ratio of the acoustic peak measurement. We attempt to improve our reconstruction method by implementing 2LPT and iterative operations: we obtain little improvement. The Fisher matrix estimates of uncertainty in the acoustic scale is tested using 5000 (Gpc/h)^3 of cosmological PM simulations from Takahashi et al. (2009). (abridged)

We introduce a new statistic omega_l for measuring and analyzing large-scale structure and particularly the baryon acoustic oscillations. omega_l is a band-filtered, configuration space statistic that is easily implemented and has advantages over the traditional power spectrum and correlation function estimators. Unlike these estimators, omega_l can localize most of the acoustic information into a single dip at the acoustic scale while also avoiding sensitivity to the poorly constrained large scale power (i.e., the integral constraint) through the use of a localized and compensated filter. It is also sensitive to anisotropic clustering through pair counting and does not require any binning. We measure the shift in the acoustic peak due to nonlinear effects using the monopole omega_0 derived from subsampled dark matter catalogues as well as from mock galaxy catalogues created via halo occupation distribution (HOD) modeling. All of these are drawn from 44 realizations of 1024^3 particle dark matter simulations in a 1h^{-1}Gpc box at z=1. We compare these shifts with those obtained from the power spectrum and conclude that the results agree. This indicates that any distance measurements obtained from omega_0 and P(k) will be consistent with each other. We also show that it is possible to extract the same amount of acoustic information using either omega_0 or P(k) from equal volume surveys.

We have carried out simulations of supersonic light jets in order to model the features observed in optical and radio images of the western hot-spot in the radio galaxy Pictor A. We have considered jets with density ratios eta=1e-2 - 1e-4, and Mach numbers ranging between 5 and 50. From each simulation, we have generated ray-traced maps of radio surface brightness at a variety of jet inclinations, in order to study the appearance of time-dependent luminous structures in the vicinity of the western hot-spot. We compare these rendered images with observed features of Pictor A. A remarkable feature of Pictor A observations is a bar-shaped "filament" inclined almost at right angles to the inferred jet direction and extending 24" (10.8 /h kpc) along its longest axis. The constraints of reproducing the appearance of this structure in simulations indicate that the jet of Pictor A lies nearly in the plane of the sky. The results of the simulation are also consistent with other features found in the radio image of Pictor A. This filament arises from the surging behaviour of the jet near the hot-spot; the surging is provoked by alternate compression and decompression of the jet by the turbulent backflow in the cocoon. We also examine the arguments for the jet in Pictor A being at a more acute angle to the line of sight and find that our preferred orientation is just consistent with the limits on the brightness ratio of the X-ray jet and counter-jet. We determine from our simulations, the structure function of hot-spot brightness and also the cumulative distribution of the ratio of intrinsic hot-spot brightnesses. The latter may be used to quantify the use of hot-spot ratios for the estimation of relativistic effects.