Monte Carlo control loops for cosmic shear cosmology with DES Year 1 data

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum C ℓ and the redshift distribution n ( z ) of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of C ℓ and n ( z ) curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of Ω m , σ 8 , and intrinsic alignment amplitude A IA , defined as S 8 D IA = σ 8 ( Ω m / 0.3 ) 0.5 D IA , where D IA = 1 − 0.11 ( A IA − 1 ) . We find S 8 D IA = 0.89 5 + 0.054 − 0.039 , where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys

Constraining radio mode feedback in galaxy clusters with the cluster radio AGNs properties to z ∼ 1

We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843 MHz radio active galactic nuclei (AGNs) population in galaxy clusters from two large catalogues created using the Dark Energy Survey (DES): ∼11 800 optically selected RM-Y3 and ∼1000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGNs are highly concentrated around cluster centres to z ∼ 1. We measure the halo occupation number for cluster radio AGNs above a threshold luminosity, finding that the number of radio AGNs per cluster increases with cluster halo mass as N ∝ M1.2 ± 0.1 (N ∝ M0.68 ± 0.34) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad-band colours and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGNs luminosity function (LF). The LF exhibits a redshift trend of the form (1 + z⁠)γ in density and luminosity, respectively, of γD = 3.0 ± 0.4 and γP = 0.21 ± 0.15 in the RM-Y3 sample, and γD = 2.6 ± 0.7 and γP = 0.31 ± 0.15 in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGNs, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core (<0.1R500) X-ray radiative losses for clusters at z < 1

Multiwavelength optical and NIR variability analysis of the Blazar PKS 0027-426

We present multiwavelength spectral and temporal variability analysis of PKS 0027-426 using optical griz observations from Dark Energy Survey between 2013 and 2018 and VEILS Optical Light curves of Extragalactic TransienT Events (VOILETTE) between 2018 and 2019 and near-infrared (NIR) JKs observations from Visible and Infrared Survey Telescope for Astronomy Extragalactic Infrared Legacy Survey (VEILS) between 2017 and 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 time-scales 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, stable when brighter, and redder when brighter trends over different time-scales and using different combinations of optical filters. Additional analysis of the optical spectra is performed to provide further understanding of this complex spectral behaviour

A DECam Search for Explosive Optical Transients Associated with IceCube Neutrino Alerts

In this work, we investigate the likelihood of association between real-time, neutrino alerts with teraelectronvolt to petaelectronvolt energy from IceCube and optical counterparts in the form of core-collapse supernovae (CC SNe). The optical follow-up of IceCube alerts requires two main instrumental capabilities: (1) deep imaging, since 73% of neutrinos would come from CC SNe at redshifts z > 0.3, and (2) a large field of view (FoV), since typical IceCube muon neutrino pointing accuracy is on the order of ~1 deg. With Blanco/DECam (gri to 24th magnitude and 2.2 deg diameter FoV), we performed a triggered optical follow-up observation of two IceCube alerts, IC170922A and IC171106A, on six nights during the three weeks following each alert. For the IC170922A (IC171106A) follow-up observations, we expect that 12.1% (9.5%) of coincident CC SNe at z lesssim 0.3 are detectable, and that, on average, 0.23 (0.07) unassociated SNe in the neutrino 90% containment regions also pass our selection criteria. We find two candidate CC SNe that are temporally coincident with the neutrino alerts in the FoV, but none in the 90% containment regions, a result that is statistically consistent with expected rates of background CC SNe for these observations. If CC SNe are the dominant source of teraelectronvolt to petaelectronvolt neutrinos, we would expect an excess of coincident CC SNe to be detectable at the 3σ confidence level using DECam observations similar to those of this work for ~60 (~200) neutrino alerts with (without) redshift information for all candidates

OzDES Reverberation Mapping Program: H beta lags from the 6-yr survey

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 to use to standardize AGN to determine cosmological distances. We present reverberation measurements made with Hβ from the 6-yr Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at 0.12 &amp;lt; z &amp;lt; 0.71, probing higher redshifts than the bulk of Hβ measurements made to date. Our fit to the R–L relation has a slope of α = 0.41 ± 0.03 and an intrinsic scatter of σ = 0.23 ± 0.02 dex. The results from our multi-object spectroscopic survey are consistent with previous measurements made by dedicated source-by-source campaigns, and with the observed dependence on accretion rate. Future surveys, including LSST, TiDES, and SDSS-V, which will be revisiting some of our observed fields, will be able to build on the results of our first-generation multi-object reverberation mapping survey

