Free Access
Issue
A&A
Volume 618, October 2018
Article Number A163
Number of page(s) 8
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201833170
Published online 23 October 2018

© ESO 2018

1 Introduction

The significance of magnetic fields on the star formation process is subject of ongoing investigations (e.g. Kandori et al. 2007; Galli 2009; Franco et al. 2010; Federrath 2015; Lizano & Galli 2015; Sadavoy et al. 2018). Recent studies reveal comparatively long lifetimes of massive (e.g. Orion) clouds (~10 Myr; e.g. Murray 2011; Meidt et al. 2015) and suggest “dynamically relevant” magnetic fields (e.g. Alves et al. 2008; Sugitani et al. 2010; Commerçon et al. 2011; Palau et al. 2013; Stutz & Gould 2016; Busquet et al. 2016). The structure of intermediate-mass clouds (e.g. Taurus), showing lifetimes of an order of magnitude shorter (~3 Myr; e.g. Hartmann et al. 2001, 2016; Ballesteros-Paredes & Hartmann 2007), however, can be at least proximately explained by turbulence (e.g. Ostriker et al. 2001; Offner et al. 2008; Kritsuk et al. 2013). In turn, recent studies on Bok globules, the lowest mass clouds known at this time, suggest that globules appear to be magnetically supported (Stutz et al. 2007; Alves et al. 2011, 2014; Bertrang et al. 2014; Das et al. 2016). These studies have been conducted on globules with rather simple density structures. In this study, we aim to test the magnetic field influence on globules with more complex density structures.

The structure of magnetic fields in the dusty envelopes around young stellar objects can be derived from the polarization of background starlight due to dichroic extinction and thermal emission by dust grains (e.g. Weintraub et al. 2000; Girart et al. 2008). In the case of polarization due to thermal emission, the polarization direction is perpendicular to the magnetic field lines, as the dust grains become partially aligned to the magnetic field, with their major axes perpendicular to the magnetic field lines, projected onto the plane-of-sky (POS; e.g. Hoang et al. 2007; Lazarian 2007; Bertrang & Wolf 2017). The light from the background stars that runs through the star-forming region also becomes polarized due to dichroic absorption by the dust grains. In this case, the polarization direction directly traces the magnetic field lines, projected onto the POS (Weintraub et al. 2000).

The influence of magnetic fields on low-mass star formation is best verified when studied independently of other phenomena such as turbulence or stellar feedback. Bok globules (Bok & Reilly 1947) arise as the ideal candidates to study the magnetic field influence on low-mass star forming regions. Bok globules are small, isolated, relative simple structured molecular clouds (e.g.; Clemens et al. 1991), and are often associated with low-mass star formation (Clemens & Barvainis 1988; Yun & Clemens 1994; Launhardt & Henning 1997). Bok globules in particular are ideal environments to study the correlation between protostellar collapse, fragmentation, and magnetic fields, as these objects are less affected by large scale turbulence and other nearby star-forming events.

The classical method to characterize magnetic fields is the linearly polarized dust emission and dichroic absorption. Thermal emission of the dust grains, observable in the sub-millimetre (sub-mm) wavelength range, traces the polarization of the densest, innermost part of the globule, due to the sensitivity of millimetre/sub-mm telescopes. The outer, less dense parts of the globule are traced by the polarization of background starlight due to dichroic absorption, observable in the near-infrared (near-IR)/optical wavelength range. Thus, by applying multi-wavelength polarimetry, it is possible to trace the magnetic field on large scales across the globule.

In the following, we describe our target selection (Sect. 2), observations (Sect. 3), and data reduction (Sect. 4). We analyse the polarization maps in Sect. 5. In Sect. 6, we discuss the magnetic field, its relation to outflows, as well as the magnetic field influence on star formation as a function of mass before we conclude in Sect. 7.

thumbnail Fig. 1

Panel a: near-IR polarization segments (red) and optical polarization segments (blue) are plotted over a 15′ × 15′ R-band DSS image of the field containing CB34. A 10% polarization segment is drawn as reference for both data sets into the lower left corner. Contours correspond to JCMT/SCUBA 850 μm dust continuum emissions from 0.1 to 0.7 Jy beam−1, with intervals of 0.1 Jy beam−1 (Das et al. 2016). Panel b: zoom into the dense centre of CB34. The near-IR polarization segments (red) and the archival sub-mm polarization segments (blue) are plotted over the total continuum intensity map at 850 (Das et al. 2016). The sub-mm segments are binned over a 10″ × 10″ grid. The solid (dashed) black (grey) contour lines correspond to blueshifted (redshifted) integrated CO emission. Contour levels are spaced at 0.1 K km s−1 intervals of 0.5 K km s−1 (black) and 0.2 K km s−1 intervals of 0.8 K km s−1 (white), respectively (from Yun & Clemens 1994; beam size of 48′′). A 5% polarization segment is drawn as reference for both data sets into the lower right corner. The well-ordered polarization pattern indicates dominant magnetic fields on scales of 103 –106 au (~ 0.0055pc).

2 Description of the sources

Our aim is to verify the large scale magnetic field influence on low-mass star formation. As stated in Sect. 1, Bok globules are ideal candidates for this kind of study. Furthermore, we must restrict ourselves to globules with available archival sub-mm or optical polarization maps that trace the magnetic field structure at scales not observable in the near-IR wavelength range. We further aim to expand our previous study of the magnetic field structure in low-mass star-forming regions (Bertrang et al. 2014)to more complex structured objects. We find that these criteria are satisfied by three Bok globules: CB34, CB56, and [OMK2002]18.

