Stellar Population Inference with Prospector Johnson, Benjamin D.; Leja, Joel; Conroy, Charlie ...
The Astrophysical journal. Supplement series,
06/2021, Letnik:
254, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Abstract
Inference of the physical properties of stellar populations from observed photometry and spectroscopy is a key goal in the study of galaxy evolution. In recent years, the quality and ...quantity of the available data have increased, and there have been corresponding efforts to increase the realism of the stellar population models used to interpret these observations. Describing the observed galaxy spectral energy distributions in detail now requires physical models with a large number of highly correlated parameters. These models do not fit easily on grids and necessitate a full exploration of the available parameter space. We present
Prospector
, a flexible code for inferring stellar population parameters from photometry and spectroscopy spanning UV through IR wavelengths. This code is based on forward modeling the data and Monte Carlo sampling the posterior parameter distribution, enabling complex models and exploration of moderate dimensional parameter spaces. We describe the key ingredients of the code and discuss the general philosophy driving the design of these ingredients. We demonstrate some capabilities of the code on several data sets, including mock and real data.
Accounting for nebular emission when modeling galaxy spectral energy distributions (SEDs) is important, as both line and continuum emissions can contribute significantly to the total observed flux. ...In this work, we present a new nebular emission model integrated within the Flexible Stellar Population Synthesis code that computes the line and continuum emission for complex stellar populations using the photoionization code Cloudy. The self-consistent coupling of the nebular emission to the matched ionizing spectrum produces emission line intensities that correctly scale with the stellar population as a function of age and metallicity. This more complete model of galaxy SEDs will improve estimates of global gas properties derived with diagnostic diagrams, star formation rates based on H , and physical properties derived from broadband photometry. Our models agree well with results from other photoionization models and are able to reproduce observed emission from H ii regions and star-forming galaxies. Our models show improved agreement with the observed H ii regions in the Ne iii/O ii plane and show satisfactory agreement with He ii emission from z = 2 galaxies, when including rotating stellar models. Models including post-asymptotic giant branch stars are able to reproduce line ratios consistent with low-ionization emission regions. The models are integrated into current versions of FSPS and include self-consistent nebular emission predictions for MIST and Padova+Geneva evolutionary tracks.
Nonparametric star formation histories (SFHs) have long promised to be the "gold standard" for galaxy spectral energy distribution (SED) modeling as they are flexible enough to describe the full ...diversity of SFH shapes, whereas parametric models rule out a significant fraction of these shapes a priori. However, this flexibility is not fully constrained even with high-quality observations, making it critical to choose a well-motivated prior. Here, we use the SED-fitting code Prospector to explore the effect of different nonparametric priors by fitting SFHs to mock UV-IR photometry generated from a diverse set of input SFHs. First, we confirm that nonparametric SFHs recover input SFHs with less bias and return more accurate errors than do parametric SFHs. We further find that, while nonparametric SFHs robustly recover the overall shape of the input SFH, the primary determinant of the size and shape of the posterior star formation rate as a function of time (SFR(t)) is the choice of prior, rather than the photometric noise. As a practical demonstration, we fit the UV-IR photometry of ∼6000 galaxies from the Galaxy and Mass Assembly survey and measure scatters between priors to be 0.1 dex in mass, 0.8 dex in SFR100 Myr, and 0.2 dex in mass-weighted ages, with the bluest star-forming galaxies showing the most sensitivity. An important distinguishing characteristic for nonparametric models is the characteristic timescale for changes in SFR(t). This difference controls whether galaxies are assembled in bursts or in steady-state star formation, corresponding respectively to (feedback-dominated/accretion-dominated) models of galaxy formation and to (larger/smaller) confidence intervals derived from SED fitting. High-quality spectroscopy has the potential to further distinguish between these proposed models of SFR(t).
ABSTRACT This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar ...evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ( ), masses ( ), and metallicities ( ). The models are self-consistently and continuously evolved from the pre-main sequence (PMS) to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the PMS to the end of core helium burning for . We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.
We explore the connection between the UV luminosity functions (LFs) of high-z galaxies and the distribution of stellar masses and star formation histories (SFHs) in their host dark matter halos. We ...provide a baseline for a redshift-independent star formation efficiency model to which observations and models can be compared. Our model assigns a star formation rate (SFR) to each dark matter halo based on the growth rate of the halo and a redshift-independent star formation efficiency. The dark matter halo accretion rate is obtained from a high-resolution N-body simulation in order to capture the stochasticity in accretion histories and to obtain spatial information for the distribution of galaxies. The halo mass dependence of the star formation efficiency is calibrated at z = 4 by requiring a match to the observed UV LF at this redshift. The model then correctly predicts the observed UV LF at z = 5-10. We present predictions for the UV luminosity and stellar mass functions, JWST number counts, and SFHs. In particular, we find a stellar-to-halo mass relation at z = 4-10 that scales with halo mass at Mh < 1011 M as M ∝ Mh2, with a normalization that is higher than the relation inferred at z = 0. The average SFRs increase as a function of time to z = 4, although there is significant scatter around the average: about 6% of the z = 4 galaxies show no significant mass growth. Using these SFHs, we present redshift-dependent UV-to-SFR conversion factors, mass return fractions, and mass-to-light ratios for different initial mass functions and metallicities, finding that current estimates of the cosmic SFR density at z ∼ 10 may be overestimated by ∼0.1-0.2 dex.
In this paper, we provide a physical model for the origin of variations in the shapes and bump strengths of dust attenuation laws in galaxies by combining a large suite of cosmological "zoom-in" ...galaxy formation simulations with 3D Monte Carlo dust radiative transfer calculations. We model galaxies over three orders of magnitude in stellar mass, ranging from Milky Way-like systems to massive galaxies at high redshift. Critically, for these calculations, we employ a constant underlying dust extinction law in all cases and examine how the role of geometry and radiative transfer effects impacts the resultant attenuation curves. Our main results follow. Despite our usage of a constant dust extinction curve, we find dramatic variations in the derived attenuation laws. The slopes of normalized attenuation laws depend primarily on the complexities of star-to-dust geometry. Increasing fractions of unobscured young stars flatten normalized curves, while increasing fractions of unobscured old stars steepen curves. Similar to the slopes of our model attenuation laws, we find dramatic variation in the 2175 ultraviolet bump strength, including a subset of curves with little to no bump. These bump strengths are primarily influenced by the fraction of unobscured O and B stars in our model, with the impact of scattered light having only a secondary effect. Taken together, these results lead to a natural relationship between the attenuation curve slope and 2175 bump strength. Finally, we apply these results to a 25 Mpc h−1 box cosmological hydrodynamic simulation in order to model the expected dispersion in attenuation laws at integer redshifts from z = 0 to 6. A significant dispersion is expected at low redshifts and decreases toward z = 6. We provide tabulated results for the best-fit median attenuation curve at all redshifts.
Abstract
The primary method for inferring the stellar mass (
M
*
) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy ...star formation history (SFH) and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here we examine the effect that the assumed SFH has on the stellar properties inferred from SED fitting by ground-truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with underreported uncertainties on results. We demonstrate that flexible nonparametric SFHs outperform traditional parametric forms in capturing variations in galaxy SFHs and, as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-
τ
model to a bias under 0.1 dex for the nonparametric model, though this is heavily dependent on the choice of prior for the nonparametric model. Similarly, using nonparametric SFHs in SED fitting results in increased accuracy in recovered galaxy star formation rates and stellar ages.
Parametric models for galaxy star formation histories (SFHs) are widely used, though they are known to impose strong priors on physical parameters. This has consequences for measurements of the ...galaxy stellar-mass function, star formation rate density (SFRD), and star-forming main sequence (SFMS). We investigate the effects of the exponentially declining, delayed exponentially declining, lognormal, and double power-law SFH models using Bagpipes. We demonstrate that each of these models imposes strong priors on specific star formation rates (SFRs), potentially biasing the SFMS, and also imposes a strong prior preference for young stellar populations. We show that stellar mass, SFR, and mass-weighted age inferences from high-quality mock photometry vary with the choice of SFH model by at least 0.1, 0.3, and 0.2 dex, respectively. However, the biases with respect to the true values depend more on the true SFH shape than the choice of model. We also demonstrate that photometric data cannot discriminate between SFH models, meaning that it is important to perform independent tests to find well-motivated priors. We finally fit a low-redshift, volume-complete sample of galaxies from the Galaxy and Mass Assembly (GAMA) Survey with each model. We demonstrate that our stellar masses and SFRs at redshift z ∼ 0.05 are consistent with other analyses. However, our inferred cosmic SFRDs peak at z ∼ 0.4, approximately 6 Gyr later than direct observations suggest, meaning that our mass-weighted ages are significantly underestimated. This makes the use of parametric SFH models for understanding mass assembly in galaxies challenging. In a companion paper, we consider nonparametric SFH models.
Forward modeling of the full galaxy SED is a powerful technique, providing self-consistent constraints on stellar ages, dust properties, and metallicities. However, the accuracy of these results is ...contingent on the accuracy of the model. One significant source of uncertainty is the contribution of obscured AGN, as they are relatively common and can produce substantial mid-IR (MIR) emission. Here we include emission from dusty AGN torii in the Prospector SED-fitting framework, and fit the UV-IR broadband photometry of 129 nearby galaxies. We find that 10% of the fitted galaxies host an AGN contributing >10% of the observed galaxy MIR luminosity. We demonstrate the necessity of this AGN component in the following ways. First, we compare observed spectral features to spectral features predicted from our model fit to the photometry. We find that the AGN component greatly improves predictions for observed H and Hβ luminosities, as well as mid-infrared Akari and Spitzer/IRS spectra. Second, we show that inclusion of the AGN component changes stellar ages and SFRs by up to a factor of 10, and dust attenuations by up to a factor of 2.5. Finally, we show that the strength of our model AGN component correlates with independent AGN indicators, suggesting that these galaxies truly host AGN. Notably, only 46% of the SED-detected AGN would be detected with a simple MIR color selection. Based on these results, we conclude that SED models which fit MIR data without AGN components are vulnerable to substantial bias in their derived parameters.
Gravitational interactions between the Large Magellanic Cloud (LMC) and the stellar and dark matter halo of the Milky Way are expected to give rise to disequilibrium phenomena in the outer Milky Way
.... A local wake is predicted to trail the orbit of the LMC, and a large-scale overdensity is predicted to exist across a large area of the northern Galactic hemisphere. Here we report the detection of both the local wake and northern overdensity (hereafter the 'collective response') in a map of the Galaxy based on 1,301 stars at Galactocentric distances between 60 and 100 kiloparsecs. The location of the wake is in good agreement with an N-body simulation that includes the dynamical effect of the LMC on the Milky Way halo. The density contrast of the wake and collective response are stronger in the data than in the simulation. The detection of a strong local wake is independent evidence that the Magellanic clouds are on their first orbit around the Milky Way. The wake traces the path of the LMC, which will provide insight into the orbit of the LMC, which in turn is a sensitive probe of the mass of the LMC and the Milky Way. These data demonstrate that the outer halo is not in dynamical equilibrium, as is often assumed. The morphology and strength of the wake could be used to test the nature of dark matter and gravity.