
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.

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.

Our goal is to characterise the dependence of the optical masstolight ratio
on galaxy colour up to z = 1.5, expanding the redshift range explored in
previous work. From the ALHAMBRA redshifts, stellar masses, and restframe
luminosities provided by the MUFFIT code, we derive the masstolight ratio vs.
colour relation (MLCR) both for quiescent and starforming galaxies. The
intrinsic relation and its physical dispersion are derived with a Bayesian
inference model. The restframe iband masstolight ratio of quiescent and
starforming galaxies presents a tight correlation with the restframe (g  i)
colour up to z = 1.5. Such MLCR is linear for quiescent galaxies and quadratic
for starforming galaxies. The intrinsic dispersion in these relations is 0.02
dex for quiescent galaxies and 0.06 dex for starforming ones. The derived
MLCRs do not present a significant redshift evolution and are compatible with
previous local results in the literature. Finally, these tight relations also
hold for g and rband luminosities. The derived MLCRs in ALHAMBRA can be used
to predict the masstolight ratio from a restframe optical colour up to z =
1.5. These tight correlations do not change with redshift, suggesting that
galaxies have evolved along the derived relations during the last 9 Gyr.

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

We present the results of a search for rapidly evolving transients in the
Dark Energy Survey Supernova Programme. These events are characterized by fast
light curve evolution (rise to peak in $\lesssim 10$ d and exponential decline
in $\lesssim30$ d after peak). We discovered 72 events, including 37 transients
with a spectroscopic redshift from host galaxy spectral features. The 37 events
increase the total number of rapid optical transients by more than factor of
two. They are found at a wide range of redshifts ($0.05<z<1.56$) and peak
brightnesses ($15.75>M_\mathrm{g}>22.25$). The multiband photometry is well
fit by a blackbody up to few weeks after peak. The events appear to be hot
($T\approx1000030000$ K) and large ($R\approx 10^{14}2\cdot10^{15}$ cm) at
peak, and generally expand and cool in time, though some events show evidence
for a receding photosphere with roughly constant temperature. Spectra taken
around peak are dominated by a blue featureless continuum consistent with hot,
optically thick ejecta. We compare our events with a previously suggested
physical scenario involving shock breakout in an optically thick wind
surrounding a corecollapse supernova (CCSNe), we conclude that current models
for such a scenario might need an additional power source to describe the
exponential decline. We find these transients tend to favor starforming host
galaxies, which could be consistent with a corecollapse origin. However, more
detailed modeling of the light curves is necessary to determine their physical
origin.

Context. Knowing the exact shape of the UV luminosity function of
highredshift galaxies is important in order to understand the star formation
history of the early universe. However, the uncertainties, especially at the
faint and bright ends of the LFs, are still significant.
Aims. In this paper, we study the UV luminosity function of redshift z = 2.5
 4.5 galaxies in 2.38 deg^2 of ALHAMBRA data with I <= 24. Thanks to the large
