Issue |
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
|
|
---|---|---|
Article Number | L69 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201014551 | |
Published online | 16 July 2010 |
Herschel: the first science highlights
LETTER TO THE EDITOR
Cool gas and dust in M 33: Results from the
HERschel M 33 Extended Survey (HERM33ES)![[*]](/icons/foot_motif.png)
J. Braine1 - P. Gratier1 - C. Kramer2 - E. M. Xilouris3 - E. Rosolowsky4 - C. Buchbender2 - M. Boquien5 - D. Calzetti5 - G. Quintana-Lacaci2 - F. Tabatabaei6 - S. Verley7 - F. Israel8 - F. van der Tak9 - S. Aalto10 - F. Combes11 - S. Garcia-Burillo12 - M. Gonzalez2 - C. Henkel6 - B. Koribalski13 - B. Mookerjea14 - M. Roellig15 - K. F. Schuster16 - M. Relaño17 - F. Bertoldi18 - P. van der Werf8 - M. Wiedner11
1 - Laboratoire d'Astrophysique de Bordeaux, Université Bordeaux 1,
Observatoire de Bordeaux, OASU, UMR 5804, CNRS/INSU, BP 89, 33270 Floirac, France
2 -
Instituto Radioastronomia Milimetrica (IRAM),
Av. Divina Pastora 7, Nucleo Central, 18012 Granada, Spain
3 -
Institute of Astronomy and Astrophysics, National Observatory of Athens,
P. Penteli, 15236 Athens, Greece
4 -
University of British Columbia, Okanagan, 3333 University Way, Kelowna BC V1V 1V7 Canada
5 -
Department of Astronomy, University of Massachusetts, Amherst,
MA 01003, USA
6 -
Max-Planck-Institut für Radioastronomie (MPIfR),
Auf dem Hügel 69, 53121 Bonn, Germany,
Argelander Institut fr Astronomie. Auf dem Hügel 71,
53121 Bonn, Germany
7 -
Dept. Física Teórica y del Cosmos, Universidad de Granada, Spain
8 -
Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
9 -
SRON Netherlands Institute for Space Research, Landleven 12,
9747 AD Groningen, The Netherlands
10 -
Onsala Space Observatory, Chalmers University of Technology,
43992 Onsala, Sweden
11 -
Observatoire de Paris, LERMA, 61 Av. de l'Observatoire,
75014 Paris, France
12 -
Observatorio Astronómico Nacional (OAN) - Observatorio de Madrid,
Alfonso XII 3, 28014 Madrid, Spain
13 -
Australia Telescope National Facility, CSIRO, PO Box 76, Epping,
NSW 1710, Australia
14 -
Department of Astronomy & Astrophysics,
Tata Institute of Fundamental Research,
Homi Bhabha Road, 400005 Mumbai, India
15 -
KOSMA, I. Physikalisches Institut, Universität zu Köln,
Zülpicher Straße 77, 50937 Köln, Germany
16 -
IRAM, 300 rue de la Piscine, 38406 St. Martin d'Hères, France
17 -
Institute of Astronomy, University of Cambridge, Madingley Road,
Cambridge CB3 0HA, UK
18 -
Argelander Institut fr Astronomie. Auf dem Hügel 71,
53121 Bonn, Germany
Received 30 March 2010 / Accepted 11 May 2010
Abstract
We present an analysis of the first space-based far-IR-submm
observations of M 33, which measure the emission from the cool
dust and resolve the giant molecular cloud complexes. With roughly
half-solar abundances, M 33 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 H2 content. In this Letter, we measure the dust emission cross-section
using SPIRE and recent CO and H I observations; a variation in
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 H I column, yields a morphology similar to that observed in CO. The H2/H I
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 M 33.
Key words: galaxies: individual: M 33 - Local Group - galaxies: evolution - galaxies: ISM - ISM: clouds - stars: formation
1 Introduction, data, and dust temperature
Understanding star formation requires studying the interplay between
the phases of the interstellar medium (ISM). Dust processes most of the
energy transiting the ISM, but the cool dust component, although
representing the vast majority of the dust mass, is difficult to
observe from the ground.
Herschel SPIRE observations (Griffin 2010; Pilbratt 2010) are the first space-based 250-500 m
data and as such provide a unique occasion to put together a global
picture of the cool gas and dust in M 33. In particular, we
compare the morphology of the Far-IR emission and that of the gas as
determined from CO and H I measurements and attempt to measure how the dust cross-section varies in M 33.
A longer term goal is to be able to use the dust emission to constrain the variation of the
factor within M 33 and elsewhere.
![]() |
Figure 1:
Dust color temperature map of M 33 calculated from the SPIRE 250 and 350 |
Open with DEXTER |
![]() |
Figure 2:
Radial distribution of the average dust temperature in the stellar disk
of M 33. The solid line joining the triangles shows the
250/350 |
Open with DEXTER |
This Letter is one of a series on the HERM33ES project on the ISM of the Local Group galaxy M 33, an overview of which is given in Kramer et al. (2010), hereafter K10. For consistency with the other M 33 papers in this volume, we adopt a distance of D = 840 kpc for M 33 (i.e., 25 arcsec = 100 pc) and orientation parameters of PA = 22.5 degrees and i = 56.
We use the recent HERM33ES SPIRE observations at 250, 350, and 500
m combined with CO(2-1) observations from the M 33CO@IRAM project, as a tracer of the
molecular component, and a high-resolution mosaic of VLA HI data (both from Gratier et al. 2010).
The SPIRE data were first processed as described in K10 and then
converted from Jy/beam to brightness units (MJy/sr).
To estimate the dust temperature, the 250
m data were convolved to the 350
m beamsize (
25'') and, assuming a single temperature grey body with an emissivity
with
,
a temperature was derived from the flux ratio.
At the temperatures of the cool component seen in M 33, Dupac et al. (2003) find
.
To minimize the effect of the uncertainty in
,
we chose to use adjacent bands to estimate temperature. The 250/350
m
ratio provides more accurate temperatures than the 350/500 micron
ratio even for temperatures below 10 K. A 15% variation (or
uncertainty) in the 250/350
m
ratios corresponds to a temperature change of 0.9, 2.1, and 4.1 K
at temperatures of 10, 15, and 20 K but the same uncertainty in
the 350/500
m
ratios yields 1.3, 3.2, and 6.3 K errors for the same dust
temperatures. A further advantage is that we obtain the dust
temperature at a resolution typical of giant molecular clouds (GMC),
100 pc.
Ideally, a multi-component fit would be used but this requires high S/N
data at many wavelengths over the whole disk. By assuming that the cool
dust component dominates the
emission beyond 250 m,
we could calculate the dust temperature out to the optical radius of
M 33. The resulting temperature map is shown in Fig. 1. Figure 2
shows a comparison with the temperature of the cool component of the
preferred two-component (warm plus cool dust) model fit to data between
24
m and 500
m from K10.
The two-component model uses
and the temperatures are higher out to 6 kpc. The temperatures in
the radial bins have been estimated by averaging the temperatures
rather than averaging the emission as in K10. The latter yields
slightly higher temperatures because the dust is usually warmer where
emission is strong (compare Fig. 1 with Fig 1. in K10).
K10 calculate the dust column density using
g-1 of dust, which is equivalent to a dust cross-section per H-atom of
for a hydrogen gas-to-dust mass ratio of 140.
In the following we use the SPIRE 250/350
m
color temperature because the data cover a greater area and agree well
with the cool dust temperature of the two-component fit, showing the
domination of the cool dust at these wavelengths.
No correction for line contamination was subtracted from the SPIRE
data. At these frequencies, the CO lines contribute very little to the
continuum flux, unlike at 1.3 mm or 850
m (e.g. Braine et al. 1997).
Figure 3 shows the SPIRE 250 m
emission with the CO(2-1) emission as contours at 25'' resolution to
show how closely the CO emission follows the dust emission peaks. Gratier et al. (2010) show a detailed comparison of the CO and H I emission with star formation tracers such as H
and Spitzer 8 and 24
m maps.
The 250
m-bright regions are detected in CO and the general morphology of the cool dust emission is seen in the H I image.
![]() |
Figure 3:
SPIRE 250 |
Open with DEXTER |
2 Dust cross-section
In order to estimate the total gas mass of a galaxy from its dust
emission, it is first necessary to measure the dust emission
cross-section per H-atom .
We do this by selecting regions with HI emission and well-constrained
dust temperatures but with little or no CO emission (in order to avoid the assumption of a
value and CO line ratio). Thereby, we can relate the observed neutral hydrogen column density,
,
to the dust optical depth,
,
and infer the cross section
.
In practice, two methods were used: (i) we selected large contiguous regions as shown on the map in Fig. 3 and (ii) all pixels with HI emission, a defined dust temperature, and no detected CO emission (
K km s-1) were taken, using maps at the resolution of the 350
m SPIRE data.
In the optically thin limit,
![]() |
(1) |
where





