We derive cosmological constraints from the probability distribution function (PDF) of evolved large-scale matter density fluctuations. We do this by splitting lines of sight by density based on their count of tracer galaxies, and by measuring both gravitational shear around and counts-in-cells in overdense and underdense lines of sight, in Dark Energy Survey (DES) First Year and Sloan Digital Sky Survey (SDSS) data. Our analysis uses a perturbation theory model (see companion paper Friedrich at al.) and is validated using N-body simulation realizations and log-normal mocks. It allows us to constrain cosmology, bias and stochasticity of galaxies w.r.t. matter density and, in addition, the skewness of the matter density field. From a Bayesian model comparison, we find that the data weakly prefer a connection of galaxies and matter that is stochastic beyond Poisson fluctuations on <=20 arcmin angular smoothing scale. The two stochasticity models we fit yield DES constraints on the matter density $\Omega_m=0.26^{+0.04}_{-0.03}$ and $\Omega_m=0.28^{+0.05}_{-0.04}$ that are consistent with each other. These values also agree with the DES analysis of galaxy and shear two-point functions (3x2pt) that only uses second moments of the PDF. Constraints on $\sigma_8$ are model dependent ($\sigma_8=0.97^{+0.07}_{-0.06}$ and $0.80^{+0.06}_{-0.07}$ for the two stochasticity models), but consistent with each other and with the 3x2pt results if stochasticity is at the low end of the posterior range. As an additional test of gravity, counts and lensing in cells allow to compare the skewness $S_3$ of the matter density PDF to its LCDM prediction. We find no evidence of excess skewness in any model or data set, with better than 25 per cent relative precision in the skewness estimate from DES alone.

We investigate the potential of weak lensing by voids to test for deviations from General Relativity. We calculate the expected lensing signal of a scalar field with derivative couplings, finding that it has the potential to boost the tangential shear both within and outside the void radius. We use voids traced by Luminous Red Galaxies in SDSS to demonstrate the methodology of testing these predictions. We find that the void central density parameter, as inferred from the lensing signal, can shift from its GR value by up to 20% in some galileon gravity models. Since this parameter can be estimated independently using the galaxy tracer profiles of voids, our method provides a consistency check of the gravity theory. Although galileon gravity is now disfavoured as a source of cosmic acceleration by other datasets, the methods we demonstrate here can be used to test for more general fifth force effects with upcoming void lensing data.

We study the ellipticity of galaxy cluster halos as characterized by the distribution of cluster galaxies and as measured with weak lensing. We use monte-carlo simulations of elliptical cluster density profiles to estimate and correct for Poisson noise bias, edge bias and projection effects. We apply our methodology to 10,428 SDSS clusters identified by the redMaPPer algorithm with richness above 20. We find a mean ellipticity $= 0.271 \pm 0.002$ (stat) $\pm 0.031$ (sys) corresponding to an axis ratio $= 0.573 \pm 0.002$ (stat) $\pm 0.039$ (sys). We compare this ellipticity of the satellites to the halo shape, through a stacked lensing measurement using optimal estimators of the lensing quadrupole based on Clampitt and Jain (2016). We find a best-fit axis ratio of $0.56 \pm 0.09$ (stat) $\pm 0.03$ (sys), consistent with the ellipticity of the satellite distribution. Thus cluster galaxies trace the shape of the dark matter halo to within our estimated uncertainties. Finally, we restack the satellite and lensing ellipticity measurements along the major axis of the cluster central galaxy's light distribution. From the lensing measurements we infer a misalignment angle with an RMS of ${30^\circ \pm 10}^\circ$ when stacking on the central galaxy. We discuss applications of halo shape measurements to test the effects of the baryonic gas and AGN feedback, as well as dark matter and gravity. The major improvements in signal-to-noise expected with the ongoing Dark Energy Survey and future surveys from LSST, Euclid and WFIRST will make halo shapes a useful probe of these effects.

