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

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

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

The luminosity distance measurement of GW170817 derived from GW analysis in
Abbott et al. 2017 (here, A17:H0) is highly correlated with the measured
inclination of the NSNS system. To improve the precision of the distance
measurement, we attempt to constrain the inclination by modeling the broadband
Xraytoradio emission from GW170817, which is dominated by the interaction of
the jet with the environment. We update our previous analysis and we consider
the radio and Xray data obtained at $t<40$ days since merger. We find that the
afterglow emission from GW170817 is consistent with an offaxis relativistic
jet with energy $10^{48}\,\rm{erg}<E_{k}\le 3\times 10^{50} \,\rm{erg}$
propagating into an environment with density $n\sim10^{2}10^{4}
\,\rm{cm^{3}}$, with preference for wider jets (opening angle $\theta_j=15$
deg). For these jets, our modeling indicates an offaxis angle $\theta_{\rm
obs}\sim2550$ deg. We combine our constraints on $\theta_{\rm obs}$ with the
joint distanceinclination constraint from LIGO. Using the same $\sim 170$
km/sec peculiar velocity uncertainty assumed in A17:H0 but with an inclination
constraint from the afterglow data, we get a value of $H_0=$$74.0 \pm
\frac{11.5}{7.5}$ $\mbox{km/s/Mpc}$, which is higher than the value of
$H_0=$$70.0 \pm \frac{12.0}{8.0}$ $\mbox{km/s/Mpc}$ found in A17:H0. Further,
using a more realistic peculiar velocity uncertainty of 250 km/sec derived from
previous work, we find $H_0=$$75.5 \pm \frac{11.6}{9.6}$ km/s/Mpc for H0 from
this system. We note that this is in modestly better agreement with the local
distance ladder than the Planck CMB, though a significant such discrimination
will require $\sim 50$ such events. Future measurements at $t>100$ days of the
Xray and radio emission will lead to tighter constraints.

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.

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.

Galaxies and their dark matter halos populate a complicated filamentary
network around large, nearly empty regions known as cosmic voids. Cosmic voids
are usually identified in spectroscopic galaxy surveys, where 3D information
about the largescale structure of the Universe is available. Although an
increasing amount of photometric data is being produced, its potential for void
studies is limited since photometric redshifts induce lineofsight position
errors of $\sim50$ Mpc/$h$ or more that can render many voids undetectable. In
this paper we present a new void finder designed for photometric surveys,
validate it using simulations, and apply it to the highquality photo$z$
redMaGiC galaxy sample of the Dark Energy Survey Science Verification (DESSV)
data. The algorithm works by projecting galaxies into 2D slices and finding
voids in the smoothed 2D galaxy density field of the slice. Fixing the
lineofsight size of the slices to be at least twice the photo$z$ scatter,
the number of voids found in these projected slices of simulated spectroscopic
and photometric galaxy catalogs is within 20% for all transverse void sizes,
and indistinguishable for the largest voids of radius $\sim 70$ Mpc/$h$ and
larger. The positions, radii, and projected galaxy profiles of photometric
voids also accurately match the spectroscopic void sample. Applying the
algorithm to the DESSV data in the redshift range $0.2<z<0.8$, we identify 87
voids with comoving radii spanning the range 18120 Mpc/$h$, and carry out a
stacked weak lensing measurement. With a significance of $4.4\sigma$, the
lensing measurement confirms the voids are truly underdense in the matter field
and hence not a product of Poisson noise, tracer density effects or systematics
in the data. It also demonstrates, for the first time in real data, the
viability of void lensing studies in photometric surveys.