The Bokglobule CB34 is located at a distance of ~1.5 kpc (Launhardt & Henning 1997). The inner region of the globule discloses three dense cores (C1, C2, C3), each extends by ~ 2.5 pc in diameter (Codella & Scappini 2003; Launhardt et al. 2010). Their masses are 52 M (C1), 94 M (C2), and 13 M (C3), respectively (Fig. 1; Codella & Scappini 2003). These cores exhibit ongoing star formation, which is based on the presence of a water maser and several collimated outflows (Fig. 1; Yun & Clemens 1994; Huard et al. 2000; Gómez et al. 2006). The chemical age of the globule is estimated to be >105 yr (Launhardt et al. 2010), with the presence of a pre-main-sequence star with an age of ~ 106 yr located in the core C1 (Alves et al. 1997). This massive globule encloses within its centre (4.05 × 2.25) × 105 au (~ 2 × 1pc) about 110–170 M (Launhardt et al. 1998, 2010; Codella & Scappini 2003).

The Bok globule CB56 is a compact, irregular shaped cloud with a major and minor axis of 4.′′ 5 and 2.′′ 5, respectively (Dutra & Bica 2002). Its major axis has a position angle of 170° with respectto the galactic plane. CB56 is associated with two infrared point sources (IRAS 07125-2503 and IRAS 07125-2507; Clemens & Barvainis 1988). The distance to CB56 has yet to be determined. The optical polarization signal of this globule is significantly stronger than the interstellar polarization in the optical wavelength range within 200 pc (Chakraborty et al. 2014; Gontcharov 2017; Gontcharov, priv. comm.). Together with the strong near-IR polarization, this indicates that CB56 is located within a distance of 200 pc and the measured polarization signal is associated with CB56.

The source [OMK2002]18 is a low-mass star-forming region at a distance of 140 pc (Maheswar et al. 2011). It consists of the two dense clouds [OMK2002]18a and [OMK2002]18b. Both clouds have a radius of ~ 4.125 × 103 au and enclose masses of 2.3 and 1.1 M, respectively (Onishi et al. 2002). OMK[2002]18 is associated with one compact infrared source (IRAS 04191+1522; Onishi et al. 2002).

3 Observations

The observations were performed in the near-IR with the instrument Son OF ISAAC (SOFI) at the 3.58 m New Technology Telescope (NTT) from November 11 to 15, 2015. SOFI/NTT is mounted at the Nasmyth A focus, equipped with a 1024 × 1024 Hawaii Rockwell array optimized for wavelengths of 1–2.5 μm. In SOFI/NTT, a single Wollaston prism is used for polarization observations. In this observing mode, the polarized flux is measured simultaneously at two different angles that differ by 90°. To derive the linear polarization degree and orientation of an object, two observations must be performed at each pointing of the telescope with different orientations of the Wollaston prism, typically 0° and 45°. This is realized by a rotation of the complete instrument. To avoid overlapping between different polarization images, an aperture mask of three alternating opaque and transmitting strips of about 40″ × 300″ for SOFI/NTT is used. We carried out Js-band polarization observations of three fields of CB34 and [OMK2002]18ab as well as one field of CB56.

4 Data reduction

For the data reduction, we applied a specialized pipeline created to work with polarization data obtained with SOFI/NTT (Bertrang et al. 2014). This pipeline performs bias correction, flat-fielding, and instrumental polarization extraction as well as the calculation of the Stokes parameters I, Q, and U via aperture photometry. Table 1summarizes the information of polarized and unpolarized standard stars used to determine the instrumental polarization.

A correction factor for the Wollaston prism was computed for the reduction of the instrumental polarization (see Appendix). This correction factor was tested for the polarized standard star, and as a successful test the polarization data obtained after bias correction correspond well to literature values.

Table 1

Polarization standard stars.

5 Polarization maps

In this section, we present the polarization maps of the three Bok globules CB34, CB56, and OMK[2002]18. We compare our near-IR data to archival polarization data in the optical and sub-mm wavelength range. As stated in Sect. 1, the polarization at these different wavelength ranges arises from different sources. While the polarization in the sub-mm range arises from thermal emission, in the near-IR and the optical it is caused by dichroic absorption. Thus, the polarization of a dust grain observed in the sub-mm is orientated perpendicular to the polarization of the very same dust grain observed in the near-IR or optical.

For CB34, in addition to our new data in the near-IR wavelength range, archival data are available in both the optical and sub-mm wavelength ranges (Das et al. 2016). This enables an analysis over the whole range of scales in this globule, from its most dense centre (sub-mm) towards the less dense outer envelope (near-IR, optical).

For CB56, the data set is restricted to archival data in the optical wavelength range (Paul et al. 2012; Chakraborty et al. 2014)in addition to the newly obtained near-IR data. However, those two wavelength ranges cover this object almost completely. For this comparison, we have made use of the optical polarization data published in Chakraborty et al. (2014), given that it is the most comprehensive data set available at this wavelength range.

For [OMK2002]18, archival data were available only in the sub-mm wavelength range (Matthews et al. 2009), in addition to the newly obtained near-IR polarization data. With these, an analysis of the polarization of the most dense centre of the globule is possible (sub-mm), as well as an analysis of the less dense outer regions (near-IR).

In the following, we discuss the polarization maps separately, first focusing on the new near-IR data, and subsequently, comparing it to the maps of the optical and sub-mm wavelength ranges.

5.1 CB34

The polarization maps of the Bok globule CB34 are shown in Fig. 1. CB34 disclose a strong polarization signal, P ≥ 1%, in the outer, less dense regions observed in the near-IR (Figs. 4 and 5). The polarization segments in the near-IR wavelength range cover this Bok globule on scales of 104 105 au (~ 0.050.5 pc). The polarization segments in the near-IR wavelength range are well ordered and orientated towards north-west with a mean orientation of 111° ± 6°.

