Search

Deuterated ions are abundant in cold (T=10 K), dense (n=10^5 cm^-3) regions, in which CO is frozen out onto dust grains. In such environments, the deuterium fractionation of such ions can exceed the elemental abundance ratio of D/H by a factor of 10^4. In this paper we use the deuterium fractionation to investigate the evolutionary state of Class 0 protostars. In a sample of 20 protostellar objects, we found a clear correlation between the N2D+/N2H+ ratio and evolutionary tracers. As expected, the coolest, i.e. the youngest, objects show the largest deuterium fractionation. Furthermore, we find that sources with a high N2D+/N2H+ ratio show clear indication for infall., Comment: 19 pages, 12 figures, accepted by A&A

We present a survey of the ortho-H2D+(1_{1,0}-1_{1,1}) line toward a sample of 10 starless cores and 6 protostellar cores, carried out at the Caltech Submillimeter Observatory. The high diagnostic power of this line is revealed for the study of the chemistry, and the evolutionary and dynamical status of low-mass dense cores. The line is detected in 7 starless cores and in 4 protostellar cores. N(ortho-H2D+) ranges between 2 and 40x10^{12} cm^{-2} in starless cores and between 2 and 9x10^{12} cm^{-2} in protostellar cores. The brightest lines are detected toward the densest and most centrally concentrated starless cores, where the CO depletion factor and the deuterium fractionation are also largest. The large scatter observed in plots of N(ortho-H2D+) vs. the observed deuterium fractionation and vs. the CO depletion factor is likely to be due to variations in the ortho-to-para (o/p) ratio of H2D+ from >0.5 for T_{kin} < 10 K gas in pre-stellar cores to ~0.03 (consistent with T_{kin} ~15 K for protostellar cores). The two Ophiuchus cores in our sample also require a relatively low o/p ratio (~0.3). Other parameters, including the cosmic-ray ionization rate, the CO depletion factor (or, more in general, the depletion factor of neutral species), the volume density, the fraction of dust grains and PAHs also largely affect the ortho-H2D+ abundance. The most deuterated and H2D+-rich objects (L429, L1544, L694-2 and L183) are reproduced by chemical models of centrally concentrated (central densties ~10^{6} cm^{-3}) cores with chemical ages between 10^4 and 10^6 yr. Upper limits of the para-H3O+ (1_1- -2_1+) and para-D2H+ (1_{1,0}-1_{0,1}) lines are also given. (Abridged), Comment: 29 pages, 6 figures. To appear in A&A

Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to emerge in the 1990s. However, detailed models of these shocks are required for the calculation of the grain sputtering contribution to gas phase abundances of species producing observed emissions. We are developing such models. Some of the molecular species introduced into the gas phase by sputtering in shocks or by thermally driven desorption in hot cores form on grain surfaces. Recently laboratory studies have begun to contribute to the understanding of surface reactions and thermally driven desorption important for the chemistry of star forming clouds. Dusty plasmas are prevalent in many evolved stars just as well as in star forming regions. Radiation pressure on dust plays a significant role in mass loss from some post-main-sequence stars. The mechanisms leading to the formation of carbonaceous dust in the stellar outflows are similar to those important for soot formation in flames. However, nucleation in oxygen-rich outflows is less well understood and remains a challenging research area. Dust is observed in supernova ejecta that have not passed through the reverse shocks that develop in the interaction of ejecta with ambient media. Dust is detected in high redshift galaxies that are sufficiently young that the only stars that could have produced the dust were so massive that they became supernovae. Consequently, the issue of the survival of dust in strong supernova shocks is of considerable interest., Comment: 4 pages, to be published in the proceedings of Fifth International Conference on Physics of Dusty Plasmas

We describe a model for the thermal and dynamical equilibrium of starless cores that includes the radiative transfer of the gas and dust and simple CO chemistry. The model shows that the structure and behavior of the cores is significantly different depending on whether the central density is either above or below about 10^5 cm-3. This density is significant as the critical density for gas cooling by gas-dust collisions and also as the critical density for dynamical stability, given the typical properties of the starless cores. The starless cores thus divide into two classes that we refer to as thermally super-critical and thermally sub-critical.This two-class distinction allows an improved interpretation of the different observational data of starless cores within a single model., Comment: ApJ in press

Aims: The production of molecular hydrogen and its deuterated forms onto carbonaceous dust grains is investigated in detail. The goal of this study is to estimate the importance of the chemistry occuring on grain surfaces for the deuteration of H2. Furthermore, we aim to find a robust and general surface chemical model which can be used in different astrophysical environments. Methods: Surface processes are described for the cases of graphitic and amorphous--carbon grains, where laboratory work is available. Langmuir--Hinshelwood as well as Eley--Rideal surface chemistries are included in the model and their relative contributions are highlighted. Analytic expressions are derived for H2, HD, and D2 formation efficiencies for both type of grains. Rate equations are tested against stochastic methods. Results: As expected, rate equations and stochastic methods diverge for grain sizes lower than a critical value acrit. For grain sizes below this critical value, D2 formation decreases to favour HD formation. The formation efficiencies of H2 and D2 can be calculated by adding a correction factor to the rate equations methods. We found that because of the presence of chemisorbed sites, which can store atoms to form molecules up to high grain temperatures, the formation efficiency of HD and D2 is very high compared to models where only physisorption sites are taken into account. When considering a realistic distribution of dust grains, we found that the formation rate of H2 and HD is enhanced by an order of magnitude if small grains are taken into account. The processes described in this paper, that allow a strong enhancement of the deuterated forms of molecular hydrogen, could explain the high degree of deuterium fractionation observed in protostellar environments., Comment: 17 pages 25 figures

We use the COMPLETE data to derive new calibrations of the X-factor and the 13CO abundance within Perseus. We divide Perseus into six sub-regions. The standard X factor, X=N(H2)/W(12CO), is derived both for the whole Perseus Complex and for each of the six sub-regions with values consistent with previous estimates. The X factor is heavily affected by the saturation of the emission above AV~4 mag, and variations are found between regions. We derive linear fits to relate W(12CO) and AV using only points below 4 mag of extinction, this yields a better estimation of the AV than the X-factor. We derive linear relations of W(13CO), N(13CO) and W(C18O) with AV . The extinction threshold above which 13CO(1-0) and C18O(1-0) are detected is about 1 mag larger than previous estimates. 12CO(1-0) and 13CO(1-0) lines saturate above 4 and 5 mag, respectively, whereas C18O(1-0) never saturates (up to 10 mag). Approximately 60% of the positions with 12CO emission have sub-thermally excited lines, and almost all positions have 12CO excitation temperatures below the dust temperature. Using the Meudon PDR code we find that 12CO and 13CO emission can be explained by uniform slab models with densities ranging between about 10^3 and 10^4 cm-3. Local variations in the volume density and non-thermal motions (linked to different star formation activity) can explain the observations. Higher densities are needed to reproduce CO data toward active star forming sites, where the larger internal motions driven by the young protostars allow more photons from the embedded high density cores to escape the cloud. In the most quiescent region, the 12CO and 13CO emission appears to arise from an almost uniform thin layer of molecular material at densities around 10^4 cm-3., Comment: 40 pages, 12 figures, accepted for publication in ApJ; version with high resolution figures available at http://www.cfa.harvard.edu/~jpineda/post/cal-co-v2.pdf

