Showing 1–50 of 70 results for author: Farahi, A

Introducing the DREAMS Project: DaRk mattEr and Astrophysics with Machine learning and Simulations

Authors: Jonah C. Rose, Paul Torrey, Francisco Villaescusa-Navarro, Mariangela Lisanti, Tri Nguyen, Sandip Roy, Kassidy E. Kollmann, Mark Vogelsberger, Francis-Yan Cyr-Racine, Mikhail V. Medvedev, Shy Genel, Daniel Anglés-Alcázar, Nitya Kallivayalil, Bonny Y. Wang, Belén Costanza, Stephanie O'Neil, Cian Roche, Soumyodipta Karmakar, Alex M. Garcia, Ryan Low, Shurui Lin, Olivia Mostow, Akaxia Cruz, Andrea Caputo, Arya Farahi , et al. (5 additional authors not shown)

Abstract: We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- f… ▽ More We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of Warm Dark Matter, each comprised of 1024 simulations generated using the Arepo code. One suite consists of uniform-box simulations covering a $(25~h^{-1}~{\rm M}_\odot)^3$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the Warm Dark Matter particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and Convolutional Neural Networks, to demonstrate how such simulation suites can be used to disentangle the effects of dark matter and baryonic physics on galactic properties. The DREAMS project can be extended further to include different dark matter models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of dark matter on these scales. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 28 pages, 8 figures, DREAMS website: https://www.dreams-project.org

Can we constrain warm dark matter masses with individual galaxies?

Authors: Shurui Lin, Francisco Villaescusa-Navarro, Jonah Rose, Paul Torrey, Arya Farahi, Kassidy E. Kollmann, Alex M. Garcia, Sandip Roy, Nitya Kallivayalil, Mark Vogelsberger, Yi-Fu Cai, Wentao Luo

Abstract: We study the impact of warm dark matter mass on the internal properties of individual galaxies using a large suite of 1,024 state-of-the-art cosmological hydrodynamic simulations from the DREAMS project. We take individual galaxies' properties from the simulations, which have different cosmologies, astrophysics, and warm dark matter masses, and train normalizing flows to learn the posterior of the… ▽ More We study the impact of warm dark matter mass on the internal properties of individual galaxies using a large suite of 1,024 state-of-the-art cosmological hydrodynamic simulations from the DREAMS project. We take individual galaxies' properties from the simulations, which have different cosmologies, astrophysics, and warm dark matter masses, and train normalizing flows to learn the posterior of the parameters. We find that our models cannot infer the value of the warm dark matter mass, even when the values of the cosmological and astrophysical parameters are given explicitly. This result holds for galaxies with stellar mass larger than $2\times10^8 M_\odot/h$ at both low and high redshifts. We calculate the mutual information and find no significant dependence between the WDM mass and galaxy properties. On the other hand, our models can infer the value of $Ω_{\rm m}$ with a $\sim10\%$ accuracy from the properties of individual galaxies while marginalizing astrophysics and warm dark matter masses. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 13 pages, 8 figures

Copacabana: A Probabilistic Membership Assignment Method for Galaxy Clusters

Authors: J. H. Esteves, M. E. S. Pereira, M. Soares-Santos, J. Annis, A. Farahi, F. Andrade-Oliveira, P. Barchi, A. Palmese, H. Lin, B. Welch, H. -Y. Wu, M. Aguena, O. Alves D. Bacon, S. Bocquet, D. Brooks, A. Carnero Rosell, J. Carretero, M. Costanzi, L. N. da Costa, J. De Vicente, P. Doel, S. Everett, B. Flaugher, J. Frieman, J. García-Bellido , et al. (30 additional authors not shown)

Abstract: Cosmological analyses using galaxy clusters in optical/NIR photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for {\it all} galaxies within a… ▽ More Cosmological analyses using galaxy clusters in optical/NIR photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for {\it all} galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89\% membership accuracy with a mild dependency on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass--cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities ($μ_{\star}$) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation. These results indicate the potential of Copacabana and $μ_{\star}$ to be used in cosmological analyses of optically selected clusters in the future. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 14 pages, 10 figures, submitted to MNRAS

Report number: FERMILAB-PUB-24-0022-PPD

U-Trustworthy Models.Reliability, Competence, and Confidence in Decision-Making

Authors: Ritwik Vashistha, Arya Farahi

Abstract: With growing concerns regarding bias and discrimination in predictive models, the AI community has increasingly focused on assessing AI system trustworthiness. Conventionally, trustworthy AI literature relies on the probabilistic framework and calibration as prerequisites for trustworthiness. In this work, we depart from this viewpoint by proposing a novel trust framework inspired by the philosoph… ▽ More With growing concerns regarding bias and discrimination in predictive models, the AI community has increasingly focused on assessing AI system trustworthiness. Conventionally, trustworthy AI literature relies on the probabilistic framework and calibration as prerequisites for trustworthiness. In this work, we depart from this viewpoint by proposing a novel trust framework inspired by the philosophy literature on trust. We present a precise mathematical definition of trustworthiness, termed $\mathcal{U}$-trustworthiness, specifically tailored for a subset of tasks aimed at maximizing a utility function. We argue that a model's $\mathcal{U}$-trustworthiness is contingent upon its ability to maximize Bayes utility within this task subset. Our first set of results challenges the probabilistic framework by demonstrating its potential to favor less trustworthy models and introduce the risk of misleading trustworthiness assessments. Within the context of $\mathcal{U}$-trustworthiness, we prove that properly-ranked models are inherently $\mathcal{U}$-trustworthy. Furthermore, we advocate for the adoption of the AUC metric as the preferred measure of trustworthiness. By offering both theoretical guarantees and experimental validation, AUC enables robust evaluation of trustworthiness, thereby enhancing model selection and hyperparameter tuning to yield more trustworthy outcomes. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Impact of Property Covariance on Cluster Weak lensing Scaling Relations

Authors: Zhuowen Zhang, Arya Farahi, Daisuke Nagai, Erwin T. Lau, Joshua Frieman, Marina Ricci, Anja von der Linden, Hao-yi Wu

Abstract: We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, $ΔΣ$, and the "true" cluster galaxy number count, $N_{\rm gal}$, as measured within a spherical volume that is void of projection effects. By quantifying the impact of thi… ▽ More We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, $ΔΣ$, and the "true" cluster galaxy number count, $N_{\rm gal}$, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the 3D vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales ($R \lesssim R_{\rm 200c}$) and a null covariance at large scales ($R \gtrsim R_{\rm 200c}$) across most mass and redshift bins. We also find that this covariance results in a $2-3\%$ bias in the halo mass estimates in most bins. Furthermore, by modeling $N_{\rm gal}$ and $ΔΣ$ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the $N_{\rm gal}$-$ΔΣ$ covariance can be explained by the secondary properties of halos that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters. △ Less

Submitted 6 June, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Dark Energy Survey Year 3 Results: Mis-centering calibration and X-ray-richness scaling relations in redMaPPer clusters

Authors: P. Kelly, J. Jobel, O. Eiger, A. Abd, T. E. Jeltema, P. Giles, D. L. Hollowood, R. D. Wilkinson, D. J. Turner, S. Bhargava, S. Everett, A. Farahi, A. K. Romer, E. S. Rykoff, F. Wang, S. Bocquet, D. Cross, R. Faridjoo, J. Franco, G. Gardner, M. Kwiecien, D. Laubner, A. McDaniel, J. H. O'Donnell, L. Sanchez , et al. (54 additional authors not shown)

Abstract: We use Dark Energy Survey Year 3 (DES Y3) clusters with archival X-ray data from XMM-Newton and Chandra to assess the centering performance of the redMaPPer cluster finder and to measure key richness observable scaling relations. In terms of centering, we find that 10-20% of redMaPPer clusters are miscentered with no significant difference in bins of low versus high richness ($20<λ<40$ and $λ>40$)… ▽ More We use Dark Energy Survey Year 3 (DES Y3) clusters with archival X-ray data from XMM-Newton and Chandra to assess the centering performance of the redMaPPer cluster finder and to measure key richness observable scaling relations. In terms of centering, we find that 10-20% of redMaPPer clusters are miscentered with no significant difference in bins of low versus high richness ($20<λ<40$ and $λ>40$) or redshift ($0.2<z<0.4$ and $0.4 <z < 0.65$). We also investigate the richness bias induced by miscentering. The dominant reasons for miscentering include masked or missing data and the presence of other bright galaxies in the cluster; for half of the miscentered clusters the correct central was one of the other possible centrals identified by redMaPPer, while for $\sim 40$% of miscentered clusters the correct central is not a redMaPPer member with most of these cases due to masking. In addition, we fit the scaling relations between X-ray temperature and richness and between X-ray luminosity and richness. We find a T$_X$-$λ$ scatter of $0.21 \pm 0.01$. While the scatter in T$_X$-$λ$ is consistent in bins of redshift, we do find modestly different slopes with high-redshift clusters displaying a somewhat shallower relation. Splitting based on richness, we find a marginally larger scatter for our lowest richness bin, $20 < λ< 40$. The X-ray properties of detected, serendipitous clusters are generally consistent with those for targeted clusters, but the depth of the X-ray data for undetected clusters is insufficient to judge whether they are X-ray underluminous in all but one case. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Fire and Smoke Digital Twin -- A computational framework for modeling fire incident outcomes

