pop-cosmos: Scaleable inference of galaxy properties and redshifts with a data-driven population model

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Thorp, Stephen, Alsing, Justin, Peiris, Hiranya V., Deger, Sinan, Mortlock, Daniel J., Leistedt, Boris, Leja, Joel, Loureiro, Arthur
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866917822951587840
author Thorp, Stephen
Alsing, Justin
Peiris, Hiranya V.
Deger, Sinan
Mortlock, Daniel J.
Leistedt, Boris
Leja, Joel
Loureiro, Arthur
author_facet Thorp, Stephen
Alsing, Justin
Peiris, Hiranya V.
Deger, Sinan
Mortlock, Daniel J.
Leistedt, Boris
Leja, Joel
Loureiro, Arthur
contents We present an efficient Bayesian method for estimating individual photometric redshifts and galaxy properties under a pre-trained population model (pop-cosmos) that was calibrated using purely photometric data. This model specifies a prior distribution over 16 stellar population synthesis (SPS) parameters using a score-based diffusion model, and includes a data model with detailed treatment of nebular emission. We use a GPU-accelerated affine invariant ensemble sampler to achieve fast posterior sampling under this model for 292,300 individual galaxies in the COSMOS2020 catalog, leveraging a neural network emulator (Speculator) to speed up the SPS calculations. We apply both the pop-cosmos population model and a baseline prior inspired by Prospector-$α$, and compare these results to published COSMOS2020 redshift estimates from the widely-used EAZY and LePhare codes. For the $\sim 12,000$ galaxies with spectroscopic redshifts, we find that pop-cosmos yields redshift estimates that have minimal bias ($\sim10^{-4}$), high accuracy ($σ_\text{MAD}=7\times10^{-3}$), and a low outlier rate ($1.6\%$). We show that the pop-cosmos population model generalizes well to galaxies fainter than its $r<25$ mag training set. The sample we have analyzed is $\gtrsim3\times$ larger than has previously been possible via posterior sampling with a full SPS model, with average throughput of 15 GPU-sec per galaxy under the pop-cosmos prior, and 0.6 GPU-sec per galaxy under the Prospector prior. This paves the way for principled modeling of the huge catalogs expected from upcoming Stage IV galaxy surveys.
format Preprint
id arxiv_https___arxiv_org_abs_2406_19437
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle pop-cosmos: Scaleable inference of galaxy properties and redshifts with a data-driven population model
Thorp, Stephen
Alsing, Justin
Peiris, Hiranya V.
Deger, Sinan
Mortlock, Daniel J.
Leistedt, Boris
Leja, Joel
Loureiro, Arthur
Cosmology and Nongalactic Astrophysics
Astrophysics of Galaxies
Instrumentation and Methods for Astrophysics
We present an efficient Bayesian method for estimating individual photometric redshifts and galaxy properties under a pre-trained population model (pop-cosmos) that was calibrated using purely photometric data. This model specifies a prior distribution over 16 stellar population synthesis (SPS) parameters using a score-based diffusion model, and includes a data model with detailed treatment of nebular emission. We use a GPU-accelerated affine invariant ensemble sampler to achieve fast posterior sampling under this model for 292,300 individual galaxies in the COSMOS2020 catalog, leveraging a neural network emulator (Speculator) to speed up the SPS calculations. We apply both the pop-cosmos population model and a baseline prior inspired by Prospector-$α$, and compare these results to published COSMOS2020 redshift estimates from the widely-used EAZY and LePhare codes. For the $\sim 12,000$ galaxies with spectroscopic redshifts, we find that pop-cosmos yields redshift estimates that have minimal bias ($\sim10^{-4}$), high accuracy ($σ_\text{MAD}=7\times10^{-3}$), and a low outlier rate ($1.6\%$). We show that the pop-cosmos population model generalizes well to galaxies fainter than its $r<25$ mag training set. The sample we have analyzed is $\gtrsim3\times$ larger than has previously been possible via posterior sampling with a full SPS model, with average throughput of 15 GPU-sec per galaxy under the pop-cosmos prior, and 0.6 GPU-sec per galaxy under the Prospector prior. This paves the way for principled modeling of the huge catalogs expected from upcoming Stage IV galaxy surveys.
title pop-cosmos: Scaleable inference of galaxy properties and redshifts with a data-driven population model
topic Cosmology and Nongalactic Astrophysics
Astrophysics of Galaxies
Instrumentation and Methods for Astrophysics
url https://arxiv.org/abs/2406.19437