Showing 1–50 of 188 results for author: Rozo, E

GANSky -- fast curved sky weak lensing simulations using Generative Adversarial Networks

Authors: Supranta S. Boruah, Pier Fiedorowicz, Rafael Garcia, William R. Coulton, Eduardo Rozo, Giulio Fabbian

Abstract: Extracting non-Gaussian information from the next generation weak lensing surveys will require fast and accurate full-sky simulations. This is difficult to achieve in practice with existing simulation methods: ray-traced $N$-body simulations are computationally expensive, and approximate simulation methods (such as lognormal mocks) are not accurate enough. Here, we present GANSky, an interpretable… ▽ More Extracting non-Gaussian information from the next generation weak lensing surveys will require fast and accurate full-sky simulations. This is difficult to achieve in practice with existing simulation methods: ray-traced $N$-body simulations are computationally expensive, and approximate simulation methods (such as lognormal mocks) are not accurate enough. Here, we present GANSky, an interpretable machine learning method that uses Generative Adversarial Networks (GANs) to produce fast and accurate full-sky tomographic weak lensing maps. The input to our GAN are lognormal maps that approximately describe the late-time convergence field of the Universe. Starting from these lognormal maps, we use GANs to learn how to locally redistribute mass to achieve simulation-quality maps. This can be achieved using remarkably small networks ($\approx 10^3$ parameters). We validate the GAN maps by computing a number of summary statistics in both simulated and GANSky maps. We show that GANSky maps correctly reproduce both the mean and $χ^2$ distribution for several statistics, specifically: the 2-pt function, 1-pt PDF, peak and void counts, and the equilateral, folded and squeezed bispectra. These successes makes GANSky an attractive tool to compute the covariances of these statistics. In addition to being useful for rapidly generating large ensembles of artificial data sets, our method can be used to extract non-Gaussian information from weak lensing data with field-level or simulation-based inference. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: 14 pages, 10 Figures, comments welcome

Dynamics-based halo model for large scale structure

Authors: Edgar M. Salazar, Eduardo Rozo, Rafael García, 1 Nickolas Kokron, Susmita Adhikari, Benedikt Diemer, Calvin Osinga

Abstract: Accurate modelling of the one-to-two halo transition has long been difficult to achieve. We demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to: 1) particles \it orbiting \rm the halo; and 2) particles \it… ▽ More Accurate modelling of the one-to-two halo transition has long been difficult to achieve. We demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to: 1) particles \it orbiting \rm the halo; and 2) particles \it infalling \rm into the halo for the first time. Motivated by this decomposition, García [R. García et. al., MNRAS 521, 2464 (2023)] advocated for defining haloes as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo--mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterised by a single physical scale $r_{\rm h}$, which we refer to as the \it halo radius. \rm The two-halo term is described using a simple yet accurate empirical model based on the Zel'dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo--mass correlation function is accurate at the $\approx 2\%$ level for $r \in [0.1,50]\ h^{-1}\ Mpc$. The Fourier transform of our best fit model describes the halo--mass power spectrum with comparable accuracy for $k\in [0.06, 6.0]\ h\ Mpc^{-1}$. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 17 pages, 14 figures

Bayesian mass mapping with weak lensing data using KARMMA -- validation with simulations and application to Dark Energy Survey Year 3 data

Authors: Supranta S. Boruah, Pier Fiedorowicz, Eduardo Rozo

Abstract: We update the field-level inference code KARMMA to enable tomographic forward-modelling of shear maps. Our code assumes a lognormal prior on the convergence field, and properly accounts for the cross-covariance in the lensing signal across tomographic source bins. We use mock weak lensing data from N-body simulations to validate our mass-mapping forward model by comparing our posterior maps to the… ▽ More We update the field-level inference code KARMMA to enable tomographic forward-modelling of shear maps. Our code assumes a lognormal prior on the convergence field, and properly accounts for the cross-covariance in the lensing signal across tomographic source bins. We use mock weak lensing data from N-body simulations to validate our mass-mapping forward model by comparing our posterior maps to the input convergence fields. We find that KARMMA produces more accurate reconstructions than traditional mass-mapping algorithms. More-over, the KARMMA posteriors reproduce all statistical properties of the input density field we tested -- one- and two-point functions, and the peak and void number counts -- with $\lesssim~10\%$ accuracy. Our posteriors exhibit a small bias that increases with decreasing source redshift, but these biases are small compared to the statistical uncertainties of current (DES) cosmic shear surveys. Finally, we apply KARMMA to Dark Energy Survey Year 3 (DES-Y3) weak lensing data, and verify that the two point shear correlation function $ξ_+$ is well fit by the correlation function of the reconstructed convergence field. This is a non-trivial test that traditional mass mapping algorithms fail. The code is publicly available at https://github.com/Supranta/KaRMMa.git. KARMMA DES-Y3 mass maps are publicly available at https://zenodo.org/records/10672062. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 14 pages, 13 figures, To be submitted to PRD, Comments welcome

Forecasting the constraints on optical selection bias and projection effects of galaxy cluster lensing with multiwavelength data

Authors: Conghao Zhou, Hao-Yi Wu, Andrés N. Salcedo, Sebastian Grandis, Tesla Jeltema, Alexie Leauthaud, Matteo Costanzi, Tomomi Sunayama, David H. Weinberg, Tianyu Zhang, Eduardo Rozo, Chun-Hao To, Sebastian Bocquet, Tamas Varga, Matthew Kwiecien

Abstract: Galaxy clusters identified with optical imaging tend to suffer from projection effects, which impact richness (the number of member galaxies in a cluster) and lensing coherently. Physically unassociated galaxies can be mistaken as cluster members due to the significant uncertainties in their line-of-sight distances, thereby changing the observed cluster richness; at the same time, projection effec… ▽ More Galaxy clusters identified with optical imaging tend to suffer from projection effects, which impact richness (the number of member galaxies in a cluster) and lensing coherently. Physically unassociated galaxies can be mistaken as cluster members due to the significant uncertainties in their line-of-sight distances, thereby changing the observed cluster richness; at the same time, projection effects alter the weak gravitational lensing signals of clusters, leading to a correlated scatter between richness and lensing at a given halo mass. As a result, the lensing signals for optically selected clusters tend to be biased high. This optical selection bias problem of cluster lensing is one of the key challenges in cluster cosmology. Fortunately, recently available multiwavelength observations of clusters provide a solution. We analyze a simulated data set mimicking the observed lensing of clusters identified by both optical photometry and gas properties, aiming to constrain this selection bias. Assuming a redMaPPer sample from the Dark Energy Survey with South Pole Telescope Sunyaev-Zeldovich effect observations, we find that an overlapping survey of 1300 square deg, 0.2 < z < 0.65, can constrain the average lensing bias to an accuracy of 5 percent. This provides an exciting opportunity for directly constraining optical selection bias from observations. We further show that our approach can remove the optical selection bias from the lensing signal, paving the way for future optical cluster cosmology analyses. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 16 pages, 5 figures. Submitted to PRD

Dark Energy Survey Year 1 Clusters are Consistent with Planck

Authors: Andrés N. Salcedo, Hao-Yi Wu, Eduardo Rozo, David H. Weinberg, Chun-Hao To, Tomomi Sunayama, Andy Lee

Abstract: The recent Dark Energy Survey Year 1 (DES-Y1) analysis of galaxy cluster abundances and weak lensing produced $Ω_{\rm m}$ and $σ_8$ constraints in 5.6$σ$ tension with Planck. It is suggested in that work that this tension is driven by unmodelled systematics in optical cluster selection. We present a novel simulation-based forward modeling framework that explicitly incorporates cluster selection in… ▽ More The recent Dark Energy Survey Year 1 (DES-Y1) analysis of galaxy cluster abundances and weak lensing produced $Ω_{\rm m}$ and $σ_8$ constraints in 5.6$σ$ tension with Planck. It is suggested in that work that this tension is driven by unmodelled systematics in optical cluster selection. We present a novel simulation-based forward modeling framework that explicitly incorporates cluster selection into its model predictions. Applying this framework to the DES-Y1 data we find consistency with Planck, resolving the tension found in the DES-Y1 analysis. An extension of this approach to the final DES data set will produce robust constraints on $Λ$CDM parameters and correspondingly strong tests of cosmological models. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 6 pages, 2 figures, 1 table, Supplemental material with 2 figures. Submitted to Physical Review Letters

Map-based cosmology inference with weak lensing -- information content and its dependence on the parameter space

Authors: Supranta S. Boruah, Eduardo Rozo

Abstract: Field-level inference is emerging as a promising technique for optimally extracting information from cosmological datasets. Indeed, previous analyses have shown field-based inference produces tighter parameter constraints than power spectrum analyses. However, estimates of the detailed quantitative gain in constraining power differ. Here, we demonstrate the gain in constraining power depends on th… ▽ More Field-level inference is emerging as a promising technique for optimally extracting information from cosmological datasets. Indeed, previous analyses have shown field-based inference produces tighter parameter constraints than power spectrum analyses. However, estimates of the detailed quantitative gain in constraining power differ. Here, we demonstrate the gain in constraining power depends on the parameter space being constrained. As a specific example, we find that field-based analysis of an LSST Y1-like mock data set only marginally improves constraints relative to a 2-point function analysis in $Λ$CDM, yet it more than doubles the constraining power of the data in the context of $w$CDM models. This effect reconciles some, but not all, of the discrepant results found in the literature. Our results demonstrate the importance of using a full systematics model when quantifying the information gain for realistic field-level analyses of future data sets. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: 5 pages, 4 figures, comments welcome

Buzzard to Cardinal: Improved Mock Catalogs for Large Galaxy Surveys

Authors: Chun-Hao To, Joseph DeRose, Risa H. Wechsler, Eli Rykoff, Hao-Yi Wu, Susmita Adhikari, Elisabeth Krause, Eduardo Rozo, David H. Weinberg

Abstract: We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including DES, DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of $z=2.35$ to a depth of $m_{\rm{r}}=27$. Compared to the Buzzard mocks, the Cardinal mocks include an update… ▽ More We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including DES, DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of $z=2.35$ to a depth of $m_{\rm{r}}=27$. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching (SHAM) model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group--galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to improve the color assignment model further and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the Dark Energy Survey data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys. The cardinal mock will be released upon publication at https://chunhaoto.com/cardinalsim. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 37 pages; 20 figures. See figure 1 for summary and appendix K for a list of main differences between Cardinal and Buzzard. To be submitted to APJ; Comments welcome. Highlights can be found at https://chunhaoto.com/cardinalsim

The Intrinsic Alignment of Red Galaxies in DES Y1 redMaPPer Galaxy Clusters

