The Dust Attenuation Law in Galaxies Salim, Samir; Narayanan, Desika
Annual review of astronomy and astrophysics,
08/2020, Letnik:
58, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Understanding the properties of dust attenuation curves in galaxies and the physical mechanisms that shape them are among the fundamental questions of extragalactic astrophysics, with great practical ...significance for deriving the physical properties of galaxies. Attenuation curves result from a combination of dust grain properties, dust content, and the spatial arrangement of dust and different populations of stars. In this review, we assess the state of the field, paying particular attention to extinction curves as the building blocks of attenuation laws. We introduce a quantitative framework to characterize extinction and attenuation curves, present a theoretical foundation for interpreting empirical results, overview an array of observational methods, and review observational results at low and high redshifts. Our main conclusions include the following:
Attenuation curves exhibit a wide range of UV-through-optical slopes, from curves with shallow (Milky Way-like) slopes to those exceeding the slope of the Small Magellanic Cloud extinction curve.
The slopes of the curves correlate strongly with the effective optical opacities, in the sense that galaxies with lower dust column density (lower visual attenuation) tend to have steeper slopes, whereas the galaxies with higher dust column density have shallower (grayer) slopes.
Galaxies exhibit a range of 2175-Å UV bump strengths, including no bump, but, on average, are suppressed compared with the average Milky Way extinction curve.
Theoretical studies indicate that both the correlation between the slope and the dust column as well as variations in bump strength may result from geometric and radiative transfer effects.
ABSTRACT
We present predictions for the evolution of the galaxy dust-to-gas ratio (DGR) and dust-to-metal ratio (DTM) from z = 0 → 6, using a model for the production, growth, and destruction of dust ...grains implemented into the simba cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z = 0, but modestly underpredicts the mass function by ∼×3 at z ∼ 1–2. The z = 0 DGR versus metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR–metallicity relationship between z = 0 and 6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion time-scale, stellar mass, and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM versus the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.
Observations of molecular gas in high-z star-forming galaxies
typically rely on emission from CO lines arising from states with rotational quantum numbers J > 1. Converting these observations to an ...estimate of the CO J = 1-0 intensity, and thus inferring H2 gas masses, requires knowledge of the CO excitation ladder or spectral line energy distribution (SLED). The few available multi-J CO observations of galaxies show a very broad range of SLEDs, even at fixed galaxy mass and star formation rate (SFR), making the conversion to J = 1-0 emission and hence molecular gas mass highly uncertain. Here, we combine numerical simulations of disc galaxies and galaxy mergers with molecular line radiative transfer calculations to develop a model for the physical parameters that drive variations in CO SLEDs in galaxies. An essential feature of our model is a fully self-consistent computation of the molecular gas temperature and excitation structure. We find that, while the shape of the SLED is ultimately determined by difficult-to-observe quantities such as the gas density, temperature and optical depth distributions, all of these quantities are well correlated with the galaxy's mean star formation rate surface density (ΣSFR), which is observable. We use this result to develop a model for the CO SLED in terms of ΣSFR, and show that this model quantitatively reproduces the SLEDs of galaxies over a dynamic range of ∼200 in SFR surface density, at redshifts from z = 0 to 6. This model should make it possible to significantly reduce the uncertainty in deducing molecular gas masses from observations of high-J CO emission.
ABSTRACT
We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via ...sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the H2 content, but for $M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $z\lesssim 1.5$. simba first quenches galaxies at $z\gtrsim 3$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$\odot$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.
In this paper, we provide a physical model for the origin of variations in the shapes and bump strengths of dust attenuation laws in galaxies by combining a large suite of cosmological "zoom-in" ...galaxy formation simulations with 3D Monte Carlo dust radiative transfer calculations. We model galaxies over three orders of magnitude in stellar mass, ranging from Milky Way-like systems to massive galaxies at high redshift. Critically, for these calculations, we employ a constant underlying dust extinction law in all cases and examine how the role of geometry and radiative transfer effects impacts the resultant attenuation curves. Our main results follow. Despite our usage of a constant dust extinction curve, we find dramatic variations in the derived attenuation laws. The slopes of normalized attenuation laws depend primarily on the complexities of star-to-dust geometry. Increasing fractions of unobscured young stars flatten normalized curves, while increasing fractions of unobscured old stars steepen curves. Similar to the slopes of our model attenuation laws, we find dramatic variation in the 2175 ultraviolet bump strength, including a subset of curves with little to no bump. These bump strengths are primarily influenced by the fraction of unobscured O and B stars in our model, with the impact of scattered light having only a secondary effect. Taken together, these results lead to a natural relationship between the attenuation curve slope and 2175 bump strength. Finally, we apply these results to a 25 Mpc h−1 box cosmological hydrodynamic simulation in order to model the expected dispersion in attenuation laws at integer redshifts from z = 0 to 6. A significant dispersion is expected at low redshifts and decreases toward z = 6. We provide tabulated results for the best-fit median attenuation curve at all redshifts.
Abstract
The primary method for inferring the stellar mass (
M
*
) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy ...star formation history (SFH) and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here we examine the effect that the assumed SFH has on the stellar properties inferred from SED fitting by ground-truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with underreported uncertainties on results. We demonstrate that flexible nonparametric SFHs outperform traditional parametric forms in capturing variations in galaxy SFHs and, as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-
τ
model to a bias under 0.1 dex for the nonparametric model, though this is heavily dependent on the choice of prior for the nonparametric model. Similarly, using nonparametric SFHs in SED fitting results in increased accuracy in recovered galaxy star formation rates and stellar ages.
ABSTRACT
We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, namely SIMBA, EAGLE, and IllustrisTNG, versus observations from z ...∼ 0 to 2. These simulations all rely on similar subresolution prescriptions to model cold interstellar gas that they cannot represent directly, and qualitatively reproduce the observed z ≈ 0 H i and H2 mass functions (HIMFs and H2MFs, respectively), CO(1–0) luminosity functions (COLFs), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density μ*, with some quantitative differences. To compare to the COLF, we apply an H2-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in αCO and significant differences between models. Using this, predicted z = 0 COLFs agree better with data than predicted H2MFs. Out to z ∼ 2, EAGLE’s and SIMBA’s HIMFs and COLFs strongly increase, while IllustrisTNG’s HIMF declines and COLF evolves slowly. EAGLE and simba reproduce high-LCO(1–0) galaxies at z ∼ 1–2 as observed, owing partly to a median αCO(z = 2) ∼ 1 versus αCO(z = 0) ∼ 3. Examining H i, H2, and CO scaling relations, their trends with M* are broadly reproduced in all models, but EAGLE yields too little H i in green valley galaxies, IllustrisTNG and SIMBA overproduce cold gas in massive galaxies, and SIMBA overproduces molecular gas in small systems. Using SIMBA variants that exclude individual active galactic nucleus (AGN) feedback modules, we find that SIMBA’s AGN jet feedback is primarily responsible by lowering cold gas contents from z ∼ 1 → 0 by suppressing cold gas in $M_*\gtrsim 10^{10}{\rm \,M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-μ* systems.
We report the serendipitous detection of two 3 mm continuum sources found in deep ALMA Band 3 observations to study intermediate-redshift galaxies in the COSMOS field. One is near a foreground galaxy ...at 1 3, but is a previously unknown dust-obscured star-forming galaxy (DSFG) at probable zCO = 3.329, illustrating the risk of misidentifying shorter wavelength counterparts. The optical-to-millimeter spectral energy distribution (SED) favors a gray λ−0.4 attenuation curve and results in significantly larger stellar mass and SFR compared to a Calzetti starburst law, suggesting caution when relating progenitors and descendants based on these quantities. The other source is missing from all previous optical/near-infrared/submillimeter/radio catalogs ("ALMA-only"), and remains undetected even in stacked ultradeep optical (>29.6 AB) and near-infrared (>27.9 AB) images. Using the ALMA position as a prior reveals faint signal-to-noise ratio ∼ 3 measurements in stacked IRAC 3.6+4.5, ultradeep SCUBA2 850 m, and VLA 3 GHz, indicating the source is real. The SED is robustly reproduced by a massive M* = 1010.8M and Mgas = 1011M , highly obscured AV ∼ 4, star-forming SFR ∼ 300 M yr−1 galaxy at redshift z = 5.5 1.1. The ultrasmall 8 arcmin2 survey area implies a large yet uncertain contribution to the cosmic star formation rate density CSFRD(z = 5) ∼ 0.9 × 10−2 M yr−1 Mpc−3, comparable to all ultraviolet-selected galaxies combined. These results indicate the existence of a prominent population of DSFGs at z > 4, below the typical detection limit of bright galaxies found in single-dish submillimeter surveys, but with larger space densities ∼3 × 10−5 Mpc−3, higher duty cycles of 50%-100%, contributing more to the CSFRD, and potentially dominating the high-mass galaxy stellar mass function.
Abstract
We present the Cosmic Sands suite of cosmological zoom-in simulations based on the
simba
galaxy formation model in order to study the buildup of the first massive and dusty galaxies in the ...early universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the buildup of stellar mass. The galaxies are already appreciably chemically evolved by
z
= 7, with modeled dust masses comparable to those inferred from observations in the same epoch, except for the most extreme systems. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250
M
⊙
yr
−1
, but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust RT, we find that, while the infrared luminosities of the most-dust-rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical
z
= 2 submillimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies.