The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $\mu$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left<z\right>\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$\mu$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work.

We present the results from a search for the electromagnetic counterpart of the LIGO/Virgo event S190510g using the Dark Energy Camera (DECam). S190510g is a binary neutron star (BNS) merger candidate of moderate significance detected at a distance of 227$\pm$92 Mpc and localized within an area of 31 (1166) square degrees at 50\% (90\%) confidence. While this event was later classified as likely non-astrophysical in nature within 30 hours of the event, our short latency search and discovery pipeline identified 11 counterpart candidates, all of which appear consistent with supernovae following offline analysis and spectroscopy by other instruments. Later reprocessing of the images enabled the recovery of 6 more candidates. Additionally, we implement our candidate selection procedure on simulated kilonovae and supernovae under DECam observing conditions (e.g., seeing, exposure time) with the intent of quantifying our search efficiency and making informed decisions on observing strategy for future similar events. This is the first BNS counterpart search to employ a comprehensive simulation-based efficiency study. We find that using the current follow-up strategy, there would need to be 19 events similar to S190510g for us to have a 99\% chance of detecting an optical counterpart, assuming a GW170817-like kilonova. We further conclude that optimization of observing plans, which should include preference for deeper images over multiple color information, could result in up to a factor of 1.5 reduction in the total number of followups needed for discovery.

We present a measurement of the Hubble constant $H_0$ using the gravitational wave (GW) event GW190814, which resulted from the coalescence of a 23 $M_\odot$ black hole with a 2.6 $M_\odot$ compact object, as a standard siren. No compelling electromagnetic counterpart with associated host galaxy has been identified for this event, thus our analysis accounts for $\sim$ 2,700 potential host galaxies within a statistical framework. The redshift information is obtained from the photometric redshift (photo-$z$) catalog from the Dark Energy Survey. The luminosity distance is provided by the gravitational wave sky map published by the LIGO/Virgo Collaboration. Since this GW event has the second-smallest sky localization area after GW170817, GW190814 is likely to provide the best constraint on cosmology from a single standard siren without identifying an electromagnetic counterpart. Our analysis uses photo-$z$ probability distribution functions and corrects for photo-$z$ biases. We also reanalyze the binary-black hole GW170814 within this updated framework. We explore how our findings impact the $H_0$ constraints from GW170817, the only GW merger associated with a unique host galaxy, and therefore the most powerful standard siren to date. From a combination of GW190814, GW170814 and GW170817, our analysis yields $H_0 = 69.0^{+ 14}_{- 7.5 }~{\rm km~s^{-1}~Mpc^{-1}}$ (68% Highest Density Interval, HDI) for a prior in $H_0$ uniform between $[20,140]~{\rm km~s^{-1}~Mpc^{-1}}$. The addition of GW190814 and GW170814 to GW170817 improves the 68% HDI from GW170817 alone by $\sim 12\%$, showing how well-localized mergers without counterparts can provide a marginal contribution to standard siren measurements, provided that a complete galaxy catalog is available at the location of the event.

We present the weak lensing mass calibration of the stellar mass based $\mu_{\star}$ mass proxy for redMaPPer galaxy clusters in the Dark Energy Survey Year 1. For the first time we are able to perform a calibration of $\mu_{\star}$ at high redshifts, $z>0.33$. In a blinded analysis, we use $\sim 6,000$ clusters split into 12 subsets spanning the ranges $0.1 \leqslant z<0.65$ and $\mu_{\star}$ up to $\sim 5.5 \times 10^{13} M_{\odot}$, and infer the average masses of these subsets through modelling of their stacked weak lensing signal. In our model we account for the following sources of systematic uncertainty: shear measurement and photometric redshift errors, miscentring, cluster-member contamination of the source sample, deviations from the NFW halo profile, halo triaxiality and projection effects. We use the inferred masses to estimate the joint mass--$\mu_{\star}$--$z$ scaling relation given by $\langle M_{200c} | \mu_{\star},z \rangle = M_0 (\mu_{\star}/5.16\times 10^{12} \mathrm{M_{\odot}})^{F_{\mu_{\star}}} ((1+z)/1.35)^{G_z}$. We find $M_0= (1.14 \pm 0.07) \times 10^{14} \mathrm{M_{\odot}}$ with $F_{\mu_{\star}}= 0.76 \pm 0.06$ and $G_z= -1.14 \pm 0.37$. We discuss the use of $\mu_{\star}$ as a complementary mass proxy to the well-studied richness $\lambda$ for: $i)$ exploring the regimes of low $z$, $\lambda<20$ and high $\lambda$, $z \sim 1$; $ii)$ testing systematics such as projection effects for applications in cluster cosmology.

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$\beta$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{{\rm 2-12\,keV}}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$\beta$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$\sigma$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