Authors: Junfeng Jiao, Ryan Hardesty Lewis, Kijin Seong, Arya Farahi, Paul Navratil, Nate Casebeer, Dev Niyogi

Abstract: Fires and burning are the chief causes of particulate matter (PM2.5), a key measurement of air quality in communities and cities worldwide. This work develops a live fire tracking platform to show active reported fires from over twenty cities in the U.S., as well as predict their smoke paths and impacts on the air quality of regions within their range. Specifically, our close to real-time tracking… ▽ More Fires and burning are the chief causes of particulate matter (PM2.5), a key measurement of air quality in communities and cities worldwide. This work develops a live fire tracking platform to show active reported fires from over twenty cities in the U.S., as well as predict their smoke paths and impacts on the air quality of regions within their range. Specifically, our close to real-time tracking and predictions culminates in a digital twin to protect public health and inform the public of fire and air quality risk. This tool tracks fire incidents in real-time, utilizes the 3D building footprints of Austin to simulate smoke outputs, and predicts fire incident smoke falloffs within the complex city environment. Results from this study include a complete fire and smoke digital twin model for Austin. We work in cooperation with the City of Austin Fire Department to ensure the accuracy of our forecast and also show that air quality sensor density within our cities cannot validate urban fire presence. We additionally release code and methodology to replicate these results for any city in the world. This work paves the path for similar digital twin models to be developed and deployed to better protect the health and safety of citizens. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 8 pages, 8 figures, conference

MSC Class: 68U20; 91D30; 92D30; 93A30 ACM Class: I.6; I.6.7; J.2; H.5.2; J.3

The XMM Cluster Survey: Exploring scaling relations and completeness of the Dark Energy Survey Year 3 redMaPPer cluster catalogue

Authors: E. W. Upsdell, P. A. Giles, A. K. Romer, R. Wilkinson, D. J. Turner, M. Hilton, E. Rykoff, A. Farahi, S. Bhargava, T. Jeltema, M. Klein, A. Bermeo, C. A. Collins, L. Ebrahimpour, D. Hollowood, R. G. Mann, M. Manolopoulou, C. J. Miller, P. J. Rooney, Martin Sahlén, J. P. Stott, P. T. P. Viana, S. Allam, O. Alves, D. Bacon , et al. (45 additional authors not shown)

Abstract: We cross-match and compare characteristics of galaxy clusters identified in observations from two sky surveys using two completely different techniques. One sample is optically selected from the analysis of three years of Dark Energy Survey observations using the redMaPPer cluster detection algorithm. The second is X-ray selected from XMM observations analysed by the XMM Cluster Survey. The sample… ▽ More We cross-match and compare characteristics of galaxy clusters identified in observations from two sky surveys using two completely different techniques. One sample is optically selected from the analysis of three years of Dark Energy Survey observations using the redMaPPer cluster detection algorithm. The second is X-ray selected from XMM observations analysed by the XMM Cluster Survey. The samples comprise a total area of 57.4 deg$^2$, bounded by the area of 4 contiguous XMM survey regions that overlap the DES footprint. We find that the X-ray selected sample is fully matched with entries in the redMaPPer catalogue, above $λ>$20 and within 0.1$< z <$0.9. Conversely, only 38\% of the redMaPPer catalogue is matched to an X-ray extended source. Next, using 120 optically clusters and 184 X-ray selected clusters, we investigate the form of the X-ray luminosity-temperature ($L_{X}-T_{X}$), luminosity-richness ($L_{X}-λ$) and temperature-richness ($T_{X}-λ$) scaling relations. We find that the fitted forms of the $L_{X}-T_{X}$ relations are consistent between the two selection methods and also with other studies in the literature. However, we find tentative evidence for a steepening of the slope of the relation for low richness systems in the X-ray selected sample. When considering the scaling of richness with X-ray properties, we again find consistency in the relations (i.e., $L_{X}-λ$ and $T_{X}-λ$) between the optical and X-ray selected samples. This is contrary to previous similar works that find a significant increase in the scatter of the luminosity scaling relation for X-ray selected samples compared to optically selected samples. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: Accepted for publication to MNRAS

Benchmarks and Explanations for Deep Learning Estimates of X-ray Galaxy Cluster Masses

Authors: Matthew Ho, John Soltis, Arya Farahi, Daisuke Nagai, August Evrard, Michelle Ntampaka

Abstract: We evaluate the effectiveness of deep learning (DL) models for reconstructing the masses of galaxy clusters using X-ray photometry data from next-generation surveys. We establish these constraints using a catalogue of realistic mock eROSITA X-ray observations which use hydrodynamical simulations to model realistic cluster morphology, background emission, telescope response, and AGN sources. Using… ▽ More We evaluate the effectiveness of deep learning (DL) models for reconstructing the masses of galaxy clusters using X-ray photometry data from next-generation surveys. We establish these constraints using a catalogue of realistic mock eROSITA X-ray observations which use hydrodynamical simulations to model realistic cluster morphology, background emission, telescope response, and AGN sources. Using bolometric X-ray photon maps as input, DL models achieve a predictive mass scatter of $σ_{\ln M_\mathrm{500c}} = 17.8\%$, a factor of two improvements on scalar observables such as richness $N_\mathrm{gal}$, 1D velocity dispersion $σ_\mathrm{v,1D}$, and photon count $N_\mathrm{phot}$ as well as a $32\%$ improvement upon idealised, volume-integrated measurements of the bolometric X-ray luminosity $L_X$. We then show that extending this model to handle multichannel X-ray photon maps, separated in low, medium, and high energy bands, further reduces the mass scatter to $16.2\%$. We also tested a multimodal DL model incorporating both dynamical and X-ray cluster probes and achieved marginal gains at a mass scatter of $15.9\%$. Finally, we conduct a quantitative interpretability study of our DL models and find that they greatly down-weight the importance of pixels in the centres of clusters and at the location of AGN sources, validating previous claims of DL modelling improvements and suggesting practical and theoretical benefits for using DL in X-ray mass inference. △ Less

Submitted 25 July, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: 14 pages, 9 figures, 3 tables, accepted in MNRAS

Journal ref: 2023 MNRAS, 524, 3, 3289-3302

Dark Energy Survey Year 3 Results: Redshift Calibration of the MagLim Lens Sample from the combination of SOMPZ and clustering and its impact on Cosmology

Authors: G. Giannini, A. Alarcon, M. Gatti, A. Porredon, M. Crocce, G. M. Bernstein, R. Cawthon, C. Sánchez, C. Doux, J. Elvin-Poole, M. Raveri, J. Myles, A. Amon, S. Allam, O. Alves, F. Andrade-Oliveira, E. Baxter, K. Bechtol, M. R. Becker, J. Blazek, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, A. Choi , et al. (89 additional authors not shown)

Abstract: We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift ca… ▽ More We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift calibration of the lens sample. We describe in detail the methodology, we validate it on simulations and discuss the main effects dominating our error budget. The new calibration is in fair agreement with the fiducial DES Y3 redshift distributions calibration, with only mild differences ($<3σ$) in the means and widths of the distributions. We study the impact of this new calibration on cosmological constraints, analysing DES Y3 galaxy clustering and galaxy-galaxy lensing measurements, assuming a $Λ$CDM cosmology. We obtain $Ω_{\rm m} = 0.30\pm 0.04$, $σ_8 = 0.81\pm 0.07 $ and $S_8 = 0.81\pm 0.04$, which implies a $\sim 0.4σ$ shift in the $Ω_{\rm}-S_8$ plane compared to the fiducial DES Y3 results, highlighting the importance of the redshift calibration of the lens sample in multi-probe cosmological analyses. △ Less

Submitted 18 October, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Journal ref: MNRAS, stad2945 (2023)

A Machine Learning Approach to Enhancing eROSITA Observations

Authors: John Soltis, Michelle Ntampaka, John Wu, John ZuHone, August Evrard, Arya Farahi, Matthew Ho, Daisuke Nagai

Abstract: The eROSITA X-ray telescope, launched in 2019, is predicted to observe roughly 100,000 galaxy clusters. Follow-up observations of these clusters from Chandra, for example, will be needed to resolve outstanding questions about galaxy cluster physics. Deep Chandra cluster observations are expensive and follow-up of every eROSITA cluster is infeasible, therefore, objects chosen for follow-up must be… ▽ More The eROSITA X-ray telescope, launched in 2019, is predicted to observe roughly 100,000 galaxy clusters. Follow-up observations of these clusters from Chandra, for example, will be needed to resolve outstanding questions about galaxy cluster physics. Deep Chandra cluster observations are expensive and follow-up of every eROSITA cluster is infeasible, therefore, objects chosen for follow-up must be chosen with care. To address this, we have developed an algorithm for predicting longer duration, background-free observations based on mock eROSITA observations. We make use of the hydrodynamic cosmological simulation Magneticum, have simulated eROSITA instrument conditions using SIXTE, and have applied a novel convolutional neural network to output a deep Chandra-like "super observation" of each cluster in our simulation sample. Any follow-up merit assessment tool should be designed with a specific use case in mind; our model produces observations that accurately and precisely reproduce the cluster morphology, which is a critical ingredient for determining cluster dynamical state and core type. Our model will advance our understanding of galaxy clusters by improving follow-up selection and demonstrates that image-to-image deep learning algorithms are a viable method for simulating realistic follow-up observations. △ Less

