Free Access
Issue
A&A
Volume 542, June 2012
Article Number A14
Number of page(s) 7
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201118710
Published online 25 May 2012

© ESO, 2012

1. Introduction

Spectropolarimetric observations of masers are a powerful tool to study the magnetic field properties in the maser pumping zone. Water masers are unique because they are found in a variety of astrophysical environments and, in particular, star-forming regions in distinct mass regimes. In contrast, methanol and OH masers are mostly found to be associated with high-mass star formation sites (e.g., Vlemmings 2008; Sanna et al. 2010a,b). The most commonly observed water-maser line is the (616 − 523) transition at 22 GHz, which is an excellent probe of molecular gas at very high volume densities (nH2 in the 108 to 1010 cm-3 range, Elitzur et al. 1989). While the Zeeman splitting is small (~10-3 Hz μG-1), the line strength is sometimes sufficient for measuring the circular polarization and therefore performing a direct measurement of the line-of-sight (LOS) component of the magnetic field. If linear polarization is also measured, then the full three-dimensional (3-D) magnetic field configuration can be derived.

There have been several recent studies of the efficiency of the maser spectropolarimetry technique. These have analyzed the magnetic field properties in very dense molecular environments around the circumstellar envelopes of evolved stars and young stars (e.g., Surcis et al. 2011a,b; Pérez-Sánchez et al. 2011). A representative illustration of how magnetic fields can be resolved on very small spatial scales was obtained by Vlemmings et al. (2006), who carried out Very Long Baseline Array (VLBA) polarization observations of H2O masers around the Cepheus A HW2 high-mass object.

Thermal dust polarization emission is difficult to observe in the most embedded portions of molecular clouds due to the low instrumental sensitivity to polarized data (and the feasibility is then limited toward a handful of objects). Submillimeter (submm) emission suffers from depolarization effects such as unresolved fields or the roundness of dust grains (Goodman et al. 1992, 1995; Lazarian et al. 1997). Water masers can be thought of as an observational tool to overcome this issue, since they are excited at very high densities.

IRAS 16293-2422 (hereafter, I16293) is a prototypical low-mass protostellar system located in the ρ Ophiuchus molecular cloud (d ≃ 120 pc, Loinard et al. 2008). This source is a well-studied binary system, usually referred to as sources A and B, that has an angular separation of 5′′ (600 AU) and is very embedded in a dense molecular core (Wootten 1989; Looney et al. 2000). Both sources have a very rich chemistry that is typical of that found in hot cores (Ceccarelli et al. 2000; Kuan et al. 2004; Bisschop et al. 2008; Jørgensen et al. 2011). Although component B appears to be a single source, high resolution interferometer data reveal that source A contains a higher degree of fragmentation. Two dust submm components were detected within this object with an angular separation of  ~0.6′′ (~72 AU, Chandler et al. 2005). I16293 has two large-scale (~1′) bipolar CO outflows, one of them associated with source A, while the powering source of the other CO outflow is a matter of debate (Walker et al. 1988; Stark et al. 2004; Yeh et al. 2008). Observations of the SiO (8–7) emission have revealed that there is also a compact molecular outflow associated with source A (Rao et al. 2009). I16293 has strong water maser emission that has been well-monitored (Wilking & Claussen 1987; Terebey et al. 1992; Claussen et al. 1996; Furuya et al. 2003). The strongest features appear at blueshifted and redshifted LSR velocities with respect to the ambient cloud velocity (~4 km s-1), typically with VLSR between  − 5 and 10 km s-1. The H2O maser emission very often has intensities of more than 100 Jy, and at times it is stronger than 300 Jy, at LSR velocities of  ~ 7 km s-1.

Tamura et al. (1993) performed observations of the 1.1 mm dust polarized emission toward I16293 at an angular resolution of 19′′. They found that the magnetic field lines are perpendicular to the major axis of the dense elongated disk-like molecular structure (~10′′ size,  ~1200 AU) observed in C18O by Mundy et al. (1990). This morphology was corroborated by Rao et al. (2009), who used the Submillimeter Array (SMA) and obtained a polarization map at an angular resolution of  ~2′′ (~240 AU), resolving the sources A and B. The mean volume density traced by the SMA maps is  ~6 × 107 cm-3. The polarization pattern around source A is compatible with a hourglass morphology for the magnetic field, whose strength was estimated to be  ~4.5 mG. The SMA maps show that there is a misalignment between the outflow direction and the magnetic field axis, that is roughly in agreement with model predictions where the magnetic energy is comparable to the centrifugal energy. In contrast, source B is associated with a uniform and apparently undisturbed magnetic field.

