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-frameU−VversusV−Jdiagram—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
IQ Collaboratory. III. The Empirical Dust Attenuation Framework—Taking Hydrodynamical Simulations with a Grain of Dust
Abstract We present the empirical dust attenuation (EDA) framework—a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color–magnitude relations, ( g − r ) − M r and (far-UV −near-UV) − M r , of the simulations to a M r < − 20 and UV complete Sloan Digital Sky Survey galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color–magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation–slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low A V since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations.
more »
« less
- Award ID(s):
- 1831412
- PAR ID:
- 10333842
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 926
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 122
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a “normal” molecular-gas to dust mass ratio, δ GDR , typically of order 100. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust have been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme δ GDR ratios. We show that simulated massive quiescent galaxies at 0 < z < 3 in the simba cosmological simulations have δ GDR values that extend >4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of >3 orders of magnitude in dust mass to return to normal values of δ GDR despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly preselect exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities.more » « less
-
ABSTRACT We present predictions for high redshift (z = 2−10) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the H α and H β + $$[\rm O \,{\small III}]$$ luminosity functions up to z = 8. The predicted H β + $$[\rm O \,{\small III}]$$ luminosity functions are consistent with present observations at z ≲ 3 with $${\lesssim} 0.1\, {\rm dex}$$ differences in luminosities. However, the predicted H α luminosity function is $${\sim }0.3\, {\rm dex}$$ dimmer than the observed one at z ≃ 2. Furthermore, we explore continuum spectral indices, the Balmer break at 4000 Å; (D4000) and the UV continuum slope β. The median D4000 versus specific star formation rate relation predicted at z = 2 is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed AUV versus β relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at z = 2−6 and investigate their dependence on galaxy colours and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colours, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high-redshift galaxies. Future James Webb Space Telecope observations will further test these predictions.more » « less
-
Abstract We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at 1.3 ≤z≤ 2.6 from the MOSFIRE Deep Evolution Field survey. We divide galaxies with a detected Hαemission line and coverage of Hβinto eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), and we then stack their spectra. From each stack, we measure the Balmer decrement and gas-phase metallicity, and then we compute the medianAVand UV continuum spectral slope (β). First, we find that none of the dust properties (Balmer decrement,AV, orβ) varies with the axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model in which attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases the SFR and lowers the metallicity, but leaves the dust column density mostly unchanged. We quantify this idea using the Kennicutt–Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation.more » « less
-
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