Daniel A. Dale

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multiwavelength library of ancillary and complementary observations, including radio continuum, H i, CO, submillimeter, BVRIJHK ,H a ,P aa, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.

The Infrared Spectral Energy Distribution of Normal Star-forming Galaxies: Calibration at Far-Infrared and Submillimeter Wavelengths

New far-infrared and submillimeter data are used to solidify and to extend to long wavelengths the empirical calibration of the infrared spectral energy distribution (SED) of normal star-forming galaxies. As was found by Dale and coworkers in 2001, a single parameter family, characterized by fν(60 μm)/fν(100 μm), is adequate to describe the range of normal galaxy SEDs observed by the Infrared Astronomical Satellite and Infrared Space Observatory from 3 to 100 μm. However, predictions based on the first-generation models at longer wavelengths (122-850 μm) are increasingly overluminous compared to the data for smaller fν(60 μm)/fν(100 μm), or alternatively, for weaker global interstellar radiation fields. After slightly modifying the far-infrared/submillimeter dust emissivity in those models as a function of the radiation field intensity to better match the long-wavelength data, a suite of SEDs from 3 μm to 20 cm in wavelength is presented. Results from relevant applications are also discussed, including submillimeter-based photometric redshift indicators, the infrared energy budget and simple formulae for recovering the bolometric infrared luminosity, and dust mass estimates in galaxies. Regarding the latter, since galaxy infrared SEDs are not well described by single blackbody curves, the usual methods of estimating dust masses can be grossly inadequate. The improved model presented herein is used to provide a more accurate relation between infrared luminosity and dust mass.

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 MID-INFRARED SPECTRUM OF STAR-FORMING GALAXIES: GLOBAL PROPERTIES OF POLYCYCLIC AROMATIC HYDROCARBON EMISSION

We present a sample of low-resolution 5-38 μm Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths that indicate that most normal galaxies undergo A_V ≲ 3 mag averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20% and is suppressed at metallicities Z ≲ Z_☉/4, as well as in low-luminosity AGN environments. Strong interband PAH feature strength variations (2-5 times) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8 μm features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGNs. The AGNs may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGNs may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled and points to more than a factor of 2 uncertainty in results that assume a fixed PAH emission spectrum, for redshifts z = 0-2.5.

Star Formation in NGC 5194 (M51a). II. The Spatially Resolved Star Formation Law

We have studied the relationship between the star formation rate (SFR), surface density, and gas surface density in the spiral galaxy M51a (NGC 5194), using multiwavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). We introduce a new SFR index based on a linear combination of Hα emission-line and 24 μm continuum luminosities, which provides reliable extinction-corrected ionizing fluxes and SFR densities over a wide range of dust attenuations. The combination of these extinction-corrected SFR densities with aperture synthesis H I and CO maps has allowed us to probe the form of the spatially resolved star formation law on scales of 0.5-2 kpc. We find that the resolved SFR versus gas surface density relation is well represented by a Schmidt power law, which is similar in form and dispersion to the disk-averaged Schmidt law. We observe a comparably strong correlation of the SFR surface density with the molecular gas surface density, but no significant correlation with the surface density of atomic gas. The best-fitting slope of the Schmidt law varies from N = 1.37 to 1.56, with zero point and slope that change systematically with the spatial sampling scale. We tentatively attribute these variations to the effects of areal sampling and averaging of a nonlinear intrinsic star formation law. Our data can also be fitted by an alternative parameterization of the SFR surface density in terms of the ratio of gas surface density to local dynamical time, but with a considerable dispersion.

The Infrared spectral energy distribution of normal star-forming galaxies

We have developed a new phenomenological model for the spectral energy distribution of normal star-forming galaxies between 3 and 1100 μm. These model spectra allow us to determine the infrared energy budget for normal galaxies, and in particular to translate far-infrared fluxes into total (bolometric) infrared fluxes. The 20 to 42 m range appears to show the most significant growth in relative terms as the activity level increases, suggesting that the 20–42 m continuum may be the best dust emission tracer of current star formation in galaxies.

Star Formation in NGC 5194 (M51a): The Panchromatic View from GALEX to Spitzer*

(Abridged) Far ultraviolet to far infrared images of the nearby galaxy NGC5194, from Spitzer, GALEX, Hubble Space Telescope and ground--based data, are used to investigate local and global star formation, and the impact of dust extinction in HII-emitting knots. In the IR/UV-UV color plane, the NGC5194 HII knots show the same trend observed for normal star-forming galaxies, having a much larger dispersion than starburst galaxies. We identify the dispersion as due to the UV emission predominantly tracing the evolved, non-ionizing stellar population, up to ages 50-100 Myr. While in starbursts the UV light traces the current SFR, in NGC5194 it traces a combination of current and recent-past SFR. Unlike the UV emission, the monochromatic 24 micron luminosity is an accurate local SFR tracer for the HII knots in NGC5194; this suggests that the 24 micron emission carriers are mainly heated by the young, ionizing stars. However, preliminary results show that the ratio of the 24 micron emission to the SFR varies by a factor of a few from galaxy to galaxy. While also correlated with star formation, the 8 micron emission is not directly proportional to the number of ionizing photons. This confirms earlier suggestions that the carriers of the 8 micron emission are heated by more than one mechanism.