Most of the outer regions of this globule are traced in both near-IR and optical polarization. Both data sets fit well in orientation (within the errors). The mean value of the polarization angle in the near-IR wavelength range in this regions is relatively consistent with the position angle of the galactic plane at the latitude of CB34, θGP ≈ 148.7°, which suggests that the magnetic field is coupled to the Galactic magnetic field in the observed POS. This is similar to the findings in the optical (Das et al. 2016).

A puzzling finding is the relation between the polarization degrees in the near-IR and the optical. As opposed to the empirical Serkowski relation (Serkowski et al. 1975), the near-IR signal is stronger polarized than the optical. Similar discrepancies have also been found, although in very few other globules, and are tentatively explained by either a limited location of origin of polarization or the presence of two different dust populations (Clemens 2012; Alves et al. 2014). Recently, it has been suggested that not only the grain size but also the chemical composition of grains alter the extinction dependency of polarization (Papoular 2018). Future theoretical work is needed to investigate this in more detail.

In the sub-mm, the polarization segments can be associated with the cores C1 and C2, while there is no detection on top of C3 (Fig. 1). The region observed in the sub-mm overlaps with our near-IR data. The polarization segments in the sub-mm wavelength range were plotted using the original data from Das et al. (2016). Thus, it has to be kept in mind that the segments must be rotated by 90° to be aligned with the magnetic field orientation. In this region, the orientation of the polarization segments obtained in the near-IR wavelength range fit very well to the sub-mm polarization segments associated with core C1, but deviate stronger from those of core C2. The polarization segments in the near-IR wavelength range are located just in between the cores C1 and C2. In this region, the mean polarization angle in the near-IR wavelength is γ¯near-IR111.9°±11.5°$\bar{\gamma}_{\rm{near-IR}}\approx 111.9^{\circ} \pm 11.5^{\circ}$, which suggests that in this case the orientation of the magnetic field in the centre of the globule is different from that of the Galactic magnetic field. Similar findings were obtained in Das et al. (2016) using the data from the sub-mm wavelength range.

The change in the polarization pattern from the outer to the central part of CB34 indicates that the magnetic field is coupled on large scales to thegalactic magnetic fields, while it disconnects in the centre where the globule is fragmented into several cores. The scale of the transition is ~104 au (~ 0.05 pc; Fig. 1).

thumbnail Fig. 2

Panel a: near-IR polarization segments plotted over a 15′× 15′ R-Band DSSimage of the field containing CB56. A 10% polarization segment is drawn for reference in the lower left corner. Only the segments with Pσp ≥ 3 are plotted. Panel b: optical polarization segments of CB56, plotted over a 10′ × 10′ R-Band DSS image (Chakraborty et al. 2014). A 1% polarization segment is drawn for reference in the lower left corner. The near-IR and optical polarization is strong and well ordered, indicating dominant magnetic fields on scales of 104 105 au (~ 0.050.5pc).

5.2 CB56

The polarization maps of the Bok Globule CB56 are shown in Fig. 2. CB56 discloses a strong polarization signal with polarization degrees of mostly 5–12% (Fig. 4). The polarization structure of this globule in the near-IR is well ordered and predominantly orientated towards north-north-west (Fig. 5).

The region of CB56 which is covered by our near-IR observations overlays with the archival optical data set (Chakraborty et al. 2014). The globule discloses a weak polarization signal in the optical wavelength range, with an average of P ≈ 1% (see Sect. 5.1 for a discussion about the Serkowsky relation).

The polarization segments obtained in both the optical and near-IR wavelength ranges are very well aligned within themselves, and fit well to each other. We find a slight deviation of the polarization orientation towards the dense centre of the globules when switching from the optical to the near-IR wavelength range.

thumbnail Fig. 3

Panel a: near-IR polarization segments plotted over a ~20′× 20′ R-band DSS image of the field containing [OMK2002]18. A 10% polarization segment is drawn for reference in the lower left corner. Only the segments with Pσp ≥ 3 are plotted. Panel b: sub-mm polarization segments of [OMK2002]18, plotted over a 20″ × 20″ grid. Contours correspond to JCMT/SCUBA 850 μm dust continuum emissions from 0.05 to 0.25 Jy beam−1, with intervals of 0.05 Jy beam−1 (Matthews et al. 2009; ©AAS, reproduced with permission). Only segments with Pσp ≥ 2 are plotted. The near-IR polarization is strong and well ordered, indicating dominant magnetic fields on scales of 104 –105 au (~ 0.05–0.5pc) in the western region of the globule.

thumbnail Fig. 4

Distribution of polarization degree, P for P ≥ 3%, counts given by N(P), of CB34, CB56, and OMK[2002]18, observed in the near-IR.

thumbnail Fig. 5

Distribution of polarization angle, γ, counts given by N(γ) of CB34, CB56, and OMK[2002]18, observed in the near-IR. The dashed lines represent the mean polarization angles, while dotted lines represent the corresponding 1σ deviation.

5.3 OMK[2002]18

The polarization maps of OMK[2002]18 is shown in Fig. 3. OMK[2002]18 shows a strong polarization signal, 2 ≥ P ≥ 10%, in all the regions observed in the near-IR (Fig. 4).

Two clearly divided regions of the globule can be observed in the near-IR wavelength range. The western region of the globule shows well-ordered polarization segments in north-north-west orientation, with strong polarization degrees of P ≈ 8%. The easternregion of the globule shows an unordered pattern with a much weaker polarization signal, P ≈ 3% indicative for depolarization effects (e.g. Brauer et al. 2016). Considering the separation between these regions, uniformly ordered polarization segments are not necessarily expected.

