2002ApJ...569..815T

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 N₂H⁺ molecule has a constant abundance, and the fraction of NH₃ 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 C¹⁸O (1-0), CS (2-1), N₂H⁺ (1-0), NH₃ (1, 1) and (2, 2), and the 1.2 mm continuum [complemented with C¹⁷O (1-0) and C³⁴S (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 NH₃, 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 C¹⁷O (1-0) and C³⁴S (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)x10⁴ cm^–3. N₂H⁺ seems unaffected by this process up to densities of several times 10⁵ cm^–3, or even 10⁶ cm^–3, while the NH₃ 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 NH₃ 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.