We present a measurement of galaxy-galaxy lensing around a magnitude-limited ($i_{AB} < 22.5$) sample of galaxies from the Dark Energy Survey Science Verification (DES-SV) data. We split these lenses into three photometric-redshift bins from 0.2 to 0.8, and determine the product of the galaxy bias $b$ and cross-correlation coefficient between the galaxy and dark matter overdensity fields $r$ in each bin, using scales above 4 Mpc/$h$ comoving, where we find the linear bias model to be valid given our current uncertainties. We compare our galaxy bias results from galaxy-galaxy lensing with those obtained from galaxy clustering (Crocce et al. 2016) and CMB lensing (Giannantonio et al. 2016) for the same sample of galaxies, and find our measurements to be in good agreement with those in Crocce et al. (2016), while, in the lowest redshift bin ($z\sim0.3$), they show some tension with the findings in Giannantonio et al. (2016). We measure $b\cdot r$ to be $0.87\pm 0.11$, $1.12 \pm 0.16$ and $1.24\pm 0.23$, respectively for the three redshift bins of width $\Delta z = 0.2$ in the range $0.2<z <0.8$, defined with the photometric-redshift algorithm BPZ. Using a different code to split the lens sample, TPZ, leads to changes in the measured biases at the 10-20\% level, but it does not alter the main conclusion of this work: when comparing with Crocce et al. (2016) we do not find strong evidence for a cross-correlation parameter significantly below one in this galaxy sample, except possibly at the lowest redshift bin ($z\sim 0.3$), where we find $r = 0.71 \pm 0.11$ when using TPZ, and $0.83 \pm 0.12$ with BPZ.

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

We use weak-lensing shear measurements to determine the mean mass of optically selected galaxy clusters in Dark Energy Survey Science Verification data. In a blinded analysis, we split the sample of more than 8,000 redMaPPer clusters into 15 subsets, spanning ranges in the richness parameter $5 \leq \lambda \leq 180$ and redshift $0.2 \leq z \leq 0.8$, and fit the averaged mass density contrast profiles with a model that accounts for seven distinct sources of systematic uncertainty: shear measurement and photometric redshift errors; cluster-member contamination; miscentering; deviations from the NFW halo profile; halo triaxiality; and line-of-sight projections. We combine the inferred cluster masses to estimate the joint scaling relation between mass, richness and redshift, $\mathcal{M}(\lambda,z) \varpropto M_0 \lambda^{F} (1+z)^{G}$. We find $M_0 \equiv \langle M_{200\mathrm{m}}\,|\,\lambda=30,z=0.5\rangle=\left[ 2.35 \pm 0.22\ \rm{(stat)} \pm 0.12\ \rm{(sys)} \right] \cdot 10^{14}\ M_\odot$, with $F = 1.12\,\pm\,0.20\ \rm{(stat)}\, \pm\, 0.06\ \rm{(sys)}$ and $G = 0.18\,\pm\, 0.75\ \rm{(stat)}\, \pm\, 0.24\ \rm{(sys)}$. The amplitude of the mass-richness relation is in excellent agreement with the weak-lensing calibration of redMaPPer clusters in SDSS by Simet et al. (2016) and with the Saro et al. (2015) calibration based on abundance matching of SPT-detected clusters. Our results extend the redshift range over which the mass-richness relation of redMaPPer clusters has been calibrated with weak lensing from $z\leq 0.3$ to $z\leq0.8$. Calibration uncertainties of shear measurements and photometric redshift estimates dominate our systematic error budget and require substantial improvements for forthcoming studies.

We measure weak lensing mass profiles of voids from a volume-limited sample of SDSS Luminous Red Galaxies (LRGs). We find voids using an algorithm designed to maximize the lensing signal by dividing the survey volume into 2D slices, and then finding holes in this 2D distribution of LRGs. We perform a stacked shear measurement on about 20,000 voids with radii between 15-55 Mpc/h and redshifts between 0.16-0.37. We measure the characteristic radial shear signal of voids with a signal-to-noise of 7. The mass profile corresponds to a fractional underdensity of about -0.4 inside the void radius and a slow approach to the mean density indicating a partially compensated void structure. We compare our measured shape and amplitude with the predictions of Krause et al 2013. Voids in the galaxy distribution have been extensively modeled using simulations and measured in the SDSS. We discuss how the addition of void mass profiles can enable studies of galaxy formation and cosmology.