Context: The detection of narrow SiO lines toward the young shocks of the L1448-mm outflow has been interpreted as a signature of the magnetic precursor of C-shocks. In contrast with the low SiO abundances (<10E-12) in the ambient gas, the narrow SiO emission at almost ambient velocities reveals enhanced SiO abundances of 10E-11. This enhancement has been proposed to be produced by the sputtering of the grain mantles at the first stages of C-shocks. However, modelling of the sputtering of grains has usually averaged the SiO abundances over the dissipation region of C-shocks, which cannot explain the recent observations. Aims: To model the evolution of the gas phase abundances of SiO, CH3OH and H2O, produced by the sputtering of grains as the shock propagates through the ambient gas. Methods: We propose a parametric model to describe the physical structure of C-shocks as a function of time. Using the known sputtering yields for water mantles (with minor constituents like silicon and CH3OH) and olivine cores by collisions with H2, He, C, O, Si, Fe and CO, we follow the evolution of the abundances of silicon, CH3OH and H2O ejected from grains. Results: The evolution of these abundances shows that CO seems to be the most efficient sputtering agent in low velocity shocks. The velocity threshold for the sputtering of silicon from the grain mantles is reduced by 5-10 km s-1 by CO compared to other models. The sputtering by CO can generate SiO abundances of 10E-11 at the early stages of low velocity shocks, consistent with those observed in the magnetic precursor of L1448-mm. Our model also satisfactorily reproduce the progressive enhancement of SiO, CH3OH and H2O observed in this outflow by the coexistence of two shocks with vs=30 and 60kms-1 within the same region., Comment: 12 pages, 7 figures, accepted for publication in A&A

We have observed the deuterated gas in the high-mass star formation region IRAS 05345+3157 at high-angular resolution, in order to determine the morphology and the nature of such gas. We have mapped the N2H+ (1-0) line with the Plateau de Bure Interferometer, and the N2D+ (3-2) and N2H+ (3-2) lines with the Submillimeter Array. The N2D+ (3-2) integrated emission is concentrated in two condensations, with masses of 2-3 and 9 M_sun and diameters of 0.05 and 0.09 pc, respectively. The high deuterium fractionation (0.1) and the line parameters in the N2D+ condensations indicate that they are likely low- to intermediate-mass pre-stellar cores, even though other scenarios are possible., Comment: 4 pages, 2 figures, accepted for publication in Astronomy and Astrophysics

We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the (sub)millimeter sources as well as at several positions with no associated continuum emission. For each detection, we have measured physical properties by fitting a simple model to every spectral line simultaneously. Where appropriate, we have refined the model by accounting for low optical depths, multiple components along the line of sight and imperfect coupling to the GBT beam. For the cores in Perseus, we find a typical kinetic temperature of T=11 K, a typical column density of N(NH3)~ 10^14.5 /cm^2 and velocity dispersions ranging from sigma_v = 0.07 km/s to 0.7 km/s. However, many cores with velocity dispersions > 0.2 km/s show evidence for multiple velocity components along the line of sight., Comment: 19 pages; Accepted to ApJS; version with high resolution figures available at http://www.cfa.harvard.edu/COMPLETE/papers/nh3-paper1.pdf ; online data at http://www.cfa.harvard.edu/COMPLETE/data_html_pages/GBT_NH3.html

We present observations with the IRAM Plateau de Bure Interferometer of three QSOs at z>5 aimed at detecting molecular gas in their host galaxies as traced by CO transitions. CO (5-4) is detected in SDSSJ033829.31+002156.3 at z=5.0267, placing it amongst the most distant sources detected in CO. The CO emission is unresolved with a beam size of ~1", implying that the molecular gas is contained within a compact region, less than ~3kpc in radius. We infer an upper limit on the dynamical mass of the CO emitting region of ~3x10^10 Msun/sin(i)^2. The comparison with the Black Hole mass inferred from near-IR data suggests that the BH-to-bulge mass ratio in this galaxy is significantly higher than in local galaxies. From the CO luminosity we infer a mass reservoir of molecular gas as high as M(H2)=2.4x10^10 Msun, implying that the molecular gas accounts for a significant fraction of the dynamical mass. When compared to the star formation rate derived from the far-IR luminosity, we infer a very short gas exhaustion timescale (~10^7 yrs), comparable to the dynamical timescale. CO is not detected in the other two QSOs (SDSSJ083643.85+005453.3 and SDSSJ163033.90+401209.6) and upper limits are given for their molecular gas content. When combined with CO observations of other type 1 AGNs, spanning a wide redshift range (05) the CO luminosity appears to saturate. We discuss the implications of these findings in terms of black hole-galaxy co-evolution., Comment: Accepted for publication in A&A Letters, 6 pages, 3 figures

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr., Comment: 29 pages, including 5 tables and 15 figures

We investigate the mm-morphology of IC1396N with unprecedented spatial resolution to analyze its dust and molecular gas properties, and draw comparisons with objects of similar mass. We have carried out sensitive observations in the most extended configurations of the IRAM Plateau de Bure interferometer, to map the thermal dust emission at 3.3 and 1.3mm, and the emission from the $J$=13$_k\to12_k$ hyperfine transitions of methyl cyanide (CH$_3$CN). We unveil the existence of a sub-cluster of hot cores in IC1396N, distributed in a direction perpendicular to the emanating outflow. The cores are embedded in a common envelope of extended and diffuse dust emission. We find striking differences in the dust properties of the cores ($\beta\simeq$ 0) and the surrounding envelope ($\beta\simeq$ 1), very likely testifying to differences in the formation and processing of dust material. The CH$_3$CN emission peaks towards the most massive hot core and is marginally extended in the outflow direction.

