We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

Deuterated interstellar molecules frequently have abundances relative to their main isotopologues much higher than the overall elemental D-to-H ratio in the cold dense interstellar medium. The H$_3^+$ and its isotopologues play a key role in the deuterium fractionation; however, the abundances of these isotopologues have not been measured empirically with respect to H$_3^+$ to date. Our aim was to constrain the relative abundances of H$_2$D$^+$ and D$_3^+$ in the cold outer envelope of the hot core CRL 2136 IRS 1. We carried out three observations targeting H$_3^+$ and its isotopologues using the spectrographs CRIRES at the VLT, iSHELL at IRTF, and EXES on board SOFIA. In addition, the CO overtone band at 2.3 $\mu$m was observed by iSHELL to characterize the gas on the line of sight. The H$_3^+$ ion was detected toward CRL 2136 IRS 1 as in previous observations. Spectroscopy of lines of H$_2$D$^+$ and D$_3^+$ resulted in non-detections. The 3$\sigma$ upper limits of $N({\rm H_2D^+})/N({\rm H_3^+})$ and $N({\rm D_3^+})/N({\rm H_3^+})$ are 0.24 and 0.13, respectively. The population diagram for CO is reproduced by two components of warm gas with the temperatures 58 K and 530 K, assuming a local thermodynamic equilibrium (LTE) distribution of the rotational levels. Cold gas ($<$20 K) makes only a minor contribution to the CO molecular column toward CRL 2136 IRS 1. The critical conditions for deuterium fractionation in a dense cloud are low temperature and CO depletion. Given the revised cloud properties, it is no surprise that H$_3^+$ isotopologues are not detected toward CRL 2136 IRS 1. The result is consistent with our current understanding of how deuterium fractionation proceeds.

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.

We present molecular line and dust continuum observations of a Planck-detected cold cloud, G074.11+00.11. The cloud consists of a system of curved filaments and a central star-forming clump. The clump is associated with several infrared sources and H2O maser emission. We aim to determine the mass distribution and gas dynamics within the clump, to investigate if the filamentary structure seen around the clump repeats itself on a smaller scale, and to estimate the fractions of mass contained in dense cores and filaments. The velocity distribution of pristine dense gas can be used to investigate the global dynamical state of the clump, the role of filamentary inflows, filament fragmentation and core accretion. We use molecular line and continuum observations from single dish observatories and interferometric facilities to study the kinematics of the region. The molecular line observations show that the central clump may have formed as a result of a large-scale filament collision. The central clump contains three compact cores. Assuming a distance of 2.3 kpc, based on Gaia observations and a three-dimensional extinction method of background stars, the mass of the central clump exceeds 700 solar masses, which is roughly 25% of the total mass of the cloud. Our virial analysis suggests that the central clump and all identified substructures are collapsing. We find no evidence for small-scale filaments associated with the cores. Our observations indicate that the clump is fragmented into three cores with masses in the range [10,50] solar masses and that all three are collapsing. The presence of an H2O maser emission suggests active star formation. However, the CO lines show only weak signs of outflows. We suggest that the region is young and any processes leading to star formation have just recently begun.

At the low temperatures ($\sim$10 K) and high densities ($\sim$100,000 H$_2$ molecules per cc) of molecular cloud cores and protostellar envelopes, a large amount of molecular species (in particular those containing C and O) freeze-out onto dust grain surfaces. It is in these regions that the deuteration of H$_3^+$ becomes very efficient, with a sharp abundance increase of H$_2$D$^+$ and D$_2$H$^+$. The multi-deuterated forms of H$_3^+$ participate in an active chemistry: (i) their collision with neutral species produces deuterated molecules such as the commonly observed N$_2$D$^+$, DCO$^+$ and multi-deuterated NH$_3$; (ii) their dissociative electronic recombination increases the D/H atomic ratio by several orders of magnitude above the D cosmic abundance, thus allowing deuteration of molecules (e.g. CH$_3$OH and H$_2$O) on the surface of dust grains. Deuterated molecules are the main diagnostic tools of dense and cold interstellar clouds, where the first steps toward star and protoplanetary disk formation take place. Recent observations of deuterated molecules are reviewed and discussed in view of astrochemical models inclusive of spin-state chemistry. We present a new comparison between models based on complete scrambling (to calculate branching ratio tables for reactions between chemical species that include protons and/or deuterons) and models based on non-scrambling (proton hop) methods, showing that the latter best agree with observations of NH$_3$ deuterated isotopologues and their different nuclear spin symmetry states.