In this work, we report spectropolarimetric Very Large Array (VLA) observations of H2O masers toward I16293. In Sect. 2, we describe the details of the observational setup. In Sect. 3, the results obtained from the maser spectroscopy are shown. In Sect. 4, we discuss a possible correlation between the masers features and the outflows. The magnetic field strength estimation and its implications for the core dynamical evolution are discussed in Sects. 5 and 6, respectively. Finally, in Sect. 7 we summarize our conclusions.

2. Observations

Our observations were performed with the VLA (NRAO1, New Mexico, USA) in its extended A-configuration, on 2007 June 25 and 27. Each observing track lasted  ~5.5 h. A total of 27 antennas were used, 10 of which had already been retrofitted with the new system, resulting in a combined VLA/EVLA (extended VLA) observation. We used the K band receivers (22 − 24 GHz) tuned to the frequency of the water maser (616 − 523) rotational transition (ν0 = 22.23508 GHz). We took advantage of the full polarization capability of the correlator, selecting a bandwidth of 0.7813 MHz (~10.5 km s-1 in velocity). The spectral setup contains a total of 128 channels, which covers most of the velocity range of the strongest water maser features (observed around the brightest feature previously reported at  ~7 km s-1) at the spectral resolution of 0.08 km s-1. The quasar J1626-298 was used as the gain calibrator. The quasar J1331+305 and the radio source J1751+096 were used as polarization calibrators in order to correct both the instrumental feed polarization and the position angle of the polarization vectors. All three calibrators were also used for bandpass corrections. We performed self-calibration using the channel with the strongest intensity. Since the detected emission was unresolved, both phase and amplitude self-calibration solutions were applied to the other channels.

The combined VLA/EVLA setup introduces about 8% closure errors into the VLA-EVLA baselines-based corrections for the bandpass. This happens because of the distinctly different bandpass response of the EVLA antennas compared to the VLA ones. These errors become larger at narrow bandwidths, when the processing of the digital signal of EVLA antennas aliases the power response at the band edges. During the data reduction, we took special care to determine the bandpass solutions by performing the particular calibration strategies recommended to the user in the NRAO webpage2. They consisted of flagging EVLA-EVLA baselines prior to bandpass calibration and unflagging them prior to applying the bandpass solutions to the data set. The bandpass corrections are then applied to both an average single-channel multi-source dataset (“channel 0” file) and the multi-channel original file.

Data reduction was performed with the Astronomical Image Processing Software package (AIPS). Imaging of Stokes parameters I, Q, and U were generated with a quasi-uniform weighting (robust = –1, Briggs 1995). Maps of polarized fraction (P) and position of polarization angles (PA) were obtained by combining the Stokes Q and U images in such a way that and PA . The resulting synthesized beam is 0.14′′ × 0.08′′, with a position angle of  − 5.7°. The rms noise for channels where no emission was detected is  ~8 mJy beam-1 and increases to 23 mJy beam-1  at the peak emission channel. A slightly lower rms is observed for polarized intensity.

3. Results

3.1. H2O maser line

The contour channel maps of the H2O emission observed with the VLA/EVLA data are shown in Fig. 1. The emission extends over a wide range of velocities (4.5 <  VLSR  < 9 km s-1), i.e., at redshifted LSR velocities with respect to the cloud systemic velocity (~4 km s-1). Our spectral setup covers only a small portion of LSR velocities blueward of it (~2 <  VLSR  < 4 km s-1) but no emission was detected within this spectral range. The water maser line has a peak intensity of 170 Jy beam-1, which is detected at VLSR  ≃ 7.4 km s-1. This emission was also extensively reported in previous surveys as the brightest feature (see, for instance, Claussen et al. 1996; Furuya et al. 2003). Some channels have emission that is extended slightly eastward (VLSR between  ≃ 6.5 and 7 km s-1), suggesting that the observed emission in these channels is marginally resolved.

thumbnail Fig. 1

Channel maps of the H2O emission toward IRAS 16293-2422. Contours are  − 50, 50, 500, 3 × 103, 1 × 104, 2 × 104 times 8 mJy beam-1, the rms noise of the maps with lower intensities (see Sect. 2). The peak flux of 167 Jy beam-1 occurs at VLSR = 7.4 km s-1. The LSR velocity of each channel is labeled.

A scheme of the distribution of dust and molecular outflows in I16293 is shown in Fig. 2. The masers detected in our observations are associated with source A. They are located  ~0.25′′  (~30 AU in projection) to the southeast of the dust condensation Aa resolved by the subarcsecond submillimeter (submm) observations of Chandler et al. (2005).

The Stokes I spectrum of the maser emission is shown in Fig. 3 (upper panel). The non-Gaussian line profile indicates that there are unresolved components in the spectrum. Apart from the peak intensity at 7.4 km s-1, unresolved emission also seems to be present at lower velocities (VLSR  ≃  5.5 km s-1) with a strong flux of  ~20 Jy beam-1. Fainter emission (~5 Jy beam-1) is also observed at higher velocity channels (VLSR  ≃  9.2 km s-1). Therefore, there are at least three unresolved components (Table 1). A two-dimensional Gaussian profile fit to each of those components provides a mean spatial separation of  ~22 milli-arcsec (mas), which is larger than the precision in the relative position determination (rmspos ≃ 2 mas, as estimated from the ratio of the width of the synthesized beam to the signal-to-noise ratio of the fainter component).

3.2. Polarized emission

The spectrum of linearly polarized intensity is very similar to the Stokes I line profile except for the flux scale. It peaks at the same systemic velocity as the Stokes I spectrum, but is weaker by a factor of  ~30. The dependence of the linear polarization intensity, polarization fraction, position angle, and Stokes I parameter with systematic velocity is shown in Fig. 3. The measured polarization fraction is 2.5 ± 0.2%. The polarization position angle θ is  −23° and shows only small changes across the maser (σθ = 2°), implying that the polarization vectors at different velocities trace basically the same region. The formal uncertainty in PA (, Wardle & Kronberg 1974) is small and ranges between  ~0.14° and 0.80° considering channels with strong and weak linear polarization emission, respectively. The linear polarization as derived from the Stokes U and Q maps is exhibited in Fig. 4, which shows the distribution of polarization vectors in the brightest velocity channel. The peak of the polarized intensity is offset by 0.05′′ from the Stokes I peak. The polarization vectors may be either parallel or perpendicular to the magnetic field orientation in the plane of the sky. In Sect. 5, we discuss any assumptions that must be considered in order to solve this ambiguity.

Table 1

H2O maser components in IRAS 16293-2422.

The line profile of the circular polarization (Stokes V) has the characteristic S-shape (see Fig. 5, lower panel). The Stokes V spectrum is proportional to both the first derivative of the Stokes I spectrum and the LOS component of the magnetic field. The red dashed line in Fig. 5 indicates the fit that most closely represents this derivation. Owing to the remaining gain differences among the two polarizations, scaled-down replicas of the Stokes I spectrum might contaminate the Stokes V data (Sarma et al. 2001; Vlemmings et al. 2002). These features were considered in the fit and subtracted from the synthetic V spectrum. The fraction of circular polarization, calculated as (Vmax − Vmin)/Imax, is  ~ 0.45% for the brightest component. The remaining maser features show only residual Zeeman profiles with amplitudes at the rms level.

4. H2O masers and the SiO/CO outflows

Figure 6 shows that the three possible maser features unresolved in our data are distributed linearly. The velocity shift observed between VLSR  ≃ 5.7 km s-1 and higher velocities is oriented in a E–W direction (PA  ≃ 110°). The proximity to source Aa and the complex outflow configuration observed in this zone suggests that the maser emission is being pumped into the dense circumstellar material around this object.

thumbnail Fig. 2

Scheme of the distribution of dust and molecular material in I16293. The plus signal indicates the position of the peak intensity of our VLA water maser data. The ellipse is the deconvolved size of the dust continuum source A as derived by Rao et al. (2009). Triangles denote the position of the submillimeter condensations observed by Chandler et al. (2005). Stars denote the VLBI water maser detections of Imai et al. (2007) and straight lines denotes the direction of the CO and SiO outflows associated with this core (Rao et al. 2009).

Rao et al. (2009) reported a strong and very young SiO outflow extending toward the northwest-southeast (NW–SE) direction, with the SE lobe detected at redshifted LSR velocities (see Fig. 2). This SE component is also traced by a CO (3-2) outflow at VLSR  ≃ 9 km s-1  (Fig. 7 in Rao et al. 2009). The CO outflow has a dominant E–W orientation and has previously been reported in the literature (Yeh et al. 2008). Rao et al. (2009) claim that the SiO and CO outflows are centered on sources Ab and Aa, respectively. Nevertheless, it is hard to identify the powering sources in their molecular maps, since the separation between Aa and Ab (~0.6′′) is much smaller than the HPBW of their SiO and CO maps (~3′′).

thumbnail Fig. 3

Stokes I (black line) and linear polarization intensity (multiplied by a factor of ten, red dashed line) spectra of the water maser emission (upper panel). The spectra of polarization fraction (black line, left scale) and position angle (dashed line, right scale) are also shown (lower panel). (A color version of this figure is provided in the online version.)

Our maser features are also detected at redshifted velocities and lie in a similar direction to the CO outflow. On the other hand, their positions south of Aa suggest that they may be excited by the SiO outflow. The SE lobes observed for both SiO and CO at redshifted velocities (Figs. 6 and 7 from Rao et al. 2009) provide a close match to the LSR of our highest velocity features (VLSR  ~  9 km s-1). In addition, the VLA 3.7 cm data from Chandler et al. (2005) show two sources, A1 and A2, close to Aa (see Fig.  6). These authors refer to A2 as a protostar that powers a radio jet because of the bipolar shape seen in their high-resolution maps. On the other hand, source A1 would be an ionized region produced by shocks between this flow and nearby dense gas. Therefore, it is possible that the water masers are generated in the interaction regions traced by the radio sources and the outflows.

Water masers detections have also been reported toward source A through very long baseline (VLBI) observations (Imai et al. 2007). The milliarcsecond resolution data revealed a spot exactly over Aa and another one to the SW of it (see Fig. 2). While the former might have been excited in the circumstellar gas of source Aa, the latter is likely to have been created by shocks in the E–W outflow. The redshifted lobe of this emission has an open shell structure, hence we expect to find masers that are offset away from the outflow main axis.

5. The line-of-sight magnetic field strength in I16293

From the Zeeman splitting formalism, the magnetic field strength can be correlated to the fraction of circular polarization by (see, for instance, Fiebig & Guesten 1989) where the AF − F coefficient depends on the maser rotational levels F and F′, the intrinsic thermal linewidth Δνth, and the maser saturation degree, while ΔVI is the full width at half maximum of the total power spectrum. For the AF − F coefficient, we adopted the range of 0.012 to 0.018, since these values are consistent with models and observations (Nedoluha & Watson 1992; Vlemmings et al. 2002). The water maser spectrum is a blend of different velocity components, hence we adopt a linewidth ranging from 0.75 km s-1, which is typical of those found in other regions (Vlemmings et al. 2006), to 1.0 km s-1, the maximum estimated value in I16293 from our spectral profile. For these values, we find that BLOS ranges from  ~  − 94 to  ~  − 188 mG. The negative signal is inferred from the Stokes V shape, meaning that the field is pointing toward the observer. For narrower linewidths (~0.75 km s-1), which would be more realistic for resolved water maser lines, the field strength ranges from  ~  − 94 to  ~  − 141 mG (see Table 2). We note, however, that since the linear polarization fraction is about 3%, this maser is likely unsaturated and the field strength determination is a fair approximation.

thumbnail Fig. 4

Distribution of H2O linear polarization vectors in the brightest emission channel (VLSR  ≃ 7.4 km s-1). Contours are  − 50,  − 30, 30, 50, 500, 3  ×  103, 1  ×  104, 2  ×  104  ×  8 mJy beam-1. Only polarization vectors whose P > 1% are plotted. Vectors are sampled as 1/2 of a beam.

Figure 4 shows the linear polarization vectors at the velocity channel with the brightest emission. There is a degeneracy between the position angle of the linear polarization vectors and the magnetic field direction projected onto the plane-of-sky (POS). The position angle of the polarization is either parallel or perpendicular to the magnetic field in the POS for θ > θcrit or θ < θcrit = 55°, respectively (Goldreich et al. 1973), where θ is the angle between B and the maser propagation direction and θcrit is the so-called Van Vleck angle. Given that the linear polarization fraction Pl depends on θ and the maser saturation level, we should be able in principle to constrain the value of θ from our data and solve for this ambiguity. We used a non-local thermodynamic equilibrium radiative transfer model in order to fit the observed intensity I and polarization fraction Pl. However, since our line is a blend of features that were spatially unresolved, we were unable to determine θ. Although observations at more extended configurations would be needed to achieve this goal, we can at least claim that the POS field topology is quite ordered in both cases, i.e., it would remain ordered if rotated by 90°.

