[abridged] Understanding how the infalling gas redistribute most of its initial angular momentum inherited from prestellar cores before reaching the stellar embryo is a key question. Disk formation has been naturally considered as a possible solution to this "angular momentum problem". However, how the initial angular momentum of protostellar cores is distributed and evolves during the main accretion phase and the beginning of disk formation has largely remained unconstrained up to now. In the framework of the IRAM CALYPSO survey, we used high dynamic range C$^{18}$O (2-1) and N$_2$H$^+$ (1-0) observations to quantify the distribution of specific angular momentum along the equatorial axis in a sample of 12 Class 0 protostellar envelopes from scales ~50 to 10000 au. The radial distributions of specific angular momentum in the CALYPSO sample suggest two distinct regimes within protostellar envelopes: the specific angular momentum decreases as $j \propto r^{1.6 \pm 0.2}$ down to ~1600 au and then tends to become relatively constant around 6 $\times$ 10$^{-4}$ km s$^{-1}$ pc down to ~50 au. The values of specific angular momentum measured in the inner Class 0 envelopes, namely that of the material directly involved in the star formation process ($<$1600 au), is on the same order of magnitude as what is inferred in small T-Tauri disks. Thus, disk formation appears to be a direct consequence of angular momentum conservation during the collapse. Our analysis reveals a dispersion of the directions of velocity gradients at envelope scales $>$1600 au, suggesting that they may not be related to rotational motions of the envelopes. We conclude that the specific angular momentum observed at these scales could find its origin in core-forming motions (infall, turbulence) or trace an imprint of the initial conditions for the formation of protostellar cores.

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]

The formation of protoplanetary disks is not well understood. To understand how and when these disks are formed, it is crucial to characterize the kinematics of the youngest protostars at a high angular resolution. Here we study a sample of 16 Class 0 protostars to measure their rotation profile at scales from 50 to 500~au and search for Keplerian rotation. We used high-angular-resolution line observations obtained with the Plateau de Bure Interferometer as part of the CALYPSO large program. From $\mathrm{^{13}CO}$ ($J=2-1$), $\mathrm{C^{18}O}$ ($J=2-1$) and SO ($N_{j}=5_{6}-4_{5}$) moment maps, we find that seven sources show rotation about the jet axis at a few hundred au scales: SerpS-MM18, L1448-C, L1448-NB, L1527, NGC1333- IRAS2A, NGC1333-IRAS4B, and SVS13-B. We analyzed the kinematics of these sources in the $uv$ plane to derive the rotation profiles down to 50~au scales. We find evidence for Keplerian rotation in only two sources, L1527 and L1448-C. Overall, this suggests that Keplerian disks larger than 50 au are uncommon around Class 0 protostars. However, in some of the sources, the line emission could be optically thick and dominated by the envelope emission. Due to the optical thickness of these envelopes, some of the disks could have remained undetected in our observations.

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.

Accurate measurements of the physical structure of protoplanetary discs are critical inputs for planet formation models. These constraints are traditionally established via complex modelling of continuum and line observations. Instead, we present an empirical framework to locate the CO isotopologue emitting surfaces from high spectral and spatial resolution ALMA observations. We apply this framework to the disc surrounding IM Lupi, where we report the first direct, i.e. model independent, measurements of the radial and vertical gradients of temperature and velocity in a protoplanetary disc. The measured disc structure is consistent with an irradiated self-similar disc structure, where the temperature increases and the velocity decreases towards the disc surface. We also directly map the vertical CO snow line, which is located at about one gas scale height at radii between 150 and 300 au, with a CO freeze-out temperature of $21\pm2$ K. In the outer disc ($> 300$ au), where the gas surface density transitions from a power law to an exponential taper, the velocity rotation field becomes significantly sub-Keplerian, in agreement with the expected steeper pressure gradient. The sub-Keplerian velocities should result in a very efficient inward migration of large dust grains, explaining the lack of millimetre continuum emission outside of 300 au. The sub-Keplerian motions may also be the signature of the base of an externally irradiated photo-evaporative wind. In the same outer region, the measured CO temperature above the snow line decreases to $\approx$ 15 K because of the reduced gas density, which can result in a lower CO freeze-out temperature, photo-desorption, or deviations from local thermodynamic equilibrium.

Aims. We wish to clarify the origin of the multiple jet features emanating from the binary protostar SVS 13A (= VLA4A/VLA4B). Methods. We used the Plateau de Bure Interferometer to map at 0.3-0.8" (~70-190 au) dust emission at 1.4 mm, CO(2-1), SiO(5-4), SO(65-54). Revised proper motions for VLA4A/4B and jet wiggling models are computed to clarify their respective contribution. Results. VLA4A shows compact dust emission suggestive of a disk < 50 au, and is the hot corino source, while CO/SiO/SO counterparts to the small-scale H2 jet originate from VLA4B and reveal the jet variable velocity structure. This jet exhibits ~ 3" wiggling consistent with orbital motion around a yet undetected ~ 20-30 au companion to VLA4B, or jet precession. Jet wiggling combined with velocity variability can explain the large apparent angular momentum in CO bullets. We also uncover a synchronicity between CO jet bullets and knots in the HH7-11 chain demonstrating that they trace two distinct jets. Their ~ 300 yr twin outburst period may be triggered by close perihelion approach of VLA4A in an eccentric orbit around VLA4B. A third jet is tentatively seen at PA ~ 0 degrees. Conclusions. SVS13 A harbors at least 2 and possibly 3 distinct jet sources. The CO and HH7-11 jets are launched from quasi-coplanar disks, separated by 20-70 au. Their synchronous major events every 300 yr favor external triggering by close binary interactions, a scenario also invoked for FU Or outbursts.

