
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 two galaxy shape catalogues from the Dark Energy Survey Year 1
data set, covering 1500 square degrees with a median redshift of $0.59$. The
catalogues cover two main fields: Stripe 82, and an area overlapping the South
Pole Telescope survey region. We describe our data analysis process and in
particular our shape measurement using two independent shear measurement
pipelines, METACALIBRATION and IM3SHAPE. The METACALIBRATION catalogue uses a
Gaussian model with an innovative internal calibration scheme, and was applied
to $riz$bands, yielding 34.8M objects. The IM3SHAPE catalogue uses a
maximumlikelihood bulge/disc model calibrated using simulations, and was
applied to $r$band data, yielding 21.9M objects. Both catalogues pass a suite
of null tests that demonstrate their fitness for use in weak lensing science.
We estimate the 1$\sigma$ uncertainties in multiplicative shear calibration to
be $0.013$ and $0.025$ for the METACALIBRATION and IM3SHAPE catalogues,
respectively.

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 describe the derivation and validation of redshift distribution estimates
and their uncertainties for the galaxies used as weak lensing sources in the
Dark Energy Survey (DES) Year 1 cosmological analyses. The Bayesian Photometric
Redshift (BPZ) code is used to assign galaxies to four redshift bins between
z=0.2 and 1.3, and to produce initial estimates of the lensingweighted
redshift distributions $n^i_{PZ}(z)$ for bin i. Accurate determination of
cosmological parameters depends critically on knowledge of $n^i$ but is
insensitive to bin assignments or redshift errors for individual galaxies. The
cosmological analyses allow for shifts $n^i(z)=n^i_{PZ}(z\Delta z^i)$ to
correct the mean redshift of $n^i(z)$ for biases in $n^i_{\rm PZ}$. The $\Delta
z^i$ are constrained by comparison of independently estimated 30band
photometric redshifts of galaxies in the COSMOS field to BPZ estimates made
from the DES griz fluxes, for a sample matched in fluxes, preseeing size, and
lensing weight to the DES weaklensing sources. In companion papers, the
$\Delta z^i$ are further constrained by the angular clustering of the source
galaxies around red galaxies with secure photometric redshifts at 0.15<z<0.9.
This paper details the BPZ and COSMOS procedures, and demonstrates that the
cosmological inference is insensitive to details of the $n^i(z)$ beyond the
choice of $\Delta z^i$. The clustering and COSMOS validation methods produce
consistent estimates of $\Delta z^i$, with combined uncertainties of
$\sigma_{\Delta z^i}=$0.015, 0.013, 0.011, and 0.022 in the four bins. We
marginalize over these in all analyses to follow, which does not diminish the
constraining power significantly. Repeating the photoz procedure using the
Directional Neighborhood Fitting (DNF) algorithm instead of BPZ, or using the
$n^i(z)$ directly estimated from COSMOS, yields no discernible difference in
cosmological inferences.