We present galaxy-galaxy lensing results from 139 square degrees of Dark Energy Survey (DES) Science Verification (SV) data. Our lens sample consists of red galaxies, known as redMaGiC, which are specifically selected to have a low photometric redshift error and outlier rate. The lensing measurement has a total signal-to-noise of 29 over scales $0.09 < R < 15$ Mpc/$h$, including all lenses over a wide redshift range $0.2 < z < 0.8$. Dividing the lenses into three redshift bins for this constant moving number density sample, we find no evidence for evolution in the halo mass with redshift. We obtain consistent results for the lensing measurement with two independent shear pipelines, ngmix and im3shape. We perform a number of null tests on the shear and photometric redshift catalogs and quantify resulting systematic uncertainties. Covariances from jackknife subsamples of the data are validated with a suite of 50 mock surveys. The results and systematics checks in this work provide a critical input for future cosmological and galaxy evolution studies with the DES data and redMaGiC galaxy samples. We fit a Halo Occupation Distribution (HOD) model, and demonstrate that our data constrains the mean halo mass of the lens galaxies, despite strong degeneracies between individual HOD parameters.

Galaxies and their dark matter halos populate a complicated filamentary network around large, nearly empty regions known as cosmic voids. Cosmic voids are usually identified in spectroscopic galaxy surveys, where 3D information about the large-scale structure of the Universe is available. Although an increasing amount of photometric data is being produced, its potential for void studies is limited since photometric redshifts induce line-of-sight position errors of $\sim50$ Mpc/$h$ or more that can render many voids undetectable. In this paper we present a new void finder designed for photometric surveys, validate it using simulations, and apply it to the high-quality photo-$z$ redMaGiC galaxy sample of the Dark Energy Survey Science Verification (DES-SV) data. The algorithm works by projecting galaxies into 2D slices and finding voids in the smoothed 2D galaxy density field of the slice. Fixing the line-of-sight size of the slices to be at least twice the photo-$z$ scatter, the number of voids found in these projected slices of simulated spectroscopic and photometric galaxy catalogs is within 20% for all transverse void sizes, and indistinguishable for the largest voids of radius $\sim 70$ Mpc/$h$ and larger. The positions, radii, and projected galaxy profiles of photometric voids also accurately match the spectroscopic void sample. Applying the algorithm to the DES-SV data in the redshift range $0.2<z<0.8$, we identify 87 voids with comoving radii spanning the range 18-120 Mpc/$h$, and carry out a stacked weak lensing measurement. With a significance of $4.4\sigma$, the lensing measurement confirms the voids are truly underdense in the matter field and hence not a product of Poisson noise, tracer density effects or systematics in the data. It also demonstrates, for the first time in real data, the viability of void lensing studies in photometric surveys.

We measure the correlation of galaxy lensing and cosmic microwave background lensing with a set of galaxies expected to trace the matter density field. The measurements are performed using pre-survey Dark Energy Survey (DES) Science Verification optical imaging data and millimeter-wave data from the 2500 square degree South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. The two lensing-galaxy correlations are jointly fit to extract constraints on cosmological parameters, constraints on the redshift distribution of the lens galaxies, and constraints on the absolute shear calibration of DES galaxy lensing measurements. We show that an attractive feature of these fits is that they are fairly insensitive to the clustering bias of the galaxies used as matter tracers. The measurement presented in this work confirms that DES and SPT data are consistent with each other and with the currently favored $\Lambda$CDM cosmological model. It also demonstrates that joint lensing-galaxy correlation measurement considered here contains a wealth of information that can be extracted using current and future surveys.

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 10<lambda<20, and find that bright (M_i<-21.5) satellites have smaller splashback radii than fainter (M_i>-20) satellites at 99% confidence.