(Abridged) Through spectrally unresolved observations of high-J CO transitions, Herschel-PACS has revealed large reservoirs of warm (300 K) and hot (700 K) molecular gas around low-mass protostars. We aim to shed light on the excitation and origin of the CO ladder observed toward protostars, and on the water abundance in different physical components using spectrally resolved Herschel-HIFI data. Observations are presented of the highly excited CO line J=16-15 with Herschel-HIFI toward 24 low-mass protostellar objects. The spectrally resolved profiles show two distinct velocity components: a broad component with an average FWHM of 20 km/s, and a narrower component with a FWHM of 5 km/s that is often offset from the source velocity. The average rotational temperature over the entire profile, as measured from comparison between CO J=16-15 and 10-9 emission, is ~300 K. A radiative-transfer analysis shows that the average H2O/CO column-density ratio is ~0.02, suggesting a total H2O abundance of ~2x10^-6. Two distinct velocity profiles observed in the HIFI line profiles suggest that the CO ladder observed with PACS consists of two excitation components. The warm component (300 K) is associated with the broad HIFI component, and the hot component (700 K) is associated with the offset HIFI component. The former originates in either outflow cavity shocks or the disk wind, and the latter in irradiated shocks. The ubiquity of the warm and hot CO components suggests that fundamental mechanisms govern the excitation of these components; we hypothesize that the warm component arises when H2 stops being the dominant coolant. In this scenario, the hot component arises in cooling molecular H2-poor gas just prior to the onset of H2 formation. High spectral resolution observations of highly excited CO transitions uniquely shed light on the origin of warm and hot gas in low-mass protostellar objects.

DH Tau is a young ($\sim$1 Myr) classical T Tauri star. It is one of the few young PMS stars known to be associated with a planetary mass companion, DH Tau b, orbiting at large separation and detected by direct imaging. DH Tau b is thought to be accreting based on copious H${\alpha}$ emission and exhibits variable Paschen Beta emission. NOEMA observations at 230 GHz allow us to place constraints on the disk dust mass for both DH Tau b and the primary in a regime where the disks will appear optically thin. We estimate a disk dust mass for the primary, DH Tau A of $17.2\pm1.7\,M_{\oplus}$, which gives a disk-to-star mass ratio of 0.014 (assuming the usual Gas-to-Dust mass ratio of 100 in the disk). We find a conservative disk dust mass upper limit of 0.42$M_{\oplus}$ for DH Tau b, assuming that the disk temperature is dominated by irradiation from DH Tau b itself. Given the environment of the circumplanetary disk, variable illumination from the primary or the equilibrium temperature of the surrounding cloud would lead to even lower disk mass estimates. A MCFOST radiative transfer model including heating of the circumplanetary disk by DH Tau b and DH Tau A suggests that a mass averaged disk temperature of 22 K is more realistic, resulting in a dust disk mass upper limit of 0.09$M_{\oplus}$ for DH Tau b. We place DH Tau b in context with similar objects and discuss the consequences for planet formation models.

Aims: In this paper we focus on the occurrence of glycolaldehyde (HCOCH2OH) in young solar analogs by performing the first homogeneous and unbiased study of this molecule in the Class 0 protostars of the nearby Perseus star forming region. Methods: We obtained sub-arcsec angular resolution maps at 1.3mm and 1.4mm of glycolaldehyde emission lines using the IRAM Plateau de Bure (PdB) interferometer in the framework of the CALYPSO IRAM large program. Results: Glycolaldehyde has been detected towards 3 Class 0 and 1 Class I protostars out of the 13 continuum sources targeted in Perseus: NGC1333-IRAS2A1, NGC1333-IRAS4A2, NGC1333-IRAS4B1, and SVS13-A. The NGC1333 star forming region looks particularly glycolaldehyde rich, with a rate of occurrence up to 60%. The glycolaldehyde spatial distribution overlaps with the continuum one, tracing the inner 100 au around the protostar. A large number of lines (up to 18), with upper-level energies Eu from 37 K up to 375 K has been detected. We derived column densities > 10^15 cm^-2 and rotational temperatures Trot between 115 K and 236 K, imaging for the first time hot-corinos around NGC1333-IRAS4B1 and SVS13-A. Conclusions: In multiple systems glycolaldehyde emission is detected only in one component. The case of the SVS13-A+B and IRAS4-A1+A2 systems support that the detection of glycolaldehyde (at least in the present Perseus sample) indicates older protostars (i.e. SVS13-A and IRAS4-A2), evolved enough to develop the hot-corino region (i.e. 100 K in the inner 100 au). However, only two systems do not allow us to firmly conclude whether the primary factor leading to the detection of glycolaldehyde emission is the environments hosting the protostars, evolution (e.g. low value of Lsubmm/Lint), or accretion luminosity (high Lint).

Fast jets are thought to be a crucial ingredient of star formation because they might extract angular momentum from the disk and thus allow mass accretion onto the star. However, it is unclear whether jets are ubiquitous, and likewise, their contribution to mass and angular momentum extraction during protostar formation remains an open question. Our aim is to investigate the ejection process in the low-mass Class 0 protostar L1157. This source is associated with a spectacular bipolar outflow, and the recent detection of high-velocity SiO suggests the occurrence of a jet. Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The jet and outflow structure were fit with a precession model. We derived the column density of CO and SiO, as well as the jet mass-loss rate and mechanical luminosity. High-velocity CO and SiO emission resolve for the first time the first 200 au of the outflow-driving molecular jet. The jet is strongly asymmetric, with the blue lobe 0.65 times slower than the red lobe. This suggests that the large-scale asymmetry of the outflow is directly linked to the jet velocity and that the asymmetry in the launching mechanism has been at work for the past 1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that the jet formation mechanism from Class 0 to Class II stages might be similar. Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow by assuming that the jet precesses counter-clockwise on a cone inclined by 73 degree to the line of sight with an opening angle of 8 degree on a period of 1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7 Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of the gravitational energy released by the forming star.

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)