Submitted 9 November, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 21 pages, 11 figures, 3 tables. Minor changes upon revision. Corrected caption of Figure 3. Added discussion of alternative asymmetry metrics. To be published in the Astrophysical Journal

Dark Energy Survey Year 3 Results: Constraints on extensions to $Λ$CDM with weak lensing and galaxy clustering

Authors: DES Collaboration, T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, J. Annis, S. Avila, D. Bacon, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, S. Birrer, J. Blazek, S. Bocquet, A. Brandao-Souza, S. L. Bridle, D. Brooks, D. L. Burke, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, J. Carretero , et al. (137 additional authors not shown)

Abstract: We constrain extensions to the $Λ$CDM model using measurements from the Dark Energy Survey's first three years of observations and external data. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data and blind analyses of real data to validate the robustness of our results. In many cases, constraini… ▽ More We constrain extensions to the $Λ$CDM model using measurements from the Dark Energy Survey's first three years of observations and external data. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data and blind analyses of real data to validate the robustness of our results. In many cases, constraining power is limited by the absence of nonlinear predictions that are reliable at our required precision. The models are: dark energy with a time-dependent equation of state, non-zero spatial curvature, sterile neutrinos, modifications of gravitational physics, and a binned $σ_8(z)$ model which serves as a probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find $(w_{\rm p}, w_a)= (-0.99^{+0.28}_{-0.17},-0.9\pm 1.2)$ at 68% confidence with $z_{\rm p}=0.24$ from the DES measurements alone, and $(w_{\rm p}, w_a)= (-1.03^{+0.04}_{-0.03},-0.4^{+0.4}_{-0.3})$ with $z_{\rm p}=0.21$ for the combination of all data considered. Curvature constraints of $Ω_k=0.0009\pm 0.0017$ and effective relativistic species $N_{\rm eff}=3.10^{+0.15}_{-0.16}$ are dominated by external data. For massive sterile neutrinos, we improve the upper bound on the mass $m_{\rm eff}$ by a factor of three compared to previous analyses, giving 95% limits of $(ΔN_{\rm eff},m_{\rm eff})\leq (0.28, 0.20\, {\rm eV})$. We also constrain changes to the lensing and Poisson equations controlled by functions $Σ(k,z) = Σ_0 Ω_Λ(z)/Ω_{Λ,0}$ and $μ(k,z)=μ_0 Ω_Λ(z)/Ω_{Λ,0}$ respectively to $Σ_0=0.6^{+0.4}_{-0.5}$ from DES alone and $(Σ_0,μ_0)=(0.04\pm 0.05,0.08^{+0.21}_{-0.19})$ for the combination of all data. Overall, we find no significant evidence for physics beyond $Λ$CDM. △ Less

Submitted 29 October, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: Updated to match published version and fix a citation reference. 46 pages, 25 figures, data available at https://dev.des.ncsa.illinois.edu/releases/y3a2/Y3key-extensions

Report number: FERMILAB-PUB-22-470-PPD

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

Correlations of Dark Matter, Gas and Stellar Profiles in Dark Matter Halos

Authors: Arya Farahi, Daisuke Nagai, Dhayaa Anbajagane

Abstract: Halos of similar mass and redshift exhibit a large degree of variability in their differential properties, such as dark matter, hot gas, and stellar mass density profiles. This variability is an indicator of diversity in the formation history of these dark matter halos that is reflected in the coupling of scatters about the mean relations. In this work, we show that the strength of this coupling d… ▽ More Halos of similar mass and redshift exhibit a large degree of variability in their differential properties, such as dark matter, hot gas, and stellar mass density profiles. This variability is an indicator of diversity in the formation history of these dark matter halos that is reflected in the coupling of scatters about the mean relations. In this work, we show that the strength of this coupling depends on the scale at which halo profiles are measured. By analyzing the outputs of the IllustrisTNG hydrodynamical cosmological simulations we report the radial- and mass-dependent couplings between the dark matter, hot gas, and stellar mass radial density profiles utilizing the population diversity in dark matter halos. We find that for the same mass halos the scatters in the density of baryons and dark matter are strongly coupled at large scales ($r>R_{200}$); but the coupling between gas and dark matter density profiles fades near the core of halos ($r < 0.3 R_{200}$). We then show that the correlation between halo profile and integrated quantities induces a radius-dependent additive bias in the profile observables of halos when halos are selected on properties other than their mass. We discuss the impact of this effect on cluster abundance and cross-correlations cosmology with multi-wavelength cosmological surveys. △ Less

Submitted 31 May, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: 8 figures, 4 appendices, data in Table 2 (Accepted)

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

The XMM Cluster Survey analysis of the SDSS DR8 redMaPPer Catalogue: Implications for scatter, selection bias, and isotropy in cluster scaling relations

Authors: P. A. Giles, A. K. Romer, R. Wilkinson, A. Bermeo, D. J. Turner, M. Hilton, E. W. Upsdell, P. J. Rooney, S. Bhargava, L. Ebrahimpour, A. Farahi, R. G. Mann, M. Manolopoulou, J. Mayers, C. Vergara, P. T. P. Viana, C. A. Collins, D. Hollowood, T. Jeltema, C. J. Miller, R. C. Nichol, R. Noorali, M. Splettstoesser, J. P. Stott

Abstract: In this paper we present the X-ray analysis of SDSS DR8 redMaPPer (SDSSRM) clusters using data products from the $XMM$ Cluster Survey (XCS). In total, 1189 SDSSRM clusters fall within the $XMM$-Newton footprint. This has yielded 456 confirmed detections accompanied by X-ray luminosity ($L_{X}$) measurements. Of the detected clusters, 382 have an associated X-ray temperature measurement ($T_{X}$).… ▽ More In this paper we present the X-ray analysis of SDSS DR8 redMaPPer (SDSSRM) clusters using data products from the $XMM$ Cluster Survey (XCS). In total, 1189 SDSSRM clusters fall within the $XMM$-Newton footprint. This has yielded 456 confirmed detections accompanied by X-ray luminosity ($L_{X}$) measurements. Of the detected clusters, 382 have an associated X-ray temperature measurement ($T_{X}$). This represents one of the largest samples of coherently derived cluster $T_{X}$ values to date. Our analysis of the X-ray observable to richness ($λ$) scaling relations has demonstrated that scatter in the $T_{X}-λ$ relation is roughly a third of that in the $L_{X}-λ$ relation, and that the $L_{X}-λ$ scatter is intrinsic, i.e. will not be significantly reduced with larger sample sizes. Our analysis of the scaling relation between $L_{X}$ and $T_{X}$ has shown that the fits are sensitive to the selection method of the sample, i.e. whether the sample is made up of clusters detected "serendipitously" compared to those deliberately targeted by $XMM$. These differences are also seen in the $L_{X}-λ$ relation and, to a lesser extent, in the $T_{X}-λ$ relation. Exclusion of the emission from the cluster core does not make a significant impact to the findings. A combination of selection biases is a likely, but as yet unproven, reason for these differences. Finally, we have also used our data to probe recent claims of anisotropy in the $L_{X}-T_{X}$ relation across the sky. We find no evidence of anistropy, but stress that this may be masked in our analysis by the incomplete declination coverage of the SDSS DR8 sample. △ Less

Submitted 22 August, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: 23 pages, 20 figures, Accepted for publication to MNRAS

KLLR: A scale-dependent, multivariate model class for regression analysis

Authors: Arya Farahi, Dhayaa Anbajagane, August Evrard

Abstract: The underlying physics of astronomical systems governs the relation between their measurable properties. Consequently, quantifying the statistical relationships between system-level observable properties of a population offers insights into the astrophysical drivers of that class of systems. While purely linear models capture behavior over a limited range of system scale, the fact that astrophysic… ▽ More The underlying physics of astronomical systems governs the relation between their measurable properties. Consequently, quantifying the statistical relationships between system-level observable properties of a population offers insights into the astrophysical drivers of that class of systems. While purely linear models capture behavior over a limited range of system scale, the fact that astrophysics is ultimately scale-dependent implies the need for a more flexible approach to describing population statistics over a wide dynamic range. For such applications, we introduce and implement a class of Kernel-Localized Linear Regression (KLLR) models. KLLR is a natural extension to the commonly-used linear models that allows the parameters of the linear model -- normalization, slope, and covariance matrix -- to be scale-dependent. KLLR performs inference in two steps: (1) it estimates the mean relation between a set of independent variables and a dependent variable and; (2) it estimates the conditional covariance of the dependent variables given a set of independent variables. We demonstrate the model's performance in a simulated setting and showcase an application of the proposed model in analyzing the baryonic content of dark matter halos. As a part of this work, we publicly release a Python implementation of the KLLR method. △ Less

Submitted 26 April, 2022; v1 submitted 20 February, 2022; originally announced February 2022.

Comments: The code is publicly available at https://github.com/afarahi/kllr and can be installed through `pip install kllr`

HSC-XXL : Baryon budget of the 136 XXL Groups and Clusters

Authors: Daichi Akino, Dominique Eckert, Nobuhiro Okabe, Mauro Sereno, Keiichi Umetsu, Masamune Oguri, Fabio Gastaldello, I-Non Chiu, Stefano Ettori, August E. Evrard, Arya Farahi, Ben Maughan, Marguerite Pierre, Marina Ricci, Ivan Valtchanov, Ian Mccarthy, Sean Mcgee, Satoshi Miyazaki, Atsushi J. Nishizawa, Masayuki Tanaka

