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).
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.
ABSTRACT
The physical properties of galaxies are encoded within their spectral energy distribution and require comparison with models to be extracted. These models must contain a synthetic stellar ...population and, where infrared data are to be used, also consider prescriptions for energy reprocessing and re-emission by dust. While many such models have been constructed, there are few analyses of the impact of stellar population model choice on derived dust parameters, or vice versa. Here, we apply a simple framework to compare the impact of these choices, combining three commonly used stellar population synthesis models and three dust emission models. We compare fits to the ultraviolet to far-infrared spectral energy distributions of a validation sample of infrared-luminous galaxies. We find that including different physics, such as binary stellar evolution, in the stellar synthesis model can introduce biases and uncertainties in the derived parameters of the dust and stellar emission models, largely due to differences in the far-ultraviolet emission available for reprocessing. This may help to reconcile the discrepancy between the cosmic star formation rate and stellar mass density histories. Notably the inclusion of a dusty stellar birth cloud component in the dust emission model provides more flexibility in accommodating the stellar population model, as its re-emission is highly sensitive to the ultraviolet radiation field spectrum and density. Binary populations favour a longer birth cloud dissipation time-scale than is found when assuming only single star population synthesis.
Abstract
We present constraints on the physical properties (including stellar mass, age, and star formation rate) of 207 6 ≲
z
≲ 8 galaxy candidates from the Reionization Lensing Cluster Survey ...(RELICS) and Spitzer-RELICS surveys. We measure photometry using T-PHOT and perform spectral energy distribution fitting using EA
z
Y and BAGPIPES. Of the 207 candidates for which we could successfully measure (or place limits on) Spitzer fluxes, 23 were demoted to likely
z
< 4. Among the high-
z
candidates, we find intrinsic stellar masses between 1 × 10
6
M
⊙
and 4 × 10
9
M
⊙
, and rest-frame UV absolute magnitudes between −22.6 and −14.5 mag. While our sample is mostly comprised of
L
m
UV
/
L
m
UV
*
<
1
galaxies, it extends to
L
m
UV
/
L
m
UV
*
∼
2
. Our sample spans ∼4 orders of magnitude in stellar mass and star formation rates, and exhibits ages that range from maximally young to maximally old. We highlight 11
z
≥ 6.5 galaxies with detections in Spitzer/IRAC imaging, several of which show evidence for some combination of evolved stellar populations, large contributions of nebular emission lines, and/or dust. Among these is PLCKG287+32-2013, one of the brightest
z
∼ 7 candidates known (AB mag 24.9 at 1.6
μ
m) with a Spitzer 3.6
μ
m flux excess suggesting strong O
iii
+ H-
β
emission (∼1000 Å rest-frame equivalent width). We discuss the possible uses and limits of our sample and present a public catalog of Hubble + Spitzer photometry along with physical property estimates for all objects in the sample. Because of their apparent brightnesses, high redshifts, and variety of stellar populations, these objects are excellent targets for follow-up with the James Webb Space Telescope.
We report on the Hubble Space Telescope (HST) detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z = 3.794 dubbed Ion1. The LyC from Ion1 is detected at 820−890 with HST ...WFC3/UVIS in the F410M band (m410 = 27.60 0.36 mAB, peak signal-to-noise ratio (S/N) = 4.17 in an r = 0 12 aperture) and 700−830 with the Very Large Telescope (VLT)/VIMOS in the U band (mU = 27.84 0.19 mAB, peak S/N = 6.7 with an r = 0 6 aperture). A 20 hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines and the P Cygni profile of C iv and Ly in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the nonionizing UV emission by 0 12 0 03, corresponding to 0.85 0.21 kpc, and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low H i column density "cavities" in the interstellar medium. We also collect the VIMOS U-band images of 107 Lyman-break galaxies at 3.40 < zspec < 3.95, i.e., sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 mAB flux limit (1 ) and an upper limit of absolute LyC escape fraction fescabs ≤ 0.63%.
Abstract
The study of galaxy evolution hinges on our ability to interpret multiwavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy ...distribution (SED) models, which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to wide and deep multiwave band galaxy surveys, the volume of high-quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED-fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply 14 of the most commonly used SED-fitting codes on samples from the CANDELS photometric catalogs at
z
∼ 1 and
z
∼ 3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust-attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust and active galactic nucleus models) on the derived stellar masses, SFRs, and
A
V
values. We then assess the difference among the codes on the SFR–stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (∼0.1 dex), SFR (∼0.3 dex), and dust attenuation (∼0.3 mag). Finally, we present some resources summarizing best practices in SED fitting.
The extremely rapid assembly of the earliest galaxies during the first billion years of cosmic history is a major challenge for our understanding of galaxy formation physics
. The advent of the James ...Webb Space Telescope (JWST) has exacerbated this issue by confirming the existence of galaxies in substantial numbers as early as the first few hundred million years
. Perhaps even more surprisingly, in some galaxies, this initial highly efficient star formation rapidly shuts down, or quenches, giving rise to massive quiescent galaxies as little as 1.5 billion years after the Big Bang
. However, due to their faintness and red colour, it has proven extremely challenging to learn about these extreme quiescent galaxies, or to confirm whether any existed at earlier times. Here we report the spectroscopic confirmation of a massive quiescent galaxy, GS-9209, at redshift, z = 4.658, just 1.25 billion years after the Big Bang, using the JWST Near-Infrared Spectrograph (NIRSpec). From these data we infer a stellar mass of M
= 3.8 ± 0.2 × 10
M
, which formed over a roughly 200 Myr period before this galaxy quenched its star-formation activity at Formula: see text, when the Universe was approximately 800 Myr old. This galaxy is both a likely descendent of the highest-redshift submillimetre galaxies and quasars, and a likely progenitor for the dense, ancient cores of the most massive local galaxies.
Abstract
We present a rest-frame UV–optical (
λ
= 2500–6400 Å) stacked spectrum representative of massive quiescent galaxies at 1.0 <
z
< 1.3 with log(
M
*
/
M
⊙
) > 10.8. The stack is constructed ...using VANDELS survey data, combined with new KMOS observations. We apply two independent full-spectral-fitting approaches, measuring a total metallicity Z/H = −0.13 ± 0.08 with
Bagpipes
and Z/H = 0.04 ± 0.14 with
Alf
, a fall of ∼0.2–0.3 dex compared with the local universe. We also measure an iron abundance Fe/H = −0.18 ± 0.08, a fall of ∼0.15 dex compared with the local universe. We measure the alpha enhancement via the magnesium abundance, obtaining Mg/Fe = 0.23 ± 0.12, consistent with galaxies of similar mass in the local universe, indicating no evolution in the average alpha enhancement of log(
M
*
/
M
⊙
) ∼ 11 quiescent galaxies over the last ∼8 Gyr. This suggests the very high alpha enhancements recently reported for several bright
z
∼ 1–2 quiescent galaxies are due to their extreme masses, log(
M
*
/
M
⊙
) ≳ 11.5, in accordance with the well-known downsizing trend, rather than being typical of the
z
≳ 1 population. The metallicity evolution we observe with redshift (falling Z/H, Fe/H, constant Mg/Fe) is consistent with recent studies. We recover a mean stellar age of
2.5
−
0.4
+
0.6
Gyr, corresponding to a formation redshift
z
form
=
2.4
−
0.3
+
0.6
. Recent studies have obtained varying average formation redshifts for
z
≳ 1 massive quiescent galaxies, and, as these studies report consistent metallicities, we identify models with different star formation histories as the most likely cause. Larger spectroscopic samples from upcoming ground-based instruments will provide precise constraints on ages and metallicities at
z
≳ 1. Combining these with precise stellar mass functions for
z
> 2 quiescent galaxies from the James Webb Space Telescope will provide an independent test of formation redshifts derived from spectral fitting.
Abstract
MACS0647–JD is a triply lensed
z
∼ 11 galaxy originally discovered with the Hubble Space Telescope. The three lensed images are magnified by factors of ∼8, 5, and 2 to AB mag 25.1, 25.6, and ...26.6 at 3.5
μ
m. The brightest is over a magnitude brighter than other galaxies recently discovered at similar redshifts
z
> 10 with JWST. Here, we report new JWST imaging that clearly resolves MACS0647–JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. The brighter larger component “A” is intrinsically very blue (
β
∼ −2.6 ± 0.1), likely due to very recent star formation and no dust, and is spatially extended with an effective radius ∼70 ± 24 pc. The smaller component “B” (
r
∼ 20
−
5
+
8
pc) appears redder (
β
∼ −2 ± 0.2), likely because it is older (100–200 Myr) with mild dust extinction (
A
V
∼ 0.1 mag). With an estimated stellar mass ratio of roughly 2:1 and physical projected separation ∼400 pc, we may be witnessing a galaxy merger 430 million years after the Big Bang. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be dissimilar, which is also suggested by the spectral energy distribution fitting, suggesting they formed further apart. We also identify a candidate companion galaxy “C” ∼3 kpc away, likely destined to merge with A and B. Upcoming JWST Near Infrared Spectrograph observations planned for 2023 January will deliver spectroscopic redshifts and more physical properties for these tiny magnified distant galaxies observed in the early universe.