Obtaining accurate distributions of galaxy redshifts is a critical aspect of weak lensing cosmology experiments. One of the methods used to estimate and validate redshift distributions is apply weights to a spectroscopic sample so that their weighted photometry distribution matches the target sample. In this work we estimate the \textit{selection bias} in redshift that is introduced in this procedure. We do so by simulating the process of assembling a spectroscopic sample (including observer-assigned confidence flags) and highlight the impacts of spectroscopic target selection and redshift failures. We use the first year (Y1) weak lensing analysis in DES as an example data set but the implications generalise to all similar weak lensing surveys. We find that using colour cuts that are not available to the weak lensing galaxies can introduce biases of $\Delta~z\sim0.015$ in the weighted mean redshift of different redshift intervals. To assess the impact of incompleteness in spectroscopic samples, we select only objects with high observer-defined confidence flags and compare the weighted mean redshift with the true mean. We find that the mean redshift of the DES Y1 weak lensing sample is typically biased at the $\Delta~z=0.005-0.05$ level after the weighting is applied. The bias we uncover can have either sign, depending on the samples and redshift interval considered. For the highest redshift bin, the bias is larger than the uncertainties in the other DES Y1 redshift calibration methods, justifying the decision of not using this method for the redshift estimations. We discuss several methods to mitigate this bias.

Cosmological studies of large-scale structure have relied on two-point statistics, not fully exploiting the rich structure of the cosmic web. In this paper we show how to capture some of this cosmic web information by using the minimum spanning tree (MST), for the first time using it to estimate cosmological parameters in simulations. Discrete tracers of dark matter such as galaxies, $N$-body particles or haloes are used as nodes to construct a unique graph, the MST, that traces skeletal structure. We study the dependence of the MST on cosmological parameters using haloes from a suite of COLA simulations with a box size of $250\ h^{-1}{\rm Mpc}$, varying the amplitude of scalar fluctuations $\left(A_{\rm s}\right)$, matter density $\left(\Omega_{\rm m}\right)$, and neutrino mass $\left(\sum m_{\nu}\right)$. The power spectrum $P$ and bispectrum $B$ are measured for wavenumbers between $0.125$ and $0.5$ $h{\rm Mpc}^{-1}$, while a corresponding lower cut of $\sim12.6$ $h^{-1}{\rm Mpc}$ is applied to the MST. The constraints from the individual methods are fairly similar but when combined we see improved $1\sigma$ constraints of $\sim 17\%$ ($\sim 12\%$) on $\Omega_{\rm m}$ and $\sim 12\%$ ($\sim 10\%$) on $A_{\rm s}$ with respect to $P$ ($P+B$) thus showing the MST is providing additional information. The MST can be applied to current and future spectroscopic surveys (BOSS, DESI, Euclid, PSF, WFIRST, and 4MOST) in 3D and photometric surveys (DES and LSST) in tomographic shells to constrain parameters and/or test systematics.

We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters. Our forecast shows that the combination of second and third moments provides a 1.5 percent constraint on $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ for DES Y3 data. This is 20 percent better than an analysis using a simulated DES Y3 shear 2-point statistics, owing to the non-Gaussian information captured by the inclusion of higher-order statistics. The methodology developed here can be applied to current and future weak lensing datasets to make use of information from non-Gaussianity in the cosmic density field and to improve constraints on cosmological parameters.

Cosmic voids gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint on degree scales. We use the simulated CMB lensing convergence map from the MICE N-body simulation to calibrate our detection strategy for a given void definition and galaxy tracer density. We then identify cosmic voids in DES Year 1 data and stack the Planck 2015 lensing convergence map on their locations, probing the consistency of simulated and observed void lensing signals. When fixing the shape of the stacked convergence profile to that calibrated from simulations, we find imprints at the $3{\sigma}$ significance level for various analysis choices. The best measurement strategies based on the MICE calibration process yield $S/N \sim 4$ for DES Y1, and the best-fit amplitude recovered from the data is consistent with expectations from MICE ($A \sim 1$). Given these results as well as the agreement between them and N-body simulations, we conclude that the previously reported excess integrated Sachs-Wolfe (ISW) signal associated with cosmic voids in DES Y1 has no counterpart in the Planck CMB lensing map.

We propose a principled Bayesian method for quantifying tension between correlated datasets with wide uninformative parameter priors. This is achieved by extending the Suspiciousness statistic, which is insensitive to priors. Our method uses global summary statistics, and as such it can be used as a diagnostic for internal consistency. We show how our approach can be combined with methods that use parameter space and data space to identify the existing internal discrepancies. As an example, we use it to test the internal consistency of the KiDS-450 data in 4 photometric redshift bins, and to recover controlled internal discrepancies in simulated KiDS data. We propose this as a diagnostic of internal consistency for present and future cosmological surveys, and as a tension metric for data sets that have non-negligible correlation, such as LSST and Euclid.