Authors: C. Zhou, A. Tong, M. A. Troxel, J. Blazek, C. Lin, D. Bacon, L. Bleem, A. Carnero Rosell, C. Chang, M. Costanzi, J. DeRose, J. P. Dietrich, A. Drlica-Wagner, D. Gruen, R. A. Gruendl, B. Hoyle, M. Jarvis, N. MacCrann, B. Mawdsley, T. McClintock, P. Melchior, J. Prat, A. Pujol, E. Rozo, E. S. Rykoff , et al. (57 additional authors not shown)

Abstract: Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We meas… ▽ More Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 redMaPPer clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshift 0.1-0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude ($A_{\textrm{IA}}$) to the measurement, finding $A_{\textrm{IA}}=0.15\pm 0.04$, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modeling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies. △ Less

Submitted 5 September, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 13 figures. Accepted to MNRAS

KaRMMa 2.0 -- Kappa Reconstruction for Mass Mapping

Authors: Pier Fiedorowicz, Eduardo Rozo, Supranta S. Boruah

Abstract: We present KaRMMa 2.0, an updated version of the mass map reconstruction code introduced in Fiedorowicz et al. (2022). KaRMMa is a full-sky Bayesian algorithm for reconstructing weak lensing mass maps from shear data. It forward-models the convergence field as a realization of a lognormal field. The corresponding shear map is calculated using the standard Kaiser-Squires transformation, and compare… ▽ More We present KaRMMa 2.0, an updated version of the mass map reconstruction code introduced in Fiedorowicz et al. (2022). KaRMMa is a full-sky Bayesian algorithm for reconstructing weak lensing mass maps from shear data. It forward-models the convergence field as a realization of a lognormal field. The corresponding shear map is calculated using the standard Kaiser-Squires transformation, and compared to observations at the field level. The posterior distribution of maps given the shear data is sampled using Hamiltonian Monte Carlo chains. Our work improves on the original algorithm by making it numerically efficient, enabling full-sky reconstructions at $\approx$ 7 arcmin resolution with modest computational resources. These gains are made with no loss in accuracy or precision relative to KaRMMa 1.0. We compare the KaRMMa 2.0 posteriors against simulations across a variety of summary statistics (one-point function, two-point functions, and peak/void counts) to demonstrate our updated algorithm provides an accurate reconstruction of the convergence field at mildly non-linear scales. Unsurprisingly, the lognormal model fails as we approach non-linear scales ($\ell \gtrsim 200$), which in turn biases the map posteriors. These biases are at the 2% level in the recovered power spectrum, and at the 5% to 15% level for other statistics, depending on the resolution. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 12 pages, 9 figures

A measurement of the mean central optical depth of galaxy clusters via the pairwise kinematic Sunyaev-Zel'dovich effect with SPT-3G and DES

Authors: E. Schiappucci, F. Bianchini, M. Aguena, M. Archipley, L. Balkenhol, L. E. Bleem, P. Chaubal, T. M. Crawford, S. Grandis, Y. Omori, C. L. Reichardt, E. Rozo, E. S. Rykoff, C. To, T. M. C. Abbott, P. A. R. Ade, O. Alves, A. J. Anderson, F. Andrade-Oliveira, J. Annis, J. S. Avva, D. Bacon, K. Benabed, A. N. Bender, B. A. Benson , et al. (117 additional authors not shown)

Abstract: We infer the mean optical depth of a sample of optically-selected galaxy clusters from the Dark Energy Survey (DES) via the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect. The pairwise kSZ signal between pairs of clusters drawn from the DES Year-3 cluster catalog is detected at $4.1 σ$ in cosmic microwave background (CMB) temperature maps from two years of observations with the SPT-3G camera o… ▽ More We infer the mean optical depth of a sample of optically-selected galaxy clusters from the Dark Energy Survey (DES) via the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect. The pairwise kSZ signal between pairs of clusters drawn from the DES Year-3 cluster catalog is detected at $4.1 σ$ in cosmic microwave background (CMB) temperature maps from two years of observations with the SPT-3G camera on the South Pole Telescope. After cuts, there are 24,580 clusters in the $\sim 1,400$ deg$^2$ of the southern sky observed by both experiments. We infer the mean optical depth of the cluster sample with two techniques. The optical depth inferred from the pairwise kSZ signal is $\barτ_e = (2.97 \pm 0.73) \times 10^{-3}$, while that inferred from the thermal SZ signal is $\barτ_e = (2.51 \pm 0.55^{\text{stat}} \pm 0.15^{\rm syst}) \times 10^{-3}$. The two measures agree at $0.6 σ$. We perform a suite of systematic checks to test the robustness of the analysis. △ Less

Submitted 16 June, 2023; v1 submitted 25 July, 2022; originally announced July 2022.

A Better Way to Define Dark Matter Haloes

Authors: Rafael Garcia, Edgar Salazar, Eduardo Rozo, Susmita Adhikari, Han Aung, Benedikt Diemer, Daisuke Nagai, Brandon Wolfe

Abstract: Dark matter haloes have long been recognized as one of the fundamental building blocks of large scale structure formation models. Despite their importance -- or perhaps because of it! -- halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, and effectively unique and parameter-free: ''A dark matter halo is co… ▽ More Dark matter haloes have long been recognized as one of the fundamental building blocks of large scale structure formation models. Despite their importance -- or perhaps because of it! -- halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, and effectively unique and parameter-free: ''A dark matter halo is comprised of the collection of particles orbiting in their own self-generated potential.'' This definition is enabled by the fact that, even with as few as $\approx 300$ particles per halo, nearly every particle in the vicinity of a halo can be uniquely classified as either orbiting or infalling based on its dynamical history. For brevity, we refer to haloes selected in this way as physical haloes. We demonstrate that: 1) the mass function of physical haloes is Press-Schechter, provided the critical threshold for collapse is allowed to vary slowly with peak height; and 2) the peak-background split prediction of the clustering amplitude of physical halos is statistically consistent with the simulation data, with an accuracy no worse than $\approx 5\%$. △ Less

Submitted 24 July, 2022; originally announced July 2022.

Comments: 16 pages, 14 figures

C2-GaMe: Classification of Cluster Galaxy Membership with Machine Learning

Authors: Daniel Farid, Han Aung, Daisuke Nagai, Arya Farahi, Eduardo Rozo

Abstract: We present Classification of Cluster GAlaxy MEmbers (C$^2$-GaMe), a classification algorithm based on a suite of machine learning models that differentiates galaxies into orbiting, infalling, and background (interloper) populations, using phase space information as input. We train and test C$^2$-GaMe with the galaxies from UniverseMachine mock catalog based on Multi-Dark Planck 2 N-body simulation… ▽ More We present Classification of Cluster GAlaxy MEmbers (C$^2$-GaMe), a classification algorithm based on a suite of machine learning models that differentiates galaxies into orbiting, infalling, and background (interloper) populations, using phase space information as input. We train and test C$^2$-GaMe with the galaxies from UniverseMachine mock catalog based on Multi-Dark Planck 2 N-body simulations. We show that probabilistic classification is superior to deterministic classification in estimating the physical properties of clusters, including density profiles and velocity dispersion. We propose a set of estimators to get an unbiased estimation of cluster properties. We demonstrate that C$^2$-GaMe can recover the distribution of orbiting and infalling galaxies' position and velocity distribution with $<1\%$ statistical error when using probabilistic predictions in the presence of interlopers in the projected phase space. Additionally, we demonstrate the robustness of trained models by applying them to a different simulation. Finally, adding a specific star formation rate and the ratio of the galaxy's halo mass to the cluster's halo mass as additional features improves the classification performance. We discuss potential applications of this technique to enhance cluster cosmology and galaxy quenching. △ Less

Submitted 22 September, 2023; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted in Astronomy and Computing. Some figures are downsampled from journal version to appear correctly in Chrome

Map-based cosmology inference with lognormal cosmic shear maps

Authors: Supranta Sarma Boruah, Eduardo Rozo, Pier Fiedorowicz

Abstract: Most cosmic shear analyses to date have relied on summary statistics (e.g. $ξ_+$ and $ξ_-$). These types of analyses are necessarily sub-optimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint po… ▽ More Most cosmic shear analyses to date have relied on summary statistics (e.g. $ξ_+$ and $ξ_-$). These types of analyses are necessarily sub-optimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative to single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps -- and to the extent that the convergence field can be described by the lognormal model -- our map posteriors allow us to reconstruct \it all \rm summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the $σ_8$--$Ω_{\rm m}$ plane by $\approx 30\%$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $S_8=σ_8Ω_{\rm m}^{1/2}$ directions, meaning map-based inference fails to significantly improve constraints on $S_8$. △ Less

Submitted 27 April, 2022; originally announced April 2022.

Comments: 12 pages, 11 figures, To be submitted to MNRAS, Comments welcome

Accurate Model of the Projected Velocity Distribution of Galaxies in Dark Matter Halos

Authors: Han Aung, Daisuke Nagai, Eduardo Rozo, Brandon Wolfe, Susmita Adhikari

Abstract: We present a percent-level accurate model of the line-of-sight velocity distribution of galaxies around dark matter halos as a function of projected radius and halo mass. The model is developed and tested using synthetic galaxy catalogs generated with the UniverseMachine run on the Multi-Dark Planck 2 N-body simulations. The model decomposes the galaxies around a cluster into three kinematically d… ▽ More We present a percent-level accurate model of the line-of-sight velocity distribution of galaxies around dark matter halos as a function of projected radius and halo mass. The model is developed and tested using synthetic galaxy catalogs generated with the UniverseMachine run on the Multi-Dark Planck 2 N-body simulations. The model decomposes the galaxies around a cluster into three kinematically distinct classes: orbiting, infalling, and interloping galaxies. We demonstrate that: 1) we can statistically distinguish between these three types of galaxies using only projected line-of-sight velocity information; 2) the halo edge radius inferred from the line-of-sight velocity dispersion is an excellent proxy for the three-dimensional halo edge radius; 3) we can accurately recover the full velocity dispersion profile for each of the three populations of galaxies. Importantly, the velocity dispersion profiles of the orbiting and infalling galaxies contain five independent parameters -- three distinct radial scales and two velocity dispersion amplitudes -- each of which is correlated with mass. Thus, the velocity dispersion profile of galaxy clusters has inherent redundancies that allow us to perform nontrivial systematics check from a single data set. We discuss several potential applications of our new model for detecting the edge radius and constraining cosmology and astrophysics using upcoming spectroscopic surveys. △ Less

Submitted 26 January, 2023; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 10 pages, 11 figures, accepted in MNRAS

LINNA: Likelihood Inference Neural Network Accelerator

Authors: Chun-Hao To, Eduardo Rozo, Elisabeth Krause, Hao-Yi Wu, Risa H. Wechsler, Andrés N. Salcedo