We perform a comparison of different approaches to stargalaxy classification
using the broadband photometric data from Year 1 of the Dark Energy Survey.
This is done by performing a wide range of tests with and without external
`truth' information, which can be ported to other similar datasets. We make a
broad evaluation of the performance of the classifiers in two science cases
with DES data that are most affected by this systematic effect: largescale
structure and Milky Way studies. In general, even though the default
morphological classifiers used for DES Y1 cosmology studies are sufficient to
maintain a low level of systematic contamination from stellar
misclassification, contamination can be reduced to the O(1%) level by using
multiepoch and infrared information from external datasets. For Milky Way
studies the stellar sample can be augmented by ~20% for a given flux limit.
Reference catalogs used in this work will be made available upon publication.

We constrain the massrichness scaling relation of redMaPPer galaxy clusters
identified in the Dark Energy Survey Year 1 data using weak gravitational
lensing. We split clusters into $4\times3$ bins of richness $\lambda$ and
redshift $z$ for $\lambda\geq20$ and $0.2 \leq z \leq 0.65$ and measure the
mean masses of these bins using their stacked weak lensing signal. By modeling
the scaling relation as $\langle M_{\rm 200m}\lambda,z\rangle = M_0
(\lambda/40)^F ((1+z)/1.35)^G$, we constrain the normalization of the scaling
relation at the 5.0 per cent level as $M_0 = [3.081 \pm 0.075 ({\rm stat}) \pm
0.133 ({\rm sys})] \cdot 10^{14}\ {\rm M}_\odot$ at $\lambda=40$ and $z=0.35$.
The richness scaling index is constrained to be $F=1.356 \pm 0.051\ ({\rm
stat})\pm 0.008\ ({\rm sys})$ and the redshift scaling index $G=0.30\pm 0.30\
({\rm stat})\pm 0.06\ ({\rm sys})$. These are the tightest measurements of the
normalization and richness scaling index made to date. We use a semianalytic
covariance matrix to characterize the statistical errors in the recovered weak
lensing profiles. Our analysis accounts for the following sources of systematic
error: shear and photometric redshift errors, cluster miscentering, cluster
member dilution of the source sample, systematic uncertainties in the modeling
of the halomass correlation function, halo triaxiality, and projection
effects. We discuss prospects for reducing this systematic error budget, which
dominates the uncertainty on $M_0$. Our result is in excellent agreement with,
but has significantly smaller uncertainties than, previous measurements in the
literature, and augurs well for the power of the DES cluster survey as a tool
for precision cosmology and upcoming galaxy surveys such as LSST, Euclid and
WFIRST.

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.

We use mock galaxy survey simulations designed to resemble the Dark Energy
Survey Year 1 (DES Y1) data to validate and inform cosmological parameter
estimation. When similar analysis tools are applied to both simulations and
real survey data, they provide powerful validation tests of the DES Y1
cosmological analyses presented in companion papers. We use two suites of
galaxy simulations produced using different methods, which therefore provide
independent tests of our cosmological parameter inference. The cosmological
analysis we aim to validate is presented in DES Collaboration et al. (2017) and
uses angular twopoint correlation functions of galaxy number counts and weak
lensing shear, as well as their crosscorrelation, in multiple redshift bins.
While our constraints depend on the specific set of simulated realisations
available, for both suites of simulations we find that the input cosmology is
consistent with the combined constraints from multiple simulated DES Y1
realizations in the $\Omega_m\sigma_8$ plane. For one of the suites, we are
able to show with high confidence that any biases in the inferred
$S_8=\sigma_8(\Omega_m/0.3)^{0.5}$ and $\Omega_m$ are smaller than the DES Y1
$1\sigma$ uncertainties. For the other suite, for which we have fewer
realizations, we are unable to be this conclusive; we infer a roughly 70%
probability that systematic biases in the recovered $\Omega_m$ and $S_8$ are
subdominant to the DES Y1 uncertainty. As cosmological analyses of this kind
become increasingly more precise, validation of parameter inference using
survey simulations will be essential to demonstrate robustness.

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 weak lensing (WL) mass constraints for a sample of massive galaxy
clusters detected by the South Pole Telescope (SPT) via the SunyaevZeldovich
effect (SZE). We use $griz$ imaging data obtained from the Science Verification
(SV) phase of the Dark Energy Survey (DES) to fit the WL shear signal of 33
clusters in the redshift range $0.25 \le z \le 0.8$ with NFW profiles and to
constrain a fourparameter SPT massobservable relation. To account for biases
in WL masses, we introduce a WL mass to true mass scaling relation described by
a mean bias and an intrinsic, lognormal scatter. We allow for correlated
scatter within the WL and SZE massobservable relations and use simulations to
calibrate priors on nuisance parameters related to bias and scatter from WL. We
constrain the normalization of the $\zetaM_{500}$ relation,
$A_\mathrm{SZ}=12.0_{6.7}^{+2.6}$ when using a prior on the mass slope
$B_\mathrm{SZ}$ from the latest SPT cluster cosmology analysis. Without this
prior, we recover $A_\mathrm{SZ}=10.8_{5.2}^{+2.3}$ and
$B_\mathrm{SZ}=1.30_{0.44}^{+0.22}$. Results in both cases imply lower cluster
masses than measured in previous work with and without WL, although the
uncertainties are large. The WL derived value of $B_\mathrm{SZ}$ is $\approx
20\%$ lower than the value preferred by the most recent SPT cluster cosmology
analysis. The method demonstrated in this work is designed to constrain cluster
masses and cosmological parameters simultaneously and will be used for
subsequent studies that employ the full SPT cluster sample together with the
DES data.

In this paper the effect of weak lensing magnification on galaxy number
counts is studied by crosscorrelating the positions of two galaxy samples,
separated by redshift, using data from the Dark Energy Survey Science
Verification dataset. The analysis is carried out for two
photometricallyselected galaxy samples, with mean photometric redshifts in the
$0.2 < z < 0.4$ and $0.7 < z < 1.0$ ranges, in the riz bands. A signal is
detected with a $3.5\sigma$ significance level in each of the bands tested, and
is compatible with the magnification predicted by the $\Lambda$CDM model. After
an extensive analysis, it cannot be attributed to any known systematic effect.
The detection of the magnification signal is robust to estimated uncertainties
in the outlier rate of the pho tometric redshifts, but this will be an
important issue for use of photometric redshifts in magnification mesurements
from larger samples. In addition to the detection of the magnification signal,
a method to select the sample with the maximum signaltonoise is proposed and
validated with data.

Mapping the underlying density field, including nonvisible dark matter,
using weak gravitational lensing measurements is now a standard tool in
cosmology. Due to its importance to the science results of current and upcoming
surveys, the quality of the convergence reconstruction methods should be well
understood. We compare three different mass map reconstruction methods:
KaiserSquires (KS), Wiener filter, and GLIMPSE. KS is a direct inversion
method, taking no account of survey masks or noise. The Wiener filter is well
motivated for Gaussian density fields in a Bayesian framework. The GLIMPSE
method uses sparsity, with the aim of reconstructing nonlinearities in the
density field. We compare these methods with a series of tests on the public
Dark Energy Survey (DES) Science Verification (SV) data and on realistic DES
simulations. The Wiener filter and GLIMPSE methods offer substantial
improvement on the standard smoothed KS with a range of metrics. For both the
Wiener filter and GLIMPSE convergence reconstructions we present a 12%
improvement in Pearson correlation with the underlying truth from simulations.
To compare the mapping methods' abilities to find mass peaks, we measure the
difference between peak counts from simulated {\Lambda}CDM shear catalogues and
catalogues with no mass fluctuations. This is a standard data vector when
inferring cosmology from peak statistics. The maximum signaltonoise value of
these peak statistic data vectors was increased by a factor of 3.5 for the
Wiener filter and by a factor of 9 using GLIMPSE. With simulations we measure
the reconstruction of the harmonic phases, showing that the concentration of
the phase residuals is improved 17% by GLIMPSE and 18% by the Wiener filter. We
show that the correlation between the reconstructions from data and the
foreground redMaPPer clusters is increased 18% by the Wiener filter and 32% by
GLIMPSE. [Abridged]

We construct the largest curvedsky galaxy weak lensing mass map to date from
the DES firstyear (DES Y1) data. The map, about 10 times larger than previous
work, is constructed over a contiguous $\approx1,500 $deg$^2$, covering a
comoving volume of $\approx10 $Gpc$^3$. The effects of masking, sampling, and
noise are tested using simulations. We generate weak lensing maps from two DES
Y1 shear catalogs, Metacalibration and Im3shape, with sources at redshift
$0.2<z<1.3,$ and in each of four bins in this range. In the highest
signaltonoise map, the ratio between the mean signaltonoise in the Emode
and the Bmode map is $\sim$1.5 ($\sim$2) when smoothed with a Gaussian filter
of $\sigma_{G}=30$ (80) arcminutes. The second and third moments of the
convergence $\kappa$ in the maps are in agreement with simulations. We also
find no significant correlation of $\kappa$ with maps of potential systematic
contaminants. Finally, we demonstrate two applications of the mass maps: (1)
crosscorrelation with different foreground tracers of mass and (2) exploration
of the largest peaks and voids in the maps.

We use a suite of simulated images based on Year 1 of the Dark Energy Survey
to explore the impact of galaxy neighbours on shape measurement and shear
cosmology. The hoopoe image simulations include realistic blending, galaxy
positions, and spatial variations in depth and PSF properties. Using the
im3shape maximumlikelihood shape measurement code, we identify four mechanisms
by which neighbours can have a nonnegligible influence on shear estimation.
These effects, if ignored, would contribute a net multiplicative bias of $m
\sim 0.03  0.09$ in the DES Y1 im3shape catalogue, though the precise impact
will be dependent on both the measurement code and the selection cuts applied.
This can be reduced to percentage level or less by removing objects with close
neighbours, at a cost to the effective number density of galaxies
$n_\mathrm{eff}$ of 30%. We use the cosmological inference pipeline of DES Y1
to explore the cosmological implications of neighbour bias and show that
omitting blending from the calibration simulation for DES Y1 would bias the
inferred clustering amplitude $S_8\equiv \sigma_8 (\Omega _\mathrm{m}
/0.3)^{0.5}$ by $2 \sigma$ towards low values. Finally, we use the hoopoe
simulations to test the effect of neighbourinduced spatial correlations in the
multiplicative bias. We find the impact on the recovered $S_8$ of ignoring such
correlations to be subdominant to statistical error at the current level of
precision.

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$.

We use weaklensing 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;
clustermember contamination; miscentering; deviations from the NFW halo
profile; halo triaxiality; and lineofsight 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 massrichness
relation is in excellent agreement with the weaklensing calibration of
redMaPPer clusters in SDSS by Simet et al. (2016) and with the Saro et al.
(2015) calibration based on abundance matching of SPTdetected clusters. Our
results extend the redshift range over which the massrichness 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.

Small temperature anisotropies in the Cosmic Microwave Background can be
sourced by density perturbations via the latetime integrated SachsWolfe
effect. Large voids and superclusters are excellent environments to make a
localized measurement of this tiny imprint. In some cases excess signals have
been reported. We probed these claims with an independent data set, using the
first year data of the Dark Energy Survey in a different footprint, and using a
different superstructure finding strategy. We identified 52 large voids and
102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation
to a priori evaluate the optimal ISW measurement configuration for our
compensated tophat filtering technique, and then performed a stacking
measurement of the CMB temperature field based on the DES data. For optimal
configurations, we detected a cumulative cold imprint of voids with $\Delta
T_{f} \approx 5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta
T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the
expected $\Delta T_{f} \approx 0.6~\mu K$ imprint of such superstructures in
$\Lambda$CDM. If we instead use an a posteriori selected filter size
($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f}
\approx 9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM
expectations and is comparable to previous measurements made using SDSS
superstructure data.

This overview article describes the legacy prospect and discovery potential
of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it
with examples from the DES early data. DES is using a widefield camera (DECam)
on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five
filters (grizY). By its completion the survey is expected to have generated a
catalogue of 300 million galaxies with photometric redshifts and 100 million
stars. In addition, a timedomain survey search over 27 sq deg is expected to
yield a sample of thousands of Type Ia supernovae and other transients. The
main goals of DES are to characterise dark energy and dark matter, and to test
alternative models of gravity; these goals will be pursued by studying large
scale structure, cluster counts, weak gravitational lensing and Type Ia
supernovae. However, DES also provides a rich data set which allows us to study
many other aspects of astrophysics. In this paper we focus on additional
science with DES, emphasizing areas where the survey makes a difference with
respect to other current surveys. The paper illustrates, using early data (from
`Science Verification', and from the first, second and third seasons of
observations), what DES can tell us about the solar system, the Milky Way,
galaxy evolution, quasars, and other topics. In addition, we show that if the
cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then
important astrophysics can be deduced from the primary DES probes. Highlights
from DES early data include the discovery of 34 Trans Neptunian Objects, 17
dwarf satellites of the Milky Way, one published z > 6 quasar (and more
confirmed) and two published superluminous supernovae (and more confirmed).

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).

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 presurvey Dark Energy Survey (DES) Science
Verification optical imaging data and millimeterwave data from the 2500 square
degree South Pole Telescope SunyaevZel'dovich (SPTSZ) survey. The two
lensinggalaxy 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 lensinggalaxy
correlation measurement considered here contains a wealth of information that
can be extracted using current and future surveys.