When dark matter halos are accreted by massive host clusters, strong gravitational tidal forces begin stripping mass from the accreted subhalos. This stripping eventually removes all mass beyond a subhalo's tidal radius, but the unbound mass remains in the vicinity of the satellite for at least a dynamical time t_dyn. The N-body subhalo study of Chamberlain et al. verified this picture and pointed out a useful observational consequence: measurements of subhalo correlations beyond the tidal radius are sensitive to the infall time, t_infall, of the subhalo onto its host. We perform this cross-correlation measurement using ~ 160,000 red satellite galaxies in SDSS redMaPPer clusters and find evidence that subhalo correlations do persist well beyond the tidal radius, suggesting that many of the observed satellites fell into their current host less than a dynamical time ago, t_infall < t_dyn. Combined with estimated dynamical times t_dyn ~ 3-5 Gyr and SED fitting results for the time at which satellites stopped forming stars, t_quench ~ 6 Gyr, we infer that for a significant fraction of the satellites, star formation quenched before those satellites entered their current hosts. The result holds for red satellites over a large range of cluster-centric distances 0.1 - 0.6 Mpc/h. We discuss the implications of this result for models of galaxy formation.

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

We present weak lensing shear catalogues for 139 square degrees of data taken during the Science Verification (SV) time for the new Dark Energy Camera (DECam) being used for the Dark Energy Survey (DES). We describe our object selection, point spread function estimation and shear measurement procedures using two independent shear pipelines, IM3SHAPE and NGMIX, which produce catalogues of 2.12 million and 3.44 million galaxies respectively. We detail a set of null tests for the shear measurements and find that they pass the requirements for systematic errors at the level necessary for weak lensing science applications using the SV data. We also discuss some of the planned algorithmic improvements that will be necessary to produce sufficiently accurate shear catalogues for the full 5-year DES, which is expected to cover 5000 square degrees.

Lensing measurements of the shapes of dark matter halos can provide tests of gravity theories and possible dark matter interactions. We measure the quadrupole weak lensing signal from the elliptical halos of 70,000 SDSS Luminous Red Galaxies. We use a new estimator that nulls the spherical halo lensing signal, isolating the shear due to anisotropy in the dark matter distribution. One of the two Cartesian components of our estimator is insensitive to the primary systematic, a spurious alignment of lens and source ellipticities, allowing us to make robust measurements of halo ellipticity. Our best-fit value for the ellipticity of the surface mass density is $0.24 \pm 0.06$, which translates to an axis ratio of 0.78. We rule out the hypothesis of no ellipticity at the $4\sigma$ confidence level, and ellipticity < 0.12 (axis ratio > 0.89) at the $2\sigma$ level. We discuss how our measurements of halo ellipticity are revised to higher values using estimates of the misalignment of mass and light from simulations. Finally, we apply the same techniques to a smaller sample of redMaPPer galaxy clusters and obtain a $3\sigma$ measurement of cluster ellipticity. We discuss how the improved signal to noise properties of our estimator can enable studies of halo shapes for different galaxy populations with upcoming surveys.

We search for the lensing signal of massive filaments between 135,000 pairs of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey. We develop a new estimator that cleanly removes the much larger shear signal of the neighboring LRG halos, relying only on the assumption of spherical symmetry. We consider two models: a thick filament model constructed from ray-tracing simulations for $\Lambda$CDM model, and a thin filament model which models the filament by a string of halos along the line connecting the two LRGs. We show that the filament lensing signal is in nice agreement with the thick simulation filament, while strongly disfavoring the thin model. The magnitude of the lensing shear due to the filament is below $10^{-4}$. Employing the likelihood ratio test, we find a 4.5$\sigma$ significance for the detection of the filament lensing signal, corresponding to a null hypothesis fluctuation probability of $3 \times 10^{-6}$. We also carried out several null tests to verify that the residual shear signal from neighboring LRGs and other shear systematics are minimized.