On 14 August 2019, the LIGO and Virgo Collaborations alerted the astronomical community of a high significance detection of gravitational waves and classified the source as a neutron star - black hole (NSBH) merger, the first event of its kind. In search of an optical counterpart, the Dark Energy Survey (DES) Gravitational Wave Search and Discovery Team performed the most thorough and accurate analysis to date, targeting the entire 90 percent confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger was detected. Objects with varying brightness were detected by the DES Search and Discovery Pipeline and we systematically reduced the list of candidate counterparts through catalog matching, light curve properties, host-galaxy photometric redshifts, SOAR spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to simulations of supernovae and kilonovae as they would appear in the DECam observations. Since there are no explicit light curve models for NSBH mergers, we characterize our sensitivity with binary NS models that are expected to have similar optical signatures as NSBH mergers. We find that if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than $10^{-8.56}$, and has a velocity between $0.18c$ and $0.21c$ are disfavored at the $2\sigma$ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the $4\sigma$ level in $95\%$ of future follow-up observations.

We present a catalog of 20,977 extended low-surface-brightness galaxies (LSBGs) identified in ~ 5000 deg$^2$ from the first three years of imaging data from the Dark Energy Survey (DES). Based on a single-component S\'ersic model fit, we define extended LSBGs as galaxies with $g$-band effective radii $R_{\scriptsize{eff}} > 2.5"$ and mean surface brightness $\bar{\mu}_{\scriptsize{eff}}(g) > 24.3$ mag arcsec$^{-2}$. We find that the distribution of LSBGs is strongly bimodal in $(g-r)$ vs. $(g-i$) color space. We divide our sample into red ($g-i \geq 0.59$) and blue ($g-i<0.59$) galaxies and study the properties of the two populations. Redder LSBGs are more clustered than their blue counterparts, and are correlated with the distribution of nearby ($z < 0.10$) bright galaxies. Red LSBGs constitute $\sim 35\%$ of our LSBG sample, and $\sim 30\%$ of these are located within 1 deg of low-redshift galaxy groups and clusters (compared to $\sim 8\%$ of the blue LSBGs). For nine of the most prominent galaxy groups and clusters, we calculate the physical properties of associated LSBGs assuming a redshift derived from the host system. In these systems, we identify 108 objects that can be classified as ultra-diffuse galaxies, defined as LSBGs with projected physical effective radii $R_{\scriptsize{eff}} > 1.5$ kpc. The wide-area sample of LSBGs in DES can be used to test the role of environment on models of LSBG formation and evolution.

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 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}$.

We use a stacking method to study the radial light profiles of luminous red galaxies (LRGs) at redshift $\sim 0.62$ and $\sim 0.25$, out to a radial range of 200 kpc. We do not find noticeable evolution of the profiles at the two redshifts. The LRG profiles appear to be well approximated by a single Sersic profile, although some excess light can be seen outside 60 kpc. We quantify the excess light by measuring the integrated flux and find that the excess is about 10\% -- a non-dominant but still nonnegligible component.

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.

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.

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.

We use imaging from the first three years of the Dark Energy Survey to characterize the dynamical state of 288 galaxy clusters at $0.1 \lesssim z \lesssim 0.9$ detected in the South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) effect survey (SPT-SZ). We examine spatial offsets between the position of the brightest cluster galaxy (BCG) and the center of the gas distribution as traced by the SPT-SZ centroid and by the X-ray centroid/peak position from Chandra and XMM data. We show that the radial distribution of offsets provides no evidence that SPT SZ-selected cluster samples include a higher fraction of mergers than X-ray-selected cluster samples. We use the offsets to classify the dynamical state of the clusters, selecting the 43 most disturbed clusters, with half of those at $z \gtrsim 0.5$, a region seldom explored previously. We find that Schechter function fits to the galaxy population in disturbed clusters and relaxed clusters differ at $z>0.55$ but not at lower redshifts. Disturbed clusters at $z>0.55$ have steeper faint-end slopes and brighter characteristic magnitudes. Within the same redshift range, we find that the BCGs in relaxed clusters tend to be brighter than the BCGs in disturbed samples, while in agreement in the lower redshift bin. Possible explanations includes a higher merger rate, and a more efficient dynamical friction at high redshift. The red-sequence population is less affected by the cluster dynamical state than the general galaxy population.

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.

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.

