Complex organic molecules (COMs) have been detected in a few Class 0 protostars but their origin is not well understood. Going beyond studies of individual objects, we want to investigate the origin of COMs in young protostars on a statistical basis. We use the CALYPSO survey performed with the IRAM PdBI to search for COMs at high angular resolution in a sample of 26 solar-type protostars, including 22 Class 0 and four Class I objects. Methanol is detected in 12 sources and tentatively in one source, which represents half of the sample. Eight sources (30%) have detections of at least three COMs. We find a strong chemical differentiation in multiple systems with five systems having one component with at least three COMs detected but the other component devoid of COM emission. The internal luminosity is found to be the source parameter impacting the most the COM chemical composition of the sources, while there is no obvious correlation between the detection of COM emission and that of a disk-like structure. A canonical hot-corino origin may explain the COM emission in four sources, an accretion-shock origin in two or possibly three sources, and an outflow origin in three sources. The CALYPSO sources with COM detections can be classified into three groups on the basis of the abundances of oxygen-bearing molecules, cyanides, and CHO-bearing molecules. These chemical groups correlate neither with the COM origin scenarii, nor with the evolutionary status of the sources if we take the ratio of envelope mass to internal luminosity as an evolutionary tracer. We find strong correlations between molecules that are a priori not related chemically (for instance methanol and methyl cyanide), implying that the existence of a correlation does not imply a chemical link. [abridged]

Methanol ice is embedded in interstellar ice mantles present in dense molecular clouds. We aim to measure the sputtering efficiencies starting from different ice mantles of varying compositions experimentally, in order to evaluate their potential impact on astrochemical models. The sputtering yields of complex organic molecules is of particular interest, since few mechanisms are efficient enough to induce a significant feedback to the gas phase. We irradiated methanol and carbon dioxide ice mixtures of varying ratios with swift heavy ions in the electronic sputtering regime. We monitored the evolution of the infrared spectra and the species released to the gas phase with a mass spectrometer. Methanol and 13C-methanol isotopologue were used to remove any ambiguity on the measured irradiation products. The sputtering of methanol embedded in carbon dioxide ice is an efficient process leading to the ejection of intact methanol in the gas phase. We establish that when methanol is embedded in a carbon-dioxide-rich mantle exposed to cosmic rays, a significant fraction is sputtered as intact molecules. The sputtered fraction follows the time-dependent bulk composition of the ice mantle, the latter evolving with time due to the radiolysis-induced evolution of the bulk. If methanol is embedded in a carbon dioxide ice matrix, as the analyses of the spectral shape of the CO2 bending mode observations in some lines of sight suggest, the overall methanol sputtering yield is higher than if embedded in a water ice mantle. The sputtering is increased by a factor close to the dominant ice matrix sputtering yield, which is about six times higher for pure carbon dioxide ice when compared to water ice. These experiments are further constraining the cosmic-ray-induced ice mantle sputtering mechanisms important role in the gas-phase release of complex organic molecules from the interstellar solid phase.