thumbnail Fig. 5

Stokes I (upper panel) and Stokes V (lower panel) spectra of the water maser emission. The red dashed line is the scaled derivative of the total power I in order to guide the eye. (A color version of this figure is provided in the online version.)

According to the submm data of Rao et al. (2009), the direction of the dust polarization vectors associated with source A averages to 95°, which implies that the POS magnetic field lines have a PA of 5° with respect to the north direction. Assuming that the 22 GHz polarization vectors in Fig. 4 are parallel to the POS magnetic field, the main direction of the field lines is  − 23° (or its supplement, 157°) with respect to the north direction, as reported in Sect. 3.2. On the other hand, if our polarization vectors were perpendicular to the POS field, their main direction would then be 67°. Although two distinct scales are being compared, it is interesting to point out that in both cases there is clearly a change in the direction of the field lines on the small scales traced by our VLA data. In the case of parallel field lines and polarization vectors, such a change could be produced by the interaction of the SiO outflow and the magnetic field, given that the field lines lie almost in the same direction as the outflow (PA  ~ 145°, Fig. 2).

6. Preshock magnetic fields and post-shock densities

The only previous estimation of the magnetic field strength toward I16293 had been performed by Claussen et al. (2003), who used VLBI observations of water masers to estimate BLOS ~  − 40 mG. Their maser spots appear to be associated with the A1 cm source (Fig. 6), which lies at a projected separation of  ~ 20 AU at VLSR  ~ 2.5 km s-1 (Claussen, priv. comm.). Since this LSR velocity is similar to the core ambient velocity, this emission is unlikely to trace the same region as our VLA data. In any case, our study is one of the first determinations of the field strength in a low-mass young stellar object at densities higher than 108 cm-3. As discussed in Sect. 4, the observed maser emission is excited in zones of compressed gas produced by shocks between the outflows and the ambient gas. In an ionized medium, shocks compress the gas by a factor mA, the Alfvénic Mach number that is defined as vs/vA, where vs is the shock velocity and vA the Alfvén speed. The magnetic field B0 in the pre-shock zone is also enhanced by a factor mA in the shocked gas so that Bs ~ mA × B0. Since vA varies with B0(n) − 1/2, then Bs is proportional to vs(n)1/2, where n is the volume density of the pre-shock region. Rearranging this expression, the pre-shock density can be written as (Kaufman & Neufeld 1996) (3)The LOS magnetic field strength calculated in Sect. 5 can be rescaled in order to obtain the total field strength Btot = Bs in the shocked zone. According to Crutcher et al. (1999), the squared LOS magnetic field  ⟨ BLOS ⟩ 2 averaged over several lines of sight in the Galaxy is  ~ . The shock velocities are retrieved from submm water emission detected in shocked zones around I16293 (Ristorcelli et al. 2005). These data were modeled as C-shocks (low velocity, non-dissociative shocks, Kaufman & Neufeld 1996) with shock velocities of 12 − 15 km s-1. Applying these values to Eq. (3), we obtain pre-shock densities of  ~  cm-3, which is consistent with the typical densities at which water masers are eventually pumped. The upper/lower limits to the number density refer to the Bs determinations for the two limiting AF − F cases (Table 2).

thumbnail Fig. 6

Possible independent maser features (stars) derived using Gaussian fits. The accuracy in the position of each maser feature is  ~2 mas, while the positions of the submm and cm sources are accurate to within  ~25 mas and  ~10 mas, respectively. The numerical labels indicate the maser velocities with respect to the local standard of rest (LSR). The 3.7 cm sources A1 and A2 (triangles) and the submm sources Aa and Ab (circles) are also shown (Chandler et al. 2005). The sizes of the symbols are not proportional to either their deconvolved sizes or their intensities.

Table 2

B-field strength calculation from Eq. (2) (assuming ΔVI = 0.75 km s-1).

If the magnetic field dominates the outflow dynamics, the magnetic pressure should be similar to the shock ram pressure, i.e., (4)where ρ0 is the pre-shock density. For the field strength Bs determined in the previous section and shock velocities of Ristorcelli et al. (2005), we find a mean pre-shock density of  ~ 2 × 108 cm-3, which is consistent with the pre-shock densities calculated from Eq. (3) and in line with a magnetically controlled outflow evolution.

