We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

We explore the relation between diffuse intracluster light (central galaxy included) and the galaxy cluster (baryonic and dark) matter distribution using a sample of 528 clusters at $0.2\leq z \leq 0.35$ found in the Dark Energy Survey (DES) Year 1 data. The surface brightness of the diffuse light shows an increasing dependence on cluster total mass at larger radius, and appears to be self-similar with a universal radial dependence after scaling by cluster radius.We also compare the diffuse light radial profiles to the cluster (baryonic and dark) matter distribution measured through weak lensing and find them to be comparable. The IllustrisTNG galaxy formation simulation offers further insight into the connection between diffuse stellar mass and cluster matter distributions -- the simulation radial profile of the diffuse stellar component does not have a similar slope with the total cluster matter content, although that of the cluster satellite galaxies does. Regardless of the radial trends, the amount of diffuse stellar mass has a low-scatter scaling relation with cluster's total mass in the simulation, out-performing the total stellar mass of cluster satellite galaxies. We conclude that there is no consistent evidence on whether or not diffuse light is a faithful radial tracer of the cluster matter distribution. Nevertheless, both observational and simulation results reveal that diffuse light is an excellent indicator of the cluster's total mass.

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultra-faint dwarf galaxy Grus~II. All stars exhibit a higher than expected $\mathrm{[Mg/Ca]}$ ratio compared to metal-poor stars in other ultra-faint dwarf galaxies and in the Milky Way halo. Nucleosynthesis in high mass ($\geqslant 20$M$_\odot$) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (3) stars suggest the chemical enrichment of Grus~II could have occurred through substantial high-mass stellar evolution and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it can not be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture ($r$-process) elements. The abundance pattern of the $r$-process elements in this star matches the scaled $r$-process pattern of the solar system and $r$-process enhanced stars in other dwarf galaxies and in the Milky Way halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of $r$-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus~II would increase the likelihood of any of these events occurring. The time delay between the $\alpha$ and $r$-process element enrichment of the galaxy favors a neutron star merger as the origin of the $r$-process elements in Grus~II.

We present deep $g$- and $r$-band Magellan/Megacam photometry of two dwarf galaxy candidates discovered in the Dark Energy Survey (DES), Grus I and Indus II (DES J2038-4609). For the case of Grus I, we resolved the main sequence turn-off (MSTO) and $\sim 2$ mags below it. The MSTO can be seen at $g_0\sim 24$ with a photometric uncertainty of $0.03$ mag. We show Grus I to be consistent with an old, metal-poor ($\sim 13.3$ Gyr, [Fe/H]$\sim-1.9$) dwarf galaxy. We derive updated distance and structural parameters for Grus I using this deep, uniform, wide-field data set. We find an azimuthally averaged half-light radius more than two times larger ($\sim 151^{+21}_{-31}$ pc; $\sim 4.^{\prime} 16^{+0.54}_{-0.74}$) and an absolute $V$-band magnitude $\sim-4.1$ that is $\sim 1$ magnitude brighter than previous studies. We obtain updated distance, ellipticity, and centroid parameters which are in agreement with other studies within uncertainties. Although our photometry of Indus II is $\sim 2-3$ magnitudes deeper than the DES Y1 Public release, we find no coherent stellar population at its reported location. The original detection was located in an incomplete region of sky in the DES Y2Q1 data set and was flagged due to potential blue horizontal branch member stars. The best fit isochrone parameters are physically inconsistent with both dwarf galaxies and globular clusters. We conclude that Indus II is likely a false-positive, flagged due to a chance alignment of stars along the line of sight.

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

We combine high-resolution spectroscopic data from APOGEE-2 Survey Data Release 16 (DR16) with broad-band photometric data from several sources, as well as parallaxes from {\it Gaia} Data Release 2 (DR2). Using the Bayesian isochrone-fitting code {\tt StarHorse}, we derive distances, extinctions and astrophysical parameters for around 388,815 APOGEE stars, achieving typical distance uncertainties of $\sim 6\%$ for APOGEE giants, $\sim 2\%$ for APOGEE dwarfs, as well as extinction uncertainties of $\sim 0.07$ mag when all photometric information is available, and $\sim 0.17$ mag if optical photometry is missing. {\tt StarHorse} uncertainties vary with the input spectroscopic catalogue, with the available photometry, and with the parallax uncertainties. To illustrate the impact of our results, we show that, thanks to {\it Gaia} DR2 and the now larger sky coverage of APOGEE-2 (including APOGEE-South), we obtain an extended map of the Galactic plane, providing unprecedented coverage of the disk close to the Galactic mid-plane ($|Z_{Gal}|<1$ kpc) from the Galactic Centre out to $R_{\rm Gal}\sim 20$ kpc. The improvements in statistics as well as distance and extinction uncertainties unveil the presence of the bar in stellar density, as well as the striking chemical duality in the innermost regions of the disk, now clearly extending to the inner bulge. We complement this paper with distances and extinctions for stars in other public released spectroscopic surveys: 324,999 in GALAH DR2, 4,928,715 in LAMOST DR5, 408,894 in RAVE DR6, and 6,095 in GES DR3

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{\Delta t}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $\Lambda$CDM cosmology and a uniform prior for $\Omega_{\rm m}$ as $\Omega_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$\sigma$ discrepancy from the cosmic microwave background measurement.

We present MWFitting, a method to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) from TRILEGAL models to real data. We apply MWFitting to photometric data from the first three years of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud and the Sagittarius Stream, our main sample spans a total area of $\sim$2,300 deg$^2$ distributed across the DES footprint. We further explore a smaller subset ($\sim$ 1,300 deg$^2$) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3\%. Based on the best-fit models, we create simulated stellar catalogues covering the whole DES footprint down to $g = 24$ magnitude. Comparisons of data and simulations provide evidence for a break in the power law index describing the stellar density of the Milky Way (MW) halo. Several previously discovered stellar over-densities are recovered in the residual stellar density map, showing the reliability of MWFitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict MW star counts for surveys such as LSST and Euclid.

We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843~MHz radio AGN population in galaxy clusters from two large catalogs created using the Dark Energy Survey (DES): $\sim$11,800 optically selected RM-Y3 and $\sim$1,000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGN are highly concentrated around cluster centers to $z\sim1$. We measure the halo occupation number for cluster radio AGN above a threshold luminosity, finding that the number of radio AGN per cluster increases with cluster halo mass as $N\propto M^{1.2\pm0.1}$ ($N\propto M^{0.68\pm0.34}$) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad band colors and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGN luminosity function (LF). The LF exhibits a redshift trend of the form $(1+z)^\gamma$ in density and luminosity, respectively, of $\gamma_{\rm D}=3.0\pm0.4$ and $\gamma_{\rm P}=0.21\pm0.15$ in the RM-Y3 sample, and $\gamma_{\rm D}=2.6\pm0.7$ and $\gamma_{\rm P}=0.31\pm0.15$ in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGN, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core ($<0.1R_{500}$) X-ray radiative losses for clusters at $z<1$.

We present Magellan/IMACS spectroscopy of three recently discovered ultra-faint Milky Way satellites, Grus II, Tucana IV, and Tucana V. We measure systemic velocities of V_hel = -110.0 +/- 0.5 km/s, V_hel = 15.9 +/- 1.8 km/s, and V_hel = -36.2 +/-2.5 km/s for the three objects, respectively. Their large relative velocity differences demonstrate that the satellites are unrelated despite their close physical proximity. We determine a velocity dispersion for Tuc IV of sigma = 4.3^+1.7_-1.0 km/s, but we cannot resolve the velocity dispersions of the other two systems. For Gru II we place an upper limit (90% confidence) on the dispersion of sigma < 1.9 km/s, and for Tuc V we do not obtain any useful limits. All three satellites have metallicities below [Fe/H] = -2.1, but none has a detectable metallicity spread. We determine proper motions for each satellite based on Gaia astrometry and compute their orbits around the Milky Way. Gru II is on a tightly bound orbit with a pericenter of 25 kpc and orbital eccentricity of 0.45. Tuc V likely has an apocenter beyond 100 kpc, and could be approaching the Milky Way for the first time. The current orbit of Tuc IV is similar to that of Gru II, with a pericenter of 25 kpc and an eccentricity of 0.36. However, a backward integration of the position of Tuc IV demonstrates that it collided with the Large Magellanic Cloud at an impact parameter of 4 kpc ~120 Myr ago, deflecting its trajectory and possibly altering its internal kinematics. Based on their sizes, masses, and metallicities, we classify Gru II and Tuc IV as likely dwarf galaxies, but the nature of Tuc V remains uncertain.