{CircumStellar Envelopes (CSEs) of stars are complex chemical objects for which theoretical models encounter difficulties in elaborating a comprehensive overview of the occurring chemical processes. Along with photodissociation, ion-neutral reactions and dissociative recombination might play an important role in controlling molecular growth in outer CSEs. The aim of this work is to provide experimental insights into pathways of photochemistry-driven molecular growth within outer CSEs to draw a more complete picture of the chemical processes occurring within these molecule-rich environments. A simplified CSE environment was therefore reproduced in the laboratory through gas-phase experiments exposing relevant gas mixtures to an Extreme UltraViolet (EUV) photon source. This photochemical reactor should ultimately allow us to investigate chemical processes and their resulting products occurring under conditions akin to outer CSEs. We used a recently developed EUV lamp coupled to the APSIS photochemical cell to irradiate CSE relevant gas mixtures of H$_2$, CO and N$_2$, at one wavelength, 73.6 nm. The detection and identification of chemical species in the photochemical reactor was achieved through in-situ mass spectrometry analysis of neutral and cationic molecules. We find that exposing CO-N$_2$-H$_2$ gas mixtures to EUV photons at 73.6 nm induces photochemical reactions that yield the formation of complex, neutral and ionic species. Our work shows that N$_2$H$^+$ can be formed through photochemistry along with highly oxygenated ion molecules like HCO$^+$ in CSE environments. We also observe neutral N-rich organic species including triazole and aromatic molecules. These results confirm the suitability of our experimental setting to investigate photochemical reactions and provide fundamental insights into the mechanisms of molecular growth in the outer CSEs.

Understanding where and when the mass of stars is determined is one of the fundamental, mostly unsolved, questions in astronomy. Here, we present the first results of GASTON, the Galactic Star Formation with NIKA2 large programme on the IRAM 30m telescope, that aims to identify new populations of low-brightness sources to tackle the question of stellar mass determination across all masses. In this paper, we focus on the high-mass star formation part of the project, for which we map a $\sim2$ deg$^2$ region of the Galactic plane around $l=24^\circ$ in both 1.2 mm and 2.0 mm continuum. Half-way through the project, we reach a sensitivity of 3.7 mJy/beam at 1.2mm. Even though larger than our target sensitivity of 2 mJy, the current sensitivity already allows the identification of a new population of cold, compact sources that remained undetected in any (sub-)mm Galactic plane survey so far. In fact, about 25% of the $\sim 1600$ compact sources identified in the 1.2 mm GASTON image are new detections. We present a preliminary analysis of the physical properties of the GASTON sources as a function of their evolutionary stage, arguing for a potential evolution of the mass distribution of these sources with time.

We have developed the first gas-grain chemical model for oxygen fractionation (also including sulphur fractionation) in dense molecular clouds, demonstrating that gas-phase chemistry generates variable oxygen fractionation levels, with a particularly strong effect for NO, SO, O2, and SO2. This large effect is due to the efficiency of the neutral 18O + NO, 18O + SO, and 18O + O2 exchange reactions. The modeling results were compared to new and existing observed isotopic ratios in a selection of cold cores. The good agreement between model and observations requires that the gas-phase abundance of neutral oxygen atoms is large in the observed regions. The S16O/S18O ratio is predicted to vary substantially over time showing that it can be used as a sensitive chemical proxy for matter evolution in dense molecular clouds.

Understanding the formation mechanisms of protoplanetary disks and multiple systems, and their pristine properties, is a key question for modern astrophysics. The properties of the youngest disks, embedded in rotating infalling protostellar envelopes, have largely remained unconstrained up to now. In the framework of the IRAM-PdBI CALYPSO survey, we have obtained sub-arcsecond observations of the dust continuum emission at 231 GHz and 94 GHz, for a sample of 16 solar-type Class 0 protostars. In an attempt to identify disk-like structures embedded at small scales in the protostellar envelopes, we model the dust continuum emission visibility profiles using both Plummer-like envelope models and envelope models including additional Gaussian disk-like components. Our analysis shows that in the CALYPSO sample, 11 of the 16 Class 0 protostars are better reproduced by models including a disk-like dust continuum component contributing to the flux at small scales, but less than 25% of these candidate protostellar disks are resolved at radii > 60 au. Including all available literature constraints on Class 0 disks at subarcsecond scales, we show that our results are representative: most (> 72% in a sample of 26 protostars) Class 0 protostellar disks are small and emerge only at radii < 60 au. Our multiplicity fraction at scales 100-5000 au is in global agreement with the multiplicity properties of Class I protostars at similar scales. We confront our observational constraints on the disk size distribution in Class 0 protostars to the typical disk properties from protostellar formation models. Because they reduce the centrifugal radius, and produce a disk size distribution peaking at radii <100 au during the main accretion phase, magnetized models of rotating protostellar collapse are favored by our observations.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

The characterization of dust properties in the interstellar medium (ISM) is key for star formation. Mass estimates are crucial to determine gravitational collapse conditions for the birth of new stellar objects in molecular clouds. However, most of these estimates rely on dust models that need further observational constraints from clouds to prestellar and protostellar cores. We present results of a study of dust emissivity changes based on mm-continuum data obtained with the NIKA camera at the IRAM-30m telescope. Observing dust emission at 1.15 mm and 2 mm allows us to constrain the dust emissivity index ($\beta$) in the Rayleigh-Jeans tail of the dust spectral energy distribution (SED) far from its peak emission, where the contribution of other parameters (i.e. dust temperature) is important. Focusing on the Taurus molecular cloud, a low-mass star-forming regions in the Gould Belt, we analyze the emission properties of several distinct objects in the B213 filament: three prestellar cores, two Class-0/I protostellar cores and one Class-II object. By means of the ratio of the two NIKA channel-maps, we show that in the Rayleigh-Jeans approximation the dust emissivity index varies among the objects. For one prestellar and two protostellar cores, we produce a robust study using Herschel data to constrain the dust temperature of the sources. By using the Abel transform inversion technique we get accurate radial $\beta$ profiles. We find systematic spatial variations of $\beta$ in the protostellar cores that is not observed in the prestellar core. While in the former case $\beta$ decreases toward the center, in the latter it remains constant. Moreover, $\beta$ appears anticorrelated with the dust temperature. We discuss the implication of these results in terms of dust grain evolution between pre- and protostellar cores.

The mass of prestellar cores is an essential ingredient to understand the onset of star formation in the core. The low level of emission from cold dust may keep parts of it hidden from observation. We aim to determine the fraction of core mass in the temperature range < 8 K that can be expected for typical low- and high-mass star formation regions. We calculate the dust temperature within standard spherically symmetric prestellar cores for a grid of density powerlaws in the outer core regions, core masses, and variation in the external multi-component radiation field. The dust is assumed to be composed of amorphous silicate and carbon, and variations of its optical properties are discussed. As measure for the distribution of cores and clumps, we use core mass functions derived for various environments. In view of the high densities in very cold central regions, dust and gas temperatures are assumed to be equal. We find that the fraction of mass with temperatures < 8 K in typical low- and high-mass cores is < 20\%. It is possible to obtain higher fractions of very cold gas by placing intermediate- or high-mass cores in a typical low-mass star formation environment. We show that the mass uncertainty arising from FIR to mm modeling of very cold dust emission is smaller than the mass uncertainty due to the unknown dust opacities.

Short version: We apply a new synergetic radiative transfer method: the derived 1D density profiles are both consistent with a cut through the Herschel PACS/SPIRE and JCMT SCUBA-2 continuum maps of L1689B and with a derived local interstellar radiation field. Choosing an appropriate cut along the filament major axis, we minimize the impact of the filament emission on the modeling. For the bulk of the core (5000-20000 au) an isothermal sphere model with a temperature of around 10 K provides the best fits. We show that the power law index of the density profile, as well as the constant temperature can be derived directly from the radial surface brightness profiles. For the inner region (< 5000 au), we find a range of densities and temperatures that are consistent with the surface brightness profiles and the local interstellar radiation field. Based on our core models, we find that pixel-by-pixel single temperature spectral energy distribution fits are incapable of determining dense core properties. We conclude that, to derive physical core properties, it is important to avoid an azimuthal average of core and filament. Correspondingly, derived core masses are too high since they include some mass of the filament, and might introduce errors when determining core mass functions. The forward radiative transfer methods also avoids the loss of information owing to smearing of all maps to the coarsest spatial resolution. We find the central core region to be colder and denser than estimated in recent inverse radiative transfer modeling, possibly indicating the start of star formation in L1689B.

Context. "Snow lines", marking regions where abundant volatiles freeze out onto the surface of dust grains, play an important role for planet growth and bulk composition in protoplanetary disks. They can already be observed in the envelopes of the much younger, low-mass Class 0 protostars that are still in their early phase of heavy accretion. Aims. We aim at using the information on the sublimation regions of different kinds of ices to understand the chemistry of the envelope, its temperature and density structure, and the history of the accretion process. Methods. As part of the CALYPSO IRAM Large Program, we have obtained observations of C$^{18}$O, N$_2$H$^+$ and CH$_3$OH towards nearby Class 0 protostars with the IRAM Plateau de Bure interferometer at sub-arcsecond resolution. For four of these sources we have modeled the emission using a chemical code coupled with a radiative transfer module. Results. We observe an anti-correlation of C$^{18}$O and N$_2$H$^+$ in NGC 1333-IRAS4A, NGC 1333-IRAS4B, L1157, and L1448C, with N$_2$H$^+$ forming a ring around the centrally peaked C$^{18}$O emission due to N$_2$H$^+$ being chemically destroyed by CO. The emission regions of models and observations match for a CO binding energy of 1200 K, which is higher than the binding energy of pure CO ices ($\sim$855 K). Furthermore, we find very low CO abundances inside the snow lines in our sources, about an order of magnitude lower than the total CO abundance observed in the gas on large scales in molecular clouds before depletion sets in. Conclusions. The high CO binding energy may hint at CO being frozen out in a polar ice environment like amorphous water ice or in non-polar CO$_2$-rich ice. The low CO abundances are comparable to values found in protoplanetary disks, which may indicate an evolutionary scenario where these low values are already established in the protostellar phase. (Abbr. Version)

Understanding the processes that could lead to enrichment of molecules in $^{15}$N atoms is of particular interest in order to shed light on the relatively large variations observed in the $^{14}$N/$^{15}$N ratio in various solar system environments. Currently, the sample of molecular clouds where $^{14}$N/$^{15}$N ratios have been measured is small and has to be enlarged in order to allow statistically significant studies. In particular, the N$_2$H$^+$ molecule currently shows the largest spread of $^{14}$N/$^{15}$N ratios in high-mass star forming regions. However, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ was obtained in only two low-mass star forming regions (L1544 and B1b). The current work extends this sample to a third dark cloud. We targeted the 16293E prestellar core, where the N$^{15}$NH$^+$ $J$=1-0 line was detected. Using a model previously developed for the physical structure of the source, we solved the molecular excitation with a non-local radiative transfer code. For that purpose, we computed specific collisional rate coefficients for the N$^{15}$NH$^+$-H$_2$ collisional system. As a first step of the analysis, the N$_2$H$^+$ abundance profile was constrained by reproducing the N$_2$H$^+$ $J$=1-0 and 3-2 maps. A scaling factor was then applied to this profile to match the N$^{15}$NH$^+$ $J$=1-0 spectrum. We derive a column density ratio N$_2$H$^+$ / N$^{15}$NH$^+$ = $330^{+170}_{-100}$. The current estimate $\sim$330 agrees with the value typical of the elemental isotopic ratio in the local ISM. It is however lower than in some other cores, where values as high as 1300 have been reported.

Complex organic molecules (COMs) are detected in many regions of the interstellar medium, including prestellar cores. However, their formation mechanisms in cold (~10 K) cores remain to this date poorly understood. The formyl radical HCO is an important candidate precursor for several O-bearing terrestrial COMs in cores, as an abundant building block of many of these molecules. Several chemical routes have been proposed to account for its formation, both on grain surfaces, as an incompletely hydrogenated product of H addition to frozen-out CO molecules, or in the gas phase, either the product of the reaction between H2CO and a radical, or as a product of dissociative recombination of protonated formaldehyde H2COH+. The detection and abundance determination of H2COH+, if present, could provide clues as to whether this latter scenario might apply. We searched for protonated formaldehyde H2COH+ in the prestellar core L1689B using the IRAM 30m telescope. The H2COH+ ion is unambiguously detected, for the first time in a cold (~10 K) source. The derived abundance agrees with a scenario in which the formation of H2COH+ results from the protonation of formaldehyde. We use this abundance value to constrain the branching ratio of the dissociative recombination of H2COH+ towards the HCO channel to ~10-30%. This value could however be smaller if HCO can be efficiently formed from gas-phase neutral-neutral reactions, and we stress the need for laboratory measurements of the rate constants of these reactions at 10 K. Given the experimental difficulties in measuring branching ratios experimentally, observations can bring valuable constraints on these values, and provide a useful input for chemical networks.

The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H and ND$_3$ isotopologues. These calculations are based on the most accurate NH$_3$--H$_2$ potential energy surface available, which has been modified to describe the geometrical changes induced by the nuclear substitutions. The dynamical calculations are performed within the close--coupling formalism and are carried out in order to provide rate coefficients up to a temperature of $T$ = 50K. For the various isotopologues/symmetries, we provide rate coefficients for the energy levels below $\sim$ 100 cm$^{-1}$. Subsequently, these new rate coefficients are used in astrophysical models aimed at reproducing the NH$_2$D, ND$_2$H and ND$_3$ observations previously reported towards the prestellar cores B1b and 16293E. We thus update the estimates of the corresponding column densities and find a reasonable agreement with the previous models. In particular, the ortho--to--para ratios of NH$_2$D and NHD$_2$ are found to be consistent with the statistical ratios.

The recent unexpected detection of terrestrial complex organic molecules in the cold (~ 10 K) gas has cast doubts on the commonly accepted formation mechanisms of these species. Standard gas-phase mechanisms are inefficient and tend to underproduce these molecules, and many of the key reactions involved are unconstrained. Grain-surface mechanisms, which were presented as a viable alternative, suffer from the fact that they rely on grain surface diffusion of heavy radicals, which is not possible thermally at very low temperatures. One of the simplest terrestrial complex organic molecules, methanol is believed to form on cold grain surfaces following from successive H atom additions on CO. Unlike heavier species, H atoms are very mobile on grain surfaces even at 10 K. Intermediate species involved in grain surface methanol formation by CO hydrogenation are the radicals HCO and CH3O, as well as the stable species formaldehyde H2CO. These radicals are thought to be precursors of complex organic molecules on grain surfaces. We present new observations of the HCO and CH3O radicals in a sample of prestellar cores and carry out an analysis of the abundances of the species HCO, H2CO, CH3O, and CH3OH, which represent the various stages of grain-surface hydrogenation of CO to CH3OH. The abundance ratios between the various intermediate species in the hydrogenation reaction of CO on grains are similar in all sources of our sample, HCO:H2CO:CH3O:CH3OH ~ 10:100:1:100. We argue that these ratios may not be representative of the primordial abundances on the grains but, rather, suggest that the radicals HCO and CH3O are gas-phase products of the precursors H2CO and CH3OH, respectively. Gas-phase pathways are considered and simple estimates of HCO and CH3O abundances are compared to the observations. Critical reaction rate constants, branching ratios, and intermediate species are finally identified.

High degrees of deuterium fractionation are commonly found in cold prestellar cores and in the envelopes around young protostars. As it brings strong constraints to chemical models, deuterium chemistry is often used to infer core history or molecule formation pathways. Whereas a large number of observations is available regarding interstellar deuterated stable molecules, relatively little is known about the deuteration of hydride radicals, as their fundamental rotational transitions are at high frequencies where the atmosphere is mostly opaque. Nitrogen hydride radicals are important species in nitrogen chemistry, as they are thought to be related to ammonia formation. Observations have shown that ammonia is strongly deuterated, with [NH2D]/[NH3] ~ 10%. Models predict similarly high [ND]/[NH] ratios, but so far only one observational determination of this ratio is available, towards the envelope of the protostar IRAS16293-2422. In order to test model predictions, we aim here at determining [ND]/[NH] in a dense, starless core. We observed NH and ND in 16293E with the HIFI spectrometer on board the Herschel Space Observatory as part of the CHESS guaranteed time key programme, and derived the abundances of these two species using a non-LTE non-local radiative transfer model. Both NH and ND are detected in the source, with ND in emission and NH in absorption against the continuum arising from the cold dust emission. Our model shows however that the ND emission and the NH absorption originate from different layers in the cloud, as further evidenced by their different velocities. In the central region of the core, we can set a lower limit to the [ND]/[NH] ratio of ~2%. This estimate is consistent with recent pure gas-phase models of nitrogen chemistry

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

Context. Using observations to deduce dust properties, grain size distribution, and physical conditions in molecular clouds is a highly degenerate problem. Aims. The coreshine phenomenon, a scattering process at 3.6 and 4.5 $\mu$m that dominates absorption, has revealed its ability to explore the densest parts of clouds. We want to use this effect to constrain the dust parameters. The goal is to investigate to what extent grain growth (at constant dust mass) inside molecular clouds is able to explain the coreshine observations. We aim to find dust models that can explain a sample of Spitzer coreshine data. We also look at the consistency with near-infrared data we obtained for a few clouds. Methods. We selected four regions with a very high occurrence of coreshine cases: Taurus-Perseus, Cepheus, Chameleon and L183/L134. We built a grid of dust models and investigated the key parameters to reproduce the general trend of surface bright- nesses and intensity ratios of both coreshine and near-infrared observations with the help of a 3D Monte-Carlo radiative transfer code. The grid parameters allow to investigate the effect of coagulation upon spherical grains up to 5 $\mu$m in size derived from the DustEm diffuse interstellar medium grains. Fluffiness (porosity or fractal degree), ices, and a handful of classical grain size distributions were also tested. We used the near- and mostly mid-infrared intensity ratios as strong discriminants between dust models. Results. The determination of the background field intensity at each wavelength is a key issue. In particular, an especially strong background field explains why we do not see coreshine in the Galactic plane at 3.6 and 4.5 $\mu$m. For starless cores, where detected, the observed 4.5 $\mu$m / 3.6 $\mu$m coreshine intensity ratio is always lower than $\sim$0.5 which is also what we find in the models for the Taurus-Perseus and L183 directions. Embedded sources can lead to higher fluxes (up to four times greater than the strongest starless core fluxes) and higher coreshine ratios (from 0.5 to 1.1 in our selected sample). Normal interstellar radiation field conditions are sufficient to find suitable grain models at all wavelengths for starless cores. The standard interstellar grains are not able to reproduce observations and, due to the multi-wavelength approach, only a few grain types meet the criteria set by the data. Porosity does not affect the flux ratios while the fractal dimension helps to explain coreshine ratios but does not seem able to reproduce near-infrared observations without a mix of other grain types. Conclusions. Combined near- and mid-infrared wavelengths confirm the potential to reveal the nature and size distribution of dust grains. Careful assessment of the environmental parameters (interstellar and background fields, embedded or nearby reddened sources) is required to validate this new diagnostic.

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

Recently discovered scattered light at 3-5 $\mu$m from low-mass cores (so-called "coreshine") reveals the presence of grains around 1 $\mu$m, which is larger than the grains found in the low-density interstellar medium. But only about half of the 100+ cores investigated so far show the effect. This prompts further studies on the origin of this detection rate. From the 3D continuum radiative transfer equation, we derive the expected scattered light intensity from a core placed in an arbitrary direction seen from Earth. We use the approximation of single scattering, consider extinction up to 2nd-order Taylor approximation, and neglect spatial gradients in the dust size distribution. The impact of the directional characteristics of the scattering on the detection of scattered light from cores is calculated for a given grain size distribution, and local effects like additional radiation field components are discussed. The surface brightness profiles of a core with a 1D density profile are calculated for various Galactic locations, and the results are compared to the approximate detection limits. We find that for optically thin radiation and a constant size distribution, a simple limit for detecting scattered light from a low-mass core can be derived that holds for grains with sizes smaller than 0.5 $\mu$m. The extinction by the core prohibits detection in bright parts of the Galactic plane, especially near the Galactic center. For scattered light received from low-mass cores with grain sizes beyond 0.5 $\mu$m, the directional characteristics of the scattering favors the detection of scattered light above and below the Galactic center, and to some extent near the Galactic anti-center. We identify the local incident radiation field as the major unknown causing deviations from this simple scheme.

Context. The largest grains (0.5-1 micron) in the interstellar size distribution are efficient in scattering near- and mid-infrared radiation. These wavelengths are therefore particularly well suited to probe the still uncertain high-end of the size distribution. Aims. We investigate the change in appearance of a typical low-mass molecular core from the Ks (2.2 micron) band to the Spitzer IRAC 3.6 and 8 micron bands, and compare with model calculations, which include variations of the grain size distribution. Methods. We combine Spitzer IRAC and ground-based near-infrared observations to characterize the scattered light observed at the near- and mid-infrared wavelengths from the core L260. Using a spherical symmetric model core, we perform radiative transfer calculations to study the impact of various dust size distributions on the intensity profiles across the core. Results. The observed scattered light patterns in the Ks and 3.6 micron bands are found to be similar. By comparison with radiative transfer models the two profiles place constraints on the relative abundance of small and large (more than 0.25 micron) dust grains. The scattered light profiles are found to be inconsistent with an interstellar silicate grain distribution extending only to 0.25 micron and large grains are needed to reach the observed fluxes and the flux ratios. The shape of the Ks band surface brightness profile limits the largest grains to 1-1.5 micron. Conclusions. In addition to observing coreshine in the Spitzer IRAC channels, the combination with ground-based near-infrared observations are suited to constrain the properties of large grains in cores.

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

Nitrogen bearing species are common tracers of the physical conditions in a wide variety of objects, and most remarkably in dark clouds. The reservoir of gaseous nitrogen is expected to be atomic or molecular, but none of the two species are observable in the dark gas. Their abundances therefore derive indirectly from those of N-bearing species through chemical modelling. The recent years have accumulated data which stress our incomplete understanding of the nitrogen chemistry in dark cloud conditions. To tackle this problem of the nitrogen chemistry in cold gas, we have revised the formation of nitrogen hydrides, which is initiated by the key reaction \ce{N+ + H2 -> NH+ + H}. We propose a new rate for this reaction which depends on the ortho-to-para ratio of H2. This new rate allows to reproduce the abundance ratios of the three nitrogen hydrides, NH, \ce{NH2}, and \ce{NH3}, observed towards IRAS16293-2422, provided that the channel leading to NH from the dissociative recombination of \ce{N2H+} is not closed at low temperature. The ortho-to-para ratio of H2 is constrained to O/P=$10^{-3}$ by the abundance ratio NH:NH2, which provides a new method to measure O/P. This work stresses the need for reaction rates at the low temperatures of dark clouds, and for branching ratios of critical dissociative recombination reactions.

Cold molecular clouds are the birthplaces of stars and planets, where dense cores of gas collapse to form protostars. The dust mixed in these clouds is thought to be made of grains of an average size of 0.1 micrometer. We report the widespread detection of the coreshine effect as a direct sign of the existence of grown, micrometer-sized dust grains. This effect is seen in half of the cores we have analyzed in our survey, spanning all Galactic longitudes, and is dominated by changes in the internal properties and local environment of the cores, implying that the coreshine effect can be used to constrain fundamental core properties such as the three-dimensional density structure and ages and also the grain characteristics themselves.

While unbiased surveys observable from ground-based telescopes have previously been obtained towards several high mass protostars, very little exists on low mass protostars. To fill up this gap, we carried out a complete spectral survey of the bands at 3, 2, 1 and 0.8 mm towards the solar type protostar IRAS16293-2422. The observations covered about 200\,GHz and were obtained with the IRAM-30m and JCMT-15m telescopes. Particular attention was devoted to the inter-calibration of the obtained spectra with previous observations. All the lines detected with more than 3 sigma and free from obvious blending effects were fitted with Gaussians to estimate their basic kinematic properties. More than 4000 lines were detected (with sigma \geq 3) and identified, yielding a line density of approximatively 20 lines per GHz, comparable to previous surveys in massive hot cores. The vast majority (~2/3) of the lines are weak and due to complex organic molecules. The analysis of the profiles of more than 1000 lines belonging 70 species firmly establishes the presence of two distinct velocity components, associated with the two objects, A and B, forming the IRAS16293-2422 binary system. In the source A, the line widths of several species increase with the upper level energy of the transition, a behavior compatible with gas infalling towards a ~1 Mo object. The source B, which does not show this effect, might have a much lower central mass of ~0.1 Mo. The difference in the rest velocities of both objects is consistent with the hypothesis that the source B rotates around the source A. This spectral survey, although obtained with single-dish telescope with a low spatial resolution, allows to separate the emission from 2 different components, thanks to the large number of lines detected. The data of the survey are public and can be retrieved on the web site http://www-laog.obs.ujf-grenoble.fr/heberges/timasss.

Nitrogen is the fifth most abundant element in the Universe, yet the gas-phase chemistry of N-bearing species remains poorly understood. Nitrogen hydrides are key molecules of nitrogen chemistry. Their abundance ratios place strong constraints on the production pathways and reaction rates of nitrogen-bearing molecules. We observed the class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI, covering most of the frequency range from 0.48 to 1.78~THz at high spectral resolution. The hyperfine structure of the amidogen radical o-NH2 is resolved and seen in absorption against the continuum of the protostar. Several transitions of ammonia from 1.2 to 1.8~THz are also seen in absorption. These lines trace the low-density envelope of the protostar. Column densities and abundances are estimated for each hydride. We find that NH:NH2:NH3=5:1:300. {Dark clouds chemical models predict steady-state abundances of NH2 and NH3 in reasonable agreement with the present observations, whilst that of NH is underpredicted by more than one order of magnitude, even using updated kinetic rates. Additional modelling of the nitrogen gas-phase chemistry in dark-cloud conditions is necessary before having recourse to heterogen processes.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

Theoretical arguments suggest that dust grains should grow in the dense cold parts of molecular clouds. Evidence of larger grains has so far been gathered in near/mid infrared extinction and millimeter observations. Interpreting the data is, however, aggravated by the complex interplay of density and dust properties (as well as temperature for thermal emission). We present new Spitzer data of L183 in bands that are sensitive and insensitive to PAHs. The visual extinction AV map derived in a former paper was fitted by a series of 3D Gaussian distributions. For different dust models, we calculate the scattered MIR radiation images of structures that agree agree with the AV map and compare them to the Spitzer data. The Spitzer data of L183 show emission in the 3.6 and 4.5 micron bands, while the 5.8 micron band shows slight absorption. The emission layer of stochastically heated particles should coincide with the layer of strongest scattering of optical interstellar radiation, which is seen as an outer surface on I band images different from the emission region seen in the Spitzer images. Moreover, PAH emission is expected to strongly increase from 4.5 to 5.8 micron, which is not seen. Hence, we interpret this emission to be MIR cloudshine. Scattered light modeling when assuming interstellar medium dust grains without growth does not reproduce flux measurable by Spitzer. In contrast, models with grains growing with density yield images with a flux and pattern comparable to the Spitzer images in the bands 3.6, 4.5, and 8.0 micron.

Context. The high degree of deuteration observed in some prestellar cores depends on the ortho-to-para H2 ratio through the H3+ fractionation. Aims. We want to constrain the ortho/para H2 ratio across the L183 prestellar core. This is mandatory to correctly describe the deuter- ation amplification phenomenon in depleted cores such as L183 and to relate the total (ortho+para) H2D+ abundance to the sole ortho-H2D+ column density measurement. Methods. To constrain this ortho/para H2 ratio and derive its profile, we make use of the N2D+ /N2H+ ratio and of the ortho-H2D+ observations performed across the prestellar core. We use two simple chemical models limited to an almost totally depleted core description. New dissociative recombination and trihydrogen cation-dihydrogen reaction rates (including all isotopologues) are presented in this paper and included in our models. Results. We estimate the H2D+ ortho/para ratio in the L183 cloud, and constrain the H2 ortho/para ratio : we show that it is varying across the prestellar core by at least an order of magnitude being still very high (~0.1) in most of the cloud. Our time-dependent model indicates that the prestellar core is presumably older than 1.5-2 x 10^5 years but that it may not be much older. We also show that it has reached its present density only recently and that its contraction from a uniform density cloud can be constrained. Conclusions. A proper understanding of deuteration chemistry cannot be attained without taking into account the whole ortho/para family of molecular hydrogen and trihydrogen cation isotopologues as their relations are of utmost importance in the global scheme. Tracing the ortho/para H2 ratio should also give useful constraints on the dynamical evolution of prestellar cores.

We study the molecular content and chemistry of a circumstellar disk surrounding the Herbig Ae star AB Aur at (sub-)millimeter wavelengths. Our aim is to reconstruct the chemical history and composition of the AB Aur disk and to compare it with disks around low-mass, cooler T Tauri stars. We observe the AB Aur disk with the IRAM Plateau de Bure Interferometer in the C- and D- configurations in rotational lines of CS, HCN, C2H, CH3OH, HCO+, and CO isotopes. Using an iterative minimization technique, observed columns densities and abundances are derived. These values are further compared with results of an advanced chemical model that is based on a steady-state flared disk structure with a vertical temperature gradient, and gas-grain chemical network with surface reactions. We firmly detect HCO+ in the 1--0 transition, tentatively detect HCN, and do not detect CS, C2H, and CH3OH. The observed HCO+ and 13CO column densities as well as the upper limits to the column densities of HCN, CS, C2H, and CH3OH are in good agreement with modeling results and those from previous studies. The AB Aur disk possesses more CO, but is less abundant in other molecular species compared to the DM Tau disk. This is primarily caused by intense UV irradiation from the central Herbig A0 star, which results in a hotter disk where CO freeze out does not occur and thus surface formation of complex CO-bearing molecules might be inhibited.

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

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

Aims: To constrain the ionization fraction in protoplanetary disks, we present new high-sensitivity interferometric observations of N$_2$H$^+$ in three disks surrounding DM Tau, LkCa 15, and MWC 480. Methods: We used the IRAM PdBI array to observe the N$_2$H$^+$ J=1-0 line and applied a $\chi^2$-minimization technique to estimate corresponding column densities. These values are compared, together with HCO$^+$ column densities, to results of a steady-state disk model with a vertical temperature gradient coupled to gas-grain chemistry. Results: We report two \ndhp detections for LkCa 15 and DM Tau at the $5 \sigma$ level and an upper limit for MWC 480. The column density derived from the data for LkCa 15 is much lower than previously reported. The [N$_2$H$^+$/HCO$^+$] ratio is on the order of 0.02--0.03. So far, HCO$^+$ remains the most abundant observed molecular ion in disks. Conclusions: All the observed values generally agree with the modelled column densities of disks at an evolutionary stage of a few million years (within the uncertainty limits), but the radial distribution of the molecules is not reproduced well. The low inferred concentration of N$_2$H$^+$ in three disks around low-mass and intermediate-mass young stars implies that this ion is not a sensitive tracer of the overall disk ionization fraction.

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

(Abriged) We present the results of millimeter observations and a suitable chemical and radiative transfer model of the AB Aur (HD 31293) circumstellar disk and surrounding envelope. The integral molecular content of this system is studied by observing CO, C$^{18}$O, CS, HCO$^+$, DCO$^+$, H$_2$CO, HCN, HNC, and SiO rotational lines with the IRAM 30-m antenna, while the disk is mapped in the HCO$^+$(1-0) transition with the Plateau de Bure interferometer. Using a flared disk model with a vertical temperature gradient and an isothermal spherical envelope model with a shadowed midplane and two unshielded cones together with a gas-grain chemical network, time-dependent abundances of observationally important molecules are calculated. Then a 2D non-LTE line radiative transfer code is applied to compute excitation temperatures of several rotational transitions of HCO$^+$, CO, C$^{18}$O, and CS molecules. We synthesize the HCO$^+$(1-0) interferometric map along with single-dish CO(2-1), C$^{18}$O(2-1), HCO$^+$(1-0), HCO$^+$(3-2), CS(2-1), and CS(5-4) spectra and compared them with the observations. Our disk model successfully reproduces observed interferometric HCO$^+$(1-0) data, thereby constraining the following disk properties: (1) the inclination angle $\iota=17^{+6}_{-3}\degr$, (2) the position angle $\phi=80\pm30\degr$, (3) the size $R_\mathrm{out}=400\pm200$ AU, (4) the mass $M_\mathrm{disk}=1.3\cdot10^{-2} M_{\sun}$ (with a factor of $\sim7$ uncertainty), and (5) that the disk is in Keplerian rotation. Furthermore, indirect evidence for a local inhomogeneity of the envelope at $\ga600$ AU is found...

Constraints on the density and thermal 3D structure of the dense molecular cloud core Rho Oph D are derived from a detailed 3D radiative transfer modeling. Two ISOCAM images at 7 and 15 micron are fitted simultaneously by representing the dust distribution in the core with a series of 3D Gaussian density profiles. Size, total density, and position of the Gaussians are optimized by simulated annealing to obtain a 2D column density map. The projected core density has a complex elongated pattern with two peaks. We propose a new method to calculate an approximate temperature in an externally illuminated complex 3D structure from a mean optical depth. This T(tau)-method is applied to a 1.3 mm map obtained with the IRAM 30m telescope to find the approximate 3D density and temperature distribution of the core Rho Oph D. The spatial 3D distribution deviates strongly from spherical symmetry. The elongated structure is in general agreement with recent gravo-turbulent collapse calculations for molecular clouds. We discuss possible ambiguities of the background determination procedure, errors of the maps, the accuracy of the T(tau)-method, and the influence of the assumed dust particle sizes and properties.

We present IRAM 30m and JCMT observations of HDO lines towards the solar-type protostar IRAS 16293-2422. Five HDO transitions have been detected on-source, and two were unfruitfully searched for towards a bright spot of the outflow of IRAS 16293-2422. We interpret the data by means of the Ceccarelli, Hollenbach and Tielens (1996) model, and derive the HDO abundance in the warm inner and cold outer parts of the envelope. The emission is well explained by a jump model, with an inner abundance of 1e-7 and an outer abundance lower than 1e-9 (3 sigma). This result is in favor of HDO enhancement due to ice evaporation from the grains in theinner envelope. The deuteration ratio HDO/H2O is found to be f_in=3% and f_out < 0.2% (3 sigma) in the inner and outer envelope respectively and therefore, the fractionation also undergoes a jump in the inner part of the envelope. These results are consistent with the formation of water in the gas phase during the cold prestellar core phase and storage of the molecules on the grains, but do not explain why observations of H2O ices consistently derive a H2O ice abundance of several 1e-5 to 1e-4, some two orders of magnitude larger than the gas phase abundance of water in the hot core around IRAS 16293-2422.

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

We present a survey of the formaldehyde emission in a sample of eight Class 0 protostars obtained with the IRAM and JCMT telescopes. The data have been analyzed with three different methods with increasing level of sophistication. We first analyze the observed emission in the LTE approximation, and derive rotational temperatures between 11 and 40 K, and column densities between 1 and 20 x 10^13 cm^-2. Second, we use a LVG code and derive larger kinetic temperatures, between 30 and 90 K, consistent with subthermally populated levels and densities from 1 to 6 x 10^5 cm^-3. The column densities from the LVG modeling are within a factor of 10 with respect to those derived in the LTE approximation. Finally, we analyze the observations based upon detailed models for the envelopes surrounding the protostars, using temperature and density profiles previously derived from continuum observations. We approximate the formaldehyde abundance across the envelope with a jump function, the jump occurring when the dust temperature reaches 100 K, the evaporation temperature of the grain mantles. The observed formaldehyde emission is well reproduced only if there is a jump, more than two orders of magnitude, in four sources. In the remaining four sources the data are consistent with a formaldehyde abundance jump, but the evidence is more marginal (~2 sigma). The inferred inner H2CO abundance varies between 1 x 10^-8 and 6 x 10^-6. We discuss the implications of these jumps for our understanding of the origin and evolution of ices in low mass star forming regions. Finally, we give predictions for the submillimeter H2CO lines, which are particularly sensitive to the abundance jumps.

This paper explores the sensitivity of line profiles of H2D+, HCO+ and N2H+, observed towards the center of L1544, to various kinematic and chemical parameters. The total width of the H2D+ line can be matched by a static model and by models invoking ambipolar diffusion and gravitational collapse. The derived turbulent line width is b=0.15 km/s for the static case and <~ 0.05 km/s for the collapse case. However, line profiles of HC18O+ and N2H+ rule out the static solution. The double-peaked H2D+ line shape requires either infall speeds in the center that are much higher than predicted by ambipolar diffusion models, or a shell-type distribution of H2D+, as is the case for HCO+ and N2H+. At an offset of ~20 arcsec from the dust peak, the H2D+ abundance drops by a factor of ~5.

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

We report the detection of D2CO in a sample of starless dense cores, in which we previously measured the degree of CO depletion. The deuterium fractionation is found extremely high, [D2CO]/[H2CO] ~ 1-10 %, similar to that reported in low-mass protostars. This provides convincing evidence that D2CO is formed in the cold pre-stellar cores, and later desorbed when the gas warms up in protostars. We find that the cores with the highest CO depletions have also the largest [D2CO]/[H2CO] ratios, supporting the theoretical prediction that deuteration increases with increasing CO depletion.

We present new results on CO depletion in a sample of nearby pre-stellar cores, based on observations of the millimeter C17O and C18O lines and the 1.3 mm dust emission with the IRAM 30m telescope. In most cases, the distribution of CO is much flatter than that of the dust, whereas other tracers, like N2H+, still probe the latter. In the centre of these objects, we estimate CO to be underabundant by a factor 4-15 depending on the cores. The CO underabundance is more pronounced in the central regions and appears to decrease with increasing distance from the core centre. This underabundance is most likely due to the freezing out of CO onto the dust grains in the cold, dense parts of the cores. We find evidence for an increase of the CO depletion degree with the core density.

We report on mid-infrared observations of the compact stellar cluster located in the proximity of SGR 1900+14, and the radio/X-ray position of this soft-gamma repeater. Observations were performed in May and June of 2001 when the bursting source was in an active state. At the known radio and X-ray position of the SGR we did not detect transient mid-IR activity, although the observations were performed only hours before and after an outburst in the high-energy band.

Observations of dark cloud cores have been carried out in the mid-infrared using ISOCAM and in the far-infrared using ISOPHOT, both aboard the Infrared Space Observatory. The cores are in most cases detected in emission at 200 and 170 micron, remain undetected at 90 micron and are seen in absorption against the diffuse mid-infrared background at 7 micron. The observations are consistent with the cores being pre-stellar and not having a central heating source, and yield core temperatures of ~ 11-13 K. We were able to determine the density structure of the cores up to radii that extend beyond the sensitivity limit of previous submillimetre studies. The column density profiles of the cores studied here flatten out in the centre, as shown by previous submillimetre continuum results, and interestingly, in a few cases, the profiles steepen with radius beyond ~ 5000-10000 AU and present sharp edges. This could indicate that these cores are decoupled from their parent molecular cloud and represent finite reservoirs of mass for subsequent accretion.

We present the results of a mid-infrared (7 micron) imaging survey of a sample of 24 starless dense cores carried out at an angular resolution of 6 arcsec with the ISOCAM camera aboard the Infrared Space Observatory (ISO). The targeted cores are believed to be pre-stellar in nature and to represent the initial conditions of low-mass, isolated star formation. In previous submillimeter dust continuum studies of such pre-stellar cores, it was found that the derived column density profiles did not follow a single power-law such as N[H2] \propto r^(-1) throughout their full extent but flattened out near their center. These submillimeter observations however could not constrain the density profiles at radii greater than ~ 10000 AU. The present absorption study uses ISOCAM's sensitivity to map these pre-stellar cores in absorption against the diffuse mid-infrared background. The goal was to determine their structure at radii that extend beyond the limits of sensitivity of the submillimeter continuum maps and at twiceas good an angular resolution. Among the 24 cores observed in our survey, a majority of them show deep absorption features. The starless cores studied here all show a column density profile that flattens in the center, which confirms the submillimeter emission results. Moreover, beyond a radius of ~ 5000-10000 AU, the typical column density profile steepens with distance from core center and gets steeper than N[H2] \propto r^(-1), until it eventually merges with the low-density ambient molecular cloud. At least three of the cores present sharp edges at R ~ 15000-30000 AU and appear to be decoupled from their parent clouds, providing finite reservoirs of mass for subsequent star formation.