skip to main content


Title: In pursuit of giants: II. Evolution of dusty quiescent galaxies over the last six billion years from the hCOSMOS survey

The physical mechanisms that link the termination of star formation in quiescent galaxies and the evolution of their baryonic components, stars, and the interstellar medium (ISM; dust, gas, and metals) are poorly constrained beyond the local Universe. In this work, we characterise the evolution of the dust content in 545 quiescent galaxies observed at 0.1 < z < 0.6 as part of the hCOSMOS spectroscopic redshift survey. This is, to date, the largest sample of quiescent galaxies at intermediate redshifts for which the dust, stellar, and metal abundances are consistently estimated. We analyse how the crucial markers of a galaxy dust life cycle, such as specific dust mass (Mdust/M), evolve with different physical parameters, namely gas-phase metallicity (Zgas), time since quenching (tquench), stellar mass (M), and stellar population age. We find morphology to be an important factor in the large scatter inMdust/M(∼2 orders of magnitude). Quiescent spirals exhibit strong evolutionary trends of specific dust mass withM, stellar age, and galaxy size, in contrast to the little to no evolution experienced by ellipticals. When transitioning from solar to super-solar metallicities (8.7 ≲ 12 + log(O/H)≲9.1), quiescent spirals undergo a reversal inMdust/M, indicative of a change in dust production efficiency. By modelling the star formation histories of our objects, we unveil a broad dynamical range of post-quenching timescales (60 Myr < tquench < 3.2 Gyr). We show thatMdust/Mis highest in recently quenched systems (tquench < 500 Myr), but its further evolution is non-monotonic, as a consequence of different pathways for dust formation, growth, or removal on various timescales. Our data are best described by simulations that include dust growth in the ISM. While this process is prevalent in the majority of galaxies, for ∼15% of objects we find evidence of additional dust content acquired externally, most likely via minor mergers. Altogether, our results strongly suggest that prolonged dust production on a timescale of 0.5 − 1 Gyr since quenching may be common in dusty quiescent galaxies at intermediate redshifts, even if their gas reservoirs are heavily exhausted (i.e. cold gas fraction < 1 − 5%).

 
more » « less
Award ID(s):
1909153
NSF-PAR ID:
10532002
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
A&A
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
678
ISSN:
0004-6361
Page Range / eLocation ID:
A35
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Single flux density measurements at observed-frame submillimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive halos from the Feedback in Realistic Environments project. Our starting point is four star-forming central galaxies at seven redshifts betweenz= 1.5 andz= 4.5. We generate modified quiescent galaxies that have been quenched for 100 Myr, 500 Myr, or 1 Gyr prior to each of the studied redshifts by reassigning stellar ages. We derive spectral energy distributions for each fiducial and modified galaxy using radiative transfer. We demonstrate that the dust mass inferred is highly dependent on the assumed dust temperature,Tdust, which is often unconstrained observationally. Motivated by recent work on quiescent galaxies that assumedTdust∼ 25 K, we show that the ratio between dust mass and 1.3 mm flux density can be higher than inferred by up to an order of magnitude, due to the considerably lower dust temperatures seen in non-star-forming galaxies. This can lead to an underestimation of dust mass (and, when submillimeter flux density is used as a proxy for molecular gas content and gas mass). This underestimation is most severe at higher redshifts, where the observed-frame 1.3 mm flux density probes rest-frame wavelengths far from the Rayleigh–Jeans regime, and hence depends superlinearly on dust temperature. We fit relations between ratios of rest-frame far-infrared flux densities and mass-weighted dust temperature that can be used to constrain dust temperatures from observations and hence derive more reliable dust and molecular gas masses.

     
    more » « less
  2. Abstract

    Over the past decade, rest-frame color–color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color–color selection tool—the rest-frameUVversusVJdiagram—by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties inUVJspace and the color–color evolution of massive galaxies as they quench at redshiftsz∼ 1–2. We find that our models broadly reproduce the observedUVJdiagram atz= 1–2, including (for the first time in a cosmological simulation) reproducing the population of extremely dust-reddened galaxies in the top right of theUVJdiagram. However, our models primarily populate this region with low-mass galaxies and do not produce as clear a bimodality between star-forming and quiescent galaxies as is seen in observations. The former issue is due to an excess of dust in low-mass galaxies and relatively gray attenuation curves in high-mass galaxies, while the latter is due to the overpopulation of the green valley insimba. When investigating the time evolution of galaxies on theUVJdiagram, we find that the quenching pathway on theUVJdiagram is independent of the quenching timescale, and instead dependent primarily on the average specific star formation rate in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies.

     
    more » « less
  3. ABSTRACT

    The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($M_{*}\lesssim 10^7\, \rm {M}_{\odot }$) are mostly quiescent and higher mass satellites ($M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$) are mostly star forming, with intermediate-mass satellites ($M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced ≥1–2 pericentre passages (≈2.5–5 Gyr) within a MW-mass halo. Most galaxies with $M_{*}\gtrsim 10^{6.5}\, \rm {M}_{\odot }$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the Satellites Around Galactic Analogs (SAGA) survey’s quiescent fraction at $M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$, they do not generally reproduce SAGA’s turnover at lower masses.

     
    more » « less
  4. Abstract

    Due to their extremely dust-obscured nature, much uncertainty still exists surrounding the stellar mass growth and content in dusty, star-forming galaxies (DSFGs) atz> 1. In this work, we present a numerical model built using empirical data on DSFGs to estimate their stellar mass contributions across the first ∼10 Gyr of cosmic time. We generate a dust-obscured stellar mass function that extends beyond the mass limit of star-forming stellar mass functions in the literature, and predict that massive DSFGs constitute as much as 50%–100% of all star-forming galaxies withM≥1011Matz> 1. We predict the number density of massive DSFGs and find general agreement with observations, although more data is needed to narrow wide observational uncertainties. We forward-model mock massive DSFGs to their quiescent descendants and find remarkable agreement with observations from the literature demonstrating that, to first order, massive DSFGs are a sufficient ancestral population to describe the prevalence of massive quiescent galaxies atz> 1. We predict that massive DSFGs and their descendants contribute as much as 25%–60% to the cosmic stellar mass density during the peak of cosmic star formation, and predict an intense epoch of population growth during the ∼1 Gyr fromz= 6 to 3 during which the majority of the most massive galaxies at high-zgrow and then quench. Future studies seeking to understand massive galaxy growth and evolution in the early universe should strategize synergies with data from the latest observatories (e.g., JWST and the Atacama Large Millimeter/submillimeter Array) to better include the heavily dust-obscured galaxy population.

     
    more » « less
  5. Abstract Using stellar population synthesis models to infer star formation histories (SFHs), we analyze photometry and spectroscopy of a large sample of quiescent galaxies that are members of Sunyaev–Zel’dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < z < 1.4 using rest-frame optical spectra and the Python-based Prospector framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 < z < 0.9) and three clusters in the SPT Hi-z cluster sample (1.25 < z < 1.4). We analyze spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our data set sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass, t 50 ) of z = 2.8 ± 0.5, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age–stellar mass relations for the full sample (age of the universe at t 50 (Gyr) = 2.52 (±0.04)–1.66 (±0.12) log 10 ( M /10 11 M ⊙ )) and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate ∼0.75 Gyr earlier than lower-mass galaxies. We also find marginally steeper age–mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass. 
    more » « less