A&A 448, 213-219 (2006)
DOI: 10.1051/0004-6361:20042440
B. Aryal - R. Weinberger
Institut für Astrophysik, Universität Innsbruck, Technikerstraße 25, 6020 Innsbruck
Received 26 November 2004 / Accepted 4 October 2005
Abstract
We present a new isolated interstellar nebula (RA =
08
,
Dec = +25
(J2000);
= 197
8, b = +31
6) found at 100
m and 60
m on IRAS maps. It has dimensions of
140
and a cone-like shape, suggesting interaction with ambient matter or external radiation. The nebula contains 2 bright condensations ("nuclei'') and several prominent filaments,
the latter being approximately parallel to each other. We carried
out preliminary studies based on IRAS data and the Palomar
Observatory Sky Survey and found that there are no hints of star
formation in the nebula. The nebula's long axis is almost parallel
to the Galactic plane. A shaping due to the nebula's motion
through the interstellar medium might be not a sufficient reason
in this case. We have discussed this possibility by applying
results of Reynolds number hydrodynamics. During a search for
possible stellar candidates for shaping this nebula we noted a remarkable position of one of the very few nearby (D
0.36 kpc) pulsars known at about this latitude: PSR B0823+26. This
pulsar is projected at the tip of the nebula and its proper motion
is approaching along the long axis (
180
to the
direction of the cone) of the nebula. The true cause for the
shaping of the nebula, which was coined Skeleton Nebula by us,
however remains unknown.
Key words: infrared: ISM - ISM: structure
Comet-like or cone-like morphologies of (large) interstellar
nebulae are usually considered to be clear hints of interaction
processes. The number of objects of this type is rather small, as
was demonstrated by Odenwald & Rickard (1987) and
Odenwald (1988), who had performed a visual survey of
all IRAS sky maps in the 100 m bandpass. They discovered 15 high Galactic latitude clouds with comet-like appearance and discussed their far-IR and optical properties. It was shown that
the shapes of some of these clouds could be explained via their
motion through the interstellar medium. In order to enlarge this
sample and to find possible other shaping mechanisms, any
additional detection of nebulae of this morphological type may be
of interest.
In this paper, we present a large new high latitude cone-like
cloud found by us at 100 m and 60
m IRAS images. A
description of this nebula as well as some physical properties are
given in Sects. 2 and 3. A discussion of some stellar objects
that might shape the nebula is presented in Sect. 4. Finally, a brief hydrodynamic discussion and our conclusions are given in Sects. 5 and 6, respectively.
![]() |
Figure 1:
A 4.5![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Figures 1a,b show 4.5
4.5
IRAS 100 and
60
m images centered at RA = 08
,
Dec = +25
(J2000), i.e. at
=
197
79, b = +31
58. A remarkable A-shaped nebula with
a skeleton-like morphology (hence dubbed as "Skeleton Nebula'' by us, hereafter SkN) with a size of
2
5 can be seen. In the optical, very faint and diffuse emission is visible on a Digital Sky Survey (DSS) image (Fig. 1c).
The two nuclei (black extended regions in Figs. 1a,b and + signs
in Figs. 1c,d,e) might be dense portions of the nebula containing
embedded young stars, i.e. the whole object might be a (nearby)
star forming region. To check this possibility, we have carried
out star counts in order to compare the stellar density at the
location of the nuclei with the stellar densities of locations
outside. In case of star formation we would expect noticeable dust
extinction in these nuclei, i.e. less stars due to the obscuration
or dimming of background stars or perhaps a concentration of newly
born stars within the nuclei. The counts were done visually using
a microscope (25
magnification) on a Second Palomar
Observatory Sky Survey blue (POSS II B) film copy. Each grid size
was 15
15
(i.e. an area comparable with
that of each of the two nuclei). A total of 17 grids were chosen
for the star counts: one each at the location of the nuclei, 5 at
the filaments of the SkN, and 10 outside of the nebula. No stellar
concentration was found within the nuclei.
The average number of stars counted at the locations of the
nuclei, of the filaments and in the regions around the SkN were
107, 109, and 103 respectively. The observational and statistical
errors were 5% and
9%, respectively. Hence, the
difference is well within the error limit. In other words, the
stellar density is the same at the locations of the nuclei and at
the other positions within and outside of the nebula. This (small)
difference also remains within the
5% error limit when the
counts were repeated by increasing and then decreasing the grid
size by 25%. These results suggest that dust extinction is not
high in the optical, since it does not influence the number of
visible background stars. Hence we found that the SkN probably
does not represent a star forming region.
The central region of the SkN including the nuclei was examined
using R and B film copies from POSS II, (Figs. 1d,e): there, very
faint diffuse extended structures having nearly constant optical
surface brightness (24 mag/arcsec2) are present at and
around the position of these nuclei. The heart-like morphology of
the southern nucleus found at 100
m and 60
m can be seen
in POSS II R too, and has about the same size there. Taken as a whole, the Skeleton Nebula is hardly brighter in R compared to B. This suggests a reflection nebula, i.e. a nebula containing some
dust, but a certain contribution of faint optical emission lines
in the red wavelength range cannot be ruled out.
The nebula belongs to the "Category II cloud'' of Odenwald
(1988). Odenwald carried out optical and far-IR
studies on 14 high Galactic latitude clouds with comet-like or
filamentary appearances. Category II clouds, appearing as faint,
diffuse reflection nebulae of nearly constant optical surface
brightness but showing a clumpy internal structure, comprise the
objects G64-26, G208-28, G225-66, G228-27, G239-15 as
well as the Draco cloud, G90+38. Odenwald (1988)
estimates that category II clouds have relatively low masses in
the range 3-40
and surface brightness typical for
clouds illuminated by the general interstellar radiation field.
![]() |
Figure 2:
100 ![]() ![]() ![]() ![]() |
Open with DEXTER |
The SkN's two bright nuclei are located at RA = 08
,
Dec = +25
(J2000) and RA = 08
,
Dec =
+26
(J2000), respectively. The angular
separation between the nuclei is
39
.
At 100
m
(Fig. 1a), the northern nucleus (N1) displays 3 bright, short
(0
), "arms'' apparently originating in its brightest, that
is central, region. One extends towards the north-east, the second
one towards the south-west, and the third is perpendicular,
extending towards the south-east. Both eastern and western arms
further sharply (by
90
)
bend southwards. Together
with the above perpendicular arm all these extensions look like an upside down small letter
.
The southern nucleus (N2) bears some resemblance to a scaled
up version of N1 and its arms (Figs. 1a,b). The longest arm
(in the west) extends from north to south. The apparent size of
this filament is 80
.
It should be noted that all
arms, i.e. all filamentary structures associated with both nuclei
are approximately parallel (within
15
)
to each other and
extend towards south or south-east. Hence, the remarkable
structure of the SkN indicates that some external force is
influencing its shape. The source of this force seems to be
located in the north-west of the nebula.
The SkN is not only discernible at 100 m, but also at 60
m. More details than in Fig. 1 can be seen in Fig. 2, which displays 100
m and 60
m HIRES contour maps (resolution
1
). In Figs. 2a,b one notices that i) both nuclei are connected to the filamentary structures mentioned above; ii) these filamentary structures deviate towards south; iii) these structures are approximately parallel to each other; iv) there
is a bridge-like connection between both nuclei; v) the western
filamentary structures are much more extended than the eastern
ones; and vi) knots can be seen in the western filament and in the northern region.
Although it is possible that the SkN's long axis is in the sky
plane, it is more probable that it is inclined by a certain angle
with respect to the plane of the sky. Possibly, N1 and N2 might, seen from the location of an assumed shaping stellar source in the north-east, be at the same distance from the
source. If so, how strongly must the SkN be inclined? The
inclination angle i (angle between the line of sight and the
normal vector of the plane of the nebula) can be estimated using
cos2i = [(b/a)2-0.22]/(1-0.22) with b/a as the
measured axial ratio (Holmberg 1946). The inclination
angle (i) is found to be 54
.
Thus, our nebula is
neither face-on (i
0
)
nor edge-on (i
90
)
object. The SkN will be edge-on (minor axis/major axis
0.20, i = 90
)
if we rotate it by
36
(i.e. 90
-
)
from the southern to the northern direction (assuming that the nebula is in the sky plane). At this angle, both nuclei would lie in a plane and any
discrete object which would be located at the head of the nebula in the north-east would almost lie above in between the nuclei. In such a case the true size of the skeleton nebula would be larger
by a factor of about 2.5
the apparent size.
There is a rapid decrease of the nebula's surface brightness
towards shorter infrared wavelengths as is evident by comparing
the nebula at 100 m, 60
m, 25
m and 12
m. All
traces of the filamentary structures as well as both nuclei have
vanished at 25
m, whereas a very faint remnant of N2remains visible at 12
m. We have obtained far-IR surface
brightness measurements at each of the four passbands for the
nebular region using IRAS images from Groningen IRAS server (SRON,
Groningen). The colour correction depends on the shape of the
intrinsic energy distribution and on the details of the wavelength
response of the system. Corrections for a number of input energy
distributions are given in Table 6, Chap. 6 of the "IRAS
Explanatory Supplement'' (1985). The ratio of flux densities
before the colour correction corresponds to 240 K of intrinsic
temperature. The correction factor corresponding to this
temperature as given in Table 6 is used for the colour correction.
A comparison of far-IR spectral distributions of our nebula with
the Draco cloud (Odenwald & Rickard 1987) and the
nebula MBM 20 (Weiland et al. 1986) is shown in Fig. 3.
The SkN (solid line in Fig. 3) does not have a strong 12 m
and 25
m spectral component. This lack of emission at 12
m and 25
m suggests that the dust population may be relatively devoid of a small grain component as compared to the
cirrus material (Puget et al. 1985), or perhaps that the
small dust grain are not being adequately excited by the ambient
radiation field (Mebold et al. 1985). The strength of the
stellar radiation field near the Draco cloud is estimated to be
about 10
lower than in the Galactic plane (Mebold et al.
1985). The Draco cloud is at 38
latitude, at
slightly higher latitude than the SkN (32
). We estimate
that the stellar radiation field for the SkN is also weaker
(
10-12
)
with respect to a position in the Galactic
plane not only because of its high latitude but also due to the
smaller number of stars in this region due to the longitude of
.
![]() |
Figure 3: Far-IR spectral distributions of the Skeleton Nebula (solid line). The major cloud component (DRACO-Head) and a segment of the plume (DRACO-Plume) of the Draco cloud are represented by dotted and dashed lines, respectively (Odenwald & Rickard 1987). A dot-dashed line is used for the MBM 20 "cirrus cloud'' (position No. 1 from Table 1, Weiland et al. 1986). All intensities have been background subtracted and colour corrected for 240 K blackbody emission. The spectra have been rescaled by dividing the intensities for the Skeleton Nebula, Draco head, Draco plume, and MBM 20 by 5.5, 45.6, 9.5, and 14.0, respectively. |
Open with DEXTER |
We used SkyView (http://skyview.gsfc.nasa.gov/) in order to
search for counterparts of the SkN at H
using the all-sky
H
map, a composite of the Virginia Tech Spectral line
Survey (VTSS) in the north and the Southern H-Alpha Sky Survey
Atlas (SHASSA) in the south, and in X-rays (ROSAT All Sky Survey
(RASS3), RXTE/PCA All Sky Slew Survey (XSS)), as well as at radio
wavelengths (Sydney University Molonglo Sky Survey (SUMSS), the
VLA FIRST Survey, and the NRAO VLA Sky Survey (NVSS)), but no
emission could be detected in any of these.
We have determined the dust colour temperature ()
of
the nebula using Henning et al. (1990). The average dust
colour temperature in the region was determined to be 27
4 K.
We also determined the dust temperature distribution throughout
the nebular region by dividing the 100
m map by the 60
m
map and comparing the resulting I(100
m)/I(60
m) values at
each map location with the values given for dust grain models by
Dwek (1986). We found that the dust colour temperature of
the southern nucleus is 34
4 K and the temperature of its
northern counterpart 32
4 K. The eastern filamentary
structures are slightly cooler (20
2 K) than the western
filamentary structures (26
3 K). Hence, the locations of the
maximum emission i.e. the nuclei of the SkN are warmer than their filaments.
The distance of the nebula is unknown. In our case, the star count method (Mebold et al. 1985) does not help because of insufficient extinction and too few stars in the region. The nebula also has no obviously embedded stars, which would be helpful for this purpose. Following Odenwald (1988), a distance of 200 pc will be assumed for the moment. The maximum physical size of the nebula at this distance then is 7.6 pc, if it lies in the sky plane, and the separation between the northern (N1) and southern (N2) nucleus is about 2 pc.
We estimated the total mass of the SkN using the dust mass
relation by Hildebrand (1983). Based on fluxes of 3.22
104 MJy/str at 100
m, and assuming a gas-to-dust mass ratio (g:d) of 150, D = 0.2 kpc, and
= 27
4 K, the total mass is 3.5
2.3
.
This
value is parameterized in terms of their assumed distance (in kpc), temperature (in K) and gas-to-dust ratio (g:d) as 3.5
2.3 (D/0.2)2(
/27
4)0.37(g:d/150)1
.
This value is at the lower end of the mass range of
Odenwald's (1988) cometary clouds of category II (3-40
).
In order to find possible candidates for the shaping of the
nebula, we used SIMBAD to locate discrete sources. There are 272 entries for a 4.5
4.5
region: 178 stars, 29 IRAS sources and 65 galaxies. Figure 4 shows their distribution as well as a sketch of the nebula. In general, all
appear to be rather randomly distributed. There were 8 IRAS sources found in the region of the nebula, but no IRAS source is found at the maxima (nuclei) in the contour maps.
![]() |
Figure 4:
Known discrete sources in the field of the nebula. a) stars (*), b) galaxies (o) and IRAS point sources (![]() ![]() ![]() ![]() |
Open with DEXTER |
The structure of the SkN suggests that - as one of the possibilities - a wind might blow from the north-west, leading to the cone-like morphology. To this end, we searched for any (stellar) objects which might be capable of shaping an interstellar cloud of small or moderate mass; such an object should be located at or around the head of the nebula and along the extrapolated major axis, i.e. towards the north or north-west. With SIMBAD, we found three possible candidates, namely an M-type emission star (RX J082605.8+262740), a carbon dwarf star (WD 0824+288) and a pulsar (PSR B0823+26). These were selected because (1) all of them are rather nearby, (2) they might emit a wind in the course of their evolution, (3) they show some peculiar properties, and (4) they are placed at a suitable apparent location with respect to the nebula. Their positions (and proper motions, if known) are shown in Fig. 4c.
RX J082605.8+262740 (RIXOS 206-517; RA = 08
,
Dec = +26
(J2000)) is an M-type emission star showing X-ray features. The source was detected in the RIXOS survey (Mittaz et al. 1999), which
registered objects down to a flux limit of 3
10-14 erg cm-2 s-1 (0.5-2 keV). The proper motion (
)
of this star is unknown. Its position looks promising (Fig. 4c) for shaping the nebula. We roughly estimated the
distance of the star: its brightness in B and V is 16.5 and 15.34,
respectively (Mason et al. 2000). Assuming negligible
extinction we found, as a rough guess, D = 17 pc. At this
distance the total mass and the physical size of the nebula would
become 0.025
0.017
(
= 27
4 K, g:d = 150)
and 0.05 pc, respectively. This would make the SkN a very small
cloudlet within the Local Bubble.
WD 0824+288 (RX J0827.0+2844; RA = 08
,
Dec = +28
(J2000)) is a dwarf
carbon (dC) star positioned at
2
north from the
northern nucleus (N1) of the nebula (Fig. 4c). The object dC
is hot (
52 kK), dim (V = 14.22) and has an average mass of
0.70
(Green et al. 2000; Green 2000;
Napiwotzki et al. 1999) and proper motion
(
= -28.2
2.7 mas yr-1,
= 0
0.9 mas yr-1) (Mauron et al.
2004, private communication). Heber et al. (1993)
noticed that this star shows variability due to heating effects
from a hot DA companion - the reason why this dwarf is visible
in the optical spectrum. In general, dwarfs are unlikely to
possess active chromospheres or undergo mass loss, but WD 0824+289
is a special case, where heating of the dwarf by a hot DA
companion causes the emission (Downes et al. 2004). This
dwarf is called a smoking gun for the existence of dust (Green et al.
2000). We roughly estimated the distance to this dwarf
using a value of the absolute magnitude estimated by Green
(2000) for all dCs, and found about 70 pc. Then, the
total mass of the nebula using the same values for the
dust-to-mass ratio and dust colour temperature as mentioned above
is 0.43
0.29
.
At D = 70 pc, the dC star would be
about 2.3 pc away from the nucleus N1 - an uncomfortably
large physical separation for shaping the SkN.
PSR B0823+26 (PSR J0826+2637; RA = 08
,
Dec = +26
(J2000)) is a rather
well studied pulsar with a proper motion (
= 61
3 mas yr-1,
= -90
2 mas yr-1) whose direction practically coincides with the long axis of the nebula
(Fig. 4c). This pulsar has a very nearly equatorial geometry, in
that both magnetic poles and the sightline all fall close to the
rotational equator of the star (Rankin & Rathnasree
1995). The transverse velocity, spin-down age and
distance of PSR B0823+26 are 194
34 km s-1, 4.92 Myr and
0.36 kpc, respectively (Gwinn et al. 1986). Very recently,
the X-ray emission properties of this pulsar was studied with the
XMM-Newton and it was found that the emission is largely dominated
by nonthermal processes (Becker et al. 2004).
The facts that both the position of PSR B0823+26 is at the head of the Skeleton Nebula and that the proper motion direction of the former coincides with the long axis of the latter is remarkable. Are these a chance superpositions? Unfortunately there is no convincing way to find out this, since a lot of structure in the IRAS maps exists at all scales, and particularly there is a multitude of variously shaped clouds, blobs, arcs, filaments, holes, etc. Nevertheless, one of us (BA) has started a study to address this problem of superpositions on IRAS maps.
To roughly estimate the probability to find such a superposition
we surveyed the IRAS dust structures around all nearby (
1 kpc) pulsars using the ATNF Pulsar database (http://www.atnf.csiro.au/research/pulsar/). To this end, we used
SkyView, selected the image scaling "Histogram
equalization'' (which permits to bring out very faint structures),
used the colour table "B-W Linear'' and selected
1
1
,
5
5
and
10
10
fields, at 12, 25, 60 and 100
m,
centered on each of the 70 nearby pulsars. We looked for all
isolated emission regions of filamentary, elliptical, and circular
shape. We found 5 cases in which pulsars could perhaps be
associated with extended dust structures, namely PSR J2124-3358,
PSR J2307+2225, PSR J1730-2304, PSR J0108-1431 and PSR J1909-3744.
J2124-3358 has an associated H-alpha nebula (Gaensler et al.
2002). However, none of these dust features shows
anything special (as to direction, location of the pulsar near the
centre of dust cavities and/or dust enhancements etc.)
Particularly, no obvious superposition with an isolated nebula,
similar to the case of the SkN, could be found. Hence, we suspect
that the association of PSR B0823+26 with the SkN might be real -
or we are misguided by a quite improbable chance superposition.
Having this strange result in mind we also looked whether there is
any other pulsar in the projected neighborhood of the SkN. The
closest (in projection) is PSR J0754+3231, which is
9
north-west of PSR B0823+26. Its distance is 3.92 kpc (Taylor et al. 1993), 11 times farther away than PSR B0823+26.
Is there anything that makes PSR B0823+26 a special case among
pulsars? Clegg et al. (1993) found that its 1400 MHz
dynamic spectrum exhibits a sharp discontinuity. The changes
occurred over time scales of 2 min. They believe that
such a discontinuity is the result of strong refraction in the
interstellar medium. We may speculate: is perhaps a part of the
SkN acting as a lens for this refraction, that is, is our nebula
at
360 pc? To get an idea of the dimensions involved: Clegg
et al. (1993) explained physical characteristics of the
lens system and derived a minimum lens size of 0.1 AU assuming
that the lens is located halfway to the pulsar. It is estimated
that a lens of size
0.7 AU located halfway (180 pc) to the
pulsar would have a diffractive time scale for fluctuation of
approximately 2 min.
If we assume that the SkN is at about the same distance as PSR B0823+26 (i.e., 360 pc), the total mass of the nebula would become 11.3
7.6
(
= 27
4 K, g:d = 150) and its maximum physical size (if in the sky plane) would amount to 16.8 pc. Anyway, despite the remarkable location and proper motion direction of the pulsar and the morphology of the nebula which led us to suspect that PSR B0823+26 might possibly shape the SkN, we
have to admit that i) up to now, no interaction of this kind - namely eroding an interstellar cloud by an approaching or closely passing pulsar - has been found; ii) we would at least expect
distinctly different physical conditions (heating) at the part of
the cloud which is closer to the pulsar - but nothing like this can be seen in the SkN; iii) it is hard to imagine how such a shaping by a neutron star could take place in principle. We intend
to study shaping mechanisms in the future, where also the possible
contribution of a pulsar nebula will be taken into account.
The Skeleton Nebula could be shaped via interaction with the
ambient medium. Oort (1969, 1970) has already
suggested that structures of clouds could be understood by
studying their collision with the ambient ISM. During a collision,
shock fronts will form that convert some of the kinetic energy of
the cloud into thermal energy, heating the gas possibly to
106 K. In hydrodynamic theory, the Reynolds number (R) is of relevance here, providing a dimensionless indication of whether the passes of a body through a fluid will result in
turbulent motion. Using Eq. (1) given by Odenwald & Rickard
(1987), we roughly estimate the R value for the SkN
assuming that it is falling towards the Galactic plane. We
emphasize that the location of the main axis of the SkN is almost
parallel to the Galactic plane (see Fig. 4), making such an interaction not very plausible.
Let us assume a distance of 200 pc (as Odenwald 1988
did for his category II clouds) and estimate the R value. The
minimum value of R is found to be 9 (assuming a coefficient of viscosity for fully ionized hydrogen gas as provided by Spitzer (1962) for T > 104 K, and a grain size, density and emissivity as suggested by Hildebrand
(1983) for 100
m) for a terminal velocity (=radial velocity/sin b) of 5 km s-1. The R value goes up to 90 if the terminal velocity of the nebula would be 10 times higher. This range of R (i.e., 9-90) suggests that the body
(i.e. our nebula) might be in the transition phase between a laminar and a turbulent flow. At this stage one can assume that vortices generate faster than that of the dissipation via
molecular diffusion.
The situation will be different if our nebula is closely
associated with the emission star RX J082605.8+262740. In this
case, R < 1 for a terminal velocity <100 km s-1 (radial
velocity 54 km s-1 at 32
latitude), indicating a laminar flow. The R value turns out to be >100 if the nebula is at the same distance as the pulsar PSR B0823+26. In this case, turbulence develops in the downstream.
A new infrared nebula (RA = 08
,
Dec =
+25
(J2000)) with an apparent size of
140
was detected on 60
m and 100
m IRAS maps. We have carried out far-IR and optical studies of this spectacular cone-like nebula. We found three discrete sources which might perhaps be responsible for shaping
this nebula: the M-type emission star RX J082605.8+262740, the
carbon dwarf star WD 0824+288 and the pulsar PSR B0823+26. Among
them, the pulsar PSR B0823+26 appears to be the prime candidate
because of its position, proper motion, distance and the presence
of a discontinuity in 1400 MHz dynamic spectrum, which suggests
refraction by interstellar or ambient matter. However, how could a pulsar erode an interstellar cloud? We cannot offer any explanation and hence assume that another shaping mechanism is at
work. A possible interaction of the nebula with the ambient ISM is discussed in terms of Reynolds number hydrodynamics. This discussion suggests that the nebula might be in a transition phase
between laminar and turbulent motion. Anyway, the process or
source responsible for shaping this nebula has not been found yet.
Acknowledgements
We acknowledge the unknown referees for their constructive criticism and numerous useful comments which have allowed to improve the paper. We thank the Austrian Science Fund (Fonds zur wissenschaftlichen Forschung) for support of the project P15316. B. Aryal acknowledges Tribhuvan University, Kathmandu, Nepal for granting study leave.