Abstract: We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry.… ▽ More We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass--total mass ($M_{500}$) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass--total mass relation is shallower than unity, $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of $M_{500}$ agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution. △ Less

Submitted 16 January, 2022; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 34 pages, 16 figures, 8 tables; accepted for the publication in PASJ

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Authors: Dhayaa Anbajagane, Han Aung, August E. Evrard, Arya Farahi, Daisuke Nagai, David J. Barnes, Weiguang Cui, Klaus Dolag, Ian G. McCarthy, Elena Rasia, Gustavo Yepes

Abstract: Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the rat… ▽ More Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available. △ Less

Submitted 21 December, 2021; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: [v1]: 8 figures, 14 pages in main text + 5 figures, 3 pages in Appendix. [v2]: Minor text edits. Accepted in MNRAS

Baryonic Imprints on DM Halos: Population Statistics from Dwarf Galaxies to Galaxy Clusters

Authors: Dhayaa Anbajagane, August E. Evrard, Arya Farahi

Abstract: In a purely cold dark matter universe, the initial matter power spectrum and its subsequent gravitational growth contain no special mass- or time-scales, and so neither do the emergent population statistics of internal dark matter (DM) halo properties. Using 1.5 million halos from three IllustrisTNG realizations of a LambdaCDM universe, we show that galaxy formation physics drives non-monotonic fe… ▽ More In a purely cold dark matter universe, the initial matter power spectrum and its subsequent gravitational growth contain no special mass- or time-scales, and so neither do the emergent population statistics of internal dark matter (DM) halo properties. Using 1.5 million halos from three IllustrisTNG realizations of a LambdaCDM universe, we show that galaxy formation physics drives non-monotonic features ("wiggles") into DM property statistics across six decades in halo mass, from dwarf galaxies to galaxy clusters. We characterize these features by extracting the halo mass-dependent statistics of five DM halo properties -- velocity dispersion, NFW concentration, density- and velocity-space shapes, and formation time -- using kernel-localized linear regression (KLLR). Comparing precise estimates of normalizations, slopes, and covariances between realizations with and without galaxy formation, we find systematic deviations across all mass-scales, with maximum deviations of 25% at the Milky-Way mass of 1e12 Msun. The mass-dependence of the wiggles is set by the interplay between different cooling and feedback mechanisms, and we discuss its observational implications. The property covariances depend strongly on halo mass and physics treatment, but the correlations are mostly robust. Using multivariate KLLR and interpretable machine learning, we show the halo concentration and velocity-space shape are principal contributors, at different mass, to the velocity dispersion variance. Statistics of mass accretion rate and DM surface pressure energy are provided in an appendix. We publicly release halo property catalogs and KLLR parameters for the TNG runs at twenty epochs up to z = 12. △ Less

Submitted 20 November, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: [v1]: 10 figures, 16 pages in Main text, and 5 figures, 5 pages in Appendix. [v2]: Minor text edit. Accepted in MNRAS

Velocity Dispersions of Clusters in the Dark Energy Survey Y3 redMaPPer Catalog

Authors: V. Wetzell, T. E. Jeltema, B. Hegland, S. Everett, P. A. Giles, R. Wilkinson, A. Farahi, M. Costanzi, D. L. Hollowood, E. Upsdell, A. Saro, J. Myles, A. Bermeo, S. Bhargava, C. A. Collins, D. Cross, O. Eiger, G. Gardner, M. Hilton, J. Jobel, P. Kelly, D. Laubner, A. R. Liddle, R. G. Mann, V. Martinez , et al. (74 additional authors not shown)

Abstract: We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlation… ▽ More We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlations between cluster velocity dispersion, richness, X-ray temperature and luminosity as well as central galaxy velocity offsets. The velocity dispersion-richness relation exhibits a bimodal distribution. The majority of clusters follow scaling relations between velocity dispersion, richness, and X-ray properties similar to those found for previous samples; however, there is a significant population of clusters with velocity dispersions which are high for their richness. These clusters account for roughly 22\% of the $λ< 70$ systems in our sample, but more than half (55\%) of $λ< 70$ clusters at $z>0.5$. A couple of these systems are hot and X-ray bright as expected for massive clusters with richnesses that appear to have been underestimated, but most appear to have high velocity dispersions for their X-ray properties likely due to line-of-sight structure. These results suggest that projection effects contribute significantly to redMaPPer selection, particularly at higher redshifts and lower richnesses. The redMaPPer determined richnesses for the velocity dispersion outliers are consistent with their X-ray properties, but several are X-ray undetected and deeper data is needed to understand their nature. △ Less

Submitted 9 June, 2022; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 22 pages, accepted to MNRAS

Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

Authors: DES Collaboration, T. M. C. Abbott, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon, F. Andrade-Oliveira, J. Annis, S. Avila, D. Bacon, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, S. Bhargava, S. Birrer, J. Blazek, A. Brandao-Souza, S. L. Bridle, D. Brooks, E. Buckley-Geer, D. L. Burke, H. Camacho, A. Campos , et al. (146 additional authors not shown)

Abstract: We present the first cosmology results from large-scale structure in the Dark Energy Survey (DES) spanning 5000 deg$^2$. We perform an analysis combining three two-point correlation functions (3$\times$2pt): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions. The analysis was designed to miti… ▽ More We present the first cosmology results from large-scale structure in the Dark Energy Survey (DES) spanning 5000 deg$^2$. We perform an analysis combining three two-point correlation functions (3$\times$2pt): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions. The analysis was designed to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results. We model the data within the flat $Λ$CDM and $w$CDM cosmological models. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude $S_8=0.776^{+0.017}_{-0.017}$ and matter density $Ω_{\mathrm{m}} = 0.339^{+0.032}_{-0.031}$ in $Λ$CDM, mean with 68% confidence limits; $S_8=0.775^{+0.026}_{-0.024}$, $Ω_{\mathrm{m}} = 0.352^{+0.035}_{-0.041}$, and dark energy equation-of-state parameter $w=-0.98^{+0.32}_{-0.20}$ in $w$CDM. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed $p=0.13$ to $0.48$. When combining DES 3$\times$2pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find $p=0.34$. Combining all of these data sets with Planck CMB lensing yields joint parameter constraints of $S_8 = 0.812^{+0.008}_{-0.008}$, $Ω_{\mathrm{m}} = 0.306^{+0.004}_{-0.005}$, $h=0.680^{+0.004}_{-0.003}$, and $\sum m_ν<0.13 \;\mathrm{eV\; (95\% \;CL)}$ in $Λ$CDM; $S_8 = 0.812^{+0.008}_{-0.008}$, $Ω_{\mathrm{m}} = 0.302^{+0.006}_{-0.006}$, $h=0.687^{+0.006}_{-0.007}$, and $w=-1.031^{+0.030}_{-0.027}$ in $w$CDM. (abridged) △ Less

Submitted 21 March, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: See https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ for the full DES Y3 3x2pt cosmology release. Matches version accepted in PRD

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: High-precision measurement and modeling of galaxy-galaxy lensing

Authors: J. Prat, J. Blazek, C. Sánchez, I. Tutusaus, S. Pandey, J. Elvin-Poole, E. Krause, M. A. Troxel, L. F. Secco, A. Amon, J. DeRose, G. Zacharegkas, C. Chang, B. Jain, N. MacCrann, Y. Park, E. Sheldon, G. Giannini, S. Bocquet, C. To, A. Alarcon, O. Alves, F. Andrade-Oliveira, E. Baxter, K. Bechtol , et al. (116 additional authors not shown)

Abstract: We present and characterize the galaxy-galaxy lensing signal measured using the first three years of data from the Dark Energy Survey (DES Y3) covering 4132 deg$^2$. These galaxy-galaxy measurements are used in the DES Y3 3$\times$2pt cosmological analysis, which combines weak lensing and galaxy clustering information. We use two lens samples: a magnitude-limited sample and the redMaGic sample, wh… ▽ More We present and characterize the galaxy-galaxy lensing signal measured using the first three years of data from the Dark Energy Survey (DES Y3) covering 4132 deg$^2$. These galaxy-galaxy measurements are used in the DES Y3 3$\times$2pt cosmological analysis, which combines weak lensing and galaxy clustering information. We use two lens samples: a magnitude-limited sample and the redMaGic sample, which span the redshift range $\sim 0.2-1$ with 10.7 M and 2.6 M galaxies respectively. For the source catalog, we use the Metacalibration shape sample, consisting of $\simeq$100 M galaxies separated into 4 tomographic bins. Our galaxy-galaxy lensing estimator is the mean tangential shear, for which we obtain a total S/N of $\sim$148 for MagLim ($\sim$120 for redMaGic), and $\sim$67 ($\sim$55) after applying the scale cuts of 6 Mpc/$h$. Thus we reach percent-level statistical precision, which requires that our modeling and systematic-error control be of comparable accuracy. The tangential shear model used in the 3$\times$2pt cosmological analysis includes lens magnification, a five-parameter intrinsic alignment model (TATT), marginalization over a point-mass to remove information from small scales and a linear galaxy bias model validated with higher-order terms. We explore the impact of these choices on the tangential shear observable and study the significance of effects not included in our model, such as reduced shear, source magnification and source clustering. We also test the robustness of our measurements to various observational and systematics effects, such as the impact of observing conditions, lens-source clustering, random-point subtraction, scale-dependent Metacalibration responses, PSF residuals, and B-modes. △ Less