It is well known that the probability distribution function (PDF) of galaxy
density contrast is approximately lognormal; whether the PDF of mass
fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is
less well established. We derive PDFs of the galaxy and projected matter
density distributions via the Counts in Cells (CiC) method. We use maps of
galaxies and weak lensing convergence produced from the Dark Energy Survey
(DES) Science Verification data over 139 deg^2. We test whether the underlying
density contrast is well described by a lognormal distribution for the
galaxies, the convergence and their joint PDF. We confirm that the galaxy
density contrast distribution is well modeled by a lognormal PDF convolved with
Poisson noise at angular scales from 1040 arcmin (corresponding to physical
scales of 310 Mpc). We note that as kappa_WL is a weighted sum of the mass
fluctuations along the line of sight, its PDF is expected to be only
approximately lognormal. We find that the kappa_WL distribution is well modeled
by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and
20 arcmin, with a bestfit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian
model, corresponding to pvalues 0.35 and 0.07 respectively, at a scale of 10
arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is
also reasonably fitted by a bivariate lognormal. As a consistency check we
compare the variances derived from the lognormal modelling with those directly
measured via CiC. Our methods are validated against maps from the MICE Grand
Challenge Nbody simulation.

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 5year DES, which is expected to cover
5000 square degrees.

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.