From the regions indicated in Fig. 3, we estimate
(using Eq. (1)) for different radii and the results are shown in the first line of Table 1. We assume that the cross section
within the inner 4 kpc is the same as at 4 kpc for two reasons: (i) there are few truly CO-free regions in the central part of M 33 and (ii) Gardan et al. (2007) showed (their Fig. 13) that the H I-H2 relation was similar within 4kpc but beyond 4 kpc the H2 fraction decreased sharply even at constant total Hydrogen column density, so we consider
kpc
as a transition in ISM properties. It is very likely that some
molecular gas is present in these regions, particularly those closer to
the center, so
is likely overestimated. In
the second approach, in order to exclude averaging in undetected molecular clouds, we take the peak of the histogram of
values for the areas without detected CO emission but with HI emission
and a constrained dust temperature, yielding the lower values for
,
as shown in Table 1.
Using this method, we find a north-south difference in M 33 at equivalent radii.
The uncertainty in the histogram method is about
;
although the distribution is sometimes broad,
the peak value is well defined. The histograms of
values are clearly different in the north and south.
The polygon method was difficult to use in the south due to the absence
of large CO-free regions with HI and dust emission. An advantage of the
histogram method is that it excludes the tail of high
which may be due to H2 without detected CO emission.
The intrinsic expectation is that because the oxygen abundance is about half solar, the dust cross-section should be as well (Draine et al. 2007), and that given
the shallow abundance gradient in M 33 (Magrini et al. 2009; Rosolowsky et al. 2007),
this should hold for the entire galaxy. However, it quickly became
apparent that the dust cross-section per H-atom varied over the galaxy,
decreasing with radius by close to a factor 2, with a higher
abundance in the south - similar to the variation found by Magrini et al. (2010) from optical H II region and PNe abundance measurements. In the solar neighborhood,
is about
cm2per H-atom (Draine & Li 2007; Draine & Lee 1984) at 250
m and varying as
.
If we apply a
of half this value to M 33 using the dust temperature map in Fig. 1, we obtain unrealistically high gas masses in the inner disk:
40
pc-2 on average, whereas the H I is about 11
pc-2 and the H2 considerably less (Gratier et al. 2010).
For a half-solar
,
the (cool) dust temperature would have to be about 30K, beyond any experimental uncertainties, so we conclude that
must be higher in the inner disk.
In the outer disk, however, the H I column density (where no CO emission is present) is roughly equal to what is derived from the dust with
,
suggesting that this value is appropriate for the outer disk.
Table 1:
Dust cross-section
at 250
m as a function of radius in M 33, expressed in units of
cm2 per H-atom.
![]() |
Figure 4:
H2 column density map estimated from the 500 |
Open with DEXTER |
3 Total and molecular gas mass
To really obtain an accurate gas mass from the cold dust emission data,
the calibration uncertainties in the SPIRE bands must be reduced
because the 15% uncertainty is not sufficient to distinguish between
different solutions of
and temperature. A second requirement is that the whole disk of
M 33 must be observed in CO so that a well-determined variation of
,
based on the second technique, can be constructed. The H I data provide a reliable picture of the H I
column density, so subtracting this from the total gas column density
derived from the dust emission yields an estimate of the molecular gas
mass which does not rely on CO emission. A further caveat is that
must not change from the atomic to molecular medium. Once the [C I] and [C II]
lines have been measured, these uncertainties can be reduced by the
additional constraint that the C/H ratio (from the sum of the C
reservoirs) should be equal to the C abundance derived by other means.
![]() |
Figure 5:
Scatter plot of H2
column density as derived from the dust emission and the CO(2-1)
intensity. The solid red line is a fit to the entire dataset and
corresponds to a
|
Open with DEXTER |
Using the information currently available, we have used Table 1 (second method) to build an image of the dust cross-section ,
erring on the low-
side, as indicated in the last line of Table 1. We then used this
,
scaled to 500
m with
,
and the dust temperature (Fig. 1) to estimate the total gas mass from the 500
m data and then subtracted the H I column density to obtain the H2 column density using the interferometric H I maps from Gratier et al. (2010), which recover more than 90% of the flux found by Putman et al. (2009)
using Arecibo. The missing 21cm flux is expected to be located in the
mid-to-outer disk where the rotation curve is fairly flat. A
morphological comparison of the dust-derived H2 column and the CO contours in Fig. 4 shows that while major regions were not missed by the CO observations, the dust-derived H2 map does not correlate perfectly with the CO.
For example, the strong CO peak at (01h34
09
4,+30
49'06'') corresponds to a red (not white) region in Fig. 4; this may of course reflect a real variation in the
ratio.
The data used in the scatter plot in Fig. 5 (see also Fig. 7 in Leroy et al. 2009) yield a formal fit with a
H2 per K km s-1 value, corresponding to a roughly galactic
H2 per K km s-1 for a line ratio CO(
)
= 0.7. The inner disk
is lower and rises to about
H2 per K km s-1 beyond 2 kpc. The offset, in principle, implies that some low column density H2
is not seen in CO. Both values depend strongly on the details of the
dust cross-section, illustrating the need for a complete map of
M 33 in CO.
The total hydrogen mass derived from the dust emission within 8 kpc is estimated to be
,
very similar to the more directly measured H I+H2 masses in Gratier et al. (2010).
The dust cross-sections we derive suggest dust-to-gas mass ratios (including He) ranging from
125 in the inner 4 kpc to 200 near R25 for Milky Way like dust (cf. K10). Comparing the total dust-derived H2 mass with the CO
luminosity (Gratier et al. 2010) yields a
of about 1.5 times the Galactic
.
Because the H I dominates, however, a small error in the total dust-derived gas mass translates into a large uncertainty in the dust-derived H2 mass. It is clearly necessary to obtain whole-galaxy measurements of
in CO-free zones to refine these estimates.
References
- Braine, J., Guélin, M., Dumke, M., et al. 1997, A&A, 326, 963 Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 [Google Scholar]
- Draine, B. T., & Li, A. 2007, ApJ, 657, 810 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, 663, 866 [NASA ADS] [CrossRef] [Google Scholar]
- Dupac, X., Bernard, J., Boudet, N., et al. 2003, A&A, 404, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gardan, E., Braine, J., Schuster, K. F., Brouillet, N., & Sievers, A. 2007, A&A, 473, 91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gratier, P., Braine, J., Rodriguez-Fernandez, N. J., et al. 2010, A&A, accepted [arXiv:1003.3222] [Google Scholar]
- Griffin, M. J., et al. 2010, A&A, 518, L3 [Google Scholar]
- Kramer, C., et al. 2010, A&A, 518, L67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Leroy, A. K., Bolatto, A., Bot, C., et al. 2009, ApJ, 702, 352 [NASA ADS] [CrossRef] [Google Scholar]
- Magrini, L., Stanghellini, L., Corbelli, E., Galli, D., & Villaver, E. 2010, A&A, 512, 63 [Google Scholar]
- Magrini, L., Stanghellini, L., & Villaver, E. 2009, ApJ, 696, 729 [NASA ADS] [CrossRef] [Google Scholar]
- Pilbratt, G. L., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
- Putman, M. E., Peek, J. E. G., Muratov, A., et al. 2009, ApJ, 703, 1486 [NASA ADS] [CrossRef] [Google Scholar]
- Rosolowsky, E., Keto, E., Matsushita, S., & Willner, S. P. 2007, ApJ, 661, 830 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... (HERM33ES)
- Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
All Tables
Table 1:
Dust cross-section
at 250
m as a function of radius in M 33, expressed in units of
cm2 per H-atom.
All Figures
![]() |
Figure 1:
Dust color temperature map of M 33 calculated from the SPIRE 250 and 350 |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Radial distribution of the average dust temperature in the stellar disk
of M 33. The solid line joining the triangles shows the
250/350 |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
SPIRE 250 |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
H2 column density map estimated from the 500 |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Scatter plot of H2
column density as derived from the dust emission and the CO(2-1)
intensity. The solid red line is a fit to the entire dataset and
corresponds to a
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.