D. Calzetti

GOODS–Herschel: an infrared main sequence for star-forming galaxies

We present the deepest 100 to 500 μm far-infrared observations obtained with the Herschel Space Observatory as part of the GOODS-Herschel key program, and examine the infrared (IR) 3–500 μm spectral energy distributions (SEDs) of galaxies at 0 < z < 2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer, and AKARI data. We determine the projected star formation densities of local galaxies from their radio and mid-IR continuum sizes. We find that the ratio of total IR luminosity to rest-frame 8 μm luminosity, IR8 (≡ L_(IR)^(tot)/L_8), follows a Gaussian distribution centered on IR8 = 4 (σ = 1.6) and defines an IR main sequence for star-forming galaxies independent of redshift and luminosity. Outliers from this main sequence produce a tail skewed toward higher values of IR8. This minority population (  3 × 10^(10) L_⊙ kpc^(-2)) and a high specific star formation rate (i.e., starbursts). The rest-frame, UV-2700 A size of these distant starbursts is typically half that of main sequence galaxies, supporting the correlation between star formation density and starburst activity that is measured for the local sample. Locally, luminous and ultraluminous IR galaxies, (U)LIRGs (L_(IR)^(tot)≥ 10^(11) L_☉), are systematically in the starburst mode, whereas most distant (U)LIRGs form stars in the “normal” main sequence mode. This confusion between two modes of star formation is the cause of the so-called “mid-IR excess” population of galaxies found at z > 1.5 by previous studies. Main sequence galaxies have strong polycyclic aromatic hydrocarbon (PAH) emission line features, a broad far-IR bump resulting from a combination of dust temperatures (T_(dust) ~ 15–50 K), and an effective T_(dust)  ~ 31 K, as derived from the peak wavelength of their infrared SED. Galaxies in the starburst regime instead exhibit weak PAH equivalent widths and a sharper far-IR bump with an effective T_(dust)~ 40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray active galactic nuclei (AGN) is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty starbursts. After correcting for the effect of starbursts on IR8, we identify new candidates for extremely obscured AGNs.

The calibration of mid-infrared star formation rate indicators

With the goal of investigating the degree to which the MIR emission traces the SFR, we analyze Spitzer 8 and 24 μm data of star-forming regions in a sample of 33 nearby galaxies with available HST NICMOS images in the Paα (1.8756 μm) emission line. The galaxies are drawn from the SINGS sample and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and LIRGs are also included in the analysis. Both the stellar continuum-subtracted 8 μm emission and the 24 μm emission correlate with the extinction-corrected Paα line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed nonlinear trend of the 24 μm emission versus number of ionizing photons, including the modest deficiency of 24 μm emission in the low-metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 μm emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 μm emission is contributed, in larger measure than the 24 μm emission, by dust heated by nonionizing stellar populations, in addition to the ionizing ones, in agreement with previous findings. Two SFR calibrations, one using the 24 μm emission and the other using a combination of the 24 μm and Hα luminosities (Kennicutt and coworkers), are presented. No calibration is presented for the 8 μm emission because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by ongoing star formation.

The Spitzer Local Volume Legacy: Survey Description and Infrared Photometry

The survey description and the near-, mid-, and far-infrared flux properties are presented for the 258 galaxies in the Local Volume Legacy (LVL). LVL is a Spitzer Space Telescope legacy program that surveys the local universe out to 11 Mpc, built upon a foundation of ultraviolet, Hα, and Hubble Space Telescope imaging from 11HUGS (11 Mpc Hα and Ultraviolet Galaxy Survey) and ANGST (ACS Nearby Galaxy Survey Treasury). LVL covers an unbiased, representative, and statistically robust sample of nearby star-forming galaxies, exploiting the highest extragalactic spatial resolution achievable with Spitzer. As a result of its approximately volume-limited nature, LVL augments previous Spitzer observations of present-day galaxies with improved sampling of the low-luminosity galaxy population. The collection of LVL galaxies shows a large spread in mid-infrared colors, likely due to the conspicuous deficiency of 8 μm polycyclic aromatic hydrocarbon emission from low-metallicity, low-luminosity galaxies. Conversely, the far-infrared emission tightly tracks the total infrared emission, with a dispersion in their flux ratio of only 0.1 dex. In terms of the relation between the infrared-to-ultraviolet ratio and the ultraviolet spectral slope, the LVL sample shows redder colors and/or lower infrared-to-ultraviolet ratios than starburst galaxies, suggesting that reprocessing by dust is less important in the lower mass systems that dominate the LVL sample. Comparisons with theoretical models suggest that the amplitude of deviations from the relation found for starburst galaxies correlates with the age of the stellar populations that dominate the ultraviolet/optical luminosities.

CALIBRATING EXTINCTION-FREE STAR FORMATION RATE DIAGNOSTICS WITH 33 GHz FREE-FREE EMISSION IN NGC 6946

