P. C. Myers

Dense cores in dark clouds. VIII: Velocity gradients

An analysis of motions consistent with uniform rotation in dense cores is presented. Twenty-nine of the 43 cores studied have a statistically significant gradient. Some gradients are spatially continuous and are consistent with uniform rotation, but other apparent gradients are caused by clump-clump motion, or sharp localized gradients, within a map. The motions in L1495, B217, L1251, L43, B361, and L1551 are discussed in detail. In L1551, the residuals of the fit to the NH3 velocity field indicate an outflow from IRS5 in the same direction as the CO outflow. Gradient orientation appears to be preserved over a range of density, as evidenced by comparing results of NH3 to fits of (C-18)O and CS maps. The specific angular momentum is found to scale roughly as F exp 3/2, where R represents the diameter of the FWHM intensity contour in a map.

Systematic Molecular Differentiation in Starless Cores

We present evidence that low-mass starless cores, the simplest units of star formation, are systematically differentiated in their chemical composition. Some molecules, including CO and CS, almost vanish near the core centers, where the abundance decreases by at least 1 or 2 orders of magnitude with respect to the value in the outer core. At the same time, the N2H+ molecule has a constant abundance, and the fraction of NH3 increases toward the core center. Our conclusions are based on a systematic study of five mostly round starless cores (L1498, L1495, L1400K, L1517B, and L1544), which we have mapped in C18O (1-0), CS (2-1), N2H+ (1-0), NH3 (1, 1) and (2, 2), and the 1.2 mm continuum [complemented with C17O (1-0) and C34S (2-1) data for some systems]. For each core we have built a spherically symmetric model in which the density is derived from the 1.2 mm continuum, the kinetic temperature is derived from NH3, and the abundance of each molecule is derived using a Monte Carlo radiative transfer code, which simultaneously fits the shape of the central spectrum and the radial profile of integrated intensity. Regarding the cores for which we have C17O (1-0) and C34S (2-1) data, the model fits these observations automatically when the standard isotopomer ratio is assumed. As a result of this modeling, we also find that the gas kinetic temperature in each core is constant at approximately 10 K. In agreement with previous work, we find that if the dust temperature is also constant, then the density profiles are centrally flattened, and we can model them with a single analytic expression. We also find that for each core the turbulent line width seems constant in the inner 0.1 pc. The very strong abundance drop of CO and CS toward the center of each core is naturally explained by the depletion of these molecules onto dust grains at densities of (2-6) × 104 cm-3. N2H+ seems unaffected by this process up to densities of several times 105 cm-3, or even 106 cm-3, while the NH3 abundance may be enhanced by its lack of depletion and by reactions triggered by the disappearance of CO from the gas phase. With the help of the Monte Carlo modeling, we show that chemical differentiation automatically explains the discrepancy between the sizes of CS and NH3 maps, a problem that has remained unexplained for more than a decade. Our models, in addition, show that a combination of radiative transfer effects can give rise to the previously observed discrepancy in the line width of these two tracers. Although this discrepancy has been traditionally interpreted as resulting from a systematic increase of the turbulent line width with radius, our models show that it can arise in conditions of constant gas turbulence.

A survey for dense cores in dark clouds

A total of 149 dark cloud positions were surveyed for evidence of dense cores in the (J,K) = (1,1) rotating inversion line of NH3. Clouds with strong emission were mapped to determine the position of the core, and maps of 41 cores are presented. The spectrum of the peak emission is fitted by least-squares analysis to determine the optical depth, velocity, intrinsic line width, and excitation temperature. Statistical equilibrium analysis is used to determine the density of the core and the kinetic temperature when possible. Most of the dense cores have temperatures of 10-15 K, densities of 2000-20,000/cu cm, and intrinsic linewidths of 0.2-0.9 km/s. The core masses range from about 0.5 solar in L1517B to 760 solar in L1031B, and their sizes range from 0.06 to 0.9 pc. The cores are not generally spherically shaped, wtih aspect ratios ranging from 1.1 to 4 4. Cores with stars have broader lines than cores without stars. 74 refs.

On the internal structure of starless cores - I. Physical conditions and the distribution of CO, CS, N

