S. Verley

DUST HEATING SOURCES IN GALAXIES: THE CASE OF M33 (HERM33ES)

Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 mu m. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 mu m to 500 mu m, along with H I, H alpha maps, and Galaxy Evolution Explorer ultraviolet data, we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination-corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations.

Star formation in M 33: multiwavelength signatures across the disk

Aims. We use different tracers, such as Hα, ultraviolet (UV), and infrared (IR) emissions at various wavelengths, to study the dust and star formation (SF) conditions throughout the disk of M 33. Methods. We derive the radial distribution of dust, of the old and young stellar population using Spitzer and GALEX data, complemented by ground-based optical data and available surveys of atomic and molecular gas. We separate the contribution of discrete sources to the IR brightness from the diffuse emission. Results. At 8 and 24 μm, discrete sources account for >40% of the IR emission in the innermost 3 kpc, and for <20% further out. We find that stochastic emission from very small grains in the diffuse interstellar medium accounts for only ∼10% of the diffuse 24 μm emission, and that dusty circumstellar shells of unresolved, evolved AGB stars (carbon stars) are a viable alternative. The 8 μ mp rof ile suggests that PAH emission declines faster with radius than the dust continuum. In annular regions 0.24 kpc wide, we find a mean extinction value for the stellar continuum AV ∼ 0.25 mag with a weak dependence on radius, consistent with the shallow metallicity gradient observed. Dust opacity derived from the 160 μm emission decreases instead by a factor of 10 from the center to edge of the star forming disk. Conclusions. Using extinction corrected UV and Hα maps we find the global SF rate in M 33 over the last 100 Myr to be

PACS and SPIRE photometer maps of M33: First results of the HERschel M33 Extended Survey (HERM33ES)

Context. Within the framework of the HERM33ES key program, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. Aims. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 mu m wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Methods. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 mu m and 500 mu m using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. Results. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 mu m). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60 K +/- 10 K. The temperature of the cold component drops significantly from similar to 24 K in the inner 2 kpc radius to 13 K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta = 1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

The AMIGA sample of isolated galaxies - IX. Molecular gas properties

Aims. We characterize the molecular gas content (ISM cold phase) using CO emission of a redshift-limited subsample of isolated galaxies from the AMIGA (Analysis of the interstellar Medium of Isolated GAlaxies) project in order to provide a comparison sample for studies of galaxies in different environments. Methods. We present the 12 CO(1–0) data for 273 AMIGA galaxies, most of them (n = 186) from our own observations with the IRAM 30 m and the FCRAO 14 m telescopes and the rest from the literature. We constructed a redshift-limited sample containing galaxies with 1500 km s −1 <v< 5000 km s −1 and excluded objects with morphological evidence of possible interaction. This sample (n = 173) is the basis for our statistical analysis. It contains galaxies with molecular gas masses, MH2 , in the range of ∼10 8 −10 10 M� .I t is dominated, both in absolute number and in detection rate, by spiral galaxies of type T = 3–5 (Sb-Sc). Most galaxies were observed with a single pointing towards their centers. Therefore, we performed an extrapolation to the total molecular gas mass expected in the entire disk based on the assumption of an exponential distribution. We then studied the relationships between MH2 and other galactic properties (LB, D 2 , LK, LFIR ,a ndMHI). Results. We find correlations between MH2 and LB, D 2 , LK ,a ndLFIR. The tightest correlation of MH2 holds with LFIR and, for T = 3– 5, with LK, and the poorest with D 2 . The correlations with LFIR and LK are very close to linearity. The correlation with LB is nonlinear so that MH2 /LB increases with LB. The molecular and the atomic gas masses of our sample show no strong correlation. We find a low mean value, log(MH2 /MHI) = −0.7 (for T = 3–5), and a strong decrease in this ratio with morphological type. The molecular gas column density and the surface density of the star formation rate (the Kennicutt-Schmidt law) show a tight correlation with a rough unity slope. We compare the relations of MH2 with LB and LK found for AMIGA galaxies to samples of interacting galaxies from the literature and find an indication for an enhancement of the molecular gas in interacting galaxies of up to 0.2–0.3 dex.

