ABSTRACT The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard ΛCDM model, only about 10–20 per cent of the observed temperature depression can be accounted for via its Integrated Sachs–Wolfe imprint. However, R ≳ 100 h−1Mpc supervoids elsewhere in the sky have shown ISW imprints AISW ≈ 5.2 ± 1.6 times stronger than expected from ΛCDM (AISW = 1), which warrants further inspection. Using the Year-3 redMaGiC catalogue of luminous red galaxies from the Dark Energy Survey, here we confirm the detection of the Eridanus supervoid as a significant underdensity in the Cold Spot’s direction at z < 0.2. We also show, with S/N ≳ 5 significance, that the Eridanus supervoid appears as the most prominent large-scale underdensity in the dark matter mass maps that we reconstructed from DES Year-3 gravitational lensing data. While we report no significant anomalies, an interesting aspect is that the amplitude of the lensing signal from the Eridanus supervoid at the Cold Spot centre is about 30 per cent lower than expected from similar peaks found in N-body simulations based on the standard ΛCDM model with parameters Ωm = 0.279 and σ8 = 0.82. Overall, our results confirm the causal relation between these individually rare structures in the cosmic web and in the CMB, motivating more detailed future surveys in the Cold Spot region.

Cosmological information from weak lensing surveys is maximised by dividing source galaxies into tomographic sub-samples for which the redshift distributions are estimated. Uncertainties on these redshift distributions must be correctly propagated into the cosmological results. We present hyperrank, a new method for marginalising over redshift distribution uncertainties in cosmological analyses, using discrete samples from the space of all possible redshift distributions. This is demonstrated in contrast to previous highly simplified parametric models of the redshift distribution uncertainty. In hyperrank the set of proposed redshift distributions is ranked according to a small (in this work between one and four) number of summary values, which are then sampled along with other nuisance parameters and cosmological parameters in the Monte Carlo chain used for inference. This can be regarded as a general method for marginalising over discrete realisations of data vector variation with nuisance parameters, which can consequently be sampled separately to the main parameters of interest, allowing for increased computational efficiency. We focus on the case of weak lensing cosmic shear analyses and demonstrate our method using simulations made for the Dark Energy Survey (DES). We show the method can correctly and efficiently marginalise over a range of models for the redshift distribution uncertainty. Finally, we compare hyperrank to the common mean-shifting method of marginalising over redshift uncertainty, validating that this simpler model is sufficient for use in the DES Year 3 cosmology results presented in companion papers.

