Cluster weak lensing is a sensitive probe of cosmology, particularly the amplitude of matter clustering $\sigma_8$ and matter density parameter $\Omega_m$. The main nuisance parameter in a cluster weak lensing cosmological analysis is the scatter between the true halo mass and the relevant cluster observable, denoted $\sigma_{\ln Mc}$. We show that combining the cluster weak lensing observable $\Delta \Sigma$ with the projected cluster-galaxy cross-correlation function $w_{p,cg}$ and galaxy auto-correlation function $w_{p,gg}$ can break the degeneracy between $\sigma_8$ and $\sigma_{\ln Mc}$ to achieve tight, percent-level constraints on $\sigma_8$. Using a grid of cosmological N-body simulations, we compute derivatives of $\Delta \Sigma$, $w_{p,cg}$, and $w_{p,gg}$ with respect to $\sigma_8$, $\Omega_m$, $\sigma_{\ln Mc}$ and halo occupation distribution (HOD) parameters describing the galaxy population. We also compute covariance matrices motivated by the properties of the Dark Energy Suvery (DES) cluster and weak lensing survey and the BOSS CMASS galaxy redshift survey. For our fiducial scenario combining $\Delta \Sigma$, $w_{p,cg}$, and $w_{p,gg}$ measured over $0.3-30.0 \; h^{-1} \; \mathrm{Mpc}$, for clusters at $z=0.35-0.55$ above a mass threshold $M_c\approx 2\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$, we forecast a $1.4\%$ constraint on $\sigma_8$ while marginalizing over $\sigma_{\ln Mc}$ and all HOD parameters. Reducing the mass threshold to $1\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$ and adding a $z=0.15-0.35$ redshift bin sharpens this constraint to $0.8\%$. The small scale $(r_p < 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``mass function'' and large scale $(r_p > 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``halo-mass cross-correlation'' regimes of $\Delta \Sigma$ have comparable constraining power, allowing internal consistency tests from such an analysis.

We present a public data release of halo catalogs from a suite of 125 cosmological $N$-body simulations from the Abacus project. The simulations span 40 $w$CDM cosmologies centered on the Planck 2015 cosmology at two mass resolutions, $4\times 10^{10}\;h^{-1}M_\odot$ and $1\times 10^{10}\;h^{-1}M_\odot$, in $1.1\;h^{-1}\mathrm{Gpc}$ and $720\;h^{-1}\mathrm{Mpc}$ boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning $z=1.5$ to $0.1$ are available for friends-of-friends and Rockstar halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos

In preparation for deep extragalactic imaging with the James Webb Space Telescope, we explore the clustering of massive halos at $z=8$ and $10$ using a large N-body simulation. We find that halos with masses $10^9$ to $10^{11}$ $h^{-1}\;M_\odot$, which are those expected to host galaxies detectable with JWST, are highly clustered with bias factors ranging from 5 and 30 depending strongly on mass, as well as on redshift and scale. This results in correlation lengths of 5--10$h^{-1}\;{\rm Mpc}$, similar to that of today's galaxies. Our results are based on a simulation of 130 billion particles in a box of $250h^{-1}\;{\rm Mpc}$ size using our new high-accuracy ABACUS simulation code, the corrections to cosmological initial conditions of (Garrison et al. 2016, 2016MNRAS.461.4125G), and the Planck 2015 cosmology. We use variations between sub-volumes to estimate the detectability of the clustering. Because of the very strong inter-halo clustering, we find that surveys of order 25$h^{-1}\;{\rm Mpc}$ comoving transverse size may be able to detect the clustering of $z=8$--$10$ galaxies with only 500-1000 survey objects if the galaxies indeed occupy the most massive dark matter halos.

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