Abstract: Bayesian posterior inference of modern multi-probe cosmological analyses incurs massive computational costs. For instance, depending on the combinations of probes, a single posterior inference for the Dark Energy Survey (DES) data had a wall-clock time that ranged from 1 to 21 days using a state-of-the-art computing cluster with 100 cores. These computational costs have severe environmental impact… ▽ More Bayesian posterior inference of modern multi-probe cosmological analyses incurs massive computational costs. For instance, depending on the combinations of probes, a single posterior inference for the Dark Energy Survey (DES) data had a wall-clock time that ranged from 1 to 21 days using a state-of-the-art computing cluster with 100 cores. These computational costs have severe environmental impacts and the long wall-clock time slows scientific productivity. To address these difficulties, we introduce LINNA: the Likelihood Inference Neural Network Accelerator. Relative to the baseline DES analyses, LINNA reduces the computational cost associated with posterior inference by a factor of 8--50. If applied to the first-year cosmological analysis of Rubin Observatory's Legacy Survey of Space and Time (LSST Y1), we conservatively estimate that LINNA will save more than US $\$300,000$ on energy costs, while simultaneously reducing $\rm{CO}_2$ emission by $2,400$ tons. To accomplish these reductions, LINNA automatically builds training data sets, creates neural network surrogate models, and produces a Markov chain that samples the posterior. We explicitly verify that LINNA accurately reproduces the first-year DES (DES Y1) cosmological constraints derived from a variety of different data vectors with our default code settings, without needing to retune the algorithm every time. Further, we find that LINNA is sufficient for enabling accurate and efficient sampling for LSST Y10 multi-probe analyses. We make LINNA publicly available at https://github.com/chto/linna, to enable others to perform fast and accurate posterior inference in contemporary cosmological analyses. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: 21 pages, 12 figures, submitted to JCAP, comments are welcome

Optical selection bias and projection effects in stacked galaxy cluster weak lensing

Authors: Hao-Yi Wu, Matteo Costanzi, Chun-Hao To, Andrés N. Salcedo, David H. Weinberg, James Annis, Sebastian Bocquet, Maria Elidaiana da Silva Pereira, Joseph DeRose, Johnny Esteves, Arya Farahi, Sebastian Grandis, Eduardo Rozo, Eli S. Rykoff, Tamás N. Varga, Risa H. Wechsler, Chenxiao Zeng, Yuanyuan Zhang, Zhuowen Zhang

Abstract: Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters ident… ▽ More Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters identified by the redMaPPer algorithm in the Buzzard simulations as a case study. We find that at a given cluster halo mass, the residuals of redMaPPer richness and weak lensing signal are positively correlated. As a result, for a given richness selection, the stacked lensing signal is biased high compared with what we would expect from the underlying halo mass probability distribution. The cluster lensing selection bias can thus lead to overestimated mean cluster mass and biased cosmology results. We show that the lensing selection bias exhibits a strong scale-dependence and is approximately 20 to 60 percent for $ΔΣ$ at large scales. This selection bias largely originates from spurious member galaxies within +/- 20 to 60 Mpc/h along the line of sight, highlighting the importance of quantifying projection effects associated with the broad redshift distribution of member galaxies in photometric cluster surveys. While our results qualitatively agree with those in the literature, accurate quantitative modelling of the selection bias is needed to achieve the goals of cluster lensing cosmology and will require synthetic catalogues covering a wide range of galaxy-halo connection models. △ Less

Submitted 9 August, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: 16 pages, 16 figures; replaced to match published version

Measuring Cosmological Parameters with Type Ia Supernovae in redMaGiC galaxies

Authors: R. Chen, D. Scolnic, E. Rozo, E. S. Rykoff, B. Popovic, R. Kessler, M. Vincenzi, T. M. Davis, P. Armstrong, D. Brout, L. Galbany, L. Kelsey, C. Lidman, A. Möller, B. Rose, M. Sako, M. Sullivan, G. Taylor, P. Wiseman, J. Asorey, A. Carr, C. Conselice, K. Kuehn, G. F. Lewis, E. Macaulay , et al. (60 additional authors not shown)

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

Submitted 21 February, 2022; originally announced February 2022.

Comments: 15 pages, 8 figures. Submitting to ApJ, comments welcome

Journal ref: ApJ 938 62 (2022)

Incorporating galaxy cluster triaxiality in stacked cluster weak lensing analyses

Authors: Zhuowen Zhang, Hao-Yi Wu, Yuanyuan Zhang, Joshua Frieman, Chun-Hao To, Joseph DeRose, Matteo Costanzi, Risa H. Wechsler, Susmita Adhikari, Eli Rykoff, Tesla Jeltema, August Evrard, Eduardo Rozo

Abstract: Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis… ▽ More Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters' shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modeling the richness-mass relation as a log-linear relation, we find that the log-richness amplitude $\ln(A)$ is boosted from the lowest to highest orientation bin with a significance of $14σ$, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark halos in different orientation bins resemble a "bottleneck" shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology -- miscentering and projection -- and find a null correlation. Analytic templates for the triaxiality bias of observed-richness and lensing profiles are mapped as corrections to the observable of richness-binned lensing profiles for redMaPPer clusters. The resulting mass bias confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey: I. Evidence for thermal energy anisotropy using oriented stacking

Authors: M. Lokken, R. Hložek, A. van Engelen, M. Madhavacheril, E. Baxter, J. DeRose, C. Doux, S. Pandey, E. S. Rykoff, G. Stein, C. To, T. M. C. Abbott, S. Adhikari, M. Aguena, S. Allam, F. Andrade-Oliveira, J. Annis, N. Battaglia, G. M. Bernstein, E. Bertin, J. R. Bond, D. Brooks, E. Calabrese, A. Carnero Rosell, M. Carrasco Kind , et al. (82 additional authors not shown)

Abstract: The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel'dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from… ▽ More The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along inter-cluster bridges, visible through their thermal Sunyaev-Zel'dovich signal in the Cosmic Microwave Background. We demonstrate a new, flexible method to analyze the hot gas signal from multi-scale extended structures. We use a Compton-$y$ map from the Atacama Cosmology Telescope (ACT) stacked on redMaPPer cluster positions from the optical Dark Energy Survey (DES). Cutout images from the $y$ map are oriented with large-scale structure information from DES galaxy data such that the superclustering signal is aligned before being overlaid. We find evidence for an extended quadrupole moment of the stacked $y$ signal at the 3.5$σ$ level, demonstrating that the large-scale thermal energy surrounding galaxy clusters is anisotropically distributed. We compare our ACT$\times$DES results with the Buzzard simulations, finding broad agreement. Using simulations, we highlight the promise of this novel technique for constraining the evolution of anisotropic, non-Gaussian structure using future combinations of microwave and optical surveys. △ Less

Submitted 18 July, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 37 pages, 23 figures, 4 tables. Added explanatory figure, table, covariance matrix equations, discussion of CIB impact. Matches the version published in ApJ

KaRMMa -- Kappa Reconstruction for Mass Mapping

Authors: Pier Fiedorowicz, Eduardo Rozo, Supranta S. Boruah, Chihway Chang, Marco Gatti

Abstract: We present KaRMMa, a novel method for performing mass map reconstruction from weak-lensing surveys. We employ a fully Bayesian approach with a physically motivated lognormal prior to sample from the posterior distribution of convergence maps. We test KaRMMa on a suite of dark matter N-body simulations with simulated DES Y1-like shear observations. We show that KaRMMa outperforms the basic Kaiser-S… ▽ More We present KaRMMa, a novel method for performing mass map reconstruction from weak-lensing surveys. We employ a fully Bayesian approach with a physically motivated lognormal prior to sample from the posterior distribution of convergence maps. We test KaRMMa on a suite of dark matter N-body simulations with simulated DES Y1-like shear observations. We show that KaRMMa outperforms the basic Kaiser-Squires mass map reconstruction in two key ways: 1) our best map point estimate has lower residuals compared to Kaiser-Squires; and 2) unlike the Kaiser-Squires reconstruction, the posterior distribution of KaRMMa maps are nearly unbiased in all summary statistics we considered, namely: one-point and two-point functions, and peak/void counts. In particular, KaRMMa successfully captures the non-Gaussian nature of the distribution of $κ$ values in the simulated maps. We further demonstrate that the KaRMMa posteriors correctly characterize the uncertainty in all summary statistics we considered. △ Less

Submitted 25 January, 2022; v1 submitted 31 May, 2021; originally announced May 2021.

Comments: 13 pages, 11 figures

Dark Energy Survey Year 3 Results: Constraints on cosmological parameters and galaxy bias models from galaxy clustering and galaxy-galaxy lensing using the redMaGiC sample

Authors: S. Pandey, E. Krause, J. DeRose, N. MacCrann, B. Jain, M. Crocce, J. Blazek, A. Choi, H. Huang, C. To, X. Fang, J. Elvin-Poole, J. Prat, A. Porredon, L. F. Secco, M. Rodriguez-Monroy, N. Weaverdyck, Y. Park, M. Raveri, E. Rozo, E. S. Rykoff, G. M. Bernstein, C. Sánchez, M. Jarvis, M. A. Troxel , et al. (116 additional authors not shown)

Abstract: We constrain cosmological and galaxy-bias parameters using the combination of galaxy clustering and galaxy-galaxy lensing measurements from the Dark Energy Survey Year-3 data. We describe our modeling framework, and choice of scales analyzed, validating their robustness to theoretical uncertainties in small-scale clustering by analyzing simulated data. Using a linear galaxy bias model and redMaGiC… ▽ More We constrain cosmological and galaxy-bias parameters using the combination of galaxy clustering and galaxy-galaxy lensing measurements from the Dark Energy Survey Year-3 data. We describe our modeling framework, and choice of scales analyzed, validating their robustness to theoretical uncertainties in small-scale clustering by analyzing simulated data. Using a linear galaxy bias model and redMaGiC galaxy sample, we obtain constraints on the matter density to be $Ω_{\rm m} = 0.325^{+0.033}_{-0.034}$. We also implement a non-linear galaxy bias model to probe smaller scales that includes parameterization based on hybrid perturbation theory and find that it leads to a 17% gain in cosmological constraining power. We perform robustness tests of our methodology pipeline and demonstrate the stability of the constraints to changes in the theoretical model. Using the redMaGiC galaxy sample as foreground lens galaxies, we find the galaxy clustering and galaxy-galaxy lensing measurements to exhibit significant signals akin to de-correlation between galaxies and mass on large scales, which is not expected in any current models. This likely systematic measurement error biases our constraints on galaxy bias and the $S_8$ parameter. We find that a scale-, redshift- and sky-area-independent phenomenological de-correlation parameter can effectively capture the impact of this systematic error. We trace the source of this de-correlation to a color-dependent photometric issue and minimize its impact on our result by changing the selection criteria of redMaGiC galaxies. Using this new sample, our constraints on the $S_8$ parameter are consistent with previous studies, and we find a small shift in the $Ω_{\rm m}$ constraints compared to the fiducial redMaGiC sample. We constrain the mean host halo mass of the redMaGiC galaxies in this new sample to be approximately $1.6 \times 10^{13} M_{\odot}/h$. △ Less

