A&A 401, 873-893 (2003)
DOI: 10.1051/0004-6361:20021837
K. Wilke - M. Stickel - M. Haas - U. Herbstmeier - U. Klaas - D. Lemke
Max-Planck-Institut für Astronomie, Königstuhl, 69117 Heidelberg, Germany
Received 30 September 2002 / Accepted 10 December 2002
Abstract
The ISOPHOT experiment onboard the
ISO satellite
generated a complete view of the Small Magellanic Cloud (SMC) at
170 m with 1.5 arcmin resolution.
The map is analysed
using an automated photometry program enabling accurate photometric
characterization of the far infrared (FIR) emitting regions.
An integrated FIR luminosity of
is obtained, leading to a star
formation rate of
.
With an
average dust temperature of
,
the total dust mass follows
to
.
In this paper, the sources detected at 170
m are compared with those obtainable from
the IRAS satellite data. For this purpose, the 12
m, 25
m, 60
m, and 100
m IRAS
high resolution (HiRes) maps of the SMC are re-examined using the same
method. In contrast to former studies,
this provides an all-band ISO/IRAS source catalog which is no
longer based on eyeball classification,
but relies on an algorithm which is capable of
automated, repeatable photometry, even for irregular sources.
In the mid infrared IRAS bands
numerous bright FIR emitting regions in the SMC are detected
and classified: 73 sources are found at 12
m, 135 at 25
m
(most of them with
).
All three FIR bands at 170
m, 100
m, and 60
m reproduce the overall
morphological structure of the SMC similarly well, in contrast to the
12
m and 25
m maps which only contain a limited number of
extended sources and do not trace the main body of the SMC.
243 sources are detected in the
ISO 170
m map, 155 of them with
.
Comparable numbers are found for the two FIR IRAS maps at 60
m
(384) and 100
m (338) with fluxes up to
.
70 of the 243
170
m sources are assigned a general SED type ("cold'', "warm'', i.e.,
,
)
for the first time.
A comparison with earlier IRAS results suggests that
many source flux densities in those studies have been under- or
overestimated
because of non-standardized fitting methods. Many
sources with flux densities up to 40 Jy listed in former catalogs
cannot be identified in our
data.
Key words: galaxies: magellanic clouds - ISM: general - ISM: dust, extinction
As the nearest dwarf irregular galaxies, the SMC and the LMC serve as excellent laboratories for studying their detailed compositions since all objects in them can be assigned the same distance (50...60 kpc, see Westerlund 1997). At the same time, radio and optical data can well be obtained down to parsec and sub-parsec scales.
Both galaxies have been found to be low in dust and heavy-element abundances in comparison to the solar neighbourhood, however, high star formation rates lead to significant infrared emission from warm and cool dust. Numerous studies focus on the structural peculiarities, e.g., the extremely rich and complex HI morphology with numerous apparently expanding shells, filaments, and arcs (for results see, e.g., Martin et al. 1989, and Staveley-Smith 1997, and references therein). The SMC also hosts several giant HII regions, the referring catalog of which was produced by Davies et al. (1976).
The first infrared study by Schwering & Israel 1989 (SI89 hereafter) was based on IRAS data. They found numerous new objects not listed in the IRAS published catalogs (PSC and SSS, see Beichmann et al. 1985; Helou & Walker 1985, and references therein) and compared the infrared results to catalogs obtained in different wavelength bands. The near-infrared (NIR) DENIS-based point source catalog (PSC) towards both Magellanic Clouds has recently been published by Cioni et al. (2000).
A comprehensive radio continuum study (e.g., Haynes et al. 1991; Filipovic et al. 1997, 1998a,b) addressed the question of the possible existence of magnetic fields in the SMC, the distribution of radio sources, and the comparison with other surveys (e.g., X-ray and IR data).
The LMC and SMC also became part of an ESO-SEST key programme which performed CO survey observations in both galaxies and focused on the molecular gas content and the physical properties of the molecular clouds (see, e.g., Israel et al. 1993; Rubio et al. 1993, and Lequeux et al. 1994). An X-ray source catalog can be found in Haberl et al. (2000). Comprehensive overviews of both Magellanic Clouds as well as a list of further references can be found in Westerlund (1990) and Westerlund (1997).
Since the IRAS mission the question has been raised how much (if any)
cold dust exists in certain objects (e.g., star forming galaxies)
which escapes an detection due to its low
temperatures. The latter would cause the FIR SEDs to keep on
rising significantly beyond the IRAS 100 m band. In order to
accurately determine the temperature of the cold dust in the SMC,
we therefore made use of ISO's long wavelength capability and obtained
a complete 170
m map.
This is the first of two papers in which the FIR properties of the SMC
are discussed. In this paper, we put special
emphasis on the overall SMC morphology at 170
m, the
statistical results of the ISO SMC observations,
the comparison with the IRAS data, and former
studies based on IRAS source catalogs. In a subsequent paper we will
discuss the quantitative and detailed
properties of the warm and colder dust, ranging from color
temperatures and integrated FIR luminosities over dust properties of single
sources to the properties of the interstellar radiation field and the
star formation rates.
The observations with the ISOPHOT (Lemke et al. 1996)
photometer onboard the ISO satellite (Kessler et al. 1996)
were performed in January 1998 in raster mode (AOT PHT22) with the
C200 detector, a
pixel array of stressed Ge:Ga with a pixel
size of
,
in conjunction with the C_160 broad band filter
(reference wavelength 170
,
equivalent width
89
m). Due to the rather large area of the SMC on the
sky, the ISOPHOT observations had to be split into a mosaic of nine
separate parts, each of
which was accompanied by two observations of the ISOPHOT Fine
Calibration Source (FCS). Adjacent parts of the whole map were
designed to be slightly overlapping, while the raster step size of each part was a
full detector size without any overlap.
Since there was little redundancy in the data, cosmic ray hits could mimic compact sources in the map. Therefore, instead of the standard ramp slopes from first-order polynomial fits, the pairwise differences of consecutive ramp readouts were used to derive the detector signals. This allowed a larger distribution to be analysed, leading to considerably more robust results.
To get rid of pairwise differences affected by cosmic ray hits, the robust outlier-insensitive myriad estimator (Kalluri & Arce 1998) was computed and 20% of the most deviant signals as measured by the absolute deviation were cut off. This outlier removal is similar to a median absolute deviation trimming, but instead of the initial median, the sample myriad is used to determine the outliers. The sample myriad value in turn is a robust estimator of the mode (most common value) of a distribution but does not require binning of the actual data set, and is easily computed by simply minimizing a particular cost function with a tuning constant set to a small value (for details see Kalluri & Arce 1998). After rejecting the outliers, the trimmed set of pairwise differences was linearly interpolated to a ten times finer grid and the value giving the minimum value for the myriad cost function accepted as the final signal for each raster point. This interpolation scheme was used as an approximation to a full numerical minimization of the myriad cost function.
![]() |
Figure 1:
False color representation of the 170 ![]() |
Open with DEXTER |
The derived detector signals at each raster position were subsequently
corrected for signal dependence on ramp integration times to be
consistent with calibration observations (Laureijs et al. 2000),
dark-current subtracted, and finally flux calibrated
with PIA Version 9.1/Cal G Version 6.0 (Gabriel et al. 1997).
For the conversion to an absolute flux level,
the observations of the ISOPHOT Fine Calibration Source (FCS) obtained at
the beginning and end of each raster in each filter were used.
The flux calibrated data streams of the detector pixels still showed differences in the overall levels of up to 20%, mostly due to inappropriately corrected pixel-to-pixel sensitivities (flat field). If not removed, these varying brightness levels would lead to striping and chessboard-like patterns in the maps. Robust morphological filtering techniques (Sternberg 1986) were used to extract the overall level of the four data streams, which were then brought to a common mean level, thereby giving the relative pixel scaling factors.
Eventually, a complete map of the whole SMC was produced from the
flatfielded flux-calibrated data streams of all pixels by using the
Drizzle Mapping Method (Hook & Fruchter 2002) within
IRAF,
which took into account the pixel sizes and inter-pixel distances of
the C200 detector and the detector roll angle. The pixel size used
for the final full map was
,
the smallest possible size
which did not produce uncovered holes. This final map was restored using
the modified Richardson-Lucy-Algorithm (Hook et al. 1994) with an
additional subsampling of two pixels and a point spread
function approximated by a Gaussian with a FWHM of
,
giving a final restored map with
pixel size. This map is shown in
Fig. 1.
To allow for a direct comparison with the shorter wavelength IRAS HiRes
maps (see below), restored ISOPHOT (sub-)maps with the same center and size as the
IRAS sub-fields were created as well as with a pixel size of
.
This is
twice the IRAS HiRes pixel size, so that IRAS maps exactly aligned
with the ISO sub-fields could be constructed by a simple
pixel
block-averaging.
Although a more quantitative analysis of the FIR properties of the SMC
will be performed in subsequent paper, Table 1 lists some preliminary
global quantities. Remarkable are the numbers found for the
global star formation rate, and the gas-to-dust ratio resulting from the additional cold dust
component entering the total dust mass.
F40-220 |
![]() |
F1-1000 |
![]() |
L1-1000 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Figure 2:
Orientation of the
![]() |
Open with DEXTER |
Although source catalogs for this wavelength range were
published in the past already (for an extensive comparison of IRAS
SMC source catalogs with other wavelength bands see, e.g., SI89),
we decided to re-analyse the IRAS
high resolution (HiRes hereafter) data in order to treat ISO results and
IRAS data in an identical and reproducible way.
Since the area in the sky covered by the SMC is comparably large
(
in the optical) and IRAS HiRes data can be
requested with maximum field sizes
of
only, we split up the
SMC field into 9 single
fields with the central one covering most of
the central body of the SMC in the FIR (for the orientation of the fields in the sky see
Fig. 2). This ensured that the borderlines between two
adjacent fields (though covered
by a small overlap) would not be located in the brightest parts of the
SMC.
The central coordinates for the requested
fields (equinox
B1950.0) are given in Table 2. Pixel sizes were 15''for all fields.
Field | RA | DEC |
---|---|---|
0 | 13.320000 | -72.937271 |
1 | 20.093601 | -72.937271 |
2 | 20.093601 | -70.949771 |
3 | 13.320000 | -70.949771 |
4 | 6.5463982 | -70.949771 |
5 | 6.5463982 | -72.937271 |
6 | 6.5493982 | -74.924771 |
7 | 13.320000 | -74.924771 |
8 | 20.093601 | -74.924771 |
The data reduction for the IRAS scans at IPAC comprises several steps:
first, the calibrated,
reconstructed detector data is deglitched (which removes spurious
non-source-like signals originating from radiation impacts on the
detector) and destriped (which corrects for different detector
responsivities during different scans, i.e., additive offsets of
certain strips). The zodiacal emission model was then subtracted from
all data scans individually. For the HiRes data fields a maximum
correlation method (MCM) for the reconstruction of the original image is
applied to the single scans which not only iteratively builds a
reliable model of the sky brightness, but also
enhances the resolution to 15'' pixel size. For a more detailed description of
IRAS data reduction routines see Assendorp et al. (1995), Bontekoe et al. (1994),
and Aumann et al. (1990).
All resulting IRAS maps are calibrated in
MJy/sr, two of them (central 60 m and 100
m) are shown in Fig. 3.
The visual inspection of all nine fields yielded the result that
only in five of them objects belonging to the SMC are located. These fields are
oriented in a cross-like pattern in the sky with numbers 0, 1, 3, 5, and 7 (see
Fig. 2). All other fields were neglected for the
subsequent analysis.
![]() |
Figure 3:
IRAS HiRes maps of the SMC. Left: 60 ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Initially, the program defines a rank list of possible detections, to which one out of six basic morphologies, ranging from Gaussian profiles to rings with different major-to-minor axis ratios, is assigned (for an extensive discussion on the construction of the source candidate list see Thilker et al. 2000). In a next step, so-called "footprint'' areas are constructed by "allowing'' the program to allocate pixel areas of the SMC input image according to the morphology assigned to the source (Fig. 4a) which may contain pixels which are not bright enough to remain in the final boundary of the source after the end of the growth procedure. For this reason, "seed'' regions are constructed by rejecting pixels falling below a certain median surface brightness limit of the initial "footprint'' region (Fig. 4b). In a third step, iterative growth starts: pixels are considered down to a limit equal to the outermost isophote, this limit is reduced by 0.02 dex in every iteration until a certain lower limit is reached (Fig. 4c). The program offers the option of making arbitrary selections for this limit where the surface brightness profile has become "sufficiently'' flat (see Fig. 4d).
The
resolution of the IRAS HiRes maps is highly asymmetric, especially for the
MIR range at 12 and 25 micron. Basically, this is the result of the
rectangular detector mask shapes and the geometry of the scans
covering the sky. Consequently, point sources appear
elongated with the narrow dimension in the scan direction and the
larger dimension determined by the cross-scan-size of the detector.
Since we use IRAS HiRes data,
it is very difficult to determine
an effective resolution, however, the resolution of unenhanced coadded IRAS images
of approximately
,
,
,
and
for the 12
m, 25
m, 60
m, and 100
m data provides us with an
impression of which resolution changes may occur in our maps when
proceeding towards longer wavelengths.
After the source catalogs were generated for every wavelength and the
classification of the detected sources (the classification scheme is
explained in detail in the next section), they were correlated with
each other, i.e., we tried to identify sources in different catalogs
within a certain correlation radius. We decided to choose a value of
due to the high spatial resolution of the ISO and IRAS
data, though the average values for
that radius were slightly larger in former studies (e.g., see Filipovic et
al. 1998b who used
for the comparison
between IRAS and radio data). This ensured the detection of all relevant source
pairs in different wavelength bands and avoided a too large number of
multiple correlations at the same time, a condition which constitutes the
corresponding upper limit for
.
As is clearly visible from
the resulting tables which are presented in Appendices A-E, we never
encountered more than 5 cross-identifications with more
than one source in one band.
![]() |
Figure 4:
Illustration of the various stages of the HIIphot procedure
for a special area. a)
IRAS 100 ![]() |
Open with DEXTER |
The general classification scheme is the following:
a source is labelled "C'' (cold) when the color temperatures
(assuming a modified black body with an
emissivity of )
derived from both, the
170
m/100
m and the 100
m/60
m ratios, are
,
it is
classified "W'' (warm) when
is obtained in both cases.
"(C)'' denotes sources where a classification as "C''
remains uncertain due to missing 60
m and/or 100
m data
points. A source is labelled "Q'' (questionable) when the precise shape of the SED cannot be derived due to
source confusion. There are a few cases where the SED seems to
"drop'' between 60
m and 170
m due to a lower
100
m flux density value. This behaviour has to be attributed to different PSF sizes in the
different wavelength bands, leading to correlations of
physically independent sources in the three wavelength
bands.
The first letter of the identification label
is followed by a down- or uparrow (
or
)
in those cases where
F170/F100-ratios <1 and >1,
respectively. Identical
170
m and 100
m flux densities are indicated by a =.
Note that for a "cold'' BB-spectrum, i.e.,
,
the
F170/F100-ratio needs not necessarily exceed 1.0.
A cross-identification with the source catalogs given by SI89 leads to the next classification letter, which may be "(C)'' or "(W)'', hinting at their cold and warm sources. Note that this type of coding could in principle lead to an ambiguity in the meaning of the single-letter label "(C)'': this could either mean that a SI89 identification is referred to without any identification in our catalog or be the label for a probably cold source which is only found in our source study. In fact, all sources classified as "(C)'' are assigned this source type in our catalog only, so this ambiguity is avoided.
ISO 170
![]() |
|||||||||
n5 |
![]() |
![]() |
F5 |
![]() |
n4 | n3 | n2 | n1 | identifications |
1 | 00 41 40 | -73 02 16 | 2.2 | 0.4 | 40 42 | 20 24 | 9 | C IRAS 00395-7317 | |
2 | 00 42 06 | -73 08 34 | 0.3 | 0.0 | |||||
3 | 00 42 06 | -72 59 58 | 4.1 | 0.4 | 50 | 32 | 14 | C![]() |
|
4 | 00 42 58 | -72 59 36 | 5.5 | 0.5 | 58 | 39 | 17 | Q(C) IRAS 00413-7316 | |
5 | 00 43 25 | -73 02 31 | 3.2 | 0.4 | 59 | 41 | 19 | Q(W) | |
6 | 00 44 37 | -72 58 06 | 0.7 | 0.3 | 71 | 20 | |||
7 | 00 45 18 | -73 05 23 | 59.8 | 1.6 | 76 | 56 | 24 | 15 | C![]() |
8 | 00 45 21 | -73 16 16 | 228.8 | 2.7 | 74 | 54 | 22 | 16 | C![]() |
9 | 00 45 28 | -73 22 07 | 259.5 | 0.0 | 75 | 55 | 23 | 17 | C![]() |
10 | 00 46 32 | -73 05 59 | 105.8 | 1.8 | 82 | 64 | 29 | 20 | C
![]() |
11 | 00 46 34 | -73 21 23 | 158.0 | 2.4 | 79 | 28 | 18 | C![]() |
|
12 | 00 46 35 | -73 15 28 | 84.6 | 1.7 | 81 | 65 | C![]() |
||
13 | 00 46 45 | -73 30 51 | 46.1 | 0.0 | 84 | 68 | 30 | C
![]() |
|
14 | 00 46 46 | -73 09 37 | 93.9 | 1.7 | 83 | (C)![]() |
|||
15 | 00 47 37 | -73 05 26 | 329.3 | 3.3 | 89 | C![]() |
|||
16 | 00 47 52 | -73 15 16 | 453.4 | 3.6 | 91 | 79 | 34 | 23 | C![]() |
17 | 00 47 55 | -73 35 21 | 1.7 | 0.0 | 90 | 81 | |||
18 | 00 47 56 | -73 23 47 | 13.1 | 0.9 | 75 84 | (C)(W) IRAS 00462-7339 | |||
19 | 00 47 58 | -72 23 13 | 4.3 | 0.4 | 92 | 83 | C![]() |
||
20 | 00 48 13 | -73 08 01 | 186.2 | 2.4 | 77 | 32 | C(C) X(HFP2/409) | ||
21 | 00 48 34 | -72 57 51 | 6.6 | 0.8 | 98 | 89 | 38 | Q(W) IRAS 00467-7314 | |
22 | 00 48 37 | -72 26 54 | 0.8 | 0.3 | 97 | 88 | |||
23 | 00 48 39 | -73 19 52 | 32.6 | 1.3 | 96 | C![]() |
|||
24 | 00 48 52 | -72 48 38 | 65.2 | 1.8 | 99 | 93 | C![]() |
||
25 | 00 48 52 | -73 08 33 | 269.4 | 3.0 | 100 | 94 | 26 | C![]() |
|
26 | 00 48 55 | -72 54 54 | 12.0 | 1.0 | 101 | 96 | 40 | Q Em*(MA93 262) | |
27 | 00 49 23 | -73 34 27 | 2.6 | 0.5 | 103 | 97 | C![]() |
||
28 | 00 49 27 | -73 26 49 | 34.3 | 1.2 | 110 | 101 | 43 | C![]() |
|
29 | 00 49 42 | -73 23 45 | 9.3 | 0.8 | 46 | (C) X(HFP2/468) Em*(LIN 139) | |||
30 | 00 49 58 | -73 12 00 | 26.5 | 1.2 | 106 | (C) X(SHP2/SMC30) Em*(MA93 335) | |||
... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
Our sources detected by HIIphot in all five bands were cross-correlated with the SIMBAD database in order to identify well-known HII regions, X-ray sources, etc. The SIMBAD database was searched using a correlation radius of 60'' for all catalogs (which is different from the coordinate correlation radius of 90'' used in our IR study). From a large list of possible identifications the most important catalogs were included:
![]() |
Figure 5:
12 ![]() |
Open with DEXTER |
The question of the underlying physics in the SMC interstellar medium and the nature of single sources will be extensively discussed in Paper II where we will also put our SMC results in the context of already existing complementary observations and catalogs, e.g., the FUV spectroscopic explorer survey of interstellar molecular hydrogen in the SMC by Tumlinson et al. (2002), HI data as analysed in Stanimirovic et al. (2000), or CO data (see, e.g., the results from the ESO-SEST key programme in the SMC given in Israel et al. 1993 and subsequent papers).
The cutoff in surface brightness
was adjusted to 1.4 MJy/sr which caused the program to limit the source
growth to the visually determined size.
The MIR emission in Fig. 5 (upper plot)
is dominated by the presence of numerous discrete sources while
nearly no diffuse emission is detected. The discrete
sources are mainly located in two different regions: most of them
cluster in a ring-like structure in the SW part of the SMC, and some
of them are located in the NE part of the SMC bar. Most
detections are found in the central IRAS field 0 (upper plot in
Fig. 5).
To summarize, at 12 m the SMC appears as a patchy pattern of
single discrete sources which are not a tracer of the
morphology of this galaxy.
As given in Table 4,
25 of our 73 sources listed in Appendix A were identified as
IRAS PSC or FSC sources using the SIMBAD database within a correlation radius of 60'', 14 were classified
as HII regions already in former studies. Finally, 38 sources
fit into the category "emission line
star/emission object''. The SIMBAD results for the identification of
carbon stars are omitted in Appendix A and in Table 4, since nearly all of our sources
(65) were associated to such a star within the correlation radius. 18
sources could be classified as type "C'' (cold) according to the
classification criteria for the 170 m catalog given above. 14 of the sources show a
170
m/100
m flux ratio greater than 1.0. Of the sources
which were not found in the 170
m ISO map but could be classified
using the 60
m and 100
m data, 2 have flux ratios
,
other two ones yield
.
![]() |
Figure 6:
25 ![]() |
Open with DEXTER |
The source counts resemble numbers and flux distribution of the 12 m
map: as can be seen from Fig. 13, the bulk of the
flux densities lies in the interval
,
although there is now a larger percentage of
sources with even higher fluxes (17% of the 135 sources). Our
five marked sources - again among the brightest ones -
in Fig. 6 with numbers 34, 80, 87,
114, 127 were assigned fluxes of 10.0 Jy, 29.0 Jy, 12.5 Jy, 6.5 Jy, and 19.2 Jy.
Four cases are found where two 100
m sources lie within the
correlation radius of 90'' around a 25
m detection.
![]() |
Figure 7:
60 ![]() ![]() |
Open with DEXTER |
The number of extracted sources
listed in Appendix C is 384, 124 of which reach flux densities of 1.0 Jy. Only these bright sources were used for further
investigations. The number
of sources identified with IRAS PSC/FSC sources or HII regions in former
studies increases considerably: more than half of the sources (71)
with flux densities
1.0 Jy are listed in the IRAS catalogs, 33
are classified as HII regions in the DEM sample. Finally, 84 of
them fall into the general category "emission line star/emission object''. 95 of the
brighter sources were classified as cold sources with
K ("C'' in our scheme), only 6 warm sources ("W'')
were found. This means that nearly all relevant source detections
above the flux limit are of cold origin. Even when compared to the 25
m
results, this number has nearly doubled (see Table 11).
For 9 sources, a precise classification was not possible
due to missing wavelength bands, but they are assumed to be
predominantly of type (C), due to their 170
m flux. Confusion with
other sources in complementary wavelength bands affected the
classification of 20 sources.
Of 101 sources without 170
m identification, 75 exhibited
,
26 were found with
.
The five landmark sources with numbers 79, 221, 271, 337, and 369
are assigned fluxes of 147.6 Jy, 283.2 Jy, 48.7 Jy, 79.2 Jy, and 56.6 Jy.
At 100 m, the SMC extends even further out than
at 60
m, the contours at the lower flux density
levels suggest an overall homogeneous bar structure with
significantly brighter regions in the NE and SW part
and a few bright sources in between. Structures in the
bridge region (Fig. 8, lower left plot) become
more extended than at 60
m (Fig. 7).
![]() |
Figure 8:
100 ![]() ![]() |
Open with DEXTER |
167 sources have flux density values 1.0 Jy and are used for further analysis. The remaining sources have
the same properties as in the 170
m map (cf. following section): most of them are located
in the outer regions of the SMC, far away from the more luminous HII
regions along the major axis of the bar. Nearly all of them have low
flux densities, too. According to Table 4, 90 of the 167
luminous sources can be classified as "C'', the classification of 18
detections remains unclear due to source confusion between
different wavelength bands.
Of the sources of type "C'', 70 are found with
, 16 with
.
These numbers are very similar to those found
for the 60
m catalog, although the percentage of cold sources
shows a small decrease compared to the 60
m case.
Only 7 sources exhibit a spectrum with
K, according to the classification criteria for our
object class "W''. Among the sources which were not detected in the
170
m band, the majority (57) exhibits flux ratios (
)
above 1.0, some of them (22) also showed decreasing
fluxes with increasing wavelengths (
). The total number of objects as
well as the distribution among the different object classes makes the
100
m data directly comparable to the 60
m and the 170
m
case. 63 of the 167 sources with
Jy are identified as
IRAS sources from the PSC or FSC catalog, 22 of them are HII regions, and 70 of them
fall into the category "emission line star/emission object''.
The final map shown in Fig. 1 reveals a wealth of structure, not only along the bar of the SMC, but also in the direction of the Magellanic bridge, connecting the SMC to the LMC. The emission of the brightest HII regions (colored blue/white in Fig. 1) mainly originates from a ring-like structure in the southwestern part of the SMC, but also from isolated sources in the northeastern part of the bar. Three further bright HII regions are obviously located in the wing region. In the outer regions of the bar, numerous fainter sources (orange/red/yellow color) were found, which contribute significantly to the total FIR emission of the SMC.
![]() |
Figure 9:
ISO 170 ![]() ![]() |
Open with DEXTER |
As given in Appendix E, a total number of 243 sources at 170 m is detected
of which 155 exhibit flux
density values above 2.0 Jy. The remaining 88 weak sources were not
taken into account for further studies, since most of them are
basically cold clouds in the outer regions of the bar which do not
constitute the morphological structure (but nevertheless contribute to
the integrated FIR emission).
The five sources labelled 16, 103, 135,
195, and 232 in Fig. 9 were assigned a growth luminosity of
453 Jy, 366 Jy, 198 Jy, 327 Jy, and 22 Jy, respectively (see table in
Appendix E).
51 of the 155 bright sources were identified
as IRAS PSC or FSC sources in the SIMBAD database, 99 were classified as
cold sources (type "C''),
7 as warm (type "W''), and 19 remained unclassified ("Q'') due to
obvious miscorrelation with other sources in the complementary
wavelength bands. 30 source SEDs could not be classified due to
missing flux density values at 100 m or 60
m, but are probably
cold due to their high 170
m value. 18 objects could be identified
as HII regions from the DEM sample, 61
contributed to the "emission line star/emission object''
category. For 74 of the sources the
ratios are
1.0, 16 have
values <1.0.
To summarize, the source extraction at 170
m yields numerous
new and formerly unknown cold FIR emitting sources which are not
identified in other catalogs. Of the 155 brighter sources, 70 were
assigned a general SED type ("W'', "C'', "(C)'').
In addition, many of the sources classified as warm ("W'') in former studies were re-classified "C'' by our study.
This puts special emphasis on the importance of the 170
m flux value
in an SED of cold (<30 K) objects for the determination of the
blackbody temperature used for the dust emission models.
The basic properties of the
170 m luminosity distribution given by the histogram
in Fig.10 are found to be the same as in the 60
m
and 100
m cases: a large
number of sources in the lowest flux bins (
)
dominates the histogram distribution. On the other hand,
single sources with very high luminosities are detected,
ranging up to the 400 Jy level.
Compared to the
100
m data, the decrease in the total number of
detections is of course caused by the considerably larger
PSF FWHM at 170
m, making the separation of small
adjacent sources difficult.
Some 100
m sources seem to be missing in the 170
m
observations.
This is caused by the fact
that the ISO 170
m
map does not cover the SW extension of the SMC with a few ten sources.
In general, the number of sources which are also detected at 170 m
increases considerably with increasing wavelength. E.g., while
45% to 47% percent of all sources found at 25
m,
60
m, and 100
m are detected at 170
m, this only
holds for 29% of the 12
m sources.
Different dust heating mechanisms play a dominant role in
producing the different appearance:
While the 12
m and 25
m emission are dominated by ultrasmall
grains which are subject to single-photon heating (Draine &
Li 2001), the 100
m through 170
m emission is
mostly produced by the emission of large grains of equilibrium
temperature. We also see different types of sources in the
MIR (12
m, 25
m) and FIR (60
m, 100
m, 170
m) bands.
In the MIR, stars account for more than half of the
entries already in the IRAS point source catalog
(Beichmann 1987),
with most of the emission coming from AGB stars with
SEDs similar to those found in the papers by Groenewegen et al. (1999)
or Fujii et al. (2002). Large numbers
of these stars have been detected in the SMC already with their spectral
properties being in good agreement with those found in our Galaxy
(Zijlstra et al. 1996). These SEDs peak
at 10
m to 30
m and there is no significant emission
beyond 60
m.
The FIR emission is most likely dominated by HII
regions with an SED peaking between 60
m and 100
m
(Ghosh 2000) and molecular clouds.
Secondly, the two FIR IRAS maps (60 m, 100
m) and the 170
m ISO map
are closely related with respect to cross-identifications: of all ISO
sources, more than 60% are identified at 60
m and 100
m as
well, in addition, more than 45% of all 60
m and 100
m
sources are found at 170
m, too (see Table 4).
![]() |
Figure 10:
Source distribution histogram for the 170 ![]() |
Open with DEXTER |
![]() |
Figure 11:
12 ![]() ![]() ![]() |
Open with DEXTER |
Figure 11 compares the flux distribution of our
reliably detected 12 m sources (upper plot) with those given in the
12
m catalog in SI89.
Both data sets show a similar luminosity
distribution with the bulk of the sources below 0.2 Jy and a few discrete sources with
considerably larger values.
Obviously, our algorithm detects less sources at very low
intensities. Most of the faint
sources found by SI89 were assigned a highly uncertain flux value
which, in most cases, is very close to or identical with the detection
limit in their IRAS maps.
![]() |
Figure 12:
Comparison of the 12 ![]() |
Open with DEXTER |
Addressing the comparability of the SI89 results with ours, two questions have to be answered: i) how many sources are commonly detected in different wavelength bands and listed in the corresponding source catalogs, and ii) whether sources with different properties (luminosity, location in the SMC, background level/structure in their vicinity, etc.) are assigned different fluxes by the competing methods. To answer these questions, we directly compared our HIIphot IRAS source catalogs to those published in SI89 (see Fig. 12).
The faint detections (close to the zero point in Fig. 12)
which are commonly identified in both 12 m catalogs
tend to cluster roughly around the equality
line, i.e., on average, our study and the SI89 method seem to assign identical source
fluxes. However, statistics are
poor due to the small total number of sources.
A number of SI89 12
m objects is not detected by our
algorithm. Faint sources from our catalog (<0.5 Jy) tend to have
slightly lower fluxes
than the SI89 ones
and therefore lie below the line
.
Looking at the few brighter sources only, the situation seems
different:
the flux ratios tend to lie closer to one, and
in addition, more
sources tend to lie above the line
,
i.e., our growth provides us with higher final fluxes
than found in the SI89 catalog.
![]() |
Figure 13:
Distribution histogram for the 25 ![]() ![]() ![]() |
Open with DEXTER |
For the 25 m source list, the situation is similar: in general, our luminosity distribution and source counts
(Fig. 13, upper plot) resemble
very much those of SI89 (lower plot). The
only difference is that we detect slightly less sources in the
luminosity bin 1.0 Jy
2.0 Jy.
For large luminosties (
6 Jy),
statistics become poor in both cases, so a reliable comparison of both
data sets in this
flux regime is not possible.
Regarding the comparison of the fluxes assigned by HIIphot and the
SI89 study (Fig. 14),
it is obvious that many more sources of low luminosity than
in the 12 m case
are commonly identified in the two different 25
m catalogs,
providing us with better statistics.
The 25
m results are very similar to the 12
m
ones: for fainter sources (<1.5 Jy), our study and the SI89 results
are in very good agreement,
brighter sources are again assigned slightly higher values in our study.
![]() |
Figure 14:
Comparison of the 25 ![]() ![]() |
Open with DEXTER |
![]() |
Figure 15:
Distribution histogram for the 60 ![]() ![]() |
Open with DEXTER |
Addressing number counts and cross-correlation with other catalogs,
major changes occur at 60 m:
the overall number of reliable detections rises considerably to
384, which is the largest
value of all catalogs/wavelength bands considered here.
But while the analysis performed in SI89 yields only
68 sources with fluxes up to 2.5 Jy, our list in
Appendix C
contains 291 discrete sources at the same
flux level. Evidently, our detection algorithm is well capable of
detecting low-flux sources even in crowded regions of the main body
of the SMC.
According to the data presented in SI89, a broad source distribution
extends up to fluxes of
20 Jy. Beyond that, only single sources are detected in
each flux bin (see Fig. 15, lower plot). In this flux
regimes, statistics become poor, so a comparison between the two
methods does not make any sense anymore.
In general, the source flux comparison of the 60 m catalog with the SI89 data set shown in
Fig. 16 (upper and lower plot) yields results which
slightly differ from those found for the
25
m data (Fig. 14):
faint sources (<10 Jy) now lie distincly below the line
and exhibit a large scatter for all flux values up to 10 Jy.
Similar to the 12
m and 25
m case, brighter sources with
fluxes up to
160 Jy tend to lie above the equality line. But
since statistics are better now for the compared 60
m catalogs, it
is evident that the HIIphot fluxes do not approach the SI89 ones, even
not for the highest flux regime.
This may be a result of the
automated growth procedure down to a surface brightness threshold of
1.4 MJy/sr which collects more flux than a visual method (depending on
background structure and source morphology).
![]() |
Figure 16:
Comparison of the 60 ![]() |
Open with DEXTER |
![]() |
Figure 17:
Distribution histogram for the 100 ![]() |
Open with DEXTER |
![]() |
Figure 18:
Comparison of the 100 ![]() ![]() ![]() ![]() |
Open with DEXTER |
To summarize, all four revisited IRAS source catalogs behave remarkably similar when compared to the corresponding results given in SI89: faint sources tend to be assigned lower flux densities in our study, and bright sources (depending on the wavelength) sometimes have considerably larger flux values than found by SI89. This is the result of mainly two interfering effects which apply to sources of low and high luminosities which lie in different regions of the SMC:
The ISO and the four IRAS HiRes maps were analysed using a new computational method of detecting sources and determining source fluxes. It turned out that the HIIphot algorithm leads to an accurate and reproducible photometric characterization of the SMC. Overlapping sources in crowded fields were equally well detected as weak extended sources in regions of an inhomogeneous background. Since the ISOPHOT data from the ISO archive and the four IRAS data sets from the IPAC IRAS archive were treated in an identical manner, this allowed for a direct comparison between the ISO and the (older) IRAS data.
data set | ISO 170 ![]() |
IRAS 12 ![]() |
IRAS
25 ![]() |
IRAS 60 ![]() |
IRAS 100 ![]() |
---|---|---|---|---|---|
Appendix | E | A | B | C | D |
![]() |
243 | 73 | 135 | 384 | 338 |
155 with
![]() |
124 with
![]() |
167 with
![]() |
|||
cross-id. 170 ![]() |
-- | 21 (29%) | 61 (45%) | 180 (47%) | 162 (48%) |
cross-id. 12 ![]() |
20 (8%) | -- | 38 (28%) | 31 (8%) | 27 (8%) |
cross-id. 25 ![]() |
56 (23%) | 36 (49%) | -- | 88 (23%) | 71 (21%) |
cross-id. 60 ![]() |
153 (63%) | 29 (40%) | 90 (67%) | -- | 130 (68%) |
cross-id. 100 ![]() |
156 (64%) | 26 (36%) | 76 (56%) | 234 (61%) | -- |
![]() |
0.3...450 | 0.1...1.0 | 0.06... 3.0 | 0.06...25 | 0.08...160 |
![]() |
453/366/198/327/22 | 1.1/1.1/1.5/0.9/1.7 | 10.0/29.0/12.5/6.5/19.2 | 148/283/48/79/57 | 271/330/84/179/58 |
![]() |
99* (64%*) | 18 (25%) | 41 (30%) | 95* (77%*) | 90* (54%*) |
![]() |
7* (5%*) | 3 (4%) | 6 (4%) | 6* (5%*) | 7* (4%*) |
![]() |
19* (12%*) | -- | 9 (7%) | 9* (7%*) | 18*(11%*) |
![]() |
74* (48%*) | 2 (3%) | 24 (18%) | -- | 70* (42%*) |
![]() |
16* (10%*) | 2 (3%) | 12 (9%) | -- | 16* (10%*) |
![]() |
-- | -- | -- | 75* (60%*) | 57* (34%*) |
![]() |
-- | -- | -- | 26* (21%*) | 22* (13%*) |
![]() |
51* (33%*) | 25 (34%) | 60 (44%) | 71* (57%*) | 63* (38%*) |
![]() |
18* (12%*) | 14 (19%) | 31 (23%) | 33* (27%*) | 22* (13%*) |
![]() |
61* (40%*) | 38 (52%) | 79 (59%) | 84* (68%*) | 70* (42%*) |
![]() |
-- | 69% | 41% | 35%* | 31%* |
Entries in column #1 are
defined as follows:
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
The following statements can be made:
Acknowledgements
This work was performed at the ISOPHOT Data Centre at MPIA which is supported by Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) with funds of Bundesministerium für Bildung und Forschung, grant 50QI0201, the Max-Planck-Gesellschaft, and ESA funds supporting the scientific data reduction. This research has made use of the Simbad Database, operated at CDS, Strasbourg, France, and NASA's Astrophysics Data System Abstract Service. The authors would like to thank D.Thilker for publishing his code in the internet and thereby making it accessible for an interested public.