We have observed several emission lines of two Nitrogen-bearing (C2H5CN and C2H3CN) and two Oxygen-bearing (CH3OCH3 and HCOOCH3) molecules towards a sample of well-known hot molecular cores (HMCs) in order to check whether the chemical differentiation seen in the Orion-HMC and W3(H_2O) between O- and N-bearing molecules is a general property of HMCs. With the IRAM-30m telescope we have observed 12 HMCs in 21 bands, centered at frequencies from 86250 to 258280 MHz. The rotational temperatures obtained range from ~100 to ~150 K in these HMCs. Single Gaussian fits performed to unblended lines show a marginal difference in the line peak velocities of the C2H5CN and CH3OCH3 lines, indicating a possible spatial separation between the region traced by the two molecules. On the other hand, neither the linewidths nor the rotational temperatures and column densities confirm such a result. By comparing the abundance ratio of the pair C2H5CN/C2H3CN with the predictions of theoretical models, we derive that the age of our cores ranges between 3.7 and 5.9x10^{4} yrs. The abundances of C2H5CN and C2H3CN are strongly correlated, as expected from theory which predicts that C2H3CN is formed through gas phase reactions involving C2H5CN. A correlation is also found between the abundances of C2H3CN and CH3OCH3, and C2H5CN and CH3OCH3. In all tracers the fractional abundances increase with the H_2 column density while they are not correlated with the gas temperature., Comment: Accepted for publication in A&A, 56 pages

Abridged: The thermal structure of a starless core is crucial for our understanding of the physics in these objects and hence for our understanding of star formation. Theory predicts a gas temperature drop in the inner 5000 AU of these objects, but there has been little observational proof of this. We performed VLA observations of the NH3 (1,1) and (2,2) transitions towards the pre-stellar core L 1544 in order to measure the temperature gradient between the high density core nucleus and the surrounding core envelope. Our VLA observation for the first time provide measurements of gas temperature in a core with a resolution smaller than 1000 AU. We have also obtained high resolution Plateau de Bure observations of the 110 GHz 111-101 para-NH2D line in order to further constrain the physical parameters of the high density nucleus. We have estimated the temperature gradient using a model of the source to fit our data in the u,v plane. We find that indeed the temperature decreases toward the core nucleus from 12 K down to 5.5 K resulting in an increase of a factor of 50% in the estimated density of the core from the dust continuum if compared with the estimates done with constant temperature of 8.75 K. We also found a remarkably high abundance of deuterated ammonia with respect to the ammonia abundance (50%+-20%), which proves the persistence of nitrogen bearing molecules at very high densities (2e6 cm-3) and shows that high-resolution observations yield higher deuteration values than single-dish observations. Our analysis of the NH3 and NH2D kinematic fields shows a decrease of specific angular momentum from the large scales to the small scales., Comment: 12 pages, 6 figures. Accepted for publication by A&A

The transition between the low density groups of T Tauri stars and the high density clusters around massive stars occurs in the intermediate-mass (IM) range (M$_*$$\sim$2--8 M$_\odot$). High spatial resolution studies of IM young stellar objects (YSO) can provide important clues to understand the clustering in massive star forming regions. Aims: Our aim is to search for clustering in IM Class 0 protostars. The high spatial resolution and sensitivity provided by the new A configuration of the Plateau de Bure Interferometer (PdBI) allow us to study the clustering in these nearby objects. Methods: We have imaged three IM Class 0 protostars (Serpens-FIRS 1, IC 1396 N, CB 3) in the continuum at 3.3 and 1.3mm using the PdBI. The sources have been selected with different luminosity to investigate the dependence of the clustering process on the luminosity of the source. Results: Only one millimeter (mm) source is detected towards the low luminosity source Serpens--FIRS 1. Towards CB 3 and IC1396 N, we detect two compact sources separated by $\sim$0.05 pc. The 1.3mm image of IC 1396 N, which provides the highest spatial resolution, reveal that one of these cores is splitted in, at least, three individual sources., Comment: 4 pages, 3 figures, accepted for publication in Astronomy and Astrophysics Letters (Special Feature IRAM/PdB)

We have measured the deuterium fractionation (through the column density ratio N(N2D+)/N(N2H+)) and the CO depletion factor (ratio between expected and observed CO abundance) in a sample of 10 high-mass protostellar candidates, in order to understand whether the earliest evolutionary stages of high-mass stars have chemical characteristics similar to those of low-mass ones. The observations were carried out with the IRAM-30m telescope and the JCMT. We have detected N2D+ emission in 7 of the 10 sources of our sample, and found an average value N(N2D+)/N(N2H+)~0.015. This value is 3 orders of magnitude larger than the interstellar D/H ratio, indicating the presence of cold and dense gas, in which the physical-chemical conditions are similar to those observed in low-mass pre-stellar cores. Also, the integrated CO depletion factors show that in the majority of the sources the expected CO abundances are larger than the observed values, with a median ratio of 3.2. In principle, the cold gas that gives origin to the N2D+ emission can be the remnant of the massive molecular core in which the high-mass (proto-)star was born, not yet heated up by the central object. If so, our results indicate that the chemical properties of the clouds in which high-mass stars are born are similar to their low-mass counterparts. Alternatively, this cold gas can be located into one (or more) starless core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot decide which is the correct scenario., Comment: 27 pages, 15 figures, accepted by Astronomy and Astrophysics

Aims. Our goals are to compare the CS, N2H+ and dust distributions in a representative sample of high-mass star forming dense cores and to determine the physical and chemical properties of these cores. Methods. We compare the results of CS(5-4) and 1.2 mm continuum mapping of twelve dense cores from the southern hemisphere presented in this work, in combination with our previous N2H+(1-0) and CS(2-1) data. We use numerical modeling of molecular excitation to estimate physical parameters of the cores. Results. Most of the maps have several emission peaks (clumps). We derive basic physical parameters of the clumps and estimate CS and N2H+ abundances. Masses calculated from LVG densities are higher than CS virial masses and masses derived from continuum data, implying small-scale clumpiness of the cores. For most of the objects, the CS and continuum peaks are close to the IRAS point source positions. The CS(5-4) intensities correlate with continuum fluxes per beam in all cases, but only in five cases with the N2H+(1-0) intensities. The study of spatial variations of molecular integrated intensity ratios to continuum fluxes reveals that I(N2H+)/F{1.2} ratios drop towards the CS peaks for most of the sources, which can be due to a N2H+ abundance decrease. For CS(5-4), the I(CS)/F{1.2} ratios show no clear trends with distance from the CS peaks, while for CS(2-1) such ratios drop towards these peaks. Possible explanations of these results are considered. The analysis of normalized velocity differences between CS and N2H+ lines has not revealed indications of systematic motions towards CS peaks., Comment: 13 pages, 5 figures, accepted by Astronomy and Astrophysics

The presence of H2D+ in dense cloud cores underlies ion-molecule reactions that strongly enhance the deuterium fractionation of many molecular species. We determine the H2D+ abundance in one starless core, Barnard 68, that has a particularly well established physical, chemical, and dynamical structure. We observed the ortho-H2D+ ground-state line 1_10-1_11, the N2H+ J=4-3 line, and the H13CO+ 4-3 line with the APEX telescope. We report the probable detection of the o-H2D+ line at an intensity Tmb=0.22 +- 0.08 K and exclusively thermal line width, and find only upper limits to the N2H+ 4-3 and H13CO+ 4-3 intensities. Within the uncertainties in the chemical reaction rates and the collisional excitation rates, chemical model calculations and excitation simulations reproduce the observed intensities and that of o-H2D+ in particular., Comment: APEX A&A special issue, accepted