Submitted 25 September, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: 34 pages, 23 figures, includes the updated results from a new redmagic sample giving S8 consistent with the cosmic shear results (and also prefers Xlens=1)

Dark Energy Survey Year 3 Results: Galaxy clustering and systematics treatment for lens galaxy samples

Authors: M. Rodríguez-Monroy, N. Weaverdyck, J. Elvin-Poole, M. Crocce, A. Carnero Rosell, F. Andrade-Oliveira, S. Avila, K. Bechtol, G. M. Bernstein, J. Blazek, H. Camacho, R. Cawthon, J. De Vicente, J. DeRose, S. Dodelson, S. Everett, X. Fang, I. Ferrero, A. Ferté, O. Friedrich, E. Gaztanaga, G. Giannini, R. A. Gruendl, W. G. Hartley, K. Herner , et al. (80 additional authors not shown)

Abstract: In this work we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, MagLim, and a sample of luminous red galaxies selected with the redMaGiC algorithm. MagLim / redMaGiC sample contains over 10 million / 2.5 million galaxies and is divided into six / five photometric redshift bins spanni… ▽ More In this work we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, MagLim, and a sample of luminous red galaxies selected with the redMaGiC algorithm. MagLim / redMaGiC sample contains over 10 million / 2.5 million galaxies and is divided into six / five photometric redshift bins spanning the range $z\in[0.20,1.05]$ / $z\in[0.15,0.90]$. Both samples cover 4143 deg$^2$ over which we perform our analysis blind, measuring the angular correlation function with a S/N $\sim 63$ for both samples. In a companion paper (DES Collaboration et al. 2021)), these measurements of galaxy clustering are combined with the correlation functions of cosmic shear and galaxy-galaxy lensing of each sample to place cosmological constraints with a 3$\times$2pt analysis. We conduct a thorough study of the mitigation of systematic effects caused by the spatially varying survey properties and we correct the measurements to remove artificial clustering signals. We employ several decontamination methods with different configurations to ensure the robustness of our corrections and to determine the systematic uncertainty that needs to be considered for the final cosmology analyses. We validate our fiducial methodology using log-normal mocks, showing that our decontamination procedure induces biases no greater than $0.5σ$ in the $(Ω_m, b)$ plane, where $b$ is galaxy bias. We demonstrate that failure to remove the artificial clustering would introduce strong biases up to $\sim 7 σ$ in $Ω_m$ and of more than $4 σ$ in galaxy bias. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: 23 pages, 19 figures, see https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ for the full DES Y3 cosmology release

The HST See Change Program: I. Survey Design, Pipeline, and Supernova Discoveries

Authors: Brian Hayden, David Rubin, Kyle Boone, Greg Aldering, Jakob Nordin, Mark Brodwin, Susana Deustua, Sam Dixon, Parker Fagrelius, Andy Fruchter, Peter Eisenhardt, Anthony Gonzalez, Ravi Gupta, Isobel Hook, Chris Lidman, Kyle Luther, Adam Muzzin, Zachary Raha, Pilar Ruiz-Lapuente, Clare Saunders, Caroline Sofiatti, Adam Stanford, Nao Suzuki, Tracy Webb, Steven C. Williams , et al. (31 additional authors not shown)

Abstract: The See Change survey was designed to make $z>1$ cosmological measurements by efficiently discovering high-redshift Type Ia supernovae (SNe Ia) and improving cluster mass measurements through weak lensing. This survey observed twelve galaxy clusters with the Hubble Space Telescope spanning the redshift range $z=1.13$ to $1.75$, discovering 57 likely transients and 27 likely SNe Ia at… ▽ More The See Change survey was designed to make $z>1$ cosmological measurements by efficiently discovering high-redshift Type Ia supernovae (SNe Ia) and improving cluster mass measurements through weak lensing. This survey observed twelve galaxy clusters with the Hubble Space Telescope spanning the redshift range $z=1.13$ to $1.75$, discovering 57 likely transients and 27 likely SNe Ia at $z\sim 0.8-2.3$. As in similar previous surveys (Dawson et al. 2009), this proved to be a highly efficient use of HST for SN observations; the See Change survey additionally tested the feasibility of maintaining, or further increasing, the efficiency at yet higher redshifts, where we have less detailed information on the expected cluster masses and star-formation rates. We find that the resulting number of SNe Ia per orbit is a factor of $\sim 8$ higher than for a field search, and 45% of our orbits contained an active SN Ia within 22 rest-frame days of peak, with one of the clusters by itself yielding 6 of the SNe Ia. We present the survey design, pipeline, and SN discoveries. Novel features include fully blinded SN searches, the first random forest candidate classifier for undersampled IR data (with a 50% detection threshold within 0.05 magnitudes of human searchers), real-time forward-modeling photometry of candidates, and semi-automated photometric classifications and follow-up forecasts. We also describe the spectroscopic follow-up, instrumental in measuring host-galaxy redshifts. The cosmology analysis of our sample will be presented in a companion paper. △ Less

Submitted 24 March, 2021; originally announced March 2021.

Comments: ApJ preprint

Synthetic Galaxy Clusters and Observations Based on Dark Energy Survey Year 3 Data

Authors: T. N. Varga, D. Gruen, S. Seitz, N. MacCrann, E. Sheldon, W. G. Hartley, A. Amon, A. Choi, A. Palmese, Y. Zhang, M. R. Becker, J. McCullough, E. Rozo, E. S. Rykoff, C. To, S. Grandis, G. M. Bernstein, S. Dodelson, K. Eckert, S. Everett, R. A. Gruendl, I. Harrison, K. Herner, R. P. Rollins, I. Sevilla-Noarbe , et al. (53 additional authors not shown)

