skip to main content


Title: DSPS: Differentiable stellar population synthesis
ABSTRACT

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
Author(s) / Creator(s):
; ; ; ;
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
  1. 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.

     
    more » « less
  2. Abstract

    Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. BelowM∼ 109M, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way (MW). Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume (D≲ 12 Mpc) to date. We use 250 confirmed and 71 candidate dwarf satellites around 27 MW-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semianalyticalSatGenmodel for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMFs), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low-mass end are populated by a relation of the formMMpeakα, with a moderate slope ofαconst=2.10±0.01and a low scatter, constant as a function of the peak halo mass, ofσconst=0.060.05+0.07. A model with a steeper slope (αgrow= 2.39 ± 0.06) and a scatter that grows with decreasingMpeakis also consistent with the observed SSMF but is not required. Our new model for the satellite–subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only MW satellites.

     
    more » « less
  3. Abstract

    We report the first star formation history study of the Milky Ways nuclear star cluster (NSC), which includes observational constraints from a large sample of stellar metallicity measurements. These metallicity measurements were obtained from recent surveys from Gemini and the Very Large Telescope of 770 late-type stars within the central 1.5 pc. These metallicity measurements, along with photometry and spectroscopically derived temperatures, are forward modeled with a Bayesian inference approach. Including metallicity measurements improves the overall fit quality, as the low-temperature red giants that were previously difficult to constrain are now accounted for, and the best fit favors a two-component model. The dominant component contains 93% ± 3% of the mass, is metal-rich ([M/H]¯0.45), and has an age of52+3Gyr, which is ∼3 Gyr younger than earlier studies with fixed (solar) metallicity; this younger age challenges coevolutionary models in which the NSC and supermassive black holes formed simultaneously at early times. The minor population component has low metallicity ([M/H]¯1.1) and contains ∼7% of the stellar mass. The age of the minor component is uncertain (0.1–5 Gyr old). Using the estimated parameters, we infer the following NSC stellar remnant population (with ∼18% uncertainty): 1.5 × 105neutron stars, 2.5 × 105stellar-mass black holes (BHs), and 2.2 × 104BH–BH binaries. These predictions result in 2–4 times fewer neutron stars compared to earlier predictions that assume solar metallicity, introducing a possible new path to understand the so-called “missing-pulsar problem”. Finally, we present updated predictions for the BH–BH merger rates (0.01–3 Gpc−3yr−1).

     
    more » « less
  4. ABSTRACT

    The relationships between stellar mass, gas-phase metallicity and star-formation rate (i.e. the mass–metallicity, MZR, and the fundamental metallicity relation, FMR) in the local Universe are revisited by fully anchoring the metallicity determination for SDSS galaxies on the Te abundance scale defined exploiting the strong-line metallicity calibrations presented by Curti et al. Self-consistent metallicity measurements allow a more unbiased assessment of the scaling relations involving M, Z and SFR, which provide powerful constraints for the chemical evolution models. We parametrize the MZR with a new functional form that allows us to better characterize the turnover mass. The slope and saturation metallicity are in good agreement with previous determinations of the MZR based on the Te method, while showing significantly lower normalization compared to those based on photoionization models. The Z–SFR dependence at fixed stellar mass is also investigated, being particularly evident for highly star-forming galaxies, where the scatter in metallicity is reduced up to a factor of ${\sim}30{{\ \rm per\ cent}}$. A new parametrization of the FMR is given by explicitly introducing the SFR dependence of the turnover mass into the MZR. The residual scatter in metallicity for the global galaxy population around the new FMR is 0.054 dex. The new FMR presented in this work represents a useful local benchmark to compare theoretical predictions and observational studies (of both local and high-redshift galaxies) whose metallicity measurements are tied to the abundance scale defined by the Te method, hence allowing proper assessment of its evolution with cosmic time.

     
    more » « less
  5. ABSTRACT

    White dwarf binaries with orbital periods below 1 h will be the most numerous sources for the space-based gravitational wave detector Laser Interferometer Space Antenna (LISA). Based on thousands of individually resolved systems, we will be able to constrain binary evolution and provide a new map of the Milky Way and its close surroundings. In this paper we predict the main properties of populations of different types of detached white dwarf binaries detected by LISA over time. For the first time, we combine a high-resolution cosmological simulation of a Milky Way-mass galaxy (taken from the FIRE project) with a binary population synthesis model for low- and intermediate-mass stars. Our Galaxy model therefore provides a cosmologically realistic star formation and metallicity history for the Galaxy and naturally produces its different components such as the thin and thick disc, the bulge, the stellar halo, and satellite galaxies and streams. Thanks to the simulation, we show how different Galactic components contribute differently to the gravitational wave signal, mostly due to their typical age and distance distributions. We find that the dominant LISA sources will be He–He double white dwarfs (DWDs) and He–CO DWDs with important contributions from the thick disc and bulge. The resulting sky map of the sources is different from previous models, with important consequences for the searches for electromagnetic counterparts and data analysis. We also emphasize that much of the science-enabling information regarding white dwarf binaries, such as the chirp mass and the sky localization, becomes increasingly rich with long observations, including an extended mission up to 8 yr.

     
    more » « less