Submitted 5 April, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: 33 pages, 18 figures, accepted by PRD

Report number: FERMILAB-PUB-21-248-AE

SHAPing the Gas: Understanding Gas Shapes in Dark Matter Haloes with Interpretable Machine Learning

Authors: Luis Fernando Machado Poletti Valle, Camille Avestruz, David J. Barnes, Arya Farahi, Erwin T. Lau, Daisuke Nagai

Abstract: The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machin… ▽ More The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machine learning approach for modelling the dependence of gas shapes on dark matter and baryonic properties. With data from the IllustrisTNG hydrodynamical cosmological simulations, we develop a machine learning pipeline that applies \pkg{XGBoost}, an implementation of gradient boosted decision trees, to predict radial profiles of gas shapes from halo properties. We show that \pkg{XGBoost} models can accurately predict gas shape profiles in dark matter haloes. We also explore model interpretability with \pkg{SHAP}, a method that identifies the most predictive properties at different halo radii. We find that baryonic properties best predict gas shapes in halo cores, whereas dark matter shapes are the main predictors in the halo outskirts. This work demonstrates the power of interpretable machine learning in modelling observable properties of dark matter haloes in the era of multi-wavelength cosmological surveys. △ Less

Submitted 2 August, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 14 pages, 8 figures, Accepted 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)

PoPE: A population-based approach to model spatial structure of astronomical systems

Authors: Arya Farahi, Daisuke Nagai, Yang Chen

Abstract: We present a novel population-based Bayesian inference approach to model the average and population variance of spatial distribution of a set of observables from ensemble analysis of low signal-to-noise ratio measurements. The method consists of (1) inferring the average profile using Gaussian Processes and (2) computing the covariance of the profile observables given a set of independent variable… ▽ More We present a novel population-based Bayesian inference approach to model the average and population variance of spatial distribution of a set of observables from ensemble analysis of low signal-to-noise ratio measurements. The method consists of (1) inferring the average profile using Gaussian Processes and (2) computing the covariance of the profile observables given a set of independent variables. Our model is computationally efficient and capable of inferring average profiles of a large population size from noisy measurements, without stacking and binning data nor parameterizing the shape of the mean profile. We demonstrate the performance of our method using dark matter, gas and stellar profiles extracted from hydrodynamical cosmological simulations of galaxy formation. Population Profile Estimator (PoPE) is publicly available in a GitHub repository. Our new method should be useful for measuring the spatial distribution and internal structure of a variety of astrophysical systems using large astronomical surveys. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 14 pages, 9 figures, 2 tables, the code is publicly available at https://github.com/afarahi/PoPE

Approximate Bayesian Uncertainties on Deep Learning Dynamical Mass Estimates of Galaxy Clusters

Authors: Matthew Ho, Arya Farahi, Markus Michael Rau, Hy Trac

Abstract: We study methods for reconstructing Bayesian uncertainties on dynamical mass estimates of galaxy clusters using convolutional neural networks (CNNs). We discuss the statistical background of approximate Bayesian neural networks and demonstrate how variational inference techniques can be used to perform computationally tractable posterior estimation for a variety of deep neural architectures. We ex… ▽ More We study methods for reconstructing Bayesian uncertainties on dynamical mass estimates of galaxy clusters using convolutional neural networks (CNNs). We discuss the statistical background of approximate Bayesian neural networks and demonstrate how variational inference techniques can be used to perform computationally tractable posterior estimation for a variety of deep neural architectures. We explore how various model designs and statistical assumptions impact prediction accuracy and uncertainty reconstruction in the context of cluster mass estimation. We measure the quality of our model posterior recovery using a mock cluster observation catalog derived from the MultiDark simulation and UniverseMachine catalog. We show that approximate Bayesian CNNs produce highly accurate dynamical cluster mass posteriors. These model posteriors are log-normal in cluster mass and recover $68\%$ and $90\%$ confidence intervals to within $1\%$ of their measured value. We note how this rigorous modeling of dynamical mass posteriors is necessary for using cluster abundance measurements to constrain cosmological parameters. △ Less

Submitted 13 March, 2021; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 15 pages, 5 figures. 3 tables, accepted for publication at ApJ

Journal ref: 2021 ApJ, 908, 204H

$μ_{\star}$ Masses: Weak Lensing Calibration of the Dark Energy Survey Year 1 redMaPPer Clusters using Stellar Masses

Authors: M. E. S. Pereira, A. Palmese, T. N. Varga, T. McClintock, M. Soares-Santos, J. Burgad, J. Annis, A. Farahi, H. Lin, A. Choi, J. DeRose, J. Esteves, M. Gatti, D. Gruen, W. G. Hartley, B. Hoyle, T. Jeltema, N. MacCrann, A. Roodman, C. Sánchez, T. Shin, A. von der Linden, J. Zuntz, T. M. C. Abbott, M. Aguena , et al. (56 additional authors not shown)

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

Submitted 17 June, 2020; originally announced June 2020.

Comments: 20 pages, 10 figures, submitted to MNRAS

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)

Driving with Data in the Motor City: Mining and Modeling Vehicle Fleet Maintenance Data

Authors: Josh Gardner, Jawad Mroueh, Natalia Jenuwine, Noah Weaverdyck, Samuel Krassenstein, Arya Farahi, Danai Koutra

Abstract: The City of Detroit maintains an active fleet of over 2500 vehicles, spending an annual average of over \$5 million on purchases and over \$7.7 million on maintenance. Modeling patterns and trends in this data is of particular importance to a variety of stakeholders, particularly as Detroit emerges from Chapter 9 bankruptcy, but the structure in such data is complex, and the city lacks dedicated r… ▽ More The City of Detroit maintains an active fleet of over 2500 vehicles, spending an annual average of over \$5 million on purchases and over \$7.7 million on maintenance. Modeling patterns and trends in this data is of particular importance to a variety of stakeholders, particularly as Detroit emerges from Chapter 9 bankruptcy, but the structure in such data is complex, and the city lacks dedicated resources for in-depth analysis. The City of Detroit's Operations and Infrastructure Group and the University of Michigan initiated a collaboration which seeks to address this unmet need by analyzing data from the City of Detroit's vehicle fleet. This work presents a case study and provides the first data-driven benchmark, demonstrating a suite of methods to aid in data understanding and prediction for large vehicle maintenance datasets. We present analyses to address three key questions raised by the stakeholders, related to discovering multivariate maintenance patterns over time; predicting maintenance; and predicting vehicle- and fleet-level costs. We present a novel algorithm, PRISM, for automating multivariate sequential data analyses using tensor decomposition. This work is a first of its kind that presents both methodologies and insights to guide future civic data research. △ Less

Submitted 21 September, 2020; v1 submitted 23 February, 2020; originally announced February 2020.

Aging Halos: Implications of the Magnitude Gap on Conditional Statistics of Stellar and Gas Properties of Massive Halos

Authors: Arya Farahi, Matthew Ho, Hy Trac

Abstract: Cold dark matter model predicts that the large-scale structure grows hierarchically. Small dark matter halos form first. Then, they grow gradually via continuous merger and accretion. These halos host the majority of baryonic matter in the Universe in the form of hot gas and cold stellar phase. Determining how baryons are partitioned into these phases requires detailed modeling of galaxy formation… ▽ More Cold dark matter model predicts that the large-scale structure grows hierarchically. Small dark matter halos form first. Then, they grow gradually via continuous merger and accretion. These halos host the majority of baryonic matter in the Universe in the form of hot gas and cold stellar phase. Determining how baryons are partitioned into these phases requires detailed modeling of galaxy formation and their assembly history. It is speculated that formation time of the same mass halos might be correlated with their baryonic content. To evaluate this hypothesis, we employ halos of mass above $10^{14}\,M_{\odot}$ realized by TNG300 solution of the IllustrisTNG project. Formation time is not directly observable. Hence, we rely on the magnitude gap between the brightest and the fourth brightest halo galaxy member, which is shown that traces formation time of the host halo. We compute the conditional statistics of the stellar and gas content of halos conditioned on their total mass and magnitude gap. We find a strong correlation between magnitude gap and gas mass, BCG stellar mass, and satellite galaxies stellar mass, but not the total stellar mass of halo. Conditioning on the magnitude gap can reduce the scatter about halo property--halo mass relation and has a significant impact on the conditional covariance. Reduction in the scatter can be as significant as 30%, which implies more accurate halo mass prediction. Incorporating the magnitude gap has a potential to improve cosmological constraints using halo abundance and allows us to gain insight into the baryon evolution within these systems. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: 14 pages, 6 Figures, 6 Tables, MNRAS

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

Authors: Dhayaa Anbajagane, August E. Evrard, Arya Farahi, David J. Barnes, Klaus Dolag, Ian G. McCarthy, Dylan Nelson, Annalisa Pillepich

Abstract: We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generate… ▽ More We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $γ= -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $γ\simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 17 pages, 15 figures

Response to NITRD, NCO, NSF Request for Information on "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan"

Authors: J. Amundson, J. Annis, C. Avestruz, D. Bowring, J. Caldeira, G. Cerati, C. Chang, S. Dodelson, D. Elvira, A. Farahi, K. Genser, L. Gray, O. Gutsche, P. Harris, J. Kinney, J. B. Kowalkowski, R. Kutschke, S. Mrenna, B. Nord, A. Para, K. Pedro, G. N. Perdue, A. Scheinker, P. Spentzouris, J. St. John , et al. (5 additional authors not shown)