Using a void catalog from the SDSS survey, we present the first measurements of void clustering and the corresponding void bias. Over the range 30-200 Mpc/h the void auto-correlation is detected at 5-sigma significance for voids of radius 15-20 Mpc/h. We also measure the void-galaxy cross-correlation at higher signal-to-noise and compare the inferred void bias with the autocorrelation results. Void bias is constant with scale for voids of a given size, but its value falls from 5.6 +/- 1.0 to below zero as the void radius increases from 15 to 30 Mpc/h. The comparison of our measurements with carefully matched galaxy mock catalogs, with no free parameters related to the voids, shows that model predictions can be reliably made for void correlations. We study the dependence of void bias on tracer density and void size with a view to future applications. In combination with our previous lensing measurements of void mass profiles, these clustering measurements provide another step towards using voids as cosmological tracers.

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. This analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. We conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9\%/4.8\% with its first-year data that covered over 1000 square degrees, and to 3.9\%/2.3\% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

We measure the weak-lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey. This pathfinder study is meant to 1) validate the DECam imager for the task of measuring weak-lensing shapes, and 2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, PSF modeling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well-behaved, but the modeling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting NFW profiles to the clusters in this study, we determine weak-lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak-lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.

We revisit the excursion set approach to calculate void abundances in chameleon-type modified gravity theories, which was previously studied by Clampitt, Cai and Li (2013). We focus on properly accounting for the void-in-cloud effect, i.e., the growth of those voids sitting in over-dense regions may be restricted by the evolution of their surroundings. This effect may change the distribution function of voids hence affect predictions on the differences between modified gravity and GR. We show that the thin-shell approximation usually used to calculate the fifth force is qualitatively good but quantitatively inaccurate. Therefore, it is necessary to numerically solve the fifth force in both over-dense and under-dense regions. We then generalise the Eulerian void assignment method of Paranjape, Lam and Sheth (2012) to our modified gravity model. We implement this method in our Monte Carlo simulations and compare its results with the original Lagrangian methods. We find that the abundances of small voids are significantly reduced in both modified gravity and GR due to the restriction of environments. However, the change in void abundances for the range of void radii of interest for both models is similar. Therefore, the difference between models remains similar to the results from the Lagrangian method, especially if correlated steps of the random walks are used. As Clampitt, Cai and Li (2013), we find that the void abundance is much more sensitive to modified gravity than halo abundances. Our method can then be a faster alternative to N-body simulations for studying the qualitative behaviour of a broad class of theories. We also discuss the limitations and other practical issues associated with its applications.

We investigate the behavior of the fifth force in voids in chameleon models using the spherical collapse method. Contrary to Newtonian gravity, we find the fifth force is repulsive in voids. The strength of the fifth force depends on the density inside and outside the void region as well as its radius. It can be many times larger than the Newtonian force and their ratio is in principle unbound. This is very different from the case in halos, where the fifth force is no more than 1/3 of gravity. The evolution of voids is governed by the Newtonian gravity, the effective dark energy force and the fifth force. While the first two forces are common in both LCDM and chameleon universes, the fifth force is unique to the latter. Driven by the outward-pointing fifth force, individual voids in chameleon models expand faster and grow larger than in a LCDM universe. The expansion velocity of the void shell can be 20% to 30% larger for voids of a few Mpc/h in radius, while their sizes can be larger by ~10%. These differences are smaller for larger voids of the same density. We compare void statistics using excursion set theory; for voids of the same size, their number density is larger in chameleon models. The fractional difference increases with void size. The chance of having voids of radius ~25 Mpc/h can be 2.5 times larger. This difference is about 10 times larger than that in the halo mass function. We find strong environmental dependence of void properties in chameleon models. The differences in size and expansion velocity with GR are both larger for small voids in high density regions. In general, the difference between chameleon models and LCDM in void properties (size, expansion velocity and distribution function) are larger than the corresponding quantities for halos. This suggests that voids might be better candidates than halos for testing gravity.

We offer predictions of symmetron modified gravity in the neighborhood of realistic dark matter halos. The predictions for the fifth force are obtained by solving the nonlinear symmetron equation of motion in the spherical NFW approximation. In addition, we compare the three major known screening mechanisms: Vainshtein, Chameleon, and Symmetron around such dark matter sources, emphasizing the significant differences between them and highlighting observational tests which exploit these differences. Finally, we demonstrate the host halo environmental screening effect ("blanket screening") on smaller satellite halos by solving for the modified forces around a density profile which is the sum of satellite and approximate host components.