Oscillating Starless Cores: The Nonlinear Regime

In a previous paper, we modeled the oscillations of a thermally-supported (Bonnor-Ebert) sphere as non-radial, linear perturbations following a standard analysis developed for stellar pulsations. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. However, the linear analysis is unable to address several questions, among them the stability, and lifetime of the perturbations. In this paper we simulate the oscillations using a three-dimensional numerical hydrodynamic code. We find that the oscillations are damped predominantly by non-linear mode-coupling, and the damping time scale is typically many oscillation periods, corresponding to a few million years, and persisting over the inferred lifetime of gobules.Comment: 7 pages, 7 figures, accepted by Ap

Coreshine in L1506C - Evidence for a primitive big-grain component or indication for a turbulent core history?

The recently discovered coreshine effect can aid in exploring the core properties and in probing the large grain population of the ISM. We discuss the implications of the coreshine detected from the molecular cloud core L1506C in the Taurus filament for the history of the core and the existence of a primitive ISM component of large grains becoming visible in cores. The coreshine surface brightness of L1506C is determined from IRAC Spitzer images at 3.6 micron. We perform grain growth calculations to estimate the grain size distribution in model cores similar in gas density, radius, and turbulent velocity to L1506C. Scattered light intensities at 3.6 micron are calculated for a variety of MRN and grain growth distributions to compare with the observed coreshine. For a core with the overall physical properties of L1506C, no detectable coreshine is predicted for an MRN size distribution. Extending the distribution to grain radii of about 0.65 $\mu$m allows to reproduce the observed surface brightness level in scattered light. Assuming the properties of L1506C to be preserved, models for the growth of grains in cores do not yield sufficient scattered light to account for the coreshine within the lifetime of the Taurus complex. Only increasing the core density and the turbulence amplifies the scattered light intensity to a level consistent with the observed coreshine brightness. The grains could be part of primitive omni-present large grain population becoming visible in the densest part of the ISM, could grow under the turbulent dense conditions of former cores, or in L1506C itself. In the later case, L1506C must have passed through a period of larger density and stronger turbulence. This would be consistent with the surprisingly strong depletion usually attributed to high column densities, and with the large-scale outward motion of the core envelope observed today.Comment: 6 pages, 6 figures, accepted for publication in Astronomy & Astrophysic

Mid-infrared observations of the SGR 1900+14 error box

We report on mid-infrared observations of the compact stellar cluster located in the proximity of SGR 1900+14, and the radio/X-ray position of this soft-gamma repeater. Observations were performed in May and June of 2001 when the bursting source was in an active state. At the known radio and X-ray position of the SGR we did not detect transient mid-IR activity, although the observations were performed only hours before and after an outburst in the high-energy band.Comment: 4 pages, 3 figures, to appear in "Gamma-Ray Burst and Afterglow Astronomy 2001", Woods Hole; 5-9 Nov, 200

Formation and Collapse of Nonaxisymmetric Protostellar Cores in Planar Magnetic Interstellar Clouds: Formulation of the Problem and Linear Analysis

We formulate the problem of the formation and collapse of nonaxisymmetric protostellar cores in weakly ionized, self-gravitating, magnetic molecular clouds. In our formulation, molecular clouds are approximated as isothermal, thin (but with finite thickness) sheets. We present the governing dynamical equations for the multifluid system of neutral gas and ions, including ambipolar diffusion, and also a self-consistent treatment of thermal pressure, gravitational, and magnetic (pressure and tension) forces. The dimensionless free parameters characterizing model clouds are discussed. The response of cloud models to linear perturbations is also examined, with particular emphasis on length and time scales for the growth of gravitational instability in magnetically subcritical and supercritical clouds. We investigate their dependence on a cloud's initial mass-to-magnetic-flux ratio (normalized to the critical value for collapse), the dimensionless initial neutral-ion collision time, and also the relative external pressure exerted on a model cloud. Among our results, we find that nearly-critical model clouds have significantly larger characteristic instability lengthscales than do more distinctly sub- or supercritical models. Another result is that the effect of a greater external pressure is to reduce the critical lengthscale for instability. Numerical simulations showing the evolution of model clouds during the linear regime of evolution are also presented, and compared to the results of the dispersion analysis. They are found to be in agreement with the dispersion results, and confirm the dependence of the characteristic length and time scales on parameters such as the initial mass-to-flux ratio and relative external pressure.Comment: 30 pages, 7 figures Accepted by Ap

Three-dimensional Continuum Radiative Transfer Images of a Molecular Cloud Core Evolution