Owing to the paucity of sub-arcsecond (sub)mm observations required to probe the innermost regions of newly forming protostars, several fundamental questions are still being debated, such as the existence and coevality of close multiple systems. We study the physical and chemical properties of the jets and protostellar sources in the NGC1333-IRAS4A proto-binary system using continuum emission and molecular tracers of shocked gas. We observed NGC1333-IRAS4A in the SiO(6-5), SO(6_5-5_4), and CO(2-1) lines and the continuum emission at 1.3, 1.4, and 3 mm using the IRAM Plateau de Bure Interferometer in the framework of the CALYPSO large program. We clearly disentangle for the first time the outflow emission from the two sources A1 and A2. The two protostellar jets have very different properties: the A1 jet is faster, has a short dynamical timescale (<10^3 yr), and is associated with H2 shocked emission, whereas the A2 jet, which dominates the large-scale emission, is associated with diffuse emission, bends, and emits at slower velocities. The observed bending of the A2 jet is consistent with the change of propagation direction observed at large scale and suggests jet precession on very short timescales (~200-600 yr). In addition, a chemically rich spectrum with emission from several COMs (e.g. HCOOH, CH3OCHO, CH3OCH3) is only detected towards A2. Finally, very high-velocity shocked emission (~50 km s^-1) is observed along the A1 jet. An LTE analysis shows that SiO, SO, and H2CO abundances in the gas phase are enhanced up to (3-4)x10^{-7}, (1.4-1.7)x10^{-6}, and (3-7.9)x10^{-7}, respectively. The intrinsic different properties of the jets and driving sources in NGC1333-IRAS4A suggest different evolutionary stages for the two protostars, with A1 being younger than A2, in a very early stage of star formation previous to the hot-corino phase.

The structure and kinematics of Class 0 protostars on scales of a few hundred AU is poorly known. Recent observations have revealed the presence of Keplerian disks with a diameter of 150-180 AU in L1527-IRS and VLA1623A, but it is not clear if such disks are common in Class 0 protostars. Here we present high-angular-resolution observations of two methanol lines in NGC1333-IRAS2A. We argue that these lines probe the inner envelope, and we use them to study the kinematics of this region. Our observations suggest the presence of a marginal velocity gradient normal to the direction of the outflow. However, the position velocity diagrams along the gradient direction appear inconsistent with a Keplerian disk. Instead, we suggest that the emission originates from the infalling and perhaps slowly rotating envelope, around a central protostar of 0.1-0.2 M$_\odot$. If a disk is present, it is smaller than the disk of L1527-IRS, perhaps suggesting that NGC1333-IRAS2A is younger.

We investigate the origin of complex organic molecules (COMs) in the gas phase around the low-mass Class~0 protostar NGC1333-IRAS2A, to determine if the COM emission lines trace an embedded disk, shocks from the protostellar jet, or the warm inner parts of the protostellar envelope. In the framework of the CALYPSO (Continuum And Lines in Young ProtoStellar Objects) IRAM Plateau de Bure survey, we obtained large bandwidth spectra at sub-arcsecond resolution towards NGC 1333-IRAS2A. We identify the emission lines towards the central protostar and perform Gaussian fits to constrain the size of the emitting region for each of these lines, tracing various physical conditions and scales. The emission of numerous COMs such as methanol, ethylene glycol, and methyl formate is spatially resolved by our observations. This allows us to measure, for the first time, the size of the COM emission inside the protostellar envelope, finding that it originates from a region of radius 40-100 AU, centered on the NGC 1333-IRAS2A protostellar object. Our analysis shows no preferential elongation of the COM emission along the jet axis, and therefore does not support the hypothesis that COM emission arises from shocked envelope material at the base of the jet. Down to similar sizes, the dust continuum emission is well reproduced with a single envelope model, and therefore does not favor the hypothesis that COM emission arises from the thermal sublimation of grains embedded in a circumstellar disk. Finally, the typical scale $\sim$60 AU observed for COM emission is consistent with the size of the inner envelope where $T_{\rm{dust}} > 100$ K is expected. Our data therefore strongly suggest that the COM emission traces the hot corino in IRAS2A, i.e., the warm inner envelope material where the icy mantles of dust grains evaporate because they are passively heated by the central protostellar object.