Assuming magnetic flux-freezing, the compression of an ordered magnetic field (B0) amplifies the field strength in proportion to the gas density according to (5)where ns is the post-shock density. If we consider that the shock compresses the well-established field geometry in the submm maps of Rao et al. (2009), which predict a (scaled) B0 ≃ 9 mG at densities n0 ≃ 4.9 × 107 cm-3, and using the post-shock field intensity Bs calculated from our data, we find post-shock densities of  ~1.3 × 109 cm-3, which is expected for effective water maser pumping. In a more conservative approach, we rescale the magnetic field strength of Rao et al. (2009) using B ∝ nκ, where κ may have distinct values depending on the dynamics of the core: κ = 0.5 for a cloud evolution controlled by the magnetic field (Crutcher et al. 1999), κ = 0.47, in the case of magnetic evolution as numerically modeled by Fiedler & Mouschovias (1993), and κ = 0.67 for a free-fall cloud collapse. For the three cases, the pre-shock magnetic field in the calculated pre-shock densities ranges between  ~24 and 36 mG. Applying again Eq. (5), we obtain post-shock densities in the range  ~2−5 × 109 cm-3, which are again consistent with typical H2O maser densities, but we are unable to discern in which dynamical regime, if magnetic or turbulent, the core evolves.

The H2O maser emission of I16293 is strong (>200 Jy) and stable over a few weeks but highly variable over months, as shown by surveys performed toward this source (Wilking & Claussen 1987; Claussen et al. 1996; Furuya et al. 2003). We therefore hope to carry out VLBI observations in the future. With these observations, the blended maser components found in this work will likely be resolved both spatially and spectroscopically. We will then be able to apply radiative transfer models similar to those in Vlemmings (2006) to derive the full 3-D magnetic field properties as has already been done in other star-forming regions (e.g., Cepheus A HW2, Vlemmings et al. 2006).

7. Conclusions

We have reported on the H2O (616 − 523) maser emission observed with the VLA/EVLA toward the low-mass source IRAS 16293-2422. This is one of the first estimations of the line-of-sight magnetic field strength in a low-mass protostar at densities higher than 108 cm-3. From this study, we have obtained insights into the dynamics of this object at such high densities. Our main conclusions are:

  • We detect strong water-maser emission associated with the submm source Aa.

  • The maser spectrum has a non-Gaussian profile, which is indicative of at least three maser features covering a velocity range of  ~3.5 km s-1 distributed in a linear configuration. These components are likely excited in zones of compressed gas produced by shocks in the outflow activity of I16293.

  • The obtained Stokes V spectrum is consistent with Zeeman emission. The mean LOS magnetic field strength is  ~ 113 mG. The POS field topology retrieved from the linear polarimetry shows an ordered pattern consistent with larger-scale field morphologies.

  • The post- and pre-shock densities calculated from our field strength estimation are consistent with the typical densities expected for H2O maser pumping.

  • The dynamics of the outflow evolution in these sources are likely regulated by the magnetic field, since the magnetic pressure in the shocked gas is similar to the pre-shock ram pressure.


1

The National Radio Astronomy Observatory – NRAO – is a facility of the National Science Foundation, operated under cooperative agreement by Associated Universities, Inc.

Acknowledgments

The authors would like to thank the anonymous referee for a constructive report that helped to improve the paper. We are also grateful to Gabriele Surcis for running the radiative transfer models on the maser data. This research was partially supported by the Deutsche Forschungsgemeinschaft (DFG) through the Emmy Noether Research grant VL 61/3-1 and through SFB 956. F.O.A., J.M.G., and J.M.T. also acknowledge support from MICINN (Spain) AYA2011-30228-C03 and AGAUR (Catalonia) 2009SGR1172 grants.