For ground-based optical imaging with current CCD technology, the Poisson fluctuations in source and sky background photon arrivals dominate the noise budget and are readily estimated. Another component of noise, however, is the signal from the undetected population of stars and galaxies. Using injection of artificial galaxies into images, we demonstrate that the measured variance of galaxy moments (used for weak gravitational lensing measurements) in Dark Energy Survey (DES) images is significantly in excess of the Poisson predictions, by up to 30\%, and that the background sky levels are overestimated by current software. By cross-correlating distinct images of "empty" sky regions, we establish that there is a significant image noise contribution from undetected static sources (US), which on average are mildly resolved at DES resolution. Treating these US as a stationary noise source, we compute a correction to the moment covariance matrix expected from Poisson noise. The corrected covariance matrix matches the moment variances measured on the injected DES images to within 5\%. Thus we have an empirical method to statistically account for US in weak lensing measurements, rather than requiring extremely deep sky simulations. We also find that local sky determinations can remove the bias in flux measurements, at a small penalty in additional, but quantifiable, noise.

In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $\kappa_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $\kappa_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $\kappa_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed.

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

We test whether the population of "extreme" trans-Neptunian objects (eTNOs) detected in the Y4 Dark Energy Survey (DES) data exhibit azimuthal asymmetries which might be evidence of gravitational perturbations from an unseen super-Earth in a distant orbit. By rotating the orbits of the detected eTNOs, we construct a synthetic population which, when subject to the DES selection function, reproduces the detected distribution of eTNOs in the orbital elements $a,e,$ and $i$ as well as absolute magnitude $H$, but has uniform distributions in mean anomaly $M$, longitude of ascending node $\Omega,$ and argument of perihelion $\omega.$ We then compare the detected distributions in each of $\Omega, \omega,$ and $\varpi\equiv\Omega+\omega$ to those expected from the isotropic population, using Kuiper's variant of the Kolmogorov-Smirnov test. The three angles are tested for each of 4 definitions of the eTNO population, choosing among $a>(150,250)$~AU and perihelion $q>(30,37)$~AU. These choices yield 3--7 eTNOs in the DES Y4 sample. Among the twelve total tests, two have the likelihood of drawing the observed angles from the isotropic population at $p<0.05.$ The 3 detections at $a>250, q>37$~AU, and the 4 detections at $a>250, q>30$~AU, have $\Omega$ distribution with $p=0.03$ of coming from the isotropic construction, but this is not strong evidence of anisotropy given the 12 different tests. The DES data taken on their own are thus consistent with azimuthal isotropy and do not require a "Planet 9" hypothesis. The limited sky coverage and object count mean, however, that the DES data by no means falsify this hypothesis.

We construct and validate the selection function of the MARD-Y3 sample. This sample was selected through optical follow-up of the 2nd ROSAT faint source catalog (2RXS) with Dark Energy Survey year 3 (DES-Y3) data. The selection function is modeled by combining an empirically constructed X-ray selection function with an incompleteness model for the optical follow-up. We validate the joint selection function by testing the consistency of the constraints on the X-ray flux--mass and richness--mass scaling relation parameters derived from different sources of mass information: (1) cross-calibration using SPT-SZ clusters, (2) calibration using number counts in X-ray, in optical and in both X-ray and optical while marginalizing over cosmological parameters, and (3) other published analyses. We find that the constraints on the scaling relation from the number counts and SPT-SZ cross-calibration agree, indicating that our modeling of the selection function is adequate. Furthermore, we apply a largely cosmology independent method to validate selection functions via the computation of the probability of finding each cluster in the SPT-SZ sample in the MARD-Y3 sample and vice-versa. This test reveals no clear evidence for MARD-Y3 contamination, SPT-SZ incompleteness or outlier fraction. Finally, we discuss the prospects of the techniques presented here to limit systematic selection effects in future cluster cosmological studies.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

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.

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.

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.

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.