A spatial gap is seen between the sub-mm and near-IR polarization segments, which is about 3.3 ×104 au (~ 0.16pc). It originates inthe sensitivity of SOFI/NTT, resp. SCUBA/JCMT (Bertrang et al. 2014). The near-IR polarization segments in the western part of this object are well aligned with the most external polarization segments in the sub-mm wavelength range.

6 Magnetic fields

In our analysis we assumed that the magnetic field is orientated perpendicular to the measured polarization segments in the sub-mm and parallel orientated to the measured polarization segments in the near-IR and in the optical. This widely applied concept is based on the finding that, independently of the alignment mechanism, charged interstellar dust grains would have a substantial magnetic moment leading to a rapid precession of the grain angular momentum around the magnetic field direction, implying a net alignment of the grains with the magnetic field (e.g. Draine & Weingartner 1997; Lazarian 2007). However, one has to keep in mind that polarization observations strongly suffer from projectional effects along the line-of-sight (LOS). Thus, for a comprehensive understanding of the magnetic field structure additional 3D radiative transfer modelling is essential but beyond the scope of this study.

In general, a high polarization signal is connected to a magnetic field strong enough to align the dust grains along the LOS. In Sect. 5, we find well-ordered polarization segments on scales of 103 au to 105 au (~ 0.005–0.5pc) in all of the three Bok globules, CB34, CB56, and OMK[2002]18. This implies the presence of dominant large scale magnetic fields in each globule.

In addition to this analysis, we were able to further estimate the magnetic field strength in the POS by using a method first proposed by Chandrasekhar & Fermi (1953; CF method) |BPOS|=α4π3ρgasυturbσγ,\begin{equation*}|B_{\rm{POS}}| = \alpha\cdot\sqrt{\frac{4\pi}{3}\rho_{\rm{gas}}}\frac{\upsilon_{\rm{turb}}}{\sigma_{\gamma}}, \end{equation*}(1)

where ρgas is the gas density in units of g cm−3, υturb is the rms turbulence velocity in units of cm s−1, and σγ the standard deviation of the polarization angles in radians. In this method, it is assumed that the magnetic field is frozen in the cloud material.

Originally, α = 1 for the CF method. Different studies (e.g, Zweibel 1990; Myers & Goodman 1991) have suggested a lower value for α. Ostriker et al. (2001) found that, for σγ ≲ 25°(0.44 radians), using α ~0.5 yields a good estimate for the magnetic field strength.

The total gas density, ρgas, is given by ρgas=1.36nH2MH2,\begin{equation*}\rho_{\rm{gas}} = 1.36 \, n_{\textrm{H}_2} \, M_{\textrm{H}_2}, \end{equation*}(2)

where MH2=2.0158$M_{\textrm{H}_2} = 2.0158$ amu is the mass of a H2 molecule.

The rms turbulence velocity, υturb, is given by (Wang et al. 1995) υturb=υ2.35,\begin{equation*} \upsilon_{\rm{turb}} = \frac{\bigtriangleup\upsilon}{2.35}, \end{equation*}(3)

where △ υ is the FWHM line-width, measured at quiescent positions located away from the emission peaks.

For the outer, less dense regions of CB34, the mean density, nH2, is 3.76 × 102 cm−3 (Das et al. 2016) on average, while the FWHM line-width, △ υ, corresponds to 1.37 km s−1 (Sepúlveda et al. 2011). Considering the globule’s morphology, we split CB34 into two regions: south-east (SE) and north-west (NW) of the centre. Our near-IR observations trace the central region itself too sparsely for a robust application of the CF method. All parameters used for each region are listed in Table 2.

We determined the dispersion of the polarization angles in the SE region to σγ = 31° and in the NW region to σγ = 29°. Since the correction found by Ostriker et al. (2001) is not applicable, the derived magnetic field strengths are upper limits. Applying Eq. (1) (α = 1), we derived magnetic field strengths of BNIRNW9.6μG$B^{\rm{NW}}_{\rm{NIR}} \approx 9.6\,\mu\rm{G}$ (north-west)and BNIRSE9.1μG$B^{\rm{SE}}_{\rm{NIR}} \approx 9.1\,\mu\rm{G}$ (south-east), respectively. The similar magnetic field strengths in both regions, along with the well-ordered polarization segments (Sect. 5), implies a uniform large scale magnetic field across the globule.

Towards the globule’s centre, the magnetic field strength increases to 34 μG for core C1 and 70 μG for core C2 (Das et al. 2016). This is consistent with the theoretical picture of magnetically supported star formation (e.g. Hennebelle & Ciardi 2009; Peters et al. 2011; Federrath 2016).

For CB56 and OMK[2002]18, necessary information on the gas densities and turbulence velocities in the regions traced by the presented polarization data is not available. Hence, it is not possible to estimate the magnetic field strength of these globules.

Table 2

Gas densities, gas velocities, polarization, and magnetic field strengths of the regions of CB34 traced in the near-IR.

6.1 Correlation between magnetic field structure and the CO outflow

Alongside the collapse of the dust in Bok globules, magnetic fields are also presumed to influence the formation of circumstellar disks and outflows (e.g. Matsumoto et al. 2006; Bertrang et al. 2017; Bertrang & Wolf 2017). Both aligned and misaligned orientation of the outflow axes along the magnetic field directions have been reported (e.g. Jones & Amini 2003; Bertrang et al. 2014; Hull et al. 2014; Zhang et al. 2014). Based on MHD simulations, Matsumoto et al. (2006) find that the alignment degree between outflow and magnetic field is directly correlated to the magnetic field strength: the stronger the magnetic field, the better the alignment. In the following, we examine the relative position of the outflow axes and the magnetic field directions for CB34, using our near-IR data and compare our result to the findings in the optical and sub-mm wavelength range by Das et al. (2016).

