A&A 462, 677-682 (2007)
DOI: 10.1051/0004-6361:20066387
D. Stamatellos1 - D. Ward-Thompson1,2 - A. P. Whitworth1 - S. Bontemps2
1 - School of Physics & Astronomy, Cardiff University, 5 The Parade, Cardiff CF24 3AA, Wales, UK
2 - Observatoire de Bordeaux, BP 89, 2 Rue de l'Observatoire, 33270 Floirac, France
Received 12 September 2006 / Accepted 6 November 2006
Abstract
Aims. Four dense cores, L1582A, L1689A, B133 and B68, classified as prestellar in terms of the absence of detectable NIR emission, are observed at radio wavelengths to investigate whether they nurture very young protostars.
Methods. We perform deep radio continuum observations at 3.6 cm and 6 cm using the VLA.
Results. No definite young protostars were discovered in any of the four cores observed. A few radio sources were discovered close to the observed cores, but these are most likely extragalactic sources or YSOs unrelated to the cores observed.
In L1582A we discovered a weak radio source near the centre of the core with radio characteristics and offset from the peak of the submillimeter emission similar to that of the newly discovered protostar in the core L1014, indicating a possible protostellar nature for this source. This needs to be confirmed with near- and/or mid-infrared observations (e.g. with Spitzer). Hence based on the current observations we are unable to confirm unequivocally that L1582A is starless.
In L1689A a possible 4.5-
radio source was discovered at the centre of the core, but needs to be confirmed with future observations.
In B133 a weak radio source, possibly a protostar, was discovered at the edge of the core on a local peak of the core submm emission, but no source was detected at the centre of the core. Thus, B133 is probably starless, but may have a protostar at its edge.
In B68 no radio sources were discovered inside or at the edge of the core, and thus B68 is indeed starless.
Four more radio sources with spectral indices characteristic of young protostars were discovered outside the cores but within the extended clouds in which these cores reside.
Conclusions. We conclude that the number of cores misclassified as prestellar is probably very small and does not significantly alter the estimated lifetime of the prestellar phase.
Key words: stars: formation - radio continuum: stars - ISM: clouds
The identification of the youngest protostars is important for understanding how dense cores in molecular clouds collapse and form stars. The initial conditions for star formation are provided by dense cores in molecular clouds in which there is no evidence that star formation has occurred (starless cores, Myers et al. 1983). Some of these cores are thought to be close to collapse or already collapsing, and they are labelled prestellar cores (e.g. Ward-Thompson et al. 1994; Ward-Thompson et al. 2002). Prestellar cores evolve to Class 0 objects, which represent the youngest protostars (André et al. 1993).
The presence of a protostar in a core is inferred by one or more of the following criteria (André et al. 2000): (i) compact centimeter radio emission; (ii) a bipolar molecular outflow; (iii) NIR or MIR emission. These criteria depend on the sensitivity of the telescope used, and how young and embedded the protostar is in its parent cloud. Additionally, the amount of radio emission from young protostars is uncertain and probably very low (e.g. Curiel et al. 1987; Rodríguez et al. 1989b; Neufeld & Hollenbach 1996; Harvey et al. 2002; Eiroa et al. 2005). Finally, in many cases the complexity of star forming regions (e.g. dense environment, multiple sources, inflow of material) does not allow a clear identification of all the outflows in a given region. Thus, a null detection with a particular telescope does not necessarily mean that a source does not exist inside a core.
The need for very sensitive observations to discover young protostars residing inside dense cores was demonstrated recently by Young et al. (2004), who, using the SPITZER space telescope, detected a NIR source embedded in L1014 (L1014-IRS), a dense core that was previously classified as starless. This suggests that L1014 is in fact a young Class 0 object. Since the Young et al. discovery, the protostellar nature of L1014-IRS has been confirmed by the detection of a very weak bipolar outflow emanating from the IR source (Bourke et al. 2005) that affects the geometry of the observed scattering light nebula around the source (Huard et al. 2006). Additionally 2 more NIR sources have been discovered inside cores previously thought to be starless: in L1521F in the Taurus molecular cloud (Terebey et al. 2005; Bourke et al. 2006), and in L1148 in the Cepheus flare (Kauffmann et al. 2005). These recently discovered low-luminosity protostars are similar to the previously known Class 0 protostars VLA1623 (André et al. 1993) and HH24MMS (Ward-Thompson et al. 1995; Bontemps et al. 1996). These discoveries of low-luminosity, young protostars, open up the possibility that other cores previously classified as prestellar may in fact contain very young protostars.
Radio observations are suitable for detecting deeply embedded protostars, as long wavelength radiation can escape the large column density of the envelope practically unattenuated. Previous studies (e.g. André et al. 1987; Leous et al. 1991; Bontemps et al. 1995) have identified a large number of young protostars exhibiting detectable radio emission. However, very young and low luminosity protostars are believed to produce very weak radio emission (e.g. Harvey et al. 2002; Neufeld & Hollenbach 1994, 1996), and thus deep observations are needed to detect them.
Young protostars exhibit a positive or almost flat (>-0.1)
spectral index (defined as ,
where
)
between e.g. 3.6 cm and 6 cm
(André et al. 1987; Anglada et al. 1998; Beltrán et al. 2001). This
behaviour of the radio spectral index is characteristic of free-free
thermal emission from ionised gas. This is most likely produced by a
partially ionised jet that propagates into the collapsing envelope
producing shocks (e.g. Curiel et al. 1987). It could also be produced by
an ionised disc wind (e.g. Martin 1996) or by the accretion shock that
develops on the surface of the protostar, and/or on the surface of the
disc that surrounds the protostar, which heats and ionises the infalling
gas (Winkler & Newman 1980; Cassen & Moosman 1981). However, the latter
two mechanisms produce radio emission that is at too low a level to be
readily detected.
A negative (<-0.1) radio spectral index indicative of non-thermal gyro-synchrotron emission has also been observed in some T Tauri stars (e.g. André et al. 1988; André 1996; Rodríguez et al. 1999), in a possible Class I object (Feigelson et al. 1998), and even in a possible Class 0 object, the triple radio continuum source in Serpens (Rodríguez et al. 1989a; Curiel et al. 1993; Curiel 1995; Raga et al. 2000). The non-thermal emission of more evolved protostars (Class II, III) is believed to be associated with the magnetic fields around these objects and the star-disc interaction region (e.g. André 1996). The non-thermal emission from younger protostars is thought to be the result of particle acceleration behind a diffuse shock wave (Rodríguez et al. 1989a; Curiel et al. 1993, and references therein).
Stamatellos et al. (2005a,b) performed 3-dimensional radiative transfer
calculations on simulations of the collapse of a star-forming core and
studied the transition from a starless core to a core with an embedded
young protostar. These simulations suggest (i) that cores with young
protostars appear to be cold (
K) but still hotter than
starless cores (T<15 K), and (ii) that cores with young protostars appear
more circular (at least in their central regions) than starless cores (see
also Goodwin et al. 2002). Based on this study, we identified 4 cores
classified as starless (L1582A, L1689A, B133, B68) that might
in fact contain very young protostars. These cores were selected in terms
of their circular appearance and their relatively high estimated
temperatures (
15 K). These
cores were observed with the VLA in the continuum at 3.6 cm and 6 cm.
In this paper we present the results of these radio observations aiming to detect very young embedded sources in dense cores. In Sect. 2 we describe the observational details and the properties of the detected radio sources, in Sect. 3 we summarize our results, and in Sect. 4 we report our conclusion.
L1582A, L1689A, B133 and B68 were observed at 3.6 cm using the VLA. The
observations were performed in July/August 2005, in the VLA C
configuration. The total integration time was about 4 h for
each core including time spent on the calibrators. These cores (except B68) were also observed at 6 cm in October/November 2005, in the DnC
configuration. L1582A was observed for
5 h, B133 for
5 h, and L1689A for
1 h. The data were calibrated using the
AIPS package of the NRAO.
Table 1: Distances and positions of the dense cores observed.
Table 2: Phase calibrators, beam sizes, rms noise, and number of radio sources detected.
Table 3: Radio sources detected at 3.6 and/or 6 cm.
In Table 1 we list the positions of the cores taken from
Kirk et al. (2005). These positions correspond to the peak of the
850
and 450
emission detected using SCUBA with a beam
size FWHM of 14.8
.
We use these positions as phase centres.
In Table 2 we list the phase calibrators used, their
bootstrapped flux densities, the VLA beam size, and the achieved rms noise for each field.
Finally, we report the number of radio sources detected and the number of
background, extra-galactic sources expected in each field (see below).
The amplitude calibrators used were
3C 286 (
Jy,
Jy)
and 3C48 (
Jy,
Jy).
The number of background sources expected in our field is determined using
the Anglada et al. (1998) calculation. Assuming that extra-galactic sources have the characteristics
described by Condon (1984),
the expected number of background sources
in a field of angular diameter
is
![]() |
(1) |
In Table 3 we list the detected radio sources, their
positions, flux densities and peak fluxes,
their spectral indices
(in cases where a source has
been detected both in 3.6 cm and 6 cm), and their offsets from the peak of
the submm emission of the core.
The fluxes were corrected for the effect of primary-beam
response. We assume a detection limit of 5-
(however we make reference to a possible 4.5-
detection).
Most of the sources do not correspond to
any sources listed in the SIMBAD astronomical database and are new
sources.
In the next subsections we describe in detail the observed radio sources and their characteristics, for each of the cores observed.
In L1582A we discovered a weak radio source (VLA1) close
(
)
to the peak of the submillimeter core emission. The
source is observed both in 3.6 cm and 6 cm
(see Figs. 1, 2). The spectral index
of the source is
.
This spectral index is not what is
expected for Class 0 protostars. For example, for the prototypical Class 0
object VLA1623 (André et al. 1993)
. The spectral index of
the radio emission from protostars is thought to be between -0.1 and 2, consistent with thermal free-free emission, whereas
for extra-galactic sources it is expected to be <-0.1, characteristic of
non-thermal emission due to optically thin synchrotron emission (e.g.
André et al. 1987; Anglada et al. 1998; Beltran et al. 2001). Thus,
L1582A-VLA1 is most probably an extragalactic source.
However we note that recent observations of the confirmed low luminosity protostar L1014-IRS (Shirley et al., VLA proposal) indicate a similar spectral index. The unusual negative spectral index may be due to non-thermal (synchrotron) emission. The similarities of L1582A-VLA1 to L1014-IRS (similar 3.6 cm and 6 cm fluxes, similar offsets from the peak of the submillimeter emission of the core) are suggestive that L1582A-VLA1 may also be a very young protostar embedded in the core, but this cannot be confirmed with the current observations.
Two more sources were detected in the field. VLA2 was detected only at 3.6 cm but a negative spectral index cannot be excluded. This source is outside the core but within the extended CO (J=1-0) line emission (Dame et al. 2001; see also Kirk et al. 2005). VLA3 was detected only at 6 cm and thus has a negative spectral index and should be a background source.
Thus, we are unable to confirm unequivocally that L1582A is starless.
![]() |
Figure 1:
L1582A: 3.6 cm VLA observations (solid thick contours)
overlaid on 850
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2:
L1582A: 6 cm VLA observations (solid thick contours) overlaid on
850
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In the L1689A field we discovered 6 sources, 5 of which were detected both at 3.6 cm and 6 cm, and 1 marginally detected only at 3.6 cm.
VLA1 has a positive spectral index, indicative of free-free thermal emission from ionised gas. It is outside the main core but within the extended CO line emission of the Ophiuchus GMC (Dame et al. 2001; Kirk et al. 2005). It coincides with a local column density peak of the cloud (see Fig. 3; and Nutter et al. 2006), although the dust emission is rather weak.
VLA2 is also outside the core but within the extended CO emission. The sign of its spectral index cannot be determined; it could be either a young star or a background source.
Sources VLA3, VLA4 and VLA5 were detected in both wavelengths and they
seem to be related, as they are aligned. They are located outside the
L1689A core but they are coincident with a submm filament adjacent to the
core. The total projected extent of the structure is
(i.e.
6000 AU if they are physically located close to L1689A). The
central (VLA3) and the southern (VLA5) source have negative spectral
indices (characteristic of non-thermal synchrotron emission),
whereas the northern source (VLA4) may have positive or marginally negative
spectral index (characteristic of optically thin thermal emission).
Thus, these sources are probably extragalactic.
![]() |
Figure 3:
L1689A: 3.6 cm VLA observations (solid thick contours)
overlaid on 850
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
L1689A:
6 cm VLA observations (solid thick contours)
overlaid on 850
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The last source in this field, VLA6, was marginally detected at 3.6 cm
(4.5-
detection). It is located at the centre of the core. It was
not detected at 6 cm but a negative spectral index cannot be ruled out.
The presence of this source in the centre of the core, if indeed it is
real, would explain the temperature increase toward the centre of the core
suggested by the observations of Ward-Thompson et al. (2002). Thus, there
is the possibility that L1689A may not be starless.
In the B133 field we discovered 4 radio sources. The strongest one,
B133-VLA1, is located outside the core. It has a negative spectral index
between 3.6 and 6 cm, indicating non-thermal emission due to optically
thin synchrotron emission characteristic from an extra-galactic source. This
source was the only one from this field that was detected at 6 cm. A
search in the SIMBAD database reveals a match with PMNJ1906-0650 of the
Parkes-MIT-NRAO (PMN) Surveys (Griffith et al. 1995) at 4.85 GHz
(
mJy). The source has been detected by the NRAO VLA Sky
Survey (Condon et al. 1998) at 1.4 GHz (
mJy), and by
the Texas Survey of Discrete Radio Sources (Douglas et al. 1996) at 365 MHz (
mJy). These fluxes are also consistent with thin
synchrotron emission from an extra-galactic source.
![]() |
Figure 5:
B133:
3.6 cm VLA observations (solid contours)
overlaid on 850
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The second strongest source in the field, B133-VLA2, is located outside
the core and it is not associated with any column density peak. The
absence of detectable 6 cm emission at the 300 Jy level indicates a
positive spectral index (
). This suggests free-free
thermal emission. B133-VLA2 may be
associated with the Aquila Rift, where B133 is thought to reside
(Kirk et al. 2005).
The B133-VLA3 source is located at the edge of the core within a local column density
peak. It has been marginally detected (at the 6-
level) and
it could be a very young protostar. Unfortunately,
although it is not detected at 6 cm a negative spectral index cannot be
excluded, so it may also be an extragalactic source.
Finally, B133-VLA4 is marginally detected at 3.6 cm. It is not detected at 6 cm but a negative spectral index is possible.
Thus, based on the VLA observations B133 appears to be a starless core with a possible protostar (B133-VLA3) at its edge.
Four radio sources were discovered in the B68 field at 3.6 cm. All of them are outside the core and they are not associated with the core but they are within the extended CO emission (Dame et al. 2001; Kirk et al. 2005) around the core. This core was not observed at 6 cm, and thus the spectral indices of these sources and their nature cannot be determined. We conclude that B68 is starless.
We conducted continuum VLA radio observations of 4 dense cores previously classified as prestellar, searching for radio emission from deeply embedded protostars. 17 radio sources were detected in the observed fields, i.e. more than the expected number of background sources. Thus at least a few of them may be associated with the star forming regions observed.
Two radio sources were discovered near the centres of the cores.
L1582A-VLA1 is located within 18
of the peak of the submm emission
of the core, and its similarity with the source embedded in L1014
suggests that this may also be a very young protostar embedded in the core,
despite its negative spectral index between 3.6 and 6 cm, and its offset from the
centre of the core.
L1689A-VLA6 is located within 1
of the peak of the submm emission,
and it may also be a very young protostar embedded in the core.
However, this detection was only at the 4.5-
level.
Future Spitzer data should confirm whether these are young protostars or not.
If these two sources indeed are associated with the observed cores then they are similar to
the recently discovered low-luminosity protostars, such as that in L1521F, which
are young Class 0 objects, e.g. like VLA1623 (André et al. 1993)
and HH24MMS (Ward-Thompson et al. 1995; Bontemps et al. 1996).
However, based on the current observations and the locations of the
radio sources, we conclude that L1582A and L1689A are most likely starless.
Two radio sources were discovered close to edges of cores. L1689A-VLA345, at the edge of the L1689A, is probably an extragalactic source. B133-VLA3, at the edge of B133, coincides with a local peak in the core submillimeter emission, and it may be a young protostar.
Four radio sources with positive spectral indices were detected outside the cores (L1689A-VLA1, L1689A-VLA2, B133-VLA2, B133-VLA4) and, thus could also be young protostars or stars associated with the extended star forming regions where these cores reside.
All of the candidate protostars reported are expected to be detectable in the NIR by the Spitzer space telescope and such observations are needed to confirm whether these sources are indeed protostars. Additionally, sensitive molecular line radio observations, with e.g. the Submillimeter Array (SMA) should be able to detect molecular outflows emanating from these objects, if they are protostars.
In this paper we presented deep radio observations of 4 prestellar cores, which are relatively warm and almost circular (in their inner regions). According to the Stamatellos et al. (2005) criteria these cores are, among cores presently classified as prestellar, the ones most likely to nurture very young protostars. However, we find no definite protostars at the centres of these cores. Two radio sources that could be very young embedded protostars were discovered near the centres of L1582A and L1689A, but further near- or mid-infrared observations (e.g. with Spitzer or the upcoming Herschel) are needed for safe conclusions on the nature of these sources. Hence, based on the current observations, similarly to Kirk et al. (2006) we conclude that the number of cores misclassified as prestellar is probably very small and does not affect significantly the estimated lifetime of the prestellar phase.
Acknowledgements
We thank C. Chandler for her help in reducing and calibrating the VLA data at NRAO Sorocco, and J. Kirk and D. Nutter for providing the SCUBA data. The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the Particle Physics and Astronomy Research Council of the United Kingdom, the Netherlands Organisation for Scientific Research, and the National Research Council of Canada. The JCMT data shown in this paper were taken during observing runs M96BU56, M97BU88, M99AU34 and M00AU32. This work was carried out while DWT was on sabbatical at l'Observatoire de Bordeaux and CEA, Saclay, and he wishes to thank both institutions for the hospitality accorded to him. We also acknowledge support by PPARC grant PPA/G/O/2002/00497.