We investigate the impact of prior models on the upper bound of the sum of neutrino masses, $\sum m_{\nu}$. We use data from Large Scale Structure of galaxies, Cosmic Microwave Background, Type Ia SuperNovae, and Big Bang Nucleosynthesis. We probe physically motivated neutrino mass models (respecting oscillation experiment constraints) and compare them to constraints using standard cosmological approximations. The former give a consistent upper bound of $\sum m_{\nu} \lesssim 0.26$ eV ($95\%$ CI) and yields a strong competitive upper bound for the lightest neutrino mass species, $m_0^{\nu} < 0.086$ eV ($95\%$ CI). By contrast one of the approximations, which is somewhat inconsistent with oscillation experiments, yields an upper bound of $\sum m_{\nu} \lesssim 0.15$ eV ($95\%$ CI), which differs substantially from the former upper bound. We, therefore, argue that cosmological neutrino mass and hierarchy determination should be pursued using physically motivated models since approximations might lead to incorrect and nonphysical upper bounds.

We constrain cosmological parameters by analysing the angular power spectra of the Baryon Oscillation Spectroscopic Survey DR12 galaxies, a spectroscopic follow-up of around 1.3 million SDSS galaxies over 9,376 deg$^2$ with an effective volume of $\sim 6.5$ (Gpc $h^{-1}$)$^3$ in the redshift range $0.15 \leq z < 0.80$. We split this sample into 13 tomographic bins ($\Delta z = 0.05$); angular power spectra were calculated using a Pseudo-$C_{\ell}$ estimator, and covariance matrices were estimated using log-normal simulated maps. Cosmological constraints obtained from these data were combined with constraints from Planck CMB experiment as well as the JLA supernovae compilation. Considering a $w$CDM cosmological model measured on scales up to $k_{max} = 0.07h$ Mpc$^{-1}$, we constrain a constant dark energy equation-of-state with a $\sim 4\%$ error at the 1-$\sigma$ level: $w_0 = -0.993^{+0.046}_{-0.043}$, together with $\Omega_m = 0.330\pm 0.012$, $\Omega_b = 0.0505 \pm 0.002$, $S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3} = 0.863 \pm 0.016$, and $h = 0.661 \pm 0.012$. For the same combination of datasets, but now considering a $\Lambda$CDM model with massive neutrinos and the same scale cut, we find: $\Omega_m = 0.328 \pm 0.009$, $\Omega_b = 0.05017^{+0.0009}_{-0.0008}$, $S_8 = 0.862 \pm 0.017$, and $h = 0.663^{+0.006}_{-0.007}$ and a 95\% credible interval (CI) upper limit of $\sum m_{\nu} < 0.14$ eV for a normal hierarchy. These results are competitive if not better than standard analyses with the same dataset, and demonstrate this should be a method of choice for future surveys, opening the door for their full exploitation in cross-correlations probes.

Mapping the underlying density field, including non-visible 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 methods: Kaiser-Squires (KS), Wiener filter, and GLIMPSE. KS is a direct inversion, not accounting for survey masks or noise. The Wiener filter is well-motivated for Gaussian density fields in a Bayesian framework. GLIMPSE uses sparsity, aiming to reconstruct non-linearities in the density field. We compare these methods with several tests using public Dark Energy Survey (DES) Science Verification (SV) data and realistic DES simulations. The Wiener filter and GLIMPSE offer substantial improvements over smoothed KS with a range of metrics. Both the Wiener filter and GLIMPSE convergence reconstructions show a 12 per cent 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 (a standard data vector when inferring cosmology from peak statistics); the maximum signal-to-noise of these peak statistics is increased by a factor of 3.5 for the Wiener filter and 9 for GLIMPSE. With simulations we measure the reconstruction of the harmonic phases; the phase residuals' concentration is improved 17 per cent by GLIMPSE and 18 per cent by the Wiener filter. The correlation between reconstructions from data and foreground redMaPPer clusters is increased 18 per cent by the Wiener filter and 32 per cent by GLIMPSE.

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

Small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe 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 super-structure 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 top-hat 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 super-structures 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 super-structure data.

We present Magellan/IMACS spectroscopy of the recently discovered Milky Way satellite Tucana III (Tuc III). We identify 26 member stars in Tuc III, from which we measure a mean radial velocity of v_hel = -102.3 +/- 0.4 (stat.) +/- 2.0 (sys.) km/s, a velocity dispersion of 0.1^+0.7_-0.1 km/s, and a mean metallicity of [Fe/H] = -2.42^+0.07_-0.08. The upper limit on the velocity dispersion is sigma < 1.5 km/s at 95.5% confidence, and the corresponding upper limit on the mass within the half-light radius of Tuc III is 9.0 x 10^4 Msun. We cannot rule out mass-to-light ratios as large as 240 Msun/Lsun for Tuc III, but much lower mass-to-light ratios that would leave the system baryon-dominated are also allowed. We measure an upper limit on the metallicity spread of the stars in Tuc III of 0.19 dex at 95.5% confidence. Tuc III has a smaller metallicity dispersion and likely a smaller velocity dispersion than any known dwarf galaxy, but a larger size and lower surface brightness than any known globular cluster. Its metallicity is also much lower than those of the clusters with similar luminosity. We therefore tentatively suggest that Tuc III is the tidally-stripped remnant of a dark matter-dominated dwarf galaxy, but additional precise velocity and metallicity measurements will be necessary for a definitive classification. If Tuc III is indeed a dwarf galaxy, it is one of the closest external galaxies to the Sun. Because of its proximity, the most luminous stars in Tuc III are quite bright, including one star at V=15.7 that is the brightest known member star of an ultra-faint satellite.

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).