In discussing the relative orientation between the CO outflow and the magnetic field, we must consider that only one component of the spatial orientation of the outflow (v sin i) is known from velocity measurements, so projectional effects have to be considered, as well as that polarization segments only trace the magnetic field structure projected on the plane of sky. We assumed that the magnetic field is orientated parallel to the polarization segments in the near-IR wavelength range, and perpendicular to the polarization segments in the sub-mm wavelength range.

The orientation of the CO outflow of CB34 is roughly parallel to the axis linking the cores C1 and C2 (Khanzadyan et al. 2002), indicating a north-east to south-west orientation.

The orientation of polarization segments in the near-IR wavelength range, tracing the less-dense outer parts of the globule, is almost perpendicular to the direction of the outflow. This fits well to the findings in the optical (Das et al. 2016). However, the polarization segments in the sub-mm polarization range, tracing the dense centre of the globule, are roughly aligned to the outflow direction, especially for those related to core C1 (Das et al. 2016).

Our new data are consistent with the finding of a correlation between the magnetic field orientation and the direction of the outflow, depending on the density and distance to the protostellar core.

6.2 Magnetic fields in star formation as a function of mass

Here, we present a sample of the lowest mass star-forming clouds that are known, located completely isolated. These low-mass clouds appear to be dominated by magnetic fields (see also, e.g., Alves et al. 2011; Bertrang et al. 2014). In larger but still relatively small clouds, such as Taurus, the cloud structure can be explained at least proximately in terms of turbulent simulations (e.g. Ostriker et al. 2001; Offner et al. 2008; Kritsuk et al. 2013). These intermediate mass clouds have comparatively short lifetimes of ~ 3 Myr (e.g. Hartmann et al. 2001, 2016; Ballesteros-Paredes & Hartmann 2007). In more massive clouds, for example Orion (mass ~ 105M), recent work seems to indicate that clouds may have lifetimes that are an order of magnitude longer (~ 30 Myr; e.g. Murray 2011; Meidt et al. 2015), and be subject to the influence of “dynamically relevant” magnetic fields (e.g. Commerçon et al. 2011; Palau et al. 2013; Stutz & Gould 2016; Busquet et al. 2016; Robitaille et al. 2018; Stutz et al. 2018). In this picture, in order to reconcile the above results including globules, we postulate a mass dependent role of magnetic fields in clouds, whereby magnetic fields appear to dominate at low and high masses, but perhaps are rather sub-dominant at intermediate (e.g. Taurus) masses. To investigate why the globules present magnetic field-dominated conditions, in the absence of more rigorous age or timescale estimates, requires more detailed modelling of the origin of globules. However, it is consistent with Fig. 1 of Stutz & Gould (2016) which strongly suggests that at least in the case of B 68 it originated as part of an undulating filamentary structure which then fragmented into cores such as B 68. This can only be explained by magnetic fields (and not gravity, due to the undulation). The fact that a reasonable fraction of such globules present “filamentary tails” (e.g. Launhardt et al. 2013) in sub-mm data is consistent with this scenario. The fact that not all globules present such “vestigial tails” may also indicate that globules have comparatively long lifetimes, which would then be consistent with both the fact that they are observed in the starless phases (indicating slow evolution that would then allow for the dissipation of the tails below detectable levels in the outer less dense regions), and the fact that the magnetic fields are observed to be dominant. This tentative scenario must be tested, but we note that proposed reasons for why magnetic fields maybe dominant in high mass clouds are different than for globules. For example, Stutz & Gould (2016) propose that massive clouds evolve under their own gravity, inevitably then concentrating the magnetic fields in an initially turbulent cloud medium (see e.g. the differences between the ISF and L1641 within Orion, filaments which both have the same large scale gravitational potential but different inner gravity and density profiles). In any case, this ansatz should be tested, but for the time being is consistent with observations. Irrespective of the details of internal evolution of individual systems, the role of magnetic fields in clouds spanning between approximately five to six orders of magnitude in mass has not been previously discussed, and thus must be tested with future theoretical and observational work.

7 Conclusions

We have obtained, for the first time, near-IR polarization data and, in conjunction with archival optical and sub-mm data, compiled multi-wavelength polarization maps which allow for the verification of the magnetic field influence on scales of 103 106 (~ 0.0055pc) in the Bok globules CB34, CB56, and [OMK2002]18, covering optically thin and optically thick regions. The major results are as follows:

  • 1.

    Our findings determine the distance of CB56 to ≤ 200 pc, for the first time, by applying a comparison with measurements of the interstellar polarization within a distance of 200 pc (see Sect. 2).

  • 2.

    We find a strong polarization signal of several percents in all three Bok globules in the near-IR.

  • 3.

    In CB34, the polarization segments in the near-IR are well aligned with those obtained in the optical tracing scales of 105–106 au (~0.55pc) as well as with those obtained in the sub-mm which trace the innermost, densest part of the globule on scales of 103 au (~0.005pc).

  • 4.

    In CB56, the polarization segments in the near-IR are well aligned with those obtained in the optical wavelength range on scales of 104–105 au (~0.05–0.5pc).

  • 5.

    In [OMK2002]18, the polarization segments in the western region are well ordered and show high polarization degrees, while we find indications for depolarizing effects in the eastern region.

  • 6.

    We find a correlation between the magnetic field structure and the CO outflow of the Bok globule CB34, depending on the density and distance to the protostellar core.

  • 7.

    For CB34, comparable magnetic field strengths are found in the regions traced in the near-IR resp. sub-mm observations by applying the CF method.

  • 8.

    We find strong indications for dominant magnetic fields in all of these three globules.

Acknowledgements