The statistical properties of galaxy clusters can only be used for cosmological purposes if observational effects related to cluster detection are accurately characterized. These effects include the selection function associated to cluster finder algorithms and survey strategy. The importance of the selection becomes apparent when different cluster finders are applied to the same galaxy catalog, producing different cluster samples. We consider parametrized functional forms for the observable-mass relation, its scatter as well as the completeness and purity of cluster samples, and study how prior knowledge on these function parameters affects dark energy constraints derived from cluster statistics. Under the assumption that completeness and purity reach 50 % at masses around 10^{13.5} Msun/h, we find that self-calibration of selection parameters in current and upcoming cluster surveys is possible, while still allowing for competitive dark energy constraints. We consider a fiducial survey with specifications similar to those of the Dark Energy Survey (DES) with 5000 deg^2, maximum redshift of zmax ~ 1.0 and threshold observed mass M_{th} ~ 10^{13.8} Msun/h, such that completeness and purity ~ 60 % - 80 % at masses around M_{th}. Perfect knowledge of all selection parameters allows for constraining a constant dark energy equation of state to sigma(w)=0.033. Employing a joint fit including self-calibration of the effective selection degrades constraints to sigma(w)=0.046. External calibrations at the level of 1 % in the parameters of the observable-mass relation and completeness/purity functions are necessary to improve the joint constraints to sigma(w)=0.041. In the lack of knowledge of selection parameters, future experiments probing larger areas and greater depths suffer from stronger relative degradations on dark energy constraints compared to current surveys.

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.

The cosmological utility of galaxy cluster catalogues is primarily limited by our ability to calibrate the relation between halo mass and observable mass proxies such as cluster richness, X-ray luminosity or the Sunyaev-Zeldovich signal. Projection effects are a particularly pernicious systematic effect that can impact observable mass proxies; structure along the line of sight can both bias and increase the scatter of the observable mass proxies used in cluster abundance studies. In this work, we develop an empirical method to characterize the impact of projection effects on redMaPPer cluster catalogues. We use numerical simulations to validate our method and illustrate its robustness. We demonstrate that modeling of projection effects is a necessary component for cluster abundance studies capable of reaching $\approx 5\%$ mass calibration uncertainties (e.g. the Dark Energy Survey Year 1 sample). Specifically, ignoring the impact of projection effects in the observable--mass relation --- i.e. marginalizing over a log-normal model only --- biases the posterior of the cluster normalization condition $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{1/2}$ by $\Delta S_8 =0.05$, more than twice the uncertainty in the posterior for such an analysis.

Super-sample covariance (SSC) is the dominant source of statistical error on large scale structure (LSS) observables for both current and future galaxy surveys. In this work, we concentrate on the SSC of cluster counts, also known as sample variance, which is particularly useful for the self-calibration of the cluster observable-mass relation; our approach can similarly be applied to other observables, such as galaxy clustering and lensing shear. We first examined the accuracy of two analytical approximations proposed in the literature for the flat sky limit, finding that they are accurate at the 15% and 30-35% level, respectively, for covariances of counts in the same redshift bin. We then developed a harmonic expansion formalism that allows for the prediction of SSC in an arbitrary survey mask geometry, such as large sky areas of current and future surveys. We show analytically and numerically that this formalism recovers the full sky and flat sky limits present in the literature. We then present an efficient numerical implementation of the formalism, which allows fast and easy runs of covariance predictions when the survey mask is modified. We applied our method to a mask that is broadly similar to the Dark Energy Survey footprint, finding a non-negligible negative cross-z covariance, i.e. redshift bins are anti-correlated. We also examined the case of data removal from holes due to, for example bright stars, quality cuts, or systematic removals, and find that this does not have noticeable effects on the structure of the SSC matrix, only rescaling its amplitude by the effective survey area. These advances enable analytical covariances of LSS observables to be computed for current and future galaxy surveys, which cover large areas of the sky where the flat sky approximation fails.

The existence of multiple subclasses of type Ia supernovae (SNeIa) has been the subject of great debate in the last decade. One major challenge inevitably met when trying to infer the existence of one or more subclasses is the time consuming, and subjective, process of subclass definition. In this work, we show how machine learning tools facilitate identification of subtypes of SNeIa through the establishment of a hierarchical group structure in the continuous space of spectral diversity formed by these objects. Using Deep Learning, we were capable of performing such identification in a 4 dimensional feature space (+1 for time evolution), while the standard Principal Component Analysis barely achieves similar results using 15 principal components. This is evidence that the progenitor system and the explosion mechanism can be described by a small number of initial physical parameters. As a proof of concept, we show that our results are in close agreement with a previously suggested classification scheme and that our proposed method can grasp the main spectral features behind the definition of such subtypes. This allows the confirmation of the velocity of lines as a first order effect in the determination of SNIa subtypes, followed by 91bg-like events. Given the expected data deluge in the forthcoming years, our proposed approach is essential to allow a quick and statistically coherent identification of SNeIa subtypes (and outliers). All tools used in this work were made publicly available in the Python package Dimensionality Reduction And Clustering for Unsupervised Learning in Astronomy (DRACULA) and can be found within COINtoolbox (https://github.com/COINtoolbox/DRACULA).