Star-formation histories of massive quiescent galaxies

This thesis presents several related analyses designed to understand the star-formation histories (SFHs) and quenching mechanisms of massive quiescent galaxies across cosmic time. More generally, it contains research directed at sophisticated modelling and Bayesian fitting of galaxy spectra. I firstly present Bayesian Analysis of Galaxies for Physical Inference and Parameter EStimation, or Bagpipes, a new, publicly available Python code that can be used to rapidly generate complex model galaxy spectra and to fit these to arbitrary combinations of spectroscopic and photometric data. I then perform a detailed analysis of the SFHs of a sample of 9289 quiescent galaxies from UltraVISTA with stellar masses, M∗> 1010M⊙ and observed redshifts from 0:25 < z < 3:75. The majority of these galaxies exhibit SFHs that rise gradually then quench relatively rapidly, over 1-2 Gyr. This behaviour is consistent with recent cosmological hydrodynamic simulations, where AGN-driven feedback in the low-accretion (jet) mode is the dominant quenching mechanism. At z > 1, I also find a class of objects with SFHs that rise and fall very rapidly, with quenching timescales of < 1 Gyr, consistent with quasar-mode AGN feedback. Finally, at z < 1, I find a population with SFHs that quench more slowly than they rise, over > 3 Gyr, consistent with other such analyses in the local Universe. I confirm the trend towards earlier formation with increasing stellar mass (downsizing) at fixed observed redshift, and a trend towards more rapid quenching at higher stellar masses. I then present a general investigation of the use of parametric SFH models in spectral fitting analyses. Parametric models for galaxy SFHs are widely used, though they are known to impose strong priors on physical parameters, with consequences for measurements of the galaxy stellar-mass function, star-formation- rate density (SFRD) and star-forming main sequence (SFMS). I investigate the effects of the exponentially declining, delayed exponentially declining, lognormal and double power law SFH models. I demonstrate that each of these models imposes strong priors on specific star-formation rates (sSFRs), potentially biasing the SFMS, and also imposes a strong prior preference for young stellar populations. I 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. I also demonstrate that photometric data cannot discriminate between SFH models, meaning it is important to perform independent tests to find well-motivated priors. In response to this I finally fit a low-redshift, volume-complete sample from the Galaxy and Mass Assembly (GAMA) Survey with each model. I demonstrate that the inferred stellar masses and SFRs at redshift, z ~ 0:05 are consistent with other analyses. However, the inferred cosmic SFRDs peak at z ~ 0:4, approximately 6 Gyr later than direct observations suggest, meaning that mass-weighted ages are significantly underestimated. This makes the use of parametric SFH models for understanding mass assembly in galaxies challenging. I finally present a Bayesian full-spectral-fitting analysis of 75 massive (M∗> 1010:3M⊙) UVJ-selected galaxies at redshifts of 1:0 < z < 1:3, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry by the use of a sophisticated physical plus systematic uncertainties model. I constrain the stellar mass vs stellar age relationship, finding a strong trend towards earlier formation with increasing stellar mass (downsizing) of 1:48+0:34 ≲0:39 Gyr per decade in mass. I show that this is consistent with other spectroscopic studies from 0 < z < 2. This places strong constraints on the AGN-feedback models used in cosmological simulations. I demonstrate that, although the relationships predicted by the Simba and IllustrisTNG simulations agree well with observations at z = 0:1, they are too shallow at z = 1, predicting an evolution of . 0:5 Gyr per decade in mass. The majority of the lowest-mass galaxies in the sample (M∗~ 1010:5M⊙) are consistent with formation in recent (z < 2), intense starburst events, with timescales of ≲ 500 Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in the sample are extreme systems: already old by z = 1, they formed at z ~ 5 and quenched by z = 3. However, I find evidence for their continued evolution through both AGN and rejuvenated star-formation activity. To understand the detailed SFHs of these objects, similar studies must be extended to the highest redshifts

Introducing a real-time interactive GUI tool for visualization of galaxy spectra