[Abridged] We present a molecular survey of the starless cores L1498 and L1517B. These cores have been selected for their relative isolation and close-to-round shape, and they have been observed in a number of lines of 13 molecular species (4 already presented in the first part of this series): CO, CS, N2H+, NH3, CH3OH, SO, C3H2, HC3N, C2S, HCN, H2CO, HCO+, and DCO+. Using a physical model of core structure and a Monte Carlo radiative transfer code, we determine for each core a self-consistent set abundances that fits simultaneously the observed radial profile of integrated intensity and the emergent spectrum towards the core center (for abundant species, optically thin isopologues are used). From this work, we find that L1498 and L1517B have similar abundance patterns, with most species suffering a significant drop toward the core center. This occurs for CO, CS, CH3OH, SO, C3H2, HC3N, C2S, HCN, H2CO, HCO+, and DCO+, which we fit with profiles having a sharp central hole. The size of this hole varies with molecule: DCO+, HCN, and HC3N have the smallest holes while SO, C2S and CO have the largest holes. Only N2H+ and NH3 are present in the gas phase at the core centers. From the different behavior of molecules, we select SO, C2S, and CH3OH as the most sensitive tracers of molecular depletion. Comparing our abundance determinations with the predictions from current chemical models we find order of magnitude discrepancies. Finally, we show how the ``contribution function'' can be used to study the formation of line profiles from the different regions of a core., Comment: 22 pages, 12 figures, A&A accepted

Deuterium enhancement of monodeuterated species has been recognized for more than 30 years as a result of the chemical fractionation that results from the difference in zero point energies of deuterated and hydrogenated molecules. The key reaction is the deuteron exchange in the reaction between HD, the reservoir of deuterium in dark interstellar clouds, and the H3+ molecular ion, leading to the production of the H2D+ molecule, and the low temperature in dark interstellar clouds favors this production. Furthermore, the presence of multiply deuterated species have incited our group to proceed further and consider the subsequent reaction of H2D+ with HD, leading to D2H+ (first detected by Vastel et al. 2004), which can further react with HD to produce D3+. In prestellar cores, where CO was found to be depleted (Bacmann et al. 2003), this production should be increased, as CO would normally destroy H3+. The first model including D2H+ and D3+ (Roberts, Herbst & Millar 2003) predicted that these molecules should be as abundant as H2D+ (see contribution by H. Roberts). The first detection of the D2H+ was made possible by the recent laboratory measurement by Hirao & Amano (2003) for the frequency of the fundamental line of the para-D2H+ (see contribution by T. Amano). Here we present observations of H2D+ and D2H+ towards a sample of dark clouds and prestellar cores and show how the distribution of ortho-H2D+ (1_1,0-1_1,1) can trace the deuterium factory in prestellar cores. We also present how future instrumentation will improve our knowledge concerning the deuterium enhancement of H3+., Comment: Physics, chemistry and astronomy of H3+, accepted in Phil. Trans. R. Soc. Lond. A (2006)

Prestellar cores are unique laboratories for studies of the chemical and physical conditions preceding star formation. We observed the prestellar core L1544 in the fundamental transition of ortho-H2D+ (1_1,0-1_1,1) at different positions over 100", and found a strong correlation between its abundance and the CO depletion factor. We also present a tentative detection of the fundamental transition of para-D2H+ (1_1,0-1_0,1) at the dust emission peak. Maps in N2H+, N2D+, HCO+ and DCO+ are used, and interpreted with the aid of a spherically symmetric chemical model that predicts the column densities and abundances of these species as a function of radius. The correlation between the observed deuterium fractionation of H3+, N2H+ and HCO+ and the observed integrated CO depletion factor across the core can be reproduced by this chemical model. In addition a simpler model is used to study the H2D+ ortho-to-para ratio. We conclude that, in order to reproduce the observed ortho-H2D+ observations, the grain radius should be larger than 0.3 microns., Comment: 24 pages, 9 figures, accepted in ApJ (to be published in July 2006)

We summarize some recent results on the metallicity and dust properties of Active Galactic Nuclei (AGN) at high redshift (110). The properties of dust in high-z QSOs are discussed within the context of the dust production mechanisms in the early universe. The dust extinction curve is observed to evolve beyond z>4, and by z~6 it is well described by the properties expected for dust produced by SNe, suggesting that the latter is the main mechanism of dust production in the early universe. We also show that the huge dust masses observed in distant QSOs can be accounted for by SN dust within the observational constraints currently available. Finally, we show that QSO winds, which have been proposed as an alternative mechanism of dust production, may also contribute significantly to the total dust budget at high redshift., Comment: 10 pages, 8 figures, invited talk at the Workshop "AGN and galaxy evolution", Specola Vaticana, Castel Gandolfo, Italy, 3-6 October 2005

Low-mass protostars form from condensations inside molecular clouds when gravity overwhelms thermal and magnetic supporting forces. The first phases of the formation of a solar-type star are characterized by dramatic changes not only in the physical structure but also in the chemical composition. Since PPIV (e.g., Langer et al.), exciting new developments have occurred in our understanding of the processes driving this chemical evolution. These developments include two new discoveries : 1) extremely enhanced molecular deuteration, which is caused by the freeze-out of heavy-element-bearing molecules onto grain mantles during the Prestellar Core and Class 0 source phases; and 2) hot corinos, which are warm and dense regions at the center of Class 0 source envelopes and which are characterized by a multitude of complex organic molecules. In this chapter we will review these two new topics, and will show how they contribute to our understanding of the first phases of solar-type stars., Comment: To appear in Protostars and Planets V

We present an overview of data available for the Ophiuchus and Perseus molecular clouds from ``Phase I'' of the COMPLETE Survey of Star-Forming Regions. This survey provides a range of data complementary to the Spitzer Legacy Program ``From Molecular Cores to Planet Forming Disks.'' Phase I includes: Extinction maps derived from 2MASS near-infrared data using the NICER algorithm; extinction and temperature maps derived from IRAS 60 and 100um emission; HI maps of atomic gas; 12CO and 13CO maps of molecular gas; and submillimetre continuum images of emission from dust in dense cores. Not unexpectedly, the morphology of the regions appears quite different depending on the column-density tracer which is used, with IRAS tracing mainly warmer dust and CO being biased by chemical, excitation and optical depth effects. Histograms of column-density distribution are presented, showing that extinction as derived from 2MASS/NICER gives the closest match to a log-normal distribution as is predicted by numerical simulations. All the data presented in this paper, and links to more detailed publications on their implications are publically available at the COMPLETE website., Comment: Accepted by AJ. Full resolution version available from: http://www.cfa.harvard.edu/COMPLETE/papers/complete_phase1.pdf