Using free-free emission measured in the Ka band (26-40 GHz) for 10 star-forming regions in the nearby galaxy NGC 6946, including its starbursting nucleus, we compare a number of star formation rate (SFR) diagnostics that are typically considered to be unaffected by interstellar extinction. These diagnostics include non-thermal radio (i.e., 1.4 GHz), total infrared (IR; 8-1000 μm), and warm dust (i.e., 24 μm) emission, along with hybrid indicators that attempt to account for obscured and unobscured emission from star-forming regions including Hα + 24 μm and UV + IR measurements. The assumption is made that the 33 GHz free-free emission provides the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24 μm, Hα + 24 μm, and UV + IR SFR calibrations are in good agreement with the 33 GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimates the nuclear SFR by a factor of ~2 relative to the 33 GHz free-free SFR. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the 33 GHz free-free SFRs for the nucleus and underestimate the SFRs from the extranuclear star-forming regions by an average factor of ~2 and ~4-5 before and after subtracting local background emission, respectively. This result likely arises from the cosmic-ray (CR) electrons decaying within the starburst region with negligible escape, whereas the transient nature of star formation in the young extranuclear star-forming complexes allows for CR electrons to diffuse significantly further than dust-heating photons, resulting in an underestimate of the true SFR. Finally, we find that the SFRs estimated using the total 33 GHz flux density appear to agree well with those estimated using free-free emission due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33 GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

The calibration of monochromatic far-infrared star formation rate indicators

Spitzer data at 24, 70, and 160 μm and ground-based Hα images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 μm. We compare recently published recipes for SFR measures using combinations of the 24 μm and observed Hα luminosities with those using 24 μm luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 μm and 160 μm luminosity are found for L(70) ≳ 1.4 × 10^(42) erg s^(–1) and L(160) ≳ 2 × 10^(42) erg s^(–1), corresponding to SFR ≳ 0.1-0.3 M_☉ yr^(–1), and calibrations of SFRs based on L(70) and L(160) are proposed. Below those two luminosity limits, the relation between SFR and 70 μm (160 μm) luminosity is nonlinear and SFR calibrations become problematic. A more important limitation is the dispersion of the data around the mean trend, which increases for increasing wavelength. The scatter of the 70 μm (160 μm) data around the mean is about 25% (factor ~2) larger than the scatter of the 24 μm data. We interpret this increasing dispersion as an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. Thus, the 70 (160) μm luminosity can be reliably used to trace SFRs in large galaxy samples, but will be of limited utility for individual objects, with the exception of infrared-dominated galaxies. The nonlinear relation between SFR and the 70 and 160 μm emission at faint galaxy luminosities suggests a variety of mechanisms affecting the infrared emission for decreasing luminosity, such as increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star-forming regions across the galaxy. In all cases, the calibrations hold for galaxies with oxygen abundance higher than roughly 12 +log(O/H) ~ 8.1. At lower metallicity, the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

The CO-to-H2 Conversion Factor and Dust-to-gas Ratio on Kiloparsec Scales in Nearby Galaxies

We present ~kiloparsec spatial resolution maps of the CO-to-H_2 conversion factor (α_(CO)) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for α_(CO) and the DGR by assuming that the DGR is approximately constant on kiloparsec scales. With this assumption, we can combine maps of dust mass surface density, CO-integrated intensity, and H I column density to solve for both αCO and the DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high-resolution far-IR maps from the Herschel key program KINGFISH, ^(12)CO J = (2-1) maps from the IRAM 30 m large program HERACLES, and H I 21 cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our α_(CO) results on the more typically used ^(12)CO J = (1-0) scale and show using literature measurements that variations in the line ratio do not affect our results. In total, we derive 782 individual solutions for α_(CO) and the DGR. On average, α_(CO) = 3.1 M_☉ pc^(–2) (K km s^(–1))^(–1) for our sample with a standard deviation of 0.3 dex. Within galaxies, we observe a generally flat profile of α_(CO) as a function of galactocentric radius. However, most galaxies exhibit a lower α_(CO) value in the central kiloparsec—a factor of ~2 below the galaxy mean, on average. In some cases, the central α_(CO) value can be factors of 5-10 below the standard Milky Way (MW) value of α_(CO,MW) = 4.4 M_☉ pc^(–2) (K km s^(–1))^(–1). While for α_(CO) we find only weak correlations with metallicity, the DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate α_(CO) for studies of nearby galaxies.

KINGFISH—Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multiwavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 μm. The spectral line imaging of the principal atomic ISM cooling lines ([O I] 63 μm, [O III] 88 μm, [N II] 122,205 μm, and [C II] 158 μm) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multiwavelength data sets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This article summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6–500 μm. All imaging and spectroscopy data products will be released to the Herschel user-generated product archives.

RADIAL DISTRIBUTION OF STARS, GAS, AND DUST IN SINGS GALAXIES. II. DERIVED DUST PROPERTIES

