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

Splashback refers to the process of matter that is accreting onto a dark
matter halo reaching its first orbital apocenter and turning around in its
orbit. The clustercentric radius at which this process occurs, r_sp, defines a
halo boundary that is connected to the dynamics of the cluster. A rapid decline
in the halo profile is expected near r_sp. We measure the galaxy number density
and weak lensing mass profiles around redMaPPer galaxy clusters in the first
year Dark Energy Survey (DES) data. For a cluster sample with mean M_200m mass
~2.5 x 10^14 M_sun, we find strong evidence of a splashbacklike steepening of
the galaxy density profile and measure r_sp=1.13 +/ 0.07 Mpc/h, consistent
with earlier SDSS measurements of More et al. (2016) and Baxter et al. (2017).
Moreover, our weak lensing measurement demonstrates for the first time the
existence of a splashbacklike steepening of the matter profile of galaxy
clusters. We measure r_sp=1.34 +/ 0.21 Mpc/h from the weak lensing data, in
good agreement with our galaxy density measurements. For different cluster and
galaxy samples, we find that consistent with LCDM simulations, r_sp scales with
R_200m and does not evolve with redshift over the redshift range of 0.30.6.
We also find that potential systematic effects associated with the redMaPPer
algorithm may impact the location of r_sp. We discuss progress needed to
understand the systematic uncertainties and fully exploit forthcoming data from
DES and future surveys, emphasizing the importance of more realistic mock
catalogs and independent cluster samples.

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 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 describe the creation, content, and validation of the Dark Energy Survey
(DES) internal yearone cosmology data set, Y1A1 GOLD, in support of upcoming
cosmological analyses. The Y1A1 GOLD data set is assembled from multiple epochs
of DES imaging and consists of calibrated photometric zeropoints, object
catalogs, and ancillary data products  e.g., maps of survey depth and
observing conditions, stargalaxy classification, and photometric redshift
estimates  that are necessary for accurate cosmological analyses. The Y1A1
GOLD widearea object catalog consists of ~137 million objects detected in
coadded images covering ~1800 deg$^2$ in the DES grizY filters. The 10{\sigma}
limiting magnitude for galaxies is g = 23.4, r = 23.2, i = 22.5, z = 21.8, and
Y = 20.1. Photometric calibration of Y1A1 GOLD was performed by combining
nightly zeropoint solutions with stellarlocus regression, and the absolute
calibration accuracy is better than 2% over the survey area. DES Y1A1 GOLD is
the largest photometric data set at the achieved depth to date, enabling
precise measurements of cosmic acceleration at z $\lesssim$ 1.

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

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.

Hubble Space Telescope observations from December 11 2015 detected the
expected fifth counter image of SN Refsdal at $z = 1.49$. In this letter, we
compare the time delay predictions from numerous models with the measured value
derived by \citet{Kelly2016a} from very early data in the light curve of the SN
Refsdal, and find a best value for $H_0 =
64^{+9}_{11}~\mathrm{km~s^{1}~Mpc^{1}}$ (68\% CL), in excellent agreement
with predictions from CMB and recent weak lensing data + BAO + BBN (from the
DES Collaboration). This is the first constraint on $H_0$ derived from time
delays between multiple lensed SN images, and the first with a galaxy cluster
lens, so subject to systematic effects different from other timedelay $H_0$
estimates. Additional time delay measurements from new multiplyimaged SNe will
allow derivation of competitive constraints on $H_0$.

The Dark Energy Survey (DES) is a fiveyear optical imaging campaign with the
goal of understanding the origin of cosmic acceleration. DES performs a 5000
square degree survey of the southern sky in five optical bands (g,r,i,z,Y) to a
depth of ~24th magnitude. Contemporaneously, DES performs a deep, timedomain
survey in four optical bands (g,r,i,z) over 27 square degrees. DES exposures
are processed nightly with an evolving data reduction pipeline and evaluated
for image quality to determine if they need to be retaken. Difference imaging
and transient source detection are also performed in the time domain component
nightly. On a biannual basis, DES exposures are reprocessed with a refined
pipeline and coadded to maximize imaging depth. Here we describe the DES image
processing pipeline in support of DES science, as a reference for users of
archival DES data, and as a guide for future astronomical surveys.

