Mapping of the_extinction_in_giant_molecular_clouds_using_optical_star_counts

A&A manuscript no.
(will be inserted by hand later) ASTRONOMY
Your thesaurus codes are: AND
08 (09.03.1 ; 09 04 1 ; 09.19.1 ) ASTROPHYSICS

Mapping of the extinction in Giant Molecular Clouds using
optical star counts
L. Cambr´sy
e
arXiv:astro-ph/9903149v1 10 Mar 1999

Observatoire de Paris, D´partement de Recherche Spatiale, F-92195 Meudon Cedex, France
e

Abstract. This paper presents large scale extinction and morphological properties of the dust grains. The star
maps of most nearby Giant Molecular Clouds of the count method was first proposed by Wolf (1923) and has
Galaxy (Lupus, ρ Ophiuchus, Scorpius, Coalsack, Taurus, been applied to Schmidt plates during several decades. It
Chamaeleon, Musca, Corona Australis, Serpens, IC 5146, consists to count the number of stars by interval of mag-
Vela, Orion, Monoceros R1 and R2, Rosette, Carina) de- nitudes (i.e. between m − 1/2 and m + 1/2) in each cell
rived from a star count method using an adaptive grid of a regular rectangular grid in an obscured area and to
and a wavelet decomposition applied to the optical data compare the result with the counts obtained in a suppos-
provided by the USNO-Precision Measuring Machine. The edly unextinguished region. In order to improve the spatial
distribution of the extinction in the clouds leads to esti- resolution, Bok (1956) proposed to make count up to the
mate their total individual masses M and their maximum completeness limiting magnitude (m ≤ mlim ). Since the
of extinction. I show that the relation between the mass number of stars counted is much larger when compared
contained within an iso–extinction contour and the extinc- to counts performed in an interval of 1 magnitude, the
tion is similar from cloud to cloud and allows the extrap- step of the grid can be reduced. Quite recently extinc-
olation of the maximum of extinction in the range 5.7 to tion maps of several southern clouds have been drawn by
25.5 magnitudes. I found that about half of the mass is Gregorio Hetem et al. (1988) using this second method.
contained in regions where the visual extinction is smaller Even more recently, Andreazza and Vilas-Boas (1996) ob-
than 1 magnitude. The star count method used on large tained extinction map of the Corona Australis and Lupus
scale (∼ 250 square degrees) is a powerful and relatively clouds. Counts were done visually using a ×30 magnifi-
straightforward method to estimate the mass of molecular cation microscope. With the digitised Schmidt plates the
complexes. A systematic study of the all sky would lead star counts method can be worked out much more eas-
to discover new clouds as I did in the Lupus complex for ily across much larger fields. The star count methods have
which I found a sixth cloud of about 104 M⊙ . also tremendously evolved thanks to the processing of dig-
ital data with high capacity computers. Cambr´sy et al.
e
Key words: ISM : clouds – ISM : dust, extinction – ISM (1997), for instance, have developed a counting method,
: structure for the DENIS data in the Chamaeleon I cloud, that takes
advantage of this new environment. The aim of this paper
is to apply this method to the optical plates digitised with
1. Introduction the USNO-PMM in a sample of Giant Molecular Clouds
to derive their extinction map (Vela, Fig. 1; Carina, Fig.
Various methods have been recently developed to estimate 2; Musca, Fig. 3; Coalsack, Fig. 4; Chamaeleon, Fig. 5;
the mass of matter contained in giant molecular clouds Corona Australis, Fig. 6; IC 5146, Fig. 7; Lupus, Fig. 8;
(GMC) using millimetric and far infrared observations ρ Ophiuchus, Fig. 9; Orion, Fig. 11; Taurus, Fig. 12; Ser-
(Boulanger et al. 1998; Mizuno et al. 1998). Neverthe- pens, Fig. 13). The method is shortly described in Sect. 2.
less, the mapping of the optical/near–infrared extinction, In Sect. 3, deduced parameters from the extinction map
based on star counts still remain the most straightforward are presented and Sect. 4 deals with individual clouds.
way to estimate the mass in form of dust grains. These
maps can be usefully compared to longer wavelength emis- 2. Star counts
sion maps in order to derive the essential physical param-
eters of the interstellar medium such as the gas to dust 2.1. Method
mass ratio, the clumpiness of the medium or the optical The star counts method is based on the comparison of lo-
Send offprint requests to: Laurent Cambr´sy,
e cal stellar densities. A drawback of the classical method
Laurent.Cambresy@obspm.fr is that it requests a grid step. If the step is too small, this