We analyze a three-dimensional smoothed particle hydrodynamics simulation of an evolving and later collapsing pre-stellar core. Using a three-dimensional continuum radiative transfer program, we generate images at 7 micron, 15 micron, 175 micron, and 1.3 mm for different evolutionary times and viewing angles. We discuss the observability of the properties of pre-stellar cores for the different wavelengths. For examples of non-symmetric fragments, it is shown that, misleadingly, the density profiles derived from a one-dimensional analysis of the corresponding images are consistent with one-dimensional core evolution models. We conclude that one-dimensional modeling based on column density interpretation of images does not produce reliable structural information and that multidimensional modeling is required.Comment: accepted by ApJL, 4 pages, 4 figure

A New Recipe for Obtaining Central Volume Densities of Prestellar Cores from Size Measurements

We propose a simple analytical method for estimating the central volume density of prestellar molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of the column density and volume density profiles of the same size indicating the existence of a uniform density inner region. The size of this region is set by the thermal pressure force which depends only on the central volume density and temperature of the core, and can provide a direct measurement of the central volume density. Thus a simple length measurement can immediately yield a central density estimate independent of any dynamical model for the core and without the need for fitting. Using the radius at which the column density is 90% of the central value as an estimate of the size of the flat inner part of the column density profile yields an estimate of the central volume density within a factor of 2 for well resolved cores.Comment: 4 pages, 4 figures, accepted for publication in ApJ Letter

Detection of 6 K gas in Ophiuchus D

Cold cores in interstellar molecular clouds represent the very first phase in star formation. The physical conditions of these objects are studied in order to understand how molecular clouds evolve and how stellar masses are determined. The purpose of this study is to probe conditions in the dense, starless clump Ophichus D (Oph D). The ground-state (1(10)-1(11)) rotational transition of ortho-H2D+ was observed with APEX towards the density peak of Oph D. The width of the H2D+ line indicates that the kinetic temperature in the core is about 6 K. So far, this is the most direct evidence of such cold gas in molecular clouds. The observed H2D+ spectrum can be reproduced with a hydrostatic model with the temperature increasing from about 6 K in the centre to almost 10 K at the surface. The model is unstable against any increase in the external pressure, and the core is likely to form a low-mass star. The results suggest that an equilibrium configuration is a feasible intermediate stage of star formation even if the larger scale structure of the cloud is thought to be determined by turbulent fragmentation. In comparison with the isothermal case, the inward decrease in the temperature makes smaller, i.e. less massive, cores susceptible to externally triggered collapse.Comment: 7 pages, 5 figures, accepted for Astronomy and Astrophysic

Lifetime of the embedded phase of low-mass star formation and the envelope depletion rates

Motivated by a considerable scatter in the observationally inferred lifetimes of the embedded phase of star formation, we study the duration of the Class 0 and Class I phases in upper-mass brown dwarfs and low-mass stars using numerical hydrodynamics simulations of the gravitational collapse of a large sample of cloud cores. We resolve the formation of a star/disk/envelope system and extend our numerical simulations to the late accretion phase when the envelope is nearly totally depleted of matter. We adopted a classification scheme of Andre et al. and calculate the lifetimes of the Class 0 and Class I phases (\tau_C0 and \tau_CI, respectively) based on the mass remaining in the envelope. When cloud cores with various rotation rates, masses, and sizes (but identical otherwise) are considered, our modeling reveals a sub-linear correlation between the Class 0 lifetimes and stellar masses in the Class 0 phase with the least-squares fit exponent m=0.8 \pm 0.05. The corresponding correlation between the Class I lifetimes and stellar masses in the Class I is super-linear with m=1.2 \pm 0.05. If a wider sample of cloud cores is considered, which includes possible variations in the initial gas temperature, cloud core truncation radii, density enhancement amplitudes, initial gas density and angular velocity profiles, and magnetic fields, then the corresponding exponents may decrease by as much as 0.3. The duration of the Class I phase is found to be longer than that of the Class~0 phase in most models, with a mean ratio \tau_CI / \tau_C0 \approx 1.5--2. A notable exception are YSOs that form from cloud cores with large initial density enhancements, in which case \tau_C0 may be greater than \tau_CI. Moreover, the upper-mass (>= 1.0 Msun) cloud cores with frozen-in magnetic fields and high cloud core rotation rates may have the \tau_CI / \tau_C0 ratios as large as 3.0--4.0. (Abdridged).Comment: Accepted for publication by The Astrophysical Journa

Grain size limits derived from 3.6 {\mu}m and 4.5 {\mu}m coreshine

Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 {\mu}m (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 {\mu}m. We aim to analyze 3.6 and 4.5 {\mu}m coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 {\mu}m bands for a sample of 10 low-mass cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes around 0.9, 0.5, 0.65, 1.5, 0.6, and > 1.5 {\mu}m, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. Considering surface brightness ratios in the 3.6 and 4.5 {\mu}m bands across a molecular cloud core is an effective method of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth.Comment: 23 pages, 18 figures, accepted for publication in A&