OVRO N2H+ Observations of Class 0 Protostars: Constraints on the Formation of Binary Stars

We present the results of an interferometric study of the N2H+(1--0) emission from nine nearby, isolated, low-mass protostellar cores, using the OVRO millimeter array. The main goal of this study is the kinematic characterization of the cores in terms of rotation, turbulence, and fragmentation. Eight of the nine objects have compact N2H+ cores with FWHM radii of 1200 -- 3500 AU, spatially coinciding with the thermal dust continuum emission. The only more evolved (Class I) object in the sample (CB 188) shows only faint and extended N2H+ emission. The mean N2H+ line width was found to be 0.37 km/s. Estimated virial masses range from 0.3 to 1.2 M_sun. We find that thermal and turbulent energy support are about equally important in these cores, while rotational support is negligible. The measured velocity gradients across the cores range from 6 to 24 km/s/pc. Assuming these gradients are produced by bulk rotation, we find that the specific angular momenta of the observed Class 0 protostellar cores are intermediate between those of dense (prestellar) molecular cloud cores and the orbital angular momenta of wide PMS binary systems. There appears to be no evolution (decrease) of angular momentum from the smallest prestellar cores via protostellar cores to wide PMS binary systems. In the context that most protostellar cores are assumed to fragment and form binary stars, this means that most of the angular momentum contained in the collapse region is transformed into orbital angular momentum of the resulting stellar binary systems.Comment: 35 pages, 9 figures (one in color), 6 tables. Accepted by ApJ (to appear in Nov. 2007

Looking into the hearts of Bok globules: MM and submm continuum images of isolated star-forming cores

We present the results of a comprehensive infrared, submillimetre, and millimetre continuum emission study of isolated low-mass star-forming cores in 32 Bok globules, with the aim to investigate the process of star formation in these regions. The submillimetre and millimetre dust continuum emission maps together with the spectral energy distributions are used to model and derive the physical properties of the star-forming cores, such as luminosities, sizes, masses, densities, etc. Comparisons with ground-based near-infrared and space-based mid and far-infrared images from Spitzer are used to reveal the stellar content of the Bok globules, association of embedded young stellar objects with the submm dust cores, and the evolutionary stages of the individual sources. Submm dust continuum emission was detected in 26 out of the 32 globule cores observed. For 18 globules with detected (sub)mm cores we derive evolutionary stages and physical parameters of the embedded sources. We identify nine starless cores, most of which are presumably prestellar, nine Class 0 protostars, and twelve Class I YSOs. Specific source properties like bolometric temperature, core size, and central densities are discussed as function of evolutionary stage. We find that at least two thirds (16 out of 24) of the star-forming globules studied here show evidence of forming multiple stars on scales between 1,000 and 50,000 AU. However, we also find that most of these small prototstar and star groups are comprised of sources with different evolutionary stages, suggesting a picture of slow and sequential star formation in isolated globulesComment: 60 pages, 28 figures, accepted by The Astrophysical Journal Supplement Serie

Dust SEDs in the era of Herschel and Planck: a Hierarchical Bayesian fitting technique

We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values.Comment: 20 pages, 9 figures, ApJ format, revised version matches ApJ-accepted versio

Dust-temperature of an isolated star-forming cloud: Herschel observations of the Bok globule CB244

We present Herschel observations of the isolated, low-mass star-forming Bok globule CB244. It contains two cold sources, a low-mass Class 0 protostar and a starless core, which is likely to be prestellar in nature, separated by 90 arcsec (~ 18000 AU). The Herschel data sample the peak of the Planck spectrum for these sources, and are therefore ideal for dust-temperature and column density modeling. With these data and a near-IR extinction map, the MIPS 70 micron mosaic, the SCUBA 850 micron map, and the IRAM 1.3 mm map, we model the dust-temperature and column density of CB244 and present the first measured dust-temperature map of an entire star-forming molecular cloud. We find that the column-averaged dust-temperature near the protostar is ~ 17.7 K, while for the starless core it is ~ 10.6K, and that the effect of external heating causes the cloud dust-temperature to rise to ~ 17 K where the hydrogen column density drops below 10^21 cm^-2. The total hydrogen mass of CB244 (assuming a distance of 200 pc) is 15 +/- 5 M_sun. The mass of the protostellar core is 1.6 +/- 0.1 M_sun and the mass of the starless core is 5 +/- 2 M_sun, indicating that ~ 45% of the mass in the globule is participating in the star-formation process.Comment: Accepted for A&A Herschel Special Issue; 5 pages, 2 figure

SMA and Spitzer Observations of Bok Glouble CB17: A Candidate First Hydrostatic Core?

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.Comment: 33 pages, 11 figures, to be published by Ap

The L723 low mass star forming protostellar system: resolving a double core

We present 1.35 mm SMA observations around the low-mass Class 0 source IRAS 19156+1906, at the the center of the L723 dark cloud. We detected emission from dust as well as emission from H2CO, DCN and CN, which arise from two cores, SMA 1 and SMA 2, separated by 2.9" (880 AU). SMA 2 is associated with VLA 2. SiO 5-4 emission is detected, possibly tracing a region of interaction between the dense envelope and the outflow. We modeled the dust and the H2CO emission from the two cores: they have similar physical properties but SMA 2 has a larger p-H2CO abundance than SMA 1. The p-H2CO abundances found are compatible with the value of the outer part of the circumstellar envelopes associated with Class 0 sources. SMA 2 is likely more evolved than SMA 1. The kinematics of the two sources show marginal evidence of infall and rotation motions. The mass detected by the SMA observation, which trace scales of ~1000 AU, is only a small fraction of the mass contained in the large scale molecular envelope, which suggests that L723 is still in a very early phase of star formation. Despite the apparent quiescent nature of the L723, fragmentation is occurring at the center of the cloud at different scales. Thus, at 1000 AU the cloud has fragmented in two cores, SMA 1 and SMA 2. At the same time, at least one of these cores, SMA 2, has undergone additional fragmentation at scales of 150 AU, forming a multiple stellar system.Comment: 35 pages, 15 figures. Accepted to the Astrophysical Journa

The Earliest Phases of Star formation (EPoS): Temperature, density, and kinematic structure of the star-forming core CB 17

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.Comment: published in A&

CB17: Inferring the dynamical history of a prestellar core with chemo-dynamical models

We present a detailed theoretical study of the isolated Bok globule CB17 (L1389) based on spectral maps of CS, HCO$^+$, C$^{18}$O, C$^{34}$S, and H$^{13}$CO$^+$ lines. A phenomenological model of prestellar core evolution, a time-dependent chemical model, and a radiative transfer simulation for molecular lines are combined to reconstruct the chemical and kinematical structure of this core. We developed a general criterion that allows to quantify the difference between observed and simulated spectral maps. By minimizing this difference, we find that very high and very low values of the effective sticking probability $S$ are not appropriate for the studied prestellar core. The most probable $S$ value for CB17 is 0.3--0.5. The spatial distribution of the intensities and self-absorption features of optically thick lines is indicative of UV irradiation of the core. By fitting simultaneously optically thin and optically thick transitions, we isolate the model that reproduces all the available spectral maps to a reasonable accuracy. The line asymmetry pattern in CB17 is reproduced by a combination of infall, rotation, and turbulent motions with velocities $\sim0.05$ km s$^{-1}$, $\sim0.1$ km s$^{-1}$, and $\sim0.1$ km s$^{-1}$, respectively. These parameters corresponds to energy ratios $E_{\rm rot}/E_{\rm grav}\approx0.03$, $E_{\rm therm}/E_{\rm grav}\approx0.8$, and $E_{\rm turb}/E_{\rm grav}\approx0.05$ (the rotation parameters are determined for $i=90^\circ$). The chemical age of the core is about 2 Myrs. In particular, this is indicated by the central depletion of CO, CS, and HCO$^+$. Based on the angular momentum value, we argue that the core is going to fragment, i.e., to form a binary (multiple) star. (abridged)Comment: ApJ, in pres

The circumstellar disc in the Bok globule CB 26: Multi-wavelength observations and modelling of the dust disc and envelope

Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step in the planet formation process. The Bok globule CB 26 harbours such a young disc. We present a detailed model of the edge-on circumstellar disc and its envelope in the Bok globule CB 26. The model is based on HST near-infrared maps in the I, J, H, and K bands, OVRO and SMA radio maps at 1.1mm, 1.3mm and 2.7mm, and the spectral energy distribution (SED) from 0.9 microns to 3mm. New photometric and spectroscopic data from the Spitzer Space Telescope and the Caltech Submilimeter Observatory have been obtained and are part of our analysis. Using the self-consistent radiative transfer code MC3D, the model we construct is able to discriminate parameter sets and dust properties of both its parts, namely envelope and disc. We find that the disc has an inner hole with a radius of 45 +/- 5 AU. Based on a dust model including silicate and graphite the maximum grain size needed to reproduce the spectral millimetre index is 2.5 microns. Features seen in the near-infrared images, dominated by scattered light, can be described as a result of a rotating envelope. Successful employment of ISM dust in both the disc and envelope hint that grain growth may not yet play a significant role for the appearance of this system. A larger inner hole gives rise to the assumption that CB 26 is a circumbinary disc.Comment: 18 pages, 15 figures, Accepted for publication in A&

Optical and submillimetre observations of Bok globules -- tracing the magnetic field from low to high density

We present optical and submillimetre polarimetry data of the Bok globule CB3 and optical polarimetry data of the Bok globule CB246. We use each set of polarimetry data to infer the B-field orientation in each of the clouds. The optical data can only be used in the low density, low extinction edge regions of clouds. The submillimetre data can only be used in the high column-density, central regions of the clouds. It has previously been found that near-infrared polarisation mapping of background stars does not accurately trace the magnetic field in dense cloud regions. This may be due to a lack of aligned grains in dense regions. We test this by comparing the field orientations measured by our two independent polarimetry methods. We find that the field orientation deduced from the optical data matches up well with the orientation estimated from the submillimetre data. We therefore claim that both methods are accurately tracing the same magnetic field in CB3. Hence, in this case, there must be significant numbers of aligned dust grains in the high density region, and they do indeed trace the magnetic field in the submillimetre. We find an offset of 40$\pm$14 degrees between the magnetic field orientation and the short axis of the globule. This is consistent with the mean value of 31$\pm$3 degrees found in our previous work on prestellar cores, even though CB3 is a protostellar core. Taken together, the six prestellar cores that we have now studied in this way show a mean offset between magnetic field orientation and core short axis of $\sim30\pm$3 degrees, in apparent contradiction with some models of magnetically dominated star formation.Comment: 8 pages, 3 figures, accepted for publication in MNRA