We describe the first public data release of the Dark Energy Survey, DES DR1,
consisting of reduced single epoch images, coadded images, coadded source
catalogs, and associated products and services assembled over the first three
years of DES science operations. DES DR1 is based on optical/nearinfrared
imaging from 345 distinct nights (August 2013 to February 2016) by the Dark
Energy Camera mounted on the 4m Blanco telescope at Cerro Tololo
InterAmerican Observatory in Chile. We release data from the DES widearea
survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad
photometric bands, grizY. DES DR1 has a median delivered pointspread function
of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a
photometric precision of < 1% in all bands, and an astrometric precision of 151
mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10
is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1
includes nearly 400M distinct astronomical objects detected in ~10,000 coadd
tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures.
Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects,
respectively, following a basic object quality selection. These data are
accessible through a range of interfaces, including query web clients, image
cutout servers, jupyter notebooks, and an interactive coadd image visualization
tool. DES DR1 constitutes the largest photometric data set to date at the
achieved depth and photometric precision.

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 characterize the variation in photometric response of the Dark Energy
Camera (DECam) across its 520~Mpix science array during 4 years of operation.
These variations are measured using high signaltonoise aperture photometry of
$>10^7$ stellar images in thousands of exposures of a few selected fields, with
the telescope dithered to move the sources around the array. A calibration
procedure based on these results brings the RMS variation in aperture
magnitudes of bright stars on cloudless nights down to 23 mmag, with <1 mmag
of correlated photometric errors for stars separated by $\ge20$". On cloudless
nights, any departures of the exposure zeropoints from a secant airmass law
exceeding >1 mmag are plausibly attributable to spatial/temporal variations in
aperture corrections. These variations can be inferred and corrected by
measuring the fraction of stellar light in an annulus between 6" and 8"
diameter. Key elements of this calibration include: correction of amplifier
nonlinearities; distinguishing pixelarea variations and stray light from
quantumefficiency variations in the flat fields; fielddependent color
corrections; and the use of an aperturecorrection proxy. The DECam response
pattern across the 2degree field drifts over months by up to $\pm7$ mmag, in a
nearlywavelengthindependent loworder pattern. We find no fundamental
barriers to pushing global photometric calibrations toward mmag accuracy.

We present the calibration of the Dark Energy Survey Year 1 (DES Y1) weak
lensing source galaxy redshift distributions from clustering measurements. By
crosscorrelating the positions of source galaxies with luminous red galaxies
selected by the redMaGiC algorithm we measure the redshift distributions of the
source galaxies as placed into different tomographic bins. These measurements
constrain any such shifts to an accuracy of $\sim0.02$ and can be computed even
when the clustering measurements do not span the full redshift range. The
highestredshift source bin is not constrained by the clustering measurements
because of the minimal redshift overlap with the redMaGiC galaxies. We compare
our constraints with those obtained from $\texttt{COSMOS}$ 30band photometry
and find that our two very different methods produce consistent constraints.

We use numerical simulations to characterize the performance of a
clusteringbased method to calibrate photometric redshift biases. In
particular, we crosscorrelate the weak lensing (WL) source galaxies from the
Dark Energy Survey Year 1 (DES Y1) sample with redMaGiC galaxies (luminous red
galaxies with secure photometric redshifts) to estimate the redshift
distribution of the former sample. The recovered redshift distributions are
used to calibrate the photometric redshift bias of standard photo$z$ methods
applied to the same source galaxy sample. We apply the method to three
photo$z$ codes run in our simulated data: Bayesian Photometric Redshift (BPZ),
Directional Neighborhood Fitting (DNF), and Random Forestbased photo$z$ (RF).
We characterize the systematic uncertainties of our calibration procedure, and
find that these systematic uncertainties dominate our error budget. The
dominant systematics are due to our assumption of unevolving bias and
clustering across each redshift bin, and to differences between the shapes of
the redshift distributions derived by clustering vs photo$z$'s. The systematic
uncertainty in the mean redshift bias of the source galaxy sample is $\Delta z
\lesssim 0.02$, though the precise value depends on the redshift bin under
consideration. We discuss possible ways to mitigate the impact of our dominant
systematics in future analyses.