2 L. Cambr´sy: Mapping of the extinction using star counts
e

may lead to empty cells in highly extinguished regions, plates (δ ≤ −35◦ ) in blue and red. Internal photometry
and if it is too large, it results in a low spatial resolution. estimators are believed to be accurate to about 0.15 mag-
My new approach consists in fixing the number of counted nitude but systematic errors can reach 0.25 magnitude in
stars per cell rather than the step of the grid. Since un- the North and 0.5 magnitude in the South. Astrometric
certainties in star counts follow a poissonian distribution, error is typically of the order of 0.25 arcsecond. This ac-
they are independent of the local extinction with an adap- curacy is an important parameter in order to count only
tive grid where the number of stars in each cell remains once those stars which are detected twice because they are
constant. Practically, I used a fixed number of 20 stars per located in the overlap of two adjacent plates.
cell and a filtering method involving a wavelet decomposi- All the extinction maps presented here have been
tion that filters the noise. This method has been described drawn in greyscale with iso–extinction contours overlaid.
in more details by Cambr´sy (1998).
e On the right side of each map, a scale indicates the corre-
spondence between colours and visual extinction, and the
value of the contours. Stars brighter than the 4th visual
Fig. 1. Extinction map of Vela from B counts (J2000 co- magnitude (Hoffleit and Jaschek 1991) are marked with a
ordinates) filled circle.

2.2. Artifacts

Fig. 2. Extinction map of Carina from R counts (J2000 Bright stars can produce artifacts in the extinction maps.
coordinates) A very bright star, actually, produces a large disc on the
plates that prevents the detection of the fainter star (for
example, α Crux in the Coalsack or Antares in ρ Ophi-
uchus, in Fig. 4 and 9, respectively). The magnitude of
Antares is mV = 0.96, and it shows up in the extinction
Fig. 3. Extinction map of Musca from R counts (J2000 map as a disc of 35′ diameter which mimics an extinction
coordinates) of 8 magnitudes. Moreover Antares is accompanied of re-
flection nebulae that prevent source extraction. In Fig.
I have applied the method to 24 GMCs. Counts and fil- 11 the well known Orion constellation is drawn over the
tering are fully automatic. Because of the wide field (∼ 250 map and the brighter stars appear. ǫ Ori (Alnilam) the
square degrees for Orion), the counts must be corrected for central star of the constellation, free of any reflection neb-
the variation of the background stellar density with galac- ula, is represented by a disc of ∼ 14′ for a magnitude of
tic latitude. Extinction and stellar density are related by: mV = 1.7. Fortunately, these artifacts can be easily iden-
tified when the bright star is isolated. When stars are in
1 Dref (b)
Aλ = log (1) the line of sight of the obscured area, the circularity of
a D a small extinguished zone is just an indication, but the
where Aλ is the extinction at the wavelength λ, D is the only straightforward way to rule out a doubt is to make a
background stellar density, Dref the density in the refer- direct visual inspection of the Schmidt plate.
ence field (depending on the galactic latitude b), and a is Reflection nebulae are a more difficult problem to iden-
defined by : tified since they are not always circular. Each time it was
possible, I chose the R plate because the reflection is much
log(Dref ) − cst.
a= (2) lower in R than in B. B plates were preferred when R
mλ plates showed obvious important defaults (e.g. edge of the
where mλ is the magnitude at the wavelength λ. plate).
Assuming an exponential law for the stellar density,
Dref (b) = D0 e−α|b| , a linear correction with the galactic
2.3. Uncertainties
latitude b must be applied to the extinction value given
by Eq. (1). The correction consists, therefore, in subtract- Extinction estimations suffer from intrinsic and system-
ing log[Dref (b)] = log(D0 ) − α|b| log (e) to the extinction atic errors. Intrinsic uncertainties result essentially from
value Aλ (b). This operation corrects the slope of Aλ (b) the star counts itself. The obtained distribution follows
which becomes close to zero, and set the zero point of ex- a Poisson law for which the parameter is precisely the
tinction. All maps are converted into visual magnitudes number of stars counted in each cell, i.e. 20. Equation (1)
assuming an extinction law of Cardelli et al. (1989) for shows that two multiplicative factors, depending on which
AB AR
which AV = 1.337 and AV = 0.751. colour the star counts is done, are needed to convert the
Here, the USNO-PMM catalogue (Monet 1996) is used stellar density into visual extinction. These factors are a
to derive the extinction map. It results from the digiti- (the slope of the luminosity function (2)), and the conver-
sation of POSS (down to −35◦ in declination) and ESO sion factor Aλ /AV . For B and R band, the derived extinc-

