When combined with infrared observations with the Spitzer telescope (3 to 160 μm), the Herschel
Space Observatory now fully samples the thermal dust emission up to 500 μm and enables us to better ...estimate the total infrared-submm energy budget (L
TIR) of nearby galaxies. We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved spectral energy distribution (SED) modelling of these objects using the Draine & Li dust models and investigate the influence of the addition of Spectral and Photometric Imaging Receiver (SPIRE) measurements in the estimation of L
TIR. We find that using data up to 250 μm leads to local L
TIR values consistent with those obtained with a complete coverage (up to 500 μm) within ±10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the 70-160 μm band. The 24-70 μm fraction decreases with increasing metallicity. The 160-1100 μmsubmillimetre band can account for up to 25 per cent of the L
TIR in metal-rich galaxies. We investigate the correlation between the total infrared (TIR) surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three Photodetector Array Camera and Spectrometer (PACS) bands can be used as reliable monochromatic estimators of the L
TIR, the 100 μm band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250 μm and L
TIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated L
TIR of nearby galaxies as well as L
TIR of z ∼ 1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.
The C ii 158 m fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, C ii balances the heating, ...including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of C ii emission remains unclear because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of C ii emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the N ii 205 m fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of C ii/N ii 122 m. Using the FIR C ii and N ii emission detected by the KINGFISH (Key Insights on Nearby Galaxies: a Far- Infrared Survey with Herschel) and Beyond the Peak Herschel programs, we show that 60%-80% of C ii emission originates from neutral gas. We find that the fraction of C ii originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and has a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total C ii emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.
The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ...interstellar medium (ISM) in many post-starburst galaxies raise important questions. Are they truly quiescent, and if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infrared spectroscopic and photometric survey of 33 E+A post-starbursts selected by the Sloan Digital Sky Survey, aimed at resolving these questions. We find compact, warm dust reservoirs with high PAH abundances and total gas and dust masses significantly higher than expected from stellar recycling alone. Both polycyclic aromatic hydrocarbon (PAH)/total infrared (TIR) and dust-to-burst stellar mass ratios are seen to decrease with post-burst age, indicative of the accumulating effects of dust destruction and an incipient transition to hot, early-type ISM properties. Their infrared spectral properties are unique, with dominant PAH emission, very weak nebular lines, unusually strong H2 rotational emission, and deep C ii deficits. There is substantial scatter among star formation rate (SFR) indicators, and both PAH and TIR luminosities provide overestimates. Even as potential upper limits, all tracers show that the SFR has typically experienced a decline of more than two orders of magnitude since the starburst and that the SFR is considerably lower than expected given both their stellar masses and molecular gas densities. These results paint a coherent picture of systems in which star formation was, indeed, rapidly truncated, but in which the ISM was not completely expelled, and is instead supported against collapse by latent or continued injection of turbulent or mechanical heating. The resulting aging burst populations provide a "high-soft" radiation field that seemingly dominates the E+A galaxies' unusual ISM energetics.
We have fit the far-ultraviolet (FUV) to sub-millimeter (850 μm) spectral energy distributions (SEDs) of the 61 galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel ...(KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multiwavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple “recipes” for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime: in the 5–10 μm region dominated by PAH emission, and also between 25 and 70 μm where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard recipes for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a “standard” constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 μm) luminosities and colors. From a principal component analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within ∼0.3 dex using only Mstar and SFR.
The CII 157.74 mu m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby ...galaxies, we investigate the relation of CII surface brightness and luminosity with SFR. We conclude that CII can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGNs). The uncertainty of the capital sigma sub(CII) - capital sigma sub(SFR) calibration is + or -0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted capital sigma sub(CII) - capital sigma sub(SFR) correlation is valid over almost five orders of magnitude in capital sigma sub(SFR), holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using CII luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to CII finds that heating efficiencies are l%-3% in normal galaxies.
Abstract
Taking advantage of the unprecedented combination of sensitivity and angular resolution afforded by the Herschel Space Observatory at far-infrared and submillimetre wavelengths, we aim to ...characterize the physical properties of cold dust within nearby galaxies, as well as the associated uncertainties, namely the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) key programme, we perform two-temperature fits of the Spitzer and Herschel photometric data (from 24 to 500 μm), with a warm and a cold component, both globally and in each resolution element. We compare the results obtained from different analysis strategies. At global scale, we observe a range of values of the modified blackbody fit parameters βc (0.8-2.5) and T
c (19.1-25.1 K). We compute maps of our modelling parameters with βc fixed or treated as a free parameter to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe steeper temperature gradients as a function of radius than when it is allowed to vary. When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities. Gathering the parameters obtained in each resolution element in a T
c-βc diagram underlines an anticorrelation between the two parameters. It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar-ends (coinciding with molecular gas density enhancements usually found in these locations). We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (β < 2.0) lead to significantly lower dust masses compared to the β < 2.0 case, with dust masses lower by up to 50 per cent if βc = 1.5, for instance. The working resolution affects our total dust mass estimates: masses increase from global fits to spatially resolved fits.
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 Delta *mm), and warm dust (i.e., 24 Delta *mm) emission, along with hybrid indicators that attempt to account for obscured and unobscured emission from star-forming regions including H Delta *a + 24 Delta *mm 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 Delta *mm, H Delta *a + 24 Delta *mm, 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.
ABSTRACT The H i and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass of nearby galaxies. Both components become too faint to be used as a tracer in ...observations of high-redshift galaxies. In those cases, the 158 m line of atomic carbon (C ii) may be the only way to derive . As the distribution and kinematics of the ISM tracer affects the determination of , it is important to quantify the relative distributions of H i, CO, and C ii. H i and CO are well-characterized observationally, however, for C ii only very few measurements exist. Here we compare observations of CO, H i, and C ii emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS, and KINGFISH surveys. We find that within R25, the average C ii exponential radial profile is slightly shallower than that of the CO, but much steeper than the H i distribution. This is also reflected in the integrated spectrum ("global profile"), where the C ii spectrum looks more like that of the CO than that of the H i. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the C ii data. Using high-spectral-resolution SOFIA C ii data of a number of star forming regions in two nearby galaxies, we find that their C ii linewidths agree better with those of the CO than the H i. As the radial extent of a given ISM tracer is a key input in deriving from spatially unresolved data, we conclude that the relevant length-scale to use in determining based on C ii data, is that of the well-characterized CO distribution. This length scale is similar to that of the optical disk.
New far-infrared and submillimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel ...(KINGFISH) sample. The spatially integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500 mu m emission shows evidence for a submillimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photodissociation regions is found to be (21 + or - 4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.
Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over ...a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields (ISRFs) responsible for heating the dust. We first fit the radial profiles with exponential functions in order to compare stellar and cool-dust disk scalelengths, as measured by 3.6 μm and 250 μm surface brightnesses. Our results show thatthe stellar and dust scalelengths are comparable, with a mean ratio of 1.04, although several galaxies show dust-to-stellar scalelength ratios of 1.5 or more. We then fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified blackbodies using both variable (MBBV) and fixed (MBBF) emissivity indices β, as well as with physically motivated dust models. The KINGFISH profiles are well suited to examining trends of dust temperature Tdust and β because they span a factor of ~200 in the ISRF intensity heating the bulk of the dust mass, Umin. Results from fitting the profile SEDs suggest that, on average, Tdust, dust optical depth τdust, and Umin decrease with radius. The emissivity index β also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average β ~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with β ≲ 1. An analysis of the wavelength variations of dust emissivities in both the data and the models shows that flatter slopes (β ≲ 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust – β degeneracy. This trend is related to variations in Umin since β and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat β ≲ 1. Both these results strongly suggest that the low apparent β values (flat slopes) in MBBV fits are caused by temperature mixing along the line of sight, rather than by intrinsic variations in grain properties. Finally, a comparison of dust models and the data show a slight ~10% excess at 500 μm for low metallicity (12 + log (O/H) ≲ 8) and low far-infrared surface brightness (Σ500).