Abstract: We present a response to the 2018 Request for Information (RFI) from the NITRD, NCO, NSF regarding the "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan." Through this document, we provide a response to the question of whether and how the National Artificial Intelligence Research and Development Strategic Plan (NAIRDSP) should be updated from the perspect… ▽ More We present a response to the 2018 Request for Information (RFI) from the NITRD, NCO, NSF regarding the "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan." Through this document, we provide a response to the question of whether and how the National Artificial Intelligence Research and Development Strategic Plan (NAIRDSP) should be updated from the perspective of Fermilab, America's premier national laboratory for High Energy Physics (HEP). We believe the NAIRDSP should be extended in light of the rapid pace of development and innovation in the field of Artificial Intelligence (AI) since 2016, and present our recommendations below. AI has profoundly impacted many areas of human life, promising to dramatically reshape society --- e.g., economy, education, science --- in the coming years. We are still early in this process. It is critical to invest now in this technology to ensure it is safe and deployed ethically. Science and society both have a strong need for accuracy, efficiency, transparency, and accountability in algorithms, making investments in scientific AI particularly valuable. Thus far the US has been a leader in AI technologies, and we believe as a national Laboratory it is crucial to help maintain and extend this leadership. Moreover, investments in AI will be important for maintaining US leadership in the physical sciences. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Report number: FERMILAB-FN-1092-SCD

Detection of anti-correlation of hot and cold baryons in galaxy clusters

Authors: Arya Farahi, Sarah L. Mulroy, August E. Evrard, Graham P. Smith, Alexis Finoguenov, Herve Bourdin, John E. Carlstrom, Chris P. Haines, Daniel P. Marrone, Rossella Martino, Pasquale Mazzotta, Christine O'Donnell, Nobuhiro Okabe

Abstract: The largest clusters of galaxies in the Universe contain vast amounts of dark matter, plus baryonic matter in two principal phases, a majority hot gas component and a minority cold stellar phase comprising stars, compact objects, and low-temperature gas. Hydrodynamic simulations indicate that the highest-mass systems retain the cosmic fraction of baryons, a natural consequence of which is anti-cor… ▽ More The largest clusters of galaxies in the Universe contain vast amounts of dark matter, plus baryonic matter in two principal phases, a majority hot gas component and a minority cold stellar phase comprising stars, compact objects, and low-temperature gas. Hydrodynamic simulations indicate that the highest-mass systems retain the cosmic fraction of baryons, a natural consequence of which is anti-correlation between the masses of hot gas and stars within dark matter halos of fixed total mass. We report observational detection of this anti-correlation based on 4 elements of a $9\times9$ element covariance matrix for nine cluster properties, measured from X-ray, optical, infrared and millimetre wavelength observations of 41 clusters from the Local Cluster Substructure Survey. These clusters were selected using explicit and quantitative selection rules that were then encoded in our hierarchical Bayesian model. Our detection of anti-correlation is consistent with predictions from contemporary hydrodynamic cosmological simulations that were not tuned to reproduce this signal. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: Published in the 2 July issue of Nature Communications. For published version, see https://rdcu.be/bICJG . The full posterior chains is available at https://doi.org/10.6084/m9.figshare.8001218

Identification of RR Lyrae stars in multiband, sparsely-sampled data from the Dark Energy Survey using template fitting and Random Forest classification

Authors: K. M. Stringer, J. P. Long, L. M. Macri, J. L. Marshall, A. Drlica-Wagner, C. E. Martínez-Vázquez, A. K. Vivas, K. Bechtol, E. Morganson, M. Carrasco Kind, A. B. Pace, A. R. Walker, C. Nielsen, T. S. Li, E. Rykoff, D. Burke, A. Carnero Rosell, E. Neilsen, P. Ferguson, S. A. Cantu, J. L. Myron, L. Strigari, A. Farahi, F. Paz-Chinchón, D. Tucker , et al. (53 additional authors not shown)

Abstract: Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multi… ▽ More Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets. △ Less

Submitted 1 May, 2019; originally announced May 2019.

Comments: 34 pages, 21 figures. Accepted for publication in AJ. Data products are available at https://des.ncsa.illinois.edu/releases/other/y3-rrl

Stellar mass as a galaxy cluster mass proxy: application to the Dark Energy Survey redMaPPer clusters

Authors: A. Palmese, J. Annis, J. Burgad, A. Farahi, M. Soares-Santos, B. Welch, M. da Silva Pereira, H. Lin, S. Bhargava, D. L. Hollowood, R. Wilkinson, P. Giles, T. Jeltema, A. K. Romer, A. E. Evrard, M. Hilton, C. Vergara Cervantes, A. Bermeo, J. Mayers, J. DeRose, D. Gruen, W. G. Hartley, O. Lahav, B. Leistedt, T. McClintock , et al. (60 additional authors not shown)

Abstract: We introduce a galaxy cluster mass observable, $μ_\star$, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 observations. Stellar masses are computed using a Bayesian Model Averaging method, and are validated for DES data using simulations and COSMOS data. We show that $μ_\star$ works as a promising mass proxy by comparing our prediction… ▽ More We introduce a galaxy cluster mass observable, $μ_\star$, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 observations. Stellar masses are computed using a Bayesian Model Averaging method, and are validated for DES data using simulations and COSMOS data. We show that $μ_\star$ works as a promising mass proxy by comparing our predictions to X-ray measurements. We measure the X-ray temperature-$μ_\star$ relation for a total of 150 clusters matched between the wide-field DES Year 1 redMaPPer catalogue, and Chandra and XMM archival observations, spanning the redshift range $0.1<z<0.7$. For a scaling relation which is linear in logarithmic space, we find a slope of $α= 0.488\pm0.043$ and a scatter in the X-ray temperature at fixed $μ_\star$ of $σ_{{\rm ln} T_X|μ_\star}=0.266^{+0.019}_{-0.020}$ for the joint sample. By using the halo mass scaling relations of the X-ray temperature from the Weighing the Giants program, we further derive the $μ_\star$-conditioned scatter in mass, finding $σ_{{\rm ln} M|μ_\star}=0.26^{+ 0.15}_{- 0.10}$. These results are competitive with well-established cluster mass proxies used for cosmological analyses, showing that $μ_\star$ can be used as a reliable and physically motivated mass proxy to derive cosmological constraints. △ Less

Submitted 18 November, 2019; v1 submitted 20 March, 2019; originally announced March 2019.

Comments: 14 pages, 7 figures, addressing MNRAS referee comments

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

Mass Variance from Archival X-ray Properties of Dark Energy Survey Year-1 Galaxy Clusters

Authors: A. Farahi, X. Chen, A. E. Evrard, D. L. Hollowood, R. Wilkinson, S. Bhargava, P. Giles, A. K. Romer, T. Jeltema, M. Hilton, A. Bermeo, J. Mayers, C. Vergara Cervantes, E. Rozo, E. S. Rykoff, C. Collins, M. Costanzi, S. Everett, A. R. Liddle, R. G. Mann, A. Mantz, P. Rooney, M. Sahlen, J. Stott, P. T. P. Viana , et al. (54 additional authors not shown)

Abstract: Using archival X-ray observations and a log-normal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behavior of X-ray temperatures, $T_X$, with optical richness, $λ_{RM}$, for clusters in the redshift range… ▽ More Using archival X-ray observations and a log-normal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behavior of X-ray temperatures, $T_X$, with optical richness, $λ_{RM}$, for clusters in the redshift range $0.2<z<0.7$. X-ray temperatures are obtained from Chandra and XMM observations for 58 and 110 redMaPPer systems, respectively. Despite non-uniform sky coverage, the $T_X$ measurements are $> 50\%$ complete for clusters with $λ_{RM} > 130$. Regression analysis on the two samples produces consistent posterior scaling parameters, from which we derive a combined constraint on the residual scatter, $σ_{\ln Tx | λ} = 0.275 \pm 0.019$. Joined with constraints for $T_X$ scaling with halo mass from the Weighing the Giants program and richness--temperature covariance estimates from the LoCuSS sample, we derive the richness-conditioned scatter in mass, $σ_{\ln M | λ} = 0.30 \pm 0.04\, _{({\rm stat})} \pm 0.09\, _{({\rm sys})}$, at an optical richness of approximately 70. Uncertainties in external parameters, particularly the slope and variance of the $T_X$--mass relation and the covariance of $T_X$ and $λ_{RM}$ at fixed mass, dominate the systematic error. The $95\%$ confidence region from joint sample analysis is relatively broad, $σ_{\ln M | λ} \in [0.14, \, 0.55]$, or a factor ten in variance. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Comments: 14 pages. Main results are Figure 3, 5, and 6, Table 2, and Equation 9. Submitted to MNRAS. Comments welcome

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

The Role of Machine Learning in the Next Decade of Cosmology

Authors: Michelle Ntampaka, Camille Avestruz, Steven Boada, Joao Caldeira, Jessi Cisewski-Kehe, Rosanne Di Stefano, Cora Dvorkin, August E. Evrard, Arya Farahi, Doug Finkbeiner, Shy Genel, Alyssa Goodman, Andy Goulding, Shirley Ho, Arthur Kosowsky, Paul La Plante, Francois Lanusse, Michelle Lochner, Rachel Mandelbaum, Daisuke Nagai, Jeffrey A. Newman, Brian Nord, J. E. G. Peek, Austin Peel, Barnabas Poczos , et al. (5 additional authors not shown)