We present the discovery of 255 binary and six multiple system candidates with wide (> 5") separation composed by L or T dwarfs companions to stars, plus nine double brown dwarf systems. The sample of brown dwarf candidates was found in the Dark Energy Survey and the possible stellar companions are from Gaia DR2 and DES data. Our search is based in a common distance criterion with no proper motion information. For the Gaia DR2 stars we estimate distances based on their parallaxes and photometry, using the StarHorse code, while for DES stars, the StarHorse distances were purely photometric for the majority of cases, with a fraction having parallax measurement from Gaia DR2. L and T dwarfs distances are based on empirical templates ranging from L0 to T9. We also compute chance alignment probabilities in order to assess the physical nature of each pair. We find 174 possible pairs with Gaia DR2 primaries with chance alignment probabilities < 5%. We also find 85 binary pair candidates with a DES star as a primary, 81 of them with chance alignment probabilities < 5%. Only nine candidate systems composed of two brown dwarfs were identified. The sample of multiple systems is made up of five triple systems and one quadruple system. We determine that the typical wide binary fraction over the L and T spectral types is 2-4%. The significant leap provided by this sample will enable constraints on the formation and evolution of L and T dwarfs.

This work presents the first search for RR Lyrae stars (RRLs) in four of the ultra-faint systems imaged by the Dark Energy Survey (DES) using SOAR/Goodman and Blanco/DECam imagers. We have detected two RRLs in the field of Grus I, none in Kim 2, one in Phoenix II, and four in Grus II. With the detection of these stars, we accurately determine the distance moduli for these ultra-faint dwarf satellite galaxies; $\mu_0$=20.51$\pm$0.10 mag (D$_{\odot}$=127$\pm$6 kpc) for Grus I and $\mu_0$=20.01$\pm$0.10 mag (D$_{\odot}$=100$\pm$5 kpc) for Phoenix II. These measurements are larger than previous estimations by Koposov et al. 2015 and Bechtol et al. 2015, implying larger physical sizes; 5\% for Grus I and 33\% for Phoenix II. For Grus II, out of the four RRLs detected, one is consistent with being a member of the galactic halo (D$_\odot$=24$\pm$1 kpc, $\mu_0$=16.86$\pm$0.10 mag), another is at D$_\odot$=55$\pm$2 kpc ($\mu_0$=18.71$\pm$0.10 mag), which we associate with Grus II, and the two remaining at D$_\odot$=43$\pm$2 kpc ($\mu_0$=18.17$\pm$0.10 mag). Moreover, the appearance of a subtle red horizontal branch in the color-magnitude diagram of Grus II at the same brightness level of the latter two RRLs, which are at the same distance and in the same region, suggests that a more metal-rich system may be located in front of Grus II. The most plausible scenario is the association of these stars with the Chenab/Orphan Stream. Finally, we performed a comprehensive and updated analysis of the number of RRLs in dwarf galaxies. This allows us to predict that the method of finding new ultra-faint dwarf galaxies by using two or more clumped RRLs will work only for systems brighter than M$_V\sim-6$ mag.

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr.

We present a chemical abundance analysis of four additional confirmed member stars of Tucana III, a Milky Way satellite galaxy candidate in the process of being tidally disrupted as it is accreted by the Galaxy. Two of these stars are centrally located in the core of the galaxy while the other two stars are located in the eastern and western tidal tails. The four stars have chemical abundance patterns consistent with the one previously studied star in Tucana III: they are moderately enhanced in $r$-process elements, i.e. they have $<$[Eu/Fe]$> \approx +$0.4 dex. The non-neutron-capture elements generally follow trends seen in other dwarf galaxies, including a metallicity range of 0.44 dex and the expected trend in $\alpha$-elements, i.e., the lower metallicity stars have higher Ca and Ti abundance. Overall, the chemical abundance patterns of these stars suggest that Tucana III was an ultra-faint dwarf galaxy, and not a globular cluster, before being tidally disturbed. As is the case for the one other galaxy dominated by $r$-process enhanced stars, Reticulum II, Tucana III's stellar chemical abundances are consistent with pollution from ejecta produced by a binary neutron star merger, although a different $r$-process element or dilution gas mass is required to explain the abundances in these two galaxies if a neutron star merger is the sole source of $r$-process enhancement.

Combining the precise parallaxes and optical photometry delivered by Gaia's second data release (Gaia DR2) with the photometric catalogues of PanSTARRS-1, 2MASS, and AllWISE, we derive Bayesian stellar parameters, distances, and extinctions for 265 million stars brighter than G=18. Because of the wide wavelength range used, our results substantially improve the accuracy and precision of previous extinction and effective temperature estimates. After cleaning our results for both unreliable input and output data, we retain 137 million stars, for which we achieve a median precision of 5% in distance, 0.20 mag in V-band extinction, and 245 K in effective temperature for G<14, degrading towards fainter magnitudes (12%, 0.20 mag, and 245 K at G=16; 16%, 0.23 mag, and 260 K at G=17, respectively). We find a very good agreement with the asteroseismic surface gravities and distances of 7000 stars in the Kepler, the K2-C3, and the K2-C6 fields, with stellar parameters from the APOGEE survey, as well as with distances to star clusters. Our results are available through the ADQL query interface of the Gaia mirror at the Leibniz-Institut f\"{u}r Astrophysik Potsdam (gaia.aip.de) and as binary tables at data.aip.de. As a first application, in this paper we provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps as a function of distance, and extensive density maps, demonstrating the potential of our value-added dataset for mapping the three-dimensional structure of our Galaxy. In particular, we see a clear manifestation of the Galactic bar in the stellar density distributions, an observation that can almost be considered a direct imaging of the Galactic bar.

We present the analysis underpinning the measurement of cosmological parameters from 207 spectroscopically classified type Ia supernovae (SNe Ia) from the first three years of the Dark Energy Survey Supernova Program (DES-SN), spanning a redshift range of 0.017<$z$<0.849. We combine the DES-SN sample with an external sample of 122 low-redshift ($z$<0.1) SNe Ia, resulting in a "DES-SN3YR" sample of 329 SNe Ia. Our cosmological analyses are blinded: after combining our DES-SN3YR distances with constraints from the Cosmic Microwave Background (CMB; Planck Collaboration 2016), our uncertainties in the measurement of the dark energy equation-of-state parameter, $w$, are .042 (stat) and .059 (stat+syst) at 68% confidence. We provide a detailed systematic uncertainty budget, which has nearly equal contributions from photometric calibration, astrophysical bias corrections, and instrumental bias corrections. We also include several new sources of systematic uncertainty. While our sample is <1/3 the size of the Pantheon sample, our constraints on $w$ are only larger by 1.4$\times$, showing the impact of the DES SN Ia light curve quality. We find that the traditional stretch and color standardization parameters of the DES SNe Ia are in agreement with earlier SN Ia samples such as Pan-STARRS1 and the Supernova Legacy Survey. However, we find smaller intrinsic scatter about the Hubble diagram (0.077 mag). Interestingly, we find no evidence for a Hubble residual step ( 0.007 $\pm$ 0.018 mag) as a function of host galaxy mass for the DES subset, in 2.4$\sigma$ tension with previous measurements. We also present novel validation methods of our sample using simulated SNe Ia inserted in DECam images and using large catalog-level simulations to test for biases in our analysis pipelines.

We report an analysis of the metal-rich tail ([Fe/H] $> -0.75$) of halo stars located at distances from the Galactic plane $z$ up to $|z| \sim 10$ kpc, observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We examine the chemistry, kinematics, and dynamics of this metal-rich halo sample using chemical abundances and radial velocities provided by the fourteenth APOGEE data release (DR14) and proper motions from the second Gaia data release (DR2). The analysis reveals three chemically different stellar populations in the [Mg/Fe] vs. [Fe/H] space -- the two distinct halo populations already reported in the literature, and a third group with intermediate [Mg/Fe] $\sim$+0.1. We derive the $U$, $V$ and $W$ velocity components with respect to the Local Standard of Rest, as well as orbits for the three stellar groups, and find that they differ also in their kinematical and dynamical properties. The high-[Mg/Fe] population exhibits a mean prograde rotation, as well as orbits that are more bound and closer to the plane, whereas the low-[Mg/Fe] population has $<V>$ closer to 0, and stars that move in less-bound orbits reaching larger distances from the centre and the Galactic plane. The intermediate-Mg stars exhibit different orbital characteristics, moving with a strong prograde rotation and low excentricity, but in less-bound orbits. This stellar population resembles the two stellar overdensities lying about $|z| \sim 5$ kpc recently reported in the literature, for which a disc origin has been claimed.

Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.

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/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread 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.

In this paper we present a catalogue of 11,745 brown dwarfs with spectral types ranging from L0 to T9, photometrically classified using data from the Dark Energy Survey (DES) year 3 release matched to the Vista Hemisphere Survey (VHS) DR3 and Wide-field Infrared Survey Explorer (WISE) data, covering approx 2,400 deg2 up to i_AB=22. The classification method follows the same photo-type method previously applied to SDSS-UKIDSS-WISE data. The most significant difference comes from the use of DES data instead of SDSS, which allow us to classify almost an order of magnitude more brown dwarfs than any previous search and reaching distances beyond 400 parsecs for the earliest types. Next, we also present and validate the GalmodBD simulation, which produces brown dwarf number counts as a function of structural parameters with realistic photometric properties of a given survey. We use this simulation to estimate the completeness and purity of our photometric LT catalogue down to i_AB=22, as well as to compare to the observed number of LT types. We put constraints on the thin disk scale height for the early L population to be around 450 parsecs, in agreement with previous findings. For completeness, we also publish in a separate table a catalogue of 20,863 M dwarfs that passed our colour cut with spectral types greater than M6. Both the LT and the late M catalogues are found at https://des.ncsa.illinois.edu/releases/other/y3-mlt

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 Anglo-Australian 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 metal-poor 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 metal-rich disk stars and metal-poor stream members. This metallicity-dependent color offers a more efficient method to recognize metal-poor targets and will increase the selection efficiency of stream members for future spectroscopic follow-up programs on stellar streams.

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 two-point 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 galaxy-shear cross-correlation 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 two-point 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 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

Since first noticed by Shapley in 1939, a faint object coincident with the Fornax dwarf spheroidal has long been discussed as a possible sixth globular cluster system. However, debate has continued over whether this overdensity is a statistical artifact or a blended galaxy group. In this Letter we demonstrate, using deep DECam imaging data, that this object is well resolved into stars and is a bona fide star cluster. The stellar overdensity of this cluster is statistically significant at the level of ~ 6 - 6.7 sigma in several different photometric catalogs including Gaia. Therefore, it is highly unlikely to be caused by random fluctuation. We show that Fornax 6 is a star cluster with a peculiarly low surface brightness and irregular shape, which may indicate a strong tidal influence from its host galaxy. The Hess diagram of Fornax 6 is largely consistent with that of Fornax field stars, but it appears to be slightly bluer. However, it is still likely more metal-rich than most of the globular clusters in the system. Faint clusters like Fornax 6 that orbit and potentially get disrupted in the centers of dwarf galaxies can prove crucial for constraining the dark matter distribution in Milky Way satellites.

We present angular diameter distance measurements obtained by locating the BAO scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1318 deg$^2$ with $0.6 < z_{\rm photo} < 1$ and a typical redshift uncertainty of $0.03(1+z)$. This sample was selected, as fully described in a companion paper, using a color/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the co-moving transverse separation, and spherical harmonics. Further, we compare results obtained from template based and machine learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, $D_A$, at the effective redshift of our sample divided by the true physical scale of the BAO feature, $r_{\rm d}$. We obtain close to a 4 per cent distance measurement of $D_A(z_{\rm eff}=0.81)/r_{\rm d} = 10.75\pm 0.43 $. These results are consistent with the flat $\Lambda$CDM concordance cosmological model supported by numerous other recent experimental results. All data products are publicly available here: https://des.ncsa.illinois.edu/releases/y1a1/bao

We perform a comparison of different approaches to star-galaxy classification using the broad-band 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: large-scale 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 mis-classification, contamination can be reduced to the O(1%) level by using multi-epoch 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 present a structural and morphological catalogue for 45 million objects selected from the first year of data from the Dark Energy Survey (DES). Single Sersic fits and non-parametric measurements are produced for g, r and i filters. The parameters from the best-fitting Sersic model (total magnitude, half-light radius, Sersic index, axis ratio and position angle) are measured with Galfit; the non-parametric coefficients (concentration, asymmetry, clumpiness, Gini, M20) are provided using the Zurich Estimator of Structural Types (ZEST+). To study the statistical uncertainties, we consider a sample of state-of-the-art image simulations with a realistic distribution in the input parameter space and then process and analyse them as we do with real data: this enables us to quantify the observational biases due to PSF blurring and magnitude effects and correct the measurements as a function of magnitude, galaxy size, Sersic index (concentration for the analysis of the non-parametric measurements) and ellipticity. We present the largest structural catalogue to date: we find that accurate and complete measurements for all the structural parameters are typically obtained for galaxies with SExtractor MAG AUTO I < 21. Indeed, the parameters in the filters i and r can be overall well recovered up to MAG AUTO < 21.5, corresponding to a fitting completeness of ~90% below this threshold, for a total of 25 million galaxies. The combination of parametric and non-parametric structural measurements makes this catalogue an important instrument to explore and understand how galaxies form and evolve. The catalogue described in this paper will be publicly released alongside the Dark Energy Survey collaboration Y1 cosmology data products at the following URL: https://des.ncsa.illinois.edu/releases/y1a1/gold/morphology.

Using deep wide-field photometry three-year data (Y3) from the Dark Energy Survey (DES), we present a panoramic study of the Fornax dwarf spheroidal galaxy. The data presented here -- a small subset of the full survey -- uniformly covers a region of 25 square degrees centered on the galaxy to a depth of g ~ 23.5. We use this data to study the structural properties of Fornax, overall stellar population, and its member stars in different evolutionary phases. We also search for possible signs of tidal disturbance. Fornax is found to be significantly more spatially extended than what early studies suggested. No statistically significant distortions or signs of tidal disturbances were found down to a surface brightness limit of ~32.1 mag/arcsec^2. However, there are hints of shell-like features located ~30' - 40' from the center of Fornax that may be stellar debris from past merger events. We also find that intermediate age and young main-sequence populations show different orientation at the galaxy center and have many substructures. The deep DES Y3 data allows us to characterize the age of those substructures with great accuracy, both those previously known and those newly discovered in this work, on the basis of their color-magnitude diagram morphology. We find that the youngest overdensities are all found on the Eastern side of Fornax, where the Fornax field population itself is slightly younger than in the West. The high quality DES Y3 data reveals that Fornax has many rich structures, and provides insights into its complex formation history.

A significant challenge facing photometric surveys for cosmological purposes is the need to produce reliable redshift estimates. The estimation of photometric redshifts (photo-zs) has been consolidated as the standard strategy to bypass the high production costs and incompleteness of spectroscopic redshift samples. Training-based photo-z methods require the preparation of a high-quality list of spectroscopic redshifts, which needs to be constantly updated. The photo-z training, validation, and estimation must be performed in a consistent and reproducible way in order to accomplish the scientific requirements. To meet this purpose, we developed an integrated web-based data interface that not only provides the framework to carry out the above steps in a systematic way, enabling the ease testing and comparison of different algorithms, but also addresses the processing requirements by parallelizing the calculation in a transparent way for the user. This framework called the Science Portal (hereafter Portal) was developed in the context the Dark Energy Survey (DES) to facilitate scientific analysis. In this paper, we show how the Portal can provide a reliable environment to access vast data sets, provide validation algorithms and metrics, even in the case of multiple photo-zs methods. It is possible to maintain the provenance between the steps of a chain of workflows while ensuring reproducibility of the results. We illustrate how the Portal can be used to provide photo-z estimates using the DES first year (Y1A1) data. While the DES collaboration is still developing techniques to obtain more precise photo-zs, having a structured framework like the one presented here is critical for the systematic vetting of DES algorithmic improvements and the consistent production of photo-zs in the future DES releases.

In the era of industrial Galactic astronomy and multi-object spectroscopic stellar surveys, the sample sizes and the number of available stellar chemical abundances have reached dimensions in which it has become difficult to process all the available information in an effective manner. In this paper we demonstrate the use of a dimensionality-reduction technique (t-distributed stochastic neighbour embedding; t-SNE) for analysing the stellar abundance-space distribution. While the non-parametric non-linear behaviour of this technique makes it difficult to estimate the significance of found abundance-space substructure, we show that our results depend little on parameter choices and are robust to abundance errors. By reanalysing the high-resolution high-signal-to-noise solar-neighbourhood HARPS-GTO sample with t-SNE, we find clearer chemical separations of the high- and low-[$\alpha$/Fe] disc sequences, hints for multiple populations in the high-[$\alpha$/Fe] population, and indications that the chemical evolution of the high-[$\alpha$/Fe] metal-rich stars is connected with the super-metal-rich stars. We also identify a number of chemically peculiar stars, among them a high-confidence s-process-enhanced abundance-ratio pair (HD91345/HD126681) with very similar ages and $v_X$ and $v_Y$ velocities, which we suggest to have a common birth origin, possibly a dwarf galaxy. Our results demonstrate the potential of abundance-space t-SNE and similar methods for chemical-tagging studies with large spectroscopic surveys.

We present a semi-empirical, largely model-independent approach for estimating Galactic birth radii, r_birth, for Milky Way disk stars. The technique relies on the justifiable assumption that a negative radial metallicity gradient in the interstellar medium (ISM) existed for most of the disk lifetime. Stars are projected back to their birth positions according to the observationally derived age and [Fe/H] with no kinematical information required. Applying our approach to the AMBRE:HARPS and HARPS-GTO local samples, we show that we can constrain the ISM metallicity evolution with Galactic radius and cosmic time, [Fe/H]_ISM(r, t), by requiring a physically meaningful r_birth distribution. We find that the data are consistent with an ISM radial metallicity gradient that flattens with time from ~-0.15 dex/kpc at the beginning of disk formation, to its measured present-day value (-0.07 dex/kpc). We present several chemo-kinematical relations in terms of mono-r_birth populations. One remarkable result is that the kinematically hottest stars would have been born locally or in the outer disk, consistent with thick disk formation from the nested flares of mono-age populations and predictions from cosmological simulations. This phenomenon can be also seen in the observed age-velocity dispersion relation, in that its upper boundary is dominated by stars born at larger radii. We also find that the flatness of the local age-metallicity relation (AMR) is the result of the superposition of the AMRs of mono-r_birth populations, each with a well-defined negative slope. The solar birth radius is estimated to be 7.3+-0.6 kpc, for a current Galactocentric radius of 8 kpc.

We present a novel approach for creating science-ready catalogs through a software infrastructure developed for the Dark Energy Survey (DES). We integrate the data products released by the DES Data Management and additional products created by the DES collaboration in an environment known as DES Science Portal. Each step involved in the creation of a science-ready catalog is recorded in a relational database and can be recovered at any time. We describe how the DES Science Portal automates the creation and characterization of lightweight catalogs for DES Year 1 Annual Release, and show its flexibility in creating multiple catalogs with different inputs and configurations. Finally, we discuss the advantages of this infrastructure for large surveys such as DES and the Large Synoptic Survey Telescope. The capability of creating science-ready catalogs efficiently and with full control of the inputs and configurations used is an important asset for supporting science analysis using data from large astronomical surveys.

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 lensing-weighted 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 30-band photometric redshifts of galaxies in the COSMOS field to BPZ estimates made from the DES griz fluxes, for a sample matched in fluxes, pre-seeing size, and lensing weight to the DES weak-lensing 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 photo-z 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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

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.

Combining stellar atmospheric parameters, such as effective temperature, surface gravity, and metallicity, with barycentric radial velocity data provides insight into the chemo-dynamics of the MilkyWay and our local Galactic environment. We analyze 3075 stars with spectroscopic data from the Sloan Digital Sky Survey III (SDSS-III) MARVELS radial velocity survey and present atmospheric parameters for 2343 dwarf stars using the spectral indices method, a modified version of the equivalent width method. We present barycentric radial velocities for 2610 stars with uncertainties of $\sim$0.3 km s$^{-1}$. We determine stellar ages using two independent methods and calculate ages for 2343 stars with a maximum-likelihood isochronal age-dating method and for 1694 stars with a Bayesian age-dating method. Using previously published parallax data we compute Galactic orbits and space velocities for 1995 stars to explore stellar populations based on kinematic and age parameters. This study combines good ages and exquisite velocities to explore local chemo-kinematics of the Milky Way, which complements many of the recent studies of giant stars with the APOGEE survey, and we find our results to be in agreement with current chemo-dynamical models of the Milky Way. Particularly, we find from our metallicity distributions and velocity-age relations of a kinematically-defined thin disk that the metal rich end has stars of all ages, even after we clean the sample of highly eccentric stars, suggesting that radial migration plays a key role in the metallicity scatter of the thin disk. All stellar parameters and kinematic data derived in this work are catalogued and published online in machine-readable form.

Understanding the formation and evolution of our Galaxy requires accurate distances, ages and chemistry for large populations of field stars. Here we present several updates to our spectro-photometric distance code, that can now also be used to estimate ages, masses, and extinctions for individual stars. Given a set of measured spectro-photometric parameters, we calculate the posterior probability distribution over a given grid of stellar evolutionary models, using flexible Galactic stellar-population priors. The code (called {\tt StarHorse}) can acommodate different observational datasets, prior options, partially missing data, and the inclusion of parallax information into the estimated probabilities. We validate the code using a variety of simulated stars as well as real stars with parameters determined from asteroseismology, eclipsing binaries, and isochrone fits to star clusters. Our main goal in this validation process is to test the applicability of the code to field stars with known {\it Gaia}-like parallaxes. The typical internal precision (obtained from realistic simulations of an APOGEE+Gaia-like sample) are $\simeq 8\%$ in distance, $\simeq 20\%$ in age,$\simeq 6\ %$ in mass, and $\simeq 0.04$ mag in $A_V$. The median external precision (derived from comparisons with earlier work for real stars) varies with the sample used, but lies in the range of $\simeq [0,2]\%$ for distances, $\simeq [12,31]\%$ for ages, $\simeq [4,12]\%$ for masses, and $\simeq 0.07$ mag for $A_V$. We provide StarHorse distances and extinctions for the APOGEE DR14, RAVE DR5, GES DR3 and GALAH DR1 catalogues.

We describe the DR14 APOGEE-TGAS catalogue, a new SDSS value-added catalogue that provides precise astrophysical parameters, chemical abundances, astro-spectro-photometric distances and extinctions, as well as orbital parameters for $\sim 30,000$ APOGEE-TGAS stars, among them $\sim5,000$ high-quality giant stars within 1 kpc.

We report the characterization of the first $62$ MaNGA Active Galactic Nuclei (AGN) hosts in the Fifth Product Launch (MPL-5) and the definition of a control sample of non-active galaxies. This control sample - comprising two galaxies for each AGN - was selected in order to match the AGN hosts in terms of stellar mass, redshift, visual morphology and inclination. The stellar masses are in the range $9.4<\log(M/M_\odot)<11.5$, and most objects have redshifts $\leq 0.08$. The AGN sample is mostly comprised of low-luminosity AGN, with only 17 nuclei with $L([OIII]\lambda 5007)\geq 3.8\times 10^{40}$ erg s$^{-1}$ (that we call "strong AGN"). The stellar population of the control sample galaxies within the inner $1$-$3$ kpc is dominated by the old ($\sim$ $4$ - $13$ Gyr) age component, with a small contribution of intermediate age ($\sim 640$-$940$ Myr) and young stars ($\leq 40$ Myr) to the total light at $5700\AA$. While the weaker AGN show a similar age distribution to that of the control galaxies, the strong AGN show an increased contribution of younger stars and a decreased contribution of older stars. Examining the relationship between the AGN stellar population properties and $L([OIII])$, we find that with increasing $L([OIII])$, the AGN exhibit a decreasing contribution from the oldest ($>4$ Gyr) stellar population relative to control galaxies, but have an increasing contribution from the younger components with ages $\sim 40$ Myr. We also find a correlation of the mean age differences (AGN - control) with $L([OIII])$, in the sense that more luminous AGN are younger than the control objects, while the low-luminosity AGN seem to be older. These results support a connection between the growth of the galaxy bulge via formation of new stars and the growth of the Supermassive Black Hole via matter accretion in the AGN phase.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We report the discovery of a new star cluster, DES 3, in the constellation of Indus, and deeper observations of the previously identified satellite DES J0222.7$-$5217 (Eridanus III). DES 3 was detected as a stellar overdensity in first-year Dark Energy Survey data, and confirmed with deeper photometry from the 4.1 metre Southern Astrophysical Research (SOAR) telescope. The new system was detected with a relatively high significance and appears in the DES images as a compact concentration of faint blue point sources. We determine that DES 3 is located at a heliocentric distance of $\sim 76\,\mathrm{kpc}$ and it is dominated by an old ($\simeq 9.8\,\mathrm{Gyr}$) and metal-poor ($\mathrm{[Fe/H]}\simeq -1.88$) population. While the age and metallicity values of DES 3 are similar to globular clusters, its half-light radius ($r_\mathrm{h}\sim 6.5\,\mathrm{pc}$) and luminosity ($M_V \sim -1.9$) are more indicative of faint star clusters. Based on the apparent angular size, DES 3, with a value of $r_\mathrm{h}\sim 0.\!^{\prime}3$, is among the smallest faint star clusters known to date. Furthermore, using deeper imaging of DES J0222.7$-$5217 taken with the SOAR telescope, we update structural parameters and perform the first isochrone modeling. Our analysis yields the first age ($\simeq 12.6\,\mathrm{Gyr}$) and metallicity ($\mathrm{[Fe/H]}\simeq -2.01$) estimates for this object. The half-light radius ($r_\mathrm{h}\sim 10.5\,\mathrm{pc}$) and luminosity ($M_V\sim -2.7$) of DES J0222.7$-$5217 suggest that it is likely a faint star cluster. The discovery of DES 3 indicates that the census of stellar systems in the Milky Way is still far from complete, and demonstrates the power of modern wide-field imaging surveys to improve our knowledge of the Galaxy's satellite population.

The outer Galactic halo is home to a number of substructures which still have an uncertain origin, but most likely are remnants of former interactions between the Galaxy and its former satellites. Triangulum-Andromeda (TriAnd) is one of these halo substructures, found as an overdensity of 2MASS M giants. We analyzed the region of Triangulum-Andromeda using photometric data from the Ninth Data Release of Sloan Digital Sky Survey (SDSS DR9). By comparing the observations with simulations from the TRILEGAL Galactic model, we were able to identify and map several scattered overdensities of main sequence stars that seem to be associated with TriAnd over a large area covering $\sim 500$ deg$^2$. One of these excesses may represent a new stellar overdensity. We also briefly discuss an alternative hypothesis, according to which TriAnd is one of the troughs of oscillation rings in the Galactic disk.

We present chemical abundance measurements of three stars in the ultra-faint dwarf galaxy Horologium I, a Milky Way satellite discovered by the Dark Energy Survey. Using high resolution spectroscopic observations we measure the metallicity of the three stars as well as abundance ratios of several $\alpha$-elements, iron-peak elements, and neutron-capture elements. The abundance pattern is relatively consistent among all three stars, which have a low average metallicity of [Fe/H] $\sim -2.6$ and are not $\alpha$-enhanced ([$\alpha$/Fe] $\sim 0.0$). This result is unexpected when compared to other low-metallicity stars in the Galactic halo and other ultra-faint dwarfs and hints at an entirely different mechanism for the enrichment of Hor I compared to other satellites. We discuss possible scenarios that could lead to this observed nucleosynthetic signature including extended star formation, a Population III supernova, and a possible association with the Large Magellanic Cloud.

We report the discovery of a stellar over-density 8$^{\circ}$ north of the center of the Small Magellanic Cloud (Small Magellanic Cloud Northern Over-Density; SMCNOD) using data from the first two years of the Dark Energy Survey (DES) and the first year of the MAGellanic SatelLITEs Survey (MagLiteS). The SMCNOD is indistinguishable in age, metallicity and distance from the nearby SMC stars, being primarly composed of intermediate-age stars (6 Gyr, Z=0.001), with a small fraction of young stars (1 Gyr, Z=0.01). The SMCNOD has an elongated shape with an ellipticity of 0.6 and a size of $\sim$ 6x2 deg. It has an absolute magnitude of $M_V \cong$ -7.7, $r_h = 2.1$ kpc, and $\mu_V(r<r_h)$ = 31.2 mag arcsec$^{-2}$. We estimate a stellar mass of $\sim 10^5$ $M_{\odot}$, following a Kroupa mass function. The SMCNOD was probably removed from the SMC disk by tidal stripping, since it is located near the head of the Magellanic Stream, and the literature indicates likely recent LMC-SMC encounters. This scenario is supported by the lack of significant HI gas. Other potential scenarios for the SMCNOD origin are a transient over-density within the SMC tidal radius or a primordial SMC satellite in advanced stage of disruption.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We describe a multi-component matched filter cluster confirmation tool (MCMF) designed for the study of large X-ray source catalogs produced by the upcoming X-ray all-sky survey mission eROSITA. We apply the method to confirm a sample of 88 clusters with redshifts $0.05<z<0.8$ in the recently published 2RXS catalog from the ROSAT all-sky survey (RASS) over the 208 deg$^2$ region overlapped by the Dark Energy Survey (DES) science verification (DES-SV) dataset. In our pilot study, we examine all X-ray sources, regardless of their extent. Our method employs a multi-color red sequence (RS) algorithm that incorporates the X-ray count rate and peak position in determining the region of interest for followup and extracts the positionally and color weighted optical richness $\lambda_{\mathrm{MCMF}}$ as a function of redshift for each source. Peaks in the $\lambda_{\mathrm{MCMF}}$-redshift distribution are identified and used to extract photometric redshifts, richness and uncertainties. The significances of all optical counterparts are characterized using the distribution of richnesses defined along random lines of sight. These significances are used to extract cluster catalogs and to estimate the contamination by random superpositions of unassociated optical systems. The delivered photometric redshift accuracy is $\delta z / (1+z)=0.010$. We find a well defined X-ray luminosity-$\lambda_{\mathrm{MCMF}}$ relation with an intrinsic scatter of $\delta \ln(\lambda_\mathrm{MCMF}| L_\mathrm{x})=0.21$. Matching our catalog with the DES-SV redMaPPer catalog yields good agreement in redshift and richness estimates; comparing our catalog with the South Pole Telescope (SPT) selected clusters shows no inconsistencies. SPT clusters in our dataset are consistent with the high mass extension of the RASS based $\lambda_{\mathrm{MCMF}}$-mass relation

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.

Chemically peculiar stars in dwarf galaxies provide a window for exploring the birth environment of stars with varying chemical enrichment. We present a chemical abundance analysis of the brightest star in the newly discovered ultra-faint dwarf galaxy candidate Tucana III. Because it is particularly bright for a star in an ultra-faint Milky Way satellite, we are able to measure the abundance of 28 elements, including 13 neutron-capture species. This star, DES J235532.66$-$593114.9 (DES J235532), shows a mild enhancement in neutron-capture elements associated with the $r$-process and can be classified as an $r$-I star. DES J235532 is the first $r$-I star to be discovered in an ultra-faint satellite, and Tuc III is the second extremely low-luminosity system found to contain $r$-process enriched material, after Reticulum II. Comparison of the abundance pattern of DES J235532 with $r$-I and $r$-II stars found in other dwarf galaxies and in the Milky Way halo suggests a common astrophysical origin for the neutron-capture elements seen in all $r$-process enhanced stars. We explore both internal and external scenarios for the $r$-process enrichment of Tuc III and show that with abundance patterns for additional stars it should be possible to distinguish between them.

We present Magellan/IMACS spectroscopy of the recently-discovered Milky Way satellite Eridanus II (Eri II). We identify 28 member stars in Eri II, from which we measure a systemic radial velocity of $v_{\rm hel} = 75.6 \pm 1.3~\mbox{(stat.)} \pm 2.0~\mbox{(sys.)}~\mathrm{km\,s^{-1}}$ and a velocity dispersion of $6.9^{+1.2}_{-0.9}~\mathrm{km\,s^{-1}}$. Assuming that Eri~II is a dispersion-supported system in dynamical equilibrium, we derive a mass within the half-light radius of Eri II is $1.2^{+0.4}_{-0.3} \times 10^{7}~\mathrm{M_\odot}$, indicating a mass-to-light ratio of $420^{+210}_{-140}~\mathrm{M_\odot}/\mathrm{L_\odot}$ and confirming that it is a dark matter-dominated dwarf galaxy. From the equivalent width measurements of the CaT lines of 16 red giant member stars, we derive a mean metallicity of ${\rm [Fe/H]} = -2.38 \pm 0.13$ and a metallicity dispersion of $\sigma_{\rm [Fe/H]} = 0.47 ^{+0.12}_{-0.09}$. The velocity of Eri II in the Galactic Standard of Rest frame is $v_{\rm GSR} = -66.6~\mathrm{km\,s^{-1}}$, indicating that either Eri II is falling into the Milky Way potential for the first time or it has passed the apocenter of its orbit on a subsequent passage. At a Galactocentric distance of $\sim$370 kpc, Eri II is one of the Milky Way's most distant satellites known. Additionally, we show that the bright blue stars previously suggested to be a young stellar population are not associated with Eri II. The lack of gas and recent star formation in Eri II is surprising given its mass and distance from the Milky Way, and may place constraints on models of quenching in dwarf galaxies and on the distribution of hot gas in the Milky Way halo. Furthermore, the large velocity dispersion of Eri II can be combined with the existence of a central star cluster to constrain MACHO dark matter with mass $\gtrsim10~\mathrm{M_\odot}$.

Planet searches using the radial velocity technique show a paucity of companions to solar-type stars within ~5 AU in the mass range of ~10 - 80 M$_{\text{Jup}}$. This deficit, known as the brown dwarf desert, currently has no conclusive explanation. New substellar companions in this region help asses the reality of the desert and provide insight to the formation and evolution of these objects. Here we present 10 new brown dwarf and two low-mass stellar companion candidates around solar-type stars from the Multi-object APO Radial-Velocity Exoplanet Large-Area Survey (MARVELS) of the Sloan Digital Sky Survey III (SDSS-III). These companions were selected from processed MARVELS data using the latest University of Florida Two Dimensional (UF2D) pipeline, which shows significant improvement and reduction of systematic errors over previous pipelines. The 10 brown dwarf companions range in mass from ~13 to 76 M$_{\text{Jup}}$ and have orbital radii of less than 1 AU. The two stellar companions have minimum masses of ~98 and 100 M$_{\text{Jup}}$. The host stars of the MARVELS brown dwarf sample have a mean metallicity of [Fe/H] = 0.03 $\pm$ 0.08 dex. Given our stellar sample we estimate the brown dwarf occurrence rate around solar-type stars with periods less than ~300 days to be ~0.56%.

Correlations between stellar chemistry and kinematics have long been used to gain insight into the evolution of the Milky Way Galaxy. Orbital angular momentum is a key physical parameter and it is often estimated from three-dimensional space motions. We here demonstrate the lower uncertainties that can be achieved in the estimation of one component of velocity through selection of stars in key directions and use of line-of-sight velocity alone (i.e. without incorporation of proper motion data). In this first paper we apply our technique to stars observed in the direction of Galactic rotation in the APOGEE survey. We first derive the distribution of azimuthal velocities, $v_\phi$, then from these and observed radial coordinates, estimate the stellar guiding centre radii, $R_g$, within $6.9\leq R \leq 10$ kpc with uncertainties smaller than (or of the order of) 1kpc. We show that there is no simple way to select a clean stellar sample based on low errors on proper motions and distances to obtain high-quality 3D velocities and hence one should pay particular attention when trying to identify kinematically peculiar stars based on velocities derived using the proper motions. Using our azimuthal velocity estimations, we investigate the joint distribution of elemental abundances and rotational kinematics free from the blurring effects of epicyclic motions, and we derive the $\partial v_\phi / \partial R_g$ trends for the thin and thick discs as a function of radius. Our analysis provides further evidence for radial migration within the thin disc and hints against radial migration playing a significant role in the evolution of the thick disc.

Using combined asteroseismic and spectroscopic observations of 418 red-giant stars close to the Galactic disc plane (6 kpc $<R_{\rm Gal}\lesssim13$ kpc, $|Z_{\rm Gal}|<0.3$ kpc), we measure the age dependence of the radial metallicity distribution in the Milky Way's thin disc over cosmic time. The slope of the radial iron gradient of the young red-giant population ($-0.058\pm0.008$ [stat.] $\pm0.003$ [syst.] dex/kpc) is consistent with recent Cepheid measurements. For stellar populations with ages of $1-4$ Gyr the gradient is slightly steeper, at a value of $-0.066\pm0.007\pm0.002$ dex/kpc, and then flattens again to reach a value of $\sim-0.03$ dex/kpc for stars with ages between 6 and 10 Gyr. Our results are in good agreement with a state-of-the-art chemo-dynamical Milky-Way model in which the evolution of the abundance gradient and its scatter can be entirely explained by a non-varying negative metallicity gradient in the interstellar medium, together with stellar radial heating and migration. We also offer an explanation for why intermediate-age open clusters in the Solar Neighbourhood can be more metal-rich, and why their radial metallicity gradient seems to be much steeper than that of the youngest clusters. Already within 2 Gyr, radial mixing can bring metal-rich clusters from the innermost regions of the disc to Galactocentric radii of 5 to 8 kpc. We suggest that these outward-migrating clusters may be less prone to tidal disruption and therefore steepen the local intermediate-age cluster metallicity gradient. Our scenario also explains why the strong steepening of the local iron gradient with age is not seen in field stars. In the near future, asteroseismic data from the K2 mission will allow for improved statistics and a better coverage of the inner-disc regions, thereby providing tighter constraints on the evolution of the central parts of the Milky Way.

We present a study of quasar selection using the DES supernova fields. We used a quasar catalog from an overlapping portion of the SDSS Stripe 82 region to quantify the completeness and efficiency of selection methods involving color, probabilistic modeling, variability, and combinations of color/probabilistic modeling with variability. We only considered objects that appear as point sources in the DES images. We examine color selection methods based on the WISE mid-IR W1-W2 color, a mixture of WISE and DES colors (g-i and i-W1) and a mixture of VHS and DES colors (g-i and i-K). For probabilistic quasar selection, we used XDQSOz, an algorithm that employs an empirical multi-wavelength flux model of quasars to assign quasar probabilities. Our variability selection uses the multi-band chi2-probability that sources are constant in the DES Year 1 griz-band light curves. The completeness and efficiency are calculated relative to an underlying sample of point sources that are detected in the required selection bands and pass our data quality and photometric error cuts. We conduct our analyses at two magnitude limits, i<19.8 mag and i<22 mag. For sources with W1 and W2 detections, the W1-W2 color or XDQSOz method combined with variability gives the highest completenesses of >85% for both i-band magnitude limits and efficiencies of >80% to the bright limit and >60% to the faint limit; however, the giW1 and giW1+variability methods give the highest quasar surface densities. The XDQSOz method and combinations of W1W2/giW1/XDQSOz with variability are among the better selection methods when both high completeness and high efficiency are desired. We also present the OzDES Quasar Catalog of 1,263 spectroscopically-confirmed quasars taken by the OzDES survey. The catalog includes quasars with redshifts up to z~4 and brighter than i=22 mag, although the catalog is not complete up this magnitude limit.

We present the discovery and preliminary characterization of a gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy. Since the images of gravitationally lensed quasars are the superposition of multiple point sources and a foreground lensing galaxy, we have developed a morphology independent multi-wavelength approach to the photometric selection of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised machine learning. Using this technique and $gi$ multicolour photometric observations from the Dark Energy Survey (DES), near IR $JK$ photometry from the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have identified a candidate system with two catalogue components with $i_{AB}=18.61$ and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point sources. Spectroscopic follow-up with NTT and the use of an archival AAT spectrum show that the point sources can be identified as a lensed quasar with an emission line redshift of $z=2.739\pm0.003$ and a foreground early type galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass $M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days. The relatively wide separation, month scale time delay duration and high redshift make this an ideal system for constraining the expansion rate beyond a redshift of 1.

We search for excess gamma-ray emission coincident with the positions of confirmed and candidate Milky Way satellite galaxies using 6 years of data from the Fermi Large Area Telescope (LAT). Our sample of 45 stellar systems includes 28 kinematically confirmed dark-matter-dominated dwarf spheroidal galaxies (dSphs) and 17 recently discovered systems that have photometric characteristics consistent with the population of known dSphs. For each of these targets, the relative predicted gamma-ray flux due to dark matter annihilation is taken from kinematic analysis if available, and estimated from a distance-based scaling relation otherwise, assuming that the stellar systems are dark-matter-dominated dSphs. LAT data coincident with four of the newly discovered targets show a slight preference (each ~$2 \sigma$ local) for gamma-ray emission in excess of the background. However, the ensemble of derived gamma-ray flux upper limits for individual targets is consistent with the expectation from analyzing random blank-sky regions, and a combined analysis of the population of stellar systems yields no globally significant excess (global significance $<1 \sigma$). Our analysis has increased sensitivity compared to the analysis of 15 confirmed dSphs by Ackermann et al. 2015. The observed constraints on the dark matter annihilation cross section are statistically consistent with the background expectation, improving by a factor of ~2 for large dark matter masses ($m_{{\rm DM},b \bar b} \gtrsim 1$ TeV and $m_{{\rm DM},\tau^{+}\tau^{-}} \gtrsim 70$ GeV) and weakening by a factor of ~1.5 at lower masses relative to previously observed limits.

The galaxy formation process in the $\Lambda$-Cold Dark Matter scenario can be constrained from the analysis of stars in the Milky Way's halo system. We examine the variation of chemical abundances in distant halo stars observed by the Apache Point Galactic Evolution Experiment (APOGEE), as a function of distance from the Galactic center ($r$) and iron abundance ([M/H]), in the range 5 $\lesssim r \lesssim$ 30 kpc and $-2.5 <$ [M/H] $<$ 0.0. We perform a statistical analysis of the abundance ratios derived by the APOGEE pipeline (ASPCAP) and distances calculated by several approaches. Our analysis reveals signatures of a different chemical enrichment between the inner and outer regions of the halo, with a transition at about 15 kpc. The derived metallicity distribution function exhibits two peaks, at [M/H] $\sim -1.5$ and $\sim -2.1$, consistent with previously reported halo metallicity distributions. We obtain a difference of $\sim 0.1$ dex for $\alpha$-element-to-iron ratios for stars at $r > 15$ kpc and [M/H] $> -1.1$ (larger in the case of O, Mg and S) with respect to the nearest halo stars. This result confirms previous claims for low-$\alpha$ stars found at larger distances. Chemical differences in elements with other nucleosynthetic origins (Ni, K, Na, and Al) are also detected. C and N do not provide reliable information about the interstellar medium from which stars formed because our sample comprises RGB and AGB stars and can experience mixing of material to their surfaces.

The correlation between active galactic nuclei (AGN) and environment provides important clues to AGN fueling and the relationship of black hole growth to galaxy evolution. In this paper, we analyze the fraction of galaxies in clusters hosting AGN as a function of redshift and cluster richness for X-ray detected AGN associated with clusters of galaxies in Dark Energy Survey (DES) Science Verification data. The present sample includes 33 AGN with L_X > 10^43 ergs s^-1 in non-central, host galaxies with luminosity greater than 0.5 L* from a total sample of 432 clusters in the redshift range of 0.1<z<0.95. Analysis of the present sample reveals that the AGN fraction in red-sequence cluster members has a strong positive correlation with redshift such that the AGN fraction increases by a factor of ~8 from low to high redshift, and the fraction of cluster galaxies hosting AGN at high redshifts is greater than the low-redshift fraction at 3.6 sigma. In particular, the AGN fraction increases steeply at the highest redshifts in our sample at z>0.7. This result is in good agreement with previous work and parallels the increase in star formation in cluster galaxies over the same redshift range. However, the AGN fraction in clusters is observed to have no significant correlation with cluster mass. Future analyses with DES Year 1 through Year 3 data will be able to clarify whether AGN activity is correlated to cluster mass and will tightly constrain the relationship between cluster AGN populations and redshift.

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

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$\alpha$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$\alpha$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

We present the discovery of a z=0.65 low-ionization broad absorption line (LoBAL) quasar in a post-starburst galaxy in data from the Dark Energy Survey (DES) and spectroscopy from the Australian Dark Energy Survey (OzDES). LoBAL quasars are a minority of all BALs, and rarer still is that this object also exhibits broad FeII (an FeLoBAL) and Balmer absorption. This is the first BAL quasar that has signatures of recently truncated star formation, which we estimate ended about 40 Myr ago. The characteristic signatures of an FeLoBAL require high column densities, which could be explained by the emergence of a young quasar from an early, dust-enshrouded phase, or by clouds compressed by a blast wave. The age of the starburst component is comparable to estimates of the lifetime of quasars, so if we assume the quasar activity is related to the truncation of the star formation, this object is better explained by the blast wave scenario.

We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains 786 clusters with richness $\lambda>20$ (roughly equivalent to $M_{\rm{500c}}\gtrsim10^{14}\,h_{70}^{-1}\,M_{\odot}$) and $0.2<z<0.9$. The DR8 catalog consists of 26311 clusters with $0.08<z<0.6$, with a sharply increasing richness threshold as a function of redshift for $z\gtrsim 0.35$. The photometric redshift performance of both catalogs is shown to be excellent, with photometric redshift uncertainties controlled at the $\sigma_z/(1+z)\sim 0.01$ level for $z\lesssim0.7$, rising to $\sim0.02$ at $z\sim0.9$ in DES SV. We make use of \emph{Chandra} and \emph{XMM} X-ray and South Pole Telescope Sunyaev-Zeldovich data to show that the centering performance and mass--richness scatter are consistent with expectations based on prior runs of \redmapper{} on SDSS data. We also show how the \redmapper{} \photoz{} and richness estimates are relatively insensitive to imperfect star/galaxy separation and small-scale star masks.

We study the galaxy populations in 74 Sunyaev Zeldovich Effect (SZE) selected clusters from the South Pole Telescope (SPT) survey that have been imaged in the science verification phase of the Dark Energy Survey (DES). The sample extends up to $z\sim 1.1$ with $4 \times 10^{14} M_{\odot}\le M_{200}\le 3\times 10^{15} M_{\odot}$. Using the band containing the 4000~\AA\ break and its redward neighbor, we study the color-magnitude distributions of cluster galaxies to $\sim m_*+2$, finding: (1) the intrinsic rest frame $g-r$ color width of the red sequence (RS) population is $\sim$0.03 out to $z\sim0.85$ with a preference for an increase to $\sim0.07$ at $z=1$ and (2) the prominence of the RS declines beyond $z\sim0.6$. The spatial distribution of cluster galaxies is well described by the NFW profile out to $4R_{200}$ with a concentration of $c_{\mathrm{g}} = 3.59^{+0.20}_{-0.18}$, $5.37^{+0.27}_{-0.24}$ and $1.38^{+0.21}_{-0.19}$ for the full, the RS and the blue non-RS populations, respectively, but with $\sim40$\% to 55\% cluster to cluster variation and no statistically significant redshift or mass trends. The number of galaxies within the virial region $N_{200}$ exhibits a mass trend indicating that the number of galaxies per unit total mass is lower in the most massive clusters, and shows no significant redshift trend. The red sequence (RS) fraction within $R_{200}$ is $(68\pm3)$\% at $z=0.46$, varies from $\sim$55\% at $z=1$ to $\sim$80\% at $z=0.1$, and exhibits intrinsic variation among clusters of $\sim14$\%. We discuss a model that suggests the observed redshift trend in RS fraction favors a transformation timescale for infalling field galaxies to become RS galaxies of 2 to 3~Gyr.

The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates.

We report the discovery of an excess of main sequence turn-off stars in the direction of the constellations of Eridanus and Phoenix from the first year data of the Dark Energy Survey (DES). The Eridanus-Phoenix (EriPhe) overdensity is centered around l~285 deg and b~-60 deg and spans at least 30 deg in longitude and 10 deg in latitude. The Poisson significance of the detection is at least 9 sigma. The stellar population in the overdense region is similar in brightness and color to that of the nearby globular cluster NGC 1261, indicating that the heliocentric distance of EriPhe is about d~16 kpc. The extent of EriPhe in projection is therefore at least ~4 kpc by ~3 kpc. On the sky, this overdensity is located between NGC 1261 and a new stellar stream discovered by DES at a similar heliocentric distance, the so-called Phoenix Stream. Given their similar distance and proximity to each other, it is possible that these three structures may be kinematically associated. Alternatively, the EriPhe overdensity is morphologically similar to the Virgo overdensity and the Hercules-Aquila cloud, which also lie at a similar Galactocentric distance. These three overdensities lie along a polar plane separated by ~120 deg and may share a common origin. Spectroscopic follow-up observations of the stars in EriPhe are required to fully understand the nature of this overdensity.

We report the discovery of a stellar stream in the Dark Energy Survey (DES) Year 1 (Y1A1) data. The discovery was made through simple color-magnitude filters and visual inspection of the Y1A1 data. We refer to this new object as the Phoenix stream, after its resident constellation. After subtraction of the background stellar population we detect a clear signal of a simple stellar population. By fitting the ridge line of the stream in color-magnitude space, we find that a stellar population with age $\tau=11.5\pm0.5$ Gyr and $[Fe/H]<-1.6$ located 17.5$\pm$0.9 kpc from the Sun gives an adequate description of the stream stellar population. The stream is detected over an extension of 8$^{\circ}.$1 (2.5 kpc) and has a width of $\sim$54 pc assuming a Gaussian profile, indicating that a globular cluster is a probable progenitor. There is no known globular cluster within 5 kpc compatible with being the progenitor of the stream, assuming that the stream traces its orbit. We examined overdensities along the stream, however no obvious counterpart bound stellar system is visible in the coadded images. We also find overdensities along the stream that appear to be symmetrically distributed - consistent with the epicyclic overdensity scenario for the formation of cold streams - as well as a misalignment between the Northern and Southern part of stream. Despite the close proximity we find no evidence that this stream and the halo cluster NGC 1261 have a common accretion origin linked to the recently found EriPhe overdensity (Li et al. 2016).

We report the discovery of 13 confirmed two-image quasar lenses from a systematic search for gravitationally lensed quasars in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We adopted a methodology similar to that used in the SDSS Quasar Lens Search (SQLS). In addition to the confirmed lenses, we report 11 quasar pairs with small angular separations ($\lesssim$2") confirmed from our spectroscopy, which are either projected pairs, physical binaries, or possibly quasar lens systems whose lens galaxies have not yet been detected. The newly discovered quasar lens system, SDSS J1452+4224 at zs$\approx$4.8 is one of the highest redshift multiply imaged quasars found to date. Furthermore, we have over 50 good lens candidates yet to be followed up. Owing to the heterogeneous selection of BOSS quasars, the lens sample presented here does not have a well-defined selection function.

We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $\Delta z = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing results using a template-based photo-$z$ algorithm (BPZ) to a machine-learning algorithm (TPZ). A companion paper (Leistedt et al 2015) presents maps of several observational variables (e.g. seeing, sky brightness) which could modulate the galaxy density. Here we characterize and mitigate systematic errors on the measured clustering which arise from these observational variables, in addition to others such as Galactic dust and stellar contamination. After correcting for systematic effects we measure galaxy bias over a broad range of linear scales relative to mass clustering predicted from the Planck $\Lambda$CDM model, finding agreement with CFHTLS measurements with $\chi^2$ of 4.0 (8.7) with 5 degrees of freedom for the TPZ (BPZ) redshifts. We test a "linear bias" model, in which the galaxy clustering is a fixed multiple of the predicted non-linear dark-matter clustering. The precision of the data allow us to determine that the linear bias model describes the observed galaxy clustering to $2.5\%$ accuracy down to scales at least $4$ to $10$ times smaller than those on which linear theory is expected to be sufficient.

We present the first results of a survey for high redshift, z $\ge$ 6, quasars using izY multi-colour photometric observations from the Dark Energy Survey (DES). Here we report the discovery and spectroscopic confirmation of the $\rm z_{AB}, Y_{AB}$ = 20.2, 20.2 (M$_{1450}$ = $-$26.5) quasar DES J0454$-$4448 with an emission line redshift of z = 6.10$\pm$0.03 and a HI near zone size of 4.6 $\pm$ 1.7 Mpc.The quasar was selected as an i-band drop out with i$-$z = 2.46 and z$_{AB} < 21.5$ from an area of $\rm \sim$300 deg$^2$. It is the brightest of our 43 candidates and was identified for follow-up spectroscopically solely based on the DES i$-$z and z$-$Y colours. The quasar is detected by WISE and has $W1_{AB} = 19.68$. The discovery of one spectroscopically confirmed quasar with 5.7 $<$ z $<$ 6.5 and z$_{AB} \leq$ 20.2 is consistent with recent determinations of the luminosity function at z $\sim$ 6. DES when completed will have imaged $\rm \sim$5000 deg$^2$ to $Y_{AB}$ = 23.0 ($5\sigma$ point source) and we expect to discover $>$ 50-100 new quasars with z $>$ 6 including 3-10 with z $>$ 7 dramatically increasing the numbers of quasars currently known that are suitable for detailed studies including determination of the neutral HI fraction of the intergalactic medium (IGM) during the epoch of Hydrogen reionization.

The Large Magellanic Cloud (LMC) harbors a rich and diverse system of star clusters, whose ages, chemical abundances, and positions provide information about the LMC history of star formation. We use Science Verification imaging data from the Dark Energy Survey to increase the census of known star clusters in the outer LMC and to derive physical parameters for a large sample of such objects using a spatially and photometrically homogeneous data set. Our sample contains 255 visually identified cluster candidates, of which 109 were not listed in any previous catalog. We quantify the crowding effect for the stellar sample produced by the DES Data Management pipeline and conclude that the stellar completeness is < 10% inside typical LMC cluster cores. We therefore develop a pipeline to sample and measure stellar magnitudes and positions around the cluster candidates using DAOPHOT. We also implement a maximum-likelihood method to fit individual density profiles and colour-magnitude diagrams. For 117 (from a total of 255) of the cluster candidates (28 uncatalogued clusters), we obtain reliable ages, metallicities, distance moduli and structural parameters, confirming their nature as physical systems. The distribution of cluster metallicities shows a radial dependence, with no clusters more metal-rich than [Fe/H] ~ -0.7 beyond 8 kpc from the LMC center. The age distribution has two peaks at ~ 1.2 Gyr and ~ 2.7 Gyr.

We present the results of the first test plates of the extended Baryon Oscillation Spectroscopic Survey. This paper focuses on the emission line galaxies (ELG) population targetted from the Dark Energy Survey (DES) photometry. We analyse the success rate, efficiency, redshift distribution, and clustering properties of the targets. From the 9000 spectroscopic redshifts targetted, 4600 have been selected from the DES photometry. The total success rate for redshifts between 0.6 and 1.2 is 71\% and 68\% respectively for a bright and faint, on average more distant, samples including redshifts measured from a single strong emission line. We find a mean redshift of 0.8 and 0.87, with 15 and 13\% of unknown redshifts respectively for the bright and faint samples. In the redshift range 0.6<z<1.2, for the most secure spectroscopic redshifts, the mean redshift for the bright and faint sample is 0.85 and 0.9 respectively. Star contamination is lower than 2\%. We measure a galaxy bias averaged on scales of 1 and 10~Mpc/h of 1.72 \pm 0.1 for the bright sample and of 1.78 \pm 0.12 for the faint sample. The error on the galaxy bias have been obtained propagating the errors in the correlation function to the fitted parameters. This redshift evolution for the galaxy bias is in agreement with theoretical expectations for a galaxy population with MB-5\log h < -21.0. We note that biasing is derived from the galaxy clustering relative to a model for the mass fluctuations. We investigate the quality of the DES photometric redshifts and find that the outlier fraction can be reduced using a comparison between template fitting and neural network, or using a random forest algorithm.

We report the discovery of eight new ultra-faint dwarf galaxy candidates in the second year of optical imaging data from the Dark Energy Survey (DES). Six of these candidates are detected at high confidence, while two lower-confidence candidates are identified in regions of non-uniform survey coverage. The new stellar systems are found by three independent automated search techniques and are identified as overdensities of stars, consistent with the isochrone and luminosity function of an old and metal-poor simple stellar population. The new systems are faint (Mv > -4.7 mag) and span a range of physical sizes (17 pc < $r_{1/2}$ < 181 pc) and heliocentric distances (25 kpc < D < 214 kpc). All of the new systems have central surface brightnesses consistent with known ultra-faint dwarf galaxies (\mu < 27.5 mag arcsec$^{-2}$). Roughly half of the DES candidates are more distant, less luminous, and/or have lower surface brightnesses than previously known Milky Way satellite galaxies. Most of the candidates are found in the southern part of the DES footprint close to the Magellanic Clouds. We find that the DES data alone exclude (p < 0.001) a spatially isotropic distribution of Milky Way satellites and that the observed distribution can be well, though not uniquely, described by an association between several of the DES satellites and the Magellanic system. Our model predicts that the full sky may hold ~100 ultra-faint galaxies with physical properties comparable to the DES satellites and that 20-30% of these would be spatially associated with the Magellanic Clouds.

We developed a code that estimates distances to stars using measured spectroscopic and photometric quantities. We employ a Bayesian approach to build the probability distribution function over stellar evolutionary models given these data, delivering estimates of model parameters for each star individually. The code was first tested on simulations, successfully recovering input distances to mock stars with <1% bias.The method-intrinsic random distance uncertainties for typical spectroscopic survey measurements amount to around 10% for dwarf stars and 20\% for giants, and are most sensitive to the quality of $\log g$ measurements. The code was validated by comparing our distance estimates to parallax measurements from the Hipparcos mission for nearby stars (< 300 pc), to asteroseismic distances of CoRoT red giant stars, and to known distances of well-studied open and globular clusters. The external comparisons confirm that our distances are subject to very small systematic biases with respect to the fundamental Hipparcos scale (+0.4 % for dwarfs, and +1.6% for giants). The typical random distance scatter is 18% for dwarfs, and 26% for giants. For the CoRoT-APOGEE sample, the typical random distance scatter is ~15%, both for the nearby and farther data. Our distances are systematically larger than the CoRoT ones by about +9%, which can mostly be attributed to the different choice of priors. The comparison to known distances of star clusters from SEGUE and APOGEE has led to significant systematic differences for many cluster stars, but with opposite signs, and with substantial scatter. Finally, we tested our distances against those previously determined for a high-quality sample of giant stars from the RAVE survey, again finding a small systematic trend of +5% and an rms scatter of 30%.

Tidal streams of globular clusters are ideal tracers of the Galactic gravitational potential. Compared to the few known, complex and diffuse dwarf-galaxy streams, they are kinematically cold, have thin morphologies and are abundant in the halo of the Milky Way. Their coldness and thinness in combination with potential epicyclic substructure in the vicinity of the stream progenitor turns them into high-precision scales. With the example of Palomar 5, we demonstrate how modeling of a globular cluster stream allows us to simultaneously measure the properties of the disrupting globular cluster, its orbital motion, and the gravitational potential of the Milky Way.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

The Dark Energy Survey (DES) is a 5000 sq. degree survey in the southern hemisphere, which is rapidly reducing the existing north-south asymmetry in the census of MW satellites and other stellar substructure. We use the first-year DES data down to previously unprobed photometric depths to search for stellar systems in the Galactic halo, therefore complementing the previous analysis of the same data carried out by our group earlier this year. Our search is based on a matched filter algorithm that produces stellar density maps consistent with stellar population models of various ages, metallicities, and distances over the survey area. The most conspicuous density peaks in these maps have been identified automatically and ranked according to their significance and recurrence for different input models. We report the discovery of one additional stellar system besides those previously found by several authors using the same first-year DES data. The object is compact, and consistent with being dominated by an old and metal-poor population. DES J0034-4902 is found at high significance and appears in the DES images as a compact concentration of faint blue point sources at ~ 87 {kpc}. Its half-light radius of r_h = 9.88 +/- 4.31 {pc} and total luminosity of M_V ~ -3.05_{-0.42}^{+0.69} are consistent with it being a low mass halo cluster. It is also found to have a very elongated shape. In addition, our deeper probe of DES 1st year data confirms the recently reported satellite galaxy candidate Horologium II as a significant stellar overdensity. We also infer its structural properties and compare them to those reported in the literature.

We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

We obtained SOAR telescope B and V photometry of 14 star clusters and 2 associations in the Bridge tidal structure connecting the LMC and SMC. These objects are used to study the formation and evolution of star clusters and associations under tidal stresses from the Clouds. Typical star clusters in the Bridge are not richly populated and have in general relatively large diameters (~30-35 pc), being larger than Galactic counterparts of similar age. Ages and other fundamental parameters are determined with field-star decontaminated photometry. A self-consistent approach is used to derive parameters for the most-populated sample cluster NGC 796 and two young CMD templates built with the remaining Bridge clusters. We find that the clusters are not coeval in the Bridge. They range from approximately a few Myr (still related to optical HII regions and WISE and Spitzer dust emission measurements) to about 100-200 Myr. The derived distance moduli for the Bridge objects suggests that the Bridge is a structure connecting the LMC far-side in the East to the foreground of the SMC to the West. Most of the present clusters are part of the tidal dwarf candidate D 1, which is associated with an H I overdensity. We find further evidence that the studied part of the Bridge is evolving into a tidal dwarf galaxy, decoupling from the Bridge.

We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based photometric redshift methods -- ANNZ2, BPZ calibrated against BCC-Ufig simulations, SkyNet, and TPZ -- are analysed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-zs. From the galaxies in the DES SV shear catalogue, which have mean redshift $0.72\pm0.01$ over the range $0.3<z<1.3$, we construct three tomographic bins with means of $z=\{0.45, 0.67, 1.00\}$. These bins each have systematic uncertainties $\delta z \lesssim 0.05$ in the mean of the fiducial SkyNet photo-z $n(z)$. We propagate the errors in the redshift distributions through to their impact on cosmological parameters estimated with cosmic shear, and find that they cause shifts in the value of $\sigma_8$ of approx. 3%. This shift is within the one sigma statistical errors on $\sigma_8$ for the DES SV shear catalog. We further study the potential impact of systematic differences on the critical surface density, $\Sigma_{\mathrm{crit}}$, finding levels of bias safely less than the statistical power of DES SV data. We recommend a final Gaussian prior for the photo-z bias in the mean of $n(z)$ of width $0.05$ for each of the three tomographic bins, and show that this is a sufficient bias model for the corresponding cosmology analysis.

Spatially-varying depth and characteristics of observing conditions, such as seeing, airmass, or sky background, are major sources of systematic uncertainties in modern galaxy survey analyses, in particular in deep multi-epoch surveys. We present a framework to extract and project these sources of systematics onto the sky, and apply it to the Dark Energy Survey (DES) to map the observing conditions of the Science Verification (SV) data. The resulting distributions and maps of sources of systematics are used in several analyses of DES SV to perform detailed null tests with the data, and also to incorporate systematics in survey simulations. We illustrate the complementarity of these two approaches by comparing the SV data with the BCC-UFig, a synthetic sky catalogue generated by forward-modelling of the DES SV images. We analyse the BCC-UFig simulation to construct galaxy samples mimicking those used in SV galaxy clustering studies. We show that the spatially-varying survey depth imprinted in the observed galaxy densities and the redshift distributions of the SV data are successfully reproduced by the simulation and well-captured by the maps of observing conditions. The combined use of the maps, the SV data and the BCC-UFig simulation allows us to quantify the impact of spatial systematics on $N(z)$, the redshift distributions inferred using photometric redshifts. We conclude that spatial systematics in the SV data are mainly due to seeing fluctuations and are under control in current clustering and weak lensing analyses. The framework presented here is relevant to all multi-epoch surveys, and will be essential for exploiting future surveys such as the Large Synoptic Survey Telescope (LSST), which will require detailed null-tests and realistic end-to-end image simulations to correctly interpret the deep, high-cadence observations of the sky.

We introduce redMaGiC, an automated algorithm for selecting Luminous Red Galaxies (LRGs). The algorithm was specifically developed to minimize photometric redshift uncertainties in photometric large-scale structure studies. redMaGiC achieves this by self-training the color-cuts necessary to produce a luminosity-thresholded LRG sample of constant comoving density. We demonstrate that redMaGiC photozs are very nearly as accurate as the best machine-learning based methods, yet they require minimal spectroscopic training, do not suffer from extrapolation biases, and are very nearly Gaussian. We apply our algorithm to Dark Energy Survey (DES) Science Verification (SV) data to produce a redMaGiC catalog sampling the redshift range $z\in[0.2,0.8]$. Our fiducial sample has a comoving space density of $10^{-3}\ (h^{-1} Mpc)^{-3}$, and a median photoz bias ($z_{spec}-z_{photo}$) and scatter $(\sigma_z/(1+z))$ of 0.005 and 0.017 respectively. The corresponding $5\sigma$ outlier fraction is 1.4%. We also test our algorithm with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) and Stripe 82 data, and discuss how spectroscopic training can be used to control photoz biases at the 0.1% level.

We report the discovery of the eighth and ninth known Trojans in stable orbits around Neptune's leading Lagrange point, L4. The objects 2014 QO$_{441}$ and 2014 QP$_{441}$ were detected in data obtained during the 2013-14 and 2014-15 observing seasons by the Dark Energy Survey, using the Dark Energy Camera (DECam) on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory. Both are in high-inclination orbits (18.8$^{\circ}$ and 19.4$^{\circ}$ respectively). With an eccentricity of 0.104, 2014 QO$_{441}$ has the most eccentric orbit of the eleven known stable Neptune Trojans. Here we describe the search procedure and investigate the objects' long-term dynamical stability and physical properties.

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.

We report the discovery of eight new Milky Way companions in ~1,800 deg^2 of optical imaging data collected during the first year of the Dark Energy Survey (DES). Each system is identified as a statistically significant over-density of individual stars consistent with the expected isochrone and luminosity function of an old and metal-poor stellar population. The objects span a wide range of absolute magnitudes (M_V from -2.2 mag to -7.4 mag), physical sizes (10 pc to 170 pc), and heliocentric distances (30 kpc to 330 kpc). Based on the low surface brightnesses, large physical sizes, and/or large Galactocentric distances of these objects, several are likely to be new ultra-faint satellite galaxies of the Milky Way and/or Magellanic Clouds. We introduce a likelihood-based algorithm to search for and characterize stellar over-densities, as well as identify stars with high satellite membership probabilities. We also present completeness estimates for detecting ultra-faint galaxies of varying luminosities, sizes, and heliocentric distances in the first-year DES data.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We address the problem of separating stars from galaxies in future large photometric surveys. We focus our analysis on simulations of the Dark Energy Survey (DES). In the first part of the paper, we derive the science requirements on star/galaxy separation, for measurement of the cosmological parameters with the Gravitational Weak Lensing and Large Scale Structure probes. These requirements are dictated by the need to control both the statistical and systematic errors on the cosmological parameters, and by Point Spread Function calibration. We formulate the requirements in terms of the completeness and purity provided by a given star/galaxy classifier. In order to achieve these requirements at faint magnitudes, we propose a new method for star/galaxy separation in the second part of the paper. We first use Principal Component Analysis to outline the correlations between the objects parameters and extract from it the most relevant information. We then use the reduced set of parameters as input to an Artificial Neural Network. This multi-parameter approach improves upon purely morphometric classifiers (such as the classifier implemented in SExtractor), especially at faint magnitudes: it increases the purity by up to 20% for stars and by up to 12% for galaxies, at i-magnitude fainter than 23.

We report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [alpha/Fe] abundance ratios (as determined from APOGEE spectra) with respect to Solar values. Their existence is not explained by standard chemical evolution models of the Milky Way, and shows that the chemical-enrichment history of the Galactic disc is more complex. We find similar stars in previously published samples for which isochrone-ages could be robustly obtained, although in smaller relative numbers, which could explain why these stars have not received prior attention. The young [alpha/Fe]-rich stars are much more numerous in the CoRoT-APOGEE (CoRoGEE) inner-field sample than in any other high-resolution sample available at present, as only CoRoGEE can explore the inner-disc regions and provide ages for its field stars. The kinematic properties of the young [$\alpha$/Fe]-rich stars are not clearly thick-disc like, despite their rather large distances from the Galactic mid-plane. Our tentative interpretation of these and previous intriguing observations in the Milky Way is that these stars were formed close to the end of the Galactic bar, near corotation -- a region where gas can be kept inert for longer times, compared to other regions shocked more frequently by the passage of spiral arms. Moreover, that is where the mass return from older inner-disc stellar generations should be maximal (according to an inside-out disc-formation scenario), further diluting the in-situ gas. Other possibilities to explain these observations (e.g., a recent gas-accretion event) are also discussed.

We present a forward-modelling simulation framework designed to model the data products from the Dark Energy Survey (DES). This forward-model process can be thought of as a transfer function -- a mapping from cosmological and astronomical signals to the final data products used by the scientists. Using output from the cosmological simulations (the Blind Cosmology Challenge), we generate simulated images (the Ultra Fast Image Simulator, Berge et al. 2013) and catalogs representative of the DES data. In this work we simulate the 244 sq. deg coadd images and catalogs in 5 bands for the DES Science Verification (SV) data. The simulation output is compared with the corresponding data to show that major characteristics of the images and catalogs can be captured. We also point out several directions of future improvements. Two practical examples, star/galaxy classification and proximity effects on object detection, are then used to demonstrate how one can use the simulations to address systematics issues in data analysis. With clear understanding of the simplifications in our model, we show that one can use the simulations side-by-side with data products to interpret the measurements. This forward modelling approach is generally applicable for other upcoming and future surveys. It provides a powerful tool for systematics studies which is sufficiently realistic and highly controllable.

In the last few decades many efforts have been made to understand the effect of spiral arms on the gas and stellar dynamics in the Milky Way disc. One of the fundamental parameters of the spiral structure is its angular velocity, or pattern speed $\Omega_p$, which determines the location of resonances in the disc and the spirals' radial extent. The most direct method for estimating the pattern speed relies on backward integration techniques, trying to locate the stellar birthplace of open clusters. Here we propose a new method based on the interaction between the spiral arms and the stars in the disc. Using a sample of around 500 open clusters from the {\it New Catalogue of Optically Visible Open Clusters and Candidates}, and a sample of 500 giant stars observed by APOGEE, we find $\Omega_p = 23.0\pm0.5$ km s$^{-1}$ kpc$^{-1}$, for a local standard of rest rotation $V_0=220$~km s$^{-1}$ and solar radius $R_0=8.0$~kpc. Exploring a range in $V_0$ and $R_0$ within the acceptable values, 200-240 km s$^{-1}$ and 7.5-8.5 kpc, respectively, results only in a small change in our estimate of $\Omega_p$, that is within the error. Our result is in close agreement with a number of studies which suggest values in the range 20-25 km s$^{-1}$ kpc$^{-1}$. An advantage of our method is that we do not need knowledge of the stellar age, unlike in the case of the birthplace method, which allows us to use data from large Galactic surveys. The precision of our method will be improved once larger samples of disk stars with spectroscopic information will become available thanks to future surveys such as 4MOST.

We measure the weak-lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey. This pathfinder study is meant to 1) validate the DECam imager for the task of measuring weak-lensing shapes, and 2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, PSF modeling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well-behaved, but the modeling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting NFW profiles to the clusters in this study, we determine weak-lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak-lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.

