Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Interstellar dust and starlight are modeled for the galaxies of the project “Key Insights on Nearby Galaxies: A FarInfrared Survey with Herschel.” The galaxies were observed by the Infrared Array Camera and the Multiband Imaging Photometer for Spitzer on Spitzer Space Telescope, and the Photodetector Array Camera and Spectrometer and the Spectral and Photometric Imaging Receiver on Herschel Space Observatory. With data from 3.6 to 500 μm, dust models are strongly constrained. Using a physical dust model, for each pixel in each galaxy we estimate (1) dust surface density, (2) dust mass fraction in polycyclic aromatic hydrocarbons (PAHs), (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in subregions with high starlight intensity. The dust models successfully reproduce the observed global and resolved spectral energy distributions. With the angular resolution of Herschel, we obtain well-resolved maps (available online) for the dust properties. As in previous studies, we find the PAH fraction qPAH to be an increasing function of metallicity, with a threshold oxygen abundance Z/Ze ≈ 0.1, but we find the data to be fitted best with qPAH increasing linearly with log O H ( ) above a threshold value of 0.15(O/H)e. We obtain total dust masses for each galaxy by summing the dust mass over the individual map pixels; these “resolved” dust masses are consistent with the masses inferred from a model fit to the global photometry. The global dust-to-gas ratios obtained from this study are found to correlate with galaxy metallicities. Systems with Z/Ze 0.5 have most of their refractory elements locked up in dust, whereas in systems with Z/Ze 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that qPAH is suppressed in regions with unusually warm dust with nL L n ( ) 70 m 0.4 m dust. With knowledge of one long-wavelength flux density ratio (e.g., f160/f500), the minimum starlight intensity heating the dust (Umin) can be estimated to within ∼50%, despite a variation in Umin of more than two orders of magnitude. For the adopted dust model, dust masses can be estimated to within ∼0.2 dex accuracy using the f160/f500 flux ratio and the integrated dust luminosity, and to ∼0.07 dex accuracy using the 500 μm luminosity nLn ( ) 500 mm alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions for estimating starlight heating intensity and dust mass may be useful for studies of high-redshift galaxies

Quantifying non-star formation associated 8um dust emission in NGC 628

Combining Ha and IRAC images of the nearby spiral galaxy NGC 628, we find that between 30-43% of its 8um dust emission is not related to recent star formation. Contributions from dust heated by young stars are separated by identifying HII regions in the Ha map and using these areas as a mask to determine the 8um dust emission that must be due to heating by older stars. Corrections are made for sub-detection-threshold HII regions, photons escaping from HII regions and for young stars not directly associated to HII regions (i.e. 10-100 Myr old stars). A simple model confirms this amount of 8um emission can be expected given dust and PAH absorption cross-sections, a realistic star-formation history, and the observed optical extinction values. A Fourier power spectrum analysis indicates that the 8um dust emission is more diffuse than the Ha emission (and similar to observed HI), supporting our analysis that much of the 8um-emitting dust is heated by older stars. The 8um dust-to-Ha emission ratio declines with galactocentric radius both within and outside of HII regions, probably due to a radial increase in disk transparency. In the course of this work, we have also found that intrinsic diffuse Ha fractions may be lower than previously thought in galaxies, if the differential extinction between HII regions and diffuse regions is taken into account.Comment: 14 pages, 11 figures, accepted in Ap

Investigating the Presence of 500 μm Submillimeter Excess Emission in Local Star Forming Galaxies

Submillimeter excess emission has been reported at 500 μm in a handful of local galaxies, and previous studies suggest that it could be correlated with metal abundance. We investigate the presence of an excess submillimeter emission at 500 μm for a sample of 20 galaxies from the Key Insights on Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) that span a range of morphologies and metallicities (12 + log (O/H) = 7.8-8.7). We probe the far-infrared (IR) emission using images from the Spitzer Space Telescope and Herschel Space Observatory in the wavelength range 24-500 μm. We model the far-IR peak of the dust emission with a two-temperature modified blackbody and measure excess of the 500 μm photometry relative to that predicted by our model. We compare the submillimeter excess, where present, with global galaxy metallicity and, where available, resolved metallicity measurements. We do not find any correlation between the 500 μm excess and metallicity. A few individual sources do show excess (10%-20%) at 500 μm; conversely, for other sources, the model overpredicts the measured 500 μm flux density by as much as 20%, creating a 500 μm "deficit." None of our sources has an excess larger than the calculated 1σ uncertainty, leading us to conclude that there is no substantial excess at submillimeter wavelengths at or shorter than 500 μm in our sample. Our results differ from previous studies detecting 500 μm excess in KINGFISH galaxies largely due to new, improved photometry used in this study