The authors thank George A. Gontcharov for providing the data on interstellar polarization in advance to publication. We thank Amelia M. Stutz for valuable discussions which improved this paper. We thank the anonymous referee for a critical and helpful report. G.H.M.B. acknowledges financial support from CONICYT through FONDECYT grant 3170657. This project has also received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 757957).

Appendix A Correction factor derivation for Wollaston prism

As stated in Sect. 4, a correction factor, Cλ for the Wollaston prism must be calculated to correct for the instrumental polarization, as the transmission ratio of the prism is not ideal. This correction factor is applied to the intensities of the upper, iu, and lower, il, beams created by the Wollaston prism in the following way (Wolf et al. 2002): iu,2=iu,2*Cλ,il,2=il,1.\begin{equation*} \begin{aligned} & i_{u,2} = i_{u,2}*C_{\lambda}, \\ & i_{l,2} = i_{l,1}. \end{aligned} \end{equation*}(A.1)

In our previous study, we find that the instrumental polarization for an instrument in Nasmyth focus has a strong dependency of the elevation angle of the telescope (ISAAC/VLT; Bertrang et al. 2014). For the here presented study, we test SOFI/NTT for this dependency as well. To determine the instrumental polarization, Cλ has to be adjusted in such a way that the measured polarization of the unpolarized standard star is minimal.

From November 11 to 15, 2015, we observed unpolarized standard stars (see Table 1) and derive Cλ dependent on the elevation angle of the telescope, with respect to the airmass of the object (see Fig. A.1). We find no correlation between the correction factor, Cλ, and the elevation angle of SOFI/NTT. Hence, we use the mean value of the correction factor to reduce our data. We observed several standard stars with SOFI/NTT over only a few nights. Observations of the same target over a longer period, as performed in Bertrang et al. (2014), might result in a different finding.

thumbnail Fig. A.1

Comparison between the correction factor that minimizes the instrumental polarization and the air mass of the observed unpolarized standard stars. The dashed line corresponds to the mean value of the computed correction factors. In contrast to ISAAC/VLT (Bertrang et al. 2014), we do not find a correlation between the instrumental polarization and the elevation angle of the telescope, respectively, the airmass of the object.