Abstract: We develop a novel data-driven method for generating synthetic optical observations of galaxy clusters. In cluster weak lensing, the interplay between analysis choices and systematic effects related to source galaxy selection, shape measurement and photometric redshift estimation can be best characterized in end-to-end tests going from mock observations to recovered cluster masses. To create such… ▽ More We develop a novel data-driven method for generating synthetic optical observations of galaxy clusters. In cluster weak lensing, the interplay between analysis choices and systematic effects related to source galaxy selection, shape measurement and photometric redshift estimation can be best characterized in end-to-end tests going from mock observations to recovered cluster masses. To create such test scenarios, we measure and model the photometric properties of galaxy clusters and their sky environments from the Dark Energy Survey Year 3 (DES Y3) data in two bins of cluster richness $λ\in[30;\,45)$, $λ\in[45;\,60)$ and three bins in cluster redshift ($z\in[0.3;\,0.35)$, $z\in[0.45;\,0.5)$ and $z\in[0.6;\,0.65)$. Using deep-field imaging data we extrapolate galaxy populations beyond the limiting magnitude of DES Y3 and calculate the properties of cluster member galaxies via statistical background subtraction. We construct mock galaxy clusters as random draws from a distribution function, and render mock clusters and line-of-sight catalogs into synthetic images in the same format as actual survey observations. Synthetic galaxy clusters are generated from real observational data, and thus are independent from the assumptions inherent to cosmological simulations. The recipe can be straightforwardly modified to incorporate extra information, and correct for survey incompleteness. New realizations of synthetic clusters can be created at minimal cost, which will allow future analyses to generate the large number of images needed to characterize systematic uncertainties in cluster mass measurements. △ Less

Submitted 31 January, 2022; v1 submitted 20 February, 2021; originally announced February 2021.

Comments: accepted to MNRAS, 22 pages, 14 figures

CODEX Weak Lensing Mass Catalogue and implications on the mass-richness relation

Authors: K. Kiiveri, D. Gruen, A. Finoguenov, T. Erben, L. van Waerbeke, E. Rykoff, L. Miller, S. Hagstotz, R. Dupke, J. Patrick Henry, J-P. Kneib, G. Gozaliasl, C. C. Kirkpatrick, N. Cibirka, N. Clerc, M. Costanzi, E. S. Cypriano, E. Rozo, H. Shan, P. Spinelli, J. Valiviita, J. Weller

Abstract: The COnstrain Dark Energy with X-ray clusters (CODEX) sample contains the largest flux limited sample of X-ray clusters at $0.35 < z < 0.65$. It was selected from ROSAT data in the 10,000 square degrees of overlap with BOSS, mapping a total number of 2770 high-z galaxy clusters. We present here the full results of the CFHT CODEX program on cluster mass measurement, including a reanalysis of CFHTLS… ▽ More The COnstrain Dark Energy with X-ray clusters (CODEX) sample contains the largest flux limited sample of X-ray clusters at $0.35 < z < 0.65$. It was selected from ROSAT data in the 10,000 square degrees of overlap with BOSS, mapping a total number of 2770 high-z galaxy clusters. We present here the full results of the CFHT CODEX program on cluster mass measurement, including a reanalysis of CFHTLS Wide data, with 25 individual lensing-constrained cluster masses. We employ $lensfit$ shape measurement and perform a conservative colour-space selection and weighting of background galaxies. Using the combination of shape noise and an analytic covariance for intrinsic variations of cluster profiles at fixed mass due to large scale structure, miscentring, and variations in concentration and ellipticity, we determine the likelihood of the observed shear signal as a function of true mass for each cluster. We combine 25 individual cluster mass likelihoods in a Bayesian hierarchical scheme with the inclusion of optical and X-ray selection functions to derive constraints on the slope $α$, normalization $β$, and scatter $σ_{\ln λ| μ}$ of our richness-mass scaling relation model in log-space: $\left<\ln λ| μ\right> = αμ+ β$, with $μ= \ln (M_{200c}/M_{\mathrm{piv}})$, and $M_{\mathrm{piv}} = 10^{14.81} M_{\odot}$. We find a slope $α= 0.49^{+0.20}_{-0.15}$, normalization $ \exp(β) = 84.0^{+9.2}_{-14.8}$ and $σ_{\ln λ| μ} = 0.17^{+0.13}_{-0.09}$ using CFHT richness estimates. In comparison to other weak lensing richness-mass relations, we find the normalization of the richness statistically agreeing with the normalization of other scaling relations from a broad redshift range ($0.0<z<0.65$) and with different cluster selection (X-ray, Sunyaev-Zeldovich, and optical). △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 37 pages, 12 figures

Dark Energy Survey Year 3 Results: Clustering Redshifts -- Calibration of the Weak Lensing Source Redshift Distributions with redMaGiC and BOSS/eBOSS

Authors: M. Gatti, G. Giannini, G. M. Bernstein, A. Alarcon, J. Myles, A. Amon, R. Cawthon, M. Troxel, J. DeRose, S. Everett, A. J. Ross, E. S. Rykoff, J. Elvin-Poole, J. Cordero, I. Harrison, C. Sanchez, J. Prat, D. Gruen, H. Lin, M. Crocce, E. Rozo, T. M. C. Abbott, M. Aguena, S. Allam, J. Annis , et al. (73 additional authors not shown)

Abstract: We present the calibration of the Dark Energy Survey Year 3 (DES Y3) weak lensing source galaxy redshift distributions $n(z)$ from clustering measurements. In particular, we cross-correlate the weak lensing (WL) source galaxies sample with redMaGiC galaxies (luminous red galaxies with secure photometric redshifts) and a spectroscopic sample from BOSS/eBOSS to estimate the redshift distribution of… ▽ More We present the calibration of the Dark Energy Survey Year 3 (DES Y3) weak lensing source galaxy redshift distributions $n(z)$ from clustering measurements. In particular, we cross-correlate the weak lensing (WL) source galaxies sample with redMaGiC galaxies (luminous red galaxies with secure photometric redshifts) and a spectroscopic sample from BOSS/eBOSS to estimate the redshift distribution of the DES sources sample. Two distinct methods for using the clustering statistics are described. The first uses the clustering information independently to estimate the mean redshift of the source galaxies within a redshift window, as done in the DES Y1 analysis. The second method establishes a likelihood of the clustering data as a function of $n(z)$, which can be incorporated into schemes for generating samples of $n(z)$ subject to combined clustering and photometric constraints. Both methods incorporate marginalisation over various astrophysical systematics, including magnification and redshift-dependent galaxy-matter bias. We characterise the uncertainties of the methods in simulations; the first method recovers the mean $z$ of tomographic bins to RMS (precision) of $\sim 0.014$. Use of the second method is shown to vastly improve the accuracy of the shape of $n(z)$ derived from photometric data. The two methods are then applied to the DES Y3 data. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 25 pages, 14 figures. See https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ for the full DES Y3 cosmological releae

Spectroscopic Quantification of Projection Effects in the SDSS redMaPPer Galaxy Cluster Catalogue

Authors: J. Myles, D. Gruen, A. B. Mantz, S. W. Allen, R. G. Morris, E. Rykoff, M. Costanzi, C. To, J. DeRose, R. H. Wechsler, E. Rozo, T. Jeltema, E. R. Carrasco, A. Kremin, R. Kron

Abstract: Projection effects, whereby galaxies along the line-of-sight to a galaxy cluster are mistakenly associated with the cluster halo, present a significant challenge for optical cluster cosmology. We use statistically representative spectral coverage of luminous galaxies to investigate how projection effects impact the low-redshift limit of the Sloan Digital Sky Survey (SDSS) redMaPPer galaxy cluster… ▽ More Projection effects, whereby galaxies along the line-of-sight to a galaxy cluster are mistakenly associated with the cluster halo, present a significant challenge for optical cluster cosmology. We use statistically representative spectral coverage of luminous galaxies to investigate how projection effects impact the low-redshift limit of the Sloan Digital Sky Survey (SDSS) redMaPPer galaxy cluster catalogue. Spectroscopic redshifts enable us to differentiate true cluster members from false positives and determine the fraction of candidate cluster members viewed in projection. Our main results can be summarized as follows: first, we show that a simple double-Gaussian model can be used to describe the distribution of line-of-sight velocities in the redMaPPer sample; second, the incidence of projection effects is substantial, accounting for $\sim 16$ per cent of the weighted richness for the lowest richness objects; third, projection effects are a strong function of richness, with the contribution in the highest richness bin being several times smaller than for low-richness objects; fourth, our measurement has a similar amplitude to state-of-the-art models, but finds a steeper dependence of projection effects on richness than these models; and fifth, the slope of the observed velocity dispersion--richness relation, corrected for projection effects, implies an approximately linear relationship between the true, three-dimensional halo mass and three-dimensional richness. Our results provide a robust, empirical description of the impact of projection effects on the SDSS redMaPPer cluster sample and exemplify the synergies between optical imaging and spectroscopic data for studies of galaxy cluster astrophysics and cosmology. △ Less

Submitted 3 June, 2021; v1 submitted 13 November, 2020; originally announced November 2020.

Comments: 12 pages, 9 figures

Measuring Cosmological Distances Using Cluster Edges as a Standard Ruler

Authors: Erika L. Wagoner, Eduardo Rozo, Han Aung, Daisuke Nagai

Abstract: The line-of-sight velocity dispersion profile of galaxy clusters exhibits a "kink" corresponding to the spatial extent of orbiting galaxies. Because the spatial extent of a cluster is correlated with the amplitude of the velocity dispersion profile, we can utilise this feature as a gravity-calibrated standard ruler. Specifically, the amplitude of the velocity dispersion data allows us to infer the… ▽ More The line-of-sight velocity dispersion profile of galaxy clusters exhibits a "kink" corresponding to the spatial extent of orbiting galaxies. Because the spatial extent of a cluster is correlated with the amplitude of the velocity dispersion profile, we can utilise this feature as a gravity-calibrated standard ruler. Specifically, the amplitude of the velocity dispersion data allows us to infer the physical cluster size. Consequently, observations of the angular scale of the "kink" in the profile can be translated into a distance measurement to the cluster. Assuming the relation between cluster radius and cluster velocity dispersion can be calibrated from simulations, we forecast that with existing data from the Sloan Digital Sky Survey (SDSS) we will be able to measure the Hubble constant with $3\%$ precision. Implementing our method with data from the Dark Energy Spectroscopic Instrument (DESI) will result in a $1.3\%$ measurement of the Hubble constant. Adding cosmological supernova data improves the uncertainty of the DESI measurement to $0.7\%$. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: 8 pages, 5 figures, to be submitted to MNRAS

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances, Weak Lensing, and Galaxy Correlations

Authors: C. To, E. Krause, E. Rozo, H. Wu, D. Gruen, R. H. Wechsler, T. F. Eifler, E. S. Rykoff, M. Costanzi, M. R. Becker, G. M. Bernstein, J. Blazek, S. Bocquet, S. L. Bridle, R. Cawthon, A. Choi, M. Crocce, C. Davis, J. DeRose, A. Drlica-Wagner, J. Elvin-Poole, X. Fang, A. Farahi, O. Friedrich, M. Gatti , et al. (83 additional authors not shown)

Abstract: Combining multiple observational probes is a powerful technique to provide robust and precise constraints on cosmological parameters. In this letter, we present the first joint analysis of cluster abundances and auto/cross correlations of three cosmic tracer fields measured from the first year data of the Dark Energy Survey: galaxy density, weak gravitational lensing shear, and cluster density spl… ▽ More Combining multiple observational probes is a powerful technique to provide robust and precise constraints on cosmological parameters. In this letter, we present the first joint analysis of cluster abundances and auto/cross correlations of three cosmic tracer fields measured from the first year data of the Dark Energy Survey: galaxy density, weak gravitational lensing shear, and cluster density split by optical richness. From a joint analysis of cluster abundances, three cluster cross-correlations, and auto correlations of galaxy density, we obtain $Ω_{\rm{m}}=0.305^{+0.055}_{-0.038}$ and $σ_8=0.783^{+0.064}_{-0.054}$. This result is consistent with constraints from the DES-Y1 galaxy clustering and weak lensing two-point correlation functions for the flat $νΛ$CDM model. We thus combine cluster abundances and all two-point correlations from three cosmic tracer fields and find improved constraints on cosmological parameters as well as on the cluster observable--mass scaling relation. This analysis is an important advance in both optical cluster cosmology and multi-probe analyses of upcoming wide imaging surveys. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 14 pages, 7 figures, 2 tables; submitted to PRL; comments are welcome

Journal ref: Phys. Rev. Lett. 126, 141301 (2021)

The Atacama Cosmology Telescope: A Catalog of > 4000 Sunyaev-Zel'dovich Galaxy Clusters

Authors: M. Hilton, C. Sifón, S. Naess, M. Madhavacheril, M. Oguri, E. Rozo, E. Rykoff, T. M. C. Abbott, S. Adhikari, M. Aguena, S. Aiola, S. Allam, S. Amodeo, A. Amon, J. Annis, B. Ansarinejad, C. Aros-Bunster, J. E. Austermann, S. Avila, D. Bacon, N. Battaglia, J. A. Beall, D. T. Becker, G. M. Bernstein, E. Bertin , et al. (124 additional authors not shown)

Abstract: We present a catalog of 4195 optically confirmed Sunyaev-Zel'dovich (SZ) selected galaxy clusters detected with signal-to-noise > 4 in 13,211 deg$^2$ of sky surveyed by the Atacama Cosmology Telescope (ACT). Cluster candidates were selected by applying a multi-frequency matched filter to 98 and 150 GHz maps constructed from ACT observations obtained from 2008-2018, and confirmed using deep, wide-a… ▽ More We present a catalog of 4195 optically confirmed Sunyaev-Zel'dovich (SZ) selected galaxy clusters detected with signal-to-noise > 4 in 13,211 deg$^2$ of sky surveyed by the Atacama Cosmology Telescope (ACT). Cluster candidates were selected by applying a multi-frequency matched filter to 98 and 150 GHz maps constructed from ACT observations obtained from 2008-2018, and confirmed using deep, wide-area optical surveys. The clusters span the redshift range 0.04 < z < 1.91 (median z = 0.52). The catalog contains 222 z > 1 clusters, and a total of 868 systems are new discoveries. Assuming an SZ-signal vs. mass scaling relation calibrated from X-ray observations, the sample has a 90% completeness mass limit of M500c > 3.8 x 10$^{14}$ MSun, evaluated at z = 0.5, for clusters detected at signal-to-noise ratio > 5 in maps filtered at an angular scale of 2.4'. The survey has a large overlap with deep optical weak-lensing surveys that are being used to calibrate the SZ-signal mass-scaling relation, such as the Dark Energy Survey (4566 deg$^2$), the Hyper Suprime-Cam Subaru Strategic Program (469 deg$^2$), and the Kilo Degree Survey (825 deg$^2$). We highlight some noteworthy objects in the sample, including potentially projected systems; clusters with strong lensing features; clusters with active central galaxies or star formation; and systems of multiple clusters that may be physically associated. The cluster catalog will be a useful resource for future cosmological analyses, and studying the evolution of the intracluster medium and galaxies in massive clusters over the past 10 Gyr. △ Less

Submitted 2 December, 2020; v1 submitted 23 September, 2020; originally announced September 2020.

Comments: 35 pages, 27 figures, accepted for publication in ApJS; v1.0 catalogs will be available from LAMBDA https://lambda.gsfc.nasa.gov/product/act/actpol_prod_table.cfm; v1.0 catalogs available from https://astro.ukzn.ac.za/~mjh/ACTDR5/v1.0/ until then

Linear Systematics Mitigation in Galaxy Clustering in the Dark Energy Survey Year 1 Data

Authors: Erika L. Wagoner, Eduardo Rozo, Xiao Fang, Martín Crocce, Jack Elvin-Poole, Noah Weaverdyck

Abstract: We implement a linear model for mitigating the effect of observing conditions and other sources of contamination in galaxy clustering analyses. Our treatment improves upon the fiducial systematics treatment of the Dark Energy Survey (DES) Year 1 (Y1) cosmology analysis in four crucial ways. Specifically, our treatment: 1) does not require decisions as to which observable systematics are significan… ▽ More We implement a linear model for mitigating the effect of observing conditions and other sources of contamination in galaxy clustering analyses. Our treatment improves upon the fiducial systematics treatment of the Dark Energy Survey (DES) Year 1 (Y1) cosmology analysis in four crucial ways. Specifically, our treatment: 1) does not require decisions as to which observable systematics are significant and which are not, allowing for the possibility of multiple maps adding coherently to give rise to significant bias even if no single map leads to a significant bias by itself; 2) characterizes both the statistical and systematic uncertainty in our mitigation procedure, allowing us to propagate said uncertainties into the reported cosmological constraints; 3) explicitly exploits the full spatial structure of the galaxy density field to differentiate between cosmology-sourced and systematics-sourced fluctuations within the galaxy density field; 4) is fully automated, and can therefore be trivially applied to any data set. The updated correlation function for the DES Y1 redMaGiC catalog minimally impacts the cosmological posteriors from that analysis. Encouragingly, our analysis does improve the goodness of fit statistic of the DES Y1 3$\times$2pt data set ($Δχ^2 = -6.5$ with no additional parameters). This improvement is due in nearly equal parts to both the change in the correlation function and the added statistical and systematic uncertainties associated with our method. We expect the difference in mitigation techniques to become more important in future work as the size of cosmological data sets grows. △ Less