Star formation in M 33: the radial and local relations with the gas

Aims. In the Local Group spiral galaxy M 33, we investigate the correlation between the star formation rate (SFR) surface density, ΣSFR, and the gas density Σgas (molecular, atomic, and total). We also explore whether there are other physical quantities, such as the hydrostatic pressure and dust optical depth, which establish a good correlation with ΣSFR. Methods. We use the Hα, far-ultraviolet (FUV), and bolometric emission maps to infer the SFR locally at different spatial scales, and in radial bins using azimuthally averaged values. Most of the local analysis is done using the highest spatial resolution allowed by gas surveys, 180 pc. The Kennicutt-Schmidt (KS) law, ΣSFR ∝ Σ n is analyzed by three statistical methods. Results. At all spatial scales, with Hα emission as a SFR tracer, the KS indices n are always steeper than those derived with the FUV and bolometric emissions. We attribute this to the lack of Hα emission in low luminosity regions where most stars form in small clusters with an incomplete initial mass function at their high mass end. For azimuthally averaged values the depletion timescale for the molecular gas is constant, and the KS index is nH2 = 1.1 ± 0.1. Locally, at a spatial resolution of 180 pc, the correlation between ΣSFR and Σgas is generally poor, even though it is tighter with the molecular and total gas than with the atomic gas alone. Considering only positions where the CO J = 1−0 line is above the 2-σ detection threshold and taking into account uncertainties in ΣH2 and ΣSFR, we obtain a steeper KS index than obtained with radial averages: nH2 = 2.22 ± 0.07 (for FUV and bolometric SFR tracers), flatter than that relative to the total gas (nHtot = 2.59 ± 0.05). The gas depletion timescale is therefore larger in regions of lower ΣSFR. Lower KS indices (nH2 = 1.46 ± 0.34 and nH2 = 1.12) are found using different fitting techniques, which do not account for individual position uncertainties. At coarser spatial resolutions these indices get slightly steeper, and the correlation improves. We find an almost linear relation and a better correlation coefficient between the local ΣSFR and the ISM hydrostatic pressure or the gas volume density. This suggests that the stellar disk, gravitationally dominant with respect to the gaseous disk in M 33, has a non-marginal role in driving the SFR. However, the tight local correlation that exists between the dust optical depth and the SFR sheds light on the alternative hypothesis that the dust column density is a good tracer of the gas that is prone to star formation.

Dust and gas power spectrum in M 33 (HERM33ES)

Power spectra of deprojected images of late-type galaxies in gas or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power laws, a shallow one on large scales (larger than 500 pc) and a steeper one on small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. The break separates the 3D behavior of the interstellar medium on small scales, controlled by star formation and feedback, from the 2D behavior on large scales, driven by density waves in the disk. The break between these two regimes depends on the thickness of the plane, which is determined by the natural self-gravitating scale of the interstellar medium. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M 33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The wide dynamical range (2–3 dex in scale) of most images allows us to clearly determine the change in slopes from −1.5 to −4, with some variations with wavelength. The break scale increases with wavelength from 100 pc at 24 and 100 μm to 350 pc at 500 μm, suggesting that the cool dust lies in a thicker disk than the warm dust, perhaps because of star formation that is more confined to the plane. The slope on small scales tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, such as M 33, are carried out to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power spectra, with two power-law slopes and breaks, that nicelybracket the data. The small-scale power-law does indeed reflect the 3D behavior of the gas layer, steepening strongly while the feedback smoothes the structures by increasing the gas turbulence. M 33 appears to correspond to a fiducial model with an SFR of ~ 0.7 M_⊙/yr, with 10% supernovae energy coupled to the gas kinematics.