Context: The earliest evolutionary stages of low-mass protostars are characterised by hot and fast jets which remove angular momentum from the circumstellar disk, thus allowing mass accretion onto the central object. However, the launch mechanism is still being debated. Aims: We would like to exploit high-angular (~ 0.8") resolution and high-sensitivity images to investigate the origin of protostellar jets using typical molecular tracers of shocked regions, such as SiO and SO. Methods: We mapped the inner 22" of the NGC1333-IRAS2A protostar in SiO(5-4), SO(65-54), and the continuum emission at 1.4 mm using the IRAM Plateau de Bure interferometer in the framework of the CALYPSO IRAM large program. Results: For the first time, we disentangle the NGC1333-IRAS2A Class 0 object into a proto-binary system revealing two protostars (MM1, MM2) separated by ~ 560 AU, each of them driving their own jet, while past work considered a single protostar with a quadrupolar outflow. We reveal (i) a clumpy, fast (up to |V-VLSR| > 50 km/s), and blueshifted jet emerging from the brightest MM1 source, and (ii) a slower redshifted jet, driven by MM2. Silicon monoxide emission is a powerful tracer of high-excitation (Tkin > 100 K; n(H2) > 10^5 cm-3) jets close to the launching region. At the highest velocities, SO appears to mimic SiO tracing the jets, whereas at velocities close to the systemic one, SO is dominated by extended emission, tracing the cavity opened by the jet. Conclusions: Both jets are intrinsically monopolar, and intermittent in time. The dynamical time of the SiO clumps is < 30-90 yr, indicating that one-sided ejections from protostars can take place on these timescales.

The aim of the present work is to perform a comprehensive analysis of the interstellar chemistry of nitrogen, focussing on the gas-phase formation of the smallest polyatomic species and in particular nitrogen hydrides. We present a new chemical network in which the kinetic rates of critical reactions have been updated based on recent experimental and theoretical studies, including nuclear spin branching ratios. Our network thus treats the different spin symmetries of the nitrogen hydrides self-consistently together with the ortho and para forms of molecular hydrogen. This new network is used to model the time evolution of the chemical abundances in dark cloud conditions. The steady-state results are analysed, with special emphasis on the influence of the overall amounts of carbon, oxygen, and sulphur. Our calculations are also compared with Herschel/HIFI observations of NH, NH$_2$, and NH$_3$ detected towards the external envelope of the protostar IRAS 16293-2422. The observed abundances and abundance ratios are reproduced for a C/O gas-phase elemental abundance ratio of $\sim0.8$, provided that the sulphur abundance is depleted by a factor larger than 2. The ortho-to-para ratio of H$_2$ in these models is $\sim10^{-3}$. Our models also provide predictions for the ortho-to-para ratios of NH$_2$ and NH$_3$ of $\sim2.3$ and $\sim0.7$ respectively. We conclude that the abundances of nitrogen hydrides in dark cloud conditions are consistent with the gas-phase synthesis predicted with our new chemical network.

Prestellar cores exhibit a strong chemical differentiation, which is mainly caused by the freeze-out of molecules onto the grain surfaces. Understanding this chemical structure is important, because molecular lines are often used as probes to constrain the core physical properties. Here we present new observations and analysis of the C18O (1-0) and H13CO+ (1-0) line emission in the L1498 and L1517B prestellar cores, located in the Taurus-Auriga molecular complex. We model these observations with a detailed chemistry network coupled to a radiative transfer code. Our model successfully reproduces the observed C18O (1-0) emission for a chemical age of a few 10^5 years. On the other hand, the observed H13CO+ (1-0) is reproduced only if cosmic-ray desorption by secondary photons is included, and if the grains have grown to a bigger size than average ISM grains in the core interior. This grain growth is consistent with the infrared scattered light ("coreshine") detected in these two objects, and is found to increase the CO abundance in the core interior by about a factor four. According to our model, CO is depleted by about 2-3 orders of magnitude in the core center.

We present a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.

We present far-IR (57-196 mu) spectra of 21 protostars in the Orion molecular clouds, obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space observatory, as part of the Herschel Orion Protostar Survey (HOPS) program. We analyzed the CO emission lines (J_up = 14-46) in the PACS spectra, extracted within a projected distance of <= 2000 AU centered on the protostar. The total luminosity of the CO lines observed with PACS (L(CO)) is found to increase with increasing L_bol. The CO rotational temperature implied by the line ratios increases with J, and at least 3-4 rotational temperature components are required to fit the observed rotational diagram. The rotational temperature components are remarkably invariant between protostars and show no dependence on L_bol, T_bol or envelope density, implying that if the emitting gas is in LTE, the CO emission must arise in multiple temperature components that remain independent of L_bol over two orders of magnitudes. The observed CO emission can also be modeled as arising from a single temperature gas component or from a medium with a power-law temperature distribution; both of these require sub-thermally excited molecular gas at low densities (n(H_2) <= 10^6 cm^-3) and high temperatures (T >= 2000 K). Our results suggest that the contribution from PDRs along the envelope cavity walls is unlikely to be the dominant component of the CO emission observed with PACS. Instead, the "universality" of the rotational temperatures and the observed correlation between L(CO) and L_bol can most easily be explained if the observed CO emission originates in shock-heated, hot (T >= 2000 K), sub-thermally excited (n(H_2) <= 10^6 cm^-3) molecular gas. Post-shock gas at these densities is more likely to be found within the outflow cavities along the molecular outflow or along the cavity walls at radii >= several 100-1000 AU.

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

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.

We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $\Delta$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $\Delta$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

GNU Data Language (GDL) is an open-source interpreted language aimed at numerical data analysis and visualisation. It is a free implementation of the Interactive Data Language (IDL) widely used in Astronomy. GDL has a full syntax compatibility with IDL, and includes a large set of library routines targeting advanced matrix manipulation, plotting, time-series and image analysis, mapping, and data input/output including numerous scientific data formats. We will present the current status of the project, the key accomplishments, and the weaknesses - areas where contributions are welcome !

The advent of large instantaneous bandwidth receivers and high spectral resolution spectrometers on (sub-)millimeter telescopes has opened up the possibilities for unbiased spectral surveys. Because of the large amount of data they contain, any analysis of these surveys requires dedicated software tools. Here we present an extension of the widely used CLASS software that we developed to that purpose. This extension, named Weeds, allows for searches in atomic and molecular lines databases (e.g. JPL or CDMS) that may be accessed over the internet using a virtual observatory (VO) compliant protocol. The package permits a quick navigation across a spectral survey to search for lines of a given species. Weeds is also capable of modeling a spectrum, as often needed for line identification. We expect that Weeds will be useful for analyzing and interpreting the spectral surveys that will be done with the HIFI instrument on board Herschel, but also observations carried-out with ground based millimeter and sub-millimeter telescopes and interferometers, such as IRAM-30m and Plateau de Bure, CARMA, SMA, eVLA, and ALMA.

We present the first Spitzer-IRS spectral maps of the Herbig-Haro flow GGD 37 detected in lines of [Ne III], [O IV], [Ar III], and [Ne V]. The detection of extended [O IV] (55 eV) and some extended emission in [Ne V] (97 eV) indicates a shock temperature in excess of 100,000 K, in agreement with X-ray observations, and a shock speed in excess of 200 km s-1. The presence of an extended pho- toionization or collisional ionization region indicates that GGD 37 is a highly unusual protostellar outflow.

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

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.

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in 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 the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

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 present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores.

We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust.

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.

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

We are carrying out a physical and chemical study of the protostellar envelopes in a representative sample of IM Class 0 protostars. In our first paper (Crimier et al. 2010), we determined the physical structure (density-temperature radial profiles) of the protostellar envelopes. Here, we study the CO depletion and N2H+ deuteration. We observed the millimeter lines of C18O, C17O, N2H+ and N2D+ toward the protostars using the IRAM 30m telescope. Based on these observations, we derived the C18O, N2H+ and N2D+ radial abundance profiles across their envelopes using a radiative transfer code. In addition, we modeled the chemistry of the protostellar envelopes. All the C18O 1-0 maps are well fit assuming that the C18O abundance decreases inwards within the protostellar envelope until the gas and dust reach the CO evaporation temperature, 20-25K, where the CO is released back to the gas phase. The N2H+ deuterium fractionation in Class 0 IMs is [N2D+]/[N2H+]=0.005-0.014, two orders of magnitude higher than the elemental [D/H] value in the interstellar medium, but a factor of 10 lower than in pre-stellar clumps. Chemical models account for the C18O and N2H+ observations if we assume the CO abundance is 2 times lower than the canonical value in the inner envelope. This could be the consequence of the CO being converted into CH3OH on the grain surfaces prior to the evaporation and/or the photodissociation of CO by the stellar UV radiation. The deuterium fractionation is not fitted by chemical models. This discrepancy is very likely caused by the simplicity of our model that assumes spherical geometry and neglects important phenomena like the effect of bipolar outflows and UV radiation from the star. More important, the deuterium fractionation is dependent on the ortho-to-para H2 ratio, which is not likely to reach the steady-state value in the dynamical time scales of these protostars.

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.

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

The NGC 1999 reflection nebula features a dark patch with a size of ~10,000 AU, which has been interpreted as a small, dense foreground globule and possible site of imminent star formation. We present Herschel PACS far-infrared 70 and 160mum maps, which reveal a flux deficit at the location of the globule. We estimate the globule mass needed to produce such an absorption feature to be a few tenths to a few Msun. Inspired by this Herschel observation, we obtained APEX LABOCA and SABOCA submillimeter continuum maps, and Magellan PANIC near-infrared images of the region. We do not detect a submillimer source at the location of the Herschel flux decrement; furthermore our observations place an upper limit on the mass of the globule of ~2.4x10^-2 Msun. Indeed, the submillimeter maps appear to show a flux depression as well. Furthermore, the near-infrared images detect faint background stars that are less affected by extinction inside the dark patch than in its surroundings. We suggest that the dark patch is in fact a hole or cavity in the material producing the NGC 1999 reflection nebula, excavated by protostellar jets from the V 380 Ori multiple system.

We present 70 and 160 micron Herschel science demonstration images of a field in the Orion A molecular cloud that contains the prototypical Herbig-Haro objects HH 1 and 2, obtained with the Photodetector Array Camera and Spectrometer (PACS). These observations demonstrate Herschel's unprecedented ability to study the rich population of protostars in the Orion molecular clouds at the wavelengths where they emit most of their luminosity. The four protostars previously identified by Spitzer 3.6-40 micron imaging and spectroscopy are detected in the 70 micron band, and three are clearly detected at 160 microns. We measure photometry of the protostars in the PACS bands and assemble their spectral energy distributions (SEDs) from 1 to 870 microns with these data, Spitzer spectra and photometry, 2MASS data, and APEX sub-mm data. The SEDs are fit to models generated with radiative transfer codes. From these fits we can constrain the fundamental properties of the protostars. We find luminosities in the range 12-84 L_sun and envelope densities spanning over two orders of magnitude. This implies that the four protostars have a wide range of envelope infall rates and evolutionary states: two have dense, infalling envelopes, while the other two have only residual envelopes. We also show the highly irregular and filamentary structure of the cold dust and gas surrounding the protostars as traced at 160 microns.

Context: Intermediate mass protostars provide a bridge between low- and high-mass protostars. Furthermore, they are an important component of the UV interstellar radiation field. Despite their relevance, little is known about their formation process. Aims: We present a systematic study of the physical structure of five intermediate mass, candidate Class 0 protostars. Our two goals are to shed light on the first phase of intermediate mass star formation and to compare these protostars with low- and high-mass sources. Methods: We derived the dust and gas temperature and density profiles of the sample. We analysed all existing continuum data on each source and modelled the resulting SED with the 1D radiative transfer code DUSTY. The gas temperature was then predicted by means of a modified version of the code CHT96. Results: We found that the density profiles of five out of six studied intermediate mass envelopes are consistent with the predictions of the "inside-out" collapse theory.We compared several physical parameters, like the power law index of the density profile, the size, the mass, the average density, the density at 1000 AU and the density at 10 K of the envelopes of low-, intermediate, and high-mass protostars. When considering these various physical parameters, the transition between the three groups appears smooth, suggesting that the formation processes and triggers do not substantially differ.

Context: The low mass protostar IRAS16293-2422 is a prototype Class 0 source with respect to the studies of the chemical structure during the initial phases of life of Solar type stars. Aims: In order to derive an accurate chemical structure, a precise determination of the source physical structure is required. The scope of the present work is the derivation of the structure of IRAS16293-2422. Methods: We have re-analyzed all available continuum data (single dish and interferometric, from millimeter to MIR) to derive accurate density and dust temperature profiles. Using ISO observations of water, we have also reconstructed the gas temperature profile. Results: Our analysis shows that the envelope surrounding IRAS16293-2422 is well described by the Shu "inside-out" collapsing envelope model or a single power-law density profile with index equal to 1.8. In contrast to some previous studies, our analysis does not show evidence of a large (>/- 800 AU in diameter) cavity. Conclusions: Although IRAS16293-2422 is a multiple system composed by two or three objects, our reconstruction will be useful to derive the chemical structure of the large cold envelope surrounding these objects and the warm component, treated here as a single source, from single-dish observations of molecular emission.

We report the results of spectroscopic mapping observations carried out toward protostellar outflows in the BHR71, L1157, L1448, NGC 2071, and VLA 1623 molecular regions using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, provide detailed maps of the 8 lowest pure rotational lines of molecular hydrogen and of the [SI] 25.25 micron and [FeII] 26.0 micron fine structure lines. The molecular hydrogen lines, believed to account for a large fraction of the radiative cooling from warm molecular gas that has been heated by a non-dissociative shock, allow the energetics of the outflows to be elucidated. Within the regions mapped towards these 5 outflow sources, total H2 luminosities ranging from 0.02 to 0.75 L(solar) were inferred for the sum of the 8 lowest pure rotational transitions. By contrast, the much weaker [FeII] 26.0 micron fine structure transition traces faster, dissociative shocks; here, only a small fraction of the fast shock luminosity emerges as line radiation that can be detected with Spitzer/IRS.

Ammonia inversion lines are often used as probes of the physical conditions in the dense ISM. The excitation temperature between the first two para metastable (rotational) levels is an excellent probe of the gas kinetic temperature. However, the calibration of this ammonia thermometer depends on the accuracy of the collisional rates with H2. Here we present new collisional rates for ortho-NH3 and para-NH3 colliding with para-H2 (J=0) and we investigate the effects of these new rates on the excitation of ammonia. Scattering calculations employ a new, high accuracy, potential energy surface computed at the coupled-cluster CCSD(T) level with a basis set extrapolation procedure. Rates are obtained for all transitions involving ammonia levels with J <= 3 and for kinetic temperatures in the range 5-100 K. We find that the calibration curve of the ammonia thermometer -- which relates the observed excitation temperature between the first two para metastable levels to the gas kinetic temperature -- does not change significantly when these new rates are used. Thus, the calibration of ammonia thermometer appears to be robust. Effects of the new rates on the excitation temperature of inversion and rotation-inversion transitions are also found to be small.

We present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities and the predictions of detailed shock models indicates that the emission arises in both slow (12 - 24 km/s) and fast (36 - 53 km/s) C-type shocks with an initial ortho-to-para ratio of ~ 1. The present H2 ortho-to-para ratio exhibits a large degree of spatial variations. In the post-shocked gas, it is usually about 2, i.e. close to the equilibrium value (~ 3). However, around at least two outflows, we observe a region with a much lower (~ 0.5) ortho-to-para ratio. This region probably corresponds to gas which has been heated-up recently by the passage of a shock front, but whose ortho-to-para has not reached equilibrium yet. This, together with the low initial ortho-to-para ratio needed to reproduce the observed emission, provide strong evidence that H2 is mostly in para form in cold molecular clouds. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and thus can be used to ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for core disruption.

We present observations of 13CO, C18O, HCO+, H13CO+, DCO+ and N2H+ line emission towards the Barnard 68 starless core. The line profiles are interpreted using a chemical network coupled with a radiative transfer code in order to reconstruct the radial velocity profile of the core. Our observations and modeling indicate the presence of complex radial motions, with the inward motions in the outer layers of the core but outward motions in the inner part, suggesting radial oscillations. The presence of such oscillation would imply that B68 is relatively old, typically one order of magnitude older than the age inferred from its chemical evolution and statistical core lifetimes. Our study demonstrates that chemistry can be used as a tool to constrain the radial velocity profiles of starless cores.

We present a detailed study of the ionization fraction of the Barnard 68 pre-stellar core, using millimeter H13CO+ and DCO+ lines observations. These observations are compared to the predictions of a radiative transfer model coupled to a chemical network that includes depletion on grains and gas phase deuterium fractionation. Together with previous observations and modelling of CO and isotopologues, our H13CO+ and DCO+ observations and modelling allow to place constraints on the metal abundance and the cosmic ionization rate. The H13CO+ emission is well reproduced for metals abundances lower than 3e-9 and a standard cosmic ray ionization rate. However, the observations are also consistent with a complete depletion of metals, i.e. with cosmic rays as the only source of ionization at visual extinctions greater than a few Av . The DCO+ emission is found to be dependent of the ortho to para H2 ratio, and indicates a ratio of ~1e-2. The derived ionization fraction is about 5e-9 with respect to H nuclei, which is about an order of magnitude lower than the one observed in the L1544 core. The corresponding ambipolar diffusion timescale is found to be an order of magnitude larger than the free fall timescale at the center of the core. The inferred metal abundance suggests that magnetically inactive regions (dead zones) are present in protostellar disks.

We present a direct comparison of a chemical/physical model to multitransitional observations of C18O and 13CO towards the Barnard 68 pre-stellar core. These observations provide a sensitive test for models of low UV field photodissociation regions and offer the best constraint on the gas temperature of a pre-stellar core. We find that the gas temperature of this object is surprisingly low (~7-8 K), and significantly below the dust temperature, in the outer layers (Av < 5 mag) that are traced by C18O and 13CO emission. As shown previously, the inner layers (Av > 5 mag) exhibit significant freeze-out of CO onto grain surfaces. Because the dust and gas are not fully coupled, depletion of key coolants in the densest layers raises the core (gas) temperature, but only by ~1 K. The gas temperature in layers not traced by C18O and 13CO emission can be probed by NH3 emission, with a previously estimated temperature of ~10-11 K. To reach these temperatures in the inner core requires an order of magnitude reduction in the gas to dust coupling rate. This potentially argues for a lack of small grains in the densest gas, presumably due to grain coagulation.

Context : Despite the low cosmic abundance of deuterium (D/H ~ 1e-5), large degrees of deuterium fractionation in molecules are observed in star forming regions with enhancements that can reach 13 orders of magnitude, which current models have difficulties to account for. Aims : Multi-isotopologue observations are a very powerful constraint for chemical models. The aim of our observations is to understand the processes forming the observed large abundances of methanol and formaldehyde in low-mass protostellar envelopes (gas-phase processes ? chemistry on the grain surfaces ?) and better constrain the chemical models. Methods : Using the IRAM 30m single-dish telescope, we observed deuterated formaldehyde (HDCO and D2CO) and methanol (CH2DOH, CH3OD, and CHD2OH) towards a sample of seven low-mass class 0 protostars. Using population diagrams, we then derive the fractionation ratios of these species (abundance ratio between the deuterated molecule and its main isotopologue) and compare them to the predictions of grain chemistry models. Results : These protostars show a similar level of deuteration as in IRAS16293-2422, where doubly-deuterated methanol -- and even triply-deuterated methanol -- were first detected. Our observations point to the formation of methanol on the grain surfaces, while formaldehyde formation cannot be fully pined down. While none of the scenarii can be excluded (gas-phase or grain chemistry formation), they both seem to require abstraction reactions to reproduce the observed fractionations.

We present observations of methanol lines in a sample of Class 0 low mass protostars. Using a 1-D radiative transfer model, we derive the abundances in the envelopes. In two sources of the sample, the observations can only be reproduced by the model if the methanol abundance is enhanced by about two order of magnitude in the inner hot region of the envelope. Two other sources show similar jumps, although at a lower confidence level. The observations for the other three sources are well reproduced with a constant abundance, but the presence of a jump cannot be ruled out. The observed methanol abundances in the warm gas around low mass protostars are orders of magnitude higher than gas phase chemistry models predict. Hence, in agreement with other evidences, this suggest that the high methanol abundance reflects recent evaporation of ices due to the heating by the newly formed star. The observed abundance ratios of CH3 OH, H2 CO, and CO are in good agreement with grain surface chemistry models. However, the absolute abundances are more difficult to reproduce and may point towards the presence of multiple ice components in these regions.

We present theoretical predictions of the rotational line emission of deuterated water in low-mass protostar collapsing envelopes. The model accounts for the density and temperature structure of the envelope, according the inside-out collapse framework. The deuterated water abundance profile is approximated by a step function, with a low value in the cold outer envelope and a higher value in the inner envelope where the grain mantles evaporate. The two abundances are the two main parameters of the modeling, along with the temperature at which the mantles evaporate. We report line flux predictions for a 30 and 5 L$_\odot$ source luminosity respectively. We show that ground based observations are capable to constrain the three parameters of the model in the case of bright low-mass protostars (L$>$10 L$_{\odot}$), and that no space based observations, like for example HSO observations, are required in this case. On the contrary, we show that the study of low-luminosity sources (L$<$10 L$_{\odot}$), assuming the same HDO abundance profile, requires too much integration time to be carried out either with available ground-based telescopes or with the HIFI instrument on board HSO. For these sources, only the large interferometer ALMA will allow to constrain the HDO abundance.

The results from a single-dish molecular line survey of a set of 18 deeply embedded young stellar objects are summarized. More than 40 lines from 16 different species were observed with the JCMT, Onsala, IRAM 30m and SEST telescopes. The multi-transition data are analyzed using a temperature and density structure derived from models of the dust continuum emission. For the outer envelope (>300 AU), the data indicate a `drop' abundance profile for many species, with normal abundances in the outer- and innermost regions and highly depleted abundances in an intermediate zone. This zone is bounded at the outer than the lifetime of the core, and at the inner edge by the evaporation temperature of the species involved. In the innermost envelope (<300 AU), all ices evaporate resulting in jumps in the abundances of complex organic molecules such as CH3OH. A key project for Herschel will be to survey gas-phase water in these objects, whose abundance shows extreme variations with temperature. ALMA wil be able to directly image the chemical variations throughout the envelope and zoom in on the inner hot core and protoplanetary disks on scales on tens of AU.

Methanol lines observations of a sample of low mass Class 0 protostars are presented. Using a 1D radiative transfer model, I show that several protostars have large abundance jumps in the inner hot and dense region of envelopes, probably because of thermal grain mantle evaporation. These abundances are compared with a grain surface chemistry model.

Complex organic molecules have previously been discovered in solar type protostars, raising the questions of where and how they form in the envelope. Possible formation mechanisms include grain mantle evaporation, interaction of the outflow with its surroundings or the impact of UV/X-rays inside the cavities. In this Letter we present the first interferometric observations of two complex molecules, CH3CN and HCOOCH3, towards the solar type protostar IRAS16293-2422. The images show that the emission originates from two compact regions centered on the two components of the binary system. We discuss how these results favor the grain mantle evaporation scenario and we investigate the implications of these observations for the chemical composition and physical and dynamical state of the two components.

We report the detection of complex molecules (HCOOCH_3, HCOOH and CH_3CN), signposts of a "hot core" like region, toward the low mass, Class 0 source NGC1333-IRAS4A. This is the second low mass protostar where such complex molecules have been searched for and reported, the other source being IRAS16293-2422. It is therefore likely that compact (few tens of AUs) regions of dense and warm gas, where the chemistry is dominated by the evaporation of grain mantles, and where complex molecules are found, are common in low mass Class 0 sources.Given that the chemical formation timescale is much shorter than the gas hot core crossing time, it is not clear whether the reported complex molecules are formed on the grain surfaces (first generation molecules) or in the warm gas by reactions involving the evaporated mantle constituents (second generation molecules). We do not find evidence for large differences in the molecular abundances, normalized to the formaldehyde abundance, between the two solar type protostars, suggesting perhaps a common origin.

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.

We present a survey of the formaldehyde emission in nine class 0 protostars obtained with the IRAM 30m and the JCMT millimeter telescopes. Using a detailed radiative transfer code of the envelopes surrounding the protostars, we show that all but one of the observed objects show an inner warm evaporation region where the formaldehyde is much more abundant (up to three orders of magnitude) than in the outer cold part. The largest inner formaldehyde abundances are associated with the sources having the lowest submillimetric to bolometric luminosity ratio, i.e. with sources closer to the class I border. These abundances are compared with predictions from recent models of hot core chemistry.

We present observations towards one of the closest regions of high mass star formation, Orion BN/KL, performed at both low resolution mode (grating mode) and high resolution mode (Fabry-Perot) with the Long Wavelength Spectrometer on board the Infrared Space Observatory. We detected the CO rotational lines from Jup = 15 to Jup = 45. While the lines with Jup<= 32 are spectrally unresolved, the higher lying lines show a broadened profile. Finally, we detected two 13CO lines, namely at Jup = 18 and 24, from which we could derive the opacities of the relative 12CO lines. The LVG analysis of the observed line spectrum allows to distinguish three main physical components with different temperatures, densities and column densities: 1) lines with Jup< 20 originate mainly in the diffuse photodissociation region surrounding the source; 2) lines with Jup between 20 and 30 originate in the high velocity outflow (plateau) emanating from IrC2; 3) lines with Jup > 32 originate in the hot and dense gas of the shocked component of the outflow. We discuss how future observations with HIFI, onboard the Far Infrared Space Telescope (FIRST) will allow to spectrally and spatially disentangle the three components, and, consequently, characterise more precisely the Orion BN/KL star forming region.

Sun like stars are born from the collapse of fragment of molecular clouds. During the first embedded phase (the so called class 0/1), the protostar is surrounded by a collapsing envelope, whose physical and chemical structure would set up the initial conditions of the proto-stellar disk which may eventually form planets. In this contribution, we show how spectroscopic far infrared to millimeter observations can be a powerful diagnostic to derive the chemical and physical structure of the innermost regions of envelopes collapsing around solar type protostars, where disk are believe to form. Our claim is based on the simultaneous modeling of the dynamics, chemistry and radiative transfer in a protostellar envelope. Here we present a study of a class 0 protostar NGC1333-IRAS4, obtained using ISO-LWS observations of water far-infrared lines. We show the derived physical and chemical structure of the protostar, namely the density and temperature profiles, as well as the water abundance across the envelope. Besides, the mass of the central protostar and its accretion rate can be constrained by these observations and modeling, yielding an alternative method to measure these two poorly know parameters.

Sun like stars are born from the collapse of fragment of molecular clouds. During the first embedded phase (the so called class 0/1), the protostar is surrounded by a collapsing envelope, whose physical and chemical structure would set up the initial conditions of the proto-stellar disk which may eventually form planets. In this contribution, we show how spectroscopic far infrared to millimeter observations can be a powerful diagnostic to derive the chemical and physical structure of the innermost regions of envelopes collapsing around solar type protostars, where disk are believe to form. Our claim is based on the simultaneous modeling of the dynamics, chemistry and radiative transfer in a protostellar envelope. Here we present a study of a class 0 protostar NGC1333-IRAS4, obtained using ISO-LWS observations of water far-infrared lines. We show the derived physical and chemical structure of the protostar, namely the density and temperature profiles, as well as the water abundance across the envelope. Besides, the mass of the central protostar and its accretion rate can be constrained by these observations and modeling, yielding an alternative method to measure these two poorly know parameters.

We report ISO-LWS far infrared observations of CO, water and oxygen lines towards the protobinary system IRAS4 in the NGC1333 cloud. We detected several water, OH, CO rotational lines, and two [OI] and [CII] fine structure lines. Given the relatively poor spectral and spatial resolution of these observations, assessing the origin of the observed emission is not straightforward. In this paper, we focus on the water line emission and explore the hypothesis that it originates in the envelopes that surround the two protostars, IRAS4 A and B, thanks to an accurate model. The model reproduces quite well the observed water line fluxes, predicting a density profile, mass accretion rate, central mass, and water abundance profile in agreement with previous works. We hence conclude that the emission from the envelopes is a viable explanation for the observed water emission, although we cannot totally rule out the alternative that the observed water emission originates in the outflow.