We have characterized the physical structure and chemical composition of two close-to-round starless cores in Taurus-Auriga, L1498 and L1517B. Our analysis is based on high angular resolution observations in at least two transitions of NH3 ,N 2H + ,C S, C 34 S, C 18 O, and C 17 O, together with maps of the 1.2 mm continuum. For both cores, we derive radial profiles of constant temperature and constant turbulence, together with density distributions close to those of non-singular isothermal spheres. Using these physical conditions and a Monte Carlo radiative transfer model, we derive abundance profiles for all species and model the strong chemical differentiation of the core interiors. According to our models, the NH3 abundance increases toward the core centers by a factor of several (≈5) while N2H + has a constant abundance over most of the cores. In contrast, both C 18 O and CS (and isotopomers) are strongly depleted in the core interiors, most likely due to their freeze out onto grains at densities of a few 10 4 cm −3 . Concerning the kinematics of the dense gas, we find (in addition to constant turbulence) a pattern of internal motions at the level of 0.1 km s −1 . These motions seem correlated with asymmetries in the pattern of molecular depletion, and we interpret them as residuals of core contraction. Their distribution and size suggest that core formation occurs in a rather irregular manner and with a time scale of a Myr. A comparison of our derived core properties with those predicted by supersonic turbulence models of core formation shows that our Taurus cores are much more quiescent than representative predictions from these models. In two appendices at the end of the paper we present a simple and accurate approximation to the density profile of an isothermal (Bonnor-Ebert) sphere, and a Monte Carlo-calibrated method to derive gas kinetic temperatures using NH3 data.

\mathsf{_2}

We present a database of 264 cores mapped in the (J,K) = (1,1) and (2,2) lines of NH3. We list the core gas properties?peak positions, total ammonia column densities, intrinsic line widths, kinetic temperatures, volume densities, core sizes, aspect ratios, and velocity gradients, as well as the properties of associated young stellar objects (YSOs)?associated IRAS sources along with their luminosities and core-YSO distances, outflow velocities, and SIMBAD and cluster associations. We also present the results of our statistical analysis and enumerate important pairwise correlations among the various gas and YSO properties. The results indicate that the association of stellar clusters with star-forming cores has a greater impact on their properties than does the presence of associated YSOs within these cores, although the latter influence is also statistically significant. In other words, the difference in core properties (nonthermal line widths, kinetic temperatures, and core sizes) between cores with and without associated YSOs is less significant when compared with the difference in these properties between cores with and without cluster associations. Furthermore, core gas and YSO properties show a significant dependence on the star-forming region in which the core is located. For instance, cores in Orion have larger line widths, higher kinetic temperatures, and larger sizes compared with cores in Taurus. Similarly, YSOs in Orion are more luminous than those in Taurus. These cluster and regional dependencies seem important enough that they ought to be accounted for in any self-consistent theory of star formation. Finally, the ratio of starless to stellar cores is too small (8:12 in Taurus, 2:41 in Orion A) to be consistent with ambipolar diffusion timescales that predict ratios as high as 3-30. This result is true even for regions that are known to be well surveyed and not to suffer from significant sample biases.

H

The maps presented in Paper I are here used to infer the variation of the column densities of HCO+, DCO+, N2H+, and N2D+ as a function of distance from the dust peak. These results are interpreted with the aid of a crude chemical model that predicts the abundances of these species as a function of radius in a spherically symmetric model with radial density distribution inferred from the observations of dust emission at millimeter wavelengths and dust absorption in the infrared. Our main observational finding is that the N(N2D+)/N(N2H+) column density ratio is of order 0.2 toward the L1544 dust peak as compared to N(DCO+)/N(HCO+) = 0.04. We conclude that this result, as well as the general finding that N2H+ and N2D+ correlate well with the dust, is caused by CO being depleted to a much higher degree than molecular nitrogen in the high-density core of L1544. Depletion also favors deuterium enhancement, and thus N2D+, which traces the dense and highly CO-depleted core nucleus, is much more enhanced than DCO+. Our models do not uniquely define the chemistry in the high-density depleted nucleus of L1544, but they do suggest that the ionization degree is a few times 10-9 and that the ambipolar diffusion timescale is locally similar to the free-fall time. It seems likely that the lower limit, which one obtains to ionization degree by summing all observable molecular ions, is not a great underestimate of the true ionization degree. We predict that atomic oxygen is abundant in the dense core and, if so, H3O+ may be the main ion in the central highly depleted region of the core.

\mathsf{^+}

We report observations of 47 candidate protostars in two optically thick lines [H2CO (212-111) and CS (2-1)] and one optically thin line [N2H+ (1-0)] using the IRAM 30 m, SEST 15 m, and Haystack 37 m radio telescopes. The sources were selected for the redness of their spectra (Tbol < 200 K) and their near distance (d < 400 pc). Most of the sources have asymmetric, optically thick lines. The observed distribution of velocity differences, ?V = (Vthick - Vthin)/?Vthin, is skewed toward negative (blueshifted) velocities for both the H2CO and CS samples. This excess is much more significant for class 0 than for class I sources, suggesting that we detect infall motions toward class 0 and not toward class I sources. This indicates a difference in the physical conditions in the circumstellar envelopes around class I and class 0 sources, but it does not rule out the presence of infall onto class I sources by, for example, lower opacity gas. Bipolar outflows alone, or rotation alone, cannot reproduce these statistics if the sample of sources has randomly oriented symmetry axes. We identify 15 spectroscopic infall candidates, six of which are new. Most of these infall candidates have primarily turbulent rather than thermal motions and are associated with clusters rather than being isolated.

