E. J. Murphy

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 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.

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.

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.

A dust-obscured massive maximum-starburst galaxy at a redshift of 6.34

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts—that is, increased rates of star formation—in the most massive dark-matter haloes at early epochs. However, it remains unknown how soon after the Big Bang massive starburst progenitors exist. The measured redshift (z) distribution of dusty, massive starbursts has long been suspected to be biased low in z owing to selection effects, as confirmed by recent findings of systems with redshifts as high as ∼5 (refs 2–4). Here we report the identification of a massive starburst galaxy at z = 6.34 through a submillimetre colour-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine-structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40 per cent of the baryonic mass. A ‘maximum starburst’ converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn in cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

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 detection of a population of submillimeter-bright, strongly lensed galaxies

Through a Lens Brightly Astronomical sources detected in the submillimeter range are generally thought to be distant, dusty galaxies undergoing a vigorous burst of star formation. They can be detected because the dust absorbs the light from stars and reemits it at longer wavelengths. Their properties are still difficult to ascertain, however, because the combination of interference from dust and the low spatial resolution of submillimeter telescopes prevents further study at other wavelengths. Using data from the Herschel Space Telescope, Negrello et al. (p. 800) showed that by searching for the brightest sources in a wide enough area in the sky it was possible to detect gravitationally lensed submillimeter galaxies with nearly full efficiency. Gravitational lensing occurs when the light of an astronomical object is deflected by a foreground mass. This phenomenon increases the apparent brightness and angular size of the lensed objects, making it easier to study sources that would be otherwise too faint to probe. Data from the Herschel Space Observatory unveils distant, dusty galaxies invisible to optical telescopes. Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty star-forming galaxies. However, the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.