Submitted 12 April, 2021; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 15 pages, 8 figures, submitted to MNRAS. Updated to be consistent with accepted version and include updated citations

Probing galaxy evolution in massive clusters using ACT and DES: splashback as a cosmic clock

Authors: Susmita Adhikari, Tae-hyeon Shin, Bhuvnesh Jain, Matt Hilton, Eric Baxter, Chihway Chang, Risa H. Wechsler, Nick Battaglia, J. Richard Bond, Sebastian Bocquet, Joseph DeRose, Steve K. Choi, Mark Devlin, Jo Dunkley, August E. Evrard, Simone Ferraro, J. Colin Hill, John P. Hughes, Patricio A. Gallardo, Martine Lokken, Amanda MacInnis, Jeffrey McMahon, Mathew S. Madhavacheril, Frederico Nati, Laura B. Newburgh , et al. (91 additional authors not shown)

Abstract: We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is loc… ▽ More We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is located at $2.4^{+0.3}_{-0.4}$ Mpc $h^{-1}$. We split the sample based on galaxy color and find significant differences in the profile shapes. Red galaxies and those in the green valley show a splashback-like minimum in their slope profile consistent with theoretical predictions, while the bluest galaxies show a weak feature that appears at a smaller radius. We develop a mapping of galaxies to subhalos in $N$-body simulations by splitting subhalos based on infall time onto the cluster halos. We find that the location of the steepest slope and differences in the shapes of the profiles can be mapped to differences in the average time of infall of galaxies of different colors. The minima of the slope in the galaxy profiles trace a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time for galaxies of different colours, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in their clusters for over $3.2 ~\rm Gyrs$, green galaxies about $2.2 ~\rm Gyrs$, while blue galaxies have been accreted most recently and have not reached apocenter. Using the information from the complete radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching in clusters occurs after a delay of about a gigayear, and that galaxies quench rapidly thereafter with an exponential timescale of $0.6$ Gyr. △ Less

Submitted 26 August, 2020; originally announced August 2020.

Comments: 23 pages, 10 figures, to be submitted to ApJ

Combination of cluster number counts and two-point correlations: Validation on Mock Dark Energy Survey

Authors: C. To, E. Krause, E. Rozo, H. Wu, D. Gruen, J. DeRose, E. S. Rykoff, R. H. Wechsler, M. R. Becker, M. Costanzi, T. Eifler, M. E. S. Pereira, N. Kokron

Abstract: We present a method of combining cluster abundances and large-scale two-point correlations, namely galaxy clustering, galaxy--cluster cross-correlations, cluster auto-correlations, and cluster lensing. This data vector yields comparable cosmological constraints to traditional analyses that rely on small-scale cluster lensing for mass calibration. We use cosmological survey simulations designed to… ▽ More We present a method of combining cluster abundances and large-scale two-point correlations, namely galaxy clustering, galaxy--cluster cross-correlations, cluster auto-correlations, and cluster lensing. This data vector yields comparable cosmological constraints to traditional analyses that rely on small-scale cluster lensing for mass calibration. We use cosmological survey simulations designed to resemble the Dark Energy Survey Year One (DES-Y1) data to validate the analytical covariance matrix and the parameter inferences. The posterior distribution from the analysis of simulations is statistically consistent with the absence of systematic biases detectable at the precision of the DES Y1 experiment. We compare the $χ^2$ values in simulations to their expectation and find no significant difference. The robustness of our results against a variety of systematic effects is verified using a simulated likelihood analysis of a Dark Energy Survey Year 1-like data vectors. This work presents the first-ever end-to-end validation of a cluster abundance cosmological analysis on galaxy catalog-level simulations. △ Less

Submitted 24 January, 2021; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: 21 pages, 6 figures, 3 tables. MNRAS accepted

Dark Energy Survey Year 1 Results: Constraining Baryonic Physics in the Universe

Authors: Hung-Jin Huang, Tim Eifler, Rachel Mandelbaum, Gary M. Bernstein, Anqi Chen, Ami Choi, Juan García-Bellido, Dragan Huterer, Elisabeth Krause, Eduardo Rozo, Sukhdeep Singh, Sarah Bridle, Joseph DeRose, Jack Elvin-Pole, Xiao Fang, Oliver Friedrich, Marco Gatti, Enrique Gaztanaga, Daniel Gruen, Will Hartley, Ben Hoyle, Mike Jarvis, Niall MacCrann, Markus Rau, Vivian Miranda , et al. (52 additional authors not shown)

Abstract: Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3$\times$2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic… ▽ More Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3$\times$2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modeling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude ($Q_1$) generally reflects the strength of baryon feedback. With the upper limit of $Q_1$ prior being bound by the Illustris feedback scenarios, we reach $\sim 20\%$ improvement in the constraint of $S_8=σ_8(Ω_{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3$\times$2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative $Q_1$ prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than $2 σ$. △ Less

Submitted 26 April, 2021; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: 22 pages, 18 figures, 2 tables. accepted to MNRAS. A brief video summary of this paper is available at https://www.youtube.com/watch?v=QbeNwk5papU

A Redefinition of the Halo Boundary Leads to a Simple yet Accurate Halo Model of Large Scale Structure

Authors: Rafael Garcia, Eduardo Rozo, Matthew R. Becker, Surhud More

Abstract: We present a model for the halo--mass correlation function that explicitly incorporates halo exclusion. We assume that halos trace mass in a way that can be described using a single scale-independent bias parameter. However, our model exhibits scale dependent biasing due to the impact of halo-exclusion, the use of a ``soft'' (i.e. not infinitely sharp) halo boundary, and differences in the one hal… ▽ More We present a model for the halo--mass correlation function that explicitly incorporates halo exclusion. We assume that halos trace mass in a way that can be described using a single scale-independent bias parameter. However, our model exhibits scale dependent biasing due to the impact of halo-exclusion, the use of a ``soft'' (i.e. not infinitely sharp) halo boundary, and differences in the one halo term contributions to $ξ_{\rm hm}$ and $ξ_{\rm mm}$. These features naturally lead us to a redefinition of the halo boundary that lies at the ``by eye'' transition radius from the one--halo to the two--halo term in the halo--mass correlation function. When adopting our proposed definition, our model succeeds in describing the halo--mass correlation function with $\approx 2\%$ residuals over the radial range $0.1\ h^{-1}{\rm Mpc} < r < 80\ h^{-1}{\rm Mpc}$, and for halo masses in the range $10^{13}\ h^{-1}{\rm M}_{\odot} < M < 10^{15}\ h^{-1}{\rm M}_{\odot}$. Our proposed halo boundary is related to the splashback radius by a roughly constant multiplicative factor. Taking the 87-percentile as reference we find $r_{\rm t}/R_{\rm sp} \approx 1.3$. Surprisingly, our proposed definition results in halo abundances that are well described by the Press-Schechter mass function with $δ_{\rm sc}=1.449\pm 0.004$. The clustering bias parameter is offset from the standard background-split prediction by $\approx 10\%-15\%$. This level of agreement is comparable to that achieved with more standard halo definitions. △ Less

Submitted 23 June, 2020; originally announced June 2020.

Comments: 12 pages, 5 figures

Localizing Transformations of the Galaxy-Galaxy Lensing Observable

Authors: Youngsoo Park, Eduardo Rozo, Elisabeth Krause

Abstract: Modern cosmological analyses of galaxy-galaxy lensing face a theoretical systematic effect arising from the non-locality of the observed galaxy-galaxy lensing signal. Because the predicted tangential shear signal at a given separation depends on the physical modeling on all scales internal to that separation, systematic uncertainties in the modeling of non-linear small scales are propagated outwar… ▽ More Modern cosmological analyses of galaxy-galaxy lensing face a theoretical systematic effect arising from the non-locality of the observed galaxy-galaxy lensing signal. Because the predicted tangential shear signal at a given separation depends on the physical modeling on all scales internal to that separation, systematic uncertainties in the modeling of non-linear small scales are propagated outwards to larger scales. Even in the absence of other limiting factors, this systematic effect alone can necessitate conservative small-scale cuts, resulting in significant losses of information in the tangential shear data vector. We construct a simple linear transformation of the observable that removes this non-locality, enabling more aggressive small-scale cuts for a given theoretical model. Our modified galaxy-galaxy lensing observable makes it possible to include observations on significantly smaller scales than those under the standard approach in cosmological analyses. More importantly, it ensures that the cosmological signal contained within the observable is exclusively drawn from well-understood physical scales. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 126, 021301 (2021)

Cosmology with the Wide-Field Infrared Survey Telescope -- Multi-Probe Strategies

Authors: Tim Eifler, Hironao Miyatake, Elisabeth Krause, Chen Heinrich, Vivian Miranda, Christopher Hirata, Jiachuan Xu, Shoubaneh Hemmati, Melanie Simet, Peter Capak, Ami Choi, Olivier Dore, Cyrille Doux, Xiao Fang, Rebekah Hounsell, Eric Huff, Hung-Jin Huang, Mike Jarvis, Dan Masters, Eduardo Rozo, Dan Scolnic, David N. Spergel, Michael Troxel, Anja von der Linden, Yun Wang , et al. (3 additional authors not shown)

Abstract: We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photome… ▽ More We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST. △ Less

Submitted 10 April, 2020; originally announced April 2020.

Comments: comments welcome

Cosmology with the Wide-Field Infrared Survey Telescope -- Synergies with the Rubin Observatory Legacy Survey of Space and Time

Authors: Tim Eifler, Melanie Simet, Elisabeth Krause, Christopher Hirata, Hung-Jin Huang, Xiao Fang, Vivian Miranda, Rachel Mandelbaum, Cyrille Doux, Chen Heinrich, Eric Huff, Hironao Miyatake, Shoubaneh Hemmati, Jiachuan Xu, Paul Rogozenski, Peter Capak, Ami Choi, Olivier Dore, Bhuvnesh Jain, Mike Jarvis, Niall MacCrann, Dan Masters, Eduardo Rozo, David N. Spergel, Michael Troxel , et al. (5 additional authors not shown)

Abstract: We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combin… ▽ More We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)). △ Less