The Spatial Distribution of Dust and Stellar Emission of the Magellanic Clouds

We study the emission by dust and stars in the Large and Small Magellanic Clouds, a pair of low-metallicity nearby galaxies, as traced by their spatially resolved spectral energy distributions (SEDs). This project combines Herschel Space Observatory PACS and SPIRE far-infrared photometry with other data at infrared and optical wavelengths. We build maps of dust and stellar luminosity and mass of both Magellanic Clouds, and analyze the spatial distribution of dust/stellar luminosity and mass ratios. These ratios vary considerably throughout the galaxies, generally between the range $0.01\leq L_{\rm dust}/L_\ast\leq 0.6$ and $10^{-4}\leq M_{\rm dust}/M_\ast\leq 4\times10^{-3}$. We observe that the dust/stellar ratios depend on the interstellar medium (ISM) environment, such as the distance from currently or previously star-forming regions, and on the intensity of the interstellar radiation field (ISRF). In addition, we construct star formation rate (SFR) maps, and find that the SFR is correlated with the dust/stellar luminosity and dust temperature in both galaxies, demonstrating the relation between star formation, dust emission and heating, though these correlations exhibit substantial scatter.Comment: 15 pages, 18 figures; ApJ, in press; version published in the journal will have higher-resolution figure

Mapping dust through emission and absorption in nearby galaxies

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and similar to 1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest H II regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A(V) measures to infer global dust properties

Modeling Dust in the Interstellar Medium

We are in a very special moment for the study of the interstellar medium (ISM). The Spitzer Space Telescope had provided, and currently Herschel Space Observatory is providing, invaluable infrared (IR) observations of a variety of astrophysical systems. These observations allow us to model several ongoing processes in the ISM, and in particular to study the physical properties of the interstellar dust. Determining the dust properties accurately is an extremely difficult task: even the overall amount of dust in other galaxies has often been very uncertain. In the current work, we develop ``state of the art'' tools for image processing and dust modeling that allows study of the interstellar dust in other galaxies using the new infrared data. We start by developing, the now ``industry-standard'', convolution kernels. They allow us to accurately combine data from several space- and ground-based telescopes, to perform multi-wavelength studies. They are a key development for doing resolved studies of astrophysical systems. We follow by analyzing the performance of ``modified blackbody'' (MBB) dust models when applied to realistic spectral energy distributions (SEDs), where we use a specific physical model, the Draine and Li (2007, DL07) dust model, to generate the synthetic SEDs. We show that MBB models can have a large bias in the inferred dust parameters, and therefore it is important to use more realistic dust models. We provide ``correction'' formulae to compensate for the MBB bias, useful when the more sophisticated dust modeling is not available. Using the DL07 dust model, which contains amorphous silicate and carbonaceous grains, we perform careful modeling of the dust properties in a large sample of well-resolved galaxies observed by the KINGFISH survey. With data from 3.6um to 500um, dust models are strongly constrained. For each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). We find no evidence for significant masses of cold dust (T<12K). For two galaxies studied in detail (NGC628 and NGC6946) the derived dust maps correlates extremely well with independent observations of emission in the HI 21cm line and CO1-0 line. The derived dust/gas mass ratio are in excellent agreement with dust/gas ratios infered from other lines of evidence

Cloud-scale ISM Structure and Star Formation in M51

International audienc