References

  1. Alves, J., Hartmann, L., Briceno, C., & Lada, C. J. 1997, AJ, 113, 1395 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alves, F. O., Franco, G. A. P., & Girart, J. M. 2008, A&A, 486, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Alves, F. O., Acosta-Pulido, J. A., Girart, J. M., Franco, G. A. P., & López, R. 2011, AJ, 142, 33 [NASA ADS] [CrossRef] [Google Scholar]
  4. Alves, F. O., Frau, P., Girart, J. M., et al. 2014, A&A, 569, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Ballesteros-Paredes, J., & Hartmann, L. 2007, Rev. Mex. Astron. Astrofis, 43, 123 [NASA ADS] [Google Scholar]
  6. Bertrang, G. H.-M., & Wolf, S. 2017, MNRAS, 469, 2869 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bertrang, G., Wolf, S., & Das, H. S. 2014, A&A, 565, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bertrang, G. H. M., Flock, M., & Wolf, S. 2017, MNRAS, 464, L61 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bok, B. J., & Reilly, E. F. 1947, ApJ, 105, 255 [NASA ADS] [CrossRef] [Google Scholar]
  10. Brauer, R., Wolf, S., & Reissl, S. 2016, A&A, 588, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Busquet, G., Estalella, R., Palau, A., et al. 2016, ApJ, 819, 139 [NASA ADS] [CrossRef] [Google Scholar]
  12. Chakraborty, A., Das, H. S., & Paul, D. 2014, MNRAS, 442, 479 [NASA ADS] [CrossRef] [Google Scholar]
  13. Chandrasekhar, S., & Fermi, E. 1953, ApJ, 118, 113 [NASA ADS] [CrossRef] [Google Scholar]
  14. Clemens, D. P. 2012, ApJ, 748, 18 [NASA ADS] [CrossRef] [Google Scholar]
  15. Clemens, D. P., & Barvainis, R. 1988, ApJS, 68, 257 [NASA ADS] [CrossRef] [Google Scholar]
  16. Clemens, D. P., Yun, J. L., & Heyer, M. H. 1991, ApJS, 75, 877 [NASA ADS] [CrossRef] [Google Scholar]
  17. Codella, C., & Scappini, F. 2003, MNRAS, 344, 1257 [NASA ADS] [CrossRef] [Google Scholar]
  18. Commerçon, B., Hennebelle, P., & Henning, T. 2011, ApJ, 742, L9 [NASA ADS] [CrossRef] [Google Scholar]
  19. Das, A., Das, H. S., Medhi, B. J., & Wolf, S. 2016, Ap&SS, 361, 381 [NASA ADS] [CrossRef] [Google Scholar]
  20. Draine, B. T., & Weingartner, J. C. 1997, ApJ, 480, 633 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dutra, C. M., & Bica, E. 2002, A&A, 383, 631 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Federrath, C. 2015, MNRAS, 450, 4035 [NASA ADS] [CrossRef] [Google Scholar]
  23. Federrath, C. 2016, J. Phys. Conf. Ser., 719, 012002 [CrossRef] [Google Scholar]
  24. Fossati, L., Bagnulo, S., Mason, E., & Landi Degl’Innocenti E. 2007, in The Future of Photometric, Spectrophotometric and Polarimetric Standardization, ed. C. Sterken, ASP Conf. Ser., 364, 503 [NASA ADS] [Google Scholar]
  25. Franco, G. A. P., Alves, F. O., & Girart, J. M. 2010, ApJ, 723, 146 [Google Scholar]
  26. Galli, D. 2009, Mem. Soc. Astron. It., 80, 54 [NASA ADS] [Google Scholar]
  27. Girart, J. M., Rao, R., & Marrone, D. P. 2008, Ap&SS, 313, 87 [NASA ADS] [CrossRef] [Google Scholar]
  28. Gómez, J. F., de Gregorio-Monsalvo, I., Suárez, O., & Kuiper, T. B. H. 2006, AJ, 132, 1322 [NASA ADS] [CrossRef] [Google Scholar]
  29. Gontcharov, G. 2017, in Stars: From Collapse to Collapse, eds. Y. Y. Balega, D. O. Kudryavtsev, I. I. Romanyuk, & I. A. Yakunin, ASP Conf. Ser., 510, 78 [NASA ADS] [Google Scholar]
  30. Hartmann, L., Ballesteros-Paredes, J., & Bergin, E. A. 2001, ApJ, 562, 852 [NASA ADS] [CrossRef] [Google Scholar]
  31. Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135 [NASA ADS] [CrossRef] [Google Scholar]
  32. Hennebelle, P., & Ciardi, A. 2009, A&A, 506, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Hoang, T., Lazarian, A., & Yan, H. 2007, BAAS, 39, 984 [NASA ADS] [Google Scholar]
  34. Huard, T. L., Weintraub, D. A., & Sandell, G. 2000, A&A, 362, 635 [NASA ADS] [Google Scholar]
  35. Hull, C. L. H., Plambeck, R. L., Kwon, W., et al. 2014, ApJS, 213, 13 [NASA ADS] [CrossRef] [Google Scholar]
  36. Jones, T. J., & Amini, H. 2003, AJ, 125, 1418 [NASA ADS] [CrossRef] [Google Scholar]
  37. Kandori, R., Tamura, M., Kusakabe, N., et al. 2007, PASJ, 59, 487 [NASA ADS] [Google Scholar]
  38. Khanzadyan, T., Smith, M. D., Gredel, R., Stanke, T., & Davis, C. J. 2002, A&A, 383, 502 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Kritsuk, A. G., Lee, C. T., & Norman, M. L. 2013, MNRAS, 436, 3247 [NASA ADS] [CrossRef] [Google Scholar]
  40. Launhardt, R., & Henning, T. 1997, A&A, 326, 329 [NASA ADS] [Google Scholar]
  41. Launhardt, R., Evans, II, N. J., Wang, Y., et al. 1998, ApJS, 119, 59 [NASA ADS] [CrossRef] [Google Scholar]
  42. Launhardt, R., Nutter, D., Ward-Thompson, D., et al. 2010, ApJS, 188, 139 [NASA ADS] [CrossRef] [Google Scholar]
  43. Launhardt, R., Stutz, A. M., Schmiedeke, A., et al. 2013, A&A, 551, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Lazarian, A. 2007, J. Quant. Spectr. Rad. Transf., 106, 225 [NASA ADS] [CrossRef] [Google Scholar]
  45. Lizano, S., & Galli, D. 2015, in Magnetic Fields in Diffuse Media, eds. A. Lazarian, E. M. de Gouveia Dal Pino, & C. Melioli, Astrophys. Space Sci. Lib., 407, 459 [NASA ADS] [CrossRef] [Google Scholar]
  46. Maheswar, G., Lee, C. W., & Dib, S. 2011, A&A, 536, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Matsumoto, T., Nakazato, T., & Tomisaka, K. 2006, ApJ, 637, L105 [NASA ADS] [CrossRef] [Google Scholar]
  48. Matthews, B. C., McPhee, C. A., Fissel, L. M., & Curran, R. L. 2009, ApJS, 182, 143 [NASA ADS] [CrossRef] [Google Scholar]
  49. Meidt, S. E., Hughes, A., Dobbs, C. L., et al. 2015, ApJ, 806, 72 [NASA ADS] [CrossRef] [Google Scholar]
  50. Murray, N. 2011, ApJ, 729, 133 [NASA ADS] [CrossRef] [Google Scholar]
  51. Myers, P. C., & Goodman, A. A. 1991, ApJ, 373, 509 [NASA ADS] [CrossRef] [Google Scholar]
  52. Offner, S. S. R., Krumholz, M. R., Klein, R. I., & McKee, C. F. 2008, AJ, 136, 404 [NASA ADS] [CrossRef] [Google Scholar]
  53. Onishi, T., Mizuno, A., Kawamura, A., Tachihara, K., & Fukui, Y. 2002, ApJ, 575, 950 [NASA ADS] [CrossRef] [Google Scholar]
  54. Ostriker, E. C., Stone, J. M., & Gammie, C. F. 2001, ApJ, 546, 980 [NASA ADS] [CrossRef] [Google Scholar]
  55. Palau, A., Fuente, A., Girart, J. M., et al. 2013, ApJ, 762, 120 [NASA ADS] [CrossRef] [Google Scholar]
  56. Papoular, R. 2018, MNRAS, 479, 1685 [NASA ADS] [CrossRef] [Google Scholar]
  57. Paul, D., Das, H. S., & Sen, K. 2012, Bull. Astron. Soc. India, 40, 113 [NASA ADS] [Google Scholar]
  58. Peters, T., Banerjee, R., Klessen, R. S., & Mac Low M.-M. 2011, ApJ, 729, 72 [NASA ADS] [CrossRef] [Google Scholar]
  59. Robitaille, J. F., Scaife, A. M. M., Carretti, E., et al. 2018, A&A, 617, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Sadavoy, S. I., Myers, P. C., Stephens, I. W., et al. 2018, ApJ, 859, 165 [NASA ADS] [CrossRef] [Google Scholar]
  61. Sepúlveda, I., Anglada, G., Estalella, R., et al. 2011, A&A, 527, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Serkowski, K., Mathewson, D. S., & Ford, V. L. 1975, ApJ, 196, 261 [NASA ADS] [CrossRef] [Google Scholar]
  63. Stutz, A. M., & Gould, A. 2016, A&A, 590, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Stutz, A. M., Bieging, J. H., Rieke, G. H., et al. 2007, ApJ, 665, 466 [NASA ADS] [CrossRef] [Google Scholar]
  65. Stutz, A. M., Gonzalez-Lobos, V. I., & Gould, A. 2018, MNRAS, submitted [arXiv:1807.11496v1] [Google Scholar]
  66. Sugitani, K., Nakamura, F., Tamura, M., et al. 2010, ApJ, 716, 299 [NASA ADS] [CrossRef] [Google Scholar]
  67. Turnshek, D. A., Bohlin, R. C., Williamson, II, R. L., et al. 1990, AJ, 99, 1243 [NASA ADS] [CrossRef] [Google Scholar]
  68. Wang, Y., Evans, II, N. J., Zhou, S., & Clemens, D. P. 1995, ApJ, 454, 217 [NASA ADS] [CrossRef] [Google Scholar]
  69. Weintraub, D. A., Goodman, A. A., & Akeson, R. L. 2000, Protostars and Planets IV, 247 [Google Scholar]
  70. Whittet, D. C. B., Martin, P. G., Hough, J. H., et al. 1992, ApJ, 386, 562 [NASA ADS] [CrossRef] [Google Scholar]
  71. Wolf, S., Vanzi, L., & Ageorges, N. 2002, How to Reduce SOFI/NTT Polarimetric Data, https://www.eso.org/sci/facilities/lasilla/instruments/sofi/tools/report.ps [Google Scholar]
  72. Yun, J. L., & Clemens, D. P. 1994, ApJS, 92, 145 [NASA ADS] [CrossRef] [Google Scholar]
  73. Zhang, Q., Qiu, K., Girart, J. M., et al. 2014, ApJ, 792, 116 [NASA ADS] [CrossRef] [Google Scholar]
  74. Zweibel, E. G. 1990, ApJ, 362, 545 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Polarization standard stars.