We develop a new ''core field structure'' (CFS) model to predict the magnetic field strength and magnetic field fluctuation profile of dense cores using gas kinematics. We use spatially resolved observations of the nonthermal velocity dispersion from the Green Bank Ammonia survey along with column density maps from SCUBA-2 to estimate the magnetic field strength across seven dense cores located in the L1688 region of Ophiuchus. The CFS model predicts the profile of the relative field fluctuation, which is related to the observable dispersion in direction of the polarization vectors. Within the context of our model we find that all the cores have a transcritical mass-to-flux ratio.

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.

The occultation of the radio galaxy 0141+268 by the asteroid (372) Palma on 2017 May 15 was observed using six antennas of the Very Long Baseline Array (VLBA). The shadow of Palma crossed the VLBA station at Brewster, Washington. Owing to the wavelength used, and the size and the distance of the asteroid, a diffraction pattern in the Fraunhofer regime was observed. The measurement retrieves both the amplitude and the phase of the diffracted electromagnetic wave. This is the first astronomical measurement of the phase shift caused by diffraction. The maximum phase shift is sensitive to the effective diameter of the asteroid. The bright spot at the shadow's center, the so called Arago--Poisson spot, is clearly detected in the amplitude time-series, and its strength is a good indicator of the closest angular distance between the center of the asteroid and the radio source. A sample of random shapes constructed using a Markov chain Monte Carlo algorithm suggests that the silhouette of Palma deviates from a perfect circle by 26+-13%. The best-fitting random shapes resemble each other, and we suggest their average approximates the shape of the silhouette at the time of the occultation. The effective diameter obtained for Palma, 192.1+-4.8 km, is in excellent agreement with recent estimates from thermal modeling of mid-infrared photometry. Finally, our computations show that because of the high positional accuracy, a single radio interferometric occultation measurement can reduce the long-term ephemeris uncertainty by an order of magnitude.

Context: The $^{15}$N fractionation has been observed to show large variations among astrophysical sources, depending both on the type of target and on the molecular tracer used. These variations cannot be reproduced by the current chemical models. Aims: Until now, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ has been accurately measured in only one prestellar source, L1544, where strong levels of fractionation, with depletion in $^{15}$N, are found ($^{14}$N/$^{15}$N$\: \approx 1000$). In this paper we extend the sample to three more bona fide prestellar cores, in order to understand if the antifractionation in N$_2$H$^+$ is a common feature of this kind of sources. Methods: We observed N$_2$H$^+$,N$^{15}$NH$^+$ and $^{15}$NNH$^+$ in L183, L429 and L694-2 with the IRAM 30m telescope. We modeled the emission with a non-local radiative transfer code in order to obtain accurate estimates of the molecular column densities, including the one for the optically thick N$_2$H$^+$. We used the most recent collisional rate coefficients available, and with these we also re-analysed the L1544 spectra previously published. Results: The obtained isotopic ratios are in the range $630-770$ and significantly differ with the value, predicted by the most recent chemical models, of $\approx 440$, close to the protosolar value. Our prestellar core sample shows high level of depletion of $^{15}$N in diazenylium, as previously found in L1544. A revision of the N chemical networks is needed in order to explain these results.

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.

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".

We report on the detection of the ground-state rotational line of ortho-D2H+ at 1.477 THz (203 micron) using the German REceiver for Astronomy at Terahertz frequencies (GREAT) onboard the Stratospheric Observatory For Infrared Astronomy (SOFIA). The line is seen in absorption against far-infrared continuum from the protostellar binary IRAS 16293-2422 in Ophiuchus. The para-D2H+ line at 691.7 GHz was not detected with the APEX telescope toward this position. These D2H+ observations complement our previous detections of para-H2D+ and ortho-H2D+ using SOFIA and APEX. By modeling chemistry and radiative transfer in the dense core surrounding the protostars, we find that the ortho-D2H+ and para-H2D+ absorption features mainly originate in the cool (T<18 K) outer envelope of the core. In contrast, the ortho-H2D+ emission from the core is significantly absorbed by the ambient molecular cloud. Analyses of the combined D2H+ and H2D+ data result in an age estimate of ~500 000 yr for the core, with an uncertainty of ~200 000 yr. The core material has probably been pre-processed for another 500 000 years in conditions corresponding to those in the ambient molecular cloud. The inferred time scale is more than ten times the age of the embedded protobinary. The D2H+ and H2D+ ions have large and nearly equal total (ortho+para) fractional abundances of ~$10^{-9}$ in the outer envelope. This confirms the central role of H3+ in the deuterium chemistry in cool, dense gas, and adds support to the prediction of chemistry models that also D3+ should be abundant in these conditions.

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