area covered by ALHAMBRA, we particularly constrain the bright end of the
luminosity function. We also calculate the cosmic variance and the
corresponding bias values for our sample and derive their host dark matter halo
masses.
Methods. We use a novel methodology based on redshift and magnitude
probability distribution functions (PDFs). This methodology robustly takes into
account the uncertainties due to redshift and magnitude errors, shot noise and
cosmic variance, and models the luminosity function in two dimensions (z; M_UV
).
Results. We find an excess of bright ~ M*_UV galaxies as compared to the
studies based on broadband photometric data. However, our results agree well
with the luminosity function of the magnitudeselected spectroscopic VVDS data.
We measure high bias values, b ~ 8  10, that are compatible with the previous
measurements considering the redshifts and magnitudes of our galaxies and
further reinforce the real highredshift nature of our bright galaxies.
Conclusions. We call into question the shape of the luminosity function at
its bright end; is it a double powerlaw as suggested by the recent broadband
photometric studies or rather a brighter Schechter function, as suggested by
our multifilter analysis and the spectroscopic VVDS data.

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 study the theoretical systematics and optimize the methodology in Baryon
Acoustic Oscillations (BAO) detections using the angular correlation function
with tomographic bins. We calibrate and optimize the pipeline for the Dark
Energy Survey Year 1 dataset using 1800 mocks. We compare the BAO fitting
results obtained with three estimators: the Maximum Likelihood Estimator (MLE),
Profile Likelihood, and Markov Chain Monte Carlo. The MLE method yields the
least bias in the fit results (bias/spread $\sim 0.02$) and the error bar
derived is the closest to the Gaussian results (1% from 68% Gaussian
expectation). When there is mismatch between the template and the data either
due to incorrect fiducial cosmology or photo$z$ error, the MLE again gives the
leastbiased results. The BAO angular shift that is estimated based on the
sound horizon and the angular diameter distance agree with the numerical fit.
Various analysis choices are further tested: the number of redshift bins,
crosscorrelations, and angular binning. We propose two methods to correct the
mock covariance when the final sample properties are slightly different from
those used to create the mock. We show that the sample changes can be
accommodated with the help of the Gaussian covariance matrix or more
effectively using the eigenmode expansion of the mock covariance. The eigenmode
expansion is significantly less susceptible to statistical fluctuations
relative to the direct measurements of the covariance matrix because the number
of free parameters is substantially reduced [$p$ parameters versus $p(p+1)/2$
from direct measurement].