Abstract: In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster an… ▽ More In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors. △ Less

Submitted 14 January, 2021; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: Submitted to the Astro2020 call for science white papers

Percent-Level Test of Isotropic Expansion Using Type Ia Supernovae

Authors: John Soltis, Arya Farahi, Dragan Huterer, C. Michael Liberato II

Abstract: We propose and implement a novel, robust, and non-parametric test of statistical isotropy of the expansion of the universe, and apply it to around one thousand type Ia supernovae from the Pantheon sample. We calculate the angular clustering of supernova magnitude residuals and compare it to the noise expected under the isotropic assumption. We also test for systematic effects and demonstrate that… ▽ More We propose and implement a novel, robust, and non-parametric test of statistical isotropy of the expansion of the universe, and apply it to around one thousand type Ia supernovae from the Pantheon sample. We calculate the angular clustering of supernova magnitude residuals and compare it to the noise expected under the isotropic assumption. We also test for systematic effects and demonstrate that their effects are negligible or are already accounted for in our procedure. We express our constraints as an upper limit on the rms spatial variation in the Hubble parameter at late times. For the sky smoothed with a Gaussian with fwhm=60 deg, less than 1% rms spatial variation in the Hubble parameter is allowed at 99.7% confidence. △ Less

Submitted 23 May, 2019; v1 submitted 19 February, 2019; originally announced February 2019.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 122, 091301 (2019)

A Robust and Efficient Deep Learning Method for Dynamical Mass Measurements of Galaxy Clusters

Authors: Matthew Ho, Markus Michael Rau, Michelle Ntampaka, Arya Farahi, Hy Trac, Barnabas Poczos

Abstract: We demonstrate the ability of convolutional neural networks (CNNs) to mitigate systematics in the virial scaling relation and produce dynamical mass estimates of galaxy clusters with remarkably low bias and scatter. We present two models, CNN$_\mathrm{1D}$ and CNN$_\mathrm{2D}$, which leverage this deep learning tool to infer cluster masses from distributions of member galaxy dynamics. Our first m… ▽ More We demonstrate the ability of convolutional neural networks (CNNs) to mitigate systematics in the virial scaling relation and produce dynamical mass estimates of galaxy clusters with remarkably low bias and scatter. We present two models, CNN$_\mathrm{1D}$ and CNN$_\mathrm{2D}$, which leverage this deep learning tool to infer cluster masses from distributions of member galaxy dynamics. Our first model, CNN$_\text{1D}$, infers cluster mass directly from the distribution of member galaxy line-of-sight velocities. Our second model, CNN$_\text{2D}$, extends the input space of CNN$_\text{1D}$ to learn on the joint distribution of galaxy line-of-sight velocities and projected radial distances. We train each model as a regression over cluster mass using a labeled catalog of realistic mock cluster observations generated from the MultiDark simulation and UniverseMachine catalog. We then evaluate the performance of each model on an independent set of mock observations selected from the same simulated catalog. The CNN models produce cluster mass predictions with lognormal residuals of scatter as low as $0.132$ dex, greater than a factor of 2 improvement over the classical $M$-$σ$ power-law estimator. Furthermore, the CNN model reduces prediction scatter relative to similar machine learning approaches by up to $17\%$ while executing in drastically shorter training and evaluation times (by a factor of 30) and producing considerably more robust mass predictions (improving prediction stability under variations in galaxy sampling rate by $30\%$). △ Less

Submitted 22 December, 2020; v1 submitted 15 February, 2019; originally announced February 2019.

Comments: 22 pages, 10 figures, 4 tables, accepted for publication at ApJ

Journal ref: 2019 ApJ, 887, 25

LoCuSS: Scaling relations between galaxy cluster mass, gas, and stellar content

Authors: Sarah L. Mulroy, Arya Farahi, August E. Evrard, Graham P. Smith, Alexis Finoguenov, Christine O'Donnell, Daniel P. Marrone, Zubair Abdulla, Hervé Bourdin, John E. Carlstrom, Jessica Démoclès, Chris P. Haines, Rossella Martino, Pasquale Mazzotta, Sean L. McGee, Nobuhiro Okabe

Abstract: We present a simultaneous analysis of galaxy cluster scaling relations between weak-lensing mass and multiple cluster observables, across a wide range of wavelengths, that probe both gas and stellar content. Our new hierarchical Bayesian model simultaneously considers the selection variable alongside all other observables in order to explicitly model intrinsic property covariance and account for s… ▽ More We present a simultaneous analysis of galaxy cluster scaling relations between weak-lensing mass and multiple cluster observables, across a wide range of wavelengths, that probe both gas and stellar content. Our new hierarchical Bayesian model simultaneously considers the selection variable alongside all other observables in order to explicitly model intrinsic property covariance and account for selection effects. We apply this method to a sample of 41 clusters at $0.15<z<0.30$, with a well-defined selection criteria based on RASS X-ray luminosity, and observations from Chandra / XMM, SZA, Planck, UKIRT, SDSS and Subaru. These clusters have well-constrained weak-lensing mass measurements based on Subaru / Suprime-Cam observations, which serve as the reference masses in our model. We present 30 scaling relation parameters for 10 properties. All relations probing the intracluster gas are slightly shallower than self-similar predictions, in moderate tension with prior measurements, and the stellar fraction decreases with mass. K-band luminosity has the lowest intrinsic scatter with a 95th percentile of 0.16, while the lowest scatter gas probe is gas mass with a fractional intrinsic scatter of $0.16 \pm 0.03$. We find no distinction between the core-excised X-ray or high-resolution Sunyaev-Zel'dovich relations of clusters of different central entropy, but find with modest significance that higher entropy clusters have higher stellar fractions than their lower entropy counterparts. We also report posterior mass estimates from our likelihood model. △ Less

Submitted 31 January, 2019; originally announced January 2019.

Comments: 21 pages, 9 figures

Journal ref: MNRAS, Volume 484, Issue 1, 2019, Pages 60-80

Dark Energy Survey Year 1 Results: Calibration of Cluster Mis-centering in the redMaPPer Catalogs

Authors: Y. Zhang, T. Jeltema, D. L. Hollowood, S. Everett, E. Rozo, A. Farahi, A. Bermeo, S. Bhargava, P. Giles, A. K. Romer, R. Wilkinson, E. S. Rykoff, A. Mantz, H. T. Diehl, A. E. Evrard, C. Stern, D. Gruen, A. von der Linden, M. Splettstoesser, X. Chen, M. Costanzi, S. Allen, C. Collins, M. Hilton, M. Klein , et al. (61 additional authors not shown)

Abstract: The center determination of a galaxy cluster from an optical cluster finding algorithm can be offset from theoretical prescriptions or $N$-body definitions of its host halo center. These offsets impact the recovered cluster statistics, affecting both richness measurements and the weak lensing shear profile around the clusters. This paper models the centering performance of the \RM~cluster finding… ▽ More The center determination of a galaxy cluster from an optical cluster finding algorithm can be offset from theoretical prescriptions or $N$-body definitions of its host halo center. These offsets impact the recovered cluster statistics, affecting both richness measurements and the weak lensing shear profile around the clusters. This paper models the centering performance of the \RM~cluster finding algorithm using archival X-ray observations of \RM-selected clusters. Assuming the X-ray emission peaks as the fiducial halo centers, and through analyzing their offsets to the \RM~centers, we find that $\sim 75\pm 8 \%$ of the \RM~clusters are well centered and the mis-centered offset follows a Gamma distribution in normalized, projected distance. These mis-centering offsets cause a systematic underestimation of cluster richness relative to the well-centered clusters, for which we propose a descriptive model. Our results enable the DES Y1 cluster cosmology analysis by characterizing the necessary corrections to both the weak lensing and richness abundance functions of the DES Y1 redMaPPer cluster catalog. △ Less

Submitted 1 June, 2019; v1 submitted 21 January, 2019; originally announced January 2019.

Comments: Code used in this analysis is available from https://github.com/yyzhang/center_modeling_y1

Dark Energy Survey Year 1 Results: Methods for Cluster Cosmology and Application to the SDSS

Authors: M. Costanzi, E. Rozo, M. Simet, Y. Zhang, A. E. Evrard, A. Mantz, E. S. Rykoff, T. Jeltema, D. Gruen, S. Allen T. McClintock, A. K. Romer, A. von der Linden, A. Farahi, J. DeRose, T. N. Varga, J. Weller, P. Giles, D. L. Hollowood, S. Bhargava, A. Bermeo-Hernandez, X. Chen, T. M. C. Abbott, F. B. Abdalla, S. Avila, K. Bechtol , et al. (61 additional authors not shown)

