Issue |
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
|
|
---|---|---|
Article Number | L70 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201014649 | |
Published online | 16 July 2010 |
Herschel: the first science highlights
LETTER TO THE EDITOR
100
m and 160
m emission as resolved star-formation rate estimators in M 33 (HERM33ES)![[*]](/icons/foot_motif.png)
M. Boquien1 -
D. Calzetti1 -
C. Kramer2 -
E. M. Xilouris3 -
F. Bertoldi4 -
J. Braine5 -
C. Buchbender2 -
F. Combes6 -
F. Israel7 -
B. Koribalski8 -
S. Lord9 -
G. Quintana-Lacaci2 -
M. Relaño10 -
M. Röllig11 -
G. Stacey12 -
F. S. Tabatabaei13 -
R. P. J. Tilanus14 -
F. van der Tak15 -
P. van der Werf7 -
S. Verley16
1 - University of Massachusetts, Department of Astronomy, LGRT-B 619E, Amherst, MA 01003, USA
2 - Instituto Radioastronomia Milimetrica, 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 - Argelander Institut für Astronomie. Auf dem Hügel 71, 53121 Bonn, Germany
5
- Laboratoire d'Astrophysique de Bordeaux, Université Bordeaux 1,
Observatoire de Bordeaux, OASU, UMR 5804, CNRS/INSU, B.P. 89, Floirac
33270, France
6 - Observatoire de Paris, LERMA, 61 Av. de l'Observatoire, 75014 Paris, France
7 - Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
8 - ATNF, CSIRO, PO Box 76, Epping, NSW 1710, Australia
9 - IPAC, MS 100-22 California Institute of Technology, Pasadena, CA 91125, USA
10 - Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, England
11 - KOSMA, I. Physikalisches Institut, Universität zu Köln, Zülpicher Straße 77, 50937 Köln, Germany
12 - Department of Astronomy, Cornell University, Ithaca, NY 14853, USA
13 - Max Planck Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
14 - JAC, 660 North A'ohoku Place, University Park, Hilo, HI 96720, USA
15 - SRON Netherlands Institute for Space Research, Landleven 12, 9747 AD Groningen, The Netherlands
16 - Dept. Física Teórica y del Cosmos, Universidad de Granada, Spain
Received 31 March 2010 / Accepted 21 April 2010
Abstract
Context. Over the past few years several studies have
provided estimates of the SFR (star-formation rate) or the total
infrared luminosity from just one infrared band. However these
relations are generally derived for entire galaxies, which are known to
contain a large scale diffuse emission that is not necessarily related
to the latest star-formation episode.
Aims. We provide new relations to estimate the SFR from resolved star-forming regions at 100 m and 160
m.
Methods. We select individual star-forming regions in the nearby
(840 kpc) galaxy M 33. We estimate the SFR combining the
emission in H
and at 24
m to calibrate the emission at 100
m and 160
m as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key program.
Results. There is less emission in the HII regions at 160 m than at 100
m. Over a dynamic range of almost 2 dex in
we find that the 100
m emission is a nearly linear estimator of the SFR, whereas that at 160
m is slightly superlinear.
Conclusions. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high
,
star formation drives the dust temperature, whereas uncertainties and
variations in radiation-transfer and dust-heated processes dominate at
low
.
Detailed modelling of both galaxies and individual star forming regions
will be needed to interpret similarities and differences between the
two and assess the fraction of diffuse emission in galaxies.
Key words: galaxies: individual: M 33 - galaxies: spiral - infrared: galaxies - galaxies: star formation
1 Introduction
Star formation is one of the main drivers of galaxy formation and
evolution and as such the accuracy of the SFR (star-formation rate)
determination is of great importance for deriving the cosmic history of
galaxies. Along with the UV and the H,
the total infrared luminosity is widely used to estimate the SFR. To
properly quantify the infrared luminosity a good sampling of the
infrared SED (spectral energy distribution) is needed (Dale & Helou 2002; Draine & Li 2007). Over the past few years, many authors have shown that the TIR (total infrared) luminosity, and by extension the SFR (Kennicutt 1998), can also be evaluated from monochromatic emission measures (Rieke et al. 2009; Boquien et al. 2010; Takeuchi et al. 2005; Calzetti et al. 2007,2010,
Li et al. 2010, in prep.). However, most relations between the
far-infrared luminosity and the SFR are established for entire
galaxies. At shorter wavelengths, some relations have been derived from
individual star-forming regions (Relaño et al. 2007; Calzetti et al. 2005; Pérez-González et al. 2006; Calzetti et al. 2007).
As new deep surveys will become available at wavelengths where most of
the energy is emitted, it is prudent to try to understand the relation
between the integral emission of a galaxy and that of the individual
star forming regions. The physical conditions, such as temperature,
abundance and emissivity, of the infrared-emitting dust can vary widely
in a galaxy, and so does the significant contribution from evolved
stars (Sauvage & Thuan 1992; Buat & Xu 1996; Lonsdale Persson & Helou 1987).
Therefore, any scaling relation established from the emission of entire
galaxies may not be appropriate when applied to resolved star-forming
regions in these same galaxies.
![]() |
Figure 1:
PACS maps of M 33 at 100 |
Open with DEXTER |
The limited resolution of far-infrared instruments onboard IRAS, ISO or even Spitzer beyond 60 m made the study of individual star-forming regions within galaxies difficult. However, the recently launched Herschel Space Observatory (Pilbratt et al. 2010)
with its unprecedented resolution provides the first opportunity to
study the spatially resolved far-infrared dust emission in exquisite
detail. Such a fine resolution is of the utmost importance to study the
emission of star-forming regions located in nearby galaxies in order to
provide insights into the fundamental properties of the dust and to
quantify the SFR.
With an inclination of 56
(Regan & Vogel 1994) and a distance of only 840 kpc (Freedman et al. 1991), M 33 is one of the closest spiral galaxies. It has been imaged by Herschel in the context of the HERM33ES open time key program (Kramer et al. 2010), providing one of the finest views Herschel will ever provide of a spiral galaxy from 100
m to 500
m (Verley et al. 2010b; Kramer et al. 2010; Braine et al. 2010).
2 Observations and data reduction
2.1 PACS
The observations provided by PACS (Poglitsch et al. 2010) at 100 m and 160
m are presented by Kramer et al. (2010)
along with a detailed description of the data processing pipeline. The
observations were carried out on 2010-01-07 in parallel mode with a 20
/s
scanning speed for a total of 6.3 h, through a single scan and a
perpendicular cross scan. The frames were first processed to
level 1 with HIPE (Ott 2010), the drifts were corrected and were deglitched with the second-order deglitcher with a 6-
threshold. The maps were produced with photproject mapmaker
using a two-step masking technique to preserve the diffuse emission
from being affected by the high-pass filter. The total flux of the
galaxy is consistent with the measures provided by IRAS and Spitzer at 100
m and 160
m (Kramer et al. 2010). We present the two maps in Fig. 1. The pixel size is 3.2
at 100
m and 6.4
at 160
m for a spatial resolution of 6.7
6.9
at 100
m and 10.7
12.1
at 160
m. The absolute calibration uncertainty is 5% at 100
m and 10% at 160
m. The total fluxes of M 33 agree to within a few percent with those from ISO (Hippelein et al. 2003) and MIPS. In addition, for all radial averages the PACS 160 flux is within 20% of the MIPS 160 one.
2.2 H
and Spitzer MIPS 24
m
We used the H
image presented by Hoopes & Walterbos (2000) that is commonly used in the literature (Verley et al. 2010a,2009; Gardan et al. 2007; Tabatabaei et al. 2007; Verley et al. 2007). The NII contamination was corrected assuming [NII]/H
in the filter bandpass. We also corrected the fluxes for Galactic foreground extinction using the Cardelli et al. (1989) law, assuming
E(B-V)=0.042 from the NASA Extragalactic Database
.
We used the 24 m MIPS data presented by Tabatabaei et al. (2007). No further processing was performed on this image.
![]() |
Figure 2:
|
Open with DEXTER |
2.3 Flux measurements
All targeted HII regions flux densities were measured in polygonal
apertures using IRAF's polyphot procedure. Each polygon was constructed
manually from the PACS 160 m image and tailored to avoid subtraction artefacts in H
and background sources in MIPS 24
m images. Each source was selected to be as compact as possible, taking into account the blending at 160
m
to avoid the mix of several star-forming regions of different ages and
properties. The background was calculated measuring the mode of the
pixels distribution in an annulus around the aperture. Annulus pixels
falling into the aperture of a source were automatically discarded. The
inner radius of the annulus ranges from 20
to 70
by steps of 10
,
which were defined to be larger than the equivalent radius of the aperture:
,
where S is the area of the aperture. A scale of 50 pc corresponds to an angular size of 11
.
The width of the annulus was set to 12
.
Aperture correction was performed for 24
m data
and for the PACS bands
. As the apertures are not circular, we applied the method presented in Boquien et al. (2007) using the equivalent radius of the aperture.
As a proxy for the SFR we applied the scaling presented by Calzetti et al. (2007):
in
yr-1, where
is the H
luminosity in W and,
is defined as
at 24
m in W, assuming a Kroupa (2001) IMF (initial mass function) with a constant SFR over 100 Myr.
For an easier use of the SFR estimator we will provide in this article, we subsequently worked in (luminosity
surface density) because it is distance-independent, in order to
facilitate a comparison with other galaxies. To do so, we divided the
luminosity by the area of the polygon measured in kpc2.
3 Results
3.1 General characteristics of HII regions
We selected a total of 179 HII regions from the 160 m
map. The physical equivalent radius of the extraction apertures ranges
from 37 pc to 256 pc, with a median of 99 pc.
ranges from
W kpc-2 to
W kpc-2 and
from
W kpc-2 to
W kpc-2. The typical 1-
uncertainties are 0.09, 0.06, 0.02, and 0.03 dex in H
,
24
m, 100
m, and 160
m respectively.
(SFR density) ranges from
kpc-2 yr-1 to
kpc-2 yr-1.
The fraction of the total flux enclosed in the 179 apertures compared
to the total flux of M 33 is 0.40, 0.43, 0.35, and 0.24 in H,
24
m, 100
m, and 160
m assuming galaxy-integrated fluxes
W m-2,
F(24)=49.4 Jy (Verley et al. 2007),
F(100)=1288 Jy, and
F(160)=1944 Jy (Kramer et al. 2010). That so little of the total 160
m
flux associated with the selected HII regions suggests that the
large-scale diffuse emission seen at this wavelength may not be
directly related to the ongoing massive star-formation, which is
consistent with the result of Hinz et al. (2004) for M 33, but may be heated by non-ionising B and A stars as observed for instance by Israel et al. (1996) in another galaxy with a similar metallicity, NGC 6822.
3.2 100
m and 160
m as SFR estimators
Being closer to the peak IR emission, the 100 
First of all
and
are very well correlated with a Spearman correlation coefficient
.
This is expected because both bands probe the grey body emission of big grains. It appears that
and
are also well correlated with
,
with a Spearman correlation coefficient
and
.
In Fig. 2 we present the fits of
and
versus the estimated
.
To estimate the relations between the SFR and the PACS emission we
fitted a linear relation in log-log using an ordinary least-square
technique taking into account uncertainties on both axes.
The best fits for M 33 HII regions respectively correspond to
![]() |
(1) | |
![]() |
(2) |
The relations have a dispersion around the best fit of 0.22 dex and 0.25 dex. We notice that the 100