Table 2

Gas densities, gas velocities, polarization, and magnetic field strengths of the regions of CB34 traced in the near-IR.

All Figures

thumbnail Fig. 1

Panel a: near-IR polarization segments (red) and optical polarization segments (blue) are plotted over a 15′ × 15′ R-band DSS image of the field containing CB34. A 10% polarization segment is drawn as reference for both data sets into the lower left corner. Contours correspond to JCMT/SCUBA 850 μm dust continuum emissions from 0.1 to 0.7 Jy beam−1, with intervals of 0.1 Jy beam−1 (Das et al. 2016). Panel b: zoom into the dense centre of CB34. The near-IR polarization segments (red) and the archival sub-mm polarization segments (blue) are plotted over the total continuum intensity map at 850 (Das et al. 2016). The sub-mm segments are binned over a 10″ × 10″ grid. The solid (dashed) black (grey) contour lines correspond to blueshifted (redshifted) integrated CO emission. Contour levels are spaced at 0.1 K km s−1 intervals of 0.5 K km s−1 (black) and 0.2 K km s−1 intervals of 0.8 K km s−1 (white), respectively (from Yun & Clemens 1994; beam size of 48′′). A 5% polarization segment is drawn as reference for both data sets into the lower right corner. The well-ordered polarization pattern indicates dominant magnetic fields on scales of 103 –106 au (~ 0.0055pc).

In the text
thumbnail Fig. 2

Panel a: near-IR polarization segments plotted over a 15′× 15′ R-Band DSSimage of the field containing CB56. A 10% polarization segment is drawn for reference in the lower left corner. Only the segments with Pσp ≥ 3 are plotted. Panel b: optical polarization segments of CB56, plotted over a 10′ × 10′ R-Band DSS image (Chakraborty et al. 2014). A 1% polarization segment is drawn for reference in the lower left corner. The near-IR and optical polarization is strong and well ordered, indicating dominant magnetic fields on scales of 104 105 au (~ 0.050.5pc).

In the text
thumbnail Fig. 3

Panel a: near-IR polarization segments plotted over a ~20′× 20′ R-band DSS image of the field containing [OMK2002]18. A 10% polarization segment is drawn for reference in the lower left corner. Only the segments with Pσp ≥ 3 are plotted. Panel b: sub-mm polarization segments of [OMK2002]18, plotted over a 20″ × 20″ grid. Contours correspond to JCMT/SCUBA 850 μm dust continuum emissions from 0.05 to 0.25 Jy beam−1, with intervals of 0.05 Jy beam−1 (Matthews et al. 2009; ©AAS, reproduced with permission). Only segments with Pσp ≥ 2 are plotted. The near-IR polarization is strong and well ordered, indicating dominant magnetic fields on scales of 104 –105 au (~ 0.05–0.5pc) in the western region of the globule.

In the text
thumbnail Fig. 4

Distribution of polarization degree, P for P ≥ 3%, counts given by N(P), of CB34, CB56, and OMK[2002]18, observed in the near-IR.

In the text
thumbnail Fig. 5

Distribution of polarization angle, γ, counts given by N(γ) of CB34, CB56, and OMK[2002]18, observed in the near-IR. The dashed lines represent the mean polarization angles, while dotted lines represent the corresponding 1σ deviation.

In the text
thumbnail Fig. A.1

Comparison between the correction factor that minimizes the instrumental polarization and the air mass of the observed unpolarized standard stars. The dashed line corresponds to the mean value of the computed correction factors. In contrast to ISAAC/VLT (Bertrang et al. 2014), we do not find a correlation between the instrumental polarization and the elevation angle of the telescope, respectively, the airmass of the object.

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.