L. Cambr´sy: Mapping of the extinction using star counts
e 3

Fig. 4. Extinction map of Coalsack from R counts (J2000 coordinates)

Fig. 5. Extinction map of the Chamaeleon complex from B counts (J2000 coordinates)

Fig. 6. Extinction map of Corona Australis from B counts (J2000 coordinates)

tion accuracies are +0.29 magnitudes and +0.5 magnitudes,
−0.23 −0.4 Table 1. Cloud properties. Distances are taken from lit-
respectively. erature, masses (expressed in solar masses) are defined by
Also, for highly obscured region, densities are esti- the regression line log M (AV ) = log(MTot ) + a × AV ,
mated from counts on large surfaces, typically larger than the maxima of extinction, Am , is measured from star
V
∼ 10′ . It is obvious that, in this case, the true peak of counts, and Ae , is extrapolated from the previous equa-
V
extinction is underestimated since we have only an av- tion assuming a fractal structure and the last column is
erage value. The resulting effect on the extinction map the value of the slope a
is similar to the saturation produced by bright stars on
Schmidt plates. This effect cannot be easily estimated and Cloud Name d (pc) AmV Ae
V MT ot Slope
is highly dependant on the cloud (about ∼ 80 magnitudes Lupus I 100(1) 5.3 7.1 104 -0.56
in ρ Ophiuchus, see Sect. 3.2). Lupus II 100(1) 3.8 5.7 80 -0.33
Lupus III 100(1) 4.9 7.6 1150 -0.40
Moreover, systematic errors due to the determination
Lupus IV 100(1) 5.3 7.0 630 -0.40
of the zero point of extinction are also present. Extinction
Lupus V 100(1) 5.2 10.6 2500 -0.32
mappings use larger areas than the cloud itself in order to Lupus VI 100 4.8 7.0 104 -0.57
estimate correctly the zero point. This systematic uncer- ρ Ophiuchus 120(1) 9.4 25.5 6600 -0.15
tainty can be neglected in most cases. Scorpius 120 6.4 7.0 6000 -0.54
Taurus 140(2) 7.5 15.7 1.1 104 -0.26
Coalsack 150(1)(3) 6.6 6.3 1.4 104 -0.63
3. The extinction maps and the derived Musca 150(1) 5.7 10.1 550 -0.27
parameters Chamaeleon III 150(4) 3.7 7.8 1300 -0.40
3.1. Fractal distribution of matter in molecular clouds Chamaeleon I 160(4) 5.2 12.9 1800 -0.25
Corona- 170(1) 5.4 10.7 1600 -0.30
Fractals in molecular clouds characterize a geometrical Australis
property which is the dilatation invariance (i.e. self– Chamaeleon II 178(4) 4.9 12.3 800 -0.22
similar fractals) of their structure. A fractal dimension in Serpens 259(5) 10.1 ??? ??? ???
cloud has been first found in Earth’s atmospheric clouds IC 5146 400(6) 6.5 11.9 2900 -0.29
by comparing the perimeter of a cloud with the area of Vela 500(7) 4.0 ??? ??? ???
rain. Then, using Viking images, a fractal structure for Orion 500(8) 7.5 20.3 3 105 -0.27
Martian clouds has also been found. Bazell and D´sert e Crossbones 830(8) 4.4 10.4 7.3 104 -0.47
(1988) obtained similar results for the interstellar cirrus Monoceros R2 830(8) 4.1 10.5 1.2 105 -0.48
discovered by IRAS. Hetem and L´pine (1993) used this
e Monoceros R1 1600(9) 5.4 20.8 2.7 105 -0.26
geometrical approach to generate clouds with some sta- Rosette 1600(9) 8.4 20.4 5 105 -0.28
tistical properties observed in real clouds. They showed Carina 2500(10) 8.3 82? 5.5 105 ? -0.07?
that classical models of spherical clouds can be improved (1) : Knude and Hog (1998), (2) : Kenyon at al.
by a fractal modelisation which depends only on one or (1994), (3) : Franco (1995), (4) : Whittet et al. (1997),
two free parameters. The mass spectrum of interstellar (5) : Straiˇys et al. (1996), (6) : Lada et al. (1994), (7) : Dun-
z
clouds can also be understood assuming a fractal struc- can et al. (1996), (8) : Maddalena et al. (1986), (9) : Turner
(1976), (10) : Feinstein (1995)
ture (Elmegreen and Falgarone 1996). Larson (1995) went
further, showing that the Taurus cloud also presents a
fractal structure in the distribution of its young stellar
objects. than 0.5 pc. Fractal in physics are defined over a number
Blitz and Williams (1997) claim, however, that clouds of decades and have always a lower limit.
are no longer fractal since they found a characteristic size Fractal structures in clouds can be characterized by a
scale in the Taurus cloud. They showed that the Taurus linear relation between the radius of a circle and the mass
cloud is not fractal for a size scale of 0.25-0.5 pc which may that it encompasses in a log-log diagram. Several defini-
correspond to a transition from a turbulent outer envelope tions of fractal exist, and this definition can be written
to an inner coherent core. This is not inconsistent with a M ∝ LD , where L is the radius of the circle and D the
fractal representation of the cloud for size scales greater fractal dimension of the cloud. The mass measured is, in

