Search

The stars of the Milky Way carry the chemical history of our Galaxy in their atmospheres as they journey through its vast expanse. Like barcodes, we can extract the chemical fingerprints of stars from high-resolution spectroscopy. The fourth data release (DR4) of the Galactic Archaeology with HERMES (GALAH) Survey, based on a decade of observations, provides the chemical abundances of up to 32 elements for 917 588 stars that also have exquisite astrometric data from the $Gaia$ satellite. For the first time, these elements include life-essential nitrogen to complement carbon, and oxygen as well as more measurements of rare-earth elements critical to modern-life electronics, offering unparalleled insights into the chemical composition of the Milky Way. For this release, we use neural networks to simultaneously fit stellar parameters and abundances across the full spectrum, leveraging synthetic grids computed with Spectroscopy Made Easy. These grids account for atomic line formation in non-local thermodynamic equilibrium for 14 elements. In a two-iteration process, we first fit stellar labels for all 1 085 520 spectra, then co-add repeated observations and refine these labels using astrometric data from $Gaia$ and 2MASS photometry, improving the accuracy and precision of stellar parameters and abundances. Our validation thoroughly assesses the reliability of spectroscopic measurements and highlights key caveats for catalogue users. GALAH DR4 represents yet another milestone in Galactic archaeology, combining detailed chemical compositions from multiple nucleosynthetic channels with kinematic information and age estimates. The resulting dataset, covering nearly a million stars, opens new avenues for understanding not only the chemical and dynamical history of the Milky Way, but also the broader questions of the origin of elements and the evolution of planets, stars, and galaxies., Comment: 43 pages, 38 figures to be submitted to PASA. Accompanying the GALAH Data Release 4, see https://www.galah-survey.org and https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR4/. All code available on http://github.com/svenbuder/GALAH_DR4/ and https://github.com/svenbuder/galah_dr4_paper. Comments welcome

