Phenotypic redshifts with self-organizing maps: A novel method to characterize redshift distributions of source galaxies for weak lensing

Wide-field imaging surveys such as the Dark Energy Survey (DES) rely on coarse measurements of spectral energy distributions in a few filters to estimate the redshift distribution of source galaxies. In this regime, sample variance, shot noise, and selection effects limit the attainable accuracy of redshift calibration and thus of cosmological constraints. We present a new method to combine wide-field, few-filter measurements with catalogs from deep fields with additional filters and sufficiently low photometric noise to break degeneracies in photometric redshifts. The multi-band deep field is used as an intermediary between wide-field observations and accurate redshifts, greatly reducing sample variance, shot noise, and selection effects. Our implementation of the method uses self-organizing maps to group galaxies into phenotypes based on their observed fluxes, and is tested using a mock DES catalog created from N-body simulations. It yields a typical uncertainty on the mean redshift in each of five tomographic bins for an idealized simulation of the DES Year 3 weak-lensing tomographic analysis of $\sigma_{\Delta z} = 0.007$, which is a 60% improvement compared to the Year 1 analysis. Although the implementation of the method is tailored to DES, its formalism can be applied to other large photometric surveys with a similar observing strategy.Comment: 24 pages, 11 figures; matches version accepted to MNRA

Testing the isotropy of the dark energy Survey's extreme trans-neptunian objects

We test whether the population of "extreme"trans-Neptunian objects (eTNOs) detected in the first four years of the Dark Energy Survey (DES Y4) data exhibit azimuthal asymmetries that 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 that, 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 Ω, and argument of perihelion ω. We then compare the detected distributions in each of Ω, ω, and the longitude of perihelion {equation presented} to those expected from the isotropic population, using Kuiper's variant of the Kolmogorov-Smirnov test. The three angles are tested for each of four 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 12 total tests, two have the likelihood of drawing the observed angles from the isotropic population at p 250 and q > 37 au and the four detections at a > 250 and q > 30 au have a Ω distribution with p ≈ 0.03 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

Testing the isotropy of the Dark Energy Survey's extreme trans-Neptunian objects

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 $p250, 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.Comment: Accepted on PS

On the relative bias of void tracers in the Dark Energy Survey

Luminous tracers of large-scale structure are not entirely representative of the distribution of mass in our Universe. As they arise from the highest peaks in the matter density field, the spatial distribution of luminous objects is biased towards those peaks. On large scales, where density fluctuations are mild, this bias simply amounts to a constant offset in the clustering amplitude of the tracer, known as linear bias. In this work we focus on the relative bias between galaxies and galaxy clusters that are located inside and in the vicinity of cosmic voids, extended regions of relatively low density in the large-scale structure of the Universe. With the help of mock data we verify that the relation between galaxy and cluster overdensity around voids remains linear. Hence, the void-centric density profiles of different tracers can be linked by a single multiplicative constant. This amounts to the same value as the relative linear bias between tracers for the largest voids in the sample. For voids of small sizes, which typically arise in higher density regions, this constant has a higher value, possibly showing an environmental dependence similar to that observed for the linear bias itself. We confirm our findings by analysing data obtained during the first year of observations by the Dark Energy Survey. As a side product, we present the first catalogue of three-dimensional voids extracted from a photometric survey with a controlled photo-z uncertainty. Our results will be relevant in forthcoming analyses that attempt to use voids as cosmological probes

Observation and Confirmation of Nine Strong Lensing Systems in Dark Energy Survey Year 1 Data

We describe the observation and confirmation of nine new strong gravitational lenses discovered in Year 1 data from the Dark Energy Survey (DES). We created candidate lists based on a) galaxy group and cluster samples and b) photometrically selected galaxy samples. We selected 46 candidates through visual inspection and then used the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope to acquire spectroscopic follow-up of 21 of these candidates. Through analysis of this spectroscopic follow-up data, we confirmed nine new lensing systems and rejected two candidates, but the analysis was inconclusive on 10 candidates. For each of the confirmed systems, we report measured spectroscopic properties, estimated source image-lens separation, and estimated enclosed masses. The sources that we targeted have an i-band surface brightness range of iSB ∼ 22 − 24 mag/arcsec2 and a spectroscopic redshift range of zspec ∼ 0.8 − 2.6. The lens galaxies have a photometric redshift range of zlens ∼ 0.3 − 0.7. The lensing systems range in source image-lens separation 2 − 9″ and in enclosed mass 1012 − 1013M⊙

Core or Cusps: The Central Dark Matter Profile of a Strong Lensing Cluster with a Bright Central Image at Redshift 1

We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at z = 1.06. The arc system is notable for the presence of a bright central image. The source is a Lyman break galaxy at z s = 2.39 and the mass enclosed within the Einstein ring of radius 14 arcsec is $\sim {10}^{14.2}\ {M}_{\odot }$. We perform a full reconstruction of the light profile of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro–Frenk–White profile—with a free parameter for the inner density slope—we find that the break radius is ${270}_{-76}^{+48}$ kpc, and that the inner density falls with radius to the power −0.38 ± 0.04 at 68% confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter-only simulations predict that the inner density should fall as ${r}^{-1}$. The tension can be alleviated if this cluster is in fact a merger; a two-halo model can also reconstruct the data, with both clumps (density varying as ${r}^{-0.8}$ and ${r}^{-1.0}$) much more consistent with predictions from dark matter-only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them

