We present
Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.
more » « less- NSF-PAR ID:
- 10401679
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 521
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 1741-1756
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract pop-cosmos : a comprehensive model characterizing the galaxy population, calibrated to 140,938 (r < 25 selected) galaxies from the Cosmic Evolution Survey (COSMOS) with photometry in 26 bands from the ultraviolet to the infrared. We construct a detailed forward model for the COSMOS data, comprising: a population model describing the joint distribution of galaxy characteristics and its evolution (parameterized by a flexible score-based diffusion model); a state-of-the-art stellar population synthesis model connecting galaxies’ intrinsic properties to their photometry; and a data model for the observation, calibration, and selection processes. By minimizing the optimal transport distance between synthetic and real data, we are able to jointly fit the population and data models, leading to robustly calibrated population-level inferences that account for parameter degeneracies, photometric noise and calibration, and selection. We present a number of key predictions from our model of interest for cosmology and galaxy evolution, including the mass function and redshift distribution; the mass–metallicity-redshift and fundamental metallicity relations; the star-forming sequence; the relation between dust attenuation and stellar mass, star formation rate, and attenuation-law index; and the relation between gas-ionization and star formation. Our model encodes a comprehensive picture of galaxy evolution that faithfully predicts galaxy colors across a broad redshift (z < 4) and wavelength range. -
more » « less
Abstract The Probabilistic Value-Added Bright Galaxy Survey (PROVABGS) catalog will provide the posterior distributions of physical properties of >10 million DESI Bright Galaxy Survey galaxies. Each posterior distribution will be inferred from joint Bayesian modeling of observed photometry and spectroscopy using Markov Chain Monte Carlo sampling and the Hahn et al. stellar population synthesis (SPS) model. To make this computationally feasible, PROVABGS will use a neural emulator for the SPS model to accelerate the posterior inference. In this work, we present how we construct the emulator using the Alsing et al. approach and verify that it can be used to accurately infer galaxy properties. We confirm that the emulator is in excellent agreement with the original SPS model with ≪1% error and is 100× faster. In addition, we demonstrate that the posteriors of galaxy properties derived using the emulator are also in excellent agreement with those inferred using the original model. The neural emulator presented in this work is essential in bypassing the computational challenge posed in constructing the PROVABGS catalog. Furthermore, it demonstrates the advantages of emulation for scaling sophisticated analyses to millions of galaxies.
-
more » « less
ABSTRACT We present Diffstar , a smooth parametric model for the in situ star formation history (SFH) of galaxies. The Diffstar model is distinct from traditional SFH models because it is parametrized directly in terms of basic features of galaxy formation physics. Diffstar includes ingredients for: the halo mass assembly history; the accretion of gas into the dark matter halo; the fraction of gas that is eventually transformed into stars, ϵms; the time-scale over which this transformation occurs, τcons; and the possibility that some galaxies will experience a quenching event at time tq, and may subsequently experience rejuvenated star formation. We show that our model is sufficiently flexible to describe the average stellar mass histories of galaxies in both the IllustrisTNG (TNG) and UniverseMachine (UM) simulations with an accuracy of ∼0.1 dex across most of cosmic time. We use Diffstar to compare TNG to UM in common physical terms, finding that: (i) star formation in UM is less efficient and burstier relative to TNG; (ii) UM galaxies have longer gas consumption time-scales, relative to TNG; (iii) rejuvenated star formation is ubiquitous in UM, whereas quenched TNG galaxies rarely experience sustained rejuvenation; and (iv) in both simulations, the distributions of ϵms, τcons, and tq share a common characteristic dependence upon halo mass, and present significant correlations with halo assembly history. We conclude with a discussion of how Diffstar can be used in future applications to fit the SEDs of individual observed galaxies, as well as in forward-modelling applications that populate cosmological simulations with synthetic galaxies.
-
ABSTRACT We present predictions for the evolution of the galaxy dust-to-gas ratio (DGR) and dust-to-metal ratio (DTM) from z = 0 → 6, using a model for the production, growth, and destruction of dust grains implemented into the simba cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z = 0, but modestly underpredicts the mass function by ∼×3 at z ∼ 1–2. The z = 0 DGR versus metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR–metallicity relationship between z = 0 and 6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion time-scale, stellar mass, and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM versus the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.more » « less
-
more » « less
Abstract Stellar kinematics and metallicity are key to exploring formation scenarios for galactic disks and halos. In this work, we characterized the relationship between kinematics and photometric metallicity along the line of sight to M31's disk. We combined optical Hubble Space Telescope/Advanced Camera for Surveys photometry, from the Panchromatic Hubble Andromeda Treasury survey, with Keck/DEIMOS spectra, from the Spectroscopic and Photometric Landscape of Andromeda’s Stellar Halo survey. The resulting sample of 3512 individual red giant branch stars spans 4–19 projected kpc, making it a useful probe of both the disk and inner halo. We separated these stars into disk and halo populations, by modeling the line-of-sight velocity distributions as a function of position across the disk region, where ∼73% stars have a high likelihood of belonging to the disk and ∼14% to the halo. Although stellar halos are typically thought to be metal-poor, the kinematically identified halo contains a significant population of stars (∼29%) with disk-like metallicity ([Fe/H]phot∼ −0.10). This metal-rich halo population lags the gaseous disk to a similar extent as the rest of the halo, indicating that it does not correspond to a canonical thick disk. Its properties are inconsistent with those of tidal debris originating from the Giant Stellar Stream merger event. Moreover, the halo is chemically distinct from the phase-mixed component previously identified along the minor axis (i.e., away from the disk), implying contributions from different formation channels. These metal-rich halo stars provide direct chemodynamical evidence in favor of the previously suggested “kicked-up” disk population in M31's inner stellar halo.
