Search

Context. Pre-stellar cores are the birthplaces of Sun-like stars and represent the initial conditions for the assembly of protoplanetary systems. Due to their short lifespans, they are rare. In recent efforts to increase the number of such sources identified in the Solar neighbourhood, we have selected a sample of 40 starless cores from the publicly available core catalogs of the Herschel Gould Belt survey. In this work, we focus on one of the sources that stands out for its high central density: Corona Australis 151. Aims. We use molecular lines that trace dense gas (n>=10^6 cm-3) to confirm the exceptionally high density of this object, to study its physical structure, and to understand its evolutionary stage. Methods. We detected the N2H+ 3-2 and 5-4 transitions, and the N2D+ 3-2, 4-3, and 6-5 lines with the APEX telescope. We use the Herschel continuum data to infer a spherically symmetric model of the core's density and temperature. This is used as input to perform non-local-thermodynamic-equilibrium radiative transfer to fit the observed five lines. Results. Our analysis confirms that this core is characterised by very high densities (a few x 10^7 cm-3 at the centre) and cold temperatures. We infer a high deuteration level of N2D+/N2H+=0.50, indicative of an advanced evolutionary stage. In the large bandwidth covered by the APEX data, we detect several other deuterated species, including CHD2OH, D2CO, and ND3. We also detect multiple sulphurated species that present broader lines with signs of high-velocity wings. Conclusions. The observation of high-velocity wings and the fact that the linewidths of N2H+ and N2D+ become larger with increasing frequency can be interpreted either as an indication of supersonic infall motions developing in the central parts of a very evolved pre-stellar core or as the signature of outflows from a very low luminosity object (VeLLO). *SHORTENED*, Comment: Accepted for publication in A&A (19 Feb. 2025)