e

Fig. 7. Extinction map of IC5146 from R counts (J2000 coordinates)

Fig. 8. Extinction map of the Lupus complex from B counts (J2000 coordinates)

Fig. 9. Extinction map of ρ Ophiuchus and Scorpius clouds from R counts (J2000 coordinates)

in which the extinction is, in fact, greater. In the Tau-
rus cloud, I found that this turn off occurs for a size of
∼ 1.7 pc. According to Blitz and Williams (1997) a turn
off toward higher masses should appear for a size scale of
∼ 0.5 pc. Obviously, our value corresponds to a limitation
of the star counts method with optical data and not to a
real characteristic size scale of the cloud.
It is therefore natural to extrapolate, down to a min-
imum mass, the linear part of the relation log [M (Av )]
versus AV to determine a maximum of extinction. This
maximum is obtained using the regression line (3) and
represents the higher extinction that can be measured,
would the cloud be fractal at all size scales. As Blitz and
Williams (1997) have shown, there is a characteristic size
scale above which the density profile becomes steeper. The
Fig. 10. Mass contained inside the iso–extinction contours extrapolation of the linear relation gives, therefore, a lower
versus the extinction (solid line), and regression line for limit for the densest core extinction. Derived values of AV
the linear part. Annotations indicate the area in square are presented in the 4th column of Table 1 and correspond
degrees contained by the iso–extinction contours AV . to a minimum mass of 1M⊙ . This minimum mass is a typi-
cal stellar mass and, an extrapolation toward lower masses
fact, contained in a cylinder of base radius L and of unde- would be meaningless.
fined height H (because the cloud depth is not constant Except for the Carina cloud, maxima are found in the
over the surface of the base of the cylinder). In our case, range from 5.7 to 25.5 magnitudes of visual extinction
we are interested in the relation between the mass and the with a median value of 10.6 magnitudes.
extinction. Since the extinction is related to the size H, We stress the point that extinction can be larger. The
which represents the depth of the cloud as defined above, ρ Ophiuchus cloud ,for example, is known to show extinc-
we seek for a relation between mass and AV (or H), L be- tion peaks of about ∼ 100 magnitudes (Casanova et al.
ing now undefined. The logarithm of the mass is found to 1995), whereas we obtain only 25.5 magnitudes. Never-
vary linearly with the extinction over a range of extinction theless, this value compared to the 9.4 magnitudes effec-
magnitudes (Fig. 10). We have: tively measured indicates that we need deeper optical ob-
servations or near–infrared data — such as those provided
log M = log MT ot + slope × AV (3)
by the DENIS survey — to investigate more deeply the
This result is compatible with a fractal structure of the cloud. The Coalsack and Scorpius are the only clouds for
cloud if AV ∝ log H, i.e. if the density of matter follows a which measured and extrapolated extinctions are similar
power law, which is, precisely, what is used in modelling (see Sect. 4).
interstellar clouds (Bernard et al. 1993).
3.3. Mass
3.2. Maximum of extinction
Assuming a gas to dust ratio, the mass of a cloud can be
Fig. 10 shows the relation between iso–extinction contours obtained using the relation (Dickman 1978) :
and the logarithm of the mass contained inside these con-
tours for the Taurus cloud (see extinction map in Fig. 12). NH
M = (αd)2 µ AV (i)
The relation is linear for AV < 5.5. For higher extinction,
∼ AV i
mass is deficient because the star count leads to under-
estimating the extinction. Indeed, for highly extinguished where α is the angular size of a pixel map, d the dis-
regions, the low number density of stars requires a larger tance to the cloud, µ the mean molecular weight cor-
area to pick up enough stars and estimate the extinction. rected for helium abundance, and i is a pixel of the extinc-
The result is therefore an average value over a large area tion map. Uncertainties on the determination of masses