We present a detailed analysis of the radial distribution of dust properties in the SINGS sample, performed on a set of ultraviolet (UV), infrared (IR), and Hi surface brightness profiles, combined with published molecular gas profiles and metallicity gradients. The internal extinction, derived from the total-IR (TIR)-to-far-UV (FUV) luminosity ratio, decreases with radius, and is larger in Sb-Sbc galaxies. The TIR-to-FUV ratio correlates with the UV spectral slope β, following a sequence shifted to redder UV colors with respect to that of starbursts. The star formation history (SFH) is identified as the main driver of this departure. Both L_TIR)/L_FUV and β correlate well with metallicity, especially in moderately face-on galaxies. The relation shifts to redder colors with increased scatter in more edge-on objects. By applying physical dust models to our radial spectral energy distributions, we have derived radial profiles of the total dust mass surface density, the fraction of the total dust mass contributed by polycyclic aromatic hydrocarbons (PAHs), and the intensity of the radiation field heating the grains. The dust profiles are exponential, their radial scale length being constant from Sb to Sd galaxies (only ~ 10% larger than the stellar scale length). Many S0/a-Sab galaxies have central depressions in their dust radial distributions. The PAH abundance increases with metallicity for 12 + log(O/H) < 9, and at larger metallicities the trend flattens and even reverses, with the SFH being a plausible underlying driver for this behavior. The dust-to-gas ratio is also well correlated with metallicity and therefore decreases with galactocentric radius. Although most of the total emitted IR power (especially in the outer regions of disks) is contributed by dust grains heated by diffuse starlight with a similar intensity as the local Milky Way radiation field, a small amount of the dust mass (~ 1%) is required to be exposed to very intense starlight in order to reproduce the observed fluxes at 24 μ m, accounting for ~ 10% of the total integrated IR power.

Radial Distribution of Stars, Gas, and Dust in SINGS Galaxies. I. Surface Photometry and Morphology

We present ultraviolet through far-infrared (FIR) surface brightness profiles for the 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). The imagery used to measure the profiles includes Galaxy Evolution Explorer UV data, optical images from Kitt Peak National Observatory, Cerro Tololo Inter-American Observatory, and Sloan Digital Sky Survey, near-IR data from Two Micron All Sky Survey, and mid- and FIR images from Spitzer. Along with the radial profiles, we also provide multi-wavelength asymptotic magnitudes and several nonparametric indicators of galaxy morphology: the concentration index (C42), the asymmetry (A), the Gini coefficient (G), and the normalized second-order moment of the brightest 20% of the galaxy’s flux (M20). In this paper, the first of a series, we describe the technical aspects regarding the surface photometry, and present a basic analysis of the global and structural properties of the SINGS galaxies at different wavelengths. The homogeneity in the acquisition, reduction, and analysis of the results presented here makes these data ideal for multiple unanticipated studies on the radial distribution of the properties of stars, dust, and gas in galaxies. Our radial profiles show a wide range of morphologies and multiple components (bulges, exponential disks, inner and outer disk truncations, etc.) that vary not only from galaxy to galaxy but also with wavelength for a given object. In the optical and near-IR, the SINGS galaxies occupy the same regions in the C42–A-G-M20 parameter space as other normal galaxies in previous studies. However, they appear much less centrally concentrated, more asymmetric, and with larger values of G when viewed in the UV (due to star-forming clumps scattered across the disk) and in the mid-IR (due to the emission of polycyclic aromatic hydrocarbons at 8.0 μm and very hot dust at 24 μm). In an accompanying paper by Mu˜noz-Mateos et al., we focus on the radial distribution of dust properties in the SINGS galaxies, providing a detailed analysis of the radial variation of the attenuation, the dust column density, the dust-to-gas ratio, the abundance of PAHs, and the intensity of the heating starlight.

Warm Molecular Hydrogen in the Spitzer SINGS Galaxy Sample

Results on the properties of warm molecular hydrogen in 57 normal galaxies are derived from measurements of H_2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of emission lines of H_2 to fainter and more common systems (L_(FIR) = 10^(7)-6 × 10^(10) L_☉). The 17 μm S(1) transition is securely detected in the nuclear regions of 86% of galaxies with stellar masses above 10^(9.5) M_☉. The derived column densities of warm H_2 (T ≥ 100 K), although averaged over kiloparsec-scale areas, are commensurate with values observed in resolved photodissociation regions. They amount to between 1% and >30% of the total H_2. The power emitted in the three lowest energy transitions is on average 30% of the power of the bright [Si II] cooling line (34.8 μm) and about 4 × 10^(-4) of the total infrared power for star-forming galaxies, which is consistent with excitation in PDRs. The fact that the H_2 line intensities scale tightly with the aromatic band emission, even though the average radiation field intensity varies by a factor of 10, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. Many of the 25 LINER/Seyfert targets strongly depart from the rest of the sample, in having warmer excited H_2 and excess H_2 rotational power with respect to the dust emission. We propose a threshold in H_2-to-aromatic band power ratios, allowing the identification of low-luminosity AGNs by an excess H_2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect in nearly half the star-forming targets nonequilibrium ortho-to-para ratios, consistent with the effects of FUV pumping combined with incomplete ortho-para thermalization, or possibly nonequilibrium photodissociation fronts.