A&A 462, 677-682 (2007)
DOI: 10.1051/0004-6361:20066387

A VLA search for young protostars embedded in dense cores

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-$\sigma $ 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

1 Introduction

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 $\alpha$, where $S_\nu\propto \nu^a$) 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 ( $T\sim15{-}30$ 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 ($\sim$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.

2 Observations and results

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 $\sim$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 $\sim$5 h, B133 for $\sim$5 h, and L1689A for $\sim$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  $\hbox {$\mu $ m}$ and 450  $\hbox {$\mu $ m}$ emission detected using SCUBA with a beam size FWHM of 14.8 $^{\prime \prime }$. 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 ( $F_{\rm 3.6~ cm}=5.26$ Jy, $F_{\rm 6~ cm}=7.34$ Jy) and 3C48 ( $F_{\rm 3.6 ~cm}=3.25$ Jy, $F_{\rm 6~ cm}=5.32$ 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 $\langle N\rangle$ in a field of angular diameter $\theta_{\rm F}$ is

\begin{displaymath}\langle N\rangle=1.4\ \left[1-{\rm e}^{-0.46\
\left(\frac{S_...
...ight)^2}\right]
\left(\frac{\nu}{5 ~{\rm GHz}}\right)^{-2.52},
\end{displaymath} (1)

where S0 is the minimum detectable flux density at the centre of the field. Each of the fields has an angular diameter of $\theta_{\rm F}\approx500\hbox{$^{\prime\prime}$ }$, so the expected number of background sources at 3.6 cm is

 \begin{displaymath}
\langle N\rangle_{\rm 3.6~cm}=2.6\ \left(\frac{S_0}{\rm 50\ \mu Jy}\right)^{-0.75} ,
\end{displaymath} (2)

and at 6 cm,

 \begin{displaymath}
\langle N\rangle_{\rm 6~cm}=3.0\ \left(\frac{S_0}{\rm 100\ \mu Jy}\right)^{-0.75} \cdot
\end{displaymath} (3)

Thus, we expect a total number of $\sim$$9\pm3$ background sources in all our fields at 3.6 cm, and $\sim$$5\pm2$ sources at 6 cm. We detected 16 sources at 3.6 cm and 8 sources at 6 cm, thus the expected number of Galactic sources detected, possibly associated with the star forming regions observed, is $7\pm3$.

In Table 3 we list the detected radio sources, their positions, flux densities and peak fluxes, their spectral indices $\alpha$ (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-$\sigma $ (however we make reference to a possible 4.5-$\sigma $ 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.

2.1 L1582A

In L1582A we discovered a weak radio source (VLA1) close ($\sim$ $18\hbox{$^{\prime\prime}$ }$) 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 $\alpha=-1.2\pm 0.5$. 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) $\alpha=0.6$ . 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.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{6387fig1.ps}
\end{figure} Figure 1: L1582A: 3.6 cm VLA observations (solid thick contours) overlaid on 850  $\hbox {$\mu $ m}$ SCUBA observations from Kirk et al. (2005) (grayscale and solid thin contours). 3.6 cm solid contours are plotted at 52 $\mu $Jy(5-$\sigma $) $\times $ (1, 1.5, 2, 2.5, 3). 850 mm dashed contours are plotted at 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 0.95 of the maximum 850 mm flux (170 mJy/beam). The VLA beam size is mentioned in Table 2. The SCUBA beam size at 850  $\hbox {$\mu $ m}$ is 14 $^{\prime \prime }$. A radio source has been discovered near the centre of the core.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{6387fig2.ps}
\end{figure} Figure 2: L1582A: 6 cm VLA observations (solid thick contours) overlaid on 850  $\hbox {$\mu $ m}$ SCUBA observations from Kirk et al. (2005) (grayscale and solid thin contours). 6 cm solid contours are plotted at 145 $\mu $Jy(5-$\sigma $$\times $ (1, 1.2, 1.4, 1.6, 1.8). 850 mm dashed contours are plotted as in Fig. 1. The VLA beam size is mentioned in Table 2. The SCUBA beam size at 850  $\hbox {$\mu $ m}$ is 14 $^{\prime \prime }$.
Open with DEXTER

2.2 L1689A

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 $\sim$ $ 46\hbox{$^{\prime\prime}$ }$(i.e. $\sim$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.


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{6387fig3.ps}
\end{figure} Figure 3: L1689A: 3.6 cm VLA observations (solid thick contours) overlaid on 850  $\hbox {$\mu $ m}$ SCUBA observations from Nutter et al. (2006) (grayscale and thin solid contours). 3.6 cm solid contours are plotted at 54 $\mu $Jy (4.5-$\sigma $) $\times $ (1, 2, 4, 6, 8, 10, 12). 850 mm dashed contours are plotted at 0.2, 0.4, 0.6, and 0.8 of the maximum 850 mm flux (170 mJy/beam).The VLA beam size is mentioned in Table 2. The SCUBA beam size at 850  $\hbox {$\mu $ m}$ is 14 $^{\prime \prime }$.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{6387fig4.ps}
\end{figure} Figure 4: L1689A: 6 cm VLA observations (solid thick contours) overlaid on 850  $\hbox {$\mu $ m}$ SCUBA observations from Nutter et al. (2006) (grayscale and thin solid contours). 6 cm solid contours are plotted at 300 $\mu $Jy(5-$\sigma $) $\times $ (1, 1.5, 2, 2.5, 3, 6, 9). 850 mm dashed contours are plotted as in Fig. 3. The VLA beam size is mentioned in Table 2. The SCUBA beam size at 850  $\hbox {$\mu $ m}$ is 14 $^{\prime \prime }$.
Open with DEXTER

The last source in this field, VLA6, was marginally detected at 3.6 cm (4.5-$\sigma $ 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.

2.3 B133

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 ( $S_\nu=65\pm11$ mJy). The source has been detected by the NRAO VLA Sky Survey (Condon et al. 1998) at 1.4 GHz ( $S_\nu=295.5\pm8.9$ mJy), and by the Texas Survey of Discrete Radio Sources (Douglas et al. 1996) at 365 MHz ( $S_\nu=1216\pm41$ mJy). These fluxes are also consistent with thin synchrotron emission from an extra-galactic source.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{6387fig5.ps}
\end{figure} Figure 5: B133: 3.6 cm VLA observations (solid contours) overlaid on 850  $\hbox {$\mu $ m}$ SCUBA observations from Kirk et al. (2006) (grayscale and dashed contours). 3.6 cm solid contours are plotted at 60 $\mu $Jy (5-$\sigma $) $\times $ (1, 1.1, 1.2, 1.3). 850 mm dashed contours are plotted at 0.2, 0.4, 0.6, and 0.8 of the maximum 850 mm flux (170 mJy/beam). The VLA beam size is mentioned in Table 2. The SCUBA beam size is 14 $^{\prime \prime }$.
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 $\mu $Jy level indicates a positive spectral index ( $\alpha>1.7\pm0.2$). 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 $\sim$6-$\sigma $ 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.

2.4 B68

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.

3 Summary

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 $^{\prime \prime }$ 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 $^{\prime \prime }$ 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-$\sigma $ 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.

4 Conclusion

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.

References

 

Copyright ESO 2007