We present $griz$ photometric light curves for the full 5 years of the Dark Energy Survey Supernova program (DES-SN), obtained with both forced Point Spread Function (PSF) photometry on Difference Images (DIFFIMG) performed during survey operations, and Scene Modelling Photometry (SMP) on search images processed after the survey. This release contains $31,636$ DIFFIMG and $19,706$ high-quality SMP light curves, the latter of which contains $1635$ photometrically-classified supernovae that pass cosmology quality cuts. This sample spans the largest redshift ($z$) range ever covered by a single SN survey ($0.1

Type Ia Supernovae (SNe Ia) are a critical tool in measuring the accelerating expansion of the universe. Recent efforts to improve these standard candles have focused on incorporating the effects of dust on distance measurements with SNe Ia. In this paper, we use the state-of-the-art Dark Energy Survey 5 year sample to evaluate two different families of dust models: empirical extinction models derived from SNe Ia data, and physical attenuation models from the spectra of galaxies. Among the SNe Ia-derived models, we find that a logistic function of the total-to-selective extinction RV best recreates the correlations between supernova distance measurements and host galaxy properties, though an additional 0.02 magnitudes of grey scatter are needed to fully explain the scatter in SNIa brightness in all cases. These empirically-derived extinction distributions are highly incompatible with the physical attenuation models from galactic spectral measurements. From these results, we conclude that SNe Ia must either preferentially select extreme ends of galactic dust distributions, or that the characterisation of dust along the SNe Ia line-of-sight is incompatible with that of galactic dust distributions.

We present a precise measurement of cosmological time dilation using the light curves of 1504 type Ia supernovae from the Dark Energy Survey spanning a redshift range $0.1\lesssim z\lesssim 1.2$. We find that the width of supernova light curves is proportional to $(1+z)$, as expected for time dilation due to the expansion of the Universe. Assuming type Ia supernovae light curves are emitted with a consistent duration $\Delta t_{\rm em}$, and parameterising the observed duration as $\Delta t_{\rm obs}=\Delta t_{\rm em}(1+z)^b$, we fit for the form of time dilation using two methods. Firstly, we find that a power of $b \approx 1$ minimises the flux scatter in stacked subsamples of light curves across different redshifts. Secondly, we fit each target supernova to a stacked light curve (stacking all supernovae with observed bandpasses matching that of the target light curve) and find $b=1.003\pm0.005$ (stat) $\pm\,0.010$ (sys). Thanks to the large number of supernovae and large redshift-range of the sample, this analysis gives the most precise measurement of cosmological time dilation to date, ruling out any non-time-dilating cosmological models at very high significance., Comment: 14 pages, 13 figures. Updated in response to reviewer feedback. Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

We measure the current expansion rate of the Universe, Hubble's constant $H_0$, by calibrating the absolute magnitudes of supernovae to distances measured by Baryon Acoustic Oscillations. This `inverse distance ladder' technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning $0.01\lt z \lt1.13$ anchored to the recent Baryon Acoustic Oscillation measurements from DESI spanning $0.30 \lt z_{\mathrm{eff}} \lt 2.33$. To trace cosmology to $z=0$, we use the third-, fourth- and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher-redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find $H_0=67.19^{+0.66}_{-0.64}\mathrm{~km} \mathrm{~s}^{-1} \mathrm{~Mpc}^{-1}$, in agreement with the value found by Planck without the need to assume Flat-$\Lambda$CDM. However the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

Reverberation mapping is the leading technique used to measure direct black hole masses outside of the local Universe. Additionally, reverberation measurements calibrate secondary mass-scaling relations used to estimate single-epoch virial black hole masses. The Australian Dark Energy Survey (OzDES) conducted one of the first multi-object reverberation mapping surveys, monitoring 735 AGN up to $z\sim4$, over 6 years. The limited temporal coverage of the OzDES data has hindered recovery of individual measurements for some classes of sources, particularly those with shorter reverberation lags or lags that fall within campaign season gaps. To alleviate this limitation, we perform a stacking analysis of the cross-correlation functions of sources with similar intrinsic properties to recover average composite reverberation lags. This analysis leads to the recovery of average lags in each redshift-luminosity bin across our sample. We present the average lags recovered for the H$\beta$, Mg II and C IV samples, as well as multi-line measurements for redshift bins where two lines are accessible. The stacking analysis is consistent with the Radius-Luminosity relations for each line. Our results for the H$\beta$ sample demonstrate that stacking has the potential to improve upon constraints on the $R-L$ relation, which have been derived only from individual source measurements until now., Comment: 20 pages, 15 figures. Accepted by MNRAS

We present cosmological constraints from the sample of Type Ia supernovae (SN Ia) discovered during the full five years of the Dark Energy Survey (DES) Supernova Program. In contrast to most previous cosmological samples, in which SN are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being a SN Ia, we find 1635 DES SNe in the redshift range $0.100.5$ SNe compared to the previous leading compilation of Pantheon+, and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints we combine the DES supernova data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning $0.025

We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50\% of their stellar mass by $z=5$ (12.6~Gyr ago), similar to known UFDs around the MW, but that 10-40\% of their stellar mass formed at later times. We uncover one remarkable UFD, \A{XIII}, which formed only 10\% of its stellar mass by $z=5$, and 75\% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\times10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs., Comment: 18 pages, 14 figures, 5 appendices, accepted for publication in ApJ

Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This $R-L$ relation is used to estimate single-epoch virial black hole masses, and has been proposed for use to standardise AGN to determine cosmological distances. We present reverberation measurements made with H$\beta$ from the six-year Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at $0.12

The correlation between the broad line region radius and continuum luminosity ($R-L$ relation) of active galactic nuclei (AGN) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At $z \sim 1-2$, where AGN activity peaks, the $R-L$ relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey (OzDES) RM project based on six years of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the re-scaled, re-sampled versions of the observed lightcurves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the $R-L$ relation to higher redshifts and luminosities. The relative iron line strength $\mathcal{R}_{\rm Fe}$ has little impact on the $R-L$ relation. The best-fit Mg II $R-L$ relation has a slope $\alpha = 0.39 \pm 0.08$ with an intrinsic scatter $\sigma_{\rm rl} = 0.15^{+0.03}_{-0.02}$. The slope is consistent with previous measurements and shallower than the H$\beta$ $R-L$ relation. The intrinsic scatter of the new $R-L$ relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H$\beta$ $R-L$ relation. Our new $R-L$ relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon., Comment: 16 pages, 8 figures; MNRAS, Volume 522, pp.4132

Recent analyses have found intriguing correlations between the colour ($c$) of type Ia supernovae (SNe Ia) and the size of their 'mass-step', the relationship between SN Ia host galaxy stellar mass ($M_\mathrm{stellar}$) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically-classified SNe Ia from the Dark Energy Survey 5-year sample, we study the differences in Hubble residual for a variety of global host galaxy and local environmental properties for SN Ia subsamples split by their colour. We find a $3\sigma$ difference in the mass-step when comparing blue ($c<0$) and red ($c>0$) SNe. We observe the lowest r.m.s. scatter ($\sim0.14$ mag) in the Hubble residual for blue SNe in low mass/blue environments, suggesting that this is the most homogeneous sample for cosmological analyses. By fitting for $c$-dependent relationships between Hubble residuals and $M_\mathrm{stellar}$, approximating existing dust models, we remove the mass-step from the data and find tentative $\sim 2\sigma$ residual steps in rest-frame galaxy $U-R$ colour. This indicates that dust modelling based on $M_\mathrm{stellar}$ may not fully explain the remaining dispersion in SN Ia luminosity. Instead, accounting for a $c$-dependent relationship between Hubble residuals and global $U-R$, results in $\leq1\sigma$ residual steps in $M_\mathrm{stellar}$ and local $U-R$, suggesting that $U-R$ provides different information about the environment of SNe Ia compared to $M_\mathrm{stellar}$, and motivating the inclusion of galaxy $U-R$ colour in SN Ia distance bias correction., Comment: 19 pages, 8 figures. Published in MNRAS

We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically-confirmed supernovae over a redshift range $0.045

We use stacked spectra of the host galaxies of photometrically identified type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilising full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features ($< 1.3\sigma$), such as stellar population age, metallicity, and mass-to-light ratio. However, we find significant trends between the Hubble residuals and the strengths of [OII] ($4.4\sigma$) and the Balmer emission lines ($3\sigma$). These trends are weaker than the well known trend between Hubble residuals and host galaxy stellar mass ($7.2\sigma$) that is derived from broad band photometry. After light curve corrections, we see fainter SNe Ia residing in galaxies with larger line strengths. We also find a trend (3$\sigma$) between Hubble residual and the Balmer decrement (a measure of reddening by dust) using H${\beta}$ and H${\gamma}$. The trend, quantified by correlation coefficients, is slightly more significant in the redder SNe Ia, suggesting that bluer SNe Ia are relatively unaffected by dust in the interstellar medium of the host and that dust contributes to current Hubble diagram scatter impacting the measurement of cosmological parameters., Comment: 15 pages, 10 figures. Accepted for publication in MNRAS

Current and future cosmological analyses with Type Ia Supernovae (SNe Ia) face three critical challenges: i) measuring redshifts from the supernova or its host galaxy; ii) classifying SNe without spectra; and iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each challenge. In the context of the Dark Energy Survey (DES), we analyze a SNIa sample with host galaxies in the redMaGiC galaxy catalog, a selection of Luminous Red Galaxies. Photo-$z$ estimates for these galaxies are expected to be accurate to $\sigma_{\Delta z/(1+z)}\sim0.02$. The DES-5YR photometrically classified SNIa sample contains approximately 1600 SNe and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns with classification uncertainties. With this subsample, we find similar Hubble scatter (to within $\sim0.01$ mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show the bias due to using photo-$z$s from redMaGiC host galaxies on the measurement of the dark energy equation-of-state $w$ is up to $\Delta w \sim 0.01-0.02$. With real data, we measure a difference in $w$ when using redMaGiC photometric redshifts versus spectroscopic redshifts of $\Delta w = 0.005$. Finally, we discuss how SNe in redMaGiC galaxies appear to be a more standardizable population due to a weaker relation between color and luminosity ($\beta$) compared to the DES-3YR population by $\sim5\sigma$; this finding is consistent with predictions that redMaGiC galaxies exhibit lower reddening ratios ($\textrm{R}_\textrm{V}$) than the general population of SN host galaxies. These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data., Comment: 15 pages, 8 figures. Submitting to ApJ, comments welcome

As part of the cosmology analysis using Type Ia Supernovae (SN Ia) in the Dark Energy Survey (DES), we present photometrically identified SN Ia samples using multi-band light-curves and host galaxy redshifts. For this analysis, we use the photometric classification framework SuperNNova (SNN; M\"oller et al. 2019) trained on realistic DES-like simulations. For reliable classification, we process the DES SN programme (DES-SN) data and introduce improvements to the classifier architecture, obtaining classification accuracies of more than 98 per cent on simulations. This is the first SN classification to make use of ensemble methods, resulting in more robust samples. Using photometry, host galaxy redshifts, and a classification probability requirement, we identify 1,863 SNe Ia from which we select 1,484 cosmology-grade SNe Ia spanning the redshift range of 0.07 < z < 1.14. We find good agreement between the light-curve properties of the photometrically-selected sample and simulations. Additionally, we create similar SN Ia samples using two types of Bayesian Neural Network classifiers that provide uncertainties on the classification probabilities. We test the feasibility of using these uncertainties as indicators for out-of-distribution candidates and model confidence. Finally, we discuss the implications of photometric samples and classification methods for future surveys such as Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST)., Comment: 18 pages, 8 figures. Accepted in MNRAS

Cosmological analyses of samples of photometrically-identified Type Ia supernovae (SNe Ia) depend on understanding the effects of 'contamination' from core-collapse and peculiar SN Ia events. We employ a rigorous analysis on state-of-the-art simulations of photometrically identified SN Ia samples and determine cosmological biases due to such 'non-Ia' contamination in the Dark Energy Survey (DES) 5-year SN sample. As part of the analysis, we test on our DES simulations the performance of SuperNNova, a photometric SN classifier based on recurrent neural networks. Depending on the choice of non-Ia SN models in both the simulated data sample and training sample, contamination ranges from 0.8-3.5 %, with the efficiency of the classification from 97.7-99.5 %. Using the Bayesian Estimation Applied to Multiple Species (BEAMS) framework and its extension 'BEAMS with Bias Correction' (BBC), we produce a redshift-binned Hubble diagram marginalised over contamination and corrected for selection effects and we use it to constrain the dark energy equation-of-state, $w$. Assuming a flat universe with Gaussian $\Omega_M$ prior of $0.311\pm0.010$, we show that biases on $w$ are $<0.008$ when using SuperNNova and accounting for a wide range of non-Ia SN models in the simulations. Systematic uncertainties associated with contamination are estimated to be at most $\sigma_{w, \mathrm{syst}}=0.004$. This compares to an expected statistical uncertainty of $\sigma_{w,\mathrm{stat}}=0.039$ for the DES-SN sample, thus showing that contamination is not a limiting uncertainty in our analysis. We also measure biases due to contamination on $w_0$ and $w_a$ (assuming a flat universe), and find these to be $<$0.009 in $w_0$ and $<$0.108 in $w_a$, hence 5 to 10 times smaller than the statistical uncertainties expected from the DES-SN sample.

We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago., Comment: We welcome comments and missing references

We use six years of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15-degree long, thin, dynamically cold, low-metallicity stellar system in the southern hemisphere. We use natural splines, a non-parametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us to improve measurements of the heliocentric distance ($17.4 \pm 0.1\,{\rm (stat.)} \pm 0.8\,{\rm (sys.)}$ kpc) and distance gradient ($-0.009 \pm 0.006$ kpc deg$^{-1}$) of Phoenix, which corresponds to a small change of $0.13 \pm 0.09$ kpc in heliocentric distance along the length of the stream. We measure linear intensity variations on degree scales, as well as deviations in the stream track on $\sim 2$-degree scales, suggesting that the stream may have been disturbed during its formation and/or evolution. We recover three peaks and one gap in linear intensity along with fluctuations in the stream track. Compared to other thin streams, the Phoenix stream shows more fluctuations and, consequently, the study of Phoenix offers a unique perspective on gravitational perturbations of stellar streams. We discuss possible sources of perturbations to Phoenix including baryonic structures in the Galaxy and dark matter subhalos., Comment: 18 Pages, 7 Figures, Matches published version

We present multi-wavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from DES (Dark Energy Survey) between 2013-2018 and VOILETTE (VEILS Optical Light curves of Extragalactic TransienT Events) between 2018-2019 and near infrared (NIR) JKs observations from VEILS (VISTAExtragalactic Infrared Legacy Survey) between 2017-2019. Multiple methods of cross-correlation of each combination of light curve provides measurements of possible lags between optical-optical, optical-NIR, and NIR-NIR emission, for each observation season and for the entire observational period. Inter-band time lag measurements consistently suggest either simultaneous emission or delays between emission regions on timescales smaller than the cadences of observations. The colour-magnitude relation between each combination of filters was also studied to determine the spectral behaviour of PKS 0027-426. Our results demonstrate complex colour behaviour that changes between bluer when brighter (BWB), stable when brighter (SWB) and redder when brighter (RWB) trends over different timescales and using different combinations of optical filters. Additional analysis of the optical spectra is performed to provide further understanding of this complex spectral behaviour., Comment: Accepted for publication in MNRAS. 33 pages, 28 figures, 11 tables

We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $\lambda$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlations between cluster velocity dispersion, richness, X-ray temperature and luminosity as well as central galaxy velocity offsets. The velocity dispersion-richness relation exhibits a bimodal distribution. The majority of clusters follow scaling relations between velocity dispersion, richness, and X-ray properties similar to those found for previous samples; however, there is a significant population of clusters with velocity dispersions which are high for their richness. These clusters account for roughly 22\% of the $\lambda < 70$ systems in our sample, but more than half (55\%) of $\lambda < 70$ clusters at $z>0.5$. A couple of these systems are hot and X-ray bright as expected for massive clusters with richnesses that appear to have been underestimated, but most appear to have high velocity dispersions for their X-ray properties likely due to line-of-sight structure. These results suggest that projection effects contribute significantly to redMaPPer selection, particularly at higher redshifts and lower richnesses. The redMaPPer determined richnesses for the velocity dispersion outliers are consistent with their X-ray properties, but several are X-ray undetected and deeper data is needed to understand their nature., Comment: 22 pages, accepted to MNRAS

We present angular diameter measurements obtained by measuring the position of Baryon Acoustic Oscillations (BAO) in an optimised sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg$^2$ with $0.6 < z_{\rm photo} < 1.1$ and a typical redshift uncertainty of $0.03(1+z)$. The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive pre-unblinding protocol. We used two different statistics: angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance $D_M$ at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain $D_M(z_{\rm eff}=0.835)/r_{\rm d} = 18.92 \pm 0.51$, which is consistent with, but smaller than, the Planck prediction assuming flat \lcdm, at the level of $2.3 \sigma$. The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at $z>0.75$. When combined with DES 3x2pt + SNIa, they lead to improvements in $H_0$ and $\Omega_m$ constraints by $\sim 20\%$, Comment: 24 pages, 10 figures, matches version accepted by PRD

We use a sample of 809 photometrically classified type Ia supernovae (SNe Ia) discovered by the Dark Energy Survey (DES) along with 40415 field galaxies to calculate the rate of SNe Ia per galaxy in the redshift range $0.2 < z <0.6$. We recover the known correlation between SN Ia rate and galaxy stellar mass across a broad range of scales $8.5 \leq \log(M_*/\mathrm{M}_{\odot}) \leq 11.25$. We find that the SN Ia rate increases with stellar mass as a power-law with index $0.63 \pm 0.02$, which is consistent with previous work. We use an empirical model of stellar mass assembly to estimate the average star-formation histories (SFHs) of galaxies across the stellar mass range of our measurement. Combining the modelled SFHs with the SN Ia rates to estimate constraints on the SN Ia delay time distribution (DTD), we find the data are fit well by a power-law DTD with slope index $\beta = -1.13 \pm 0.05$ and normalisation $A = 2.11 \pm0.05 \times 10^{-13}~\mathrm{SNe}~{\mathrm{M}_{\odot}}^{-1}~\mathrm{yr}^{-1}$, which corresponds to an overall SN Ia production efficiency $N_{\mathrm{Ia}}/M_* = 0.9~_{-0.7}^{+4.0} \times 10^{-3}~\mathrm{SNe}~\mathrm{M}_{\odot}^{-1}$. Upon splitting the SN sample by properties of the light curves, we find a strong dependence on DTD slope with the SN decline rate, with slower-declining SNe exhibiting a steeper DTD slope. We interpret this as a result of a relationship between intrinsic luminosity and progenitor age, and explore the implications of the result in the context of SN Ia progenitors., Comment: 19 pages, 18 figures. Accepted for publication in MNRAS

The recently discovered Indus stellar stream exhibits a diverse chemical signature compared to what is found for most other streams due to the abundances of two outlier stars, Indus$\_$0 and Indus$\_$13. Indus$\_$13, exhibits an extreme enhancement in rapid neutron-capture ($r$-)process elements with $\mathrm{[Eu/Fe]} = +1.81$. It thus provides direct evidence of the accreted nature of $r$-process enhanced stars. In this paper we present a detailed chemical analysis of the neutron-capture elements in Indus$\_$13, revealing the star to be slightly actinide poor. The other outlier, Indus$\_0$, displays a globular cluster-like signature with high N, Na, and Al abundances, while the rest of the Indus stars show abundances compatible with a dwarf galaxy origin. Hence, Indus$\_0$ provides the first chemical evidence of a fully disrupted dwarf containing a globular cluster. We use the chemical signature of the Indus stars to discuss the nature of the stream progenitor which was likely a chemically evolved system, with a mass somewhere in the range from Ursa Minor to Fornax., Comment: 14 pages, 4 figures, and 5 tables, accepted for publication in ApJ

We have identified a first group of 33 new candidates for symbiotic stars (SySt) of the accreting-only variety among the 600,255 stars so far observed by the GALAH high-resolution spectroscopic survey of the Southern Hemisphere, more than doubling the number of those previously known. GALAH aims to high latitudes and this offers the possibility to sound the Galaxy for new SySt away from the usual Plane and Bulge hunting regions. In this paper we focus on SySt of the M spectral type, showing an Halpha emission with a peak in excess of 0.5 above the adjacent continuum level, and not affected by coherent radial pulsations. These constraints will be relaxed in future studies. The 33 new candidate SySt were subjected to a vast array of follow-up confirmatory observations (X-ray/UV observations with the Swift satellite, search for optical flickering, presence of a near-UV upturn in ground-based photometric and spectroscopic data, radial velocity changes suggestive of orbital motion, variability of the emission line profiles). According to Gaia eDR3 parallaxes, the new SySt are located at the tip of the Giant Branch, sharing the same distribution in M(Ks) of the well established SySt. The accretion luminosities of the new SySt are in the range 1-10 Lsun, corresponding to mass-accretion rates of 0.1-1x10(-9) Msun/yr for WDs of 1 Msun. The M giant of one of the new SySt presents a large Lithium over-abundance., Comment: MNRAS, revised version

We present DES14X2fna, a high-luminosity, fast-declining type IIb supernova (SN IIb) at redshift $z=0.0453$, detected by the Dark Energy Survey (DES). DES14X2fna is an unusual member of its class, with a light curve showing a broad, luminous peak reaching $M_r\simeq-19.3$ mag 20 days after explosion. This object does not show a linear decline tail in the light curve until $\simeq$60 days after explosion, after which it declines very rapidly (4.38$\pm$0.10 mag 100 d$^{-1}$ in $r$-band). By fitting semi-analytic models to the photometry of DES14X2fna, we find that its light curve cannot be explained by a standard $^{56}$Ni decay model as this is unable to fit the peak and fast tail decline observed. Inclusion of either interaction with surrounding circumstellar material or a rapidly-rotating neutron star (magnetar) significantly increases the quality of the model fit. We also investigate the possibility for an object similar to DES14X2fna to act as a contaminant in photometric samples of SNe Ia for cosmology, finding that a similar simulated object is misclassified by a recurrent neural network (RNN)-based photometric classifier as a SN Ia in $\sim$1.1-2.4 per cent of cases in DES, depending on the probability threshold used for a positive classification.

Reverberation mapping is a robust method to measure the masses of supermassive black holes (SMBHs) outside of the local Universe. Measurements of the radius -- luminosity ($R-L$) relation using the Mg II emission line are critical for determining these masses near the peak of quasar activity at $z \approx 1 - 2$, and for calibrating secondary mass estimators based on Mg II that can be applied to large samples with only single-epoch spectroscopy. We present the first nine Mg II lags from our five-year Australian Dark Energy Survey (OzDES) reverberation mapping program, which substantially improves the number and quality of Mg II lag measurements. As the Mg II feature is somewhat blended with iron emission, we model and subtract both the continuum and iron contamination from the multi-epoch spectra before analyzing the Mg II line. We also develop a new method of quantifying correlated spectroscopic calibration errors based on our numerous, contemporaneous observations of F-stars. The lag measurements for seven of our nine sources are consistent with both the H$\beta$ and Mg II $R-L$ relations reported by previous studies. Our simulations verify the lag reliability of our nine measurements, and we estimate that the median false positive rate of the lag measurements is $4\%$., Comment: 18 pages, 12 figures. MNRAS, Volume 507, Issue 3, November 2021, Pages 3771-3788

The Solo (Solitary Local) Dwarf Galaxy survey is a volume limited, wide-field g- and i- band survey of all known nearby (<3 Mpc) and isolated (>300 kpc from the Milky Way or M31) dwarf galaxies. This set of 44 dwarfs are homogeneously analysed for quantitative comparisons to the satellite dwarf populations of the Milky Way and M31. In this paper, an analysis of the 12 closest Solo dwarf galaxies accessible from the northern hemisphere is presented, including derivation of their distances, spatial distributions, morphology, and extended structures, including their inner integrated light properties and their outer resolved star distributions. All 12 galaxies are found to be reasonably well described by two-dimensional Sersic functions, although UGC 4879 in particular shows tentative evidence of two distinct components. No prominent extended stellar substructures, that could be signs of either faint satellites or recent mergers, are identified in the outer regions of any of the systems examined., Comment: 25 pages, 18 figures, Accepted to MNRAS

The analysis of current and future cosmological surveys of type Ia supernovae (SNe Ia) at high-redshift depends on the accurate photometric classification of the SN events detected. Generating realistic simulations of photometric SN surveys constitutes an essential step for training and testing photometric classification algorithms, and for correcting biases introduced by selection effects and contamination arising from core collapse SNe in the photometric SN Ia samples. We use published SN time-series spectrophotometric templates, rates, luminosity functions and empirical relationships between SNe and their host galaxies to construct a framework for simulating photometric SN surveys. We present this framework in the context of the Dark Energy Survey (DES) 5-year photometric SN sample, comparing our simulations of DES with the observed DES transient populations. We demonstrate excellent agreement in many distributions, including Hubble residuals, between our simulations and data. We estimate the core collapse fraction expected in the DES SN sample after selection requirements are applied and before photometric classification. After testing different modelling choices and astrophysical assumptions underlying our simulation, we find that the predicted contamination varies from 5.8 to 9.3 per cent, with an average of 7.0 per cent and r.m.s. of 1.1 per cent. Our simulations are the first to reproduce the observed photometric SN and host galaxy properties in high-redshift surveys without fine-tuning the input parameters. The simulation methods presented here will be a critical component of the cosmology analysis of the DES photometric SN Ia sample: correcting for biases arising from contamination, and evaluating the associated systematic uncertainty.

In order to accurately determine stellar properties, knowledge of the effective temperature of stars is vital. We implement Gaia and 2MASS photometry in the InfraRed Flux Method and apply it to over 360,000 stars across different evolutionary stages in the GALAH DR3 survey. We derive colour-effective temperature relations that take into account the effect of metallicity and surface gravity over the range 4000 to 8000 kelvin, from very metal-poor stars to super solar metallicities. The internal uncertainty of these calibrations is of order 40-80 kelvin depending on the colour combination used. Comparison against solar-twins, Gaia benchmark stars and the latest interferometric measurements validates the precision and accuracy of these calibrations from F to early M spectral types. We assess the impact of various sources of uncertainties, including the assumed extinction law, and provide guidelines to use our relations. Robust solar colours are also derived., Comment: MNRAS accepted. Include Gaia solar colours. Colour-effective temperature routines for Gaia DR2 and DR3 system available at https://github.com/casaluca/colte

Massive sets of stellar spectroscopic observations are rapidly becoming available and these can be used to determine the chemical composition and evolution of the Galaxy with unprecedented precision. One of the major challenges in this endeavour involves constructing realistic models of stellar spectra with which to reliably determine stellar abundances. At present, large stellar surveys commonly use simplified models that assume that the stellar atmospheres are approximately in local thermodynamic equilibrium (LTE). To test and ultimately relax this assumption, we have performed non-LTE calculations for $13$ different elements (H, Li, C, N, O, Na, Mg, Al, Si, K, Ca, Mn, and Ba), using recent model atoms that have physically-motivated descriptions for the inelastic collisions with neutral hydrogen, across a grid of $3756$ 1D MARCS model atmospheres that spans $3000\leq T_{\mathrm{eff}}/\mathrm{K}\leq8000$, $-0.5\leq\log{g/\mathrm{cm\,s^{-2}}}\leq5.5$, and $-5\leq\mathrm{[Fe/H]}\leq1$. We present the grids of departure coefficients that have been implemented into the GALAH DR3 analysis pipeline in order to complement the extant non-LTE grid for iron. We also present a detailed line-by-line re-analysis of $50126$ stars from GALAH DR3. We found that relaxing LTE can change the abundances by between $-0.7\,\mathrm{dex}$ and $+0.2\,\mathrm{dex}$ for different lines and stars. Taking departures from LTE into account can reduce the dispersion in the $\mathrm{[A/Fe]}$ versus $\mathrm{[Fe/H]}$ plane by up to $0.1\,\mathrm{dex}$, and it can remove spurious differences between the dwarfs and giants by up to $0.2\,\mathrm{dex}$. The resulting abundance slopes can thus be qualitatively different in non-LTE, possibly with important implications for the chemical evolution of our Galaxy., Comment: 19 pages, 25 figures, 2 tables, arXiv abstract abridged; accepted for publication in A&A

We consider the effects of weak gravitational lensing on observations of 196 spectroscopically confirmed Type Ia Supernovae (SNe Ia) from years 1 to 3 of the Dark Energy Survey (DES). We simultaneously measure both the angular correlation function and the non-Gaussian skewness caused by weak lensing. This approach has the advantage of being insensitive to the intrinsic dispersion of SNe Ia magnitudes. We model the amplitude of both effects as a function of $\sigma_8$, and find $\sigma_8 = 1.2^{+0.9}_{-0.8}$. We also apply our method to a subsample of 488 SNe from the Joint Light-curve Analysis (JLA) (chosen to match the redshift range we use for this work), and find $\sigma_8 = 0.8^{+1.1}_{-0.7}$. The comparable uncertainty in $\sigma_8$ between DES-SN and the larger number of SNe from JLA highlights the benefits of homogeneity of the DES-SN sample, and improvements in the calibration and data analysis., Comment: 10 pages, 7 figures, MNRAS

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., Comment: 20 pages, 16 figures. Accepted for publication by MNRAS. Data release available at https://datacentral.org.au

Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lema\^itre law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects., Comment: 39 pages, 22 figures, 10 tables, Accepted for publication in MNRAS

Rapidly evolving transients (RETs), also termed fast blue optical transients (FBOTs), are a distinct class of astrophysical event. They are characterised by lightcurves that decline much faster than common classes supernovae (SNe), span vast ranges in peak luminosity and can be seen to redshifts greater than 1. Their evolution on fast timescales has hindered high quality follow-up observations, and thus their origin and explosion/emission mechanism remains unexplained. In this paper we define the largest sample of RETs to date, comprising 106 objects from the Dark Energy Survey, and perform the most comprehensive analysis of RET host galaxies. Using deep-stacked photometry and emission-lines from OzDES spectroscopy, we derive stellar masses and star-formation rates (SFRs) for 49 host galaxies, and metallicities for 37. We find that RETs explode exclusively in star-forming galaxies and are thus likely associated with massive stars. Comparing RET hosts to samples of host galaxies of other explosive transients as well as field galaxies, we find that RETs prefer galaxies with high specific SFRs, indicating a link to young stellar populations, similar to stripped-envelope SNe. RET hosts appear to show a lack of chemical enrichment, their metallicities akin to long duration gamma-ray bursts and superluminous SN host galaxies. There are no clear relationships between properties of the host galaxies and the peak magnitudes or decline rates of the transients themselves., Comment: 16 pages, 15 figures, submitted to MNRAS

We present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, $\Omega_{\rm M}$, and the dark energy equation-of-state parameter, $w(\equiv p/\rho)$. We build a sample of 20 cosmologically useful SLSNe~I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak decline SLSN~I standardization relation with a larger dataset and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardisation via minimisation of the $\chi^2$ computed from a covariance matrix which includes statistical and systematic uncertainties. For a spatially flat $\Lambda$CDM cosmological model, we find $\Omega_{\rm M}=0.38^{+0.24}_{-0.19}$, with a rms of 0.27 mag for the residuals of the distance moduli. For a $w_0w_a$CDM cosmological model, the addition of SLSNe~I to a `baseline' measurement consisting of Planck temperature together with type Ia supernovae, results in a small improvement in the constraints of $w_0$ and $w_a$ of 4\%. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat $\Lambda$CDM model with the same precision (considering only statistical uncertainties) as current surveys that use type Ia supernovae, while providing a factor 2-3 improvement in the precision of the constraints on the time variation of dark energy, $w_0$ and $w_a$. This paper represents the proof-of-concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift ($z>1$) universe., Comment: 16 pages, 9 figures. MNRAS accepted

While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent' galaxies that host multiple SNe Ia 'siblings'. Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find 8 galaxies that hosted two likely Type Ia SNe. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at $2.8\sigma$ significance that at least 1/2 of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find $100\times$ more pairs of galaxies, and pave a new line of study on the consistency of Type Ia supernovae in the same parent galaxies. Finally, we argue that it is unlikely some of these SNe are actually single, lensed SN with multiple images., Comment: Submitted to ApJL. Comments welcome. Fermilab id included

We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$., Comment: Accepted in MNRAS. 17 pages, 6 figures

We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. Fitting spectral energy distributions to the $griz$ photometric measurements of the DES-SN host galaxies, we derive stellar masses and star-formation rates. For the DES-SN sample, when considering a 5D ($z$, $x_1$, $c$, $\alpha$, $\beta$) bias correction, we find evidence of a Hubble residual `mass step', where SNe Ia in high mass galaxies ($>10^{10} \textrm{M}_{\odot}$) are intrinsically more luminous (after correction) than their low mass counterparts by $\gamma=0.040\pm0.019$mag. This value is larger by $0.031$mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $\gamma=0.066\pm0.020$mag. The 1D-5D $\gamma$ difference for DES-SN is $0.026\pm0.009$mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the $x_1$ component of the 5D distance-bias correction. To better understand this effect, we include an intrinsic correlation between light-curve width and stellar mass in simulated SN Ia samples. We show that a 5D fit recovers $\gamma$ with $-9$mmag bias compared to a $+2$mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modeling correlations between galaxy properties and SN is necessary to determine the implications for $\gamma$ and ensure unbiased precision estimates of the dark energy equation-of-state as we enter the era of LSST., Comment: 27 pages, 13 figures; Submitted to MNRAS

The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create depth-optimised coadds using single-epoch DES-SN images that are selected based on sky and atmospheric conditions. For each of the five DES-SN seasons, a separate coadd is made from the other 4 seasons such that each SN has a corresponding deep coadd with no contaminating SN emission. The coadds reach limiting magnitudes of order $\sim 27$ in $g$-band, and have a much smaller magnitude uncertainty than the previous DES-SN host templates, particularly for faint objects. We present the resulting multi-band photometry of host galaxies for samples of spectroscopically confirmed type Ia (SNe Ia), core-collapse (CCSNe), and superluminous (SLSNe) as well as rapidly evolving transients (RETs) discovered by DES-SN. We derive host galaxy stellar masses and probabilistically compare stellar-mass distributions to samples from other surveys. We find that the DES spectroscopically confirmed sample of SNe Ia selects preferentially fewer high mass hosts at high redshift compared to other surveys, while at low redshift the distributions are consistent. DES CCSNe and SLSNe hosts are similar to other samples, while RET hosts are unlike the hosts of any other transients, although these differences have not been disentangled from selection effects., Comment: 18 pages, 17 figures. Submitted to MNRAS

We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ravi, Wambelong, Willka Yaku, and Turbio. We find evidence for a misalignment between stream tracks and the systemic proper motion of streams that may suggest a systematic gravitational influence from the Large Magellanic Cloud. These proper motions, when combined with radial velocity measurements, will allow for detailed orbit modeling which can be used to constrain properties of the LMC and its affect on nearby streams, as well as global properties of the Milky Way's gravitational potential., Comment: 22 pages, 9 figures, 6 tables; submitted to ApJ

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}., Comment: 25 pages, 14 figures (1 in appendix), 3 tables (1 in appendix). Published on MNRAS. See also paper from Shipp et al. 2019, which measures the proper motion of the DES streams

Black hole mass measurements outside the local universe are critically important to derive the growth of supermassive black holes over cosmic time, and to study the interplay between black hole growth and galaxy evolution. In this paper we present two measurements of supermassive black hole masses from reverberation mapping (RM) of the broad CIV emission line. These measurements are based on multi-year photometry and spectroscopy from the Dark Energy Survey Supernova Program (DES-SN) and the Australian Dark Energy Survey (OzDES), which together constitute the OzDES RM Program. The observed reverberation lag between the DES continuum photometry and the OzDES emission-line fluxes is measured to be $358^{+126}_{-123}$ and $343^{+58}_{-84}$ days for two quasars at redshifts of $1.905$ and $2.593$ respectively. The corresponding masses of the two supermassive black holes are $4.4 \times 10^{9}$ and $3.3 \times 10^{9}$ M$_\odot$, which are among the highest-redshift and highest-mass black holes measured to date with RM studies. We use these new measurements to better determine the CIV radius$-$luminosity relationship for high-luminosity quasars, which is fundamental to many quasar black hole mass estimates and demographic studies., Comment: 15 pages, 12 figures Updated with minor revisions to match version accepted for publication by MNRAS. Results remain the same

We present new cosmological hydrodynamic simulations that incorporate Dark Matter Annihilation Feedback (DMAF), whereby energy released from the annihilation of dark matter particles through decay channels such as photon or positron-electron pairs provide additional heating sources for local baryonic material. For annihilation rates comparable to WIMP-like particles, we find that the key influence of DMAF is to inhibit gas accretion onto halos. Such diminished gas accretion early in the lifetimes of halos results in reduced gas fractions in smaller halos, and the delayed halo formation times of larger structures, suggesting that DMAF could impact the stellar age distribution in galaxies, and morphology of dwarfs. For a dark matter particle mass of $m_\chi\sim10$~MeV, there is a `critical halo mass' of $\sim10^{13}$ M$_{\odot}$ at $z=0$, below which there are large differences when compared to $\Lambda$CDM, such as a reduction in the abundance of halo structures as large as 25 percent, reduced gas content by 50 percent and central gas densities reduced down to 10 percent within halos of mass $\sim10^{12}$ M$_{\odot}$ but with increasing effects in smaller halos. Higher dark matter particle mass models have a smaller `critical halo mass'. For a $m_\chi\sim100$~MeV model, we find differences start appearing below halo masses of $\sim10^{12}$ M$_\odot$ and a $m_\chi\gtrsim 1$~GeV model, this mass scale lies below the resolution of our simulations, though we still observe changes in the morphology of dwarf galaxies., Comment: 16 pages, 12 figures

We present a detailed analysis of the inner mass structure of the Cosmic Horseshoe (J1148+1930) strong gravitational lens system observed with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3). In addition to the spectacular Einstein ring, this systems shows a radial arc. We obtained the redshift of the radial arc counter image $z_\text{s,r} = 1.961 \pm 0.001$ from Gemini observations. To disentangle the dark and luminous matter, we consider three different profiles for the dark matter distribution: a power-law profile, the NFW, and a generalized version of the NFW profile. For the luminous matter distribution, we base it on the observed light distribution that is fitted with three components: a point mass for the central light component resembling an active galactic nucleus, and the remaining two extended light components scaled by a constant M/L. To constrain the model further, we include published velocity dispersion measurements of the lens galaxy and perform a self-consistent lensing and axisymmetric Jeans dynamical modeling. Our model fits well to the observations including the radial arc, independent of the dark matter profile. Depending on the dark matter profile, we get a dark matter fraction between 60 % and 70 %. With our composite mass model we find that the radial arc helps to constrain the inner dark matter distribution of the Cosmic Hoseshoe independently of the dark matter profile., Comment: 19 pages, 14 figures, 8 tables, submitted to A&A

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a)., Comment: 25 pages, 14 figures, 8 tables. Submitted to AJ