COMPARISON OF Hα AND UV STAR FORMATION RATES IN THE LOCAL VOLUME: SYSTEMATIC DISCREPANCIES FOR DWARF GALAXIES

Using a complete sample of ~300 star-forming galaxies within 11 Mpc of the Milky Way, we evaluate the consistency between star formation rates (SFRs) inferred from the far ultraviolet (FUV) non-ionizing continuum and Hα nebular emission, assuming standard conversion recipes in which the SFR scales linearly with luminosity at a given wavelength. Our analysis probes SFRs over 5 orders of magnitude, down to ultra-low activities on the order of ~10^–4 M_☉ yr^–1. The data are drawn from the 11 Mpc Hα and Ultraviolet Galaxy Survey (11HUGS), which has obtained Hα fluxes from ground-based narrowband imaging, and UV fluxes from imaging with GALEX. For normal spiral galaxies (SFR ~ 1 M_☉ yr^–1), our results are consistent with previous work which has shown that FUV SFRs tend to be lower than Hα SFRs before accounting for internal dust attenuation, but that there is relative consistency between the two tracers after proper corrections are applied. However, a puzzle is encountered at the faint end of the luminosity function. As lower luminosity dwarf galaxies, roughly less active than the Small Magellanic Cloud, are examined, Hα tends to increasingly underpredict the total SFR relative to the FUV. The trend is evident prior to corrections for dust attenuation, which affects the FUV more than the nebular Hα emission, so this general conclusion is robust to the effects of dust. Although past studies have suggested similar trends, this is the first time this effect is probed with a statistical sample for galaxies with SFR ≤0.1 M_☉ yr^–1. By SFR ~ 0.003 M_☉ yr–1, the average Hα-to-FUV flux ratio is lower than expected by a factor of two, and at the lowest SFRs probed, the ratio exhibits an order of magnitude discrepancy for the handful of galaxies that remain in the sample. A range of standard explanations does not appear to be able to fully account for the magnitude of the systematic. Some recent work has argued for a stellar initial mass function which is deficient in high-mass stars in dwarf and low surface brightness galaxies, and we also consider this scenario. Under the assumption that the FUV traces the SFR in dwarf galaxies more robustly, the prescription relating Hα luminosity to SFR is re-calibrated for use in the low SFR regime when FUV data are not available.

An Ultraviolet-to-Radio Broadband Spectral Atlas of Nearby Galaxies

The ultraviolet-to-radio continuum spectral energy distributions are presented for all 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). A principal component analysis of the sample shows that most of the samples spectral variations stem from two underlying components, one representative of a galaxy with a low infrared-to-ultraviolet ratio and one representative of a galaxy with a high infrared-to-ultraviolet ratio. The influence of several parameters on the infrared-to-ultraviolet ratio is studied (e.g., optical morphology, disk inclination, far-infrared color, ultraviolet spectral slope, and star formation history). Consistent with our understanding of normal star-forming galaxies, the SINGS sample of galaxies in comparison to more actively star-forming galaxies exhibits a larger dispersion in the infrared-to-ultraviolet versus ultraviolet spectral slope correlation. Early-type galaxies, exhibiting low star formation rates and high optical surface brightnesses, have the most discrepant infrared-to-ultraviolet correlation. These results suggest that the star formation history may be the dominant regulator of the broadband spectral variations between galaxies. Finally, a new discovery shows that the 24 μm morphology can be a useful tool for parameterizing the global dust temperature and ultraviolet extinction in nearby galaxies. The dust emission in dwarf/irregular galaxies is clumpy and warm accompanied by low ultraviolet extinction, while in spiral galaxies there is typically a much larger diffuse component of cooler dust and average ultraviolet extinction. For galaxies with nuclear 24 μm emission, the dust temperature and ultraviolet extinction are relatively high compared to disk galaxies.

OPTICAL SPECTROSCOPY AND NEBULAR OXYGEN ABUNDANCES OF THE SPITZER/SINGS GALAXIES

We present intermediate-resolution optical spectrophotometry of 65 galaxies obtained in support of the Spitzer Infrared Nearby Galaxies Survey (SINGS). For each galaxy we obtain a nuclear, circumnuclear, and semi-integrated optical spectrum designed to coincide spatially with mid- and far-infrared spectroscopy from the Spitzer Space Telescope. We make the reduced, spectrophotometrically calibrated one-dimensional spectra, as well as measurements of the fluxes and equivalent widths of the strong nebular emission lines, publically available. We use optical emission-line ratios measured on all three spatial scales to classify the sample into star-forming, active galactic nuclei (AGNs), and galaxies with a mixture of star formation and nuclear activity. We find that the relative fraction of the sample classified as star forming versus AGN is a strong function of the integrated light enclosed by the spectroscopic aperture. We supplement our observations with a large database of nebular emission-line measurements of individual H II regions in the SINGS galaxies culled from the literature. We use these ancillary data to conduct a detailed analysis of the radial abundance gradients and average H II-region abundances of a large fraction of the sample. We combine these results with our new integrated spectra to estimate the central and characteristic (globally averaged) gas-phase oxygen abundances of all 75 SINGS galaxies. We conclude with an in-depth discussion of the absolute uncertainty in the nebular oxygen abundance scale.