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.
We study the dust attenuation curves of 230,000 individual galaxies in the local universe, ranging from quiescent to intensely star-forming systems, using GALEX, SDSS, and WISE photometry calibrated ...on the Herschel ATLAS. We use a new method of constraining SED fits with infrared luminosity (SED+LIR fitting), and parameterized attenuation curves determined with the CIGALE SED-fitting code. Attenuation curve slopes and UV bump strengths are reasonably well constrained independently from one another. We find that attenuation curves exhibit a very wide range of slopes that are on average as steep as the curve slope of the Small Magellanic Cloud (SMC). The slope is a strong function of optical opacity. Opaque galaxies have shallower curves-in agreement with recent radiative transfer models. The dependence of slopes on the opacity produces an apparent dependence on stellar mass: more massive galaxies have shallower slopes. Attenuation curves exhibit a wide range of UV bump amplitudes, from none to Milky Way (MW)-like, with an average strength one-third that of the MW bump. Notably, local analogs of high-redshift galaxies have an average curve that is somewhat steeper than the SMC curve, with a modest UV bump that can be, to first order, ignored, as its effect on the near-UV magnitude is 0.1 mag. Neither the slopes nor the strengths of the UV bump depend on gas-phase metallicity. Functional forms for attenuation laws are presented for normal star-forming galaxies, high-z analogs, and quiescent galaxies. We release the catalog of associated star formation rates and stellar masses (GALEX-SDSS-WISE Legacy Catalog 2).
Abstract
Robust estimation of star formation rates (SFRs) at higher redshifts (
z
≳ 1) using UV–optical–near-infrared (NIR) photometry is contingent on the ability of spectral energy distribution ...(SED) fitting to constrain the dust attenuation, stellar metallicity, and star formation history (SFH) simultaneously. IR-derived dust luminosities can help break the degeneracy between these parameters, but IR data are often not available. Here, we explore strategies for SED fitting at
z
≳ 1 in the absence of IR data using a sample of log
M
*
> 10.2 star-forming galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) for which 24
μ
m data are available. We adopt the total IR luminosity (
L
TIR
) obtained from 24
μ
m as the “ground truth,” which allows us to assess how well it can be recovered (as
L
dust
) from UV–optical–NIR SED fitting. We test a variety of dust attenuation models, stellar population synthesis models, metallicity assumptions, and SFHs separately to identify which assumptions maximize the agreement (correlation and linearity) between
L
TIR
and
L
dust
. We find that a flexible dust attenuation law performs best. For stellar populations, we find that Bruzual & Charlot models are favored over those of Eldridge et al. Fixing the stellar metallicity at solar value is preferred to other fixed values or leaving it as a free parameter. For SFHs, we find that minimizing the variability in the recent (<100 Myr) SFH improves the agreement with
L
TIR
. Finally, we provide a catalog of galaxy parameters (including
M
*
and SFR) for CANDELS galaxies with
log
M
*
>
8
and 0.7 <
z
< 1.3, obtained using the models we found to be the most robust.
Abstract
The Sloan Digital Sky Survey (SDSS) was foundational to the study of galaxy evolution, having revealed the bimodality of galaxies and the relationship between their structure and ...star-forming activity. However, ground-based optical surveys like SDSS are limited in resolution and depth, which may lead to biases or poor quality in the derived morphological properties, potentially impacting our understanding of how and why galaxies cease their star formation (quench). We use archival Hubble Space Telescope (HST) imaging of ∼2000 SDSS objects to assess the reliability of SDSS-derived morphologies, taking advantage of both SDSS statistical samples and HST’s superior resolution and sensitivity. Single Sérsic fitting and bulge-disk decomposition is performed on HST images for direct comparison with SDSS results. Of the three catalogs of SDSS-derived morphologies considered, none is significantly more accurate or precise than the others. For disk-dominated galaxies (
n
< 2.5), global Sérsic indices (
n
) from Meert et al. (hereafter, M15) are preferred. For bulge-dominated galaxies (
n
> 2.5), Simard et al. (hereafter, S11) and M15 overestimate
n
by ∼20%, and derived global
n
from Blanton et al. are preferred. We provide an empirical relation to correct S11 global Sérsic indices. Global
R
eff
from S11 overestimates
R
eff
for the largest galaxies by 0.1 dex. Despite these moderate biases, SDSS-derived single-component parameters are generally significantly more robust than SDSS-derived two-component parameters. The bulge Sérsic index (
n
bulge
) cannot be reliably constrained from SDSS imaging at all, whereas the bulge-to-total (
B
/
T
) ratio can be inferred from SDSS (provided that
n
bulge
= 4 is enforced) but has a large random error of ∼0.2.
In this work, we investigate the reliability of the BPT diagram for excluding galaxies that host an active galactic nucleus (AGN). We determine the prevalence of X-ray AGNs in the star-forming region ...of the BPT diagram and discuss the reasons behind this apparent misclassification, focusing primarily on relatively massive ( ) galaxies. X-ray AGNs are selected from deep XMM observations using a new method that results in greater samples with a wider range of X-ray luminosities, complete to for z < 0.3. Taking X-ray detectability into account, we find that the average fraction of X-ray AGNs in the BPT star-forming branch is 2%, suggesting the BPT diagram can provide a reasonably clean sample of star-forming galaxies. However, the X-ray selection is itself rather incomplete. At the tip of the AGN branch of the BPT diagram, the X-ray AGN fraction is only 14%, which may have implications for studies that exclude AGNs based only on X-ray observations. Interestingly, the X-ray AGN fractions are similar for Seyfert and LINER populations, consistent with LINERs being true AGNs. We find that neither the star formation dilution nor the hidden broad-line components can satisfactorily explain the apparent misclassification of X-ray AGNs. On the other hand, ∼40% of all X-ray AGNs have weak emission lines such that they cannot be placed on the BPT diagram at all and often have low specific SFRs. Therefore, the most likely explanation for "misclassified" X-ray AGNs is that they have intrinsically weak AGN lines, and are only placeable on the BPT diagram when they tend to have high specific SFRs.
Abstract
The Milky Way dust extinction curve in the near-infrared (NIR) follows a power-law form, but the value of the slope,
β
NIR
, is debated. Systematic variations in the slope of the Milky Way ...UV extinction curve are known to be correlated with variations in the optical slope (through
R
V
), but whether such a dependence extends to the NIR is unclear. Finally, because of low dust column densities, the NIR extinction law is poorly understood at high Galactic latitudes where most extragalactic work takes place. In this paper, we construct extinction curves from 56,649 stars with Sloan Digital Sky Survey (SDSS) and Two Micron All Sky Survey photometry, based on stellar parameters from SDSS spectra. We use dust maps to identify dust-free stars, from which we calibrate the relation between stellar parameters and intrinsic colors. Furthermore, to probe the low-dust regime at high latitudes, we use aggregate curves based on many stars. We find no systematic variation of
β
NIR
across low-to-moderate dust columns (0.02 <
E
(
B
−
V
) ≲ 1), and report average
β
NIR
= 1.85 ± 0.01, in agreement with the law in the 2019 Fitzpatrick et al. study, but steeper than the Cardelli et al. and 1999 Fitzpatrick laws. Star-to-star scatter in
β
NIR
is relatively small (σ(
β
NIR
) = 0.13). We also find no intrinsic correlation between
β
NIR
and
R
V
(there is an apparent correlation that is the result of the correlated uncertainties in the two values). These results hold for typical sightlines; we do not probe very dusty regions near the Galactic Center, nor rare sightlines with
R
V
> 4. Finally, we find
R
H
= 0.345 ± 0.007 and comment on its bearing on Cepheid calibrations and the determination of
H
0
.
ABSTRACT In this paper, we present the GALEX-SDSS-WISE Legacy Catalog (GSWLC), a catalog of physical properties (stellar masses, dust attenuations, and star formation rates SFRs) for ∼700,000 ...galaxies with Sloan Digital Sky Survey (SDSS) redshifts below 0.3. GSWLC contains galaxies within the Galaxy Evolution Explorer footprint, regardless of a UV detection, covering 90% of SDSS. The physical properties were obtained from UV/optical spectral energy distribution (SED) fitting following Bayesian methodology of Salim et al., with improvements such as blending corrections for low-resolution UV photometry, flexible dust attenuation laws, and emission-line corrections. GSWLC also includes mid-IR SFRs derived from IR templates based on 22 Wide-field Infrared Survey Explorer observations. These estimates are independent of UV/optical SED fitting, in order to separate possible systematics. The paper argues that the comparison of specific SFRs (sSFRs) is more informative and physically motivated than the comparison of SFRs. The sSFRs resulting from the UV/optical SED fitting are compared to the mid-IR sSFRs and to sSFRs from three published catalogs. For "main-sequence" galaxies with no active galactic nucleus (AGN) all sSFRs are in very good agreement (within 0.1 dex on average). In particular, the widely used aperture-corrected SFRs from the MPA/JHU catalog show no systematic offsets, in contrast to some integral field spectroscopy results. For galaxies below the main sequence (log sSFR ), mid-IR (s)SFRs based on fixed luminosity-SFR conversion are severely biased (up to 2 dex) because the dust is primarily heated by old stars. Furthermore, mid-IR (s)SFRs are overestimated by up to 0.6 dex for galaxies with AGNs, presumably due to nonstellar dust heating. UV/optical (s)SFRs are thus preferred to IR-based (s)SFRs for quenched galaxies and those that host AGNs.
We present an atlas of ultraviolet and infrared images of ∼15,750 local (d 50 Mpc) galaxies, as observed by NASA's Wide-field Infrared Survey Explorer (WISE) and Galaxy Evolution Explorer (GALEX) ...missions. These maps have matched resolution (FWHM 7 5 and 15″), matched astrometry, and a common procedure for background removal. We demonstrate that they agree well with resolved intensity measurements and integrated photometry from previous surveys. This atlas represents the first part of a program (the z = 0 Multiwavelength Galaxy Synthesis) to create a large, uniform database of resolved measurements of gas and dust in nearby galaxies. The images and associated catalogs will be publicly available at the NASA/IPAC Infrared Science Archive. This atlas allows us estimate local and integrated star formation rates (SFRs) and stellar masses (M ) across the local galaxy population in a uniform way. In the appendix, we use the population synthesis fits of Salim et al. to calibrate integrated M and SFR estimators based on GALEX and WISE. Because they leverage a Sloan Digital Sky Survey (SDSS)-based training set of >100,000 galaxies, these calibrations have high precision and allow us to rigorously compare local galaxies to SDSS results. We provide these SFR and M estimates for all galaxies in our sample and show that our results yield a "main sequence" of star-forming galaxies comparable to previous work. We also show the distribution of intensities from resolved galaxies in NUV-to-WISE1 versus WISE1-to-WISE3 space, which captures much of the key physics accessed by these bands.
It has been proposed that the (stellar) mass-(gas) metallicity relation of galaxies exhibits a secondary dependence on star formation rate (SFR), and that the resulting Mlow *-Z-SFR relation may be ...redshift-invariant, i.e., "fundamental." However, conflicting results on the character of the SFR dependence, and whether it exists, have been reported. To gain insight into the origins of the conflicting results, we (1) devise a non-parametric, astrophysically motivated analysis framework based on the offset from the star-forming ("main") sequence at a given Mlow * (relative specific SFR); (2) apply this methodology and perform a comprehensive re-analysis of the local Mlow *-Z-SFR relation, based on SDSS, GALEX, and WISE data; and (3) study the impact of sample selection and of using different metallicity and SFR indicators. We show that metallicity is anti-correlated with specific SFR regardless of the indicators used. We do not find that the relation is spurious due to correlations arising from biased metallicity measurements or fiber aperture effects. We emphasize that the dependence is weak/absent for massive galaxies (log Mlow * > 10.5), and that the overall scatter in the Mlow *-Z-SFR relation does not greatly decrease from the Mlow *-Z relation. We find that the dependence is stronger for the highest SSFR galaxies above the star-forming sequence. This two-mode behavior can be described with a broken linear fit in 12+log(O/H) versus log (SFR/Mlow *), at a given Mlow *. Previous parameterizations used for comparative analysis with higher redshift samples that do not account for the more detailed behavior of the local Mlow *-Z-SFR relation may incorrectly lead to the conclusion that those samples follow a different relationship.
Abstract
The Galaxy Evolution Explorer (GALEX) satellite performed the first and only large-area UV survey, which in tandem with the Sloan Digital Sky Survey (SDSS) has facilitated modeling of the ...spectral energy distributions of low-redshift galaxies and the determination of various galaxy properties, in particular the star formation rate. However, the relatively crude angular resolution of GALEX (5″) made its images susceptible to blending of sources, resulting in potentially biased far-UV and near-UV (NUV) pipeline photometry. To remedy this issue and take advantage of model-fit photometry, we use the EMphot software to obtain forced GALEX photometry for ∼700,000 SDSS galaxies at
z
< 0.3. Positional priors of target galaxies and potentially contaminating neighbors were taken from SDSS. New photometry is based on the best-fitting of three model profiles: optical-like, exponential, and flat. New photometry mitigates blending present in the original pipeline catalogs, which affected 16% of galaxies at a level of >0.2 mag and 2% at a level of >1 mag. Pipeline NUV magnitudes are severely affected (≳1 mag) when the neighbor is brighter than the target galaxy and within 10″, or when the neighbor is fainter and within ∼3″ of the target. New photometry fixes edge-of-detector bias, which affected pipeline photometry by up to 0.1 mag in NUV. We present catalogs with new photometry for GALEX observations of different depths, corresponding to the all-sky imaging survey (AIS), medium imaging survey, and deep imaging survey. Catalogs feature combined magnitudes for multiple detections of the same galaxy in a survey.