Low-mass dense cores represent the state of molecular gas associated with the earliest phases of low-mass star formation. Such cores are called "protostellar" or "starless," depending on whether they do or do not contain compact sources of luminosity. In this chapter, the first half of the review of low-mass dense cores, we describe the numerous inferences made about the nature of starless cores as a result of recent observations, since these reveal the initial conditions of star formation. We focus on the identification of isolated starless cores and their internal physical and chemical properties, including morphologies, densities, temperatures, kinematics, and molecular abundances. These objects display a wide range of properties since they are each at different points on evolutionary paths from ambient molecular cloud material to cold, contracting, and centrally concentrated configurations with significant molecular depletions and, in rare cases, enhancements., Comment: 16 pages, to appear in "Protostars and Planets V"

We found that in regions of high mass star formation the CS emission correlates well with the dust continuum emission and is therefore a good tracer of the total mass while the N$_2$H$^+$ distribution is frequently very different. This is opposite to their typical behavior in low-mass cores where freeze-out plays a crucial role in the chemistry. The behavior of other high density tracers varies from source to source but most of them are closer to CS. Radial density profiles in massive cores are fitted by power laws with indices about -1.6, as derived from the dust continuum emission. The radial temperature dependence on intermediate scales is close to the theoretically expected one for a centrally heated optically thin cloud. The velocity dispersion either remains constant or decreases from the core center to the edge. Several cores including those without known embedded IR sources show signs of infall motions. They can represent the earliest phases of massive protostars. There are implicit arguments in favor of small-scale clumpiness in the cores., Comment: 6 pages, 4 figures, to be published in "Massive star birth: A crossroads of Astrophysics", IAU Symposium Proceedings of the international Astronomical Union 227, Held 16-20 May, Italy, edited by Cesaroni, R.; Felli, M.; Churchwell, E.; Walmsley, M. Cambridge: Cambridge University Press, 2005

It is not common to consider the role of uncertainties in the rate coefficients used in interstellar gas-phase chemical models. In this paper, we report a new method to determine both the uncertainties in calculated molecular abundances and their sensitivities to underlying uncertainties in the kinetic data utilized. The method is used in hot core models to determine if previous analyses of the age and the applicable cosmic-ray ionization rate are valid. We conclude that for young hot cores ($\le 10^4$ yr), the modeling uncertainties related to rate coefficients are reasonable so that comparisons with observations make sense. On the contrary, the modeling of older hot cores is characterized by strong uncertainties for some of the important species. In both cases, it is crucial to take into account these uncertainties to draw conclusions from the comparison of observations with chemical models., Comment: Accepted for publication in A&A

We report the first detection of the ground transition of the deuterated water at 464 GHz in the young proto-planetary disk surrounding the solar type protostar DM Tau. The line is observed in absorption against the continuum from the cold dust in the disk midplane, with a line to continuum ratio close to unity. The observation implies that deuterated gaseous water is present, with a relatively large abundance ($\sim 3\times10^{-9}$), in the outer disk above the midplane, where the density is, within a factor ten, $\sim 10^6$ cm$^{-3}$ and the temperature is lower than about 25 K. In these conditions, the H$_2$O condensation timescale is much smaller than the DM Tau disk age, and, therefore, water should be fully frozen onto the grain mantles. We suggest that UV photons and/or X-rays sublimate part of the mantles re-injecting the ices into the gas phase. Even though there is currently no measurement of H$_2$O, we provide arguments that the HDO/H$_2$O ratio should be about 0.01 or larger, which would be hundreds of times larger than the values measured in Solar System objects. This suggests the need of strong caution in comparing and linking the HDO/H$_2$O in Solar System and star forming environments., Comment: Accepted to be published in ApJ Letters

We report the detection of the 2P_3/2 -> 2P_1/2 fine-structure line of C+ at 157.74 micron in SDSSJ114816.64+525150.3 (hereafter J1148+5251), the most distant known quasar, at z=6.42, using the IRAM 30-meter telescope. This is the first detection of the [CII] line at high redshift, and also the first detection in a Hyperluminous Infrared Galaxy (L_FIR > 10^13 Lsun). The [CII] line is detected at a significance level of 8 sigma and has a luminosity of 4.4 x 10^9 Lsun. The L_[CII]/L_FIR ratio is 2 x 10^-4, about an order of magnitude smaller than observed in local normal galaxies and similar to the ratio observed in local Ultraluminous Infrared Galaxies. The [CII] line luminosity indicates that the host galaxy of this quasar is undergoing an intense burst of star formation with an estimated rate of ~3000 Msun/yr. The detection of C+ in SDSS J1148+5251 suggests a significant enrichment of metals at z ~ 6 (age of the universe ~870 Myr), although the data are consistent with a reduced carbon to oxygen ratio as expected from chemical evolutionary models of the early phases of galaxy formation., Comment: 5 pages, 2 figures, accepted by A&A Letters

We report high angular resolution (~0.6''x0.5'' at 1.3mm) observations of the Class 0 intermediate-mass (IM) protostar NGC 7129--FIRS 2 using the PdBI. Our observations show the existence of a hot core at the position of the Class 0 object. This is, up to our knowledge, the first IM hot core detected so far. Interferometric maps of the region in the the continuum at 3mm and 1.3mm and in the CH3CN 5_k-->4_k, CH3OH 5_kk'-->4_kk', and D2CO 4_04-->3_03 lines are presented in this paper. Enhanced CH3CN and CH3OH abundances are measured towards the hot core (X(CH3CN)~7.010^{-9}, X(CH3OH)~3 10^{-8}-- a few 10^{-7}). While intense D2CO emission is detected towards the hot core, the N2D+ line has not been detected in our interferometric observations. The different behaviors of D2CO and N2D+ emissions suggest different formation mechanisms for the two species and different deuteration processes for H2CO and N2H+ (surface and gas-phase chemistry, respectively). Finally, the spectrum of the large bandwidth correlator show a forest of lines at the hot core position revealing that this object is extraordinarily rich in complex molecules., Comment: 14 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