Dense cloud cores present chemical differentiation due to the different distribution of C-bearing and N-bearing molecules, the latter being less affected by freeze-out onto dust grains. In this letter we show that two C-bearing molecules, CH$_3$OH and $c$-C$_3$H$_2$, present a strikingly different (complementary) morphology while showing the same kinematics toward the prestellar core L1544. After comparing their distribution with large scale H$_2$ column density N(H$_2$) map from the Herschel satellite, we find that these two molecules trace different environmental conditions in the surrounding of L1544: the $c$-C$_3$H$_2$ distribution peaks close to the southern part of the core, where the surrounding molecular cloud has a N(H$_2$) sharp edge, while CH$_3$OH mainly traces the northern part of the core, where N(H$_2$) presents a shallower tail. We conclude that this is evidence of chemical differentiation driven by different amount of illumination from the interstellar radiation field: in the South, photochemistry maintains more C atoms in the gas phase allowing carbon chain (such as $c$-C$_3$H$_2$) production; in the North, C is mainly locked in CO and methanol traces the zone where CO starts to freeze out significantly. During the process of cloud contraction, different gas and ice compositions are thus expected to mix toward the central regions of the core, where a potential Solar-type system will form. An alternative view on carbon-chain chemistry in star-forming regions is also provided.

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.

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.

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.

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.

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.

(abridged) We correlated near-infrared stellar H-Ks colour excesses of background stars from NTT/SOFI with the far-IR optical depth map, tauFIR, derived from Herschel 160, 250, 350, and 500 um data. The Herschel maps were also used to construct a model for the cloud to examine the effect of temperature gradients on the estimated optical depths and dust absorption cross-sections. A linear correlation is seen between the colour H-Ks and tauFIR up to high extinctions (AV ~ 25). The correlation translates to the average extinction ratio A250um/AJ = 0.0014 +/- 0.0002, assuming a standard near-infrared extinction law and a dust emissivity index beta=2. Using an empirical NH/AJ ratio we obtain an average absorption cross-section per H nucleus of sigmaH(250um) = (1.8 +/- 0.3) * 10^(-25) cm^2 / H-atom, corresponding to a cross-section per unit mass of gas kappaG(250 um) = 0.08 +/- 0.01 cm^2 / g. The cloud model however suggests that owing to the bias caused by temperature changes along the line-of-sight these values underestimate the true cross-sections by up to 40% near the centre of the core. Assuming that the model describes the effect of the temperature variation on tauFIR correctly, we find that the relationship between H-Ks and tauFIR agrees with the recently determined relationship between sigmaH and NH in Orion A. The derived far-IR cross-section agrees with previous determinations in molecular clouds with moderate column densities, and is not particularly large compared with some other cold cores. We suggest that this is connected to the core not beng very dense (the central density is likely to be ~10^5 cm^-3) and judging from previous molecular line data, it appears to be at an early stage of chemical evolution.

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.

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.

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.

We aim to further constrain the properties and evolutionary stages of dense cores in Orion B9. The central part of Orion B9 was mapped at 350 micron with APEX/SABOCA. A sample of nine cores in the region were observed in C17O(2-1), H13CO+(4-3) (towards 3 sources), DCO+(4-3), N2H+(3-2), and N2D+(3-2) with APEX/SHFI. These data are used in conjunction with our previous APEX/LABOCA 870-micron dust continuum data. Many of the LABOCA cores show evidence of substructure in the higher-resolution SABOCA image. In particular, we report on the discovery of multiple very low-mass condensations in the prestellar core SMM 6. Based on the 350-to-870 micron flux density ratios, we determine dust temperatures of ~7.9-10.8 K, and dust emissivity indices of ~0.5-1.8. The CO depletion factors are in the range ~1.6-10.8. The degree of deuteration in N2H+ is ~0.04-0.99, where the highest value (seen towards the prestellar core SMM 1) is, to our knowledge, the most extreme level of N2H+ deuteration reported so far. The level of HCO+ deuteration is about 1-2%. We also detected D2CO towards two sources. The detection of subcondensations within SMM 6 shows that core fragmentation can already take place during the prestellar phase. The origin of this substructure is likely caused by thermal Jeans fragmentation of the elongated parent core. A low depletion factor and the presence of gas-phase D2CO in SMM 1 suggest that the core chemistry is affected by the nearby outflow. The very high N2H+ deuteration in SMM 1 is likely to be remnant of the earlier CO-depleted phase.