, and NH

We present an unbiased census of starless cores in Perseus, Serpens, and Ophiuchus, assembled by comparing large-scale Bolocam 1.1 mm continuum emission maps with Spitzer c2d surveys. We use the c2d catalogs to separate 108 starless from 92 protostellar cores in the 1.1 mm core samples from Enoch and Young and their coworkers. A comparison of these populations reveals the initial conditions of the starless cores. Starless cores in Perseus have similar masses but larger sizes and lower densities on average than protostellar cores, with sizes that suggest density profiles substantially flatter than ρ∝r^-2. By contrast, starless cores in Serpens are compact and have lower masses than protostellar cores; future star formation will likely result in lower mass objects than the currently forming protostars. Comparison to dynamical masses estimated from the NH3 survey of Perseus cores by Rosolowsky and coworkers suggests that most of the starless cores are likely to be gravitationally bound, and thus prestellar. The combined prestellar core mass distribution includes 108 cores and has a slope of α = -2.3 ± 0.4 for M > 0.8 M☉. This slope is consistent with recent measurements of the stellar initial mass function, providing further evidence that stellar masses are directly linked to the core formation process. We place a lower limit on the core-to-star efficiency of 25%. There are approximately equal numbers of prestellar and protostellar cores in each cloud; thus the dense prestellar core lifetime must be similar to the lifetime of embedded protostars, or 4.5 x 10^5 yr, with a total uncertainty of a factor of 2. Such a short lifetime suggests a dynamic, rather than quasi-static, core evolution scenario, at least at the relatively high mean densities (n > 2 x 10^4 cm^-3) to which we are sensitive.

\mathsf{_3}

We present a multiline study of the dense core L1544 in the Taurus molecular complex. Although L1544 does not harbor an embedded star, it presents several characteristics of cores that have already undergone star formation, suggesting that it may be rather advanced in its evolution toward becoming a star-forming core. The spectral lines from L1544 present an interesting dichotomy, with the thick dense gas tracers su†ering very strong self absorption while CO and its isotopes are not being absorbed at all. The presence of the self absorptions allows us to study both the density structure and kinematics of the gas in detail. A simple analysis shows that the core is almost isothermal and that the self absorptions are due to very subthermal excitation of the dense gas tracers in the outer layers. The density has to decrease outward rapidly, and a detailed radiative transfer calculation that simultaneously -ts three iso- topes of CO and two of CS shows that the density approximately follows a r~1.5 power law. The self absorptions, in addition, allow us to measure the relative velocity between the inner and outer layers of the core, and we -nd that there is a global pattern of inward motions (background and foreground approaching each other). The relative speed between the foreground and background changes with posi- tion, and we use a simple two-layer model to deduce that while the foreground gas has a constant veloc- ity, the background material presents systematic velocity changes that we interpret as arising from two velocity components. We explore the origin of the inward motions by comparing our observations with models of gravitational collapse. A model in which the infall starts at the center and propagates outward (as in the inside-out collapse of Shu) is inconsistent with the large extension of the absorption (that sug- gests an advanced age) and the lack of a star at the core center (that suggests extreme youth). Ambipolar di†usion seems also ruled out because of the large amount of the inward speed (up to 0.1 km s~1) and the fact that ionized species move with speeds similar to those of the neutrals. Other infall models seem also to have problems -tting the data, so if L1544 is infalling, it seems to be doing so in a manner not contemplated by the standard theories of star formation. Our study of L1544 illustrates how little is still known about the physical conditions that precede star formation and how detailed studies of starless cores are urgently needed. Subject headings: ISM: individual (L1544) E ISM: kinematics and dynamics E stars: formation

in L1498 and L1517B

This paper reports observations of 34 IRAS sources associated with dense cores in dark clouds at wavelengths 0.4-20 microns. The stars near cores tend to be visible T Tauri stars, while stars in cores tend to have circumstellar extinction 30-90 mag and luminosity about one solar, similar to that of T Tauri stars. The typical highly obscured star is probably accompanied by a luminous structure of substellar temperature, such as a circumstellar disk. In Taurus-Auriga, stars in cores probably become visible T Tauri stars less than 100,000 yr after they become detectable by IRAS; i.e., after they attain luminosity greater than about 0.1 solar. This implies that they are extremely young and may still be accreting. 54 references.