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

Clusters of galaxies gravitationally lens the cosmic microwave background
(CMB) radiation, resulting in a distinct imprint in the CMB on arcminute
scales. Measurement of this effect offers a promising way to constrain the
masses of galaxy clusters, particularly those at high redshift. We use CMB maps
from the South Pole Telescope SunyaevZel'dovich (SZ) survey to measure the CMB
lensing signal around galaxy clusters identified in optical imaging from first
year observations of the Dark Energy Survey. The cluster catalog used in this
analysis contains 3697 members with mean redshift of $\bar{z} = 0.45$. We
detect lensing of the CMB by the galaxy clusters at $8.1\sigma$ significance.
Using the measured lensing signal, we constrain the amplitude of the relation
between cluster mass and optical richness to roughly $17\%$ precision, finding
good agreement with recent constraints obtained with galaxy lensing. The error
budget is dominated by statistical noise but includes significant contributions
from systematic biases due to the thermal SZ effect and cluster miscentering.

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.

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.

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

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.

The detection of GW170817 in both gravitational waves and electromagnetic
waves heralds the age of gravitationalwave multimessenger astronomy. On 17
August 2017 the Advanced LIGO and Virgo detectors observed GW170817, a strong
signal from the merger of a binary neutronstar system. Less than 2 seconds
after the merger, a gammaray burst (GRB 170817A) was detected within a region
of the sky consistent with the LIGOVirgoderived location of the
gravitationalwave source. This sky region was subsequently observed by optical
astronomy facilities, resulting in the identification of an optical transient
signal within $\sim 10$ arcsec of the galaxy NGC 4993. These multimessenger
observations allow us to use GW170817 as a standard siren, the
gravitationalwave analog of an astronomical standard candle, to measure the
Hubble constant. This quantity, which represents the local expansion rate of
the Universe, sets the overall scale of the Universe and is of fundamental
importance to cosmology. Our measurement combines the distance to the source
inferred purely from the gravitationalwave signal with the recession velocity
inferred from measurements of the redshift using electromagnetic data. This
approach does not require any form of cosmic "distance ladder;" the
gravitational wave analysis can be used to estimate the luminosity distance out
to cosmological scales directly, without the use of intermediate astronomical
distance measurements. We determine the Hubble constant to be
$70.0^{+12.0}_{8.0} \, \mathrm{km} \, \mathrm{s}^{1} \, \mathrm{Mpc}^{1}$
(maximum a posteriori and 68% credible interval). This is consistent with
existing measurements, while being completely independent of them. Additional
standardsiren measurements from future gravitationalwave sources will provide
precision constraints of this important cosmological parameter.

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

Galaxy crosscorrelations with highfidelity redshift samples hold the
potential to precisely calibrate systematic photometric redshift uncertainties
arising from the unavailability of complete and representative training and
validation samples of galaxies. However, application of this technique in the
Dark Energy Survey (DES) is hampered by the relatively low number density,
small area, and modest redshift overlap between photometric and spectroscopic
samples. We propose instead using photometric catalogs with reliable
photometric redshifts for photoz calibration via crosscorrelations. We verify
the viability of our proposal using redMaPPer clusters from the Sloan Digital
Sky Survey (SDSS) to successfully recover the redshift distribution of SDSS
spectroscopic galaxies. We demonstrate how to combine photoz with
crosscorrelation data to calibrate photometric redshift biases while
marginalizing over possible clustering bias evolution in either the calibration
or unknown photometric samples. We apply our method to DES Science Verification
(DES SV) data in order to constrain the photometric redshift distribution of a
galaxy sample selected for weak lensing studies, constraining the mean of the
tomographic redshift distributions to a statistical uncertainty of $\Delta z
\sim \pm 0.01$. We forecast that our proposal can in principle control
photometric redshift uncertainties in DES weak lensing experiments at a level
near the intrinsic statistical noise of the experiment over the range of
redshifts where redMaPPer clusters are available. Our results provide strong
motivation to launch a program to fully characterize the systematic errors from
bias evolution and photoz shapes in our calibration procedure.

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.

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.

We report the results of a Dark Energy Camera (DECam) optical followup of
the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO
detectors. Our observations cover 28.8 deg$^2$ of the localization region in
the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability),
starting 10 hours after the event was announced and spanning four epochs at
$224$ days after the GW detection. We achieve $5\sigma$ pointsource limiting
magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag,
in our difference images. Given the two day delay, we search this area for a
rapidly declining optical counterpart with $\gtrsim 3\sigma$ significance
steady decline between the first and final observations. We recover four
sources that pass our selection criteria, of which three are cataloged AGN. The
fourth source is offset by $5.8$ arcsec from the center of a galaxy at a
distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and
has a red color of $iz\approx 0.3$ mag. These properties roughly match the
expectations for a kilonova. However, this source was detected several times,
starting $94$ days prior to GW151226, in the PanSTARRS Survey for Transients
(dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its
longterm behavior, PS15cdi is likely a Type IIP supernova that transitioned
out of its plateau phase during our observations, mimicking a kilonovalike
behavior. We comment on the implications of this detection for contamination in
future optical followup observations.