e 5

come essentially from the distance which is always diffi- imum of extinction was estimated to be ∼ 10 magnitudes.
cult to evaluate. Assuming a correct distance estimation, Using B star counts we find here 5.2 magnitudes. This dif-
error resulting from magnitude uncertainties can be eval- ference is normal since J is less sensitive to extinction than
uated: an underestimation of 0.5 magnitude of visual ex- B. The important remark is that the extrapolated value
tinction implies a reduction of the total mass of a fac- for the extinction derived from B star counts is 12.9, con-
tor ∼ 2. According to Savage and Mathis (1979), the gas sistent with the value obtained with J star counts. We ob-
to dust ratio is NH = 1.87 × 1021 cm−2 .mag−1 where
AV tain the same result for the Chamaeleon II cloud for which
NH = NHI + 2NH2 . Kim and Martin (1996) show that J star counts lead also to a maximum of about ∼ 10 mag-
this value depends on the total to selective extinction ra- nitudes whereas the extrapolated value from B counts is
tio RV = AV /EB−V . For RV = 5.3, the gas to dust ratio 12.3 magnitudes. Infrared data are definitely necessary to
would be divided by a factor 1.2. The value of RV is sup- investigate the cores of the clouds.
posed to be larger in molecular clouds than in the general Besides, the shape of the whole Musca–Chamaeleon
interstellar medium but the variation with the extinction extinction and the IRAS 100µm maps are very similar.
is not clearly established. So, I used the general value of Disregarding the far–infrared gradient produced by the
3.1 and the gas to dust ratio of Savage and Mathis. Mass heating by the galactic plane, there is a good match of
of the cores of the clouds may, therefore, be overestimated the far–infrared emission and extinction contours and fila-
by a factor ∼ 1.2. mentary connections. This region is well adapted to study
In Fig. 10, the relation log [M (Av )] vs. AV extrapo- the correlation between the extinction and the far–infrared
lated toward the zero extinction gives an estimation of emission because there is no massive stars which heat the
the total mass of the cloud using Eq. (3). Masses obtained dust. The 100µm flux can, therefore, be converted in a
are shown in the Table 1. The median mass is 2900M⊙ relatively straightforward way into a column density unit
and the range is from 80 to 5 105 M⊙ . Using expression (3), using the 60/100µm colour temperature (Boulanger et al.
we also remark that half of the total mass is located out- 1998).
side the iso–extinction curve 1.0 magnitudes. This value
is remarkably stable from cloud to cloud with a standard
4.4. Coalsack and Scorpius
deviation of 0.3.
The extinction map of the Coalsack is displayed in Fig.
4. Remarks on individual clouds 4. The edge of the cloud contains the brightest star of
the Southern Cross (α = 12h 26m 36s , δ = −63◦ 05′ 57′′ ).
4.1. Vela and Serpens This cloud is known to be a conglomerate of dust mate-
In the Vela and the Serpens clouds (Fig. 1 and 13, respec- rial, its distance is therefore difficult to estimate. Franco
tively), there is no linear relation between log [M (Av )] and (1995) using Str¨mgren photometry gives a distance of
o
AV . Extrapolation for the maximum of extinction or for 150-200 pc. More recently Knude and Hog (1998) estimate
the total mass estimations is not possible. However, mass a distance of 100-150 pc using Hipparcos data. Finally,
lower limits can be obtained using the extinction map di- I adopted an intermediate value of 150 pc to derive the
rectly: 5.7 104 M⊙ and 1.1 105 M⊙ for Vela and Serpens, mass of the Coalsack. Maximum extinction estimations
respectively. It is difficult to understand why there is no for this cloud are 6.6 and 6.3 for the measured and the
linear relation for these two clouds, even for low values of extrapolated values, respectively. The lowest limit for the
extinction. maximum of extinction, as defined in Sect. 3.2, is reached,
but no characteristic size scale has been found by studying
the shape of log [M (Av )] at high extinction. These values
4.2. Carina are too close, regarding their uncertainties, to reflect any
The study of the Carina (Fig. 2) presents aberrant values evidence of a clumpy structure.
for the slope of the line log [M (Av )] (see Table 1). Conse- Nyman et al. (1989) have made a CO survey of the
quently, the maximum of the extrapolated extinction, 82 Coalsack cloud. They have divided the cloud into 4 re-
magnitudes, cannot be trusted. The important reflection gions. Regions I and II which correspond to the northern
in the Carina region is probably responsible for the shape part of the area near α Crux, are well correlated with
of the extinction map. Stars cannot be detected because of the extinction map. Nyman et al. have defined two other
the reflection and thus, extinction cannot be derived from regions which have no obvious counterpart in extinction.
R star counts. Infrared data are requested to eliminate the Region III below -64◦ of declination in the western part
contribution of the nebulae. of the cloud is not seen in the extinction map and this
may result of a scanning defect of the plates. The same
problem exists for the region IV which is a filament in the
4.3. Musca–Chamaeleon
eastern part of the cloud at δ ≃ −64◦ .
The Chamaeleon I has already been mapped with DENIS The Scorpius cloud is located near the ρ Ophiuchus
star counts in J band (Cambr´sy et al. 1997) and the max-
e cloud (Fig. 9). The 120 pc distance used to derive its mass