The Dark Energy Camera has captured a large set of images as part of Science Verification (SV) for the Dark Energy Survey. The SV footprint covers a lar ge portion of the outer Large Magellanic Cloud (LMC), providing photometry 1.5 magnitudes fainter than the main sequence turn-off of the oldest LMC stel lar population. We derive geometrical and structural parameters for various stellar populations in the LMC disk. For the distribution of all LMC stars, we find an inclination of $i=-38.14^{\circ}\pm0.08^{\circ}$ (near side in the North) and a position angle for the line of nodes of $\theta_0=129.51^{\circ}\pm0.17^{\circ}$. We find that stars younger than $\sim 4$ Gyr are more centrally concentrated than older stars. Fitting a projected exponential disk shows that the scale radius of the old populations is $R_{>4 Gyr}=1.41\pm0.01$ kpc, while the younger population has $R_{<4 Gyr}=0.72\pm0.01$ kpc. Howe ver, the spatial distribution of the younger population deviates significantly from the projected exponential disk model. The distribution of old stars suggests a large truncation radius of $R_{t}=13.5\pm0.8$ kpc. If this truncation is dominated by the tidal field of the Galaxy, we find that the LMC is $\simeq 24^{+9}_{-6}$ times less massive than the encircled Galactic mass. By measuring the Red Clump peak magnitude and comparing with the best-fit LM C disk model, we find that the LMC disk is warped and thicker in the outer regions north of the LMC centre. Our findings may either be interpreted as a warped and flared disk in the LMC outskirts, or as evidence of a spheroidal halo component

Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models, and inference with MCMC. By generating $\approx10^7$ model streams, we show that the geometry of the Pal 5 debris yields powerful constraints on the solar position and motion, the Milky Way and Pal 5 itself. All 10 model parameters were allowed to vary over large ranges without additional prior information. Using only SDSS data and a few radial velocities from the literature, we find that the distance of the Sun from the Galactic Center is $8.30\pm0.25$ kpc, and the transverse velocity is $253\pm16$ km/s. Both estimates are in excellent agreement with independent measurements of these quantities. Assuming a standard disk and bulge model, we determine the Galactic mass within Pal 5's apogalactic radius of 19 kpc to be $(2.1\pm0.4)\times10^{11}$ M$_\odot$. Moreover, we find the potential of the dark halo with a flattening of $q_z = 0.95^{+0.16}_{-0.12}$ to be essentially spherical within the radial range that is effectively probed by Pal 5. We also determine Pal 5's mass, distance and proper motion independently from other methods, which enables us to perform vital cross-checks. We conclude that with more observational data and by using additional prior information, the precision of this method can be significantly increased.

We present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova (SLSN) from the Dark Energy Survey (DES). We briefly discuss the data and search algorithm used to find this event in the first year of DES operations, and outline the spectroscopic data obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm its redshift (z = 0.663 +/- 0.001 based on the host-galaxy emission lines) and likely spectral type (type I). Using this redshift, we find M_U_peak = -21.05 +0.10 -0.09 for the peak, rest-frame U-band absolute magnitude, and find DES13S2cmm to be located in a faint, low metallicity (sub-solar), low stellar-mass host galaxy (log(M/M_sun) = 9.3 +/- 0.3); consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest declining tails (beyond +30 days rest frame past peak), and is the faintest at peak. Moreover, we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of only 0.2-0.3 magnitudes between +25 and +30 days after peak (rest frame) depending on redshift range studied; this could be important for 'standardising' such supernovae, as is done with the more common type Ia. We fit the bolometric light curve of DES13S2cmm with two competing models for SLSNe-I - the radioactive decay of 56Ni, and a magnetar - and find that while the magnetar is formally a better fit, neither model provides a compelling match to the data. Although we are unable to conclusively differentiate between these two physical models for this particular SLSN-I, further DES observations of more SLSNe-I should break this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 days in the rest frame of the supernova.