We present the results of a multiline survey at mm-wavelengths of the Cepheus A star forming region. Four main flows have been identified: three pointing in the SW, NE, and SE directions and accelerating high density CS clumps. The fourth outflow, revealed by high-sensitivity HDO observations, is pointing towards South and is associated with conditions particularly favourable to a chemical enrichment. The analysis of the line profiles shows that the SiO molecule dominates at the highest velocities and at the highest excitation conditions, confirming its close association with shocks. H2S, SO2, and SO preferentially trace more quiescent regions than SiO, and in particular a lack of bright H2S emission at the highest velocities is found. OCS and H2CS emit at quite high velocities, where the abundances of three shock tracers like SiO, CH3OH, and HDO are higher. These results may indicate that H2S is not the only major sulphur carrier in the grain mantles, and that OCS and H2CS may probably play an important role on the grains; or that alternatively they rapidly form once the mantle is evaporated after the passage of a shock. Finally, the outflow peak emission has been compared with recent time-dependent sulphur chemistry models: the results indicate that, if associated with accurate measurements of the physical conditions, the CH3OH/H2CS column density ratio can be used as an effective chemical clock to date the age of shocked gas., Comment: Accepted by MNRAS

Velocity profiles of ground state lines of H2D+, HC18O+ and N2H+, observed previously with the CSO and IRAM 30m telescopes, are modeled with a Monte Carlo radiative transfer program to study the temperature, density and velocity structure of the pre-stellar core LDN 1544. The H2D+ line is double-peaked like that of the other ions, but previous models that fit the HC18O+ and N2H+ profiles are found not to fit the H2D+ data. Matching the H2D+ observations requires at least three modifications to the model at small radii: (1) the density profile must continue to rise inward and not flatten off toward the center; (2) the gas temperature must be nearly constant and not drop inwards significantly; (3) the infall velocity must increase inward, in a fashion intermediate between `quasi-static' (ambipolar diffusion) and `fully dynamic' (Larson-Penston) collapse. The C18O emission indicates a chemical age of <~0.1 Myr. The effects of a flattened structure and rotation on the line profiles are shown to be unimportant, at least on the scales probed by single-dish telescopes. Alternatively, the H2D+ profile is affected by absorption in the outer layers of the core, if gas motions in these layers are sufficiently small., Comment: Accepted by A&A; 10 pages, 10 colour figures

This paper will review the basic concepts of gas-phase and grain surface chemistry of dense molecular clouds, where low mass and high mass stars form. The chemistry of cold pre-stellar cloud cores, where molecular freeze-out and deuterium fractionation dominate, will be presented. Then, following cloud evolution after protostellar birth, hot core and shock chemistry will be discussed in view of recent observations. A brief summary of the chemistry in protoplanetary disks will also be furnished. The aim is to identify important gas tracers in the various steps of the star formation process, pointing out the main problems still open in the field of astrochemistry., Comment: Review article to appear in proceedings of the "Cores to Clusters" meeting held in Porto, Portugal, October 2004

We reconsider H2 formation on grain surfaces. We develop a rate equation model which takes into account the presence of both physisorbed and chemisorbed sites on the surface, including quantum mechanical tunnelling and thermal diffusion. In this study, we took into consideration the uncertainties on the characteristics of graphitic surfaces. We calculate the H2 formation efficiency with the Langmuir Hinshelwood and Eley Rideal mechanisms, and discuss the importance of these mechanisms for a wide range of grain and gas temperatures. We also develop a Monte Carlo simulation to calculate the H2 formation efficiency and compare the results to our rate equation models. Our results are the following: (1) Depending on the barrier against chemisorption, we predict the efficiency of H2 formation for a wide range of grain and gas temperatures. (2) The Eley-Rideal mechanism has an impact on the H2 formation efficiency at high grain and gas temperatures. (3) The fact that we consider chemisorption in our model makes the rate equation and Monte Carlo approaches equivalent., Comment: in "Light, dust and chemical evolution", Journal of Physics: Conference Series

The rotational spectra of DCO+ and its two isotopomers D{13}^CO+ and DC{18}^O+, produced in a negative glow discharge cell, have been recorded in the 137-792 GHz region, which includes lines from J=2-1 up to J=11-10. The determined rotational and centrifugal distortion constants allow to predict the DCO+ rotational spectrum up to 1000 GHz with an accuracy of 1 part in 10^8 or better. This is important for kinematic studies of dense molecular cloud cores and for future far-infrared observations. We also report on the first detection of the hyperfine structure of the DCO+(1-0) line made at the IRAM 30m antenna, toward the quiescent starless cloud core L1512, in the Taurus Molecular Cloud. We point out that this is the first observation of the hyperfine splitting due to the deuteron. This allowed us to quantify the effects of the hyperfine splitting on the line width determination; if the hyperfine structure is not taken into account in the line fit, the DCO+(1-0) line width is overestimated by a significant factor (> 2)., Comment: 14 pages, 3 figures, accepted by A&A

We have carried out a molecular survey of the Class 0 IM protostar NGC 7129 -- FIRS 2 (hereafter FIRS 2) and the Herbig Be star LkH$\alpha$ 234 with the aim of studying the chemical evolution of the envelopes of intermediate-mass (IM) young stellar objects (YSOs). Both objects have similar luminosities (~500 Lsun) and are located in the same molecular cloud which minimizes the chemical differences due to different stellar masses or initial cloud conditions. Moreover, since they are located at the same distance, we have the same spatial resolution in both objects. A total of 17 molecular species (including rarer isotopes) have been observed in both objects and the structure of their envelopes and outflows is determined with unprecedent detail. Our results show that the protostellar envelopes are dispersed and warmed up during the evolution to become a pre-main sequence star. In fact, the envelope mass decreases by a factor >5 from FIRS 2 to LkH$\alpha$234, while the kinetic temperature increases from ~13K to 28K. On the other hand, there is no molecular outflow associated with LkH$\alpha$234. The molecular outflow seems to stop before the star becomes visible. These physical changes strongly affect the chemistry of their envelopes. Based on our results in FIRS2 and LkH$\alpha$ 234, we propose some abundance ratios that can be used as chemical clocks for the envelopes of IM YSOs. The SiO/CS, CN/N2H+, HCN/N2H+, DCO+/HCO+ and D2CO/DCO+ ratios are good diagnostics of the protostellar evolutionary stage., Comment: 24 pages, 17 figures

Recent observations suggest that the behaviour of tracer species such as N_2H+ and CS is significantly different in regions of high and low mass star formation. In the latter, N_2H+ is a good tracer of mass, while CS is not. Observations show the reverse to be true in high-mass star formation regions. We use a computational chemical model to show that the abundances of these and other species may be significantly altered by a period of accelerated collapse in high mass star forming regions. We suggest these results provide a potential explanation of the observations, and make predictions for the behaviour of other species., Comment: 5 pages, 3 figures To be published in the Astrophysical Journal