References

  1. Bisschop, S. E., Jørgensen, J. K., Bourke, T. L., Bottinelli, S., & van Dishoeck, E. F. 2008, A&A, 488, 959 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Briggs, D. 1995, Ph.D. Thesis, The New Mexico Institute of Mining and Technology [Google Scholar]
  3. Ceccarelli, C., Loinard, L., Castets, A., Tielens, A. G. G. M., & Caux, E. 2000, A&A, 357, L9 [NASA ADS] [Google Scholar]
  4. Chandler, C. J., Brogan, C. L., Shirley, Y. L., & Loinard, L. 2005, ApJ, 632, 371 [NASA ADS] [CrossRef] [Google Scholar]
  5. Claussen, M. J., Wilking, B. A., Benson, P. J., et al. 1996, ApJS, 106, 111 [NASA ADS] [CrossRef] [Google Scholar]
  6. Claussen, M., Sarma, A. P., Wootten, A., Marvel, K. B., & Wilking, B. 2003, in Meeting Abstracts, BAAS, 35, 1362 [NASA ADS] [Google Scholar]
  7. Crutcher, R. M., Roberts, D. A., Troland, T. H., & Goss, W. M. 1999, ApJ, 515, 275 [NASA ADS] [CrossRef] [Google Scholar]
  8. Elitzur, M., Hollenbach, D. J., & McKee, C. F. 1989, ApJ, 346, 983 [NASA ADS] [CrossRef] [Google Scholar]
  9. Fiebig, D., & Guesten, R. 1989, A&A, 214, 333 [NASA ADS] [Google Scholar]
  10. Fiedler, R. A., & Mouschovias, T. C. 1993, ApJ, 415, 680 [NASA ADS] [CrossRef] [Google Scholar]
  11. Furuya, R. S., Kitamura, Y., Wootten, A., Claussen, M. J., & Kawabe, R. 2003, ApJS, 144, 71 [NASA ADS] [CrossRef] [Google Scholar]
  12. Goldreich, P., Keeley, D. A., & Kwan, J. Y. 1973, ApJ, 179, 111 [NASA ADS] [CrossRef] [Google Scholar]
  13. Goodman, A. A., Jones, T. J., Lada, E. A., & Myers, P. C. 1992, ApJ, 399, 108 [Google Scholar]
  14. Goodman, A. A., Jones, T. J., Lada, E. A., & Myers, P. C. 1995, ApJ, 448, 748 [NASA ADS] [CrossRef] [Google Scholar]
  15. Imai, H., Nakashima, K., Bushimata, T., et al. 2007, PASJ, 59, 1107 [NASA ADS] [Google Scholar]
  16. Jørgensen, J. K., Bourke, T. L., Nguyen Luong, Q., & Takakuwa, S. 2011, A&A, 534, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Kaufman, M. J., & Neufeld, D. A. 1996, ApJ, 456, 250 [NASA ADS] [CrossRef] [Google Scholar]
  18. Kuan, Y., Huang, H., Charnley, S. B., et al. 2004, ApJ, 616, L27 [NASA ADS] [CrossRef] [Google Scholar]
  19. Lazarian, A., Goodman, A. A., & Myers, P. C. 1997, ApJ, 490, 273 [NASA ADS] [CrossRef] [Google Scholar]
  20. Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJ, 675, L29 [NASA ADS] [CrossRef] [Google Scholar]
  21. Looney, L. W., Mundy, L. G., & Welch, W. J. 2000, ApJ, 529, 477 [NASA ADS] [CrossRef] [Google Scholar]
  22. Mundy, L. G., Wootten, H. A., & Wilking, B. A. 1990, ApJ, 352, 159 [NASA ADS] [CrossRef] [Google Scholar]
  23. Nedoluha, G. E., & Watson, W. D. 1992, ApJ, 384, 185 [NASA ADS] [CrossRef] [Google Scholar]
  24. Pérez-Sánchez, A. F., Vlemmings, W. H. T., & Chapman, J. M. 2011, MNRAS, 418, 1402 [NASA ADS] [CrossRef] [Google Scholar]
  25. Rao, R., Girart, J. M., Marrone, D. P., Lai, S., & Schnee, S. 2009, ApJ, 707, 921 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ristorcelli, I., Falgarone, E., Schöier, F., et al. 2005, in IAU Symp., 235, 227P [Google Scholar]
  27. Sanna, A., Moscadelli, L., Cesaroni, R., et al. 2010a, A&A, 517, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Sanna, A., Moscadelli, L., Cesaroni, R., et al. 2010b, A&A, 517, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Sarma, A. P., Troland, T. H., & Romney, J. D. 2001, ApJ, 554, L217 [NASA ADS] [CrossRef] [Google Scholar]
  30. Stark, R., Sandell, G., Beck, S. C., et al. 2004, ApJ, 608, 341 [NASA ADS] [CrossRef] [Google Scholar]
  31. Surcis, G., Vlemmings, W. H. T., Curiel, S., et al. 2011a, A&A, 527, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Surcis, G., Vlemmings, W. H. T., Torres, R. M., van Langevelde, H. J., & Hutawarakorn Kramer, B. 2011b, A&A, 533, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Tamura, M., Hayashi, S. S., Yamashita, T., Duncan, W. D., & Hough, J. H. 1993, ApJ, 404, L21 [NASA ADS] [CrossRef] [Google Scholar]
  34. Terebey, S., Vogel, S. N., & Myers, P. C. 1992, ApJ, 390, 181 [NASA ADS] [CrossRef] [Google Scholar]
  35. Vlemmings, W. H. T. 2006, A&A, 445, 1031 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Vlemmings, W. H. T. 2008, A&A, 484, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Vlemmings, W. H. T., Diamond, P. J., & van Langevelde, H. J. 2002, A&A, 394, 589 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Vlemmings, W. H. T., Diamond, P. J., van Langevelde, H. J., & Torrelles, J. M. 2006, A&A, 448, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Walker, C. K., Lada, C. J., Young, E. T., & Margulis, M. 1988, ApJ, 332, 335 [NASA ADS] [CrossRef] [Google Scholar]
  40. Wardle, J. F. C., & Kronberg, P. P. 1974, ApJ, 194, 249 [NASA ADS] [CrossRef] [Google Scholar]
  41. Wilking, B. A., & Claussen, M. J. 1987, ApJ, 320, L133 [NASA ADS] [CrossRef] [Google Scholar]
  42. Wootten, A. 1989, ApJ, 337, 858 [NASA ADS] [CrossRef] [Google Scholar]
  43. Yeh, S. C. C., Hirano, N., Bourke, T. L., et al. 2008, ApJ, 675, 454 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