Synchronous Rotation in the (136199) Eris–Dysnomia System

We combine photometry of Eris from a 6 month campaign on the Palomar 60 inch telescope in 2015, a 1 month Hubble Space Telescope WFC3 campaign in 2018, and Dark Energy Survey data spanning 2013–2018 to determine a light curve of definitive period 15.771 ± 0.008 days (1σ formal uncertainties), with nearly sinusoidal shape and peak-to-peak flux variation of 3%. This is consistent at part-per-thousand precision with the P = 15.785 90 ± 0.00005 day sidereal period of Dysnomia's orbit around Eris, strengthening the recent detection of synchronous rotation of Eris by Szakáts et al. with independent data. Photometry from Gaia are consistent with the same light curve. We detect a slope of 0.05 ± 0.01 mag per degree of Eris's brightness with respect to illumination phase averaged across g, V, and r bands, intermediate between Pluto's and Charon's values. Variations of 0.3 mag are detected in Dysnomia's brightness, plausibly consistent with a double-peaked light curve at the synchronous period. The synchronous rotation of Eris is consistent with simple tidal models initiated with a giant-impact origin of the binary, but is difficult to reconcile with gravitational capture of Dysnomia by Eris. The high albedo contrast between Eris and Dysnomia remains unexplained in the giant-impact scenario

Dark Energy Survey Year 1 results: validation of weak lensing cluster member contamination estimates from P(z) decomposition

Weak lensing source galaxy catalogues used in estimating the masses of galaxy clusters can be heavily contaminated by cluster members, prohibiting accurate mass calibration. In this study, we test the performance of an estimator for the extent of cluster member contamination based on decomposing the photometric redshift P(z) of source galaxies into contaminating and background components. We perform a full scale mock analysis on a simulated sky survey approximately mirroring the observational properties of the Dark Energy Survey Year One observations (DES Y1), and find excellent agreement between the true number profile of contaminating cluster member galaxies in the simulation and the estimated one. We further apply the method to estimate the cluster member contamination for the DES Y1 redMaPPer cluster mass calibration analysis, and compare the results to an alternative approach based on the angular correlation of weak lensing source galaxies. We find indications that the correlation based estimates are biased by the selection of the weak lensing sources in the cluster vicinity, which does not strongly impact the P(z) decomposition method. Collectively, these benchmarks demonstrate the strength of the P(z) decomposition method in alleviating membership contamination and enabling highly accurate cluster weak lensing studies without broad exclusion of source galaxies, thereby improving the total constraining power of cluster mass calibration via weak lensing

Dark Energy Survey Year 1 results: validation of weak lensing cluster member contamination estimates from P(z) decomposition

Weak lensing source galaxy catalogues used in estimating the masses of galaxy clusters can be heavily contaminated by cluster members, prohibiting accurate mass calibration. In this study, we test the performance of an estimator for the extent of cluster member contamination based on decomposing the photometric redshift P(z) of source galaxies into contaminating and background components. We perform a full scale mock analysis on a simulated sky survey approximately mirroring the observational properties of the Dark Energy Survey Year One observations (DES Y1), and find excellent agreement between the true number profile of contaminating cluster member galaxies in the simulation and the estimated one. We further apply the method to estimate the cluster member contamination for the DES Y1 redMaPPer cluster mass calibration analysis, and compare the results to an alternative approach based on the angular correlation of weak lensing source galaxies. We find indications that the correlation based estimates are biased by the selection of the weak lensing sources in the cluster vicinity, which does not strongly impact the P(z) decomposition method. Collectively, these benchmarks demonstrate the strength of the P(z) decomposition method in alleviating membership contamination and enabling highly accurate cluster weak lensing studies without broad exclusion of source galaxies, thereby improving the total constraining power of cluster mass calibration via weak lensing

Characterizing the intracluster light over the redshift range 0.2 < z < 0.8 in the DES-ACT overlap

We characterize the properties and evolution of bright central galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT) overlapping regions, covering the redshift range 0.20 14.4. We also measure the stellar mass–halo mass (SMHM) relation for the BCG+ICL system and find that the slope, β, which characterizes the dependence of M200m,SZ on the BCG+ICL stellar mass, increases with radius. The outskirts are more strongly correlated with the halo than the core, which supports that the BCG+ICL system follows a two-phase growth, where recent growth (z < 2) occurs beyond the BCG’s core. Additionally, we compare our observed SMHM relation results to the IllustrisTNG300-1 cosmological hydrodynamic simulations and find moderate qualitative agreement in the amount of diffuse light. However, the SMHM relation’s slope is steeper in TNG300-1 and the intrinsic scatter is lower, likely from the absence of projection effects in TNG300-1. Additionally, we find that the ICL exhibits a colour gradient such that the outskirts are bluer than the core. Moreover, for the lower halo mass clusters (log10(M200m,SZ/M⊙) < 14.59), we detect a modest change in the colour gradient’s slope with lookback time, which combined with the absence of stellar mass growth may suggest that lower mass clusters have been involved in growth via tidal stripping more recently than their higher mass counterparts