Abstract: We perform the first blind analysis of cluster abundance data. Specifically, we derive cosmological constraints from the abundance and weak-lensing signal of \redmapper\ clusters of richness $λ\geq 20$ in the redshift range $z\in[0.1,0.3]$ as measured in the Sloan Digital Sky Survey (SDSS). We simultaneously fit for cosmological parameters and the richness--mass relation of the clusters. For a fla… ▽ More We perform the first blind analysis of cluster abundance data. Specifically, we derive cosmological constraints from the abundance and weak-lensing signal of \redmapper\ clusters of richness $λ\geq 20$ in the redshift range $z\in[0.1,0.3]$ as measured in the Sloan Digital Sky Survey (SDSS). We simultaneously fit for cosmological parameters and the richness--mass relation of the clusters. For a flat $Λ$CDM cosmological model with massive neutrinos, we find $S_8 \equiv σ_{8}(Ω_m/0.3)^{0.5}=0.79^{+0.05}_{-0.04}$. This value is both consistent and competitive with that derived from cluster catalogues selected in different wavelengths. Our result is also consistent with the combined probes analyses by the Dark Energy Survey (DES) and the Kilo-Degree Survey (KiDS), and with the Cosmic Microwave Background (CMB) anisotropies as measured by \planck. We demonstrate that the cosmological posteriors are robust against variation of the richness--mass relation model and to systematics associated with the calibration of the selection function. In combination with Baryon Acoustic Oscillation (BAO) data and Big-Bang Nucleosynthesis (BBN) data, we constrain the Hubble rate to be $h=0.66\pm 0.02$, independent of the CMB. Future work aimed at improving our understanding of the scatter of the richness--mass relation has the potential to significantly improve the precision of our cosmological posteriors. The methods described in this work were developed for use in the forthcoming analysis of cluster abundances in the DES. Our SDSS analysis constitutes the first part of a staged-unblinding analysis of the full DES data set. △ Less

Submitted 22 October, 2018; originally announced October 2018.

Comments: 25 pages, 14 figures

Gaia DR2 unravels incompleteness of nearby cluster population: New open clusters in the direction of Perseus

Authors: T. Cantat-Gaudin, A. Krone-Martins, N. Sedaghat, A. Farahi, R. S. de Souza, R. Skalidis, A. I. Malz, S. Macêdo, B. Moews, C. Jordi, A. Moitinho, A. Castro-Ginard, E. E. O. Ishida, C. Heneka, A. Boucaud, A. M. M. Trindade

Abstract: Open clusters (OCs) are popular tracers of the structure and evolutionary history of the Galactic disk. The OC population is often considered to be complete within 1.8 kpc of the Sun. The recent Gaia Data Release 2 (DR2) allows the latter claim to be challenged. We perform a systematic search for new OCs in the direction of Perseus using precise and accurate astrometry from Gaia DR2. We implement… ▽ More Open clusters (OCs) are popular tracers of the structure and evolutionary history of the Galactic disk. The OC population is often considered to be complete within 1.8 kpc of the Sun. The recent Gaia Data Release 2 (DR2) allows the latter claim to be challenged. We perform a systematic search for new OCs in the direction of Perseus using precise and accurate astrometry from Gaia DR2. We implement a coarse-to-fine search method. First, we exploit spatial proximity using a fast density-aware partitioning of the sky via a k-d tree in the spatial domain of Galactic coordinates, (l, b). Secondly, we employ a Gaussian mixture model in the proper motion space to quickly tag fields around OC candidates. Thirdly, we apply an unsupervised membership assignment method, UPMASK, to scrutinise the candidates. We visually inspect colour-magnitude diagrams to validate the detected objects. Finally, we perform a diagnostic to quantify the significance of each identified overdensity in proper motion and in parallax space We report the discovery of 41 new stellar clusters. This represents an increment of at least 20% of the previously known OC population in this volume of the Milky Way. We also report on the clear identification of NGC 886, an object previously considered an asterism. This letter challenges the previous claim of a near-complete sample of open clusters up to 1.8 kpc. Our results reveal that this claim requires revision, and a complete census of nearby open clusters is yet to be found. △ Less

Submitted 21 March, 2019; v1 submitted 12 October, 2018; originally announced October 2018.

Comments: accepted for publication in A&A

Journal ref: A&A 624, A126 (2019)

Chandra Follow-Up of the SDSS DR8 redMaPPer Catalog Using the MATCha Pipeline

Authors: Devon L. Hollowood, Tesla Jeltema, Xinyi Chen, Arya Farahi, August Evrard, Spencer Everett, Eduardo Rozo, Eli Rykoff, Rebecca Bernstein, Alberto Bermeo, Lena Eiger, Paul Giles, Holger Israel, Renee Michel, Raziq Noorali, Kathy Romer, Philip Rooney, Megan Splettstoesser

Abstract: In order to place constraints on cosmology through optical surveys of galaxy clusters, one must first understand the properties of those clusters. To this end, we introduce the Mass Analysis Tool for Chandra (MATCha), a pipeline which uses a parallellized algorithm to analyze archival Chandra data. MATCha simultaneously calculates X-ray temperatures and luminosities and performs centering measurem… ▽ More In order to place constraints on cosmology through optical surveys of galaxy clusters, one must first understand the properties of those clusters. To this end, we introduce the Mass Analysis Tool for Chandra (MATCha), a pipeline which uses a parallellized algorithm to analyze archival Chandra data. MATCha simultaneously calculates X-ray temperatures and luminosities and performs centering measurements for hundreds of potential galaxy clusters using archival X-ray exposures. We run MATCha on the redMaPPer SDSS DR8 cluster catalog and use MATCha's output X-ray temperatures and luminosities to analyze the galaxy cluster temperature-richness, luminosity-richness, luminosity-temperature, and temperature-luminosity scaling relations. We detect 447 clusters and determine 246 r2500 temperatures across all redshifts. Within 0.1 < z < 0.35 we find that r2500 Tx scales with optical richness as ln(kB Tx / 1.0 keV) = (0.52 \pm 0.05) ln(λ/70) + (1.85 \pm 0.03) with intrinsic scatter of 0.27 \pm 0.02 (1 σ). We investigate the distribution of offsets between the X-ray center and redMaPPer center within 0.1 < z < 0.35, finding that 68.3 \pm 6.5% of clusters are well-centered. However, we find a broad tail of large offsets in this distribution, and we explore some of the causes of redMaPPer miscentering. △ Less

Submitted 13 August, 2019; v1 submitted 20 August, 2018; originally announced August 2018.

Comments: 36 pages, 17 figures, 7 tables. Accepted to ApJS

Modeling projection effects in optically-selected cluster catalogues

Authors: M. Costanzi, E. Rozo, E. S. Rykoff, A. Farahi, T. Jeltema, A. E. Evrard, A. Mantz, D. Gruen, R. Mandelbaum, J. DeRose, T. McClintock, T. N. Varga, Y. Zhang, J. Weller, R. H. Wechsler, M. Aguena

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

Submitted 18 July, 2018; originally announced July 2018.

Comments: 16 pages, 11 figures

ActiveRemediation: The Search for Lead Pipes in Flint, Michigan

Authors: Jacob Abernethy, Alex Chojnacki, Arya Farahi, Eric Schwartz, Jared Webb

Abstract: We detail our ongoing work in Flint, Michigan to detect pipes made of lead and other hazardous metals. After elevated levels of lead were detected in residents' drinking water, followed by an increase in blood lead levels in area children, the state and federal governments directed over $125 million to replace water service lines, the pipes connecting each home to the water system. In the absence… ▽ More We detail our ongoing work in Flint, Michigan to detect pipes made of lead and other hazardous metals. After elevated levels of lead were detected in residents' drinking water, followed by an increase in blood lead levels in area children, the state and federal governments directed over $125 million to replace water service lines, the pipes connecting each home to the water system. In the absence of accurate records, and with the high cost of determining buried pipe materials, we put forth a number of predictive and procedural tools to aid in the search and removal of lead infrastructure. Alongside these statistical and machine learning approaches, we describe our interactions with government officials in recommending homes for both inspection and replacement, with a focus on the statistical model that adapts to incoming information. Finally, in light of discussions about increased spending on infrastructure development by the federal government, we explore how our approach generalizes beyond Flint to other municipalities nationwide. △ Less

Submitted 17 August, 2018; v1 submitted 10 June, 2018; originally announced June 2018.

Comments: 10 pages, 10 figures, To appear in KDD 2018, For associated promotional video, see https://www.youtube.com/watch?v=YbIn_axYu9E

Astronomers' and Physicists' Attitudes Toward Education & Public Outreach: A Programmatic Study with The Dark Energy Survey

Authors: A. Farahi, R. R. Gupta, C. Krawiec, A. A. Plazas, R. C. Wolf

Abstract: We present a programmatic study of physicists' and astronomers' attitudes toward education and public outreach (EPO) using 131 survey responses from members of the Dark Energy Survey. We find a disparity between the types of EPO activities researchers deem valuable and those in which they participate. Most respondents are motivated to engage in EPO by a desire to educate the public. Barriers to en… ▽ More We present a programmatic study of physicists' and astronomers' attitudes toward education and public outreach (EPO) using 131 survey responses from members of the Dark Energy Survey. We find a disparity between the types of EPO activities researchers deem valuable and those in which they participate. Most respondents are motivated to engage in EPO by a desire to educate the public. Barriers to engagement include career- and skill-related concerns, but lack of time is the main deterrent. We explore the value of centralized EPO efforts and conclude with a list of recommendations for increasing researchers' engagement. △ Less

Submitted 18 June, 2019; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 16 pages, 5 figures; metadata updated 18 June 2019

Journal ref: Published in Journal of STEM Outreach, 2019 6-11. pp 1-15. https://tinyurl.com/y6d7mg36