H2O maser components in IRAS 16293-2422.

Table 2

B-field strength calculation from Eq. (2) (assuming ΔVI = 0.75 km s-1).

All Figures

thumbnail Fig. 1

Channel maps of the H2O emission toward IRAS 16293-2422. Contours are  − 50, 50, 500, 3 × 103, 1 × 104, 2 × 104 times 8 mJy beam-1, the rms noise of the maps with lower intensities (see Sect. 2). The peak flux of 167 Jy beam-1 occurs at VLSR = 7.4 km s-1. The LSR velocity of each channel is labeled.

In the text
thumbnail Fig. 2

Scheme of the distribution of dust and molecular material in I16293. The plus signal indicates the position of the peak intensity of our VLA water maser data. The ellipse is the deconvolved size of the dust continuum source A as derived by Rao et al. (2009). Triangles denote the position of the submillimeter condensations observed by Chandler et al. (2005). Stars denote the VLBI water maser detections of Imai et al. (2007) and straight lines denotes the direction of the CO and SiO outflows associated with this core (Rao et al. 2009).

In the text
thumbnail Fig. 3

Stokes I (black line) and linear polarization intensity (multiplied by a factor of ten, red dashed line) spectra of the water maser emission (upper panel). The spectra of polarization fraction (black line, left scale) and position angle (dashed line, right scale) are also shown (lower panel). (A color version of this figure is provided in the online version.)

In the text
thumbnail Fig. 4

Distribution of H2O linear polarization vectors in the brightest emission channel (VLSR  ≃ 7.4 km s-1). Contours are  − 50,  − 30, 30, 50, 500, 3  ×  103, 1  ×  104, 2  ×  104  ×  8 mJy beam-1. Only polarization vectors whose P > 1% are plotted. Vectors are sampled as 1/2 of a beam.

In the text
thumbnail Fig. 5

Stokes I (upper panel) and Stokes V (lower panel) spectra of the water maser emission. The red dashed line is the scaled derivative of the total power I in order to guide the eye. (A color version of this figure is provided in the online version.)

In the text
thumbnail Fig. 6

Possible independent maser features (stars) derived using Gaussian fits. The accuracy in the position of each maser feature is  ~2 mas, while the positions of the submm and cm sources are accurate to within  ~25 mas and  ~10 mas, respectively. The numerical labels indicate the maser velocities with respect to the local standard of rest (LSR). The 3.7 cm sources A1 and A2 (triangles) and the submm sources Aa and Ab (circles) are also shown (Chandler et al. 2005). The sizes of the symbols are not proportional to either their deconvolved sizes or their intensities.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.