The gravitational collapse of a spherical cloud core is investigated by numerical calculations. The initial conditions of the core lie close to the critical Bonnor-Ebert sphere with a central density of \sim 10^4 cm^{-3} in one model (alpha=1.1), while gravity overwhelms pressure in the other (alpha=4.0), where alpha is the internal gravity-to-pressure ratio. The alpha=1.1 model shows reasonable agreement with the observed velocity field in prestellar cores. Molecular distributions in cores are calculated by solving a chemical reaction network that includes both gas-phase and grain-surface reactions. When the central density of the core reaches 10^5 cm^{-3}, carbon-bearing species are significantly depleted in the central region of the alpha=1.1 model, while the depletion is only marginal in the other model. The two different approaches encompass the observed variations of molecular distributions in different prestellar cores, suggesting that molecular distributions can be probes of contraction or accumulation time scales of cores. The central enhancement of the NH3/N2H+ ratio, which is observed in some prestellar cores, can be reproduced under certain conditions by adopting recently measured branching fractions for N2H+ recombination. Various molecular species, such as CH3OH and CO2, are produced by grain-surface reactions. The ice composition depends sensitively on the assumed temperature. Multi-deuterated species are included in our most recent gas-grain chemical network. The deuterated isotopomers of H3+ are useful as probes of the central regions of evolved cores, in which gas-phase species with heavy elements are strongly depleted. At 10 K, our model can reproduce the observed abundance ratio of ND3/NH3, but underestimates the isotopic ratios of deuterated to normal methanol., Comment: 38 pages, 9 figures

We have undertaken a survey of N2H+ and N2D+ towards 31 low-mass starless cores using the IRAM 30m telescope. Our main objective has been to determine the abundance ratio of N2D+ and N2H+ towards the nuclei of these cores and thus to obtain estimates of the degree of deuterium enrichment, a symptom of advanced chemical evolution according to current models. We find that the N(N2D+)/N(N2H+) ratio is larger in more "centrally concentrated cores" with larger peak H2 and N2H+ column density than the sample mean. The deuterium enrichment in starless cores is presently ascribed to depletion of CO in the high density (> 3*10^4 cm-3) core nucleus. To substantiate this picture, we compare our results with observations in dust emission at 1.2 mm and in two transitions of C18O. We find a good correlation between deuterium fractionation and N(C18O)/N(H2) for the nuclei of 14 starless cores. We, thus, identified a set of properties that characterize the most evolved, or "pre-stellar", starless cores. These are: higher N2H+ and N2D+ column densities, higher N(N2D+)/N(N2H+), more pronounced CO depletion, broader N2H+ lines with infall asymmetry, higher central H2 column densities and a more compact density profile than in the average core. We conclude that this combination of properties gives a reliable indication of the evolutionary state of the core. Seven cores in our sample (L1521F, OphD, L429, L694, L183, L1544 and TMC2) show the majority of these features and thus are believed to be closer to forming a protostar than are the other members of our sample. Finally, we note that the subsample of Taurus cores behaves more homogeneously than the total sample, an indication that the external environment could play an important role in the core evolution., Comment: 42 pages, 17 figures, 12 tables, published in ApJ. Figures 8 to 15 available in gif format; full resolution version can be downloaded at http://cfa-www.harvard.edu/~acrapsi/publications/ Replaced with final version as published on the journal

We report a theoretical study of sulphur chemistry, as applied to hot cores, where S-bearing molecular ratios have been previously proposed and used as chemical clocks. As in previous models, we follow the S-bearing molecular composition after the injection of grain mantle components into the gas phase. For this study, we developed a time-dependent chemical model with up-to-date reaction rate coefficients. We ran several cases, using different realistic chemical compositions for the grain mantles and for the gas prior to mantle evaporation. The modeling shows that S-bearing molecular ratios depend very critically on the gas temperature and density, the abundance of atomic oxygen, and, most importantly, on the form of sulphur injected in the gas phase, which is very poorly known. Consequently, ratios of S-bearing molecules cannot be easily used as chemical clocks. However, detailed observations and careful modeling of both physical and chemical structure can give hints on the source age and constrain the mantle composition (i.e. the form of sulphur in cold molecular clouds) and, thus, help to solve the mystery of the sulphur depletion. We analyse in detail the cases of Orion and IRAS16293-2422. The comparison of the available observations with our model suggests that the majority of sulphur released from the mantles is mainly in, or soon converted into, atomic form., Comment: Accepted for publication in A&A

