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

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

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

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

Ray-tracing for complex astrophysical high-opacity structures

We present a ray-tracing technique for radiative transfer modeling of complex three-dimensional (3D) structures which include dense regions of high optical depth like in dense molecular clouds, circumstellar disks, envelopes of evolved stars, and dust tori around active galactic nuclei. The corresponding continuum radiative transfer problem is described and the numerical requirements for inverse 3D density and temperature modeling are defined. We introduce a relative intensity and transform the radiative transfer equation along the rays to solve machine precision problems and to relax strong gradients in the source term. For the optically thick regions where common ray-tracers are forced to perform small trace steps, we give two criteria for making use of a simple approximative solver crossing the optically thick region quickly. Using an example of a density structure with optical depth changes of 6 orders of magnitude and sharp temperature variations, we demonstrate the accuracy of the proposed scheme using a common 5th-order Runge-Kutta ray-tracer with adaptive step size control. In our test case, the gain in computational speed is about a factor of 870. The method is applied to calculate the temperature distribution within a massive molecular cloud core for different boundary conditions for the radiation field.Comment: 21 pages, 5 figures to appear in Astrophysical Journa

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

Depletion and low gas temperature in the L183 prestellar core: the N2H+ - N2D+ tool

Context. The study of pre-stellar cores (PSCs) suffers from a lack of undepleted species to trace the gas physical properties in their very dense inner parts. Aims. We want to carry out detailed modelling of N2H+ and N2D+ cuts across the L183 main core to evaluate the depletion of these species and their usefulness as a probe of physical conditions in PSCs. Methods. We have developed a non-LTE (NLTE) Monte-Carlo code treating the 1D radiative transfer of both N2H+ and N2D+, making use of recently published collisional coefficients with He between individual hyperfine levels. The code includes line overlap between hyperfine transitions. An extensive set of core models is calculated and compared with observations. Special attention is paid to the issue of source coupling to the antenna beam. Results. The best fitting models indicate that i) gas in the core center is very cold (7$\pm$ 1 K) and thermalized with dust, ii) depletion of N2H+ does occur, starting at densities 5-7E5 cm−3 and reaching a factor of 6 (+13/−3) in abundance, iii) deuterium fractionation reaches ∼70% at the core center, and iv) the density profile is proportional to r^-1 out to ∼4000 AU, and to r^−2 beyond. Conclusions. Our NLTE code could be used to (re-)interpret recent and upcoming observations of N2H+ and N2D+ in many pre-stellar cores of interest, to obtain better temperature and abundance profiles

The Intrinsic Shapes of Molecular Cloud Fragments over a Range of Length Scales

We decipher intrinsic three-dimensional shape distributions of molecular clouds, cloud cores, Bok globules, and condensations using recently compiled catalogues of observed axis ratios for these objects mapped in carbon monoxide, ammonia, through optical selection, or in continuum dust emission. We apply statistical techniques to compare assumed intrinsic axis ratio distributions with observed projected axis ratio distributions. Intrinsically triaxial shapes produce projected distributions which agree with observations. Molecular clouds mapped in $^{12}$CO are intrinsically triaxial but more nearly prolate than oblate, while the smaller cloud cores, Bok globules, and condensations are also intrinsically triaxial but more nearly oblate than prolate.Comment: 12 pages, 11 figures. Version with color figures can be found at http://www.astro.uwo.ca/~cjones/ or http://www.astro.uwo.ca/~basu/. To appear in ApJ, 10 April 2002, v. 569, no.

Upper limit for the D2H+ ortho-to-para ratio in the prestellar core 16293E (CHESS)

The H3+ ion plays a key role in the chemistry of dense interstellar gas clouds where stars and planets are forming. The low temperatures and high extinctions of such clouds make direct observations of H3+ impossible, but lead to large abundances of H2D+ and D2H+, which are very useful probes of the early stages of star and planet formation. The ground-state rotational ortho-D2H+ 111-000 transition at 1476.6 GHz in the prestellar core 16293E has been searched for with the Herschel/HIFI instrument, within the CHESS (Chemical HErschel Surveys of Star forming regions) Key Program. The line has not been detected at the 21 mK km/s level (3 sigma integrated line intensity). We used the ortho-H2D+ 110-111 transition and para-D2H+ 110-101 transition detected in this source to determine an upper limit on the ortho-to-para D2H+ ratio as well as the para-D2H+/ortho-H2D+ ratio from a non-LTE analysis. The comparison between our chemical modeling and the observations suggests that the CO depletion must be high (larger than 100), with a density between 5e5 and 1e6 cm-3. Also the upper limit on the ortho-D2H+ line is consistent with a low gas temperature (~ 11 K) with a ortho-to-para ratio of 6 to 9, i.e. 2 to 3 times higher than the value estimated from the chemical modeling, making it impossible to detect this high frequency transition with the present state of the art receivers.Comment: Accepted in A&