
We present cosmological results from a combined analysis of galaxy clustering
and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from
the first year of the Dark Energy Survey (DES Y1). We combine three twopoint
functions: (i) the cosmic shear correlation function of 26 million source
galaxies in four redshift bins, (ii) the galaxy angular autocorrelation
function of 650,000 luminous red galaxies in five redshift bins, and (iii) the
galaxyshear crosscorrelation of luminous red galaxy positions and source
galaxy shears. To demonstrate the robustness of these results, we use
independent pairs of galaxy shape, photometric redshift estimation and
validation, and likelihood analysis pipelines. To prevent confirmation bias,
the bulk of the analysis was carried out while blind to the true results; we
describe an extensive suite of systematics checks performed and passed during
this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM
cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for
$\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino
mass density and including the 457 $\times$ 457 element analytic covariance
matrix. We find consistent cosmological results from these three twopoint
functions, and from their combination obtain $S_8 \equiv \sigma_8
(\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{0.025}$ and $\Omega_m =
0.264^{+0.032}_{0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 =
0.794^{+0.029}_{0.027}$, $\Omega_m = 0.279^{+0.043}_{0.022}$, and
$w=0.80^{+0.20}_{0.22}$ at 68% CL. The precision of these DES Y1 results
rivals that from the Planck cosmic microwave background measurements, allowing
a comparison of structure in the very early and late Universe on equal terms.
Although the DES Y1 bestfit values for $S_8$ and $\Omega_m$ are lower than the
central values from Planck ...

We present galaxygalaxy lensing measurements from 1321 sq. deg. of the Dark
Energy Survey (DES) Year 1 (Y1) data. The lens sample consists of a selection
of 660,000 red galaxies with highprecision photometric redshifts, known as
redMaGiC, split into five tomographic bins in the redshift range $0.15 < z <
0.9$. We use two different source samples, obtained from the Metacalibration
(26 million galaxies) and Im3shape (18 million galaxies) shear estimation
codes, which are split into four photometric redshift bins in the range $0.2 <
z < 1.3$. We perform extensive testing of potential systematic effects that can
bias the galaxygalaxy lensing signal, including those from shear estimation,
photometric redshifts, and observational properties. Covariances are obtained
from jackknife subsamples of the data and validated with a suite of lognormal
simulations. We use the shearratio geometric test to obtain independent
constraints on the mean of the source redshift distributions, providing
validation of those obtained from other photo$z$ studies with the same data.
We find consistency between the galaxy bias estimates obtained from our
galaxygalaxy lensing measurements and from galaxy clustering, therefore
showing the galaxymatter crosscorrelation coefficient $r$ to be consistent
with one, measured over the scales used for the cosmological analysis. The
results in this work present one of the three twopoint correlation functions,
along with galaxy clustering and cosmic shear, used in the DES cosmological
analysis of Y1 data, and hence the methodology and the systematics tests
presented here provide a critical input for that study as well as for future
cosmological analyses in DES and other photometric galaxy surveys.

We measure the clustering of DES Year 1 galaxies that are intended to be
combined with weak lensing samples in order to produce precise cosmological
constraints from the joint analysis of largescale structure and lensing
correlations. Twopoint correlation functions are measured for a sample of $6.6
\times 10^{5}$ luminous red galaxies selected using the \textsc{redMaGiC}
algorithm over an area of $1321$ square degrees, in the redshift range $0.15 <
z < 0.9$, split into five tomographic redshift bins. The sample has a mean
redshift uncertainty of $\sigma_{z}/(1+z) = 0.017$. We quantify and correct
spurious correlations induced by spatially variable survey properties, testing
their impact on the clustering measurements and covariance. We demonstrate the
sample's robustness by testing for stellar contamination, for potential biases
that could arise from the systematic correction, and for the consistency
between the twopoint auto and crosscorrelation functions. We show that the
corrections we apply have a significant impact on the resultant measurement of
cosmological parameters, but that the results are robust against arbitrary
choices in the correction method. We find the linear galaxy bias in each
redshift bin in a fiducial cosmology to be $b(z$=$0.24)=1.40 \pm 0.08$,
$b(z$=$0.38)=1.61 \pm 0.05$, $b(z$=$0.53)=1.60 \pm 0.04$ for galaxies with
luminosities $L/L_*>$$0.5$, $b(z$=$0.68)=1.93 \pm 0.05$ for $L/L_*>$$1$ and
$b(z$=$0.83)=1.99 \pm 0.07$ for $L/L_*$$>1.5$, broadly consistent with
expectations for the redshift and luminosity dependence of the bias of red
galaxies. We show these measurements to be consistent with the linear bias
obtained from tangential shear measurements.

We derive cosmological constraints from the probability distribution function
(PDF) of evolved largescale 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 countsincells 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
Nbody simulation realizations and lognormal 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 twopoint 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 use 26 million galaxies from the Dark Energy Survey (DES) Year 1 shape
catalogs over 1321 deg$^2$ of the sky to produce the most significant
measurement of cosmic shear in a galaxy survey to date. We constrain
cosmological parameters in both the flat $\Lambda$CDM and $w$CDM models, while
also varying the neutrino mass density. These results are shown to be robust
using two independent shape catalogs, two independent \photoz\ calibration
methods, and two independent analysis pipelines in a blind analysis. We find a
3.5\% fractional uncertainty on $\sigma_8(\Omega_m/0.3)^{0.5} =
0.782^{+0.027}_{0.027}$ at 68\% CL, which is a factor of 2.5 improvement over
the fractional constraining power of our DES Science Verification results. In
$w$CDM, we find a 4.8\% fractional uncertainty on $\sigma_8(\Omega_m/0.3)^{0.5}
= 0.777^{+0.036}_{0.038}$ and a dark energy equationofstate
$w=0.95^{+0.33}_{0.39}$. We find results that are consistent with previous
cosmic shear constraints in $\sigma_8$  $\Omega_m$, and see no evidence for
disagreement of our weak lensing data with data from the CMB. Finally, we find
no evidence preferring a $w$CDM model allowing $w\ne 1$. We expect further
significant improvements with subsequent years of DES data, which will more
than triple the sky coverage of our shape catalogs and double the effective
integrated exposure time per galaxy.

We explore the impact of an update to the typical approximation for the shape
noise term in the analytic covariance matrix for cosmic shear experiments that
assumes the absence of survey boundary and mask effects. We present an exact
expression for the number of galaxy pairs in this term based on the the survey
mask, which leads to more than a factor of three increase in the shape noise on
the largest measured scales for the KiloDegree Survey (KIDS450) realspace
cosmic shear data. We compare the result of this analytic expression to several
alternative methods for measuring the shape noise from the data and find
excellent agreement. This update to the covariance resolves any internal model
tension evidenced by the previously large cosmological bestfit $\chi^2$ for
the KiDS450 cosmic shear data. The bestfit $\chi^2$ is reduced from 161 to
121 for 118 degrees of freedom. We also apply a correction to how the
multiplicative shear calibration uncertainty is included in the covariance.
This change, along with a previously known update to the reported effective
angular values of the data vector, jointly shift the inferred amplitude of the
correlation function to higher values. We find that this improves agreement of
the KiDS450 cosmic shear results with Dark Energy Survey Year 1 and Planck
results.

We present density split statistics, a framework that studies lensing and
countsincells as a function of foreground galaxy density, thereby providing a
largescale measurement of both 2point and 3point statistics. Our method
extends our earlier work on trough lensing and is summarized as follows: given
a foreground (low redshift) population of galaxies, we divide the sky into
subareas of equal size but distinct galaxy density. We then measure lensing
around uniformly spaced points separately in each of these subareas, as well as
countsincells statistics (CiC). The lensing signals trace the matter density
contrast around regions of fixed galaxy density. Through the CiC measurements
this can be related to the density profile around regions of fixed matter
density. Together, these measurements constitute a powerful probe of cosmology,
the skewness of the density field and the connection of galaxies and matter.
In this paper we show how to model both the density split lensing signal and
CiC from basic ingredients: a nonlinear power spectrum, clustering hierarchy
coefficients from perturbation theory and a parametric model for galaxy bias
and shotnoise. Using Nbody simulations, we demonstrate that this model is
sufficiently accurate for a cosmological analysis on year 1 data from the Dark
Energy Survey.

Optical imaging surveys measure both the galaxy density and the gravitational
lensinginduced shear fields across the sky. Recently, the Dark Energy Survey
(DES) collaboration used a joint fit to twopoint correlations between these
observables to place tight constraints on cosmology (DES Collaboration et al.
2017). In this work, we develop the methodology to extend the DES Collaboration
et al. (2017) analysis to include crosscorrelations of the optical survey
observables with gravitational lensing of the cosmic microwave background (CMB)
as measured by the South Pole Telescope (SPT) and Planck. Using simulated
analyses, we show how the resulting set of five twopoint functions increases
the robustness of the cosmological constraints to systematic errors in galaxy
lensing shear calibration. Additionally, we show that contamination of the
SPT+Planck CMB lensing map by the thermal SunyaevZel'dovich effect is a
potentially large source of systematic error for twopoint function analyses,
but show that it can be reduced to acceptable levels in our analysis by masking
clusters of galaxies and imposing angular scale cuts on the twopoint
functions. The methodology developed here will be applied to the analysis of
data from the DES, the SPT, and Planck in a companion work.

We present the methodology for and detail the implementation of the Dark
Energy Survey (DES) 3x2pt DES Year 1 (Y1) analysis, which combines
configurationspace twopoint statistics from three different cosmological
probes: cosmic shear, galaxygalaxy lensing, and galaxy clustering, using data
from the first year of DES observations. We have developed two independent
modeling pipelines and describe the code validation process. We derive
expressions for analytical realspace multiprobe covariances, and describe
their validation with numerical simulations. We stresstest the inference
pipelines in simulated likelihood analyses that vary 67 cosmology parameters
plus 20 nuisance parameters and precisely resemble the analysis to be presented
in the DES 3x2pt analysis paper, using a variety of simulated input data
vectors with varying assumptions.
We find that any disagreement between pipelines leads to changes in assigned
likelihood $\Delta \chi^2 \le 0.045$ with respect to the statistical error of
the DES Y1 data vector. We also find that angular binning and survey mask do
not impact our analytic covariance at a significant level. We determine lower
bounds on scales used for analysis of galaxy clustering (8 Mpc$~h^{1}$) and
galaxygalaxy lensing (12 Mpc$~h^{1}$) such that the impact of modeling
uncertainties in the nonlinear regime is well below statistical errors, and
show that our analysis choices are robust against a variety of systematics.
These tests demonstrate that we have a robust analysis pipeline that yields
unbiased cosmological parameter inferences for the flagship 3x2pt DES Y1
analysis. We emphasize that the level of independent code development and
subsequent code comparison as demonstrated in this paper is necessary to
produce credible constraints from increasingly complex multiprobe analyses of
current data.

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 2point 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 stateoftheart
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$.

Data resampling methods such as the deleteone jackknife are a common tool
for estimating the covariance of large scale structure probes. In this paper we
investigate the concepts of internal covariance estimation in the context of
cosmic shear twopoint statistics. We demonstrate how to use lognormal
simulations of the convergence field and the corresponding shear field to carry
out realistic tests of internal covariance estimators and find that most
estimators such as jackknife or subsample covariance can reach a satisfactory
compromise between bias and variance of the estimated covariance.
In a forecast for the complete, 5year DES survey we show that internally
estimated covariance matrices can provide a large fraction of the true
uncertainties on cosmological parameters in a 2D cosmic shear analysis. The
volume inside contours of constant likelihood in the $\Omega_m$$\sigma_8$
plane as measured with internally estimated covariance matrices is on average
$\gtrsim 85\%$ of the volume derived from the true covariance matrix. The
uncertainty on the parameter combination $\Sigma_8 \sim \sigma_8
\Omega_m^{0.5}$ derived from internally estimated covariances is $\sim 90\%$ of
the true uncertainty.

We present measurements of weak gravitational lensing cosmic shear twopoint
statistics using Dark Energy Survey Science Verification data. We demonstrate
that our results are robust to the choice of shear measurement pipeline, either
ngmix or im3shape, and robust to the choice of twopoint statistic, including
both real and Fourierspace statistics. Our results pass a suite of null tests
including tests for Bmode contamination and direct tests for any dependence of
the twopoint functions on a set of 16 observing conditions and galaxy
properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We
furthermore use a large suite of simulations to compute the covariance matrix
of the cosmic shear measurements and assign statistical significance to our
null tests. We find that our covariance matrix is consistent with the halo
model prediction, indicating that it has the appropriate level of halo sample
variance. We compare the same jackknife procedure applied to the data and the
simulations in order to search for additional sources of noise not captured by
the simulations. We find no statistically significant extra sources of noise in
the data. The overall detection significance with tomography for our highest
source density catalog is 9.7sigma. Cosmological constraints from the
measurements in this work are presented in a companion paper (DES et al. 2015).

Shear peak statistics has gained a lot of attention recently as a practical
alternative to the two point statistics for constraining cosmological
parameters. We perform a shear peak statistics analysis of the Dark Energy
Survey (DES) Science Verification (SV) data, using weak gravitational lensing
measurements from a 139 deg$^2$ field. We measure the abundance of peaks
identified in aperture mass maps, as a function of their signaltonoise ratio,
in the signaltonoise range $0<\mathcal S / \mathcal N<4$. To predict the peak
counts as a function of cosmological parameters we use a suite of $N$body
simulations spanning 158 models with varying $\Omega_{\rm m}$ and $\sigma_8$,
fixing $w = 1$, $\Omega_{\rm b} = 0.04$, $h = 0.7$ and $n_s=1$, to which we
have applied the DES SV mask and redshift distribution. In our fiducial
analysis we measure $\sigma_{8}(\Omega_{\rm m}/0.3)^{0.6}=0.77 \pm 0.07$, after
marginalising over the shear multiplicative bias and the error on the mean
redshift of the galaxy sample. We introduce models of intrinsic alignments,
blending, and source contamination by cluster members. These models indicate
that peaks with $\mathcal S / \mathcal N>4$ would require significant
corrections, which is why we do not include them in our analysis. We compare
our results to the cosmological constraints from the two point analysis on the
SV field and find them to be in good agreement in both the central value and
its uncertainty. We discuss prospects for future peak statistics analysis with
upcoming DES data.

We measure the weak lensing shear around galaxy troughs, i.e. the radial
alignment of background galaxies relative to underdensities in projections of
the foreground galaxy field over a wide range of redshift in Science
Verification data from the Dark Energy Survey. Our detection of the shear
signal is highly significant (10 to 15$\sigma$ for the smallest angular scales)
for troughs with the redshift range z in [0.2,0.5] of the projected galaxy
field and angular diameters of 10 arcmin...1{\deg}. These measurements probe
the connection between the galaxy, matter density, and convergence fields. By
assuming galaxies are biased tracers of the matter density with Poissonian
noise, we find agreement of our measurements with predictions in a fiducial
Lambda cold dark matter model. The prediction for the lensing signal on large
trough scales is virtually independent of the details of the underlying model
for the connection of galaxies and matter. Our comparison of the shear around
troughs with that around cylinders with large galaxy counts is consistent with
a symmetry between galaxy and matter over and underdensities. In addition, we
measure the twopoint angular correlation of troughs with galaxies which, in
contrast to the lensing signal, is sensitive to galaxy bias on all scales. The
lensing signal of troughs and their clustering with galaxies is therefore a
promising probe of the statistical properties of matter underdensities and
their connection to the galaxy field.

Intrinsic variations of the projected density profiles of clusters of
galaxies at fixed mass are a source of uncertainty for cluster weak lensing. We
present a semianalytical model to account for this effect, based on a
combination of variations in halo concentration, ellipticity and orientation,
and the presence of correlated haloes. We calibrate the parameters of our model
at the 10 per cent level to match the empirical cosmic variance of cluster
profiles at M_200m=10^14...10^15 h^1 M_sol, z=0.25...0.5 in a cosmological
simulation. We show that weak lensing measurements of clusters significantly
underestimate mass uncertainties if intrinsic profile variations are ignored,
and that our model can be used to provide correct mass likelihoods. Effects on
the achievable accuracy of weak lensing cluster mass measurements are
particularly strong for the most massive clusters and deep observations (with
~20 per cent uncertainty from cosmic variance alone at M_200m=10^15 h^1 M_sol
and z=0.25), but significant also under typical groundbased conditions. We
show that neglecting intrinsic profile variations leads to biases in the
massobservable relation constrained with weak lensing, both for intrinsic
scatter and overall scale (the latter at the 15 per cent level). These biases
are in excess of the statistical errors of upcoming surveys and can be avoided
if the cosmic variance of cluster profiles is accounted for.