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

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

We present density split statistics, a framework that studies lensing and
countsincells as a function of foreground galaxy density, thereby providing a
largescale measurement of both 2point and 3point statistics. Our method
extends our earlier work on trough lensing and is summarized as follows: given
a foreground (low redshift) population of galaxies, we divide the sky into
subareas of equal size but distinct galaxy density. We then measure lensing
around uniformly spaced points separately in each of these subareas, as well as
countsincells statistics (CiC). The lensing signals trace the matter density
contrast around regions of fixed galaxy density. Through the CiC measurements
this can be related to the density profile around regions of fixed matter
density. Together, these measurements constitute a powerful probe of cosmology,
the skewness of the density field and the connection of galaxies and matter.
In this paper we show how to model both the density split lensing signal and
CiC from basic ingredients: a nonlinear power spectrum, clustering hierarchy
coefficients from perturbation theory and a parametric model for galaxy bias
and shotnoise. Using Nbody simulations, we demonstrate that this model is
sufficiently accurate for a cosmological analysis on year 1 data from the Dark
Energy Survey.

Aims. We develop, test and characterise of a new statistical tool
(intelligent system) for the sifting and analysis of nearby young open cluster
(NYOC) populations. Methods. Using a Bayesian formalism, this statistical tool
is able to obtain the posterior distributions of parameters governing the
cluster model. It also uses hierarchical bayesian models to establish weakly
informative priors, and incorporates the treatment of missing values and
nonhomogeneous (heteroscedastic) observational uncertainties. Results. From
simulations, we estimate that this statistical tool renders kinematic (proper
motion) and photometric (luminosity) distributions of the cluster population
with a contamination rate of $5.8 \pm 0.2$ %. The luminosity distributions and
present day mass function agree with the ones found by Bouy et al. (2015b) on
the completeness interval of the survey. At the probability threshold of
maximum accuracy, the classifier recovers $\sim$ 90% of Bouy et al. (2015b)
candidate members and finds 10% of new ones. Conclusions. A new statistical
tool for the analysis of NYOC is introduced, tested and characterised. Its
comprehensive modelling of the data properties allows it to get rid of the
biases present in previous works. In particular, those resulting from the use
of only completely observed (nonmissing) data and the assumption of
homoskedastic uncertainties. Also, its Bayesian framework allows it to properly
propagate observational uncertainties into membership probabilities and cluster
velocity and luminosity distributions. Our results are in a general agreement
with those from the literature, although we provide the most uptodate and
extended list of candidate members of the Pleiades cluster.

Methods. We compute Bayesian evidences and Bayes Factors for a set of
variations of the classical radial models by King (1962), Elson et al. (1987)
and Lauer et al. (1995). The variations incorporate different degrees of model
freedom and complexity, amongst which we include biaxial (elliptical) symmetry,
and luminosity segregation. As a byproduct of the model comparison, we obtain
posterior distributions and maximum a posteriori estimates for each set of
model parameters. Results. We find that the model comparison results depend on
the spatial extent of the region used for the analysis. For a circle of 11.5
parsecs around the cluster centre (the most homogeneous and complete region),
we find no compelling reason to abandon Kings model, although the Generalised
King model, introduced in this work, has slightly better fitting properties.
Furthermore, we find strong evidence against radially symmetric models when
compared to the elliptic extensions. Finally, we find that including mass
segregation in the form of luminosity segregation in the J band, is strongly
supported in all our models. Conclusions. We have put the question of the
projected spatial distribution of the Pleiades cluster on a solid probabilistic
framework, and inferred its properties using the most exhaustive and least
contaminated list of Pleiades candidate members available to date. Our results
suggest however that this sample may still lack about 20% of the expected
number of cluster members. Therefore, this study should be revised when the
completeness and homogeneity of the data can be extended beyond the 11.5
parsecs limit. Such study will allow a more precise determination of the
Pleiades spatial distribution, its tidal radius, ellipticity, number of objects
and total mass.

We present a spectroscopic study of the tidal tails and core of the Milky Way
satellite Tucana III, collectively referred to as the Tucana III stream, using
the 2dF+AAOmega spectrograph on the AngloAustralian Telescope and the IMACS
spectrograph on the Magellan/Baade Telescope. In addition to recovering the
brightest 9 previously known member stars in the Tucana III core, we identify
22 members in the tidal tails. We observe strong evidence for a velocity
gradient of 8.0 km/s/deg (or 18.3 km/s/kpc over at least 3$^{\circ}$ (or 1.3
kpc) on the sky. Based on the continuity in velocity we confirm that the Tucana
III tails are real tidal extensions of Tucana III. The large velocity gradient
of the stream implies that Tucana III is likely on a radial orbit. We
successfully obtain metallicities for 4 members in the core and 12 members in
the tails. We find that members close to the ends of the stream tend to be more
metalpoor than members in the core, indicating a possible metallicity gradient
between the center of the progenitor halo and its edge. The spread in
metallicity suggests that the progenitor of the Tucana III stream is a dwarf
galaxy rather than a star cluster. Furthermore, we find that with the precise
photometry of the Dark Energy Survey data, there is a discernible color offset
between metalrich disk stars and metalpoor stream members. This
metallicitydependent color offers a more efficient method to recognize
metalpoor targets and will increase the selection efficiency of stream members
for future spectroscopic followup programs on stellar streams.

We present results of the first dynamical stream fits to the recently
discovered Tucana III stream. These fits assume a fixed Milky Way potential and
give proper motion predictions, which can be tested with the upcoming Gaia Data
Release 2. These fits reveal that Tucana III is on an eccentric orbit around
the Milky Way and, more interestingly, that Tucana III passed within 15 kpc of
the Large Magellanic Cloud (LMC) approximately 75 Myr ago. Given this close
passage, we fit the Tucana III stream in the combined presence of the Milky Way
and the LMC. We find that the predicted proper motions depend on the assumed
mass of the LMC and that the LMC can induce a substantial proper motion
perpendicular to the stream track. A detection of this misalignment will
directly probe the extent of the LMC's influence on our Galaxy, and has
implications for nearly all methods which attempt to constraint the Milky Way
potential. Such a measurement will be possible with the upcoming Gaia DR2,
allowing for a measurement of the LMC's mass.

We study the effect of introducing altruistic agents in a Schellinglike
model of residential segregation. We find that even an infinitesimal proportion
of altruists has dramatic catalytic effects on the collective utility of the
system. Altruists provide pathways that move the system away from the
suboptimal equilibrium it would reach if the system included only egoist
agents, allowing it to reach the optimal steady state.

We introduce and study in two dimensions a new class of dry, aligning, active
matter that exhibits a direct transition to orientational order, without the
phaseseparation phenomenology usually observed in this context. Characterized
by selfpropelled particles with velocity reversals and ferromagnetic alignment
of polarities, systems in this class display quasilongrange polar order with
continuouslyvarying scaling exponents and yet a numerical study of the
transition leads to conclude that it does not belong to the
BerezinskiiKosterlitzThouless universality class, but is best described as a
standard critical point with algebraic divergence of correlations. We
rationalize these findings by showing that the interplay between order and
density changes the role of defects.

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.

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

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

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