e

Fig. 11. Extinction map of Orion, Monoceros I, Rosette, Monoceros II, Crossbones from R counts (J2000 coordinates)

Fig. 12. Extinction map of Taurus from R counts (J2000 coordinates)

Fig. 13. Extinction map of Serpens from R counts (J2000 coordinates)

is the ρ Ophiuchus distance (Knude and Hog 1998). As for in these regions. The complex has been first separated into
the Coalsack cloud, measured and extrapolated extinction 4 clouds (Schwartz 1977), and then, a fifth cloud has been
are similar : 6.4 and 7.0, respectively. But, no evidence of recently discovered using 13 CO survey (Tachihara et al.
a characteristic size scale can be found. For both clouds, 1996). I have discovered, here, a sixth cloud which hap-
this result is not surprising since the maximum of mea- pens to be as massive as the Lupus I cloud, ∼ 104 M⊙
sured extinction reaches the extrapolated value. Would (assuming it is also located at 100 pc). The measured ex-
dense cores with steeper extinction profile be found, the tinction for this cloud reaches 4.8 magnitudes.
measured extinction would have been significantly greater The comparison between the 13 CO and the extinction
than the extrapolated value. map is striking, especially for the Lupus I cloud for which
each core detected in the molecular observations has a
counterpart in extinction. Mass estimations can be com-
4.5. Corona Australis
pared on condition that the same field is used for both
The extinction in this cloud has recently been derived us- maps. Moreover, 13 CO observations are less sensitive than
ing star counts on B plates by Andreazza and Vilas-Boas B star counts for low extinction. The lower contour in the
13
(1996). Our methods are very similar and we obtain, there- CO map corresponds to a visual extinction of ∼ 2 magni-
fore, comparable results. The main difference comes from tudes. Using this iso–extinction contour to define the edge
the most extinguished core. Since they use a regular grid, of the cloud and the same distance (150 pc) as Tachihara
they cannot investigate cores where the mean distance et al. (1996), I find a mass of ∼ 1300M⊙ in agreement
between two stars is greater than their grid step. Con- with their estimation of 1200M⊙. Murphy et al. (1986)
sequently, they obtained a plateau where I find 4 distinct estimate the mass of the whole complex to be ∼ 3 104 M⊙
cores. using 12 CO observations and a distance of 130 pc. Using
the same distance, I would obtain 4 104M⊙ (2.3 104 M⊙
for a 100 pc distance).
4.6. IC5146
CO observations are presented in Lada et al. (1994). The 4.8. ρ Ophiuchus
two eastern cores in the 13 CO map have only one coun-
terpart in the extinction map because the bright nebula, Because of the important star formation activity of the
Lynds 424, prevents the star detections in that region. inner part of the cloud, the IRAS flux at 100µm shows dif-
Lada et al. (1994) also present an extinction map of a ferent structures of those seen in the extinction map. On
part of the cloud derived from H − K colour excess obser- the other hand, the 3 large filaments are present in both
vations. This colour is particularly well adapted for such maps. Even if these regions are complex because several
investigations because, in one hand, infrared wavelengths stars heat them, the comparison between the far–infrared
allow deeper studies and, on the other hand, H −K colour emission and the extinction should allow to derive the
has a small dispersion versus the spectral type of stars. I dust temperature and a 3-dimensional representation of
obtain similar low iso–extinction contours but they reach the cloud and of the stars involved in the heating. Unfor-
a much greater maximum of visual extinction of about 20 tunately, uncertainties about the distances for these stars
magnitudes. are too large (about 15%) and corresponds roughly to the
cloud size.
4.7. Lupus
4.9. Orion
New estimations of distance using Hipparcos data (Knude
and Hog 1998) have led to locate the Lupus complex (Fig. Maddalena et al (1986) have published a large scale CO
8) at only 100 pc from the Sun. Lupus turns out to be the map of Orion and Monoceros R2. Masses derived from CO
most nearby star–forming cloud. This distance is used to emission are consistent with those I obtain : 1.9 105 M⊙
estimate the complex mass but, it is important to remark and 0.86 105M⊙ for Orion and Monoceros R2, respec-
that evidence of reddening suggests that dust material is tively, from CO data and 3 105M⊙ and 1.2 105M⊙ from
present up to a distance of about 170 pc (Franco 1990). the extinction maps. The Orion B maps looks very alike.
Masses could therefore be underestimated by a factor < 3
∼ The Orion A maps show a significant difference near the