Measurements of the physical properties of accretion disks in active galactic nuclei are important for better understanding the growth and evolution of supermassive black holes. We present the accretion disk sizes of 22 quasars from continuum reverberation mapping with data from the Dark Energy Survey (DES) standard star fields and the supernova C fields. We construct continuum lightcurves with the \textit{griz} photometry that span five seasons of DES observations. These data sample the time variability of the quasars with a cadence as short as one day, which corresponds to a rest frame cadence that is a factor of a few higher than most previous work. We derive time lags between bands with both JAVELIN and the interpolated cross-correlation function method, and fit for accretion disk sizes using the JAVELIN Thin Disk model. These new measurements include disks around black holes with masses as small as $\sim10^7$ $M_{\odot}$, which have equivalent sizes at 2500\AA \, as small as $\sim 0.1$ light days in the rest frame. We find that most objects have accretion disk sizes consistent with the prediction of the standard thin disk model when we take disk variability into account. We have also simulated the expected yield of accretion disk measurements under various observational scenarios for the Large Synoptic Survey Telescope Deep Drilling Fields. We find that the number of disk measurements would increase significantly if the default cadence is changed from three days to two days or one day., Comment: 33 pages, 24 figures

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., Comment: 30 Pages, 18 Figures, 12 Tables. Submitted to ApJ. Comments welcome

We present an improved measurement of the Hubble constant (H_0) using the 'inverse distance ladder' method, which adds the information from 207 Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) at redshift 0.018 < z < 0.85 to existing distance measurements of 122 low redshift (z < 0.07) SNe Ia (Low-z) and measurements of Baryon Acoustic Oscillations (BAOs). Whereas traditional measurements of H_0 with SNe Ia use a distance ladder of parallax and Cepheid variable stars, the inverse distance ladder relies on absolute distance measurements from the BAOs to calibrate the intrinsic magnitude of the SNe Ia. We find H_0 = 67.8 +/- 1.3 km s-1 Mpc-1 (statistical and systematic uncertainties, 68% confidence). Our measurement makes minimal assumptions about the underlying cosmological model, and our analysis was blinded to reduce confirmation bias. We examine possible systematic uncertainties and all are below the statistical uncertainties. Our H_0 value is consistent with estimates derived from the Cosmic Microwave Background assuming a LCDM universe (Planck Collaboration et al. 2018)., Comment: 15 pages, 5 figures, updated to match accepted version

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.