This field FC-01863+00035 at
is close to the Galactic
mid-plane where the interstellar extinction is serious. The measurement of
interstellar extinction becomes then impossible in optical at large
distances. Based on the optical observations of O-type to F-type
stars, the
value at the direction of
and
which is close to the ISOGAL field, is estimated to be about 3 mag and approximately constant from a distance of
1 kpc up
to
4 kpc (Neckel & Klare 1980). This may be a good reference for short
distances although the patchy distribution of interstellar matter in
the Galactic plane (Schultheis et al. 1999) may lead to some significant
difference from the extinction in the ISOGAL field
FC-01863+00035. However, the sources detected in the mid-infrared LW2 and LW3
bands extend to a much larger distance because of the much smaller
extinction at these wavelengths than in optical.
Our calculation of interstellar extinction is based on two
assumptions. One is that most of the objects in the ISOGAL field are
luminous RGB stars or AGB stars with moderate mass loss rate. These
objects are very bright with a luminosity of about one thousand to
several thousand solar luminosity, and cold with an effective
temperature lower than 4000 K or so. In addition, stellar winds from
the photosphere have formed a circumstellar envelope for many of
them. So they are strong emitters in mid-infrared. The
analysis of the ISOGAL fields at (l=0.0
,
b=1.0
)
and in Baade's Window,
has actually found most of the sources there are RGB and AGB stars
whose colors are mildly reddened by interstellar extinction
(Omont et al. 1999; Glass et al. 1999). The other assumption is that the intrinsic
colors (
)0 of all these sources are approximately similar as well as
(
-[7])0, which is not true for the minority of objects with
large mass-loss. According to the compilation of Wainscoat et al. (1992),
the value of (
J-K)0 of an M 0III star is about 1.0 and of an M 5
early AGB star is about 1.3 (see also van Loon et al. 2003 and references therein). In one ISOGAL field Sgr I that is inside
the Baade's Window (Glass et al. 1999), the concentration of (
)0and (
-[7])0 is very evident; for the objects with (
-[7]) < 1.0,
the average (
)0 is
1.08 with standard deviation of 0.17 mag and the average (
-[7])0being 0.20 with standard deviation of 0.23 mag (with the revised photometry of
the ISOGAL PSC, see Sect. 3.1) when the interstellar extinction
is subtracted as 0.2 mag in the
-band.
The assumption to start to extract the extinction law in mid-infrared
is that the intrinsic color index (
)0 of most stars in the field is
the same so that their observed color
represents the interstellar
extinction. Another presumption is the linear correlation between two
color indexes if they are both originated from interstellar
extinction, e.g. a linear correlation between the interstellar
reddening at
and
-[7], and at
and
-[15]. For the first
assumption, we had to remove some objects with apparently different
intrinsic color. In addition to RGB stars and AGB stars with moderate
mass loss, YSOs and AGB stars with high mass loss may be important
members of the mid-infrared objects. Their cold and thick
circumstellar envelopes absorb near-infrared radiation and radiate in
mid-infrared strongly and make them easily detected by the ISOGAL
survey. Their intrinsic color indexes (
-[15])0 and
(
-[7])0 [and even (
)0 for some of them] are redder
than RGB stars or early AGB
stars, and their color index [7]-[15] is also redder. Because the observed
values of
,
-[7] and
-[15] are significantly influenced by
interstellar extinction, they can not be used as a criterion to pick
up the late AGB stars or YSOs. On the other hand, the interstellar
extinction in the LW2 and LW3 bands is much smaller, the color index
[7]-[15] should be affected little by interstellar extinction. From
the ISOGAL results in the Baade's Window (Glass et al. 1999) and from
modeling (Ojha et al. 2003) we estimate that, for
[7]-[15]<0.4,
the effect of circumstellar dust of AGB stars is negligible on
(
-[7])0, that it remains small (
0.3) on (
-[15])0,
and that it is similarly negligible on (
)0 for
[7]-[15]<1.0. Therefore we exclude the stars with
or
in the discussion of the relation of
the colors
-[7] and
,
respectively, with interstellar
extinction. Applying the same criteria for young stars also warrants
that the effect of circumstellar dust is negligible on their intrinsic
colors (see e.g. Bontemps et al. 2001).
If we set the intrinsic color (
)0 of most ISOGAL stars as
and of (
-[7])0 as
,
and the
observed color indexes as
and
-[7]
accordingly,
the linear relationship is expected to be
(
-
.
Based on model calculation
(e.g. Bertelli et al. 1994) and the observation of the field Sgr I, the
intrinsic color
of RGB stars or early AGB stars is
taken to be 1.2 on average, while the intrinsic color index
of these stars is not well known (however, see van Loon et al. 2003). By leaving
as a variable and assuming that the observed
of sources with
[7]-[15] < 0.4 is mainly caused by
interstellar extinction, a robust linear fitting method, which
minimizes absolute deviation and is insensitive to large departures
for a small number of points, was used to fit the data and the result
is displayed in the color-color diagram (Fig. 4, left).
The resulted values for this linear fitting to the points decoded by
pluses and triangles in Fig. 4 are: k=0.35 and
(assuming
).
The processing of
-[15] is similar to that of
-[7]. The results
are then k=0.39 and
(assuming
)
(Fig. 4, right). In
Fig. 4, the sources with
[7]-[15]>0.4 are displayed by
diamonds though they were abandoned during the fitting process to both
the
-[7] and
-[15] with
.
Their large deviation to the right
side from the fit line further proves that they are intrinsically
redder. However, there could be a residual effect of circumstellar
dust on the value of
-[15] even for sources with
[7]-[15]<0.4,
especially for sources with large values of
which are more
luminous; it could lead to an overestimation of k by maybe
10-20%.
For a blackbody radiation at temperature 4000 K typical of an M 0 RGB
star, the Planck function yields
and
when the flux density for
is taken as
W/m2/Hz (Schuller et al. 2003),
which is consistent with the fitting results
and
.
Our values of
and
are in reasonable agreement with those from Glass et al. (1999) in Baade's
Windows
when one takes into account the correction of -0.45 mag to be
applied to the preliminary ISOGAL photometry used by Glass et al. (1999)
(Schuller et al. 2003).
The linear coefficients between
and
or
define the extinction values at 7 and 15 microns which
depend on the extinction values in near infrared. With the extinction
values in near infrared given by van de Hulst (1946), Glass (1999) and Rieke & Lebofsky (1989), the corresponding values at 7 and 15 microns are
calculated and listed in Table 2 by adopting the
coefficients k=0.35 and k=0.39 between
and
,
and
and
,
respectively.
While the extinction values in the near-infrared from Rieke &
Lebofsky (1989) are higher, the ones preferred by Glass (1999) are
just an update of those of van de Hulst (1946) and thus the two latter
sets of values are quite close to each other. Similarly, the
extinction values at 7 and 15
m derived from Rieke & Lebofsky
are higher than those
inferred from Glass-van de Hulst, i.e.
and
AV (see Table 2). As generally in
other ISOGAL papers, we prefer Glass' values for the reasons exposed
in Glass (1999), except possibly for the fields close to the
Galactic Center and in star forming regions. To compare the values we
have derived for the mid-infrared extinction with previous ones, we
will distinguish the case of 7
m where the fitting of
Fig. 4 (left) is very robust, from the one of 15
m
which is more uncertain. In both cases one should note that the
averages of the extinction on the wavelength range of the LW2 and LW3
broad bands yield slightly larger values than at 7 and 15
m,
respectively, which are both close to a minimum of the extinction
curve.
At 7
m the most widely used reference for the extinction value is probably Mathis (1990) which is similar to Draine & Lee (1984).
They used similar extinction law in near infrared to Rieke & Lebofsky (1989),
and estimated the infrared extinction from an extrapolation of the optical
extinction law and near-infrared observational data, yielding
.
This value is slightly lower than the one we derived, 0.03
.
Such a difference may arise
from the uncertainty in the ISOGAL photometry or in the fitting
procedure in Fig. 4. It is probably included as well within the
uncertainty of the extrapolation by Mathis (1990). We will use our value, 0.03
,
in the following of this paper, keeping in view its uncertainty. Let however note that the higher value we could derive from the near-infrared extinction law of Rieke & Lebofsky (1989) (Table 2), is also consistent with the results derived by Lutz (1999) (
)
who used the hydrogen recombination lines and the same near-infrared
extinction law towards the Galactic Center.
![]() |
Figure 4:
Color-color diagram of
|
The extinction ratio
A7/A15 is given as
1.3 by Mathis (1990). However, the value of Draine & Lee (1984) is almost twice smaller. The best value for this ratio has probably been derived by Hennebelle et al. (2001) from an analysis of infrared dark clouds from the ISOGAL survey. They give
for the clouds away from the Galactic Center, which is similar to Draine & Lee (1984).
The value we derived for this ratio in Table 2, Col. 2, 1.2, is more compatible with Mathis (1990) than with Hennebelle et al. (2001).
In Fig. 4 (right),
the 15
m value deduced from the Hennebelle et al. ratio is shown by a dashed line
derived from the extinction values at J,
and 7
m taken from
Col. 2
of Table 2. It is evident that their result implies more
extinction at 15
m. However, let us stress that our direct
fitting is somewhat
uncertain at 15
m because of the smaller number of sources and of the possibility of residual effects
of circumstellar reddening on
-[15], especially for the distant
luminous sources with large extinction. Therefore, we consider that
the 15
m value deduced from Hennebelle et al. (2001) (dashed line in
Fig. 4 (right)) is still compatible with the data of
our fitting in Fig. 4 (right, full line), and that our
data cannot provide a really accurate value of the extinction at 15
m.
Nevertheless, we will use the values
derived from our fitting
in the subsequent sections, since such differences on the small
extinction at 15
m are practically negligible in the following
discussions. One has to keep in mind that the extinction law may vary
with the directions due to the inhomogeneous distribution of
the interstellar matter in the Galactic plane.
| Band |
|
||
| vdH1 | ISG2 | R & L3 | |
| J | 0.245 | 0.256 | 0.281 |
| 0.087 | 0.089 | 0.112 | |
| 7 | 0.032 | 0.031 | 0.0534 |
| 15 | 0.026 | 0.025 | 0.046 |
|
1 Values for J and 2 Values for J and 3 Values for J and 4 7 and 15 (see text): |
By assuming the intrinsic color index
is the same
for all the objects which were detected in both J and
bands and
with
[7]-[15]<1.0 if detected in both LW2 and LW3 bands as well, the
interstellar extinction to individual objects can then be calculated
from the observed
.
Theoretical calculation shows that
doesn't differ much for late-type RGB stars and
early AGB stars which are the major components of the ISOGAL
sources. Though there may be some foreground main-sequence stars
and early-type AGB stars
with intrinsic bluer color than late-type RGB stars that would bring about the
underestimation of
,
the number of such sources should be
small and they should be mostly at small distance to be detectable by
ISOGAL as seen in Fig. 4 for those with
< 1.0.
While many of the ISOGAL sources are detected in the
-band and not
detected in the J-band, the value of
may then be inferred from
according to the extinction law derived above. This
method is right only if the non-detection in the J-band is caused by large
interstellar extinction. Because some of the non-detections in the J-band
may come from serious absorption by circumstellar dust of AGB
stars or YSOs, the estimation from
may overestimate the
real interstellar extinction, in particular when the index
is large, e.g. larger than 2 that results in
.
So
when inferring the extinction value from
-[7], we have excluded the
sources with
[7]-[15]>0.4 because of the risk they present the
occurrence of large mass loss. For the same reason, we have excluded
the sources with
[7]-[15]>1.0 when inferring the extinction value
from
.
However, we have thus altogether discarded a large proportion of
sources with the largest extinction.
We adopted the value of intrinsic color index
as 1.2
and the extinction values listed in Col. 2 of
Table 2. The global distribution of
is shown in
Fig. 5d, where values of
from
or
-[7]
are distinguished by the dash or dot lines respectively, while the
summation of the two types is represented by the solid line. It can
be seen, as expected, that the sources not detected in the J-band
experience higher extinction than those detected in the J-band on
average. All the sources not detected in the J-band have
and peak at about 28. In Fig. 5, we have added to the bin
all the sources with
.
Such sources should not be
AGB or RGB-tip stars, but nearby earlier stars, mostly K giants with
a few A-B stars. Most of them must have
AV<2-3. We have checked
that the corresponding values of I-J are consistent with J-K for all these
19 sources, except for one where there is a problem with the DENIS
I-J associations. As expected, all of them are associated with a
visible GSC2-2 star.
The combined visible/near-infrared colors are consistent with
K-early or M giants
for most of them, with, however, an identified A2 star.
In addition to the concentration at
,
the sources are
distributed along
unevenly. After a dip at
,
there is first a progressive increase of their
number up to
,
and then a steeper rise until a
relatively sharp maximum at
.
After that,
their number decreases rather regularly up to
.
We think that the uneven distribution of
reflects the
inhomogeneities in the distribution of the interstellar medium,
partly along the line of sight with the crossing of the molecular
ring and of several arms, but probably mainly perpendicular to the
line of sight
across the observed ISOGAL field. In order to see the influence of
the spiral arms, we estimated the extinction values to the arms from
the distribution of the emission in the radio lines of CO and HI. The
kinematical distance to the interstellar clouds can be inferred from
their radial velocity. The line of sight in the
direction l=-18.63
and b=0.35
touches the outer edges of
the Sagittarius,
the Centaurus and the Norma arm at respectively 1, 4, and 6 kpc, runs
through the bulge and then reaches the far side of the Norma arm at
about 13 kpc.
Indeed, most of the line of sight, between 3 to 13 kpc, is
in the molecular ring, including the tangential point at
8 kpc
(
km s-1).
As for the whole Galactic disk, there exist 12CO data for our
field taken from the Milky Way survey at the coarse resolution of
0.125
(7.5
)
(Dame et al. 2001). Our
field
thus implies nine pointings of this survey at
![]()
and
.
A rapid look at these CO emission data
shows that it corresponds mainly to radial velocities characteristic
of the molecular rings and more precisely, of the Centaurus and Norma
arms. However, there is an ambiguity for the emission between -70--80 km s-1 which could be attributed either to the near side or to the far
side of Norma. Especially in the Centaurus arm, the CO emission
displays a strong gradient with respect to b across this field
(Table 3). In order to take into account this gradient in the
discussion of the distribution of AV, we will consider
separately the three b ranges of the pointings of the CO survey
(Table 3). For each b range of the pointings we computed
the CO integrated intensity, W
,
averaged over the three lpointings, for each velocity interval roughly corresponding to the
different spiral arms (Table 3). Then we estimated the
corresponding AV shown in Table 3, adding a contribution
from the HI regions from the HI survey of Bloemen et al. (1990).
![]() |
Figure 5:
Histograms of the interstellar extinction
|
We estimated
shown in Table 3 by adopting the
conversion factor
from the CO integrated line intensity W
to H2 column density (Dame et al. 2001) and
the factor
from H2column density to
.
This estimation, in particular of the
HI column density, suffers some uncertainty from integrating the
velocity on unclear contours of the paper that could be about 30%.
The W
to N
conversion factor is also known
to be rather uncertain (Bachiller & Cernicharo 1986; Harjunpaa & Mattila
1996).
Similarly, we split our sample of ISOGAL sources in three unequal parts corresponding to these b ranges (Table 3). For each subsample, we have represented in Fig. 5 the histogram of the distribution of AV. We note, as expected, important differences in the total ranges of these three distributions reflecting the gradient of the extinction with b.
In view of discussing the correspondence between the
distribution from the ISOGAL sources and the determination from the
interstellar gas, let us stress the difficulty that any spatial
inhomogeneity in the CO intensity smaller than the large CO beam is
smoothed with the present CO data. Therefore, it is impossible to
estimate the actual spatial dispersion of
within a spiral
arm. It is certainly significantly larger than the dispersion,
20%, between the three l pointings for the same b range. In
each of the three histograms of Fig. 5 the average value of
the accumulated extinction expected from the interstellar gas just
beyond the Centaurus arm is,
,
18 and 12, respectively
(Table 3). Looking at
the values of the sharp rise of the
distribution
due the Centaurus arm, with
,
14 and 10, respectively,
one sees that
these are 20-50% lower than
.
A natural interpretation
is that this difference essentially represents the dispersion of the
extinction through the Centaurus arm.
The case of the subfield with
b>0.4375
is special because the b sample observed by ISOGAL,
0.4375
<b< 0.495
is about half the total range of the
CO pointing, 0.4375
<b< 0.5625
.
Because of the strong
gradient with b, it is likely that the average extinction in the
region observed by ISOGAL is larger than the one in the total CO
beam which yields
(see Table 3). Indeed
the main feature of Fig. 5c is approximately symmetric with respect
to
,
which should be close to the actual value of
the average of
in the region observed by ISOGAL.
However, it is unclear in this case whether the main
contribution of Norma is absent because it comes from the far side
with very few sources behind it, or it is included in the pedestal
between
and
,
together with the extinction of
a few sources behind the tangential point region.
The situation is less clear for the other two subsamples with
smaller b and larger extinction. It is true that the distribution
again extends beyond the value
,
symmetric of
,
with respect to
,
as expected from such a dispersion of the
distribution through the Centaurus clouds and the additional extinction beyond Centaurus. But there is a large decline
of the number of the sources towards large
,
making the
distribution very asymmetric with respect to
.
The loss of
sources may be explained by the very large value of
which
prevents the detection at 7
m or even in the
-band
of many sources.
Finally, we note that for the three subsamples, the total range of
values of
inferred from the ISOGAL sources is consistent
with the maximum accumulated extinction along the line of sight
expected from the interstellar gas, including the dispersion of AV values within clouds, the Norma ambiguity and the loss of sources with very large
.
Despite the many uncertainties in the relation between W
and the ISOGAL colours and in the interpretation of the
distribution, one feature remains intriguing. Before the sharp rise at
-14 that we associate with the Centaurus arm, there
is a distribution of a smaller number of sources from
to
12. The local and Sagittarius arm emission of CO is quite
unable to account for such extinctions. The only explanation is that
such sources are located within the Centaurus arm, behind the first
dust layers. The subsequent very sharp rise suggests a narrow
distribution of the bulk of the Centaurus extinction. The appreciable
number of sources with
to
12 suggests that
the first Centaurus layers of dust are distributed over a rather large
distance with a substantial source density. This implies a large
width for the Centaurus "arm'', which is consistent with the broad CO
velocity distribution, and/or an over-density of ISOGAL sources inside
it, which could be mostly young/massive K giants.
Copyright ESO 2003