We perform a search for stellar streams around the Milky Way using the first
three years of multiband optical imaging data from the Dark Energy Survey
(DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of $g > 23.5$
with a relative photometric calibration uncertainty of $< 1 \%$. This data set
yields unprecedented sensitivity to the stellar density field in the southern
celestial hemisphere, enabling the detection of faint stellar streams to a
heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a
matchedfilter in colormagnitude space derived from a synthetic isochrone of
an old, metalpoor stellar population. Our detection technique recovers four
previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a
possible extension of Molonglo. In addition, we report the discovery of eleven
new stellar streams. In general, the new streams detected by DES are fainter,
more distant, and lower surface brightness than streams detected by similar
techniques in previous photometric surveys. As a byproduct of our stellar
stream search, we find evidence for extratidal stellar structure associated
with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The
evergrowing sample of stellar streams will provide insight into the formation
of the Galactic stellar halo, the Milky Way gravitational potential, as well as
the large and smallscale distribution of dark matter around the Milky Way.

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 present observations of DES16C2nm, the first spectroscopically confirmed
hydrogenfree superluminous supernova (SLSNI) at redshift z~2. DES16C2nm was
discovered by the Dark Energy Survey (DES) Supernova Program, with followup
photometric data from the Hubble Space Telescope, Gemini, and the European
Southern Observatory Very Large Telescope supplementing the DES data.
Spectroscopic observations confirm DES16C2nm to be at z=1.998, and
spectroscopically similar to Gaia16apd (a SLSNI at z=0.102), with a peak
absolute magnitude of U=22.26$\pm$0.06. The high redshift of DES16C2nm
provides a unique opportunity to study the ultraviolet (UV) properties of
SLSNeI. Combining DES16C2nm with ten similar events from the literature, we
show that there exists a homogeneous class of SLSNeI in the UV (~2500A), with
peak luminosities in the (restframe) U band, and increasing absorption to
shorter wavelengths. There is no evidence that the mean photometric and
spectroscopic properties of SLSNeI differ between low (z<1) and high redshift
(z>1), but there is clear evidence of diversity in the spectrum at <2000A,
possibly caused by the variations in temperature between events. No significant
correlations are observed between spectral line velocities and photometric
luminosity. Using these data, we estimate that SLSNeI can be discovered to
z=3.8 by DES. While SLSNeI are typically identified from their blue observed
colors at low redshift (z<1), we highlight that at z>2 these events appear
optically red, peaking in the observerframe zband. Such characteristics are
critical to identify these objects with future facilities such as the Large
Synoptic Survey Telescope, Euclid, and the WideField Infrared Survey
Telescope, which should detect such SLSNeI to z=3.5, 3.7, and 6.6,
respectively.

We present a study of NGC 4993, the host galaxy of the GW170817 gravitational
wave event, the GRB170817A short gammaray burst (sGRB) and the AT2017gfo
kilonova. We use Dark Energy Camera imaging, AAT spectra and publicly available
data, relating our findings to binary neutron star (BNS) formation scenarios
and merger delay timescales. NGC4993 is a nearby (40 Mpc) earlytype galaxy,
with $i$band S\'ersic index $n=4.0$ and low asymmetry ($A=0.04\pm 0.01$).
These properties are unusual for sGRB hosts. However, NGC4993 presents
shelllike structures and dust lanes indicative of a recent galaxy merger, with
the optical transient located close to a shell. We constrain the star formation
history (SFH) of the galaxy assuming that the galaxy merger produced a star
formation burst, but find little to no ongoing star formation in either
spatiallyresolved broadband SED or spectral fitting. We use the bestfit SFH
to estimate the BNS merger rate in this type of galaxy, as $R_{NSM}^{gal}=
5.7^{+0.57}_{3.3} \times 10^{6} {\rm yr}^{1}$. If star formation is the only
considered BNS formation scenario, the expected number of BNS mergers from
earlytype galaxies detectable with LIGO during its first two observing seasons
is $0.038^{+0.004}_{0.022}$, as opposed to $\sim 0.5$ from all galaxy types.
Hypothesizing that the binary system formed due to dynamical interactions
during the galaxy merger, the subsequent time elapsed can constrain the delay
time of the BNS coalescence. By using velocity dispersion estimates and the
position of the shells, we find that the galaxy merger occurred $t_{\rm
mer}\lesssim 200~{\rm Myr}$ prior to the BNS coalescence.

The discovery of a kilonova (KN) associated with the Advanced LIGO
(aLIGO)/Virgo event GW170817 opens up new avenues of multimessenger
astrophysics. Here, using realistic simulations, we provide estimates of the
number of KNe that could be found in data from past, present and future surveys
without a gravitationalwave trigger. For the simulation, we construct a
spectral timeseries model based on the DESGW multiband lightcurve from the
single known KN event, and we use an average of BNS rates from past studies of
$10^3 \rm{Gpc}^{3}/\rm{year}$, consistent with the $1$ event found so far.
Examining past and current datasets from transient surveys, the number of KNe
we expect to find for ASASSN, SDSS, PS1, SNLS, DES, and SMT is between 0 and
$0.3$. We predict the number of detections per future survey to be: 8.3 from
ATLAS, 10.6 from ZTF, 5.5/69 from LSST (the Deep Drilling / Wide Fast Deep),
and 16.0 from WFIRST. The maximum redshift of KNe discovered for each survey is
z = 0.8 for WFIRST, z = 0.25 for LSST and z = 0.04 for ZTF and ATLAS. For the
LSST survey, we also provide contamination estimates from Type Ia and
Corecollapse supernovae: after lightcurve and templatematching requirements,
we estimate a background of just 2 events. More broadly, we stress that future
transient surveys should consider how to optimize their search strategies to
improve their detection efficiency, and to consider similar analyses for GW
followup programs.

We present the Dark Energy Camera (DECam) discovery of the optical
counterpart of the first binary neutron star merger detected through
gravitational wave emission, GW170817. Our observations commenced 10.5 hours
postmerger, as soon as the localization region became accessible from Chile.
We imaged 70 deg$^2$ in the $i$ and $z$ bands, covering 93\% of the initial
integrated localization probability, to a depth necessary to identify likely
optical counterparts (e.g., a kilonova). At 11.4 hours postmerger we detected
a bright optical transient located $10.6''$ from the nucleus of NGC\,4993 at
redshift $z=0.0098$, consistent (for $H_0 = 70$\, km s$^{1}$ Mpc$^{1}$) with
the distance of $40 \pm 8$\, Mpc reported by the LIGO Scientific Collaboration
and the Virgo Collaboration (LVC). At detection the transient had magnitudes
$i\approx 17.30$ and $z\approx 17.45$, and thus an absolute magnitude of $M_i =
15.7$, in the luminosity range expected for a kilonova. We identified 1,500
potential transient candidates. Applying simple selection criteria aimed at
rejecting background events such as supernovae, we find the transient
associated with NGC\,4993 as the only remaining plausible counterpart, and
reject chance coincidence at the 99.5\% confidence level. We therefore conclude
that the optical counterpart we have identified near NGC\,4993 is associated
with GW170817. This discovery ushers in the era of multimessenger astronomy
with gravitational waves, and demonstrates the power of DECam to identify the
optical counterparts of gravitationalwave sources.

We present UV, optical, and NIR photometry of the first electromagnetic
counterpart to a gravitational wave source from Advanced LIGO/Virgo, the binary
neutron star merger GW170817. Our data set extends from the discovery of the
optical counterpart at $0.47$ days to $18.5$ days postmerger, and includes
observations with the Dark Energy Camera (DECam), GeminiSouth/FLAMINGOS2
(GS/F2), and the {\it Hubble Space Telescope} ({\it HST}). The spectral energy
distribution (SED) inferred from this photometry at $0.6$ days is well
described by a blackbody model with $T\approx 8300$ K, a radius of $R\approx
4.5\times 10^{14}$ cm (corresponding to an expansion velocity of $v\approx
0.3c$), and a bolometric luminosity of $L_{\rm bol}\approx 5\times10^{41}$ erg
s$^{1}$. At $1.5$ days we find a multicomponent SED across the optical and
NIR, and subsequently we observe rapid fading in the UV and blue optical bands
and significant reddening of the optical/NIR colors. Modeling the entire data
set we find that models with heating from radioactive decay of $^{56}$Ni, or
those with only a single component of opacity from $r$process elements, fail
to capture the rapid optical decline and red optical/NIR colors. Instead,
models with two components consistent with lanthanidepoor and lanthaniderich
ejecta provide a good fit to the data, the resulting "blue" component has
$M_\mathrm{ej}^\mathrm{blue}\approx 0.01$ M$_\odot$ and
$v_\mathrm{ej}^\mathrm{blue}\approx 0.3$c, and the "red" component has
$M_\mathrm{ej}^\mathrm{red}\approx 0.04$ M$_\odot$ and
$v_\mathrm{ej}^\mathrm{red}\approx 0.1$c. These ejecta masses are broadly
consistent with the estimated $r$process production rate required to explain
the Milky Way $r$process abundances, providing the first evidence that BNS
mergers can be a dominant site of $r$process enrichment.

Using $\sim 100$ Xray selected clusters in the Dark Energy Survey Science
Verification data, we constrain the luminosity function (LF) of cluster red
sequence galaxies as a function of redshift. This is the first homogeneous
optical/Xray sample large enough to constrain the evolution of the luminosity
function simultaneously in redshift ($0.1<z<1.05$) and cluster mass ($13.5 \le
\rm{log_{10}}(M_{200crit}) \sim< 15.0$). We pay particular attention to
completeness issues and the detection limit of the galaxy sample. We then apply
a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schecter
function, including its characteristic break ($m^*$) to a faint end powerlaw
slope ($\alpha$). Our method enables us to avoid known issues in similar
analyses based on stacking or binning the clusters. We find weak and
statistically insignificant ($\sim 1.9 \sigma$) evolution in the faint end
slope $\alpha$ versus redshift. We also find no dependence in $\alpha$ or $m^*$
with the Xray inferred cluster masses. However, the amplitude of the LF as a
function of cluster mass is constrained to $\sim 20\%$ precision. As a
byproduct of our algorithm, we utilize the correlation between the LF and
cluster mass to provide an improved estimate of the individual cluster masses
as well as the scatter in true mass given the Xray inferred masses. This
technique can be applied to a larger sample of Xray or optically selected
clusters from the Dark Energy Survey, significantly improving the sensitivity
of the analysis.