Cool and warm dust emission from M 33 (HerM33es)

In the framework of the open-time key program “Herschel M 33 extended survey (HerM33es)”, we study the far-infrared emission from the nearby spiral galaxy M 33 in order to investigate the physical properties of the dust such as its temperature and luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350, and 500 μm) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 μm and Spitzer-MIPS 24 and 70 μm data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the “cool” dust grains are heated to temperatures of between 11 K and 28 K, with the lowest temperatures being found in the outskirts of the galaxy and the highest ones both at the center and in the bright HII regions. The infrared/submillimeter total luminosity (5–1000 μm) is estimated to be 1.9 × 10^9 _(-4.4×10)^8^(+4.0×10)^8L_⊙. Fifty-nine percent of the total infrared/submillimeter luminosity of the galaxy is produced by the “cool” dust grains (~15 K), while the remaining 41% is produced by “warm” dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly higher at the center than the outer parts of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms), we find that the fraction of the emission from the disk in the mid-infrared (24 μm) is 21%, while it gradually rises up to 57% in the submillimeter (500 μm). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The “cool” dust inside the disk is heated to temperatures in a narrow range between 18 K and 15 K (going from the center to the outer parts of the galaxy).

Variation in the dust emissivity index across M 33 with Herschel and Spitzer (HerM 33es)

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M 33 using Spitzer and Herschel photometric data. M 33 is a Local Group spiral with slightly subsolar metallicity, which makes it an ideal stepping-stone to less regular and lower-metallicity objects such as dwarf galaxies and, probably, young-universe objects. Expressing the emissivity of the dust as a power law, the power-law exponent (β) was estimated from two independent approaches designed to properly treat the degeneracy between β and the dust temperature (T ). Both β and T are higher in the inner than in the outer disk, contrary to reported β − T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Hα) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate way. Both cold- and warm-dust column densities are high in star-forming regions and reach their maxima toward the giant star-forming complexes NGC 604 and NGC 595. β declines from close to 2 in the center to about 1.3 in the outer disk. β is positively correlated with star formation and with the molecular gas column, as traced by the Hα and CO emission. The lower dust-emissivity index in the outer parts of M 33 is most likely related to the reduced metallicity (different grain composition) and possibly to a different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the far-infrared-bright regions in the outer disk also have a low β. Like most spirals, M 33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Hα or CO tend to trace the inner disk, which makes it difficult to distinguish between their effects on the dust. The assumption of a constant emissivity index β is obviously not appropriate.

Cool gas and dust in M 33: Results from the HERschel M 33 Extended Survey (HERM33ES)

We present an analysis of the first space-based far-IR-submm observations of M33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H-2 content. In this Letter, we measure the dust emission cross-section sigma using SPIRE and recent CO and HI observations; a variation in s is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the HI column, yields a morphology similar to that observed in CO. The H-2/HI mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M33.

Properties of compact 250 μm emission and H II regions in M 33 (HERM33ES)

Aims. Within the framework of the HERM33ES key program, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250 μm in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Hα shells in the northern part of the galaxy. Methods. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250 μm SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24 μm and Hα for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Hα, far-ultraviolet from GALEX, and infrared from both Spitzer IRAC 8 μm and MIPS 24 μm in order to make comparisons. Results. For the 159 compact sources selected at 250 μm, we find a very strong Pearson correlation coefficient with the MIPS 24 μm emission (r_(24) = 0.94) and a rather strong correlation with the Hα emission, although with more scatter (r_(Hα) = 0.83). The morphological study of the Hα shells shows a displacement between far-ultraviolet, Hα, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Hα shell structures. Conclusions. The very strong link between the 250 μm compact emission and the 24 μm and Hα emissions, by recovering the star formation rate from standard recipes for H II regions, allows us to provide star formation rate calibrations based on the 250 μm compact emission alone. The different locations of the Hα and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Hα shells and the associated cool dust.