We constrain the matter density $\Omega_{\mathrm{m}}$ and the amplitude of density fluctuations $\sigma_8$ within the $\Lambda$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv~\sigma_8\sqrt{\Omega_{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68% confidence limits), with a precision of 1.8%, an improvement of ~38% compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in $S_8$ is recorded. We find a mild tension of $1.5 \thinspace \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower $S_8$. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment $A_{\mathrm{IA}}$ and $S_8$, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

In this paper we present and validate the galaxy sample used for the analysis of the Baryon Acoustic Oscillation signal (BAO) in the Dark Energy Survey (DES) Y3 data. The definition is based on a colour and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.5, while ensuring a high quality photometric redshift determination. The sample covers $\approx 4100$ square degrees to a depth of $i = 22.3 \ (AB)$ at $10\sigma$. It contains 7,031,993 galaxies in the redshift range from $z$= 0.6 to 1.1, with a mean effective redshift of 0.835. Photometric redshifts are estimated with the machine learning algorithm DNF, and are validated using the VIPERS PDR2 sample. We find a mean redshift bias of $z_{\mathrm{bias}} \approx 0.01$ and a mean uncertainty, in units of $1+z$, of $\sigma_{68} \approx 0.03$. We evaluate the galaxy population of the sample, showing it is mostly built upon Elliptical to Sbc types. Furthermore, we find a low level of stellar contamination of $\lesssim 4\%$. We present the method used to mitigate the effect of spurious clustering coming from observing conditions and other large-scale systematics. We apply it to the DES Y3 BAO sample and calculate sample weights that are used to get a robust estimate of the galaxy clustering signal. This paper is one of a series dedicated to the analysis of the BAO signal in the DES Y3 data. In the companion papers, Ferrero et al. (2021) and DES Collaboration (2021), we present the galaxy mock catalogues used to calibrate the analysis and the angular diameter distance constraints obtained through the fitting to the BAO scale, respectively. The galaxy sample, masks and additional material will be released in the public DES data repository upon acceptance.

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultrafaint dwarf (UFD) galaxy Grus II. All stars exhibit a higher than expected [Mg/Ca] ratio compared to metal-poor stars in other UFD galaxies and in the Milky Way (MW) halo. Nucleosynthesis in high-mass (M) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (three) stars suggests 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 cannot 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 ther-process elements in this star matches the scaledr-process pattern of the solar system andr-process enhanced stars in other dwarf galaxies and in the MW halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites ofr-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 andr-process element enrichment of the galaxy favors a neutron star merger as the origin of ther-process elements in Grus II.

We perform a comparison of different approaches to star-galaxy classification using the broad-band photometric data from Year 1 of the Dark Energy Survey. This is done by performing a wide range of tests with and without external `truth' information, which can be ported to other similar datasets. We make a broad evaluation of the performance of the classifiers in two science cases with DES data that are most affected by this systematic effect: large-scale structure and Milky Way studies. In general, even though the default morphological classifiers used for DES Y1 cosmology studies are sufficient to maintain a low level of systematic contamination from stellar mis-classification, contamination can be reduced to the O(1%) level by using multi-epoch and infrared information from external datasets. For Milky Way studies the stellar sample can be augmented by ~20% for a given flux limit. Reference catalogs used in this work will be made available upon publication.

For decades optical time-domain searches have been tuned to find ordinary supernovae, which rise and fall in brightness over a period of weeks. Recently, supernova searches have improved their cadences and a handful of fast-evolving luminous transients (FELTs) have been identified. FELTs have peak luminosities comparable to Type Ia supernovae, but rise to maximum in <10 days and fade from view in <month. Here we present the most extreme example of this class thus far, KSN2015K, with a rise time of only 2.2 days and a time above half-maximum (t_1/2) of only 6.8 days. Here we show that, unlike Type Ia supernovae, the light curve of KSN2015K was not powered by the decay of radioactive elements. We further argue that it is unlikely that it was powered by continuing energy deposition from a central remnant (a magnetar or black hole). Using numerical radiation hydrodynamical models, we show that the light curve of KSN2015K is well fit by a model where the supernova runs into external material presumably expelled in a pre-supernova mass loss episode. The rapid rise of KSN2015K therefore probes the venting of photons when a hypersonic shock wave breaks out of a dense extended medium.

We report the discovery of a stellar overdensity 8 degrees north of the centre of the Small Magellanic Cloud (SMC; Small Magellanic Cloud Northern Over-Density; SMCNOD), using data from the first 2 yr of the Dark Energy Survey (DES) and the first year of the MAGellanic SatelLITEs Survey (MagLiteS). The SMCNOD is indistinguishable in age, metallicity and distance from the nearby SMC stars, being primarily composed of intermediate-age stars (6 Gyr, Z=0.001), with a small fraction of young stars (1 Gyr, Z=0.01). The SMCNOD has an elongated shape with an ellipticity of 0.6 and a size of similar to 6 degrees x 2 degrees. It has an absolute magnitude of M-V congruent to -7.7, r(h) = 2.1 kpc, and mu v(r < r(h)) = 31.2 mag arcsec(-2). We estimate a stellar mass of similar to 10(5) M-circle dot, following a Kroupa mass function. The SMCNOD was probably removed from the SMC disc by tidal stripping, since it is located near the head of the Magellanic Stream, and the literature indicates likely recent Large Magellanic Cloud-SMC encounters. This scenario is supported by the lack of significant H-1 gas. Other potential scenarios for the SMCNOD origin are a transient overdensity within the SMC tidal radius or a primordial SMC satellite in advanced stage of disruption.

We present a study of NGC 4993, the host galaxy of the GW170817 gravitational-wave event, the GRB 170817A short gamma-ray burst (sGRB), and the AT 2017gfo kilonova. We use Dark Energy Camera imaging, AAT spectra, and publicly available data, relating our findings to binary neutron star (BNS) formation scenarios and merger delay timescales. NGC 4993 is a nearby early-type galaxy, with an i-band Sérsic index n=4.0 and low asymmetry (A= 0.04± 0.01). These properties are unusual for sGRB hosts. However, NGC 4993 presents shell-like structures and dust lanes indicative of a recent galaxy merger, with the optical transient located close to a shell. We constrain the star formation history (SFH) of the galaxy assuming that the galaxy merger produced a star formation burst, but find little to no ongoing star formation in either spatially resolved broadband SED or spectral fitting. We use the best-fit SFH to estimate the BNS merger rate in this type of galaxy, as = ́ + R 5.7 10 yr NSM gal 3.3 0.57 6 1. If star formation is the only considered BNS formation scenario, the expected number of BNS mergers from early-type galaxies detectable with LIGO during its first two observing seasons is + 0.038 0.022 , as opposed to ∼0.5 from all galaxy types. Hypothesizing that the binary formed due to dynamical interactions during the galaxy merger, the subsequent time elapsed can constrain the delay time of the BNS coalescence. By using velocity dispersion estimates and the position of the shells, we find that the galaxy merger occurred tmer200Myr prior to the BNS coalescence.

We report the joint WASP/KELT discovery of WASP-167b/KELT-13b, a transiting hot Jupiter with a 2.02-d orbit around a $V$ = 10.5, F1V star with [Fe/H] = 0.1 $\pm$ 0.1. The 1.5 R$_{\rm Jup}$ planet was confirmed by Doppler tomography of the stellar line profiles during transit. We place a limit of $<$ 8 M$_{\rm Jup}$ on its mass. The planet is in a retrograde orbit with a sky-projected spin-orbit angle of $\lambda = -165^{\circ} \pm 5^{\circ}$. This is in agreement with the known tendency for orbits around hotter stars to be more likely to be misaligned. WASP-167/KELT-13 is one of the few systems where the stellar rotation period is less than the planetary orbital period. We find evidence of non-radial stellar pulsations in the host star, making it a $\delta$-Scuti or $\gamma$-Dor variable. The similarity to WASP-33, a previously known hot-Jupiter host with pulsations, adds to the suggestion that close-in planets might be able to excite stellar pulsations.

We present the first study of the evolution of the galaxy luminosity and stellar-mass functions (GLF and GSMF) carried out by the Dark Energy Survey (DES). We describe the COMMODORE galaxy catalogue selected from Science Verification images. This catalogue is made of ∼ 4× 10^6 galaxies at 0<z≲1.3 over a sky area of ∼155 sq. deg with i-band limiting magnitude i=23 mag. Such characteristics are unprecedented for galaxy catalogues and they enable us to study the evolution of GLF and GSMF at 0<z<1 homogeneously with the same statistically-rich data-set and free of cosmic variance effects. The aim of this study is twofold: i) we want to test our method based on the use of photometric-redshift probability density functions against literature results obtained with spectroscopic redshifts; ii) we want to shed light on the way galaxies build up their masses over cosmic time. We find that both the i-band galaxy luminosity and stellar mass functions are characterised by a double-Schechter shape at z<0.2. Both functions agree well with those based on spectroscopic redshifts. The DES GSMF agrees especially with those measured for the GAlaxy Mass Assembly and the PRism MUlti-object Survey out to z∼1. At 0.2<z<1, we find the i-band luminosity and stellar-mass densities respectively to be constant (ρ_ L∝ (1+z)^-0.12±0.11) and decreasing (ρ_ Mstar∝ (1+z)^-0.5±0.1) with z. This indicates that, while at higher redshift galaxies have less stellar mass, their luminosities do not change substantially because of their younger and brighter stellar populations. Finally, we also find evidence for a top-down mass-dependent evolution of the GSMF.

We analyzed data from the first year of a survey for Near Earth Objects (NEOs) that we are carrying out with the Dark Energy Camera (DECam) on the 4-meter Blanco telescope at the Cerro Tololo Inter-American Observatory. We implanted synthetic NEOs into the data stream to derive our nightly detection efficiency as a function of magnitude and rate of motion. Using these measured efficiencies and the Solar System absolute magnitudes derived by the Minor Planet Center for the 1377 measurements of 235 unique NEOs detected, we directly derive, for the first time from a single observational data set, the NEO size distribution from 1 km down to 10 meters. We find that there are 10^6.6 NEOs larger than 10 meters. This result implies a factor of ten fewer small NEOs than some previous results, though our derived size distribution is in good agreement with several other estimates.

We present a multi-band photometric catalog of ≈1.7 million cluster members for a field of view of ≈2° × 2° across ω Cen. Photometry is based on images collected with the Dark Energy Camera on the 4 m Blanco telescope and the Advanced Camera for Surveys on the Hubble Space Telescope. The unprecedented photometric accuracy and field coverage allowed us, for the first time, to investigate the spatial distribution of ω Cen multiple populations from the core to the tidal radius, confirming its very complex structure. We found that the frequency of blue main-sequence stars is increasing compared to red main-sequence stars starting from a distance of ≈25′ from the cluster center. Blue main-sequence stars also show a clumpy spatial distribution, with an excess in the northeast quadrant of the cluster pointing toward the direction of the Galactic center. Stars belonging to the reddest and faintest red-giant branch also show a more extended spatial distribution in the outskirts of ω Cen, a region never explored before. Both these stellar sub-populations, according to spectroscopic measurements, are more metal-rich compared to the cluster main stellar population. These findings, once confirmed, make ω Cen the only stellar system currently known where metal-rich stars have a more extended spatial distribution compared to metal-poor stars. Kinematic and chemical abundance measurements are now needed for stars in the external regions of ω Cen to better characterize the properties of these sub-populations.

Abstract The Dark Energy Survey is undertaking an observational programme imaging 1/4 of the southern hemisphere sky with unprecedented photometric accuracy. In the process of observing millions of faint stars and galaxies to constrain the parameters of the dark energy equation of state, the Dark Energy Survey will obtain pre-discovery images of the regions surrounding an estimated 100 gamma-ray bursts over 5 yr. Once gamma-ray bursts are detected by, e.g., the Swift satellite, the DES data will be extremely useful for follow-up observations by the transient astronomy community. We describe a recently-commissioned suite of software that listens continuously for automated notices of gamma-ray burst activity, collates information from archival DES data, and disseminates relevant data products back to the community in near-real-time. Of particular importance are the opportunities that non-public DES data provide for relative photometry of the optical counterparts of gamma-ray bursts, as well as for identifying key characteristics (e.g., photometric redshifts) of potential gamma-ray burst host galaxies. We provide the functional details of the DESAlert software, and its data products, and we show sample results from the application of DESAlert to numerous previously detected gamma-ray bursts, including the possible identification of several heretofore unknown gamma-ray burst hosts.

Using the science verification data of the Dark Energy Survey for a new sample of 106 X-ray selected clusters and groups, we study the stellar mass growth of bright central galaxies (BCGs) since redshift z ~ 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into an analysis of a redshift-dependent BCG-cluster mass relation, m * ∝ M 200 1.5 × 10 14 M ⊙ 0.24 ± 0.08 ( 1 + z ) − 0.19 ± 0.34 ?> , and compare the observed relation to the model prediction. We estimate the average growth rate since z = 1.0 for BCGs hosted by clusters of M200,z = 1013.8 M⊙; at z = 1.0: m*,BCG appears to have grown by 0.13 ± 0.11 dex, in tension at the ∼2.5σ significance level with the 0.40 dex growth rate expected from the semi-analytic model. We show that the build-up of extended intracluster light after z = 1.0 may alleviate this tension in BCG growth rates.

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