To derive LTE column densities from the 12CO and 13CO(1-0) data,
we follow the method outlined by e.g.
Brand & Wouterloot ([1998]), and assume for the 13CO abundance
the value of
derived by Dickman ([1978]).
![]() |
Figure 9:
a) Distribution of the IRAS 60 ![]() ![]() |
The highest value for the H2 column density in MBM 32 that we derive
is
cm-2 at offset (8
, 20
).
The ratio 12CO/13CO(1-0) increases
with
(12CO(1-0)) from 5.2 (>5 K) to 6.9 (2 - 3 K),
suggesting a decrease in optical depth from the cloud center to the edge.
Due to the smaller 13CO line
widths the integrated
ratio increases from 6.9 to 11.1 in this
interval.
These ratios are higher than those found in Giant Molecular Clouds, where
they are typically 3 - 5. We cannot directly derive column densities towards
all positions because 13CO(1-0) was not measured at all
positions where we have 12CO(1-0) HRS data, and therefore we
first compared for those positions where this is possible, N(H
and
(see Fig. 11). It is seen that the ratio
=
N(H
is independent of
(though with a large spread) at a value of 7.0
2.5 (s.d.)
1019 cm-2 (K km s-1)-1. To derive cloud masses we
use this value, rather than the value of
cm-2 (K km s-1)-1
(Strong & Mattox [1996]) or
cm-2 (K km s-1)-1 (Hunter et al.
[1997])
found for inner Galaxy clouds. But it is close to the number derived from
gamma-ray data by Digel et al. ([1996]) for the Polaris Flare,
0.92
0.14 1020 cm-2 (K km s-1)-1.
Similarly we find from the 13CO(1-0) data a ratio
N(H
CO) of 5.4
cm-2 (K km s-1)-1.
This number will be used to derive masses
of clumps embedded in MBM 32.
The resulting masses in the three velocity intervals in Fig. 1a are
16.9
(2 - 7 km s-1), 4.1
(-5 - 0 km s-1), and
0.46
(0 - 2 km s-1), including a factor 1.36 for He, and
assuming a distance to the cloud of 100 pc. The cloud-averaged H2densities for these three regions are then 70, 78, and 135 cm-3
respectively. It was assumed that the clouds are spherical, which is not
very likely, and this may cause the lower mean density for the main cloud
component, compared to the other two values.
If the depth along the line of sight is smaller, the density for
this component is larger than 56 cm-3.
Maximum derived 13CO column densities are about
cm-2.
MBM 32 has a low average
extinction of AB=0.6 mag over an area of
(Heithausen & Mebold [1989]; see also Magnani & de Vries
[1986]). The maximum value derived for the
extinction (at a resolution of about 8
)
is only slightly higher.
Assuming a peak extinction AV of 0.8 mag and using
N(H) =
AV (Bertoldi & McKee [1992]),
the 13CO
abundance would be
,
equal to the Dickman value. However
the assumption of a constant 13CO abundance in a HLC such as
MBM 32 is probably not correct. Van Dishoeck & Black ([1988]) show
that near
N(H2) = 1021 cm-2 there is a strong increase of CO abundance
with N(H2). In addition part of the extinction can be related to
H I gas rather than to H2 since the cloud is clearly detected in
H I (see Figs. 5-7).
The derived J=1-0 excitation temperatures are between 5 and 9 K,
with the highest value of 9.6 K at offset (4,
24
). The
kinetic temperature of 24 K derived by Schreiber et al. ([1993])
from NH3observations therefore suggests that CO is subthermally excited. Also the
dust temperature of 20 K is higher than
.
However also beam filling
(the presence of many tiny clumps) could explain lower than expected line
temparatures.
12CO(3-2) has been observed in a much smaller region due to the
relative weakness of the lines. We compare it in the same way with the
12CO(2-1) emission as above. However, because the (3-2) transition has
not been observed to the edge of the cloud, we did not add spectra without
lines around the observed region before convolving to the lower resolution.
This transition was only mapped in two small
regions of the main cloud component (2 - 7 km s-1). The ratio (2-1)/(3-2)
of the peak
is 2.17
0.45 (N = 37) in the northern region, and
2.74
0.51 (N = 22) in the southern region. The difference in the two ratios
might be real, but part of it might also be explained by the lower efficiency
assumed for part of the measurements in the southern region (those made in
1990/1991). However for observations in this region we found no systematic
differences in intensity after efficiency correction between observations
made in 1990/1991 and in 1991/1992. The ratio (1-0)/(2-1) in both regions
are equal within the uncertainty. As for this ratio (Fig. 14), the
ratio
(2-1)/(3-2) is not dependent on the
(2-1) (see Sect. 4.3).
The ratio of the 12CO(2-1) and 13CO(2-1) emission at the
positions where the latter line was detected is shown in Fig. 14a.
We show both the values determined from the peak intensities and those from
the integrated intensity. It is seen that at the
positions of strongest 13CO emission at offsets of 5
to 8
,
the
ratio is about 11. At offsets west of this range the ration increases only
slightly to values between 11 and 13 (with the exception of the position at
offset 3
). However there is not a systematic increase in this range, as
one might expect if the 12CO optical depth decreases. This should be
confirmed by deeper integrations in the most western part of the cloud. As
in the case for the (1-0) ratios, the ratios of the peak
are
slightly higher (median 11.1) than those of the integrated
(median
14.2). Both ratios are higher than those of the (1-0) ratios for the
main cloud component (about 6 and 9 for peak and integrated
,
respectively; see Sect. 4.1) and higher than in GMCs (e.g. Kramer et al.
[1996] found a 13CO/12CO(2-1) ratio of 1.5 to 8 for
Orion B). We can
compare the (2-1) ratios with those for the negative velocity component at
offsets
in the range -1
to 4
.
Adding the emission in this range
(omitting offset 3
)
we find ratios of 19.4
4.3 and 19.0
10.1
respectively
for the peak and integrated intensities of the negative velocity component (the
peak 13CO(2-1)
is 0.028
0.015 K). For the positive
velocity
component in this interval these values are 13.4
0.4 and 13.3
1.1. If
the intrinsic ratio 12CO/13CO is equal to the typical value in the
solar neighbourhood of 76
7 (Wilson & Rood [1994]),
the optical depths in both cloud components are
very similar. Because of fractionation, the intrinsic
ratio in the outer parts of a HLC such as MBM 32, will decrease slightly
(see e.g. Fig. 6 of Turner et al. [1992]), and possibly significantly
(van Dishoeck et al. [1991] suggest a factor 2 to 3), while both the
12CO and
13CO abundance decrease to values much lower than the Dickman
([1978]) one
(for 13CO). Deriving LTE column densities from the (2-1) data along
this line, we find values that are a factor 2-4 lower than those from the
(1-0) HRS data at the same positions, while the derived excitation temperatures
are almost equal for (1-0) and (2-1). This suggests that the assumptions
used for the LTE column densities are not correct, as was found also for
other clouds, such as Orion (see e.g. Kramer et al. [1996]). Since for
deriving the used abundances similar assumptions were used, it has little
influence on the derived cloud masses.
![]() |
Figure 11:
Plot of factor
![]() ![]() ![]() ![]() ![]() |
The distribution of the 12CO(3-2) emission on a 40
raster along
the same line as in Fig. 13 for (2-1) is show in Fig. 15.
The negative velocity emission (Fig. 15a) was observed with lower
sensitivity than the positive velocity emission (Fig. 15b). An
indication for a minimum near offset +3
in the latter panel
is also seen in
Fig. 13c. In all transitions the clouds edge appears not very sharp,
but the decrease in intensity is steeper than exponential.
We compare in Fig. 14b the 12CO(2-1) and (3-2) emission after
convolution of the latter data to the (2-1) angular resolution. Both cloud
components show a line ratio which is independent of position. However the
average values for both clouds are different: 2.77
0.23 for the
positive velocity emission at
offsets larger than -5
,
and 2.09
0.21 for the negative velocity
emission. This difference is in the same direction as found for the two mapped
regions (see Sect. 4.2). Both regions were again observed in different
periods for which different telescope efficiencies were assumed.
![]() |
Figure 14:
a) The line ratios 12CO(
2-1)/13CO(2-1)
along the cut through MBM 32 in Fig. 14b. The filled squares indicate
the ratios for
the peak
![]() ![]() |
![]() |
0<
![]() |
![]() |
|
12CO(1-0)/(2-1) | 2.40 | 2.72 | 1.93 |
12CO(2-1)/(3-2) | 2.1 | 2.7 | |
12CO/13CO(1-0) (
![]() |
5.2 - 6.9 | ||
12CO/13CO(1-0) (
![]() |
7 - 11 | ||
12CO/13CO(2-1) | 19.2 | 13.4 | |
X'(lte) (for 12CO(1-0)) | 0.7 1020 a | ||
X'(lte) (for 12CO(2-1)) | 1.4 1020 a | ||
X'(lte) (for 13CO(1-0)) | 5.4 1020 a | ||
X(fir) (for 12CO(1-0)) | >0.2 1020 a | ||
a Units K km s-1cm-2. |
In Table 2 we summarize the measured average line ratios in the
different parts of the cloud as well as working values for
for the
different transitions which result in column densities equal to the LTE
column densities (for 12CO(2-1) we assumed a (1-0)/(2-1) ratio of 2
to obtain this
). The ratios can be compared with average
dv ratios listed by
Ingalls et al. ([2000]) for translucent clouds (obtained from data
by van Dishoeck et al. ([1991]): (
1-0)/(2-1)=1.30+0.59-0.29,
(2-1)/(3-2)=1.82+0.74-0.41. Their 12CO/13CO(1-0) ratio is
7.7
+3.4-1.8. Our ratios appear to be larger, but might be
consistent if calibration uncertainties are taken into account.
We have tried to analyze these ratios with an escape probability model
(Stutzki & Winnewisser [1985]). It appears that the 12CO line
ratios are
consistent with low density (n<100 cm-3) gas with
about 12 K.
Also the strongest 12CO(1-0) lines of about 6 K can be explained by such
gas. However the weaker 12CO lines with the same line ratios cannot
be explained. It is possible that this is caused by similar gas with
a smaller beam filling factor. But then the relatively high
12CO/13CO(1-0) ratios still cannot be explained.
The results of the analysis of the line intensities towards one central position
in MBM 32 by Schreiber et al. ([1993]), which indicate a flat
density distribution and the constancy of the line ratios within the whole
cloud, may suggest that the cloud consist of small clumps (much smaller than
the beam size) with average properties that do not change within the mapped
area. However the number of such clumps decreases towards the cloud edges.
The presence of small clumps within a beam was also suggested by Tauber et al.
([1991]) from high spectral resolution 12CO and 13CO(1-0)
measurements towards Orion. Ingalls et al. ([2000]) conclude from
the constant line rations in translucent clouds (including their own
12CO(4-3) data) that the clouds consist of
high density (n(H
2)=104.5 cm-3), low temperature (8 K) cells.
However we note that they only used the observed line ratios to compare with
their LVG models, no line temperatures. We found that in that case solutions
are easier to find (see above).
![]() |
![]() |
![]() |
|
![]() ![]() |
2.6 K | 1.3 K | 4.8 K |
![]() |
1.92 km s-1 | 2.1 km s-1 | 1.12 - 1.49 km s-1 |
![]() |
4.35 km s-1 | 4.75 km s-1 | 4.95 km s-1 |
![]() |
0.57 pc | 0.23 pc | 0.95 pc |
![]() |
3.04 ![]() |
0.34 ![]() |
12.4 ![]() |
![]() |
2.61 ![]() |
1.23 ![]() |
8.38 ![]() |
![]() |
0.020 ![]() |
0.053 ![]() |
Pound & Blitz ([1993]) list 3 clumps found in a still unpublished 13CO(1-0) map of MBM 32. Two of those clumps are located close to clumps listed in Table 6. Our radii are larger because of the angular resolution of the observations.
We have also tried to analyse the cloud data with more detailed algorithms, investigating for instance the power law distribution of clump masses (see e.g. Stutzki et al. [1998]). However the data cubes are not large enough to allow us to obtain reliable results.
Similar to the CO maps, we have investigated the structure of the H I
gas. We used a data cube where the broad velocity component had been
subtracted. The algorithm finds about 60 clumps larger than the 9
2
resolution. The residual map shows fairly uniform emission with lines of
several K. The largest clump found represents the 24 K emission at
3.8 km
s-1associated with the main part of the molecular cloud. However, most clumps
are located at the edge of the mapped region and therefore the data are
of no use for a detailed analysis.
In Fig. 16 we compare some of the derived parameters (size, linewidth
and mass) of the clumps listed in Tables 3 to 6.
We distinguish between the three velocity components and between the clumps
which are visible in the original maps before subtracting other Gaussian clumps
and those which are not visible in those maps.
It is seen in Fig. 16a that the 12CO(1-0) linewidths found from
the MRS data are somewhat larger than those derived from the HRS data, even
after corection for the resolution. There is no significant difference in
distribution between the filled and open symbols (+ and -, respectively,
in Col. 10 of Tables 3 to 6).
We can investigate the linear correlation between the parameters for those
velocity components and transitions where the number of clumps found is large
enough. We only use the clumps indicated by the filled symbols.
These results can be compared with correlations found for integral parameters
of HLC. For the main component (12CO(1-0) HRS data) we obtained
(correlation coefficient 0.86):
The data points in Fig. 16b show a larger correlation which is
visible in all four panels. We obtain
from the same two data sets as above the following results:
The mean slope in Eqs. (6) and (7) of about +2.3 is somewhat steeper than the one found for GMC's of 2.0 by Brand & Wouterloot ([1995]) although with a much smaller constant term, which reflects the lower densities of the clumps. The slope is equal to the one found by Heithausen et al. ([1998b]) towards the Polaris Flare from the combination of different data sets.
Copyright The European Southern Observatory (ESO)