The gas kinetic temperature in the centres of starless, high-density cores is predicted to fall as low as 5-6 K. The aim of this study was to determine the kinetic temperature distribution in the low-mass prestellar core Oph D where previous observations suggest a very low central temperature. The densest part of the Oph D core was mapped in the NH3(1,1) and (2,2) inversion lines using the Very Large Array (VLA). The physical quantities were derived from the observed spectra by fitting the hyperfine structure of the lines, and subsequently the temperature distribution of Oph D was calculated using the standard rotational temperature techniques. A physical model of the cores was constructed, and the simulated spectra after radiative transfer calculations with a 3D Monte Carlo code were compared with the observed spectra. Temperature, density, and ammonia abundance of the core model were tuned until a satisfactory match with the observation was obtained. The high resolution of the interferometric data reveals that the southern part of Oph D comprises of two small cores. The gas kinetic temperatures, as derived from ammonia towards the centres of the southern and northern core are 7.4 and 8.9 K, respectively. The observed masses of the cores are only 0.2 M_Sun. Their potential collapse could lead to formation of brown dwarfs or low-mass stars.

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.

We observed the 9-cm Lambda-doubling lines of CH along the dense filament of TMC-1. The CH column densities were compared with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation of the observed OH/CH abundance ratios. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0e-8, to the southeast where it reaches a value of 2.0e-8. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios of ~ 16 - 20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase provides a natural explanation for the observation of "early-type" molecules in this region.

We aim to determine the physical and chemical properties of dense cores in Orion B9. We observed the NH3(1,1) and (2,2), and the N2H+(3-2) lines towards the submm peak positions. These data are used in conjunction with our LABOCA 870 micron dust continuum data. The gas kinetic temperature in the cores is between ~9.4-13.9 K. The non-thermal velocity dispersion is subsonic in most of the cores. The non-thermal linewidth in protostellar cores appears to increase with increasing bolometric luminosity. The core masses are very likely drawn from the same parent distribution as the core masses in Orion B North. Starless cores in the region are likely to be gravitationally bound, and thus prestellar. Some of the cores have a lower radial velocity than the systemic velocity of the region, suggesting that they are members of the "low-velocity part" of Orion B. The observed core-separation distances deviate from the corresponding random-like model distributions. The distances between the nearest-neighbours are comparable to the thermal Jeans length. The fractional abundances of NH3 and N2H+ in the cores are ~1.5-9.8x10^{-8} and ~0.2-5.9x10^{-10}, respectively. The NH3 abundance appears to decrease with increasing H2 column and number densities. The NH3/N2H+ column density ratio is larger in starless cores than in cores with embedded protostars. The core population in Orion B9 is comparable in physical properties to those in nearby low-mass star-forming regions. It is unclear if the origin of cores could be explained by turbulent fragmentation. On the other hand, many of the core properties conform with the picture of dynamic core evolution. The Orion B9 region has probably been influenced by the feedback from the nearby Ori OB 1b group, and the fragmentation of the parental cloud into cores could be caused by gravitational instability.

Infrared dark clouds (IRDCs) likely represent very early stages of high-mass star/star cluster formation. In this study, we aim to determine the physical properties and spatial distribution of dense clumps in the IRDC MSXDC G304.74+01.32 (G304.74), and bring these characteristics into relation to theories concerning the origin of IRDCs and their fragmentation into clumps and star-forming cores. G304.74 was mapped in the 870 $\mu$m dust continuum with the LABOCA bolometer on APEX. Archival MSX and IRAS infrared data were used to study the nature and properties of the submillimetre clumps within the cloud. There are 8 clumps within G304.74 which are not associated with mid-infrared (MIR) emission. Some of them are candidates for being/harbouring high-mass starless cores (HMSCs). We compared the clump masses and their spatial distribution in G304.74 with those in several other recently studied IRDCs. There is a high likelihood that the clump mass distributions in G304.74 and in several other IRDCs represent the samples of the same parent distribution. In most cases the spatial distributions of clumps in IRDCs do not deviate significantly from random distributions. This is consistent with the idea that the origin of IRDCs, and their further sub-fragmentation down to scales of clumps is caused by supersonic turbulence in accordance with results from giant molecular clouds.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.