Submitted 11 April, 2020; v1 submitted 9 April, 2020; originally announced April 2020.

Comments: added references, comments welcome

The Phase Space Structure of Dark Matter Halos

Authors: Han Aung, Daisuke Nagai, Eduardo Rozo, Rafael Garcia

Abstract: The phase space structure of dark matter halos can be used to measure the mass of the halo, infer mass accretion rates, and probe the effects of modified gravity. Previous studies showed that the splashback radius can be measured in position space using the slope of the density profile. Using N-body simulations, we show that the phase space structure of the dark matter halo does not end at this sp… ▽ More The phase space structure of dark matter halos can be used to measure the mass of the halo, infer mass accretion rates, and probe the effects of modified gravity. Previous studies showed that the splashback radius can be measured in position space using the slope of the density profile. Using N-body simulations, we show that the phase space structure of the dark matter halo does not end at this splashback radius. Instead, there exists a region where infalling, splashback, and virialized halos are mixed spatially. We model the distribution of the three kinematically distinct populations and show that there exists an "edge radius" beyond which a dark matter halo has no orbiting substructures. This radius is a fixed multiple of the splashback radius as defined in previous works, and can be interpreted as a radius which contains a fixed fraction of the apocenters of dark matter particles. Our results provide a firm theoretical foundation to the satellite galaxy model adopted in the companion paper by Tomooka et al., where we analyzed the phase space distribution of SDSS redMaPPer clusters. △ Less

Submitted 19 February, 2021; v1 submitted 25 March, 2020; originally announced March 2020.

Comments: 7 pages, 5 figures, accepted in MNRAS, companion paper to Tomooka et al. 2020

Journal ref: MNRAS, 502, 1041 (2021)

Clusters Have Edges: The Projected Phase SpaceStructure of SDSS redMaPPer Clusters

Authors: Paxton Tomooka, Eduardo Rozo, Erika L. Wagoner, Han Aung, Daisuke Nagai, Sasha Safonova

Abstract: We study the distribution of line-of-sight velocities of galaxies in the vicinity of SDSS redMaPPer galaxy clusters. Based on their velocities, galaxies can be split into two categories: galaxies that are dynamically associated with the cluster, and random line-of-sight projections. Both the fraction of galaxies associated with the galaxy clusters, and the velocity dispersion of the same, exhibit… ▽ More We study the distribution of line-of-sight velocities of galaxies in the vicinity of SDSS redMaPPer galaxy clusters. Based on their velocities, galaxies can be split into two categories: galaxies that are dynamically associated with the cluster, and random line-of-sight projections. Both the fraction of galaxies associated with the galaxy clusters, and the velocity dispersion of the same, exhibit a sharp feature as a function of radius. The feature occurs at a radial scale $R_{\rm edge} \approx 2.2R_{\rmλ}$, where $R_{\rmλ}$ is the cluster radius assigned by redMaPPer. We refer to $R_{\rm edge}$ as the "edge radius." These results are naturally explained by a model that further splits the galaxies dynamically associated with a galaxy cluster into a component of galaxies orbiting the halo and an infalling galaxy component. The edge radius $R_{\rm edge}$ constitutes a true "cluster edge", in the sense that no orbiting structures exist past this radius. A companion paper (Aung et al. 2020) tests whether the "halo edge" hypothesis holds when investigating the full three-dimensional phase space distribution of dark matter substructures in numerical simulations, and demonstrates that this radius coincides with a suitably defined splashback radius. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: 10 pages, 5 figures, submitted to MNRAS, companion paper to Aung et al. 2020

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

Authors: DES Collaboration, Tim Abbott, Michel Aguena, Alex Alarcon, Sahar Allam, Steve Allen, James Annis, Santiago Avila, David Bacon, Alberto Bermeo, Gary Bernstein, Emmanuel Bertin, Sunayana Bhargava, Sebastian Bocquet, David Brooks, Dillon Brout, Elizabeth Buckley-Geer, David Burke, Aurelio Carnero Rosell, Matias Carrasco Kind, Jorge Carretero, Francisco Javier Castander, Ross Cawthon, Chihway Chang, Xinyi Chen , et al. (107 additional authors not shown)

Abstract: 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 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter densi… ▽ More 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 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ 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 ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 35 pages, 20 figures, submitted to Physical Review D

Journal ref: Phys. Rev. D 102, 023509 (2020)

Milky Way Satellite Census. I. The Observational Selection Function for Milky Way Satellites in DES Y3 and Pan-STARRS DR1

Authors: A. Drlica-Wagner, K. Bechtol, S. Mau, M. McNanna, E. O. Nadler, A. B. Pace, T. S. Li, A. Pieres, E. Rozo, J. D. Simon, A. R. Walker, R. H. Wechsler, T. M. C. Abbott, S. Allam, J. Annis, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, M. Costanzi, L. N. da Costa, J. De Vicente, S. Desai , et al. (33 additional authors not shown)

Abstract: We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ~80% of the sky. Our search for satellite galaxies targets ~25,000 deg$^2$ of the high-Galactic-latitude sky reaching a 10$σ$ point-source depth of… ▽ More We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ~80% of the sky. Our search for satellite galaxies targets ~25,000 deg$^2$ of the high-Galactic-latitude sky reaching a 10$σ$ point-source depth of $\gtrsim$ 22.5 mag in the $g$ and $r$ bands. While satellite galaxy searches have been performed independently on DES and PS1 before, this is the first time that a self-consistent search is performed across both data sets. We do not detect any new high-significance satellite galaxy candidates, while recovering the majority of satellites previously detected in surveys of comparable depth. We characterize the sensitivity of our search using a large set of simulated satellites injected into the survey data. We use these simulations to derive both analytic and machine-learning models that accurately predict the detectability of Milky Way satellites as a function of their distance, size, luminosity, and location on the sky. To demonstrate the utility of this observational selection function, we calculate the luminosity function of Milky Way satellite galaxies, assuming that the known population of satellite galaxies is representative of the underlying distribution. We provide access to our observational selection function to facilitate comparisons with cosmological models of galaxy formation and evolution. △ Less

Submitted 16 April, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 34 pages, 11 figures, 6 tables; updated to match published version. Selection functions available at: https://github.com/des-science/mw-sats

Report number: FERMILAB-PUB-19-604-AE

Journal ref: ApJ 893, 1 (2020)

CODEX clusters. The Survey, the Catalog, and Cosmology of the X-ray Luminosity Function

Authors: A. Finoguenov, E. Rykoff, N. Clerc, M. Costanzi, S. Hagstotz, J. Ider Chitham, K. Kiiveri, C. C. Kirkpatrick, R. Capasso, J. Comparat, S. Damsted, R. Dupke, G. Erfanianfar, J. Patrick Henry, F. Kaefer, J-P. Kneib, V. Lindholm, E. Rozo, L. van Waerbeke, J. Weller

Abstract: Large area catalogs of galaxy clusters constructed from ROSAT All Sky Survey provide the base for our knowledge on the population of clusters thanks to the long-term multiwavelength efforts on their follow-up. Advent of large area photometric surveys superseding in depth previous all-sky data allows us to revisit the construction of X-ray cluster catalogs, extending the study to lower cluster mass… ▽ More Large area catalogs of galaxy clusters constructed from ROSAT All Sky Survey provide the base for our knowledge on the population of clusters thanks to the long-term multiwavelength efforts on their follow-up. Advent of large area photometric surveys superseding in depth previous all-sky data allows us to revisit the construction of X-ray cluster catalogs, extending the study to lower cluster masses and to higher redshifts and to provide the modelling of the selection function. We perform a wavelet detection of X-ray sources and make extensive simulations of the detection of clusters in the RASS data. We assign an optical richness to each of the 24,788 detected X-ray sources in the 10,382 square degrees of SDSS BOSS area, using redMaPPer version 5.2. We name this survey COnstrain Dark Energy with X-ray (CODEX) clusters. We show that there is no obvious separation of sources on galaxy clusters and AGN, based on distribution of systems on their richness. This is a combination of increasing number of galaxy groups and their selection as identification of an X-ray sources either by chance or due to groups hosting an AGN. To clean the sample, we use a cut on the optical richness at the level corresponding to the 10\% completeness of the survey and include it into the modelling of cluster selection function. We present the X-ray catalog extending to a redshift of 0.6 down to X-ray fluxes of $10^{-13}$ ergs s$^{-1}$ cm$^{-2}$. We provide the modelling of the sample selection and discuss the redshift evolution of the high end of the X-ray luminosity function (XLF). Our results on $z<0.3$ XLF are in agreement with previous studies, while we provide new constraints on the $0.3<z<0.6$ XLF. We find a lack of strong redshift evolution of the XLF, provide exact modeling of the effect of low number statistics and AGN contamination and present the resulting constraints on the flat $Λ$CDM. △ Less

Submitted 29 April, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 10 pages, 11 figures, A&A in press, part of SDSS DR16 data release

Journal ref: A&A 638, A114 (2020)