e 7

Trapezium (α = 5h 35m , δ = −5◦ 23′ ) where the young optical band so this problem may be neglected for opti-
stars pollute the star counts. The correlation between CO cal counts. This is no longer true with near–infrared data
and extinction maps for Monoceros R2 is less striking, for which young objects must be removed before the star
because of the star forming activity. It is clear that star counts. We obtain extinction map with a spatial resolution
clusters involve an underestimation of the extinction, but always adapted to the local densities which are typically
I would like to stress the point that the heating by a star about ∼ 1′ for the outer part of cloud and ∼ 10′ for the
cluster may also destroy the CO molecules. Estimation of most extinguished regions where the stellar density be-
the column density from star counts and from CO obser- <
comes very low, i.e. ∼ 1000 stars.deg−2 . These maps allow
vations may, therefore, be substantially underestimated in the estimation of the total mass of the cloud by extrap-
regions such as the Trapezium. olation of the distribution of matter with the extinction.
It appears that mass concentrated in the regions of low
extinction represents an important part of the total mass
4.10. Taurus
of a cloud (1/2 is contained in regions of extinction lower
Onishi et al. (1996) have studied the cores in the Taurus than 1 magnitude).
cloud using a C18 O survey. All of the 40 cores identified Extrapolation of the distribution of matter in highly
in their survey are also detected in the extinction map. extinguished areas is more risky. Star counts method give
Moreover, Abergel et al. (1994) have shown a strong cor- a relation for which masses are underestimated in the
relation between the far–infrared and the 13 CO emission cores of the cloud for a well understood reason: estima-
in that region. Despite the complexity of the Taurus struc- tion of the local density requires to pick up enough stars
ture (filaments, cores), it is a region, like the Chamaeleon and thus, to use larger area because of the low number
complex, located at high galactic latitude (b ≃ −16◦), density. Moreover a characteristic size scale in the distri-
without complex stellar radiation field. This situation is bution of matter (Blitz and Williams 1997) indicates the
highly favourable to a large scale comparison of CO, far– presence of the lower limit of the fractal cloud structure.
infrared and extinction maps. For this size scale the linear extrapolation used in Fig. 10
also underestimates the real mass and extinction. Despite
5. Conclusion this difficulty, the extrapolated extinction is useful to es-
timate the saturation level in the extinction map, but it
Star counts technique is used since the beginning of the is important to keep in mind that extinction can be much
century and is still a very powerful way to investigate the larger in small cores.
distribution of solid matter in molecular clouds. Now, with Finally, the examination of the relation between mass
the development of digital data, this technique become and extinction is useful to check what we are measuring.
easy to use and can probe much larger areas. For all re- The Carina cloud show an aberrant slope (Table 1) which
gions, we have assumed that all stars were background is a strong indication that the map cannot be directly
stars. The error resulting from this hypothesis can easily interpreted as an extinction map. In that case, the elim-
be estimated. Equation (1) can be written : ination of reflection nebulae is probably a solution and
1 S therefore, near-infrared data are requested to investigate
Aλ = log + cste(Dref ) the extinction. For the Vela and the Serpens cloud, the
a nb
absence of linear part is not understood and reflection is
where S is the surface which contains nb stars. If 50% not the solution in these regions.
of the stars are foreground stars, the difference between
the real extinction and the extinction which assumes all Acknowledgements. I warmly thank N. Epchtein for initiat-
background stars is : ing this study and for his critical reading of the manuscript
which helped to clarify this paper. The Centre de Donnés as-
e
1 tronomiques de Strasbourg (CDS) is also thanked for accessing
∆Aλ = log 2
a to the USNO data.
That corresponds to ∼ 0.6 or ∼ 1.1 magnitudes of vi-
sual extinction whether star counts are done using B or R
References
band, respectively. Fortunately, most of the clouds are lo-
cated at small or intermediate distances to the Sun (except Abergel A., Boulanger F., Mizuno A., et al., 1994, ApJ
the Carina at 2500 pc) and this effect is probably small, 423, L59
at least for low extinction. The good agreement between Andreazza C.M., Vilas-Boas J.W.S., 1996, A&AS 116, 21
mass derived from extinction or from CO data argues in Bazell D., D´sert F.X., 1988, ApJ 333, 353
e
that favour. Bernard J.P., Boulanger F., Puget J.L., 1993, A&A 277,
Stars physically associated to the clouds are more prob- 609
lematic because they are located precisely close to the ex- Blitz L., Williams J.P., 1997, ApJ 488, L145
tinction cores. Young objects are generally faint in the Bok B.J., 1956, AJ 61, 309