We describe the model for the mapping from sky brightness to the digital
output of the Dark Energy Camera, and the algorithms adopted by the Dark Energy
Survey (DES) for inverting this model to obtain photometric measures of
celestial objects from the raw camera output. The calibration aims for fluxes
that are uniform across the camera field of view and across the full angular
and temporal span of the DES observations, approaching the accuracy limits set
by shot noise for the full dynamic range of DES observations. The DES pipeline
incorporates several substantive advances over standard detrending techniques,
including: principalcomponentsbased sky and fringe subtraction; correction of
the "brighterfatter" nonlinearity; use of internal consistency in onsky
observations to disentangle the influences of quantum efficiency, pixelsize
variations, and scattered light in the dome flats; and pixelbypixel
characterization of instrument spectral response, through combination of
internalconsistency constraints with auxiliary calibration data. This article
provides conceptual derivations of the detrending/calibration steps, and the
procedures for obtaining the necessary calibration data. Other publications
will describe the implementation of these concepts for the DES operational
pipeline, the detailed methods, and the validation that the techniques can
bring DECam photometry and astrometry within ~2 mmag and ~3 mas, respectively,
of fundamental atmospheric and statistical limits. The DES techniques should be
broadly applicable to widefield imagers.

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.

Many scientific goals for the Dark Energy Survey (DES) require calibration of
optical/NIR broadband $b = grizY$ photometry that is stable in time and uniform
over the celestial sky to one percent or better. It is also necessary to limit
to similar accuracy systematic uncertainty in the calibrated broadband
magnitudes due to uncertainty in the spectrum of the source. Here we present a
"Forward Global Calibration Method (FGCM)" for photometric calibration of the
DES, and we present results of its application to the first three years of the
survey (Y3A1). The FGCM combines data taken with auxiliary instrumentation at
the observatory with data from the broadband survey imaging itself and models
of the instrument and atmosphere to estimate the spatial and timedependence
of the passbands of individual DES survey exposures. "Standard" passbands are
chosen that are typical of the passbands encountered during the survey. The
passband of any individual observation is combined with an estimate of the
source spectral shape to yield a magnitude $m_b^{\mathrm{std}}$ in the standard
system. This "chromatic correction" to the standard system is necessary to
achieve subpercent calibrations. The FGCM achieves reproducible and stable
photometric calibration of standard magnitudes $m_b^{\mathrm{std}}$ of stellar
sources over the multiyear Y3A1 data sample with residual random calibration
errors of $\sigma=56\,\mathrm{mmag}$ per exposure. The accuracy of the
calibration is uniform across the $5000\,\mathrm{deg}^2$ DES footprint to
within $\sigma=7\,\mathrm{mmag}$. The systematic uncertainties of magnitudes in
the standard system due to the spectra of sources are less than
$5\,\mathrm{mmag}$ for main sequence stars with $0.5<gi<3.0$.

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.

We report the discovery of two new candidate stellar systems in the
constellation of Cetus using the data from the first two years of the Dark
Energy Survey (DES). The objects, DES J01111341 and DES J0225+0304, are
located at a heliocentric distance of ~ 25 kpc and appear to have old and
metalpoor populations. Their distances to the Sagittarius orbital plane, ~
1.73 kpc (DES J01111341) and ~ 0.50 kpc (DES J0225+0304), indicate that they
are possibly associated with the Sagittarius dwarf stream. The halflight
radius (r_h ~ 4.55 pc) and luminosity (M_V ~ +0.3) of DES J01111341 are
consistent with it being an ultrafaint stellar cluster, while the halflight
radius (r_h ~ 18.55 pc) and luminosity (M_V ~ 1.1) of DES J0225+0304 place it
in an ambiguous region of sizeluminosity space between stellar clusters and
dwarf galaxies. Determinations of the characteristic parameters of the
Sagittarius stream, metallicity spread (2.18 < [Fe/H] < 0.95) and distance
gradient (23 kpc < D_sun < 29 kpc), within the DES footprint in the Southern
hemisphere, using the same DES data, also indicate a possible association
between these systems. If these objects are confirmed through spectroscopic
followup to be gravitationally bound systems and to share a Galactic
trajectory with the Sagittarius stream, DES J01111341 and DES J0225+0304 would
be the first ultrafaint stellar systems associated with the Sagittarius stream.
Furthermore, DES J0225+0304 would also be the first confirmed case of an
ultrafaint satellite of a satellite.