Candidate periodically variable quasars from the Dark Energy Survey and the Sloan Digital Sky Survey

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-yr long multicolour (griz) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r ∼23.5 mag) powered by less massive black holes (with masses ≳ 108.5M⊙) at high redshift for the first time. We find five candidates with significant (at >99.74 per cent single-frequency significance in at least two bands with a global p-value of ∼7 × 10−4–3 × 10−3 accounting for the look-elsewhere effect) periodicity with observed periods of ∼3–5 yr (i.e. 1–2 yr in rest frame) having ∼4–6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ∼0.8+0.5−0.3 per cent or ∼1.1+0.7−0.5 quasar per deg2. Our detection rate is 4–80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins

The Observed Evolution of the Stellar Mass-Halo Mass Relation for Brightest Central Galaxies

We quantify evolution in the cluster-scale stellar mass–halo mass (SMHM) relation's parameters using 2323 clusters and brightest central galaxies (BCGs) over the redshift range 0.03 ≤ z ≤ 0.60. The precision on the inferred SMHM parameters is improved by including the magnitude gap (mgap) between the BCG and fourth-brightest cluster member (M14) as a third parameter in the SMHM relation. At fixed halo mass, accounting for mgap, through a stretch parameter, reduces the SMHM relation's intrinsic scatter. To explore this redshift range, we use clusters, BCGs, and cluster members identified using the Sloan Digital Sky Survey C4 and redMaPPer cluster catalogs and the Dark Energy Survey redMaPPer catalog. Through this joint analysis, we detect no systematic differences in BCG stellar mass, mgap, and cluster mass (inferred from richness) between the data sets. We utilize the Pareto function to quantify each parameter's evolution. We confirm prior findings of negative evolution in the SMHM relation's slope (3.5σ), and detect negative evolution in the stretch parameter (4.0σ) and positive evolution in the offset parameter (5.8σ). This observed evolution, combined with the absence of BCG growth, when stellar mass is measured within 50 kpc, suggests that this evolution results from changes in the cluster's mgap. For this to occur, late-term growth must be in the intracluster light surrounding the BCG. We also compare the observed results to IllustrisTNG 300-1 cosmological hydrodynamic simulations and find modest qualitative agreement. However, the simulations lack the evolutionary features detected in the real data

Search for RR lyrae stars in DES ultrafaint systems: grus I, kim 2, phoenix II, and grus