Improved Cosmological Constraints from SDSS redMaPPer Clusters via X-ray Follow-up of a Complete Subsample of Systems

Authors: Matthew Kirby, Eduardo Rozo, R. Glenn Morris, Steven W. Allen, Matteo Costanzi, Tesla E. Jeltema, Adam B. Mantz, A. Kathy Romer, E. S. Rykoff, Anja von der Linden

Abstract: We improve upon the cosmological constraints derived from the abundance and weak-lensing data of redMaPPer clusters detected in the Sloan Digital Sky Survey (SDSS). Specifically, we derive gas mass data using Chandra X-ray follow-up of a complete sample of the 30 richest SDSS redMaPPer clusters with $z\in[0.1,0.3]$, and use these additional data to improve upon the original analysis by Costanzi et… ▽ More We improve upon the cosmological constraints derived from the abundance and weak-lensing data of redMaPPer clusters detected in the Sloan Digital Sky Survey (SDSS). Specifically, we derive gas mass data using Chandra X-ray follow-up of a complete sample of the 30 richest SDSS redMaPPer clusters with $z\in[0.1,0.3]$, and use these additional data to improve upon the original analysis by Costanzi et al. (2019b). We simultaneously fit for the parameters of the richness-mass relation, the cluster gas mass-mass relation, and cosmology. By including our X-ray cluster sample in the SDSS cluster cosmology analysis, we measure $Ω_{\rm m} = 0.25 \pm 0.04$ and $σ_8 = 0.85^{+0.06}_{-0.08}$. These constraints represent a 25.5% and 29.8% reduction in the size of the 68% confidence intervals of $Ω_{\rm m}$ and $σ_8$ respectively, relative to the constraints published in Costanzi et al. (2019b). Our cosmological constraints are in agreement with early universe results from Planck. As a byproduct of our analysis, we also perform an independent calibration of the amplitude of the $\langle M_{\rm gas}^{\rm true}|M_{\rm 500c}\rangle$ scaling relation. Our calibration is consistent with and of comparable precision to that of Mantz et al. (2016b). △ Less

Submitted 14 November, 2019; v1 submitted 29 October, 2019; originally announced October 2019.

RedMaPPer: Evolution and Mass Dependence of the Conditional Luminosity Functions of Red Galaxies in Galaxy Clusters

Authors: Chun-Hao To, Rachel M. Reddick, Eduardo Rozo, Eli Rykoff, Risa H. Wechsler

Abstract: We characterize the luminosity distribution, halo mass dependence, and redshift evolution of red galaxies in galaxy clusters using the SDSS Data Release 8 RedMaPPer cluster sample. We propose a simple prescription for the relationship between the luminosity of both central and satellite galaxies and the mass of their host halos, and show that this model is well-fit by the data. Using a larger gala… ▽ More We characterize the luminosity distribution, halo mass dependence, and redshift evolution of red galaxies in galaxy clusters using the SDSS Data Release 8 RedMaPPer cluster sample. We propose a simple prescription for the relationship between the luminosity of both central and satellite galaxies and the mass of their host halos, and show that this model is well-fit by the data. Using a larger galaxy cluster sample than previously employed in the literature, we find that the luminosities of central galaxies scale as $\langle \log L \rangle \propto A_L \log (M_{200b})$, with $A_L=0.39\pm0.04$, and that the scatter of the central--galaxy luminosity at fixed $M_{200b}$ ( $σ_{\log L|M}$) is $0.23 ^{+0.05}_{-0.04}$ dex, with the error bar including systematics due to miscentering of the cluster finder, photometry, and photometric redshift estimation. Our data prefers a positive correlation between the luminosity of central galaxies and the observed richness of clusters at a fixed halo mass, with an effective correlation coefficient $d_{\rm{eff}}=0.36^{+0.17}_{-0.16}$. The characteristic luminosity of satellites becomes dimmer from $z=0.3$ to $z=0.1$ by $\sim 20\%$ after accounting for passive evolution. We estimate the fraction of galaxy clusters where the brightest galaxy is not the central to be $P_{\rm{BNC}} \sim 20\%$. We discuss implications of these findings in the context of galaxy evolution and the galaxy--halo connection. △ Less

Submitted 29 October, 2020; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: 29 pages, 17 figures, 5 tables. Accepted by APJ

Journal ref: APJ, 897, 15, 2020

The Aemulus Project IV: Emulating Halo Bias

Authors: Thomas McClintock, Eduardo Rozo, Arka Banerjee, Matthew R. Becker, Joseph DeRose, Sean McLaughlin, Jeremy L. Tinker, Risa H. Wechsler, Zhongxu Zhai

Abstract: Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substanti… ▽ More Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substantial improvements in accuracy compared to the widely used Tinker et al. (2010) model. Halos in this work are defined using an overdensity criteria of 200 relative to the mean background density. Halo catalogs are produced for 40 N-body simulations as part of the Aemulus project at snapshots from z=3 to z=0. The emulator is trained over the mass range $6\times10^{12}-7\times10^{15}\ h^{-1}M_{\odot}$. Using an additional suite of 35 simulations, we determine that the precision of the emulator is redshift dependent, achieving sub-percent levels for a majority of the redshift range. Two additional simulation suites are used to test the ability of the emulator to extrapolate to higher and lower masses. Our high-resolution simulation suite is used to develop an extrapolation scheme in which the emulator asymptotes to the Tinker et al. (2010) model at low mass, achieving ~3% accuracy down to $10^{11}\ h^{-1}M_{\odot}$. Finally, we present a method to propagate emulator modeling uncertainty into an error budget. Our emulator is made publicly available at \url{https://github.com/AemulusProject/bias_emulator}. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 14 pages, 7 figures, Submitted to ApJ

A Detection of CMB-Cluster Lensing using Polarization Data from SPTpol

Authors: S. Raghunathan, S. Patil, E. Baxter, B. A. Benson, L. E. Bleem, T. M. Crawford, G. P. Holder, T. McClintock, C. L. Reichardt, T. N. Varga, N. Whitehorn, P. A. R. Ade, S. Allam, A. J. Anderson, J. E. Austermann, S. Avila, J. S. Avva, D. Bacon, J. A. Beall, A. N. Bender, F. Bianchini, S. Bocquet, D. Brooks, D. L. Burke, J. E. Carlstrom , et al. (104 additional authors not shown)

Abstract: We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization… ▽ More We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $λ\ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8σ$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements. △ Less

Submitted 24 September, 2019; v1 submitted 19 July, 2019; originally announced July 2019.

Comments: 10 pages, 3 figures, 1 table; typos fixed; accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 123, 181301 (2019)

Concordance Cosmology?

Authors: Youngsoo Park, Eduardo Rozo

Abstract: We propose a new intuitive metric for evaluating the tension between two experiments, and apply it to several data sets. While our metric is non-optimal, if evidence of tension is detected, this evidence is robust and easy to interpret. Assuming a flat $Λ$CDM cosmological model, we find that there is a modest $2.2σ$ tension between the DES Year 1 results and the ${\it Planck}$ measurements of the… ▽ More We propose a new intuitive metric for evaluating the tension between two experiments, and apply it to several data sets. While our metric is non-optimal, if evidence of tension is detected, this evidence is robust and easy to interpret. Assuming a flat $Λ$CDM cosmological model, we find that there is a modest $2.2σ$ tension between the DES Year 1 results and the ${\it Planck}$ measurements of the Cosmic Microwave Background (CMB). This tension is driven by the difference between the amount of structure observed in the late-time Universe and that predicted from fitting the ${\it Planck}$ data, and appears to be unrelated to the tension between ${\it Planck}$ and local esitmates of the Hubble rate. In particular, combining DES, Baryon Acoustic Oscillations (BAO), Big-Bang Nucleosynthesis (BBN), and supernovae (SNe) measurements recovers a Hubble constant and sound horizon consistent with ${\it Planck}$, and in tension with local distance-ladder measurements. If the tension between these various data sets persists, it is likely that reconciling ${\it all}$ current data will require breaking the flat $Λ$CDM model in at least two different ways: one involving new physics in the early Universe, and one involving new late-time Universe physics. △ Less

Submitted 12 July, 2019; originally announced July 2019.

Comments: 8 pages. 5 figures

Dark Energy Survey Year 1 Results: Wide field mass maps via forward fitting in harmonic space

Authors: B. Mawdsley, D. Bacon, C. Chang, P. Melchior, E. Rozo, S. Seitz, N. Jeffrey, M. Gatti, E. Gaztanaga, D. Gruen, W. G. Hartley, B. Hoyle, S. Samuroff, E. Sheldon, M. A. Troxel, J. Zuntz, T. M. C. Abbott, J. Annis, E. Bertin, S. L. Bridle, D. Brooks, E. Buckley-Geer, D. L. Burke, A. Carnero Rosell, M. Car- rasco Kind , et al. (41 additional authors not shown)

Abstract: We present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser-Squires method, which applies a dire… ▽ More We present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser-Squires method, which applies a direct inversion on the data; our approach is in good agreement with the previous direct approach in the central regions of the footprint. The mapping technique is assessed and verified with tests on simulations; together with the Kaiser-Squires method, the technique is then applied to data from the Dark Energy Survey Year 1 data and the differences between the two methods are compared. We also produce the first DES measurements of the convergence Minkowski functionals and compare them to those measured in simulations. △ Less

Submitted 29 May, 2019; originally announced May 2019.

Comments: 19 pages, 16 figures. Submitted to MNRAS

Reconstructing Probability Distributions with Gaussian Processes

Authors: Thomas McClintock, Eduardo Rozo

Abstract: Modern cosmological analyses constrain physical parameters using Markov Chain Monte Carlo (MCMC) or similar sampling techniques. Oftentimes, these techniques are computationally expensive to run and require up to thousands of CPU hours to complete. Here we present a method for reconstructing the log-probability distributions of completed experiments from an existing MCMC chain (or any set of poste… ▽ More Modern cosmological analyses constrain physical parameters using Markov Chain Monte Carlo (MCMC) or similar sampling techniques. Oftentimes, these techniques are computationally expensive to run and require up to thousands of CPU hours to complete. Here we present a method for reconstructing the log-probability distributions of completed experiments from an existing MCMC chain (or any set of posterior samples). The reconstruction is performed using Gaussian process regression for interpolating the log-probability. This allows for easy resampling, importance sampling, marginalization, testing different samplers, investigating chain convergence, and other operations. As an example use-case, we reconstruct the posterior distribution of the most recent Planck 2018 analysis. We then resample the posterior, and generate a new MCMC chain with forty times as many points in only thirty minutes. Our likelihood reconstruction tool can be found online at https://github.com/tmcclintock/AReconstructionTool. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Comments: 7 pages, 4 figures, repository at https://github.com/tmcclintock/AReconstructionTool