Cosmic shear is sensitive to fluctuations in the cosmological matter density
field, including on small physical scales, where matter clustering is affected
by baryonic physics in galaxies and galaxy clusters, such as star formation,
supernovae feedback and AGN feedback. While muddying any cosmological
information that is contained in small scale cosmic shear measurements, this
does mean that cosmic shear has the potential to constrain baryonic physics and
galaxy formation. We perform an analysis of the Dark Energy Survey (DES)
Science Verification (SV) cosmic shear measurements, now extended to smaller
scales, and using the Mead et al. 2015 halo model to account for baryonic
feedback. While the SV data has limited statistical power, we demonstrate using
a simulated likelihood analysis that the final DES data will have the
statistical power to differentiate among baryonic feedback scenarios. We also
explore some of the difficulties in interpreting the small scales in cosmic
shear measurements, presenting estimates of the size of several other
systematic effects that make inference from small scales difficult, including
uncertainty in the modelling of intrinsic alignment on nonlinear scales,
`lensing bias', and shape measurement selection effects. For the latter two, we
make use of novel image simulations. While future cosmic shear datasets have
the statistical power to constrain baryonic feedback scenarios, there are
several systematic effects that require improved treatments, in order to make
robust conclusions about baryonic feedback.

We measure the correlation of galaxy lensing and cosmic microwave background
lensing with a set of galaxies expected to trace the matter density field. The
measurements are performed using presurvey Dark Energy Survey (DES) Science
Verification optical imaging data and millimeterwave data from the 2500 square
degree South Pole Telescope SunyaevZel'dovich (SPTSZ) survey. The two
lensinggalaxy correlations are jointly fit to extract constraints on
cosmological parameters, constraints on the redshift distribution of the lens
galaxies, and constraints on the absolute shear calibration of DES galaxy
lensing measurements. We show that an attractive feature of these fits is that
they are fairly insensitive to the clustering bias of the galaxies used as
matter tracers. The measurement presented in this work confirms that DES and
SPT data are consistent with each other and with the currently favored
$\Lambda$CDM cosmological model. It also demonstrates that joint lensinggalaxy
correlation measurement considered here contains a wealth of information that
can be extracted using current and future surveys.

We measure the crosscorrelation between the galaxy density in the Dark
Energy Survey (DES) Science Verification data and the lensing of the cosmic
microwave background (CMB) as reconstructed with the Planck satellite and the
South Pole Telescope (SPT). When using the DES main galaxy sample over the full
redshift range $0.2 < z < 1.2$, a crosscorrelation signal is detected at $6
\sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES
galaxies into five photometric redshift bins, finding significant ($>$$2
\sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we
find the redshift evolution of the signal matches expectations, although the
amplitude is consistently lower than predicted across redshift bins. We test
for possible systematics that could affect our result and find no evidence for
significant contamination. Finally, we demonstrate how these measurements can
be used to constrain the growth of structure across cosmic time. We find the
data are fit by a model in which the amplitude of structure in the $z<1.2$
universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck
cosmology, a $1.7\sigma$ deviation.

The American Physical Society's Division of Particles and Fields initiated a
longterm planning exercise over 201213, with the goal of developing the
community's long term aspirations. The subgroup "Dark Energy and CMB" prepared
a series of papers explaining and highlighting the physics that will be studied
with large galaxy surveys and cosmic microwave background experiments. This
paper summarizes the findings of the other papers, all of which have been
submitted jointly to the arXiv.

Large sets of objects with spectroscopic redshift measurements will be needed
for imaging dark energy experiments to achieve their full potential, serving
two goals:_training_, i.e., the use of objects with known redshift to develop
and optimize photometric redshift algorithms; and_calibration_, i.e., the
characterization of moments of redshift (or photoz error) distributions.
Better training makes cosmological constraints from a given experiment
stronger, while highlyaccurate calibration is needed for photoz systematics
not to dominate errors. In this white paper, we investigate the required scope
of spectroscopic datasets which can serve both these purposes for ongoing and
nextgeneration dark energy experiments, as well as the time required to obtain
such data with instruments available in the next decade. Large time allocations
on kiloobject spectrographs will be necessary, ideally augmented by infrared
spectroscopy from space. Alternatively, precision calibrations could be
obtained by measuring crosscorrelation statistics using samples of bright
objects from a large baryon acoustic oscillation experiment such as DESI. We
also summarize the additional work on photometric redshift methods needed to
prepare for ongoing and future dark energy experiments.

The anomalous lack of large angle temperature correlations has been a
surprising feature of the cosmic microwave background (CMB) since first
observed by COBEDMR and subsequently confirmed and strengthened by the
Wilkinson Microwave Anisotropy Probe. This anomaly may point to the need for
modifications of the standard model of cosmology or may indicate that our
Universe is a rare statistical fluctuation within that model. Further
observations of the temperature autocorrelation function will not elucidate
the issue; sufficiently high precision statistical observations already exist.
Instead, alternative probes are required. In this work we explore the
expectations for forthcoming polarization observations. We define a
prescription to test the hypothesis that the largeangle CMB temperature
perturbations in our Universe represent a rare statistical fluctuation within
the standard cosmological model. These tests are based on the temperature$Q$
Stokes parameter correlation. Unfortunately these tests cannot be expected to
be definitive. However, we do show that if this $TQ$correlation is observed to
be sufficiently large over an appropriately chosen angular range, then the
hypothesis can be rejected at a high confidence level. We quantify these
statements and optimize the statistics we have constructed to apply to the
anticipated polarization data. We find that we can construct a statistic that
has a 25 per cent chance of excluding the hypothesis that we live in a rare
realization of LCDM at the 99.9 per cent confidence level.

BigBOSS is a Stage IV groundbased dark energy experiment to study baryon
acoustic oscillations (BAO) and the growth of structure with a widearea galaxy
and quasar redshift survey over 14,000 square degrees. It has been
conditionally accepted by NOAO in response to a call for major new
instrumentation and a highimpact science program for the 4m Mayall telescope
at Kitt Peak. The BigBOSS instrument is a roboticallyactuated, fiberfed
spectrograph capable of taking 5000 simultaneous spectra over a wavelength
range from 340 nm to 1060 nm, with a resolution R = 30004800.
Using data from imaging surveys that are already underway, spectroscopic
targets are selected that trace the underlying dark matter distribution. In
particular, targets include luminous red galaxies (LRGs) up to z = 1.0,
extending the BOSS LRG survey in both redshift and survey area. To probe the
universe out to even higher redshift, BigBOSS will target bright [OII] emission
line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are
obtained to measure the BAO feature, trace the matter power spectrum at smaller
scales, and detect redshift space distortions. BigBOSS will provide additional
constraints on early dark energy and on the curvature of the universe by
measuring the Lyalpha forest in the spectra of over 600,000 2.2 < z < 3.5
quasars.
BigBOSS galaxy BAO measurements combined with an analysis of the broadband
power, including the Lyalpha forest in BigBOSS quasar spectra, achieves a FOM
of 395 with Planck plus Stage III priors. This FOM is based on conservative
assumptions for the analysis of broad band power (kmax = 0.15), and could grow
to over 600 if current work allows us to push the analysis to higher wave
numbers (kmax = 0.3). BigBOSS will also place constraints on theories of
modified gravity and inflation, and will measure the sum of neutrino masses to
0.024 eV accuracy.

We extend the Figure of Merit formalism usually adopted to quantify the
statistical performance of future dark energy probes to assess the robustness
of a future mission to plausible systematic bias. We introduce a new robustness
Figure of Merit which can be computed in the Fisher Matrix formalism given
arbitrary systematic biases in the observable quantities. We argue that
robustness to systematics is an important new quantity that should be taken
into account when optimizing future surveys. We illustrate our formalism with
toy examples, and apply it to future type Ia supernova (SNIa) and baryonic
acoustic oscillation (BAO) surveys. For the simplified systematic biases that
we consider, we find that SNIa are a somewhat more robust probe of dark energy
parameters than the BAO. We trace this back to a geometrical alignement of
systematic bias direction with statistical degeneracy directions in the dark
energy parameter space.

We propose a novel representation of cosmic microwave anisotropy maps, where
each multipole order l is represented by l unit vectors pointing in directions
on the sky and an overall magnitude. These "multipole vectors and scalars"
transform as vectors under rotations. Like the usual spherical harmonics,
multipole vectors form an irreducible representation of the proper rotation
group SO(3). However, they are related to the familiar spherical harmonic
coefficients, alm, in a nonlinear way, and are therefore sensitive to different
aspects of the CMB anisotropy. Nevertheless, it is straightforward to determine
the multipole vectors for a given CMB map and we present an algorithm to
compute them. Using the WMAP fullsky maps, we perform several tests of the
hypothesis that the CMB anisotropy is statistically isotropic and Gaussian
random. We find that the result from comparing the oriented area of planes
defined by these vectors between multipole pairs 2<=l1!=l2<=8 is inconsistent
with the isotropic Gaussian hypothesis at the 99.4% level for the ILC map and
at 98.9% level for the cleaned map of Tegmark et al. A particular correlation
is suggested between the l=3 and l=8 multipoles, as well as several other
pairs. This effect is entirely different from the now familiar planarity and
alignment of the quadrupole and octupole: while the aforementioned is fairly
unlikely, the multipole vectors indicate correlations not expected in Gaussian
random skies that make them unusually likely. The result persists after
accounting for pixel noise and after assuming a residual 10% dust contamination
in the cleaned WMAP map. While the definitive analysis of these results will
require more work, we hope that multipole vectors will become a valuable tool
for various cosmological tests, in particular those of cosmic isotropy.

The Supernova / Acceleration Probe (SNAP) is a proposed spacebased
experiment designed to study the dark energy and alternative explanations of
the acceleration of the Universe's expansion by performing a series of
complementary systematicscontrolled measurements. We describe a
selfconsistent reference mission design for building a Type Ia supernova
Hubble diagram and for performing a widearea weak gravitational lensing study.
A 2m widefield telescope feeds a focal plane consisting of a 0.7
squaredegree imager tiled with equal areas of optical CCDs and near infrared
sensors, and a highefficiency lowresolution integral field spectrograph. The
SNAP mission will obtain highsignaltonoise calibrated lightcurves and
spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A
widefield survey covering one thousand square degrees resolves ~100 galaxies
per square arcminute. If we assume we live in a cosmologicalconstantdominated
Universe, the matter density, dark energy density, and flatness of space can
all be measured with SNAP supernova and weaklensing measurements to a
systematicslimited accuracy of 1%. For a flat universe, the
densitytopressure ratio of dark energy can be similarly measured to 5% for
the present value w0 and ~0.1 for the time variation w'. The large survey area,
depth, spatial resolution, timesampling, and nineband optical to NIR
photometry will support additional independent and/or complementary darkenergy
measurement approaches as well as a broad range of auxiliary science programs.
(Abridged)