FINEP - FINANCIADORA DE ESTUDOS E PROJETOSFAPERJ - FUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DO RIO DE JANEIROCNPQ - CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICOMCTIC - MINISTÉRIO DA CIÊNCIA, TECNOLOGIA, INOVAÇÕES E COMUNICAÇÕESThis work presents the first search for RR Lyrae stars (RRLs) in four he ultrafaint systems imaged by the Dark Energy Survey using SOAR/Goodman and Blanco/DECam imagers, We have detected two RRLs in the field of Grus I, none in Kim 2, one in Phoenix II, and four in Grus II. With the detection of these stars, we accurately determine the distance moduli for these ultrafaint dwarf satellite galaxies; mu(0) = 20.51 +/- 0.10 mag (D-circle dot = 127 +/- 6 kpc) for Grus I and mu(0) = 20.01 +/- 0.10 mag (D-circle dot = 100 +/- 5 kpc) for Phoenix IL These measurements are larger than previous estimations by Koposov et al. and Bechtol et al., implying larger physical sizes; S per cent for Grus I and 33 per cent for Phoenix II. For Grus II, of the four RRLs detected, one is consistent with being a member of the galactic halo (D-circle dot = 24 +/- 1 kpc, mu(0) = 16.86 +/- 0.10 mag), another is at D-circle dot = 55 +/- 2 kpc (mu(0) = 18.71 +/- 0.10 mag), which we associate with Grus II, and the two remaining at D-circle dot = 43 +/- 2 kpc (mu(0) = 18.17 +/- 0.10 mag). Moreover, the appearance of a subtle red horizontal branch in the colour magnitude diagram of Grus II at the same brightness level of the latter two RRLs, which are at the same distance and in the same region, suggests that a more metal-rich system may be located in front of Grus II. The most plausible scenario is the association of these stars with the Chenab/Orphan Stream. Finally, we performed a comprehensive and updated analysis of the number of RRLs in dwarf galaxies. This allows us to predict that the method of finding new ultrafaint dwarf galaxies using two or more clumped RRLs will work only for systems brighter than M-V similar to 6 mag.490221832199FINEP - FINANCIADORA DE ESTUDOS E PROJETOSFAPERJ - FUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DO RIO DE JANEIROCNPQ - CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICOMCTIC - MINISTÉRIO DA CIÊNCIA, TECNOLOGIA, INOVAÇÕES E COMUNICAÇÕESFINEP - FINANCIADORA DE ESTUDOS E PROJETOSFAPERJ - FUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DO RIO DE JANEIROCNPQ - CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICOMCTIC - MINISTÉRIO DA CIÊNCIA, TECNOLOGIA, INOVAÇÕES E COMUNICAÇÕESSem informaçãoSem informação465376/2014-2Sem informaçãoWe thank the anonymous referee for the useful comments that helped to improve the manuscript. CEMV thanks M. Monelli for photometry advices and helpful conversations. RH is partially supported by NASA grant NNH15ZDA001N-WFIRST. Based on observations obtained at the Southern Astrophysical Research (SOAR) telescope (NOAO Prop. ID 2016A-0196; PI: Vivas), which is a joint project of the Ministerio da Ciencia, Tecnologia, Inovacoes e Comunicacoes (MCTIC) do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill (UNC), and Michigan State University (MSU). Based on observations at CTIO, National Optical Astronomy Observatory (NOAO Prop. ID: 2016A-0196, PI: Vivas; NOAO Prop. ID 2012B-0001; PI: J. Frieman), which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation (NSF). Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, FINEP - FINANCIADORA DE ESTUDOS E PROJETOS, Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and the Ministerio da Ciencia, Tecnologia e Inovacao, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the DES. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciencies de l'Espai (IEEC/CSIC), the Institut de Fisica d'Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universitat Munchen and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium. The DES data management system is supported by the National Science Foundation under grants AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-66861, FPA2015-68048, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. Institut de Fisica d'Altes Energies (IFAE) is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Program (FP7/2007-2013) including ERC grants 240672, 291329, and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020, and the Brazilian Instituto Nacional de Ciencia e Tecnologia (INCT) e-Universe (CNPq grant 465376/2014-2). This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes

Consistency of cosmic shear analyses in harmonic and real space

Recent cosmic shear studies have reported discrepancies of up to 1σ on the parameter S8=σ8Ωm/0.3‾‾‾‾‾‾‾√S8=σ8Ωm/0.3 between the analysis of shear power spectra and two-point correlation functions, derived from the same shear catalogues. It is not a priori clear whether the measured discrepancies are consistent with statistical fluctuations. In this paper, we investigate this issue in the context of the forthcoming analyses from the third year data of the Dark Energy Survey (DES Y3). We analyse DES Y3 mock catalogues from Gaussian simulations with a fast and accurate importance sampling pipeline. We show that the methodology for determining matching scale cuts in harmonic and real space is the key factor that contributes to the scatter between constraints derived from the two statistics. We compare the published scales cuts of the KiDS, Subaru-HSC, and DES surveys, and find that the correlation coefficients of posterior means range from over 80 per cent for our proposed cuts, down to 10 per cent for cuts used in the literature. We then study the interaction between scale cuts and systematic uncertainties arising from multiple sources: non-linear power spectrum, baryonic feedback, intrinsic alignments, uncertainties in the point spread function, and redshift distributions. We find that, given DES Y3 characteristics and proposed cuts, these uncertainties affect the two statistics similarly; the differential biases are below a third of the statistical uncertainty, with the largest biases arising from intrinsic alignment and baryonic feedback. While this work is aimed at DES Y3, the tools developed can be applied to Stage-IV surveys where statistical errors will be much smaller