Issue |
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
|
|
---|---|---|
Article Number | L89 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201014541 | |
Published online | 16 July 2010 |
Herschel: the first science highlights
LETTER TO THE EDITOR
Determining dust temperatures and masses in the Herschel
era: The importance of observations longward of 200 micron![[*]](/icons/foot_motif.png)
K. D. Gordon1 - F. Galliano2 - S. Hony2 - J.-P. Bernard3 - A. Bolatto4 - C. Bot5 - C. Engelbracht6 - A. Hughes7,8 - F. P. Israel9 - F. Kemper10 - S. Kim11 - A. Li12 - S. C. Madden2 - M. Matsuura13,14 - M. Meixner1,15 - K. Misselt6 - K. Okumura2 - P. Panuzzo2 - M. Rubio16 - W. T. Reach17,18 - J. Roman-Duval1 - M. Sauvage2 - R. Skibba6 - A. G. G. M. Tielens9
1 - Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA
2 - CEA, Laboratoire AIM, Irfu/SAp, Orme des Merisiers, 91191 Gif-sur-Yvette, France
3 - Centre d'Étude Spatiale des Rayonnements, CNRS, 9 av. du Colonel Roche, BP 4346, 31028 Toulouse, France
4 - Department of Astronomy, University of Maryland. College Park, MD 20742, USA
5 - Observatoire Astronomique de Strasbourg, 11 rue de l'université, 67000 STRASBOURG, France
6 - Steward Observatory, University of Arizona, 933 North Cherry Ave., Tucson, AZ 85721, USA
7 - Centre for Supercomputing and Astrophysics, Swinburne University of Technology, Hawthorn VIC 3122, Australia
8 - CSIRO Australia Telescope National Facility, PO Box 76, Epping NSW 1710, Australia
9 - Sterrewacht Leiden, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
10 - Jodrell Bank Centre for Astrophysics, Alan Turing Building,
School of Physics & Astronomy, University of Manchester, Oxford
Road, Manchester M13 9PL, UK - Astronomy & Space Science, Sejong
University, 143-747, Seoul, South Korea
11 - 314 Physics Building, Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA
12 - Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK
13 - MSSL, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK
14 - Visiting Scientist at Smithsonian Astrophysical Observatory, Harvard-CfA, 60 Garden St., Cambridge, MA, 02138, USA
15 - Departamento de Astronomia, Universidad de Chile, Casilla 36-D, Santiago, Chile
16 - Spitzer Science Center, California Institute of Technology, MS 220-6, Pasadena, CA 91125, USA
17 - Stratospheric Observatory for Infrared Astronomy, Universities
Space Research Association, Mail Stop 211-3, Moffett Field, CA 94035,
USA
Received 30 March 2010 / Accepted 1 May 2010
Abstract
Context. The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (100
m). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500
m that is beyond the peak (
160
m) of dust emission in most galaxies.
Aims. We investigate the differences in the fitted dust temperatures and masses determined using only <200 m data and then also including >200
m data (new SPIRE observations) to determine how important having >200
m data is for deriving these dust properties.
Methods. We fit the 100 to 350 m
observations of the Large Magellanic Cloud (LMC) point-by-point with a
model that consists of a single temperature and fixed emissivity law.
The data used are existing observations at 100 and 160
m (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE key project as part of the Herschel science demonstration phase.
Results. The dust temperatures and masses computed using only 100 and 160 m data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350
m data. We find that an emissivity law proportional to
minimizes the 100-350
m fractional residuals. We find that the emission at 500
m is
10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500
m excess is weakly anti-correlated with MIPS 24
m flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500
m
excess. Our results do not allow us to distinguish between a systematic
variation in the wavelength dependent emissivity law or a population of
very cold dust only detectable at
for the origin of the 500
m excess.
Key words: ISM: general - galaxies: individual: LMC - Magellanic Clouds - infrared: ISM
1 Introduction
Among nearby galaxies, the Large Magellanic Cloud (LMC) and Small
Magellanic Cloud (SMC) represent unique astrophysical laboratories for
interstellar medium (ISM) studies. Both Clouds are relatively nearby,
the LMC at 50 kpc (Schaefer 2008) and the SMC at
60 kpc
(Hilditch et al. 2005), and provide ISM measurements that are relatively
unconfused by line-of-sight uncertainties when compared to the Milky
Way. The two Clouds span an interesting metallicity range with the
LMC at
1/2
(Russell & Dopita 1992) being above the
threshold of 1/3-1/4
where the properties of the ISM
change as traced by the rapid reduction in the PAH dust mass fractions
and possible dust-to-gas ratios (Draine et al. 2007) and
the SMC at
1/5
(Russell & Dopita 1992) below this threshold.
Finally, the dust in the LMC and SMC shows strong variations in its
ultraviolet characteristics (Gordon et al. 2003).
![]() |
Figure 1:
The best fit (assuming an emissivity law with
|
Open with DEXTER |
The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) in
the Magellanic Clouds Herschel key program will map both Clouds using
the PACS/SPIRE Parallel observing mode providing observations at 100,
160, 250, 350, and 500 m (Meixner et al. 2010).
The HERITAGE wavelength coverage (100-500
m) and spatial
resolution (
10 pc at 500
m) is well suited to measuring
the spatial variations of dust temperatures and masses. The infrared
dust emission in most galaxies peaks between 100-200
m (Dale et al. 2005) and observations >200
m are important for
accurate dust temperature and masses (Willmer et al. 2009). Ground-based
submilimeter observations do provide the needed >200
m observations, but they have been seen to be in excess of that expected
from extrapolating fits to the <200
m data for sub-solar
metallicity galaxies (Galliano et al. 2005). This excess could be due to
very cold dust that only emits at submilimeter wavelengths or
variations in the wavelength dependent dust emissivity law
(Paradis et al. 2009a; Reach et al. 1995). As part of the science demonstration
program (SDP), two HERITAGE AORs centered on the LMC were executed.
These observations are used in this paper to explore the impact SPIRE
observations have on the measurement of dust temperatures and masses
including the behavior of any submilimeter excess.
2 Data
The observation and data reduction for the HERITAGE SDP data are given
in Meixner et al. (2010). For this paper, we use high
quality IRAS 100 m and MIPS 160
m observations instead
of the PACS observations which display large residual instrumental
signatures (expected to be eliminated with the full HERITAGE dataset).
We extracted the HERITAGE SDP region from the existing IRAS/IRIS
100
m (Miville-Deschênes & Lagache 2005) and MIPS 160
m (Meixner et al. 2006; Bernard et al. 2008) mosaics. We have used custom convolution
kernels created using the technique of Gordon et al. (2008) to convolve the images to a common resolution of
of the IRAS 100
m data. We also created a 2nd set of images (excluding the IRAS 100
m data) at the common resolution of the SPIRE 500
m and MIPS 160
m data of
38
.
Emission from Milky Way (MW) foreground cirrus clouds contributes to the far-IR emission seen in the LMC. We use the HI column density map created by integrating the MW velocities in the full HI cube (Staveley-Smith et al. 2003) to correct all the images for the MW infrared cirrus emission. The HI column densities were transformed to IR surface brightnesses using the model of the MW emission used by Bernard et al. (2008). Finally, any residual emission was removed by fitting a gradient across the SDP region using the regions in the strip beyond the IR edge of the LMC.
3 Results
For each point in the image, we determined the dust temperature by
fitting the observed far-IR SED to a modified black body of the form
![]() |
(1) |
The dust mass is computed from the measured 160

![]() |
(2) |
where








3.1 Dust temperatures and masses
The best fit dust temperature and mass values
were determined by for fits using only the pre-Herschel data
(IRAS 100 m and MIPS 160
m) and fits including
the Herschel SPIRE data (IRAS 100
m, MIPS 160
m, and SPIRE 250/350
m). Given the inclusion of the IRAS 100
m, we
used the
resolution images. The SPIRE 500
m data
are not included in these fits as it is usually systematically high
(see Sect. 3.2 and Meixner et al. 2010)
and including it in the fits only causes the residuals at the other
wavelengths to increase without significantly improving the fit. The
value of
used in the
fits was set to 1, 1.5, or 2 as this range encompasses realistic dust
grains (amorphous to crystalline grains) and is what has been used in
the past (Dunne et al. 2000). The dust temperature and
mass maps and fractional residual images for the
case
are shown in Fig. 1.
The differences between the best fit dust temperatures and masses (with and
without the SPIRE data) depends on the value of
used. For
,
the with SPIRE to without SPIRE temperature ratio is
and mass ratio is
.
For
,
the with/without SPIRE temperature ratio is
and mass
ratio is
.
For
,
the with/without
temperature ratio is
and mass ratio is
.
Thus, the inclusion of >200
data in the fits can
change the derived dust temperature by up to 8% and mass by up to
23% depending on the assumed value of
.
Prior to the Herschel observations, it was not possible to constrain
the best value of
given that there were only two infrared maps
of the LMC with
.
With the Herschel
observations, the behavior of the residuals as a function of
can be used to determine the optimal
value. Histograms of the
fractional residuals at different
values are shown in
Fig. 2.
A value of
clearly minimizes
the fractional residuals at all wavelengths with most of the pixels
having residuals of less 10% at all wavelengths except 500
m.
This result implies that either the characteristics of the dust grains
are intermediate between the two extremes or that a more complex dust
emission model including a distribution of dust temperatures and grain
sizes is needed (Paradis et al. 2009b; Draine et al. 2007). Assuming a
for the pre-Herschel fits (a common assumption) and using the
best fit
for the fits including the SPIRE data, we find
the with/without temperature ratio is
and mass ratio
is
.
This decrease in dust masses reduces the
magnitude of the ``FIR excess'' found by Bernard et al. (2008) for the
LMC. Roman-Duval et al. (2010) explore this issue in
detail for two specific LMC molecular clouds.
3.2 500
m excess
In the previous section, we have not included the 500 m observations in the analysis as it was seen not to improve the quality
of the fits and previous studies (Galametz et al. 2009; Galliano et al. 2005) have
observed submm fluxes in excess of that expected from fits to the
far-IR fluxes. At
resolution, the average fractional
500
m fit residual is 0.25, 0.10, and -0.05 for
values
of 2, 1.5, and 1 (Fig. 2). As a
is
strongly favored as it minimizes the residuals at all other wavelengths, we find a 500
m excess of approximately 10%. We find the same level of 500
m excess for fits done at both the
and
resolutions.
![]() |
Figure 2:
The histograms of the fractional residuals at different wavelengths are shown for
|
Open with DEXTER |
There are four possible origins of the 500 m
excess: 1) systematics due to our assumptions on our fitting;
2) a flux calibration error; 3) variations in the wavelength
dependent emissivity law (Agladze et al. 1996; Reach et al. 1995); and 4) very cold dust that mostly emits at
(Finkbeiner et al. 1999; Galliano et al. 2005). Whatever the the physical process responsible for the 500
m excess, the HERITAGE SDP SPIRE data of the
LMC allow us to probe the origin of the 500
m excess at high
spatial resolution in an external galaxy for the first time. We
tested the systematics of our fitting algorithm and searched for
correlations of the 500
m excess with different tracers of the
ISM conditions (dust temperature, dust mass, HI mass, and MIPS
24
m flux) in an attempt to determine the origin of the
500
m excess. The two strongest correlations are seen for MIPS
24
m flux (probing the ISM conditions for small grains) and the
total gas mass (probing the ISM density) are are shown in
Fig. 3.
To test 1), we performed Monte Carlo simulations where the
observations were simulated both with and without an excess at 500 m and with different
values. These simulations were fit with
varying
laws and realistic uncertainties. A 500
m excess was found in the simulations if it was part of the
simulation or if the fitting
was smaller than the simulation
.
Given that we empirically determine
from the <
data, our conclusion is that the excess we find is not a
result of our fitting method.
![]() |
Figure 3:
The 500 |
Open with DEXTER |
For 2), it is possible that there is a systematic 500 m flux
calibration error on the order of 10%. The official maximal possible
flux calibration error for SPIRE is 15% at any wavelength
(Griffin, et al. 2010). Given that we are including the SPIRE 250 and
350
m measurements in our fitting, the 500
m flux
calibration error would have to be relative to the other two SPIRE
bands and so is likely smaller than 15%. In addition, we would
expect to see no correlation between the excess and ISM condition
tracers, yet we see weak correlations
(Fig. 3).
For 3), a wavelength dependent increase in the dust emissivity law at
500 m on the order of 10% is possible (Paradis et al. 2009a). This
variation may be attributed to the dust grains amorphous/crystalline nature, size distribution, temperature, and material (Meny et al. 2007; Henning et al. 1995). For example, if the 500
m excess is due to small grains having a different
than large grains, we would expect the excess to be correlated with the MIPS 24
m emission (Reach et al. 1995). Yet the excess is weakly anti-correlated with MIPS 24
m flux (Fig. 3).
For 4), very cold dust that only emits at 500
m is
physically possible. The very coldest dust would necessarily be the
dust that is best shielded and, thus, we would then expect the excess
to be strongest in the highest density regions and regions with the
lowest radiation fields. Figure 3 gives a
conflicting answer as we see the largest excesses in the faintest
24
m regions (as expected) and least dense regions (not as
expected).
4 Conclusions
We investigate the importance of >200 m data in determining
dust temperatures and masses using new Herschel SPIRE observations of
the LMC (taken for the HERITAGE key project as part of the
Herschel science demonstration phase) combined with existing IRAS
100
m and Spitzer MIPS 160
m images. We fit the
observations with a model consisting of dust emitting as a single
temperature blackbody modified with an emissivity law proportional to
.
For fixed values of
,
fits using only the
100-160
m data give dust temperatures and masses that are on
average up to 8% and 23% different from fits using the same
and the 100-350
m data. The new SPIRE observations allowed us
to determine that
minimizes the residuals from 100 to
350
m. Using a
for the 100-160
m and a
for the 100-350
m fits results in an increase of
10% for the dust temperature and a decrease in the dust mass by 36%.
On average, there is a fractional excess at 500 m of
10%. The origin of the fractional excess is unlikely to be due
to our fitting algorithm or a flux calibration error, but it could be
due to either very cold dust that emits only
500
m or a
variation in the wavelength dependent change in the dust emissivity.
Planned HERITAGE observations of the LMC and SMC will allow for a more
detailed investigation of including >
data (mainly the
500
m excess) due to better quality PACS and SPIRE images
(optimized observations and cross-scans).
We acknowledge financial support from the NASA Herschel Science Center, J.P.L. contracts # 1381522 & 1381650. M.R. is supported by FONDECYT No1080335 and FONDAP No15010003. We thank the contributions and support from the European Space Agency (ESA), the PACS and SPIRE teams, the Herschel Science Center and the NASA Herschel Science Center (esp. A. Barbar and K. Xu) and the PACS and SPIRE data center at CEA-Saclay, without which none of this work would be possible. We thank the referee, Thomas Henning, for suggestions that improved the clarity of this paper.
References
- Agladze, N. I., Sievers, A. J., Jones, S. A., Burlitch, J. M., & Beckwith, S. V. W. 1996, ApJ, 462, 1026 Bernard, J.-P., Reach, W. T., Paradis, D., et al. 2008, AJ, 136, 919 [NASA ADS] [CrossRef] [Google Scholar]
- Dale, D. A., Bendo, G. J., Engelbracht, C. W., et al. 2005, ApJ, 633, 857 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, 663, 866 [NASA ADS] [CrossRef] [Google Scholar]
- Dunne, L., Eales, S., Edmunds, M., et al. 2000, MNRAS, 315, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 [NASA ADS] [CrossRef] [Google Scholar]
- Fukui, Y., Kawamura, A., Minamidani, T., et al. 2008, ApJS, 178, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Galametz, M., Madden, S., Galliano, F., et al. 2009, A&A, 508, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Galliano, F., Madden, S. C., Jones, A. P., Wilson, C. D., & Bernard, J.-P. 2005, A&A, 434, 867 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Galliano, F., Dwek, E., & Chanial, P. 2008, ApJ, 672, 214 [NASA ADS] [CrossRef] [Google Scholar]
- Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 594, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gordon, K. D., Engelbracht, C. W., Rieke, G. H., et al. 2008, ApJ, 682, 336 [NASA ADS] [CrossRef] [Google Scholar]
- Griffin, M. J., et al. 2010, A&A, 518, L3 [Google Scholar]
- Henning, T., Michel, B., & Stognienko, R. 1995, Planet. Space Sci., 43, 1333 [NASA ADS] [CrossRef] [Google Scholar]
- Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, S., Staveley-Smith, L., Dopita, M. A., et al. 2003, ApJS, 148, 473 [NASA ADS] [CrossRef] [Google Scholar]
- Laor, A., & Draine, B. T. 1993, ApJ, 402, 441 [NASA ADS] [CrossRef] [Google Scholar]
- Meixner, M., Gordon, K. D., Indebetouw, R., et al. 2006, AJ, 132, 2268 [NASA ADS] [CrossRef] [Google Scholar]
- Meixner, M., et al. 2010, A&A, 518, L71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meny, C., Gromov, V., Boudet, N., et al. 2007, A&A, 468, 171 [Google Scholar]
- Miville-Deschênes, M.-A., & Lagache, G. 2005, ApJS, 157, 302 [NASA ADS] [CrossRef] [Google Scholar]
- Paradis, D., Bernard, J., & Mény, C. 2009a, A&A, 506, 745 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Paradis, D., Reach, W. T., Bernard, J., et al. 2009b, AJ, 138, 196 [NASA ADS] [CrossRef] [Google Scholar]
- Reach, W. T., Dwek, E., Fixsen, D. J., et al. 1995, ApJ, 451, 188 [NASA ADS] [CrossRef] [Google Scholar]
- Roman-Duval, J., et al. 2010, A&A, 518, L74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Russell, S. C., & Dopita, M. A. 1992, ApJ, 384, 508 [NASA ADS] [CrossRef] [Google Scholar]
- Schaefer, B. E. 2008, AJ, 135, 112 [NASA ADS] [CrossRef] [Google Scholar]
- Stansberry, J. A., Gordon, K. D., Bhattacharya, B., et al. 2007, PASP, 119, 1038 [NASA ADS] [CrossRef] [Google Scholar]
- Staveley-Smith, L., Kim, S., Calabretta, M. R., Haynes, R. F., & Kesteven, M. J. 2003, MNRAS, 339, 87 [NASA ADS] [CrossRef] [Google Scholar]
- Swinyard, B. M., et al. 2010, A&A, 518, L4 [Google Scholar]
- Willmer, C. N. A., Rieke, G. H., Le Floc'h, E., et al. 2009, AJ, 138, 146 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... 200 micron
- Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
All Figures
![]() |
Figure 1:
The best fit (assuming an emissivity law with
|
Open with DEXTER | |
In the text |
![]() |
Figure 2:
The histograms of the fractional residuals at different wavelengths are shown for
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
The 500 |
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.