We observed the pre-stellar core L1521F in dust emission at 1.2mm and in two transitions each of N2H+, N2D+, C18O, and C17O in order to increase the sample of well studied centrally concentrated and chemically evolved starless cores, likely on the verge of star formation, and to determine the initial conditions for low--mass star formation in the Taurus Molecular Cloud. We derived in this object a molecular hydrogen number density n(H2) ~ 10^6 cm-3 and a CO depletion factor, integrated along the line of sight, fD ~ 15 in the central 20'', similar to the pre-stellar core L1544. However, the N(N2D+)/N(N2H+) column density ratio is ~0.1, a factor of about 2 lower than that found in L1544. The observed relation between the deuterium fractionation and the integrated CO depletion factor across the core can be reproduced by chemical models if N2H+ is slightly (factor of ~2 in fractional abundance) depleted in the central 3000 AU. The N2H+ and N2D+ linewidths in the core center are ~0.3 km/s, significantly larger than in other more quiescent Taurus starless cores but similar to those observed in the center of L1544. The kinematical behaviour of L1521F is more complex than seen in L1544, and a model of contraction due to ambipolar diffusion is only marginally consistent with the present data. Other velocity fields, perhaps produced by unresolved substructure, are present. Both chemical and kinematical analyses suggest that L1521F is less evolved than L1544, but, in analogy with L1544, it is approaching the ``critical'' state., Comment: 15 pages, 17 figures, accepted by A&A (higher resolution version of images at http://www.arcetri.astro.it/~starform/publ2004.htm)

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, N2H+, CS, C34S, C18O, and C17O, 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 while N2H+ has a constant abundance over most of the cores. In contrast, both C18O 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., Comment: A&A accepted. 23 pages, 15 figures. v2: Minor typos corrected. v3: additional clarifications on turbulent models; final journal version

This paper explores the sensitivity of line profiles of H2D+, HCO+ and N2H+, observed towards the center of L1544, to various kinematic and chemical parameters. The total width of the H2D+ line can be matched by a static model and by models invoking ambipolar diffusion and gravitational collapse. The derived turbulent line width is b=0.15 km/s for the static case and <~ 0.05 km/s for the collapse case. However, line profiles of HC18O+ and N2H+ rule out the static solution. The double-peaked H2D+ line shape requires either infall speeds in the center that are much higher than predicted by ambipolar diffusion models, or a shell-type distribution of H2D+, as is the case for HCO+ and N2H+. At an offset of ~20 arcsec from the dust peak, the H2D+ abundance drops by a factor of ~5., Comment: four pages, two colour figures; to appear in The Dense Interstellar Medium in Galaxies, proceedings of the fourth Cologne-Bonn-Zermatt Symposium, Sept 22-26, 2003

We present optical spectra of eight candidate brown dwarfs and a previously known T Tauri star (Sz 33) of the Chamaeleon I dark cloud. We derived spectral types based on the strength of the TiO or VO absorption bands present in the spectra of these objects as well as on the PC3 index of Martin et al. (1999). Photometric data from the literature are used to estimate the bolometric luminosities for these sources. We apply D'Antona & Mazzitelli (1997) pre-main sequence evolutionary tracks and isochrones to derive masses and ages. Based on the presence of Halpha in emission, we confirm that most of the candidates are young objects. Our sample however includes two sources for which we can only provide upper limits for the emission in Halpha; whereas these two objects are most likely foreground/background stars, higher resolution spectra are required to confirm their true nature. Among the likely cloud members, we detect one new sub-stellar object and three transition stellar/sub-stellar sources., Comment: 22 pages - manuscript format

We present the results of N2H+(1-0) observations of 35 dense molecular cloud cores from the northern and southern hemispheres where massive stars and star clusters are formed. Line emission has been detected in 33 sources, for 28 sources detailed maps have been obtained. The optical depth of (23-12) component toward peak intensity positions of 10 sources is ~ 0.2-1. In total, 47 clumps have been revealed in 26 sources. Integrated intensity maps with aspect ratios < 2 have been fitted with a power-law radial distribution $r^{-p}$ convolved with the telescope beam. Mean power-law index is close to unity corresponding to the $\sim r^{-2}$ density profile provided N2H+ excitation conditions do not vary inside these regions. Line widths of the cores either decrease or stay constant with distance from the center. The ratio of rotational to gravitational energy is too low for rotation to play a significant role in the dynamics of the cores. A correlation between mean line widths and sizes of clumps has been found., Comment: 17 pages, Latex

We have detected the 372 GHz line of ortho-H2D+ towards the pre-stellar core L1544. The strongest emission (T_mb ~ 1 K) occurs at the peak of the millimeter continuum emission, while measurements at offset positions indicate that H2D+ is confined within ~ 20 arcsec, where CO is highly depleted. The derived H2D+ abundance of ~ 10^{-9} is comparable with previous estimates of the electron abundance in the core, which suggests that H2D+ is the main molecular ion in the central 20 arcsec (2800 AU) of L1544. This confirms the expectations that H2D+ is dramatically enhanced in gas depleted of molecules other than H2. The measured abundance even exceeds the present model predictions by about a factor ten. One possibility is that all CNO-bearing neutral species, including atomic oxygen, are almost completely (> 98%) frozen within a radius of ~2800 AU., Comment: A&A Letters, in press

We have observed the dark cloud Barnard 68 with SCUBA at 850 um and with SIMBA at 1.2 mm. The submillimetre and millimetre dust emission correlate well with the extinction map of Alves, Lada and Lada (2001).The A_V/850um correlation is clearly not linear and suggests lower temperatures for the dust in the inner core of the cloud. Assuming a model for the temperature gradient, we derive the cloud-averaged dust emissivities (normalised to the V-Band extinction efficiency) at 850 um and 1.2 mm. We find k_850um/k_V = 4.0 +/- 1.0 x 10^-5 and k_1.2mm/k_V = 9.0 +/- 3.0 x 10^-6. These values are compared with other determinations in this wavelength regime and with expectations for models of diffuse dust and grain growth in dense clouds., Comment: 8 pages, 4 figures, A&A accepted (Letter), referee format

We have used the IRAM 30-m and FCRAO 14-m telescopes to observe the molecular clumps associated with 12 ultracompact (UC) HII regions in the J=6-5, 8-7 and 13-12 rotational transitions of methyl-acetylene (CH3C2H). Under the assumption of LTE and optically thin emission, we have derived temperature estimates ranging from 30 to 56 K. We estimate that the clumps have diameters of 0.2-1.6 pc, H_2 densities of 10^5-10^6 {cm^{-3}}, and masses of 10^2-2 10^4 M_\odot. We compare these values with those obtained by other authors from different molecular tracers and find that the H_2 density and the temperature inside the clumps vary respectively like n_{H_2} ~ R^{-2.6} and T ~ R^{-0.5}, with R distance from the centre. We also find that the virial masses of the clumps are ~3 times less than those derived from the CH3C2H column densities: we show that a plausible explanation is that magnetic fields play an important role to stabilise the clumps, which are on the verge of gravitational collapse. Finally, we show that the CH3C2H line width increases for decreasing distance from the clump centre: this effect is consistent with infall in the inner regions of the clumps. We conclude that the clumps around UC HII regions are likely to be transient (~10^(5) yr) entities, remnants of isothermal spheres currently undergoing gravitational collapse: the high mass accretion rates (~10^{-2} M_\odot yr^{-1}) lead to massive star formation at the centre of such clumps., Comment: 15 pages, 11 figures, A & A in press

Deuterium fractionation in molecular ions, in particular HCO+, has been extensively used to estimate the degree of ionization in molecular clouds. This paper reviews recent work on ionization degree in homogeneous clouds. We will show that the N(DCO+)/N(HCO+) column density ratio furnishes a measurement of x(e) only in regions where CO is not significantly depleted, thus in the outer skirts of dense cloud cores. To probe x(e) deep inside the clouds, one has to gauge deuterium enhancement in molecular ions with parent species not affected by depletion (e.g. N2H+), and rely on chemical models which take into account the cloud density structure. Unlike N(DCO+)/N(HCO+), the N(N2D+)/N(N2H+) column density ratio is predicted to considerably increase with core evolution (and/or the amount of CO depletion), reaching large values (> 0.2) in cloud cores on the verge of forming a star., Comment: 25 pages, 5 figures, accepted for publication in the special issue on "Deuterium in The Universe" of Planetary and Space Science

The formation of singly and doubly deuterated isotopomers of formaldehyde and of singly, doubly, and multiply deuterated isotopomers of methanol on interstellar grain surfaces has been studied with a semi-empirical modified rate approach and a Monte Carlo method in the temperature range 10-20 K. Agreement between the results of the two methods is satisfactory for all major and many minor species throughout this range. If gas-phase fractionation can produce a high abundance of atomic deuterium, which then accretes onto grain surfaces, diffusive surface chemistry can produce large abundances of deuterated species, especially at low temperatures and high gas densities. Warming temperatures will then permit these surface species to evaporate into the gas, where they will remain abundant for a considerable period. We calculate that the doubly deuterated molecules CHD2OH and CH2DOD are particularly abundant and should be searched for in the gas phase of protostellar sources. For example, at 10 K and high density, these species can achieve up to 10-20% of the abundance of methanol., Comment: 27 pages, 3 figures, Planetary and Space Science, in press