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
How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-dust Geometry Model in Spectral Energy Distribution Fitting
Abstract One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the simba cosmological simulation with the prospector SED fitting code. We consider the impact of the commonly assumed uniform screen model and introduce a new nonuniform screen model parameterized by the fraction of unobscured stellar light. This nonuniform screen model allows for a nonzero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a nonuniform screen model, with the median offset in UV attenuation decreasing from −0.30 dex with a uniform screen model to −0.17 dex with the nonuniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the nonuniform screen model, decreasing the SFR offset on average by 0.12 dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.
more »
« less
- Award ID(s):
- 1909153
- NSF-PAR ID:
- 10347523
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 931
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 14
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
more » « less
ABSTRACT We present an analysis of the dust attenuation of star-forming galaxies at z = 2.5–4.0 through the relationship between the UV spectral slope (β), stellar mass (M*), and the infrared excess (IRX = LIR/LUV) based on far-infrared continuum observations from the Atacama Large Millimeter/sub-millimeter Array (ALMA). Our study exploits the full ALMA archive over the COSMOS field processed by the A3COSMOS team, which includes an unprecedented sample of ∼1500 galaxies at z ∼ 3 as primary or secondary targets in ALMA band 6 or 7 observations with a median continuum sensitivity of 126 $\rm {\mu Jy\, beam}^{-1}$ (1σ). The detection rate is highly mass dependent, decreasing drastically below log (M*/M⊙) = 10.5. The detected galaxies show that the IRX–β relationship of massive (log M*/M⊙ > 10) main-sequence galaxies at z = 2.5–4.0 is consistent with that of local galaxies, while starbursts are generally offset by $\sim 0.5\, {\rm dex}$ to larger IRX values. At the low-mass end, we derive upper limits on the infrared luminosities through stacking of the ALMA data. The combined IRX–M* relation at $\rm {log\, ({\it M}_{\ast }/\mathrm{M}_{\odot })\gt 9}$ exhibits a significantly steeper slope than reported in previous studies at similar redshifts, implying little dust obscuration at log M*/M⊙ < 10. However, our results are consistent with earlier measurements at z ∼ 5.5, indicating a potential redshift evolution between z ∼ 2 and z ∼ 6. Deeper observations targeting low-mass galaxies will be required to confirm this finding.
-
more » « less
ABSTRACT We present 65 extremely dust-obscured galaxies from the UltraVISTA DR3 survey of the COSMOS field at 1 < z < 4. In contrast to other studies of dusty galaxies, we select our sample based on dust attenuation measured by UV–MIR spectral energy distribution (SED) modelling that allows for extreme attenuation levels. We construct our sample by making cuts at 1 ≤ z ≤ 4, AV ≥ 3, and log(M*/M⊙) ≥ 10.5. This method reliably selects galaxies exhibiting independent indicators of significant dust content, including far-infrared detection rates. We perform panchromatic SED modelling with matched Herschel photometry and find stellar and dust properties that differ from typical submillimetre galaxy (SMG) samples as well as Herschel sources matched in redshift and stellar mass. Our sources have lower star formation rates and higher AV than SMGs, but comparable total infrared luminosities. Most of our sample falls on or near the star-forming main sequence for this redshift range. Finally, we perform a morphological analysis with galfit using the KS-band images and Hubble F814W and F160W imaging when available. Typical axial ratios of ∼0.4 suggest disc-like morphology for the majority of our sources, and we note only three apparent merging systems. Our sample generally agrees with the size–mass relation for star-forming galaxies, with a tail extending to smaller sizes. We conclude that the most heavily obscured galaxies in this redshift range share many characteristics with typical star-forming galaxies, forming a population of dusty galaxies that overlaps, but is not encompassed by, those selected through dust emission.
-
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 medianA Vand UV continuum spectral slope (β ). First, we find that none of the dust properties (Balmer decrement,A V, 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.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4357/ac6959