Pre-stellar cores are the first steps in the process of star and planet formation. However, the dynamical and chemical evolution of pre-stellar cores is still not well understood. We aim at estimating the central density of the pre-stellar core IRAS16293E and at carrying out an inventory of molecular species towards the density peak of the core. We observed high-$J$ rotational transitions of N$_2$H$^+$ and N$_2$D$^+$, and several other molecular lines towards the dust emission peak using the Atacama Pathfinder EXperiment (APEX) telescope, and derived the density and temperature profiles of the core using far-infrared surface brightness maps from $Herschel$. The N$_2$H$^+$ and N$_2$D$^+$ lines were analysed by non-LTE radiative transfer modelling. Our best-fit core model consists in a static inner region, embedded in an infalling envelope with an inner radius of approximately 3000 au (21" at 141 pc). The observed high-J lines of N$_2$H$^+$ and N$_2$D$^+$ (with critical densities greater than 10$^6$ cm$^{-3}$) turn out to be very sensitive to depletion; the present single-dish observations are best explained with no depletion of N$_2$H$^+$ and N$_2$D$^+$ in the inner core. The N$_2$D$^+$/N$_2$H$^+$ ratio that best reproduces our observations is 0.44, one of the largest observed to date in pre-stellar cores. Additionally, half of the molecules that we observed are deuterated isotopologues, confirming the high-level of deuteration towards this source. Non-LTE radiative transfer modelling of N$_2$H$^+$ and N$_2$D$^+$ lines proved to be an excellent diagnostic of the chemical structure and dynamics of a pre-stellar core. Probing the physical conditions immediately before the protostellar collapse is a necessary reference for theoretical studies and simulations with the aim of understanding the earliest stages of star and planet formation and the time scale of this process., Comment: Accepted for publication in A&A

We constructed two new models for deuterium and spin-state chemistry for the purpose of modeling the low-temperature environment prevailing in starless and pre-stellar cores. The fundamental difference between the two models is in the treatment of ion-molecule proton-donation reactions of the form $\rm XH^+ + Y \longrightarrow X + YH^+$, which are allowed to proceed either via full scrambling or via direct proton hop, i.e., disregarding proton exchange. The choice of the reaction mechanism affects both deuterium and spin-state chemistry, and in this work our main interest is on the effect on deuterated ammonia. We applied the new models to the starless core H-MM1, where several deuterated forms of ammonia have been observed. Our investigation slightly favors the proton hop mechanism over full scrambling because the ammonia D/H ratios are better fit by the former model, although neither model can reproduce the observed $\rm NH_2D$ ortho-to-para ratio of 3 (the models predict a value of $\sim$2). Extending the proton hop scenario to hydrogen atom abstraction reactions yields a good agreement for the spin-state abundance ratios, but greatly overestimates the deuterium fractions of ammonia. However, one can find a reasonably good agreement with the observations with this model by increasing the cosmic-ray ionization rate over the commonly-adopted value of $\sim$$10^{-17}\,\rm s^{-1}$. We also find that the deuterium fractions of several other species, such as $\rm H_2CO$, $\rm H_2O$, and $\rm CH_3$, are sensitive to the adopted proton-donation reaction mechanism. Whether the full scrambling or proton hop mechanism dominates may be dependent on the reacting system, and new laboratory and theoretical studies for various reacting systems are needed to constrain chemical models., Comment: Accepted to A&A; 14 pages, 13 figures

Context: The study of dust emission at millimeter wavelengths is important to shed light on the dust properties and physical structure of pre-stellar cores, the initial conditions in the process of star and planet formation. Aims: Using two new continuum facilities, AzTEC at the LMT and MUSTANG-2 at the GBO, we aim to detect changes in the optical properties of dust grains as a function of radius for the well-known pre-stellar core L1544. Methods: We determine the emission profiles at 1.1 and 3.3 mm and examine whether they can be reproduced in terms of the current best physical models for L1544. We also make use of various tools to determine the radial distributions of the density, temperature, and the dust opacity in a self-consistent manner. Results: We find that our observations cannot be reproduced without invoking opacity variations. With the new data, new temperature and density profiles, as well as opacity variations across the core, have been derived. The opacity changes are consistent with the expected variations between uncoagulated bare grains, toward the outer regions of the core, and grains with thick ice mantles, toward the core center. A simple analytical grain growth model predicts the presence of grains of ~3-4 um within the central 2000 au for the new density profile., Comment: Accepted for publication in A&A

Aims. We study whether rotational excitation makes a difference to the abundances of the $\rm H_3^+$ isotopologs, including spin states, in physical conditions corresponding to starless cores and protostellar envelopes. Methods. We developed a new rate coefficient set for the $\rm H_3^+$ isotopologs, allowing for rotational excitation, using the state-to-state rate coefficients from Hugo et al. These new so-called species-to-species rate coefficients are compared with previously-used ground state-to-species rate coefficients. Results. The species-to-species and ground state-to-species model results differ at high density and toward increasing temperatures ($T > 10$ K). The species-to-species model predicts a lower $\rm H_3^+$ deuteration degree at high density owing to an increase of the rate coefficients of endothermic reactions that decrease deuteration. At 20 K the ground state-to-species model overestimates the abundance of $\rm H_2D^+$ by a factor of about two while the abundance of $\rm D_3^+$ can differ by an order of magnitude between the models. Spin-state abundance ratios are also affected, and the new model better reproduces recent observations of ortho and para $\rm H_2D^+$ and $\rm D_2H^+$. The applicability regime of the new rate coefficients depends on the critical densities of the various rotational transitions. Conclusions. The difference in the abundances of the $\rm H_3^+$ isotopologs predicted by the two models is negligible at 10 K but excited states are very important in studies of deuteration at higher temperatures, for example in protostellar envelopes. The species-to-species rate coefficients provide a more realistic approach to the chemistry of the $\rm H_3^+$ isotopologs than the ground state-to-species rate coefficients do, and so the former should be adopted in chemical models describing the chemistry of the $\rm H_3^+ + H_2$ reacting system., Comment: accepted to A&A; 22 pages, 8 figures, 2 appendices; arXiv abstract heavily redacted and somewhat modified

In dense and cold molecular clouds dust grains are surrounded by thick icy mantles. It is however not clear if dust growth and coagulation take place before the switch-on of a protostar. This is an important issue, as the presence of large grains may affect the chemical structure of dense cloud cores, including the dynamically important ionization fraction, and the future evolution of solids in protoplanetary disks. To study this further, we focus on L1544, one of the most centrally concentrated pre-stellar cores on the verge of star formation, and with a well-known physical structure. We observed L1544 at 1.2 and 2 mm using NIKA, a new receiver at the IRAM 30 m telescope, and we used data from the Herschel Space Observatory archive. We find no evidence of grain growth towards the center of L1544 at the available angular resolution. Therefore, we conclude that single dish observations do not allow us to investigate grain growth toward the pre-stellar core L1544 and high sensitivity interferometer observations are needed. We predict that dust grains can grow to 200 $\mu$m in size toward the central ~300 au of L1544. This will imply a dust opacity change by a factor of ~2.5 at 1.2 mm, which can be detected using the Atacama Large Millimeter and submillimeter Array (ALMA) at different wavelengths and with an angular resolution of 2"., Comment: 12 pages, 14 figures. Accepted for publication in A&A

The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density estimates. We investigate how the inclusion of the hyperfine splitting caused by the deuterium (D) nucleus affects the analysis of the rotational lines of NH$_2$D. We present 30m IRAM observations of the above mentioned lines, as well as APEX o/p-NH$_2$D observations of the 1$_{01}$-0$_{00}$ lines at 333 GHz. The hyperfine spectra are first analyzed with a line list that only includes the hyperfine splitting due to the $^{14}$N nucleus. We find inconsistencies between the line widths of the 1$_{01}$-0$_{00}$ and 1$_{11}$-1$_{01}$ lines, the latter being larger by a factor of $\sim$1.6$\pm0.3$. Such a large difference is unexpected given the two sets of lines are likely to originate from the same region. We next employ a newly computed line list for the o/p-NH$_2$D transitions, where the hyperfine structure induced by both nitrogen and deuterium nuclei is included. With this new line list, the analysis of the previous spectra leads to linewidths which are compatible. Neglecting the hyperfine structure owing to D leads to overestimate the linewidths of the o/p-NH$_2$D lines at 3 mm. The error for a cold molecular core is about 50\%. This error propagates directly to the column density estimate. It is therefore recommended to take into account the hyperfine splittings caused by both the $^{14}$N and D nuclei in any analysis relying on these lines.

Aims. We investigate the stability of non-isothermal Bonnor-Ebert spheres with a model that includes a self-consistent calculation of the gas temperature. This way we can discard the assumption of equality between the dust and gas temperatures, and study the stability as the gas temperature changes with chemical evolution of the gas. Methods. We use a gas-grain chemical model including a time-dependent treatment of depletion onto grain surfaces, which strongly influences the gas temperature as the main coolant, CO, depletes from the gas. Dust and gas temperatures are solved with radiative transfer. For comparison with previous work, we assume that the cores are deeply embedded in a larger external structure, corresponding to visual extinction $A_{\rm V}^{\rm ext}=10$ mag. Results. We find that the critical non-dimensional radius $\xi_1$ derived here is similar to our previous work where we assumed $T_{\rm dust}=T_{\rm gas}$; the $\xi_1$ values lie below the isothermal critical value $\xi_0\sim6.45$, but the difference is less than 10%. Chemical evolution does not affect notably the stability condition of low-mass cores (<0.75 $M_\odot$). For higher masses the decrease of cooling owing to CO depletion causes substantial temporal changes in the temperature and density profiles of the cores. In the mass range 1-2 $M_\odot$ , $\xi_1$ decreases with chemical evolution, whereas above 3 $M_\odot$ , $\xi_1$ instead increases. We also find that decreasing $A_{\rm V}^{\rm ext}$ increases the gas temperature especially when the gas is chemically old, causing $\xi_1$ to increase with respect to models with higher $A_{\rm V}^{\rm ext}$. The derived $\xi_1$ values are close to $\xi_0$. The density contrast between the core center and edge varies between 8 to 16 depending on core mass and the chemical age of the gas, compared to the constant value $\sim$ 14.1 for the isothermal BES., Comment: 7 pages, 5 figures; accepted for publication in A&A; abstract (heavily) abridged for arXiv

Aims. We aim to develop a chemical model that contains a consistent description of spin-state chemistry in reactions involving chemical species with multiple deuterons. We apply the model to the specific case of deuterated ammonia, to derive values for the various spin-state ratios. Methods. We apply symmetry rules in the complete scrambling assumption to calculate branching ratio tables for reactions between chemical species that include multiple protons and/or deuterons. Reaction sets for both gas-phase and grain-surface chemistry are generated using an automated routine that forms all possible spin-state variants of any given reaction with up to six H/D atoms. Single-point and modified Bonnor-Ebert models are used to study the density and temperature dependence of ammonia and its isotopologs, and the associated spin-state ratios. Results. We find that the spin-state ratios of the ammonia isotopologs are, at late times, very different from their statistical values. The ratios are rather insensitive to variations in the density, but present strong temperature dependence. We derive high peak values ($\sim$ 0.1) for the deuterium fraction in ammonia, in agreement with previous (gas-phase) models. The deuterium fractionation is strongest at high density, corresponding to a high degree of depletion, and also presents temperature dependence. We find that in the temperature range 5 to 20 K, the deuterium fractionation peaks at $\sim$ 15 K while most of the ortho/para (and meta/para for $\rm ND_3$) ratios present a minimum at 10 K (ortho/para $\rm NH_2D$ has instead a maximum at this temperature). Conclusions. Owing to the density and temperature dependence found in the abundances and spin-state ratios of ammonia and its isotopologs, it is evident that observations of ammonia and its deuterated forms can provide important constraints on the physical structure of molecular clouds., Comment: 13 pages, 13 figures, 6 tables; accepted for publication in A&A; arXiv abstract abridged

Context. he study of prestellar cores is essential to understanding the initial stages of star formation. With $Herschel$ more cold clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 $Herschel$ fields observed as a follow-up on original $Planck$ detections. We have observed these clumps in $^{13}$CO (1-0), C$^{18}$O (1-0), and N$_2$H$^+$ (1-0) lines. Aims. Our aim is to find out if these cold clumps are prestellar. We have examined to what extent independent analysis of the dust and the molecular lines lead to similar conclusions about the masses of these objects. Methods. We calculate the clump masses and densities from the dust continuum and molecular line observations and compare these to each other and to the virial and Bonnor-Ebert masses calculated for each clump. Finally we examine two of the fields with radiative transfer models to estimate CO abundances. Results. When excitation temperatures could be estimated, the column densities derived from molecular line observations were comparable to those from dust continuum data. The median column density estimates are 4.2$\times 10^{21}$cm$^{-2}$ and 5.5$\times 10^{21}$cm$^{-2}$ for the line and dust emission data, respectively. The calculated abundances, column densities, volume densities, and masses all have large uncertainties and one must be careful when drawing conclusions. Abundance of $^{13}$CO was found in modeling the two clumps in the field G131.65$+$9.75 to be close to the usual value of 10$^{-6}$. The abundance ratio of $^{13}$CO and C$^{18}$O was $\sim$10. Molecular abundances could only be estimated with modeling, relying on dust column density data. Conclusions. The results indicate that most cold clumps, even those with dust color temperatures close to 11 K, are not necessarily prestellar., Comment: Accepted to A&A

Aims. We aim to present simulated chemical abundance profiles for a variety of important species, with special attention given to spin-state chemistry, in order to provide reference results against which present and future models can be compared. Methods. We employ gas-phase and gas-grain models to investigate chemical abundances in physical conditions corresponding to starless cores. To this end, we have developed new chemical reaction sets for both gas-phase and grain-surface chemistry, including the deuterated forms of species with up to six atoms and the spin-state chemistry of light ions and of the species involved in the ammonia and water formation networks. The physical model is kept simple in order to facilitate straightforward benchmarking of other models against the results of this paper. Results. We find that the ortho/para ratios of ammonia and water are similar in both gas-phase and gas-grain models, at late times in particular, implying that the ratios are determined by gas-phase processes. We derive late-time ortho/para ratios of ~0.5 and ~1.6 for ammonia and water, respectively. We find that including or excluding deuterium in the calculations has little effect on the abundances of non-deuterated species and on the ortho/para ratios of ammonia and water, especially in gas-phase models where deuteration is naturally hindered owing to the presence of abundant heavy elements. Although we study a rather narrow temperature range (10-20 K), we find strong temperature dependence in, e.g., deuteration and nitrogen chemistry. For example, the depletion timescale of ammonia is significantly reduced when the temperature is increased from 10 to 20 K; this is because the increase in temperature translates into increased accretion rates, while the very high binding energy of ammonia prevents it from being desorbed at 20 K., Comment: 24 pages, 18 figures; accepted for publication in A&A; abstract abridged

Cha-MMS1 was mapped in the NH_3(1,1) line and the 1.2 cm continuum using the Australia Telescope Compact Array, ATCA. The angular resolution of the ATCA observations is 7'' (~ 1000 AU), and the velocity resolution is 50 m s^{-1}. The core was also mapped with the 64-m Parkes telesope in the NH_3(1,1) and (2,2) lines. Observations from Herschel Space Observatory and Spitzer Space telescope were used to help interpretation. A compact high column density core with a steep velocity gradient is detected in ammonia, with a fractional ammonia abundance compatible with determinations towards other dense cores. The direction of the velocity gradient agrees with previous single-dish observations, and the overall velocity distribution can be interpreted as rotation. The rotation axis goes through the position of a compact far-infrared source detected by Spitzer and Herschel. The specific angular momentum of the core is typical for protostellar envelopes. A string of 1.2 cm continuum sources is tentatively detected near the rotation axis. The ammonia spectra suggest the presence of warm embedded gas in its vicinity. An hourglass-shaped structure is seen in ammonia at the cloud's average LSR velocity, also aligned with the rotation axis. Although this structure resembles a pair of outflow lobes the ammonia spectra show no indications of shocked gas. The observed ammonia structure mainly delineates the inner envelope around the central source. The velocity gradient is likely to originate in the angular momentum of the contracting core, although influence of the outflow from the neighbouring young star IRS4 is possibly visible on one side of the core. The tentative continuum detection and the indications of a warm background component near the rotation axis suggest that the core contains a deeply embedded outflow which may have been missed in previous single-dish CO surveys owing to beam dilution., Comment: 15 pages, 12 figures. Submitted to Astronomy & Astrophysics

Aims: We aim to investigate the abundances of light deuterium-bearing species such as HD, H2D+ and D2H+ in a gas-grain chemical model including an extensive description of deuterium and spin state chemistry, in physical conditions appropriate to the very centers of starless cores. Methods: We combine a gas-grain chemical model with radiative transfer calculations to simulate density and temperature structure in starless cores. The chemical model includes deuterated forms of species with up to 4 atoms and the spin states of the light species H2, H2+ and H3+ and their deuterated forms. Results: We find that HD eventually depletes from the gas phase because deuterium is efficiently incorporated to grain-surface HDO, resulting in inefficient HD production on grains. HD depletion has consequences not only on the abundances of e.g. H2D+ and D2H+, whose production depends on the abundance of HD, but also on the spin state abundance ratios of the various light species, when compared with the complete depletion model where heavy elements do not influence the chemistry. Conclusions: While the eventual HD depletion leads to the disappearance of light deuterium-bearing species from the gas phase in a relatively short timescale at high density, we find that at late stages of core evolution the abundances of H2D+ and D2H+ increase toward the core edge and the disributions become extended. The HD depletion timescale increases if less oxygen is initially present in the gas phase, owing to chemical interaction between the gas and the dust predecing the starless core phase. Our results are greatly affected if H2 is allowed to tunnel on grain surfaces, and therefore more experimental data not only on tunneling but also on the O + H2 surface reaction in particular is needed., Comment: 14 pages, 12 figures, abstract abridged; accepted for publication in A & A

Chemical reactions in starless molecular clouds are heavily dependent on interactions between gas phase material and solid phase dust and ices. We have observed the abundance and distribution of molecular gases in the cold, starless core DC 000.4-19.5 (SL42) in Corona Australis using data from the Swedish ESO Submillimeter Telescope. We present column density maps determined from measurements of C18O(J=2-1,1-0) and N2H+(J=1-0) emission features. Herschel data of the same region allow a direct comparison to the dust component of the cloud core and provide evidence for gas phase depletion of CO at the highest extinctions. The dust color emperature in the core calculated from Herschel maps ranges from roughly 10.7 to 14.0 K. This range agrees with the previous determinations from Infrared Space Observatory and Planck observations. The column density profile of the core can be fitted with a Plummer-like density distribution approaching n(r) ~ r^{-2} at large distances. The core structure deviates clearly from a critical Bonnor-Ebert sphere. Instead, the core appears to be gravitationally bound and to lack thermal and turbulent support against the pressure of the surrounding low-density material: it may therefore be in the process of slow contraction. We test two chemical models and find that a steady-state depletion model agrees with the observed C18O column density profile and the observed N(C18O) versus AV relationship., Comment: 29 pages, 9 figures

The temperature is a central parameter affecting the chemical and physical properties of dense cores of interstellar clouds and their evolution to star formation. The chemistry and the dust properties are temperature dependent and the interpretation of observation requires the knowledge of the temperature and its variations. Measurement of the gas kinetic temperature is possible with molecular line spectroscopy, the ammonia molecule, NH3, being the most commonly used tracer. We want to determine the accuracy of the temperature estimates derived from ammonia spectra. The normal interpretation of NH3 observations assumes that all the hyperfine line components are tracing the same gas volume. In the case of temperature gradients they may be sensitive to different layers and cause errors in the optical depth and gas temperature estimates. We examine a series of spherical cloud models, 1.0 and 0.5 M_Sun Bonnor-Ebert spheres, with different radial temperature profiles. We calculate synthetic NH3 spectra and compare the derived column densities and temperatures to the true values. For high signal-to-noise observations, the estimated gas kinetic temperatures are within ~0.3 K of the real mass averaged temperature and the column densities are correct to within ~10%. When the S/N ratio of the (2,2) spectrum decreases below 10, the temperature errors are of the order of 1K but without a significant bias. When the density of the models is increased by a factor of a few, the results begin to show significant bias because of the saturation of the (1,1) main group. The ammonia spectra are found to be a reliable tracer of the mass averaged gas temperature. Because the radial temperature profiles of the cores are not well constrained, the central temperature could still differ from this value. If the cores are optically very thick, there are no guarantees of the accuracy., Comment: 7 pages, accepted to A&A

Aims: We aim to derive a stability condition for non-isothermal Bonnor-Ebert spheres and compare the physical properties of critical non-isothermal and isothermal gas spheres. These configurations can serve as models for prestellar cores before gravitational collapse. Methods: A stability condition for non-isothermal spheres is derived by constructing an expression for the derivative of boundary pressure with respect to core volume. The temperature distribution is determined by means of radiative transfer calculations. Based on the stability analysis, we derive the physical parameters of critical cores for the mass range 0.1 - 5.0 M_sun. In addition, the properties of roughly Jupiter-mass cores are briefly examined. Results: At the low-mass end the critical non-isothermal sphere has lower central density and a slightly larger physical radius than the corresponding isothermal sphere (i.e. one with the same mass and average temperature). The temperature decrease towards the core centre becomes steeper towards smaller masses as the central density becomes higher. The slope depends on the adopted dust model. We find that the critical dimensionless radius increases above the isothermal value xi_0 = 6.45 for very low-mass cores (< 0.2 M_sun). However, in the mass-range studied here the changes are within 5% from the isothermal value. Conclusions: The density structures of non-isothermal and isothermal Bonnor-Ebert spheres for a given mass are fairly similar. However, the present models predict clear differences in the average temperatures for the same physical radius. Especially for low-mass cores, the temperature gradient probably has implications on the chemistry and the observed line emission. We also find that hydrostatic Jupiter-mass cores with radii less than 100 AU would have very high boundary pressures compared with typical pressures in the interstellar space., Comment: Accepted for publication in A&A

Context: The depletion of heavy elements in cold cores of interstellar molecular clouds can lead to a situation where deuterated forms of H_3^+ are the most useful spectroscopic probes of the physical conditions. Aims: The aim is to predict the observability of the rotational lines of H_2D^+ and D_2H^+ from prestellar cores. Methods: Recently derived rate coefficients for the H_3^+ + H_2 isotopic system were applied to the "complete depletion" reaction scheme to calculate abundance profiles in hydrostatic core models. The ground-state lines of H_2D^+(o) (372 GHz) and D_2H^+(p) (692 GHz) arising from these cores were simulated. The excitation of the rotational levels of these molecules was approximated by using the state-to-state coefficients for collisions with H_2. We also predicted line profiles from cores with a power-law density distribution advocated in some previous studies. Results: The new rate coefficients introduce some changes to the complete depletion model, but do not alter the general tendencies. One of the modifications with respect to the previous results is the increase of the D_3^+ abundance at the cost of other isotopologues. Furthermore, the present model predicts a lower H_2D^+ (o/p) ratio, and a slightly higher D_2H^+ (p/o) ratio in very cold, dense cores, as compared with previous modelling results. These nuclear spin ratios affect the detectability of the submm lines of H_2D^+(o) and D_2H^+(p). The previously detected H_2D^+ and D_2H^+ lines towards the core I16293E, and the H_2D^+ line observed towards Oph D can be reproduced using the present excitation model and the physical models suggested in the original papers., Comment: 10 pages, 11 Figures; ver2: updated some of the Figures, added some references, added an entry to acknowledgements

The aims of this study are to determine the properties and spatial distribution of dense cores in Ori B9, and to estimate their ages and dynamical timescales. The cloud was mapped in the 870 micron continuum with APEX/LABOCA, and selected positions were observed in the lines of N2H+ and N2D+ using IRAM-30m. These were used together with our previous H2D+ observations. Moreover, archival FIR Spitzer/MIPS maps were combined with the LABOCA map to distinguish between pre- and protostellar cores, and to estimate the evolutionary stages of protostars. Twelve dense cores were detected at 870 micron in the Ori B9 cloud. The submm cores constitute ~4% of the total mass of the Ori B9 region. There is an equal number of pre- and protostellar cores. Two of the submm sources, which we call SMM 3 and SMM 4, are previously unknown Class 0 candidates. We found a moderate degree of deuteration in N2H+ (0.03-0.04). There is, furthermore, evidence for N2H+ depletion in the core SMM 4. We derive a relatively high degree of ionization (~10^-7) in the clump associated with IRAS 05405-0117. The ambipolar diffusion timescales for two of the cores are ~70-100 times longer than the free-fall time. The distribution and masses of dense cores in Ori B9 are similar to those observed in more active regions in Orion, where the statistical core properties have been explained by turbulent fragmentation. The 50/50 proportions of pre- and protostellar cores suggest that duration of the prestellar phase is comparable to the free-fall time. This timescale can be questioned, however, on the basis of chemical data on the IRAS 05405-0117 region. A possible explanation is that this survey samples only the densest, i.e., dynamically most advanced cores., Comment: 17 pages, 8 figures, 9 tables. Revised version accepted for publication in A&A (Figs. 1, 3, and 7 revised, minor revisions in Sect. 5.2, updated reference information)

The aim of this study is to investigate the nature of radio sources associated with young stellar objects (YSOs) belonging to the R CrA cluster. By combining the centimetre radio data with the wealth of shorter wavelength data accumulated recently we wish to refine estimates of the evolutionary stages of the YSOs. Fluxes and spectral indices for the brightest radio sources were derived from the observations at 3, 6, and 20 cm using the ATCA. Seven of detected sources can be assigned to YSOs, which have counterparts in the X-rays, infrared or submm. One of the YSOs, Radio Source 9, is a Class 0 candidate, and another, IRS 7B, is suggested to be in the Class 0/I transition stage. IRS 7B is associated with extended radio lobes at 6 and 20 cm. The lobes may have a gyrosynchrotron emission component, which could be understood in terms of Fermi accleration in shocks. The Class I objects detected here seem to be a mixed lot. One of these, the wide binary IRS 5, shows a negative spectral index, rapid variability, and a high degree of circular polarisation with $V/I\approx33$ % on one of the days of observation. These signs of magnetic activity suggest that at least one of the binary components has advanced beyond the Class I stage. The radio source without YSO assigment, Radio Source 5, has been suggested to be a brown dwarf. The radio properties, in particular its persistent strong emission, do not support this classification. The radio characteristics of the detected YSOs roughly agree with the scheme where the dominant emission mechanism changes with age. The heterogeneity of the Class I group can possibly be explained by a drastic decline in the jet activity during this stage, which also changes the efficiency of free-free absorption around the protostar., Comment: Accepted for publication in A&A (8 pages, 4 figures, 4 tables)

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 Astrophysics

Molecular oxygen, O2 has been expected historically to be an abundant component of the chemical species in molecular clouds and, as such, an important coolant of the dense interstellar medium. However, a number of attempts from both ground and from space have failed to detect O2 emission. The work described here uses heterodyne spectroscopy from space to search for molecular oxygen in the interstellar medium. The Odin satellite carries a 1.1 m sub-millimeter dish and a dedicated 119 GHz receiver for the ground state line of O2. Starting in 2002, the star forming molecular cloud core rho Oph A was observed with Odin for 34 days during several observing runs. We detect a spectral line at v(LSR) = 3.5 km/s with dv(FWHM) = 1.5 km/s, parameters which are also common to other species associated with rho Ohp A. This feature is identified as the O2 (N_J = 1_1 - 1_0) transition at 118 750.343 MHz. The abundance of molecular oxygen, relative to H2,, is 5E-8 averaged over the Odin beam. This abundance is consistently lower than previously reported upper limits., Comment: 10 pages, 11 figures, accepted for publication in A&A

Aims. The main goal of the present study is to determine the fractional SiO abundance in high-mass star-forming cores, and to investigate its dependence on the physical conditions. In this way we wish to provide constraints on the chemistry models concerning the formation of SiO in the gas phase or via grain mantle evaporation. The work addresses also CH3CCH chemistry as the kinetic temperature is determined using this molecule. Methods. We estimate the physical conditions of 15 high-mass star-forming cores and derive the fractional SiO and CH3CCH abundances in them by using spectral line and dust continuum observations with the SEST. Results. The kinetic temperatures as derived from CH3CCH range from 25 to 39 K. The SiO emission regions are extended and typically half of the integrated line emission comes from the velocity range traced out by CH3CCH emission. The upper limit of SiO abundance in this 'quiescent' gas component is ~10^-10. The average CH3CCH abundance is about 7 x 10^-9. It shows a shallow, positive correlation with the temperature, whereas SiO shows the opposite tendency. Conclusions. We suggest that the high CH3CCH abundance and its possible increase when the clouds get warmer is related to the intensified desorption of the chemical precursors of the molecule from grain surfaces. In contrast, the observed tendency of SiO does not support the idea that the evaporation of Si-containing species from the grain mantles would be important, and it contradicts with the models where neutral reactions with activation barriers dominate the SiO production. A possible explanation for the decrease is that warmer cores represent more evolved stages of core evolution with fewer high-velocity shocks and thus less efficient SiO replenishment., Comment: 19 pages with 6 figures, accepted for publication in Astronomy & Astrophysics

Aims. The purpose of this study is to investigate the distributions of the isomeric molecules HCN and HNC and estimate their abundance ratio in the protostellar core Cha-MMS1 located in Chamaeleon {\sc i}. Methods. The core was mapped in the J=1-0 rotational lines of HCN, HNC, and HN13C. The column densities of H13CN, HN13C, H15NC and NH3 were estimated towards the centre of the core. Results. The core is well delineated in all three maps. The kinetic temperature in the core, derived from the NH3 (1,1) and (2,2) inversion lines, is 12.1+/-0.1 K. The HN13C/H13CN column density ratio is between 3 and 4, i.e. similar to values found in several other cold cores. The HN13C/H15NC column density ratio is about 7. In case no 15N fractionation occurs in HNC (as suggested by recent modelling results), the HNC/HN13C abundance ratio is in the range 30-40, which indicates a high degree of 13C fractionation in HNC. Assuming no differential 13C fractionation the HCN and HNC abundances are estimated to be about 7E-10 and about 2E-9, respectively, the former being nearly two orders of magnitude smaller than that of NH3. Using also previously determined column densities in Cha-MMS1, we can put the most commonly observed nitrogenous molecules in the following order according to their fractional abundances: X(NH3) > X(HC3N) > X(HNC) > X(HCN) > X(N2H+). Conclusions. The relationships between molecular abundances suggest that Cha-MMS1 represents an evolved chemical stage, experiencing at present the 'late-time' cyanopolyyne peak. The possibility that the relatively high HNC/HCN ratio derived here is only valid for the $^{13}$C isotopic substitutes cannot be excluded on the basis of the present and other available data., Comment: 9 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Aims. The purpose of this study is to examine the prediction that the deuterated H3+ ion, H2D+, can be found exclusively in the coldest regions of molecular cloud cores. This is also a feasibility study for the detection of the ground-state line of ortho-H2D+ at 372 GHz with APEX. Methods. The 1(10)-1(11) transition of H2D+ at 372 GHz was searched towards selected positions in the massive star forming cloud OriB9, in the dark cloud L183, and in the low- to intermediate mass star-forming cloud R CrA. Results. The line was detected in cold, prestellar cores in the regions of OriB9 and L183, but only upper limits were obtained towards other locations which either have elevated temperatures or contain a newly born star. The H2D+ detection towards OriB9 is the first one in a massive star-forming region. The fractional ortho-H2D+ abundances (relative to H2) are estimated to be about 1.0E-10 in two cold cores in OriB9, and 3.0E-10 in the cold core of L183. Conclusions. The H2D+ detection in OriB9 shows that also massive star forming regions contain very cold prestellar cores which probably have reached matured chemical composition characterized, e.g., by a high degree of deuterium fractionation. Besides as a tracer of the interior parts of prestellar cores, H2D+ may therefore be used to put contraints on the timescales related to massive star formation., Comment: 5 pages, 1 figure, accepted for publication in Astronomy and Astrophysics

The feasibility of observing the C18O (3-2) spectral line in cold clouds with the APEX telescope has been tested. As the line at 329.330 GHz lies in the wing of a strong atmospheric H2O absorption it can be observed only at high altitude observatories. Using the three lowest rotational levels instead of only two helps to narrow down the physical properties of dark clouds and globules. The centres of two C18O maxima in the high latitude low mass star forming region CG 12 were mapped in C18O (3-2) and the data were analyzed together with spectral line data from the SEST. The T_MB(3-2)/T_MB(2-1) ratio in the northern C18O maximum, CG 12 N, is 0.8, and in the southern maximum, CG 12 S, ~2. CG 12 N is modelled as a 120'' diameter (0.4pc) cold core with a mass of 27 Msun. A small size maximum with a narrow, 0.8 kms-1, C18O (3-2) spectral line with a peak temperature of T_MB ~11 K was detected in CG 12 S. This maximum is modelled as a 60'' to 80'' diameter (~0.2pc) hot (80 K < Tex < 200 K) ~1.6 Msun clump. The source lies on the axis of a highly collimated bipolar molecular outflow near its driving source. This is the first detection of such a compact, warm object in a low mass star forming region., Comment: APEX A&A special issue, accepted

We investigate the relationship between the C18O column density and the visual extinction in Chamaeleon I and in a part of the Chamaeleon III molecular cloud. The C18O column densities, N(C18O), are calculated from J=1-0 rotational line data observed with the SEST telescope. The visual extinctions, Av, are derived using JHK photometry from the 2MASS survey and the NICER color excess technique. In contrast with the previous results of Hayakawa et al. (2001), we find that the average N(C18O) / Av ratios are similar in Cha I and Cha III, and lie close to values derived for other clouds, i.e. N(C18O) ~ 2 x 10^14 cm^-2 (Av - 2). We find, however, clear deviations from this average relationship towards individual clumps. Larger than average N(C18O) / Av ratios can be found in clumps associated with the active star forming region in the northern part of Ch a I. On the other hand, some regions in the relatively quiescent southern part of Cha I show smaller than average N(C18O) / Av ratios and also very shallow proportionality between N(C18O) and Av. The shallow proportionality suggests that C18O is heavily depleted in these regions. As the degree of depletion is proportional to the gas density, these regions probably contain very dense, cold cores, which do not stand out in CO mappings. A comparison with the dust temperature map derived from the ISO data shows that the most prominent of the potentially depleted cores indeed coincides with a dust temperature minimum. It seems therefore feasible to use N(C18O) and Av data together for identifying old, dense cores in large scale mappings., Comment: 15 pages, 8 figures. Accepted for publication in A&A. Replaced with a version with corrected illustrations

The Chamaeleon I dark cloud (Cha I} has been mapped in C18O J=1-0 with an angular resolution of 1 arcmin using the SEST telescope. The large scale structures previously observed with lower spatial resolution in the cloud turn into a network of clumpy filaments. The automatic Clumpfind routine developed by Williams et al. 1994 is used to identify individual clumps in a consistent way. Altogether 71 clumps were found and the total mass of these clumps is 230 Mo. The dense 'cores' detected with the NANTEN telescope (Mizuno et al. 1999) and the very cold cores detected in the ISOPHOT serendipity survey (Toth et al. 2000) form parts of these filaments but decompose into numerous 'clumps' The filaments are preferentially oriented at right angles to the large-scale magnetic field in the region. We discuss the cloud structure, the physical characteristics of the clumps and the distribution of young stars. The observed clump mass spectrum is compared with the predictions of the turbulent fragmentation model of Padoan & Nordlund 2002. An agreement is found if fragmentation has been driven by very large-scale hypersonic turbulence, and if it has had time to dissipate into modestly supersonic turbulence in the interclump gas by the present time., Comment: 15 pages, 11 figures, 1 table, to be published in A&A (Accepted)

Odin has successfully observed the molecular core rho Oph A in the 572.5 GHz rotational ground state line of ammonia, NH3 (J,K = 1,0 - 0,0). The interpretation of this result makes use of complementary molecular line data obtained from the ground (C17O and CH3OH) as part of the Odin preparatory work. Comparison of these observations with theoretical model calculations of line excitation and transfer yields a quite ordinary abundance of methanol, X(CH3OH) = 3e-9. Unless NH3 is not entirely segregated from C17O and CH3OH, ammonia is found to be significantly underabundant with respect to typical dense core values, viz. X(NH3) = 8e-10., Comment: 4 pages, 2 figures, 2 tables, to appear in Astron. Astrophys. Letters

We have mapped the submillimeter emission from the 1(10)-1(01) transition of ortho-water in the W3 star-forming region. A 5'x5' map of the W3 IRS4 and W3 IRS5 region reveals strong water lines at half the positions in the map. The relative strength of the Odin lines compared to previous observations by SWAS suggests that we are seeing water emission from an extended region. Across much of the map the lines are double-peaked, with an absorption feature at -39 km/s; however, some positions in the map show a single strong line at -43 km/s. We interpret the double-peaked lines as arising from optically thick, self-absorbed water emission near the W3 IRS5, while the narrower blue-shifted lines originate in emission near W3 IRS4. In this model, the unusual appearance of the spectral lines across the map results from a coincidental agreement in velocity between the emission near W3 IRS4 and the blue peak of the more complex lines near W3 IRS5. The strength of the water lines near W3 IRS4 suggests we may be seeing water emission enhanced in a photon-dominated region., Comment: Accepted to A&A Letters as part of the special Odin issue; 4 pages

Odin has successfully observed three regions in the Orion A cloud, i.e. Ori KL, Ori S and the Orion Bar, in the 572.5 GHz rotational ground state line of ammonia, ortho-NH3 (J,K) = (1,0) -> (0,0), and the result for the Orion Bar represents the first detection in an ammonia line. Several velocity components are present in the data. Specifically, the observed line profile from the Orion Bar can be decomposed into two components, which are in agreement with observations in high-J CO lines by Wilson et al. 2001. Using the source model for the Orion Bar by these authors, our Odin observation implies a total ammonia abundance of NH3/H2 = 5E-9., Comment: 4 pages with 2 figures, accepted for publication in A&A

The low-mass star formation region associated with the reflection nebula Cederblad 110 in the Chamaeleon I cloud was mapped with the Australian Telescope Compact Array (ATCA) at 6 and 3.5cm. Altogether 11 sources were detected, three of which are previously known low mass young stellar objects associated with the nebula: the illuminating star IRS2 (Class III, Einstein X-ray source CHX7), the brightest far-infrared source IRS4 (Class I), and the weak X-ray source CHX10a (Class III). The other young stellar objects in the region, including the Class 0 protostar candidate Cha-MMS1, were not detected. The radio spectral index of IRS4 (alpha = 1.7 +/- 0.3) is consistent with optically thick free-free emission arising from a dense ionized region, probably a jet-induced shock occurring in the circumstellar material. As the only Class I protostar with a 'thermal jet' IRS4 is the strongest candidate for the central source of the molecular outflow found previously in the region. The emission from IRS2 has a flat spectrum (alpha = 0.05 +/- 0.05) but shows no sign of polarization, and therefore its origin is likely to be optically thin free-free emission either from ionized wind or a collimated jet. The strongest source detected in this survey is a new compact object with a steep negative spectral index (-1.1) and a weak linear polarization (about 2 %), which probably represents a background radio galaxy., Comment: 7 pages, 2 Postscript figures. Accepted for publication in Astronomy & Astrophysics

We have observed the nearby isolated globule Barnard 68 (B68) in the (J,K)=(1,1) and (2,2) inversion lines of ammonia. The gas kinetic temperature derived from these is T=10+-1.2 K. The observed line-widths are almost thermal: DV=0.181+-0.003 km/s (DVtherm=0.164+-0.010 km/s), supporting the earlier hypothesis that B68 is in hydrostatic equilibrium. The kinetic temperature is an input parameter to the physical cloud model put forward recently, and we discuss the impact of the new value in this context., Comment: 4 pages, 2 figures, accepted for publication in A&A (Letters)

Two dense dark cloud cores representing different stages of dynamical evolution were observed in a number of molecular spectral lines. One of the cores, Cha- MMS1 in the Chamaeleon cloud I contains a Class 0 protostar, whereas the other, CrA C in the R Coronae Australis cloud, is pre-stellar. We find that the cores have very different chemical compositions. Cha-MMS1 exhibits characteristics of so-called `early-type' chemistry with high abundan- ces of carbon-chain molecules such as HC3N, but it also has a large N2H+ abundance, which is expected only to build up at later stages. In contrast, none of the carbon-chain molecules were detected in CrA C. CrA C has a higher SO abundance than Cha-MMS1, which implies that it is chemically `older' than Cha- MMS1. The most striking difference between the two cores is seen in the HC3N/SO abundance ratio, which is at least three orders of magnitude higher in Cha-MMS than in CrA C. This result is surprising since starless cores are usually thought to be chemically younger than star-forming cores. Because of the high N2H+ abundance, we suggest that Cha-MMS1 represents the `late-time cyanopolyyne peak'. This would also be a more natural explanation for the carbon-chain molecules than the `early-time' picture in view of the fact that the core is presently collapsing to form a star. Thus, the dynamical evolution in CrA C seems to have been very slow compared with that of Cha-MMS1, and we discuss possible reasons for this., Comment: To appear in A&A. 19 pages, 4 figures