To aid the understanding of the nonlinear relationship between galaxy properties and predicted spectral energy distributions (SED), we present a new interactive graphical user interface tool pipes_vis based on Bagpipes. It allows for real-time manipulation of a model galaxy's star formation history, dust and other relevant properties through sliders and text boxes, with each change's effect on the predicted SED reflected instantaneously. We hope the tool will assist in building intuition about what affects the SED of galaxies, potentially helping to speed up fitting stages such as prior construction, and aid in undergraduate and graduate teaching. pipes_vis is available online (pipes_vis is maintained and documented online at https://github.com/HinLeung622/pipes_vis, or version 0.4.1 is archived in Zenodo (Leung 2021) and also available for installation through pip install pipes_vis).PostprintNon peer reviewe

How to Measure Galaxy Star Formation Histories. II. Nonparametric Models

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 {\it 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 \texttt{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 (SFR) as a function of time is the choice of prior, rather than the photometric noise. As a practical demonstration, we fit the UV-IR photometry of $\sim$6000 galaxies from the GAMA survey and measure inter-prior scatters in mass (0.1 dex), SFR$_{100\; \mathrm{Myr}}$ (0.8 dex), and mass-weighted ages (0.2 dex), 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).Comment: replacing with ApJ accepted versio

Chemical evolution of local post-starburst galaxies: Implications for the mass-metallicity relation

We use the stellar fossil record to constrain the stellar metallicity evolution and star-formation histories of the post-starburst regions within 45 local post-starburst galaxies from the MaNGA survey. The direct measurement of the regions' stellar metallicity evolution is achieved by a new two-step metallicity model that allows for stellar metallicity to change at the peak of the starburst. We also employ a Gaussian process noise model that accounts for correlated errors introduced by the observational data reduction or inaccuracies in the models. We find that a majority of post-starburst regions (69% at $>1\sigma$ significance) increased in stellar metallicity during the recent starburst, with an average increase of 0.8 dex and a standard deviation of 0.4 dex. A much smaller fraction of PSBs are found to have remained constant (22%) or declined in metallicity (9%, average decrease 0.4 dex, standard deviation 0.3 dex). The pre-burst metallicities of the post-starburst galaxies are in good agreement with the mass-metallicity relation of local star-forming galaxies. These results are consistent with hydrodynamic simulations, which suggest that mergers between gas-rich galaxies are the primary formation mechanism of local PSBs, and rapid metal recycling during the starburst outweighs the impact of dilution by any gas inflows. The final mass-weighted metallicities of the post-starburst galaxies are consistent with the mass-metallicity relation of local passive galaxies. Our results suggest that rapid quenching following a merger-driven starburst is entirely consistent with the observed gap between the stellar mass-metallicity relations of local star-forming and passive galaxies.Comment: 18+4 pages, 8+2 figures, submitted to MNRA

Chemical evolution of local post-starburst galaxies : implications for the mass-metallicity relation

We use the stellar fossil record to constrain the stellar metallicity evolution and star-formation histories of the post-starburst (PSB) regions within 45 local PSB galaxies from the MaNGA survey. The direct measurement of the regions’ stellar metallicity evolution is achieved by a new two-step metallicity model that allows for stellar metallicity to change at the peak of the starburst. We also employ a Gaussian process noise model that accounts for correlated errors introduced by the observational data reduction or inaccuracies in the models. We find that a majority of PSB regions (69 per cent at >1σ significance) increased in stellar metallicity during the recent starburst, with an average increase of 0.8 dex and a standard deviation of 0.4 dex. A much smaller fraction of PSBs are found to have remained constant (22 per cent) or declined in metallicity (9 per cent, average decrease 0.4 dex, standard deviation 0.3 dex). The pre-burst metallicities of the PSB galaxies are in good agreement with the mass–metallicity (MZ) relation of local star-forming galaxies. These results are consistent with hydrodynamic simulations, which suggest that mergers between gas-rich galaxies are the primary formation mechanism of local PSBs, and rapid metal recycling during the starburst outweighs the impact of dilution by any gas inflows. The final mass-weighted metallicities of the PSB galaxies are consistent with the MZ relation of local passive galaxies. Our results suggest that rapid quenching following a merger-driven starburst is entirely consistent with the observed gap between the stellar mass–metallicity relations of local star-forming and passive galaxies.Peer reviewe

A census of star formation histories of massive galaxies at 0.6 < z < 1 from spectro-photometric modeling using Bagpipes and Prospector

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($\tau_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe.Comment: 25 pages, 12 figures, 4 tables, submitted to Ap

How to Measure Galaxy Star Formation Histories. I. Parametric Models

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 (GSMF), 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 (sSFRs), 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 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\sim0.05$ are consistent with other analyses. However, our inferred cosmic SFRDs peak at $z\sim0.4$, approximately 6 Gyr later than direct observations suggest, meaning 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 non-parametric SFH models.Comment: 20 pages, 12 figures, ApJ accepte

The star formation histories of z ∼ 1 post-starburst galaxies

Funding: LTA acknowledges support from the Ministry of Higher Education and Scientific Research (MOHESR), Iraq. AW acknowledges financial support from the Royal Society Newton Fund (grant NAF/R1/180403, PI Natalia Vale Asari) and Fundação de à Amparo Pesquisa do Estado de São Paulo (FAPESP) process number 2019/01768-6.We present the star formation histories of 39 galaxies with high-quality rest-frame optical spectra at 0.5 <z <1.3 selected to have strong Balmer absorption lines and/or Balmer break, and compare to a sample of spectroscopically selected quiescent galaxies at the same redshift. Photometric selection identifies a majority of objects that have clear evidence for a recent short-lived burst of star formation within the last 1.5 Gyr, i.e. 'post-starburst' galaxies, however we show that good quality continuum spectra are required to obtain physical parameters such as burst mass fraction and burst age. Dust attenuation appears to be the primary cause for misidentification of post-starburst galaxies, leading to contamination in spectroscopic samples where only the [O II] emission line is available, as well as a small fraction of objects lost from photometric samples. The 31 confirmed post-starburst galaxies have formed 40-90 per cent of their stellar mass in the last1-1.5 Gyr. We use the derived star formation histories to find that the post-starburst galaxies are visible photometrically for 0.5-1 Gyr. This allows us to update a previous analysis to suggest that 25-50 per cent of the growth of the red sequence at z ∼ 1 could be caused by a starburst followed by rapid quenching. We use the inferred maximum historical star formation rates of several 100-1000 M⊙yr-1 and updated visibility times to confirm that sub-mm galaxies are likely progenitors of post-starburst galaxies. The short quenching time-scales of 100-200 Myr are consistent with cosmological hydrodynamic models in which rapid quenching is caused by the mechanical expulsion of gas due to an acive galactic neucleus.Publisher PDFPeer reviewe

Remarkably Compact Quiescent Candidates at 3<z<53<z<5 in JWST-CEERS

In this letter, we measure the rest-frame optical and near-infrared sizes of ten quiescent candidates at $3<z<5$, first reported by Carnall et al. (2023a). We use James Webb Space Telescope (JWST) Near-Infrared Camera (NIRCam) F277W and F444W imaging obtained through the public CEERS Early Release Science (ERS) program and imcascade, an astronomical fitting code that utilizes Multi-Gaussian Expansion, to carry out our size measurements. When compared to the extrapolation of rest-optical size-mass relations for quiescent galaxies at lower redshift, eight out of ten candidates in our sample (80%) are on average more compact by $\sim$40%. Seven out of ten candidates (70%) exhibit rest-frame infrared sizes $\sim$10% smaller than rest-frame optical sizes, indicative of negative color gradients. Two candidates (20%) have rest-frame infrared sizes $\sim$1.4$\times$ larger than rest-frame optical sizes; one of these candidates exhibits signs of ongoing or residual star formation, suggesting this galaxy may not be fully quenched. The remaining candidate is unresolved in both filters, which may indicate an Active Galactic Nuclei (AGN). Strikingly, we observe three of the most massive galaxies in the sample (log(M$_{\star}$/M$_{\odot}$) = 10.74 - 10.95) are extremely compact, with effective radii ${\sim}$0.7 kpc. Our findings provide no indication that the size evolution relation flattens out, and may indicate that the size evolution of quiescent galaxies is steeper than previously anticipated beyond $z>3$.Comment: Accepted for publication in ApJL. 11 pages, 4 figures, 1 tabl

The star formation histories of z~1 post-starburst galaxies

We present the star formation histories of 39 galaxies with high quality rest-frame optical spectra at 0.5<z<1.3 selected to have strong Balmer absorption lines and/or Balmer break, and compare to a sample of spectroscopically selected quiescent galaxies at the same redshift. Photometric selection identifies a majority of objects that have clear evidence for a recent short-lived burst of star formation within the last 1.5 Gyr, i.e. "post-starburst" galaxies, however we show that good quality continuum spectra are required to obtain physical parameters such as burst mass fraction and burst age. Dust attenuation appears to be the primary cause for misidentification of post-starburst galaxies, leading to contamination in spectroscopic samples where only the [OII] emission line is available, as well as a small fraction of objects lost from photometric samples. The 31 confirmed post-starburst galaxies have formed 40-90% of their stellar mass in the last 1-1.5 Gyr. We use the derived star formation histories to find that the post-starburst galaxies are visible photometrically for 0.5-1 Gyr. This allows us to update a previous analysis to suggest that 25-50% of the growth of the red sequence at z~1 could be caused by a starburst followed by rapid quenching. We use the inferred maximum historical star formation rates of several 100-1000 Msun/yr and updated visibility times to confirm that sub-mm galaxies are likely progenitors of post-starburst galaxies. The short quenching timescales of 100-200 Myr are consistent with cosmological hydrodynamic models in which rapid quenching is caused by the mechanical expulsion of gas due to an AGN.Comment: 24 pages, 16 figures, MNRAS accepted 06/03/202