4 Discussion
4.1 Modeling
To model the individual HII regions we used the Calzetti et al. (2007) model as a baseline. The ionising flux is determined using S TARBURST99 (Leitherer et al. 1999) with a Kroupa (2001) IMF, an instantaneous burst and solar metallicity. The extinction is assumed to follow the Calzetti (2001) law. As the
ratio does not show a significant correlation with
for our selected regions in M 33, we assumed a constant
E(B-V)=0.25 mag, set to reproduce the observed mean
.
For the dust emissivity we assumed the Draine & Li (2007) model prescriptions for Spitzer MIPS 160, pending updated dust emissivities for Herschel bands. The model is plotted in Fig. 2.
4.2 Comparison with entire galaxies
Calzetti et al. (2010) showed that for entire galaxies the emission at 160







4.3 Dust temperature
The slight non-linearity in the
relation hints at a higher dust temperature with increasing
,
which would have different effects on the two PACS bands because the
peak of the emission passes through the filter bandpasses as the
temperature increases. Indeed, as dust gets warmer, an increasing
fraction will be emitted at shorter wavelengths. In Fig. 3 we plot
,
a proxy for the dust temperature of the warm component, versus the
.
![]() |
Figure 3:
|
Open with DEXTER |
We see a clear trend with higher
(
)
leading to a higher dust temperature. The spanned range is compatible with the emissivity values published by Draine & Li (2007) for
,
U being the interstellar radiation field normalised to that of
the solar neighbourhood. This means that an increasing fraction of the
total dust emission is detected in the 100
m band compared to the 160
m
band. Interestingly we also notice that the trend in entire galaxies is
very similar to the trend in individual HII regions in M 33
despite the fact that individual regions should have little
contamination from the diffuse large scale emission. One possible
interpretation is that in both entire galaxies and individual HII
regions, star-formation dominates at higher
and creates the trend whereas, at lower
the trend is influenced by the uncertainties and variations of
conditions in the radiation-transfer and dust-heating processes such as
the opacity of the star-formation region, the clumpiness of the media,
the relative locations of stars and dust clouds, etc.
5 Conclusions
We used the high-resolution Herschel 100 m and 160
m observations of a nearby star-forming galaxy, M 33. We combined Herschel PACS data with Spitzer MIPS 24
m and ground-based H
to provide new calibrations of the 100
m and 160
m to estimate
from individual star-forming regions.For the selected star-forming regions in M 33, the 100
m
luminosity is a linear SFR estimator over a factor 100 in surface
brightness, whereas the 160 micron luminosity is slightly superlinear.
It appears that individual star forming regions exhibit a similar
behaviour as entire galaxies taken from the LVL, SINGS, starburst
galaxies from Engelbracht et al. (2008) samples when estimating
from the 100
m and 160
m bands emission. In a similar fashion, the dust temperature - as measured by the ratio of the 160
m to 100
m emission - increases as a function of
,
suggesting that at high
the star formation drives the trend for both systems, while at lower
uncertainties and variations of conditions in the radiation-transfer
and dust-heating processes contribute to the scatter. In other words,
the fairly wide dust temperature distribution at low
becomes increasingly biased towards higher temperatures at higher
in both HII regions and entire galaxies.
We thank Herschel scientists for their valuable help with the PACS data reduction, in particular Babar Ali, Bruno Altieri, Bidushi Bhattacharya, Nicolas Billot and Marc Sauvage. We also thank the NHSC for providing the computing architecture used in the reduction of the data. We also thank our referee, C. K. Xu, for useful comments that helped improve the quality of this article.
References
- Boquien, M., Duc, P.-A., Braine, J., et al. 2007, A&A, 467, 93 Boquien, M., Bendo, G., Calzetti, D., et al. 2010, ApJ, 713, 626 [NASA ADS] [CrossRef] [Google Scholar]
- Braine, J., et al. 2010, A&A, 518, L69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Buat, V., & Xu, C. 1996, A&A, 306, 61 [NASA ADS] [Google Scholar]
- Calzetti, D. 2001, PASP, 113, 1449 [NASA ADS] [CrossRef] [Google Scholar]
- Calzetti, D., Kennicutt, R. C., Bianchi, L., et al. 2005, ApJ, 633, 871 [NASA ADS] [CrossRef] [Google Scholar]
- Calzetti, D., Kennicutt, R. C., Engelbracht, C. W., et al. 2007, ApJ, 666, 870 [NASA ADS] [CrossRef] [Google Scholar]
- Calzetti, D., Wu, S., Hong, S., et al. 2010, ApJ, 714, 1256 [NASA ADS] [CrossRef] [Google Scholar]
- Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Dale, D. A., & Helou, G. 2002, ApJ, 576, 159 [NASA ADS] [CrossRef] [Google Scholar]
- Dale, D. A., Gil de Paz, A., Gordon, K. D., et al. 2007, ApJ, 655, 863 [NASA ADS] [CrossRef] [Google Scholar]
- Dale, D. A., Cohen, S. A., Johnson, L. C., et al. 2009, ApJ, 703, 517 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T., & Li, A. 2007, ApJ, 657, 810 [NASA ADS] [CrossRef] [Google Scholar]
- Engelbracht, C. W., Rieke, G. H., Gordon, K. D., et al. 2008, ApJ, 678, 804 [NASA ADS] [CrossRef] [Google Scholar]
- Freedman, W. L., Wilson, C. D., & Madore, B. F. 1991, ApJ, 372, 455 [NASA ADS] [CrossRef] [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]
- Hinz, J. L., Rieke, G. H., Gordon, K. D., et al. 2004, ApJS, 154, 259 [NASA ADS] [CrossRef] [Google Scholar]
- Hippelein, H., Haas, M., Tuffs, R. J., et al. 2003, A&A, 407, 137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hoopes, C. G., & Walterbos, R. A. M. 2000, ApJ, 541, 597 [NASA ADS] [CrossRef] [Google Scholar]
- Israel, F. P., Bontekoe, T. R., & Kester, D. J. M. 1996, A&A, 308, 723 [NASA ADS] [Google Scholar]
- Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 [Google Scholar]
- Kramer, C., et al. 2010, A&A, 518, L67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
- Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Lonsdale Persson, C. J., & Helou, G. 1987, ApJ, 314, 513 [NASA ADS] [CrossRef] [Google Scholar]
- Magrini, L., Stanghellini, L., Corbelli, E., Galli, D., & Villaver, E. 2010, A&A, 512, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ott, S. 2010, in Astronomical Data Analysis Software and Systems XIX, ed. Y. Mizumoto, K.-I. Morita, & M. Ohishi, ASP Conf. Ser. [Google Scholar]
- Pérez-González, P. G., Kennicutt, Jr., R. C., Gordon, K. D., et al. 2006, ApJ, 648, 987 [NASA ADS] [CrossRef] [Google Scholar]
- Pilbratt, G. L., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
- Poglitsch, A., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Regan, M. W., & Vogel, S. N. 1994, ApJ, 434, 536 [NASA ADS] [CrossRef] [Google Scholar]
- Relaño, M., Lisenfeld, U., Pérez-González, P. G., Vílchez, J. M., & Battaner, E. 2007, ApJ, 667, L141 [NASA ADS] [CrossRef] [Google Scholar]
- Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., et al. 2009, ApJ, 692, 556 [NASA ADS] [CrossRef] [Google Scholar]
- Sauvage, M. & Thuan, T. X. 1992, ApJ, 396, L69 [NASA ADS] [CrossRef] [Google Scholar]
- Tabatabaei, F. S., Beck, R., Krause, M., et al. 2007, A&A, 466, 509 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Takeuchi, T. T., Buat, V., Iglesias-Páramo, J., Boselli, A., & Burgarella, D. 2005, A&A, 432, 423 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Verley, S., Hunt, L. K., Corbelli, E., & Giovanardi, C. 2007, A&A, 476, 1161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Verley, S., Corbelli, E., Giovanardi, C., & Hunt, L. K. 2009, A&A, 493, 453 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Verley, S., Corbelli, E., Giovanardi, C., & Hunt, L. K. 2010a, A&A, 510, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Verley, S., et al. 2010b, A&A, 518, L68 [NASA ADS] [CrossRef] [EDP Sciences] [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.
- ... mapmaker
- PACS photometer - Prime and Parallel scan mode release note. V.1.2, 23 February 2010.
- ... Database
- The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.
- ... data
- Following the formula provided in the MIPS instrument handbook.
- ... bands
- Following the correction provided in PACS photometer - Prime and Parallel scan mode release note. V.1.2, 23 February 2010.
All Figures
![]() |
Figure 1:
PACS maps of M 33 at 100 |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
|
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.