e

Boulanger F., Bronfman L., Dame T., et al., 1998, A&A
332, 273
Cambr´sy L., 1998, The Impact of Near Infrared Sky
e
Surveys on Galactic and Extragalactic Astronomy, ed.
N. Epchtein, vol. 230 of ASSL series, 157–163, Kluwer
Academic Publishers
Cambr´sy L., Epchtein N., Copet E., et al., 1997, A&A
e
324, L5
Cardelli J., Geoffrey C., Mathis J., 1989, ApJ 345, 245
Casanova S., Montmerle T., Feigelson E.D., et al., 1995,
ApJ 439, 752
Dickman R.L., 1978, AJ 83, 363
Duncan A.R., Stewart R.T., Haynes R.F., et al., 1996,
MNRAS 280, 252
Elmegreen B.G., Falgarone E., 1996, ApJ 471, 816
Feinstein A., 1995, Revista Mexicana de Astronomia y
Astrofisica Conference Series 2, 57
Franco G.A.P., 1990, A&A 227, 499
Franco G.A.P., 1995, A&AS 114, 105
Gregorio Hetem J., Sanzovo G., L´pine J., 1988, A&AS
e
76, 347
Hetem A. J., L´pine J.R.D., 1993, A&A 270, 451
e
Hoffleit D., Jaschek C., 1991, New Haven, Conn.: Yale
University Observatory, 5th rev.ed., edited by Hoffleit,
D.; Jaschek, C. (coll.)
Kenyon S.J., Dobrzycka D., Hartmann L., 1994, AJ 108,
1872
Kim S.H., Martin P., 1996, ApJ 462, 296
Knude J., Hog E., 1998, A&A 338, 897
Lada C., Lada E., Clemens D., et al., 1994, ApJ 429, 694
Larson R.B., 1995, MNRAS 272, 213
Maddalena R.J., Morris M., Moscowitz J., et al., 1986,
ApJ 303, 375
Mizuno A., Hayakawa T., Yamaguchi N., et al., 1998, ApJ
507, L83
Monet D., 1996, BAAS 188, 5404
Murphy D.C., Cohen R., May J., 1986, A&A 167, 234
Nyman L.A., Bronfman L., Thaddeus P., 1989, A&A 216,
185
Onishi T., Mizuno A., Kawamura A., et al., 1996, ApJ
465, 815
Savage B.D., Mathis J.S., 1979, ARA&A 17, 73
Schwartz R.D., 1977, ApJS 35, 161
z ˇ
Straiˇys V., Cernis K., Bartaˇi¯ t˙ S., 1996, Baltic Astron-
su e
omy 5, 125
Tachihara K., Dobashi K., Mizuno A., et al., 1996, PASJ
48, 489
Turner D.G., 1976, ApJ 210, 65
Whittet D., Prusti T., Franco G., et al., 1997, A&A 327,
1194
Wolf M., 1923, Astron. Nachr. 219, 109

This figure "fig1.gif" is available in "gif" format from:

http://arXiv.org/ps/astro-ph/9903149v1

Mapping of the_extinction_in_giant_molecular_clouds_using_optical_star_counts

More Related Content

What's hot (20)

Viewers also liked (16)

Similar to Mapping of the_extinction_in_giant_molecular_clouds_using_optical_star_counts (20)

More from Sérgio Sacani (20)

Mapping of the_extinction_in_giant_molecular_clouds_using_optical_star_counts