Large scale spectral maps of star forming regions enable the comparative study of the gas excitation around an ensemble of sources at a common frame of reference, providing direct insights in the multitude of processes involved. In this paper we employ spectral-line maps to decipher the excitation, the kinematical and dynamical processes in NGC 1333 as revealed by a number of different emission lines, aiming to set a reference for the applicability of tracers in constraining diverse physical processes. We reconstruct line maps for H$_2$ , CO, H$_2$O and C$^+$ using data obtained with the Spitzer-IRS and Herschel HIFI-SPIRE. We compare the morphological features of the maps and derive the gas excitation for regions of interest employing LTE and non-LTE methods. We also calculate the kinematical and dynamical properties for each outflow tracer consistently for all outflows in NGC 1333. We finally measure the water abundance in outflows with respect to carbon monoxide and molecular hydrogen. CO and H$_2$ are highly excited around B-stars and at lower levels trace protostellar outflows. H$_2$O emission is dominated by a moderately fast component associated with outflows. Intermediate J CO lines appear brightest at the locations traced by a narrow H$_2$O component, indicating that beyond the dominating collisional processes, a secondary, radiative excitation component can also be active. The morphology, kinematics, excitation and abundance variations of water are consistent with its excitation and partial dissociation in shocks. Water abundance ranges between 5 x 10$^{-7}$ and 10$^{-5}$, with the lower values being more representative. Water is brightest and most abundant around IRAS 4A which is consistent with the latter hosting a hot corino source. Finally, the outflow mass flux is found highest for CO and decreases by one and two orders of magnitude for H$_2$ and H$_2$O, respectively.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S $(J=2-1)$ and N$_2$H$^+$ $(J=1-0)$ with the NRO 45-m telescope, doubling the number of cores where the CS depletion was probed using C$^{34}$S. In most of our targets, the distribution of C$^{34}$S emission shows features that suggest that the CS molecule is generally depleted in the center of the prestellar cores. The radial profile of the CS abundance with respect to H$_2$ directly measured from the CS emission and the Herschel dust emission indicates that the CS molecule is depleted by a factor of $\sim$3 toward the central regions of the cores with respect to their outer regions. The degree of the depletion is found to be even more enhanced by an order of magnitude when the contaminating effect introduced by the presence of CS molecules in the surrounding envelope that lie along the line-of-sight is removed. Except for L1197 which is classified as relatively the least evolved core in our targets based on its observed physical parameters, we found that the remaining four prestellar cores are suffering from significant CS depletion at their central region regardless of the relative difference in their evolutionary status.

L1521E is a dense starless core in Taurus that was found to have relatively low molecular depletion by earlier studies, thus suggesting a recent formation. We aim to characterize the chemical structure of L1521E and compare it to the more evolved L1544 pre-stellar core. We have obtained $\sim$2.5$\times$2.5 arcminute maps toward L1521E using the IRAM-30m telescope in transitions of various species. We derived abundances for the species and compared them to those obtained toward L1544. We estimated CO depletion factors. Similarly to L1544, $c$-C$_3$H$_2$ and CH$_3$OH peak at different positions. Most species peak toward the $c$-C$_3$H$_2$ peak. The CO depletion factor derived toward the $Herschel$ dust peak is 4.3$\pm$1.6, which is about a factor of three lower than that toward L1544. The abundances of sulfur-bearing molecules are higher toward L1521E than toward L1544 by factors of $\sim$2-20. The abundance of methanol is similar toward the two cores. The higher abundances of sulfur-bearing species toward L1521E than toward L1544 suggest that significant sulfur depletion takes place during the dynamical evolution of dense cores, from the starless to pre-stellar stage. The CO depletion factor measured toward L1521E suggests that CO is more depleted than previously found. Similar CH$_3$OH abundances between L1521E and L1544 hint that methanol is forming at specific physical conditions in Taurus, characterized by densities of a few $\times$10$^4$ cm$^{-3}$ and $N$(H$_2$)$\gtrsim$10$^{22}$ cm$^{-2}$, when CO starts to catastrophically freeze-out, while water can still be significantly photodissociated, so that the surfaces of dust grains become rich in solid CO and CH$_3$OH, as already found toward L1544. Methanol can thus provide selective crucial information about the transition region between dense cores and the surrounding parent cloud.

We present an IRAM 30m line survey of the dense core L483 in the 3 mm band. We detected 71 molecules (140 including different isotopologs), most of which are present in the cold and quiescent ambient cloud. Of particular interest among the detected molecules are the cis isomer of HCOOH, the complex organic molecules HCOOCH3, CH3OCH3, and C2H5OH, a wide variety of carbon chains, nitrogen oxides like N2O, and saturated molecules like CH3SH. In general, fractional molecular abundances in L483 are systematically lower than in TMC-1 (especially for carbon chains), tend to be higher than in L1544 and B1-b, and are similar to those in L1527. Apart from the overabundance of carbon chains in TMC-1, we find that L483 does not have a marked chemical differentiation with respect to starless/prestellar cores like TMC-1 and L1544, although it does chemically differentiate from Class 0 hot corino sources like IRAS 16293-2422. This fact suggests that the chemical composition of the ambient cloud of some Class 0 sources could be largely inherited from the dark cloud starless/prestellar phase. We also derived isotopic ratios for a variety of molecules, some of which show isotopic anomalies like an extreme depletion in 13C for one of the two isotopologs of c-C3H2, a drastic enrichment in 18O for SO and HNCO (SO being also largely enriched in 17O), and different abundances for the two 13C substituted species of C2H and the two 15N substituted species of N2H+. We report the first detection in space of some minor isotopologs like c-C3D. The exhaustive chemical characterization of L483 presented here, together with similar studies of other prestellar and protostellar sources, should allow us to identify the main factors that regulate the chemical composition of cores along the process of formation of low-mass protostars.

The main carrier of sulfur in dense clouds, where it is depleted from the gas phase, remains a mystery. Shock waves in young molecular outflows disrupt the ice mantles and allow us to directly probe the material that is ejected into the gas phase. A comprehensive study of sulfur-bearing species towards L1157-B1, a shocked region along a protostellar outflow, has been carried out as part of the IRAM-30m large program ASAI. The dataset contains over 100 lines of CCS, H$_2$CS, OCS, SO, SO$_2$ and isotopologues. The results of these observations are presented, complementing previous studies of sulfur-bearing species in the region. The column densities and fractional abundances of these species are measured and together these species account for 10\% of the cosmic sulfur abundance in the region. The gas properties derived from the observations are also presented, demonstrating that sulfur bearing species trace a wide range of different gas conditions in the region.

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.

Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble and fragmentation may take place, with the consequent formation of binaries and multiple systems. We present ALMA Band 6 observations (ACA and 12m array) of the dust continuum emission of the 8 Msun pre-stellar core L1544, with angular resolution of 2'' x 1.6'' (linear resolution 270 au x 216 au). Within the primary beam, a compact region of 0.1 Msun, which we call a "kernel", has been unveiled. The kernel is elongated, with a central flat zone with radius Rker ~ 10'' (~ 1400 au). The average number density within Rker is ~1 x 10^6 cm^{-3}, with possible local density enhancements. The region within Rker appears to have fragmented, but detailed analysis shows that similar substructure can be reproduced by synthetic interferometric observations of a smooth centrally concentrated dense core with a similar central flat zone. The presence of a smooth kernel within a dense core is in agreement with non-ideal magneto-hydro-dynamical simulations of a contracting cloud core with a peak number density of 1 x 10^7 cm^{-3}. Dense cores with lower central densities are completely filtered out when simulated 12m-array observations are carried out. These observations demonstrate that the kernel of dynamically evolved dense cores can be investigated at high angular resolution with ALMA.

The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region.

We study the kinematics of the dense gas of starless and protostellar cores traced by the N2D+(2-1), N2H+(1-0), DCO+(2-1), and H13CO+(1-0) transitions along the L1495 filament and the kinematic links between the cores and the surrounding molecular cloud. We measure velocity dispersions, local and total velocity gradients and estimate the specific angular momenta of 13 dense cores in the four transitions using the on-the-fly observations with the IRAM 30 m antenna. To study a possible connection to the filament gas, we use the fit results of the C18O(1-0) survey performed by Hacar et al. (2013). All cores show similar properties along the 10 pc-long filament. N2D+(2-1) shows the most centrally concentrated structure, followed by N2H+(1-0) and DCO+(2-1), which show similar spatial extent, and H13CO+(1-0). The non-thermal contribution to the velocity dispersion increases from higher to lower density tracers. The change of magnitude and direction of the total velocity gradients depending on the tracer used indicates that internal motions change at different depths within the cloud. N2D+ and N2H+ show smaller gradients than the lower density tracers DCO+ and H13CO+, implying a loss of specific angular momentum at small scales. At the level of cloud-core transition, the core's external envelope traced by DCO+ and H13CO+ is spinning up, consistent with conservation of angular momentum during core contraction. C18O traces the more extended cloud material whose kinematics is not affected by the presence of dense cores. The decrease in specific angular momentum towards the centres of the cores shows the importance of local magnetic fields to the small scale dynamics of the cores. The random distributions of angles between the total velocity gradient and large scale magnetic field suggests that the magnetic fields may become important only in the high density gas within dense cores.

The isocyanate radical (NCO) is the simplest molecule containing the backbone of the peptide bond, C(=O)-N. This bond has a prebiotic interest since is the one linking two amino acids to form large chains of proteins. It is also present in some organic molecules observed in space such as HNCO, NH2CHO and CH3NCO. In this letter we report the first detection in space of NCO towards the dense core L483. We also report the identification of the ion H2NCO+, definitively confirming its presence in space, and observations of HNCO, HOCN, and HCNO in the same source. For NCO, we derive a column density of 2.2e12 cm-2, which means that it is only about 5 times less abundant than HNCO. We find that H2NCO+, HOCN and HCNO have abundances relative to HNCO of 1/400, 1/80, and 1/160, respectively. Both NCO and H2NCO+ are involved in the production of HNCO and several of its isomers. We have updated our previous chemical models involving NCO and the production of the CHNO isomers. Taking into account the uncertainties in the model, the observed abundances are reproduced relatively well. Indeed, the detection of NCO and H2NCO+ in L483 supports the chemical pathways to the formation of the detected CHNO isomers. Sensitive observations of NCO in sources where other molecules containing the C(=O)-N subunit have been detected could help in elucidating its role in prebiotic chemistry in space.

Evidence is mounting that the small bodies of our Solar System, such as comets and asteroids, have at least partially inherited their chemical composition from the first phases of the Solar System formation. It then appears that the molecular complexity of these small bodies is most likely related to the earliest stages of star formation. It is therefore important to characterize and to understand how the chemical evolution changes with solar-type protostellar evolution. We present here the Large Program "Astrochemical Surveys At IRAM" (ASAI). Its goal is to carry out unbiased millimeter line surveys between 80 and 272 GHz of a sample of ten template sources, which fully cover the first stages of the formation process of solar-type stars, from prestellar cores to the late protostellar phase. In this article, we present an overview of the surveys and results obtained from the analysis of the 3 mm band observations. The number of detected main isotopic species barely varies with the evolutionary stage and is found to be very similar to that of massive star-forming regions. The molecular content in O- and C- bearing species allows us to define two chemical classes of envelopes, whose composition is dominated by either a) a rich content in O-rich complex organic molecules, associated with hot corino sources, or b) a rich content in hydrocarbons, typical of Warm Carbon Chain Chemistry sources. Overall, a high chemical richness is found to be present already in the initial phases of solar-type star formation.

We present the first identification in interstellar space of the thioformyl radical (HCS) and its metastable isomer HSC. These species were detected toward the molecular cloud L483 thanks to observations carried out with the IRAM 30m telescope in the 3 mm band. We derive beam-averaged column densities of 7e12 cm-2 for HCS and 1.8e11 cm-2 for HSC, which translate to fractional abundances relative to H2 of 2e-10 and 6e-12, respectively. Although the amount of sulfur locked by these radicals is low, their detection allows to put interesting constraints on the chemistry of sulfur in dark clouds. Interestingly, the H2CS/HCS abundance ratio is found to be quite low, around 1, in contrast with the oxygen analogue case, in which the H2CO/HCO abundance ratio is around 10 in dark clouds. Moreover, the radical HCS is found to be more abundant than its oxygen analogue, HCO. The metastable species HOC, the oxygen analogue of HSC, has not been yet observed in space. These observational constraints are confronted with the outcome of a recent model of the chemistry of sulfur in dark clouds. The model underestimates the fractional abundance of HCS by at least one order of magnitude, overestimates the H2CS/HCS abundance ratio, and does not provide an abundance prediction for the metastable isomer HSC. These observations should prompt a revision of the chemistry of sulfur in interstellar clouds.

Abridged. Are all filaments bundles of fibers? To address this question, we have investigated the gas organization within the paradigmatic Integral Shape Filament (ISF). We combined two new ALMA Cycle 3 mosaics with previous IRAM 30m observations to produce a high-dynamic range N$_2$H$^+$(1-0) emission map of the ISF tracing its high-density material and velocity structure down to scales of 0.009 pc. From the analysis of the gas kinematics, we identify a total of 55 dense fibers in the central region of the ISF. Independently of their location, these fibers are characterized by transonic internal motions, lengths of ~0.15 pc, and masses per-unit-length close to those expected in hydrostatic equilibrium. The ISF fibers are spatially organized forming a dense bundle with multiple hub-like associations likely shaped by the local gravitational potential. Within this complex network, the ISF fibers show a compact radial emission profile with a median FWHM of 0.035 pc systematically narrower than the previously proposed universal 0.1 pc filament width. Our ALMA observations reveal complex bundles of fibers in the ISF, suggesting strong similarities between the internal substructure of this massive filament and previously studied lower-mass objects. The fibers show identical dynamic properties in both low- and high-mass regions, and their widespread detection suggests a preferred organizational mechanism of gas in which the physical fiber dimensions (width and length) are self-regulated depending on their intrinsic gas density. Combined with previous works, we identify a systematic increase of the surface density of fibers as a function of the total mass per-unit-length in filamentary clouds. Based on this empirical correlation, we propose a unified star-formation scenario where the observed differences between low- and high-mass clouds emerge naturally from the initial concentration of fibers.

The high-sensitivity of the IRAM 30-m ASAI unbiased spectral survey in the mm-window allows us to detect NO emission towards both the Class I object SVS13-A and the protostellar outflow shock L1157-B1. We detect the hyperfine components of the $^2\Pi_{\rm 1/2}$ $J$ = 3/2 $\to$ 1/2 (at 151 GHz) and the $^2\Pi_{\rm 1/2}$ $J$ = 5/2 $\to$ 3/2 (250 GHz) spectral pattern. The two objects show different NO profiles: (i) SVS13-A emits through narrow (1.5 km s$^{-1}$) lines at the systemic velocity, while (ii) L1157-B1 shows broad ($\sim$ 5 km s$^{-1}$) blue-shifted emission. For SVS13-A the analysis leads to $T_{\rm ex}$ $\geq$ 4 K, $N(\rm NO)$ $\leq$ 3 $\times$ 10$^{15}$ cm$^{-2}$, and indicates the association of NO with the protostellar envelope. In L1157-B1, NO is tracing the extended outflow cavity: $T_{\rm ex}$ $\simeq$ 4--5 K, and $N(\rm NO)$ = 5.5$\pm$1.5 $\times$ 10$^{15}$ cm$^{-2}$. Using C$^{18}$O, $^{13}$C$^{18}$O, C$^{17}$O, and $^{13}$C$^{17}$O ASAI observations we derive an NO fractional abundance less than $\sim$ 10$^{-7}$ for the SVS13-A envelope, in agreement with previous measurements towards extended PDRs and prestellar objects. Conversely, a definite $X(NO)$ enhancement is measured towards L1157-B1, $\sim$ 6 $\times$ 10$^{-6}$, showing that the NO production increases in shocks. The public code UCLCHEM was used to interpret the NO observations, confirming that the abundance observed in SVS13-A can be attained in an envelope with a gas density of 10$^5$ cm$^{-3}$ and a kinetic temperature of 40 K. The NO abundance in L1157-B1 is reproduced with pre-shock densities of 10$^5$ cm$^{-3}$ subjected to a $\sim$ 45 km s$^{-1}$ shock.

Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the $^{13}$CH$_3$OH and CH$_2$DOH distribution in the inner regions ($\sim$100 au) of the HH212 system in Orion B. To this end, we used ALMA Cycle 1 and Cycle 4 observations in Band 7 with angular resolution down to $\sim$0.15$"$. Results: We detected 6 lines of $^{13}$CH$_3$OH and 13 lines of CH$_2$DOH with upper level energies up to 438 K in temperature units. We derived a rotational temperature of (171 $\pm$ 52) K and column densities of 7$\times$10$^{16}$ cm$^{-2}$ ($^{13}$CH$_3$OH) and 1$\times$10$^{17}$ cm$^{-2}$ (CH$_2$DOH), respectively. Consequently, the D/H ratio is (2.4 $\pm$ 0.4)$\times$10$^{-2}$, a value lower by an order of magnitude with respect to what was previously measured using single dish telescopes toward protostars located in Perseus. Our findings are consistent with the higher dust temperatures in Orion B with respect to that derived for the Perseus cloud. The emission is tracing a rotating structure extending up to 45 au from the jet axis and elongated by 90 au along the jet axis. So far, the origin of the observed emission appears to be related with the accretion disk. Only higher spatial resolution measurements however, will be able to disentangle between different possible scenarios: disk wind, disk atmosphere, or accretion shocks.

Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations.

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

We have investigated the global dynamical state of the Integral Shaped Filament in the Orion A cloud using new N$_2$H$^+$ (1-0) large-scale, IRAM30m observations. Our analysis of its internal gas dynamics reveals the presence of accelerated motions towards the Orion Nebula Cluster, showing a characteristic blue-shifted profile centred at the position of the OMC-1 South region. The properties of these observed gas motions (profile, extension, and magnitude) are consistent with the expected accelerations for the gravitational collapse of the OMC-1 region and explain both the physical and kinematic structure of this cloud.

Are the initial conditions for clustered star formation the same as for non-clustered star formation? To investigate the initial gas properties in young proto-clusters we carried out a comprehensive and high-sensitivity study of the internal structure, density, temperature, and kinematics of the dense gas content of the NGC1333 region in Perseus, one of the nearest and best studied embedded clusters. The analysis of the gas velocities in the Position-Position-Velocity space reveals an intricate underlying gas organization both in space and velocity. We identified a total of 14 velocity-coherent, (tran-)sonic structures within NGC1333, with similar physical and kinematic properties than those quiescent, star-forming (aka fertile) fibers previously identified in low-mass star-forming clouds. These fibers are arranged in a complex spatial network, build-up the observed total column density, and contain the dense cores and protostars in this cloud. Our results demonstrate that the presence of fibers is not restricted to low-mass clouds but can be extended to regions of increasing mass and complexity. We propose that the observational dichotomy between clustered and non-clustered star-forming regions might be naturally explained by the distinct spatial density of fertile fibers in these environments.

[Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms{}) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,{\sc i}] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,{\sc i}] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,{\sc i}] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources.

Although surveys of infall motions in dense cores have been carried out for years, few surveys have focused on mapping infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30-m Telescope maps of N2H+(1-0), DCO+(2-1), DCO+(3-2), and HCO+(3-2) emission towards two prestellar cores (L492 and L694-2) and one protostellar core (L1521F). We find that the measured infall velocity varies with position across each core and choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Line-of-sight infall speeds estimated from DCO+(2-1) line profiles can decrease by 40-50 m/s when observing at a radial offset >= 0.04 pc from the core's dust continuum emission peak. Median infall speeds calculated from all observed positions across a core can also vary by as much as 65 m/s depending on the transition. These results show that while single-pointing, single-transition surveys of core infall velocities may be good indicators of whether a core is either contracting or expanding, the magnitude of the velocities they measure are significantly impacted by the choice of molecular line, proximity to the core center, and core evolutionary state.

The mechanisms leading to the formation of disks around young stellar objects (YSOs) and to the launching of the associated jets are crucial to the understanding of the earliest stages of star and planet formation. HH 212 is a privileged laboratory to study a pristine jet-disk system. Therefore we investigate the innermost region ($<100$ AU) around the HH 212-MM1 protostar through ALMA band\,7 observations of methanol. The 8 GHz bandwidth spectrum towards the peak of the continuum emission of the HH 212 system reveals at least 19 transitions of methanol. Several of these lines (among which several vibrationally excited lines in the v$_{\rm t}=1,2$ states) have upper energies above 500 K. They originate from a compact ($<135$ AU in diameter), hot ($\sim 295$ K) region elongated along the direction of the SiO jet. We performed a fit in the $uv$ plane of various velocity channels of the strongest high-excitation lines. The blue- and red-shifted velocity centroids are shifted roughly symmetrically on either side of the jet axis, indicating that the line-of-sight velocity beyond 0.7 km s$^{-1}$ from systemic is dominated by rotational motions. The velocity increases moving away from the protostar further indicating that the emission of methanol is not associated with a Keplerian disk or rotating-infalling cavity, and it is more likely associated with outflowing gas. We speculate that CH$_3$OH traces a disk wind gas accelerated at the base. The launching region would be at a radius of a few astronomical units from the YSO.

$Aims.$ We study the relation between the jet and the outflow in the IRAS 04166+2706 protostar. This Taurus protostar drives a molecular jet that contains multiple emission peaks symmetrically located from the central source. The protostar also drives a wide-angle outflow consisting of two conical shells. $Methods.$ We have used the Atacama Large Millimeter/submillimeter Array (ALMA) interferometer to observe two fields along the IRAS 04166+2706 jet. The fields were centered on a pair of emission peaks that correspond to the same ejection event, and were observed in CO(2-1), SiO(5-4), and SO(65-54). $ Results.$ Both ALMA fields present spatial distributions that are approximately elliptical and have their minor axes aligned with the jet direction. As the velocity increases, the emission in each field moves gradually across the elliptical region. This systematic pattern indicates that the emitting gas in each field lies in a disk-like structure that is perpendicular to the jet axis and is expanding away from the jet. A small degree of curvature in the first-moment maps indicates that the disks are slightly curved in the manner expected for bow shocks moving away from the IRAS source. A simple geometrical model confirms that this scenario fits the main emission features. $Conclusions.$ The emission peaks in the IRAS 04166+2706 jet likely represent internal bow shocks where material is being ejected laterally away from the jet axis. While the linear momentum of the ejected gas is dominated by the component in the jet direction, the sideways component is not negligible, and can potentially affect the distribution of gas in the surrounding outflow and core.

Sulfur-bearing molecules are highly reactive in the gas phase of the ISM. However, the form in which most of the sulfur is locked onto interstellar dust grains is unknown. By taking advantage of the short time-scales of shocks in young molecular outflows, one could track back the main form of sulfur in the ices. In this paper, six transitions of H$_2$S and its isotopologues in the L1157-B1 bowshock have been detected using data from the Herschel-CHESS survey and the IRAM-30m ASAI large program. These detections are used to calculate the properties of H$_2$S gas in L1157-B1 through use of a rotation diagram and to explore the possible carriers of sulfur on the grains. The isotopologue detections allow the first calculation of the H$_2$S deuteration fraction in an outflow from a low mass protostar. The fractional abundance of H$_2$S in the region is found to be 6.0$\times$10$^{-7}$ and the deuteration fraction is 2$\times$10$^{-2}$. In order to investigate the form of sulfur on the grains, a chemical model is run with four different networks, each with different branching ratios for the freeze out of sulfur bearing species into molecules such as OCS and H$_2$S. It is found that the model best fits the data when at least half of each sulfur bearing species hydrogenates when freezing. We therefore conclude that a significant fraction of sulfur in L1157-B1 is likely to be locked in H$_2$S on the grains.

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar.

Aims: Using the unprecedented combination of high resolution and sensitivity offered by ALMA, we aim to investigate whether and how hot corinos, circumstellar disks, and ejected gas are related in young solar-mass protostars. Methods: We observed CH$_3$CHO and deuterated water (HDO) high-excitation ($E_{\rm u}$ up to 335 K) lines towards the Sun-like protostar HH212--MM1. Results: For the first time, we have obtained images of CH$_3$CHO and HDO emission in the inner $\simeq$ 100 AU of HH212. The multifrequency line analysis allows us to contrain the density ($\geq$ 10$^{7}$ cm$^{-3}$), temperature ($\simeq$ 100 K), and CH$_3$CHO abundance ($\simeq$ 0.2--2 $\times$ 10$^{-9}$) of the emitting region. The HDO profile is asymmetric at low velocities ($\leq$ 2 km s$^{-1}$ from $V_{\rm sys}$). If the HDO line is optically thick, this points to an extremely small ($\sim$ 20--40 AU) and dense ($\ge$ 10$^{9}$ cm$^{-3}$) emitting region. Conclusions: We report the first detection of a hot corino in Orion. The HDO asymmetric profile indicates a contribution of outflowing gas from the compact central region, possibly associated with a dense disk wind.

Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its $\sim$6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e., $\sigma_{NT}/c_s\lesssim1$) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the ISM. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson's velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales.

Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecular cloud that is an ideal candidate for a filament at an early stage of fragmentation. We employ dust extinction mapping in conjunction with near-infrared data from the NEWFIRM instrument, and 870 um dust continuum emission data from the LABOCA instrument, to estimate column densities. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud's morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, 1.6 pc long filament in its Center region. The line mass of the filament is 21-31 Ms pc^-1 and FWHM 0.07 pc. Its radial profile can be fitted with a Plummer profile that has the power-index of 2.6 \pm 11%, flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of 5-10 x 10^4 cm^-3. The fragments in the cloud have a mean separation of 0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the cloud ends, leaving its center yet relatively non-fragmented, possibly because of gravitational focusing in a finite geometry.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

To investigate the disk formation and jet launch in protostars is crucial to comprehend the earliest stages of star and planet formation. We aim to constrain the properties of the molecular jet and the disk of the HH 212 protostellar system at unprecedented angular scales through ALMA observations of sulfur-bearing molecules, SO 9(8)-8(7), SO 10(11)-10(10), SO2 8(2,6)-7(1,7). SO 9(8)-8(7) and SO2 8(2,6)-7(1,7) show broad velocity profiles. At systemic velocity they probe the circumstellar gas and the cavity walls. Going from low to high blue-/red-shifted velocities the emission traces the wide-angle outflow and the fast (~100-200 km/s) and collimated (~90 AU) molecular jet revealing the inner knots with timescales <50 years. The jet transports a mass loss rate >0.2-2e-6 Msun/yr, implying high ejection efficiency (>0.03-0.3). The SO and SO2 abundances in the jet are ~1e-7-1e-6. SO 10(11)-10(10) emission is compact and shows small-scale velocity gradients indicating that it originates partly from the rotating disk previously seen in HCO+ and C17O, and partly from the base of the jet. The disk mass is >0.002-0.013 Msun, and the SO abundance in the disk is ~1e-8-1e-7. SO and SO2 are effective tracers of the molecular jet in the inner few hundreds AU from the protostar. Their abundances indicate that 1% - 40% of sulfur is in SO and SO2 due to shocks in the jet/outflow and/or to ambipolar diffusion at the wind base. The SO abundance in the disk is 3-4 orders of magnitude larger than in evolved protoplanetary disks. This may be due to an SO enhancement in the accretion shock at the envelope-disk interface or in spiral shocks if the disk is partly gravitationally unstable.

Formamide (NH2CHO) has been proposed as a pre-biotic precursor with a key role in the emergence of life on Earth. While this molecule has been observed in space, most of its detections correspond to high-mass star-forming regions. Motivated by this lack of investigation in the low-mass regime, we searched for formamide, as well as isocyanic acid (HNCO), in 10 low- and intermediate-mass pre-stellar and protostellar objects. The present work is part of the IRAM Large Programme ASAI (Astrochemical Surveys At IRAM), which makes use of unbiased broadband spectral surveys at millimetre wavelengths. We detected HNCO in all the sources and NH2CHO in five of them. We derived their abundances and analysed them together with those reported in the literature for high-mass sources. For those sources with formamide detection, we found a tight and almost linear correlation between HNCO and NH2CHO abundances, with their ratio being roughly constant -between 3 and 10- across 6 orders of magnitude in luminosity. This suggests the two species are chemically related. The sources without formamide detection, which are also the coldest and devoid of hot corinos, fall well off the correlation, displaying a much larger amount of HNCO relative to NH2CHO. Our results suggest that, while HNCO can be formed in the gas phase during the cold stages of star formation, NH2CHO forms most efficiently on the mantles of dust grains at these temperatures, where it remains frozen until the temperature rises enough to sublimate the icy grain mantles. We propose hydrogenation of HNCO as a likely formation route leading to NH2CHO.

Context: Outflows are an important part of the star formation process as both the result of ongoing active accretion and one of the main sources of mechanical feedback on small scales. Water is the ideal tracer of these effects because it is present in high abundance in various parts of the protostar. Method: We present \textit{Herschel} HIFI spectra of multiple water-transitions towards 29 nearby Class 0/I protostars as part of the WISH Survey. These are decomposed into different Gaussian components, with each related to one of three parts of the protostellar system; quiescent envelope, cavity shock and spot shocks in the jet and at the base of the outflow. We then constrain the excitation conditions present in the two outflow-related components. Results: Water emission is optically thick but effectively thin, with line ratios that do not vary with velocity, in contrast to CO. The physical conditions of the cavity and spot shocks are similar, with post-shock H$_{2}$ densities of order 10$^{5}-$10$^{8}$\,cm$^{-3}$ and H$_{2}$O column densities of order 10$^{16}-$10$^{18}$\,cm$^{-2}$. H$_{2}$O emission originates in compact emitting regions: for the spot shocks these correspond to point sources with radii of order 10-200\,AU, while for the cavity shocks these come from a thin layer along the outflow cavity wall with thickness of order 1-30\,AU. Conclusions: Water emission at the source position traces two distinct kinematic components in the outflow; J shocks at the base of the outflow or in the jet, and C shocks in a thin layer in the cavity wall. Class I sources have similar excitation conditions to Class 0 sources, but generally smaller line-widths and emitting region sizes. We suggest that it is the velocity of the wind driving the outflow, rather than the decrease in envelope density or mass, that is the cause of the decrease in H$_{2}$O intensity between Class 0 and I.

L1157-mm powers a molecular outflow that is well-known for its shock-induced chemical activity in several hot-spots. We have studied the molecular emission toward L1157-mm searching for a jet component responsible for these spots. We used the IRAM 30m telescope to observe the vicinity of L1157-mm in several lines of SiO. The SiO(5-4) and SiO(6-5) spectra toward L1157-mm present blue and red detached components about 45 km/s away from the ambient cloud. These extremely high-velocity (EHV) components are similar to those found in the L1448 and IRAS 04166+2706 outflows, and probably arise from a molecular jet driven by L1157-mm. Observations of off-center positions indicate that the jet is unresolved in SiO(5-4) (<11"). The EHV jet seen in SiO probably excites L1157-B1 and the other chemically active spots of the L1157 outflow.

(Abridged) We study the kinematics of the dense gas in the Taurus L1495/B213 filamentary region to investigate the mechanism of core formation. We use observations of N2H+(1-0) and C18O(2-1) carried out with the IRAM 30m telescope. We find that the dense cores in L1495/B213 are significantly clustered in linear chain-like groups about 0.5pc long. The internal motions in these chains are mostly subsonic and the velocity is continuous, indicating that turbulence dissipation in the cloud has occurred at the scale of the chains and not at the smaller scale of the individual cores. The chains also present an approximately constant abundance of N2H+ and radial intensity profiles that can be modeled with a density law that follows a softened power law. A simple analysis of the spacing between the cores using an isothermal cylinder model indicates that the cores have likely formed by gravitational fragmentation of velocity-coherent filaments. Combining our analysis of the cores with our previous study of the large-scale C18O emission from the cloud, we propose a two-step scenario of core formation in L1495/B213. In this scenario, named "fray and fragment," L1495/B213 originated from the supersonic collision of two flows. The collision produced a network of intertwined subsonic filaments or fibers ("fray" step). Some of these fibers accumulated enough mass to become gravitationally unstable and fragment into chains of closely-spaced cores. This scenario may also apply to other regions of star formation.

The HH54 shock is a Herbig-Haro object, located in the nearby Chamaeleon II cloud. Observed CO line profiles are due to a complex distribution in density, temperature, velocity, and geometry. Resolving the HH54 shock wave in the far-infrared cooling lines of CO constrain the kinematics, morphology, and physical conditions of the shocked region. We used the PACS and SPIRE instruments on board the Herschel space observatory to map the full FIR spectrum in a region covering the HH54 shock wave. Complementary Herschel-HIFI, APEX, and Spitzer data are used in the analysis as well. The observed features in the line profiles are reproduced using a 3D radiative transfer model of a bow-shock, constructed with the Line Modeling Engine code (LIME). The FIR emission is confined to the HH54 region and a coherent displacement of the location of the emission maximum of CO with increasing J is observed. The peak positions of the high-J CO lines are shifted upstream from the lower J CO lines and coincide with the position of the spectral feature identified previously in CO(10-9) profiles with HIFI. This indicates a hotter molecular component in the upstream gas with distinct dynamics. The coherent displacement with increasing J for CO is consistent with a scenario where IRAS12500-7658 is the exciting source of the flow, and the 180 K bow-shock is accompanied by a hot (800 K) molecular component located upstream from the apex of the shock and blueshifted by -7 km s$^{-1}$. The spatial proximity of this knot to the peaks of the atomic fine-structure emission lines observed with Spitzer and PACS ([OI]63, 145 $\mu$m) suggests that it may be associated with the dissociative shock as the jet impacts slower moving gas in the HH54 bow-shock.

Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario.

Context: Herschel observations suggest that the H$_2$O distribution in outflows from low-mass stars resembles the H$_2$ emission. It is still unclear which of the different excitation components that characterise the mid- and near-IR H$_2$ distribution is associated with H$_2$O. Aim: The aim is to spectrally resolve the different excitation components observed in the H$_2$ emission. This will allow us to identify the H$_2$ counterpart associated with H$_2$O and finally derive directly an H$_2$O abundance estimate with respect to H$_2$. Methods: We present new high spectral resolution observations of H$_2$ 0-0 S(4), 0-0 S(9), and 1-0 S(1) towards HH 54, a bright nearby shock region in the southern sky. In addition, new Herschel-HIFI H$_2$O (2$_{12}$$-$1$_{01}$) observations at 1670~GHz are presented. Results: Our observations show for the first time a clear separation in velocity of the different H$_2$ lines: the 0-0 S(4) line at the lowest excitation peaks at $-$7~km~s$^{-1}$, while the more excited 0-0 S(9) and 1-0 S(1) lines peak at $-$15~km~s$^{-1}$. H$_2$O and high-$J$ CO appear to be associated with the H$_2$ 0-0 S(4) emission, which traces a gas component with a temperature of 700$-$1000 K. The H$_2$O abundance with respect to H$_2$ 0-0 S(4) is estimated to be $X$(H$_2$O)$<$1.4$\times$10$^{-5}$ in the shocked gas over an area of 13$^{\prime\prime}$. Conclusions: We resolve two distinct gas components associated with the HH 54 shock region at different velocities and excitations. This allows us to constrain the temperature of the H$_2$O emitting gas ($\leq$1000 K) and to derive correct estimates of H$_2$O abundance in the shocked gas, which is lower than what is expected from shock model predictions.

Context: Because it is viewed simply edge-on, the HH212 protostellar system is an ideal laboratory for studying the interplay of infall, outflow, and rotation in the earliest stages of low-mass star formation. Aims: We wish to exploit the unmatched combination of high angular resolution, high sensitivity, high-imaging fidelity, and spectral coverage provided by ALMA to shed light on the complex kinematics of the innermost central regions of HH212. Methods: We mapped the inner 10" (4500 AU) of the HH212 system at about 0.5 arcsec resolution in several molecular tracers and in the 850 $\mu$m dust continuum using the ALMA interferometer in band 7 in the extended configuration of the Early Science Cycle 0 operations. Results: Within a single ALMA spectral set-up, we simultaneously identify all the crucial ingredients known to be involved in the star formation recipe namely: (i) the fast, collimated bipolar SiO jet driven by the protostar, (ii) the large-scale swept-up CO outflow, (iii) the flattened rotating and infalling envelope, with bipolar cavities carved by the outflow (in C$^{17}$O(3--2)), and (iv) a rotating wide-angle flow that fills the cavities and surrounding the axial jet (in C$^{34}$S(7--6)). In addition, the compact high-velocity C$^{17}$O emission ($\pm$ 1.9--3.5 km s$^{-1}$ from systemic) shows a velocity gradient along the equatorial plane consistent with a rotating disk of about 0.2 arcsec = 90 AU around a $\simeq 0.3 \pm 0.1 M_{\rm \odot}$ source. The rotating disk is possibly Keplerian. Conclusions: HH212 is the third Class 0 protostar with possible signatures of a Keplerian disk of radius $\geq 30 AU$. The warped geometry in our CS data suggests that this large keplerian disk might result from misaligned magnetic and rotation axes during the collapse phase. The wide-angle CS flow suggests that disk winds may be present in this source.

We present the study of the H2O spatial distribution at two bright shocked regions along IRAS4A, one of the strongest H2O emitters among the Class 0 outflows. We obtained Herschel-PACS maps of the IRAS4A outflow and HIFI observations of two shocked positions. The largest HIFI beam of 38 arcsec at 557 GHz was mapped in several key water lines with different upper energy levels, to reveal possible spatial variations of the line profiles. We detect four H2O lines and CO (16-15) at the two selected positions. In addition, transitions from related outflow and envelope tracers are detected. Different gas components associated with the shock are identified in the H2O emission. In particular, at the head of the red lobe of the outflow, two distinct gas components with different excitation conditions are distinguished in the HIFI emission maps: a compact component, detected in the ground-state water lines, and a more extended one. Assuming that these two components correspond to two different temperature components observed in previous H2O and CO studies, the excitation analysis of the H2O emission suggests that the compact (about 3 arcsec) component is associated with a hot (T~1000 K) gas with densities ~(1-4)x10^5 cm^{-3}, whereas the extended one (10-17 arcsec) traces a warm (T~300-500 K) and dense gas (~(3-5)x10^7 cm^{-3}). Finally, using the CO (16-15) emission observed at R2, we estimate the H2O/H2 abundance of the warm and hot components to be (7-10)x10^{-7} and (3-7)x10^{-5}. Our data allowed us, for the first time, to resolve spatially the two temperature components previously observed with HIFI and PACS. We propose that the compact hot component may be associated with the jet that impacts the surrounding material, whereas the warm, dense, and extended component originates from the compression of the ambient gas by the propagating flow.

We have investigated the water abundance in shock-heated molecular gas, making use of $Herschel$ measurements of far-infrared CO and H$_2$O line emissions in combination with $Spitzer$ measurements of mid-IR H$_2$ rotational emissions. We present far-infrared line spectra obtained with $Herschel$'s PACS instrument in range spectroscopy mode towards two positions in the protostellar outflow NGC 2071 and one position each in the supernova remnants W28 and 3C391. These spectra provide unequivocal detections, at one or more positions, of 12 rotational lines of water, 14 rotational lines of CO, 8 rotational lines of OH (4 lambda doublets), and 7 fine-structure transitions of atoms or atomic ions. We first used a simultaneous fit to the CO line fluxes, along with H$_2$ rotational line fluxes measured previously by $Spitzer$, to constrain the temperature and density distribution within the emitting gas; and we then investigated the water abundances implied by the observed H$_2$O line fluxes. The water line fluxes are in acceptable agreement with standard theoretical models for nondissociative shocks that predict the complete vaporization of grain mantles in shocks of velocity $v \sim 25$ km/s, behind which the characteristic gas temperature is $\sim 1300$ K and the H$_2$O/CO ratio is 1.2

The molecular outflow from IRAS 04166+2706 was mapped with the Submillimeter Array (SMA) at 350 GHz continuum and CO J = 3$-$2 at an angular resolution of ~1 arcsec. The field of view covers the central arc-minute, which contains the inner four pairs of knots of the molecular jet. On the channel map, conical structures are clearly present in the low velocity range (|V$-$V$_0$|$<$10 km s$^{-1}$), and the highly collimated knots appear in the Extremely High Velocity range (EHV, 50$>$|V$-$V$_0$|$>$30 km s$^{-1}$). The higher angular resolution of ~1 arcsec reveals the first blue-shifted knot (B1) that was missing in previous PdBI observation of Sant\'iago-Garc\'ia et al. (2009) at an offset of ~6 arcsec to the North-East of the central source. This identification completes the symmetric sequence of knots in both the blue- and red-shifted lobes of the outflow. The innermost knots R1 and B1 have the highest velocities within the sequence. Although the general features appear to be similar to previous CO J = 2$-$1 images in Sant\'iago-Garc\'ia et al. (2009), the emission in CO J = 3$-$2 almost always peaks further away from the central source than that of CO J = 2$-$1 in the red-shifted lobe of the channel maps. This gives rise to a gradient in the line-ratio map of CO J = 3$-$2/J = 2$-$1 from head to tail within a knot. A large velocity gradient (LVG) analysis suggests that the differences may reflect a higher gas kinetic temperature at the head. We also explore possible constraints imposed by the non-detection of SiO J = 8$-$7.

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.

In the framework of the WISH key program, several H2O (E_u>190 K), high-J CO, [OI], and OH transitions are mapped with PACS in two shock positions along the two prototypical low-luminosity outflows L1448 and L1157. Previous HIFI H2O observations (E_u=53-249 K) and complementary Spitzer mid-IR H2 data are also used, with the aim of deriving a complete picture of the excitation conditions. At all selected spots a close spatial association between H2O, mid-IR H2, and high-J CO emission is found, whereas the low-J CO emission traces either entrained ambient gas or a remnant of an older shock. The excitation analysis at L1448-B2 suggests that a two-component model is needed to reproduce the H2O, CO, and mid-IR H2 lines: an extended warm component (T~450 K) is traced by the H2O emission with E_u =53-137 K and by the CO lines up to J=22-21, and a compact hot component (T=1100 K) is traced by the H2O emission with E_u>190 K and by the higher-J CO lines. At L1448-B2 we obtain an H2O abundance (3-4)x10^{-6} for the warm component and (0.3-1.3)x10^{-5} for the hot component; we also detect OH and blue-shifted [OI] emission, spatially coincident with the other molecular lines and with [FeII] emission. This suggests a dissociative shock for these species, related to the embedded atomic jet. On the other hand, a non-dissociative shock at the point of impact of the jet on the cloud is responsible for the H2O and CO emission. The other examined shock positions show an H2O excitation similar to L1448-B2, but a slightly higher H2O abundance (a factor of 4). The two gas components may represent a gas stratification in the post-shock region. The extended and low-abundance warm component traces the post-shocked gas that has already cooled down to a few hundred Kelvin, whereas the compact and possibly higher-abundance hot component is associated with the gas that is currently undergoing a shock episode.

The evolution of deeply embedded young stellar objects (YSOs) has traditionally been traced through dust continuum spectral energy distributions (SEDs), but the use of CO excitation as an evolutionary probe has not yet been explored due to lack of high-J CO observations. The aim is to constrain the physical characteristics (excitation, kinematics, column density) of the warm gas toward 26 low-mass Class 0 and I YSOs using spectrally-resolved Herschel Space Observatory data of high-J lines of CO. Data are complemented by ground-based observations from APEX and the JCMT to compare those with the colder gas traced by lower-J CO lines. This is the first large spectrally resolved high-J CO survey conducted for these types of sources. The median excitation temperatures for 12CO, 13CO and C18O derived from single-temperature fits to the J_u=2-10 integrated intensities are ~70K, 48K and 37K, respectively, with no significant difference between Class 0 and I sources and no trend with M_env or L_bol. Thus, in contrast with the continuum SEDs, the spectral line energy distributions (SLEDs) do not show an evolution during the embedded stage. In contrast, the integrated line intensities of all CO isotopologs show a clear decrease with evolutionary stage as the envelope is dispersed. The H2O 110-101/CO 10-9 intensity ratio does not change significantly with velocity, in contrast with the H2O/CO 3-2 ratio, indicating that CO 10-9 is the lowest transition for which the line wings probe the same warm shocked gas as H2O. Modeling of the full suite of C18O lines indicates an abundance profile for Class 0 sources that is consistent with a freeze-out zone below 25 K and evaporation at higher temperatures, with but some fraction of the CO transformed into other species in the cold phase. In contrast, observations for two Class I sources in Ophiuchus are consistent with a constant high CO abundance profile.

Context. A small group of bipolar protostellar outflows display strong emission from shock-tracer molecules such as SiO and CH3OH, and are generally referred to as "chemically active." The best-studied outflow from this group is the one in L 1157. Aims. We study the molecular emission from the bipolar outflow powered by the very young stellar object HH 114 MMS and compare its chemical composition with that of the L1157 outflow. Methods. We have used the IRAM 30m radio telescope to observe a number of transitions from CO, SiO, CH3OH, SO, CS, HCN, and HCO+ toward the HH 114 MMS outflow. The observations consist of maps and a two-position molecular survey. Results. The HH 114 MMS outflow presents strong emission from a number of shock-tracer molecules that dominate the appearance of the maps around the central source. The abundance of these molecules is comparable to the abundance in L 1157. Conclusions. The outflow from HH 114 MMS is a spectacular new case of a chemically active outflow.

Context. The observed physical properties of outflows from low-mass sources put constraints on possible ejection mechanisms. Historically, these quantities have been derived from CO using ground-based observations. It is thus important to investigate whether parameters such as momentum rate (thrust) and mechanical luminosity (power) are the same when different molecular tracers are used. Aims. We aim at determining the outflow momentum, dynamical time-scale, thrust, energy and power using CO and H2O as tracers of outflow activity. Methods. Within the framework of the WISH key program, three molecular outflows from Class 0 sources have been mapped using the HIFI instrument aboard Herschel. We use these observations together with previously published H2 data to infer the physical properties of the outflows. We compare the physical properties derived here with previous estimates based on CO observations. Results. Inspection of the spatial distribution of H2O and H2 confirms that these molecules are co-spatial. The most prominent emission peaks in H2 coincide with strong H2O emission peaks and the estimated widths of the flows when using the two tracers are comparable. Conclusions. For the momentum rate and the mechanical luminosity, inferred values are independent of which tracer that is used, i.e., the values agree to within a factor of 4 and 3 respectively.

(Abridged) Context. Core condensation is a critical step in the star-formation process, but is still poorly characterized observationally. Aims. We have studied the 10 pc-long L1495/B213 complex in Taurus to investigate how dense cores have condensed out of the lower-density cloud material. Results. From the N$_2$H$^+$ emission, we identify 19 dense cores, some starless and some protostellar. They are not distributed uniformly, but tend to cluster with relative separations on the order of 0.25 pc. From the C$^{18}$O emission, we identify multiple velocity components in the gas. We have characterized them by fitting gaussians to the spectra, and by studying the distribution of the fits in position-position-velocity space. In this space, the C$^{18}$O components appear as velocity-coherent structures, and we have identified them automatically using a dedicated algorithm (FIVe: Friends In Velocity). Using this algorithm, we have identified 35 filamentary components with typical lengths of 0.5 pc, sonic internal velocity dispersions, and mass-per-unit-length close to the stability threshold of isothermal cylinders at 10 K. Core formation seems to have occurred inside the filamentary components via fragmentation, with a small number of fertile components with larger mass-per-unit-length being responsible for most cores in the cloud. At large scales, the filamentary components appear grouped into families, which we refer to as bundles. Conclusions. Core formation in L1495/B213 has proceeded by hierarchical fragmentation. The cloud fragmented first into several pc-scale regions. Each of these regions later fragmented into velocity-coherent filaments of about 0.5 pc in length. Finally, a small number of these filaments fragmented quasi-statically and produced the individual dense cores we see today.

(Abridged) We present a survey of the water emission in a sample of more than 20 outflows from low mass young stellar objects with the goal of characterizing the physical and chemical conditions of the emitting gas. We have used the HIFI and PACS instruments on board the Herschel Space Observatory to observe the two fundamental lines of ortho-water at 557 and 1670 GHz. These observations were part of the "Water In Star-forming regions with Herschel" (WISH) key program, and have been complemented with CO and H2 data. We find that the emission from water has a different spatial and velocity distribution from that of the J=1-0 and 2-1 transitions of CO, but it has a similar spatial distribution to H2, and its intensity follows the H2 intensity derived from IRAC images. This suggests that water traces the outflow gas at hundreds of kelvins responsible for the H2 emission, and not the component at tens of kelvins typical of low-J CO emission. A warm origin of the water emission is confirmed by a remarkable correlation between the intensities of the 557 and 1670 GHz lines, which also indicates the emitting gas has a narrow range of excitations. A non-LTE radiative transfer analysis shows that while there is some ambiguity on the exact combination of density and temperature values, the gas thermal pressure nT is constrained within less than a factor of 2. The typical nT over the sample is 4 10^{9} cm^{-3}K, which represents an increase of 10^4 with respect to the ambient value. The data also constrain within a factor of 2 the water column density. When this quantity is combined with H2 column densities, the typical water abundance is only 3 10^{-7}, with an uncertainty of a factor of 3. Our data challenge current C-shock models of water production due to a combination of wing-line profiles, high gas compressions, and low abundances.

(Abridged) Far-infrared Herschel-PACS spectra of 18 low-mass protostars of various luminosities and evolutionary stages are studied. We quantify their far-infrared line emission and the contribution of different atomic and molecular species to the gas cooling budget during protostellar evolution. We also determine the spatial extent of the emission and investigate the underlying excitation conditions. Most of the protostars in our sample show strong atomic and molecular far-infrared emission. Water is detected in 17 objects, including 5 Class I sources. The high-excitation H2O line at 63.3 micron is detected in 7 sources. CO transitions from J=14-13 up to 49-48 are found and show two distinct temperature components on Boltzmann diagrams with rotational temperatures of ~350 K and ~700 K. H2O has typical excitation temperatures of ~150 K. Emission from both Class 0 and I sources is usually spatially extended along the outflow direction but with a pattern depending on the species and the transition. The H2O line fluxes correlate strongly with those of the high-J CO lines, as well as with the bolometric luminosity and envelope mass. They correlate less strongly with OH and not with [OI] fluxes. The PACS data probe at least two physical components. The H2O and CO emission likely arises in non-dissociative (irradiated) shocks along the outflow walls with a range of pre-shock densities. Some OH is also associated with this component, likely resulting from H2O photodissociation. UV-heated gas contributes only a minor fraction to the CO emission observed by PACS, based on the strong correlation between the shock-dominated CO 24-23 line and the CO 14-13 line. [OI] and some of the OH emission probe dissociative shocks in the inner envelope. The total far-infrared cooling is dominated by H2O and CO, with [OI] increasing for Class I sources.

We investigate on the spatial and velocity distribution of H2O along the L1448 outflow, its relationship with other tracers, and its abundance variations, using maps of the o-H2O 1_{10}-1_{01} and 2_{12}-1_{01} transitions taken with the Herschel-HIFI and PACS instruments, respectively. Water emission appears clumpy, with individual peaks corresponding to shock spots along the outflow. The bulk of the 557 GHz line is confined to radial velocities in the range \pm 10-50 km/s but extended emission associated with the L1448-C extreme high velocity (EHV) jet is also detected. The H2O 1_{10}-1_{01}/CO(3-2) ratio shows strong variations as a function of velocity that likely reflect different and changing physical conditions in the gas responsible for the emissions from the two species. In the EHV jet, a low H2O/SiO abundance ratio is inferred, that could indicate molecular formation from dust free gas directly ejected from the proto-stellar wind. We derive averaged Tkin and n(H2) values of about 300-500 K and 5 10^6 cm-3 respectively, while a water abundance with respect to H2 of the order of 0.5-1 10^{-6} along the outflow is estimated. The fairly constant conditions found all along the outflow implies that evolutionary effects on the timescales of outflow propagation do not play a major role in the H2O chemistry. The results of our analysis show that the bulk of the observed H2O lines comes from post-shocked regions where the gas, after being heated to high temperatures, has been already cooled down to a few hundred K. The relatively low derived abundances, however, call for some mechanism to diminish the H2O gas in the post-shock region. Among the possible scenarios, we favor H2O photodissociation, which requires the superposition of a low velocity non-dissociative shock with a fast dissociative shock able to produce a FUV field of sufficient strength.

Apart from being an important coolant, H2O is known to be a tracer of high-velocity molecular gas. Recent models predict relatively high abundances behind interstellar shockwaves. The dynamical and physical conditions of the H2O emitting gas, however, are not fully understood yet. We aim to determine the abundance and distribution of H2O, its kinematics and the physical conditions of the gas responsible for the H2O emission. The observed line profile shapes help us understand the dynamics in molecular outflows. We mapped the VLA1623 outflow, in the ground-state transitions of o-H2O, with the HIFI and PACS instruments. We also present observations of higher energy transitions of o-H2O and p-H2O obtained with HIFI and PACS towards selected outflow positions. From comparison with non-LTE radiative transfer calculations, we estimate the physical parameters of the water emitting regions. The observed water emission line profiles vary over the mapped area. Spectral features and components, tracing gas in different excitation conditions, allow us to constrain the density and temperature of the gas. The H2O emission originates in a region where temperatures are comparable to that of the warm H2 gas (T\gtrsim200K). Thus, the H2O emission traces a gas component significantly warmer than the gas responsible for the low-J CO emission. The H2O column densities at the CO peak positions are low, i.e. N(H2O) \simeq (0.03-10)x10e14 cm-2. The H2O abundance with respect to H2 in the extended outflow is estimated at X(H2O)<1x10e-6, significantly lower than what would be expected from most recent shock models. The H2O emission traces a gas component moving at relatively high velocity compared to the low-J CO emitting gas. However, other dynamical quantities such as the momentum rate, energy and mechanical luminosity are estimated to be the same, independent of the molecular tracer used, CO or H2O.

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.

(Abridged) Water is a key tracer of dynamics and chemistry in low-mass protostars, but spectrally resolved observations have so far been limited in sensitivity and angular resolution. In this first systematic survey of spectrally resolved water emission in low-mass protostellar objects, H2O was observed in the ground-state transition at 557 GHz with HIFI on Herschel in 29 embedded Class 0 and I protostars. Complementary far-IR and sub-mm continuum data (including PACS data from our program) are used to constrain the spectral energy distribution of each source. H2O intensities are compared to inferred envelope and outflow properties and CO 3-2 emission. H2O emission is detected in all objects except one. The line profiles are complex and consist of several kinematic components. The profiles are typically dominated by a broad Gaussian emission feature, indicating that the bulk of the water emission arises in outflows, not the quiescent envelope. Several sources show multiple shock components in either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles in 7 sources (6 Class 0, 1 Class I) indicative of infalling envelopes, and regular P-Cygni profiles in 4 sources (3 Class I, 1 Class 0) indicative of expanding envelopes. "Bullets" moving at >50 km/s are seen in 4 Class 0 sources; 3 of these are new detections. In the outflow, the H2O/CO abundance ratio as a function of velocity is nearly the same for all sources, increasing from 10^-3 at <5 km/s to >10^-1 at >10 km/s. The H2O abundance in the outer envelope is low, ~10^-10. The different H2O profile components show a clear evolutionary trend: in the Class 0 sources, emission is dominated by outflow components originating inside an infalling envelope. When the infall diminishes during the Class I phase, the outflow weakens and H2O emission disappears.

Bipolar outflows constitute some of the best laboratories to study shock chemistry in the interstellar medium. A number of molecular species have their abundance enhanced by several orders of magnitude in the outflow gas, likely as a combined result of dust mantle disruption and high temperature gas chemistry, and therefore become sensitive indicators of the physical changes taking place in the shock. Identifying these species and understanding their chemical behavior is therefore of high interest both to chemical studies and to our understanding of the star-formation process. Here we review some of the recent progress in the study of the molecular composition of bipolar outflows, with emphasis in the tracers most relevant for shock chemistry. As we discuss, there has been rapid progress both in characterizing the molecular composition of certain outflows as well as in modeling the chemical processes likely involved. However, a number of limitations still affect our understanding of outflow chemistry. These include a very limited statistical approach in the observations and a dependence of the models on plane-parallel shocks, which cannot reproduce the observed wing morphology of the lines. We finish our contribution by discussing the chemistry of the so-called extremely high velocity component, which seems different from the rest of the outflow and may originate in the wind from the very vicinity of the protostar.

As part of the WISH (Water In Star-forming regions with Herschel) key project, we report on the observations of several ortho- and para-H2O lines performed with the HIFI instrument towards two bright shock spots (R4 and B2) along the outflow driven by the L1448 low-mass proto-stellar system, located in the Perseus cloud. These data are used to identify the physical conditions giving rise to the H2O emission and infer any dependence with velocity. These observations provide evidence that the observed water lines probe a warm (T_kin~400-600 K) and very dense (n 10^6 - 10^7 cm^-3) gas, not traced by other molecules, such as low-J CO and SiO, but rather traced by mid-IR H2 emission. In particular, H2O shows strong differences with SiO in the excitation conditions and in the line profiles in the two observed shocked positions, pointing to chemical variations across the various velocity regimes and chemical evolution in the different shock spots. Physical and kinematical differences can be seen at the two shocked positions. At the R4 position, two velocity components with different excitation can be distinguished, with the component at higher velocity (R4-HV) being less extended and less dense than the low velocity component (R4-LV). H2O column densities of about 2 10^13 and 4 10^14 cm^-2 have been derived for the R4-LV and the R4-HV components, respectively. The conditions inferred for the B2 position are similar to those of the R4-HV component, with H2O column density in the range 10^14 - 5 10^14 cm^-2, corresponding to H2O/H2 abundances in the range 0.5 - 1 10^-5. The observed line ratios and the derived physical conditions seem to be more consistent with excitation in a low velocity J-type shock with large compression rather than in a stationary C-shock, although none of these stationary models seems able to reproduce all the characteristics of the observed emission.

In the framework of the Herschel-WISH key program, several ortho-H2O and para-H2O emission lines, in the frequency range from 500 to 1700 GHz, were observed with the HIFI instrument in two bow-shock regions (B2 and R) of the L1157 cloud. The primary aim is to analyse water emission lines as a diagnostic of the physical conditions in the blue (B2) and red-shifted (R) lobes to compare the excitation conditions. A total of 5 ortho- and para-H216O plus one o-H218O transitions were observed in B2 and R with a wide range of excitation energies (27 K<=Eu<=215 K). The H2O spectra, observed in the two shocked regions, show that the H2O profiles are markedly different in the two regions. In particular, at the bow-shock R, we observed broad (~30 km s-1 with respect to the ambient velocity) red-shifted wings where lines at different excitation peak at different red-shifted velocities. The B2 spectra are associated with a narrower velocity range (~6 km s-1), peaking at the systemic velocity. The excitation analysis suggests, for B2, low values of column density NH2O <=5{\times}1013 cm-2, a density range of 105 <=nH2 <=107 cm-3, and warm temperatures (>=300 K). The presence of the broad red-shifted wings and multiple peaks in the spectra of the R region, prompted the modelling of two components. High velocities are associated with relatively low temperatures (~100K),NH2O{\simeq}5{\times}1012-5{\times}1013 cm-2 and densities nH2{\simeq}106-108 cm-3.Lower velocities are associated with higher excitation conditions with Tkin>=300 K, very dense gas (nH2 ~108 cm-3) and low column density (NH2O<5{\times}1013 cm-2).

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and 12CO(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper limit for the dense core bolometric luminosity of 0.05 Lsun is obtained. By modeling the broadband SED and the continuum interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between these two models, which produce similar fits. We also detect 12CO(2-1) emission with red- and blue-shifted emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. Whatever the true nature of the central object in L1451-mm, this core presents an excellent laboratory to study the earliest phases of low-mass star formation.

We studied the abundance of HCN, H13CN, and HN13C in a sample of prestellar cores, in order to search for species associated with high density gas. We used the IRAM 30m radiotelescope to observe along the major and the minor axes of L1498, L1521E, and TMC 2, three cores chosen on the basis of their CO depletion properties. We mapped the J=1-0 transition of HCN, H13CN, and HN13C towards the source sample plus the J=1-0 transition of N2H+ and the J=2-1 transition of C18O in TMC 2. We used two different radiative transfer codes, making use of recent collisional rate calculations, in order to determine more accurately the excitation temperature, leading to a more exact evaluation of the column densities and abundances. We find that the optical depths of both H13CN(1-0) and HN13C(1-0) are non-negligible, allowing us to estimate excitation temperatures for these transitions in many positions in the three sources. The observed excitation temperatures are consistent with recent computations of the collisional rates for these species and they correlate with hydrogen column density inferred from dust emission. We conclude that HCN and HNC are relatively abundant in the high density zone, n(H2) about 10^5 cm-3, where CO is depleted. The relative abundance [HNC]/[HCN] differs from unity by at most 30 per cent consistent with chemical expectations. The three hyperfine satellites of HCN(1-0) are optically thick in the regions mapped, but the profiles become increasingly skewed to the blue (L1498 and TMC 2) or red (L1521E) with increasing optical depth suggesting absorption by foreground layers.

Context. Low-mass star-forming cores differ from their surrounding molecular cloud in turbulence, shape, and density structure. Aims. We aim to understand how dense cores form out of the less dense cloud material by studying the connection between these two regimes. Methods. We observed the L1517 dark cloud in C18O(1-0), N2H+(1-0), and SO(JN=32-21) with the FCRAO 14m telescope, and in the 1.2mm dust continuum with the IRAM 30m telescope. Results. Most of the gas in the cloud lies in four filaments that have typical lengths of 0.5 pc. Five starless cores are embedded in these filaments and have chemical compositions indicative of different evolutionary stages. The filaments have radial profiles of C18O(1-0) emission with a central flattened region and a power-law tail, and can be fitted approximately as isothermal, pressure-supported cylinders. The filaments, in addition, are extremely quiescent. They have subsonic internal motions and are coherent in velocity over their whole length. The large-scale motions in the filaments can be used to predict the velocity inside the cores, indicating that core formation has not decoupled the dense gas kinematically from its parental material. In two filaments, these large-scale motions consist of oscillations in the velocity centroid, and a simple kinematic model suggests that they may be related to core-forming flows. Conclusions. Core formation in L1517 seems to have occurred in two steps. First, the subsonic, velocity-coherent filaments have condensed out of the more turbulent ambient cloud. Then, the cores fragmented quasi-statically and inherited the kinematics of the filaments. Turbulence dissipation has therefore occurred mostly on scales on the order of 0.5 pc or larger, and seems to have played a small role in the formation of the individual cores.

We are aiming at the observational estimation of the relative contribution to the cooling by CO and H2O, as this provides decisive information for the understanding of the oxygen chemistry behind interstellar shock waves. Methods. The high sensitivity of HIFI, in combination with its high spectral resolution capability, allows us to trace the H2O outflow wings at unprecedented signal-to-noise. From the observation of spectrally resolved H2O and CO lines in the HH52-54 system, both from space and from ground, we arrive at the spatial and velocity distribution of the molecular outflow gas. Solving the statistical equilibrium and non-LTE radiative transfer equations provides us with estimates of the physical parameters of this gas, including the cooling rate ratios of the species. The radiative transfer is based on an ALI code, where we use the fact that variable shock strengths, distributed along the front, are naturally implied by a curved surface. Based on observations of CO and H2O spectral lines, we conclude that the emission is confined to the HH54 region. The quantitative analysis of our observations favours a ratio of the CO-to-H2O-cooling-rate >> 1. From the best-fit model to the CO emission, we arrive at an H2O abundance close to 1e-5. The line profiles exhibit two components, one of which is triangular and another, which is a superposed, additional feature. This additional feature likely originates from a region smaller than the beam where the ortho-water abundance is smaller than in the quiescent gas. Comparison with recent shock models indicate that a planar shock can not easily explain the observed line strengths and triangular line profiles.We conclude that the geometry can play an important role. Although abundances support a scenario where J-type shocks are present, higher cooling rate ratios than predicted by these type of shocks are derived.

Herschel-HIFI observations of water in the low-mass star-forming object L1448-MM, known for its prominent outflow, are presented, as obtained within the `Water in star-forming regions with Herschel' (WISH) key programme. Six H2-16O lines are targeted and detected (E_up/k_B ~ 50-250 K), as is CO J= 10-9 (E_up/k_B ~ 305 K), and tentatively H2-18O 110-101 at 548 GHz. All lines show strong emission in the "bullets" at |v| > 50 km/s from the source velocity, in addition to a broad, central component and narrow absorption. The bullets are seen much more prominently in H$_2$O than in CO with respect to the central component, and show little variation with excitation in H2O profile shape. Excitation conditions in the bullets derived from CO lines imply a temperature >150 K and density >10^5 cm^-3, similar to that of the broad component. The H2O/CO abundance ratio is similar in the "bullets" and the broad component, ~ 0.05-1.0, in spite of their different origins in the molecular jet and the interaction between the outflow and the envelope. The high H2O abundance indicates that the bullets are H2 rich. The H2O cooling in the "bullets" and the broad component is similar and higher than the CO cooling in the same components. These data illustrate the power of Herschel-HIFI to disentangle different dynamical components in low-mass star-forming objects and determine their excitation and chemical conditions.

`Water In Star-forming regions with Herschel' (WISH) is a key program on the Herschel Space Observatory designed to probe the physical and chemical structure of young stellar objects using water and related molecules and to follow the water abundance from collapsing clouds to planet-forming disks. About 80 sources are targeted covering a wide range of luminosities and evolutionary stages, from cold pre-stellar cores to warm protostellar envelopes and outflows to disks around young stars. Both the HIFI and PACS instruments are used to observe a variety of lines of H2O, H218O and chemically related species. An overview of the scientific motivation and observational strategy of the program is given together with the modeling approach and analysis tools that have been developed. Initial science results are presented. These include a lack of water in cold gas at abundances that are lower than most predictions, strong water emission from shocks in protostellar environments, the importance of UV radiation in heating the gas along outflow walls across the full range of luminosities, and surprisingly widespread detection of the chemically related hydrides OH+ and H2O+ in outflows and foreground gas. Quantitative estimates of the energy budget indicate that H2O is generally not the dominant coolant in the warm dense gas associated with protostars. Very deep limits on the cold gaseous water reservoir in the outer regions of protoplanetary disks are obtained which have profound implications for our understanding of grain growth and mixing in disks.

We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

(Abridged) We present a molecular survey of the outflows powered by L1448-mm and IRAS 04166+2706, two sources with prominent wing and extremely high velocity (EHV) components in their CO spectra. The molecular composition of the two outflows presents systematic changes with velocity that we analyze by dividing the outflow in three chemical regimes, two of them associated with the wing component and the other the EHV gas. The analysis of the two wing regimes shows that species like H2CO and CH3OH favor the low-velocity gas, while SiO and HCN are more abundant in the fastest gas. We also find that the EHV regime is relatively rich in O-bearing species, as is not only detected in CO and SiO (already reported elsewhere), but also in SO, CH3OH, and H2CO (newly reported here), with a tentative detection in HCO+. At the same time, the EHV regime is relatively poor in C-bearing molecules like CS and HCN. We suggest that this difference in composition arises from a lower C/O ratio in the EHV gas. The different chemical compositions of the wing and EHV regimes suggest that these two outflow components have different physical origins. The wing component is better explained by shocked ambient gas, although none of the existing shock models explains all observed features. The peculiar composition of the EHV gas may reflect its origin as a dense wind from the protostar or its surrounding disk.

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

The Heterodyne Instrument for the Far Infrared (HIFI) onboard the Herschel Space Observatory allows the first observations of light diatomic molecules at high spectral resolution and in multiple transitions. Here, we report deep integrations using HIFI in different lines of hydrides towards the high-mass star forming region AFGL 2591. Detected are CH, CH+, NH, OH+, H2O+, while NH+ and SH+ have not been detected. All molecules except for CH and CH+ are seen in absorption with low excitation temperatures and at velocities different from the systemic velocity of the protostellar envelope. Surprisingly, the CH(JF,P = 3/2_2,- - 1/2_1,+) and CH+(J = 1 - 0, J = 2 - 1) lines are detected in emission at the systemic velocity. We can assign the absorption features to a foreground cloud and an outflow lobe, while the CH and CH+ emission stems from the envelope. The observed abundance and excitation of CH and CH+ can be explained in the scenario of FUV irradiated outflow walls, where a cavity etched out by the outflow allows protostellar FUV photons to irradiate and heat the envelope at larger distances driving the chemical reactions that produce these molecules.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

Water in Star-forming regions with Herschel (WISH) is a Herschel Key Program investigating the water chemistry in young stellar objects (YSOs) during protostellar evolution. Hydroxyl (OH) is one of the reactants in the chemical network most closely linked to the formation and destruction of H2O. High-temperature chemistry connects OH and H2O through the OH + H2 <-> H2O + H reactions. Formation of H2O from OH is efficient in the high-temperature regime found in shocks and the innermost part of protostellar envelopes. Moreover, in the presence of UV photons, OH can be produced from the photo-dissociation of H2O. High-resolution spectroscopy of the OH 163.12 micron triplet towards HH 46 and NGC 1333 IRAS 2A was carried out with the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel. The low- and intermediate-mass YSOs HH 46, TMR 1, IRAS 15398-3359, DK Cha, NGC 7129 FIRS 2, and NGC 1333 IRAS 2A were observed with the Photodetector Array Camera and Spectrometer (PACS) in four transitions of OH and two [OI] lines. The OH transitions at 79, 84, 119, and 163 micron and [OI] emission at 63 and 145 micron were detected with PACS towards the class I low-mass YSOs as well as the intermediate-mass and class I Herbig Ae sources. No OH emission was detected from the class 0 YSO NGC 1333 IRAS 2A, though the 119 micron was detected in absorption. With HIFI, the 163.12 micron was not detected from HH 46 and only tentatively detected from NGC 1333 IRAS 2A. The combination of the PACS and HIFI results for HH 46 constrains the line width (FWHM > 11 km/s) and indicates that the OH emission likely originates from shocked gas. This scenario is supported by trends of the OH flux increasing with the [OI] flux and the bolometric luminosity. Similar OH line ratios for most sources suggest that OH has comparable excitation temperatures despite the different physical properties of the sources.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

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.

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

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.

The isocyanic acid (HNCO) presents an extended distribution in the centers of the Milky Way and the spiral galaxy IC342. Based on the morphology of the emission and the HNCO abundance with respect to H2, several authors made the hypothesis that HNCO could be a good tracer of interstellar shocks. Here we test this hypothesis by observing a well-known Galactic source where the chemistry is dominated by shocks. We have observed several transitions of HNCO towards L1157-mm and two positions (B1 and B2) in the blue lobe of the molecular outflow. The HNCO line profiles exhibit the same characteristics of other well-known shock tracers like CH3OH, H2CO, SO or SO2. HNCO, together with SO2 and OCS, are the only three molecules detected so far whose emission is much more intense in B2 than in B1, making these species valuable probes of chemical differences along the outflow. The HNCO abundance with respect to H2 is 0.4-1.8 10^-8 in B1 and 0.3-1 10^-7 in B2. These abundances are the highest ever measured, and imply an increment with respect to L1157-mm of a factor up to 83, demonstrating that this molecule is actually a good shock tracer. Our results probe that shocks can actually produce the HNCO abundance measured in galactic nuclei and even higher ones. We propose that the gas phase abundance of HNCO is due both to grain mantles erosion by the shock waves and by neutral-neutral reactions in gas phase involving CN and O2. The observed anticorrelation of CN and HNCO fluxes supports this scenario. The observed similarities of the HNCO emission and the sulfured molecules may arise due to formation pathways involving also O2.

We investigate 35 pre-stellar cores and 36 proto-stellar cores in the Perseus molecular cloud. We find a very tight correlation between the physical parameters describing the N2H+ and NH3 gas. Both the velocity centroids and the line widths of N2H+ and NH3 correlate much better than either species correlates with CO, as expected if the nitrogen-bearing species are probing primarily the dense core gas where the CO has been depleted. We also find a tight correlation in the inferred abundance ratio between N2H+ and para-NH3 across all cores, with N(p-NH3)/N(N2H+)= 22 +/- 10. We find a mild correlation between NH3 (and N2H+) column density and the (sub)millimeter dust continuum derived H2 column density for pre-stellar cores, N(p-NH3)/N(H2) ~ 10e-8, but do not find a fixed ratio for proto-stellar cores. The observations suggest that in the Perseus molecular cloud the formation and destruction mechanisms for the two nitrogen-bearing species are similar, regardless of the physical conditions in the dense core gas. While the equivalence of N2H+ and NH3 as powerful tracers of dense gas is validated, the lack of correspondence between these species and the (sub)millimeter dust continuum observations for proto-stellar cores is disconcerting and presently unexplained.

We study the abundance of CCH in prestellar cores both because of its role in the chemistry and because it is a potential probe of the magnetic field. We also consider the non-LTE behaviour of the N=1-0 and N=2-1 transitions of CCH and improve current estimates of the spectroscopic constants of CCH. We used the IRAM 30m radiotelescope to map the N=1-0 and N=2-1 transitions of CCH towards the prestellar cores L1498 and CB246. Towards CB246, we also mapped the 1.3 mm dust emission, the J=1-0 transition of N2H+ and the J=2-1 transition of C18O. We used a Monte Carlo radiative transfer program to analyse the CCH observations of L1498. We derived the distribution of CCH column densities and compared with the H2 column densities inferred from dust emission. We find that while non-LTE intensity ratios of different components of the N=1-0 and N=2-1 lines are present, they are of minor importance and do not impede CCH column density determinations based upon LTE analysis. Moreover, the comparison of our Monte-Carlo calculations with observations suggest that the non-LTE deviations can be qualitatively understood. For L1498, our observations in conjunction with the Monte Carlo code imply a CCH depletion hole of radius 9 x 10^{16} cm similar to that found for other C-containing species. We briefly discuss the significance of the observed CCH abundance distribution. Finally, we used our observations to provide improved estimates for the rest frequencies of all six components of the CCH(1-0) line and seven components of CCH(2-1). Based on these results, we compute improved spectroscopic constants for CCH. We also give a brief discussion of the prospects for measuring magnetic field strengths using CCH.

We have performed observations of water maser emission towards a sample of low-mass protostars, in order to investigate the properties of jets associated with the earliest stages of star formation and their interaction with the surrounding medium. The main aim is to measure the absolute positions and proper motions of the H_2O spots in order to investigate the kinematics of the region from where the jet is launched. We imaged the protostars in the nearby region NGC 1333-IRAS 4 in the water maser line at 22.2 GHz by using the VLBA in phase-reference mode at the milliarcsecond scale over four epochs, spaced by one month to measure proper motions. Two protostars (A2 and B) were detected in a highly variable H_2O maser emission, with an active phase shorter than four weeks. The H_2O maps allow us to trace the fast jet driven by the B protostar: we observed both the red- and blue-shifted lobes very close to the protostar, =< 35 AU, moving away with projected velocities of ~10-50 km/s. The comparison with the molecular outflow observed at larger scale suggests a jet precession with a 18'/yr rate. By measuring the positional spread of the H_2O spots we estimate a jet width of ~2 AU at a distance of ~12 AU from the driving protostar.

We imaged the protostars of the nearby region NGC1333 IRAS 4 in the water maser line at 22.2 GHz by using the VLBA in phase referencing at milliarcsecond scale over four epochs spaced by one month to measure proper motions. We measure the absolute positions and proper motions of the H2O spots to investigate the kinematics of the region from where the jet is launched. Two protostars (A2 and B) have been detected in a highly variable H2O maser emission, with an active phase shorter than four weeks. A 70 AU chain of several maser spots, very well aligned, has been observed close to the B protostar. The apparent proper motions have been derived, finding that the H2O spots are moving along the N-NW direction with projected velocities between 10 and 50 km/s. We conclude that in IRAS 4B, water maser trace a highly collimated bipolar jet clearly associated with the protostar.

Context: IRAS 04166+2706 in Taurus is one of the most nearby young stellar objects whose molecular outflow contains a highly collimated fast component. Methods: We have observed the IRAS 04166+2706 outflow with the IRAM Plateau de Bure interferometer in CO(J=2-1) and SiO(J=2-1) achieving angular resolutions between 2'' and 4''. To improve the quality of the CO(2-1) images, we have added single dish data to the interferometer visibilities. Results: The outflow consists of two distinct components. At velocities <10 km/s, the gas forms two opposed, approximately conical shells that have the YSO at their vertex. These shells coincide with the walls of evacuated cavities and seem to result from the acceleration of the ambient gas by a wide-angle wind. At velocities >30 km/s, the gas forms two opposed jets that travel along the center of the cavities and whose emission is dominated by a symmetric collection of at least 7 pairs of peaks. The velocity field of this component presents a sawtooth pattern with the gas in the tail of each peak moving faster than the gas in the head. This pattern, together with a systematic widening of the peaks with distance to the central source, is consistent with the emission arising from internal working surfaces traveling along the jet and resulting from variations in the velocity field of ejection. We interpret this component as the true protostellar wind, and we find its composition consistent with a chemical model of such type of wind. Conclusions: Our results support outflow wind models that have simultaneously wide-angle and narrow components, and suggest that the EHV peaks seen in a number of outflows consist of internally-shocked wind material.

We report the detection of complex organic molecules in the young protostellar outflow L1157. We identify lines from HCOOCH3, CH3CN, HCOOH and C2H5OH at the position of the B1 shock in the blueshifted lobe, making it the first time that complex species have been detected towards a molecular outflow powered by a young low-mass protostar. The time scales associated with the warm outflow gas (< 2,000 yr) are too short for the complex molecules to have formed in the gas phase after the shock-induced sputtering of the grain mantles. It is more likely that the complex species formed in the surface of grains and were then ejected from the grain mantles by the shock. The formation of complex molecules in the grains of low-mass star forming regions must be relatively efficient, and our results show the importance of considering the impact of outflows when studying complex molecules around protostars. The relative abundance with respect to methanol of most of the detected complex molecules is similar to that of hot cores and molecular clouds in the galactic center region, which suggests that the mantle composition of the dust in the L1157 dark cloud is similar to dust in those regions.

We survey the kinematics of over one hundred and fifty candidate (and potentially star-forming) dense cores in the Perseus molecular cloud with pointed N2H+(1-0) and simultaneous C18O(2-1) observations. Our detection rate of N2H+ is 62%, rising to 84% for JCMT SCUBA-selected targets. In agreement with previous observations, we find that the dense N2H+ targets tend to display nearly thermal linewidths, particularly those which appear to be starless (using Spitzer data), indicating turbulent support on the small scales of molecular clouds is minimal. For those N2H+ targets which have an associated SCUBA dense core, we find their internal motions are more than sufficient to provide support against the gravitational force on the cores. Comparison of the N2H+ integrated intensity and SCUBA flux reveals fractional N2H+ abundances between 10^-10 and 10^-9. We demonstrate that the relative motion of the dense N2H+ gas and the surrounding C18O gas is less than the sound speed in the vast majority of cases (~90%). The point-to-point motions we observe within larger extinction regions appear to be insufficient to provide support against gravity, although we sparsely sample these regions.

Cold dark clouds are nearby members of the densest and coldest phase in the galactic interstellar medium, and represent the most accessible sites where stars like our Sun are currently being born. In this review we discuss recent progress in their study, including the newly discovered infrared dark clouds that are likely precursors to stellar clusters. At large scales, dark clouds present filamentary mass distributions with motions dominated by supersonic turbulence. At small, sub-parsec scales, a population of subsonic starless cores provides a unique glimpse of the conditions prior to stellar birth. Recent studies of starless cores reveal a combination of simple physical properties together with a complex chemical structure dominated by the freeze-out of molecules onto cold dust grains. Elucidating this combined structure is both an observational and theoretical challenge whose solution will bring us closer to understanding how molecular gas condenses to form stars.

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

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

Dense cores are the simplest star-forming sites that we know, but despite their simplicity, they still hold a number of mysteries that limit our understanding of how solar-type stars form. ALMA promises to revolutionize our knowledge of every stage in the life of a core, from the pre-stellar phase to the final disruption by the newly born star. This contribution presents a brief review of the evolution of dense cores and illustrates particular questions that will greatly benefit from the increase in resolution and sensitivity expected from ALMA

In this paper we present the first SiO multiline analysis (from J=2-1 to J=11-10) of the molecular bullets along the outflows of the Class 0 sources L1448-mm and L1157-mm, obtained through observations with IRAM and JCMT. We have computed the main physical parameters in each bullet and compared them with other tracers of warm and dense gas and with models for the SiO excitation in shocks. We find that the bullets close to L1448--mm, associated with high velocity gas, have higher excitation conditions (n(H2) ~ 10^{6} cm^{-3}, T > 500 K) with respect to the L1157 bullets (n(H2) ~1-5 10^{5} cm^{-3}, T ~ 100-300 K). In both the sources, there is a clear evidence of the presence of velocity components having different excitation conditions, with the denser and/or warmer gas associated with the gas at the higher speed. In L1448 the bulk of the emission is due to the high-excitation and high velocity gas, while in L1157 most of the emission comes from the low excitation gas at ambient velocity. The observed velocity-averaged line ratios are well reproduced by shocks with speeds v_s larger than ~ 30 km/s and densities ~ 10^{5} - 10^{6} cm^{-3}. Plane-parallel shock models, however, fail to predict all the observed line profiles and in particular the very similar profiles shown by both low and high excitation lines. The overall observations support the idea that the L1157 clumps are shock interaction events older than the L1448 bullets close to the driving source. In the latter objects, the velocity structure and the variations of physical parameters with the velocity resemble very closely those found in optical/IR jets near the protostar, suggesting that similar launching and excitation mechanisms are also at the origin of collimated jets seen at millimetre wavelengths.

Aims. We present molecular line observations of the southwestern part of the IC 348 young cluster, and we use them together with NIR and mm continuum data to determine the distribution of dense gas, search for molecular outflows, and analyze the ongoing star formation activity in the region. Methods. Our molecular line data consists of C18O(1--0) and N2H+(1--0) maps obtained with the FCRAO telescope at a resolution of about 50'' and CO(2--1) data obtained with the IRAM 30m telescope at a resolution of 11''. Results. The dense gas southwest of IC 348 is concentrated in two groups of dense cores, each of them with a few solar masses of material and indications of CO depletion at high density. One of the core groups is actively forming stars, while the other seems starless. There is evidence for at least three bipolar molecular outflows in the region, two of them powered by previously identified Class 0 sources, while the other one is powered by a still not well characterized low-luminosity object. The ongoing star formation activity is producing a small stellar subgroup in the cluster. Using the observed core characteristics and the star formation rate in the cluster we propose that that similar episodes of stellar birth may have produced the subclustering observed in the halo of IC 348.

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

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

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

Our understanding of the physical and chemical structure of pre-stellar cores, the simplest star-forming sites, has significantly improved since the last IAU Symposium on Astrochemistry (South Korea, 1999). Research done over these years has revealed that major molecular species like CO and CS systematically deplete onto dust grains at the interior of pre-stellar cores, while species like N2H+ and NH3 survive in the gas phase and can usually be detected towards the core centers. Such a selective behaviour of molecular species gives rise to a differentiated (onion-like) chemical composition, and manifests itself in molecular maps as a dichotomy between centrally peaked and ring-shaped distributions. From the point of view of star-formation studies, the identification of molecular inhomogeneities in cores helps to resolve past discrepancies between observations made using different tracers, and brings the possibility of self-consistent modelling of the core internal structure. Here I present recent work on determining the physical and chemical structure of two pre-stellar cores, L1498 and L1517B, using observations in a large number of molecules and Monte Carlo radiative transfer analysis. These two cores are typical examples of the pre-stellar core population, and their chemical composition is characterized by the presence of large freeze out holes in most molecular species. In contrast with these chemically processed objects, a new population of chemically young cores has started to emerge. The characteristics of its most extreme representative, L1521E, are briefly reviewed.

Dense low mass cores in nearby clouds like Taurus and Auriga are some of the simplest sites currently forming stars like our Sun. Because of their simplicity and proximity, dense cores offer the clearest view of the different phases of star formation, in particular the conditions prior to the onset of gravitational collapse. Thanks to the combined analysis of the emission from molecular lines and the emission/absorption from dust grains, the last several years have seen a very rapid progress in our understanding of the structure and chemical composition of starless cores. Previous contradictions between molecular tracers are now understood to arise from core chemical inhomogeneities, which are caused by the selective freeze out of molecules onto cold dust grains. The analysis of the dust emission and absorption, in addition, has allowed us to derive accurate density profiles, and has made finally possible to carry out self consistent modeling of the internal structure of starless cores. In this paper I briefly review the evolution of core studies previous to the current golden age, and show how multi-tracer emission can now be modeled in a systematic manner. Finally I show how we can start to reconstruct the early history of core formation taking advantage of the chemical changes in the gas.

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

We present the first case of a highly collimated, extremely high velocity bipolar outflow in Taurus. It is powered by the low-luminosity (0.4 L_sun) source IRAS 04166+2706 and contains gas accelerated up to 50 km/s with respect to the ambient cloud both toward the blue and the red (uncorrected for projection). At the highest velocities, the outflow collimation factor exceeds 20, and the gas displays a very high degree of spatial symmetry. This very fast gas presents multiple maxima, and most likely arises from the acceleration of ambient material by a time-variable jet-like stellar wind. When scaled for luminosity, the outflow parameters of IRAS 04166 are comparable to those of other extremely high velocity outflows like L1448, indicating that even the very quiescent star-formation mode of Taurus can produce objects powering very high energy flows (L_mec/L_* > 0.15).

We present new high-sensitivity high-resolution mm-wave observations of the Onsala 1 ultra-compact HII region that bring to light the internal structure of this massive star forming cloud. The 1.2 mm continuum map obtained with the IRAM 30-m radiotelescope (~11" resolution) shows a centrally peaked condensation of 1' size (~0.5 pc at the assumed distance of 1.8 kpc) which has been further investigated at higher resolution in the 3 mm continuum and in the emission lines of H13CO+ J=1--0 and SiO J=2--1 with the IRAM Plateau de Bure interferometer. The 3 mm data, with a resolution of ~5" X 4", displays a unresolved continuum source at the peak of the extended 1.2 mm emission and closely associated with the ultra-compact HII region. The H13CO+ map traces the central condensation in good agreement with previous NH_3 maps of Zheng et al. (1985). However, the velocity field of this central condensation, which was previously thought to arise in a rapidly rotating structure, is better explained in terms of the dense and compact component of a bipolar outflow. This interpretation is confirmed by SiO and CO observations of the full region. In fact, our new SiO data unveils the presence of multiple (at least 4) outflows in the region. In particular, there is an important center of outflow activity in the region at about 1' north of the UCHII region. Indeed the different outflows are related to different members of the Onsala 1 cluster. The data presented here support a scenario in which the phases of massive star formation begin much later in the evolution of a cluster and/or UCHII region last for much longer than 10^5 yrs.

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

We have characterized the physical structure and chemical composition of two close-to-round starless cores in Taurus-Auriga, L1498 and L1517B. Our analysis is based on high angular resolution observations in at least two transitions of NH3, N2H+, CS, C34S, C18O, and C17O, together with maps of the 1.2 mm continuum. For both cores, we derive radial profiles of constant temperature and constant turbulence, together with density distributions close to those of non-singular isothermal spheres. Using these physical conditions and a Monte Carlo radiative transfer model, we derive abundance profiles for all species and model the strong chemical differentiation of the core interiors. According to our models, the NH3 abundance increases toward the core centers by a factor of several while N2H+ has a constant abundance over most of the cores. In contrast, both C18O and CS (and isotopomers) are strongly depleted in the core interiors, most likely due to their freeze out onto grains at densities of a few 10^4 cm^{-3}. Concerning the kinematics of the dense gas, we find (in addition to constant turbulence) a pattern of internal motions at the level of 0.1 km s^{-1}. These motions seem correlated with asymmetries in the pattern of molecular depletion, and we interpret them as residuals of core contraction. Their distribution and size suggest that core formation occurs in a rather irregular manner and with a time scale of a Myr. A comparison of our derived core properties with those predicted by supersonic turbulence models of core formation shows that our Taurus cores are much more quiescent than representative predictions from these models.

L1521E seems unique among starless cores. It stands out in a distribution of a ratio (R) that we define to asses core evolution, and which compares the emission of the easily-depleted C18O molecule with that of the hard to deplete, late-time species N2H+. While all cores we have studied so far have R ratio lower than 1, L1521E has an R value of 3.4, which is 8 times the mean of the other cores. To understand this difference, we have modeled the C18O and N2H+ abundance profiles in L1521E using a density distribution derived from 1.2mm continuum data. Our model shows that the C18O emission in this core is consistent with constant abundance, and this makes L1521E the first core with no C18O depletion. Our model also derives an unusually low N2H+ abundance. These two chemical peculiarities suggest that L1521E has contracted to its present density very recently, and it is therefore an extremely young starless core. Comparing our derived abundances with a chemical model, we estimate a tentative age of 1.5 x 10^5 yr, which is too short for ambipolar diffusion models.

We study the density structure of candidate contracting starless core L694-2 using 1.3 mm dust continuum observations from the IRAM Plateau de Bure Interferometer and the Berkeley-Illinois-Maryland Array, probing spatial scales from 10000-500 AU. The long baseline PdBI observations detect no emission, and limit the contamination from a compact component F_c < 2.7 mJy. This limit corresponds to a very small disk mass, M_disk < 5e-4 M_sun x (60 K / T_disk), and bolsters the ``starless'' interpretation of the core. The shorter baseline BIMA data are compared to density models using a physically motivated temperature distribution with a central minimum. This analysis provides clear evidence for a turn-over from the steep outer density profile observed in dust extinction to much more shallow behavior in the inner regions (<7500 AU). The best fit Bonnor-Ebert, Plummer-like, broken power law, and end-on cylinder models produce very similar flattened profiles and cannot be distinguished. We quantify the sensitivity of the inferred structure to various uncertainties, including the temperature distribution, the central position, and the presence of a weak unresolved central component. The largest uncertainty comes from the temperature assumption; an isothermal model modifies the best fit parameters by \~2 sigma, with the inferred profiles more shallow. Dust emission and extinction profiles are reproduced by an embedded cylinder with scale height H=13.5'' inclined at a small angle to the line of sight. The flat central density profile suggests that pressure forces still support the core, despite the extended inward motions inferred from molecular line observations (Lee, Myers, & Tafalla 2001). In the context of the cylindrical density model, these inward motions may represent the contraction of a prolate core along its major axis.

We present subarcsecond resolution dust continuum observations of the protostellar collapse candidate B335 made with the IRAM Plateau de Bure Interferometer at wavelengths of 1.2 and 3.0 mm. These observations probe to < 100 AU size scales and reveal a compact source component that we identify with a circumstellar disk. We analyze these data in concert with previous lower resolution interferometer observations and find a best fit density structure for B335 that consists of a power law envelope with index p=1.55 +/- 0.04 (r < 5000 AU) together with a disk (r < 100 AU) of flux F_{1.2 mm}=21 +/-2 mJy. We estimate a systematic uncertainty in the power law index delta(p) < 0.15, where the largest error comes from the assumed form of the dust temperature falloff with radius. This determination of the inner density structure of B335 has a precision unique amongst protostellar cores, and it is consistent with the r^{-1.5} profile of gravitational free-fall, in accord with basic expectations for the formation of a star. The flux (and implied mass) of the compact component in B335 is typical of the disks around T Tauri stars.

We present new high-sensitivity 1.3 mm bolometer observations of the young outflow L1157. These data show that the continuum emission arises from four distinct components: a circumstellar disk, a protostellar envelope, an extended flattened envelope --the dense remnant of the molecular cloud in which the protostar was formed--, and the outflow itself, which represents ~20% of the total flux. The outflow emission exhibits two peaks that are coincident with the two strong shocks in the southern lobe of L1157. We show that the mm continuum is dominated by thermal dust emission arising in the high velocity material. The spectral index derived from the new 1.3 mm data and 850 mu observations from Shirley et al. (2000), is ~5 in the outflow, significantly higher than in the protostellar envelope (~3.5). This can be explained by an important line contamination of the 850 mu map, and/or by different dust characteristics in the two regions, possibly smaller grains in the post-shocks regions of the outflow. Our observations show that bipolar outflows can present compact emission peaks which must not be misinterpreted as protostellar condensations when mapping star forming regions.

We present a study of the density structure of the protostellar collapse candidate B335 using continuum observations from the IRAM Plateau de Bure Interferometer made at wavelengths of 1.2mm and 3.0mm. We analyze these data, which probe spatial scales from 5000 AU to 500 AU, directly in the visibility domain by comparison to synthetic observations constructed from models that assume different physical conditions. This approach allows for much more stringent constraints to be derived from the data than from analysis of images. A single radial power law in density provides a good description of the data, with best fit power law index p=1.65+/-0.05. Through simulations, we quantify the sensitivity of this result to various model uncertainties, including assumptions of temperature distribution, outer boundary, dust opacity spectral index, and an unresolved central component. The largest uncertainty comes from the unknown presence of a centralized point source. A point source with 1.2mm flux of F=12+/-7 mJy reduces the density index to p=1.47+/-0.07. The remaining sources of systematic uncertainty, the most important of which is the temperature distribution, likely contribute a total uncertainty of < 0.2. We therefore find strong evidence that the power law index of the density distribution within 5000 AU is significantly less than the value at larger radii, close to 2.0 from previous studies of dust emission and extinction. These results conform well to the generic paradigm of isolated, low-mass star formation which predicts a power law density index close to p=1.5 for an inner region of gravitational free fall onto the protostar.

We present results of an extensive mapping survey of N2H+(1-0) in about 60 low mass cloud cores already mapped in the NH3(1,1) inversion transition line. The survey has been carried out at the FCRAO antenna with an angular resolution about 1.5 times finer than the previous ammonia observations. Cores with stars typically have map sizes about a factor of two smaller for N2H+ than for NH3, indicating the presence of denser and more centrally concentrated gas compared to starless cores. Significant correlations are found between NH3 and N2H+ column densities and excitation temperatures in starless cores, but not in cores with stars, suggesting a different chemical evolution of the two species. Velocity gradients range between 0.5 and 6 km/s/pc, similar to what has been found with NH3 data. ``Local'' velocity gradients show significant variation in both magnitude and direction, suggesting the presence of complexmotions not interpretable as simple solid body rotation. Integrated intensity profiles of starless cores present a ``central flattening'' and are consistent with a spherically symmetric density law n ~ r^{-1.2} for r < ~0.03 pc and n ~ r^{-2} at larger r. Cores with stars are better modelled with single density power laws with n ~ r^{-2}. Line widths change across the core but we did not find a general trend. The deviation in line width correlates with the mean line width, suggesting that the line of sight contains ~ 10 coherence lengths. The corresponding value of the coherence length, ~ 0.01 pc, is similar to the expected cutoff wavelength for MHD waves. This similarity may account for the increased ``coherence'' of line widths on small scales. Despite of the finer angular resolution, the majority of N2H+ and NH3 maps show a similar ``simple'' structure, with single peaks and no elongation.

(Abridged) We present evidence that low-mass starless cores, the simplest units of star formation, are systematically differentiated in their chemical composition. Molecules including CO and CS almost vanish near the core centers, where the abundance decreases by one or two orders of magnitude. At the same time, N2H+ has a constant abundance, and the fraction of NH3 increases toward the core center. Our conclusions are based on a study of 5 mostly-round starless cores (L1498, L1495, L1400K, L1517B, and L1544), which we have mappedin C18O(1-0), C17O(1-0), CS(2-1), C34S(2-1), N2H+(1-0), NH3(1,1) and (2,2), and the 1.2 mm continuum. For each core we have built a model that fits simultaneously the radial profile of all observed emission and the central spectrum for the molecular lines. The observed abundance drops of CO and CS are naturally explained by the depletion of these molecules onto dust grains at densities of 2-6 10^4 cm-3. N2H+ seems unaffected by this process up to densities of several 10^5, while the NH3 abundance may be enhanced by reactions triggered by the disappearance of CO from the gas phase. With the help of our models, we show that chemical differentiation automatically explains the discrepancy between the sizes of CS and NH3 maps, a problem which has remained unexplained for more than a decade. Our models, in addition, show that a combination of radiative transfer effects can give rise to the previously observed discrepancy in the linewidth of these two tracers. Although this discrepancy has been traditionally interpreted as resulting from a systematic increase of the turbulent linewidth with radius, our models show that it can arise in conditions of constant gas turbulence.

We have mapped the dense dark core L1544 in H13CO+(1-0), DCO+(2-1), DCO+(3-2), N2H+(1-0), NTH+(3-2), N2D+(2-1), N2D+(3-2), C18O(1-0), and C17O(1-0) using the IRAM 30-m telescope. We have obtained supplementary observations of HC18O+(1-0), HC17O+(1-0), and D13CO+(2-1). Many of the observed maps show a general correlation with the distribution of dust continuum emission in contrast to C18O(1-0) and C17O(1-0) which give clear evidence for depletion of CO at positions close to the continuum peak. In particular N2D+(2-1) and (3-2) and to a lesser extent N2H+(1-0) appear to be excellent tracers of the dust continuum. We find that the tracers of high density gas (in particular N2D+) show a velocity gradient along the minor axis of the L1544 core and that there is evidence for larger linewidths close to the dust emission peak. We interpret this using the model of the L1544 proposed by Ciolek & Basu (2000) and by comparing the observed velocities with those expected on the basis of their model. The results show reasonable agreement between observations and model in that the velocity gradient along the minor axis and the line broadening toward the center of L1544 are predicted by the model. This is evidence in favour of the idea that amipolar diffusion across field lines is one of the basic processes leading to gravitational collapse. However, line widths are significantly narrower than observed and are better reproduced by the Myers & Zweibel (2001) model which considers the quasistatic vertical contraction of a layer due to dissipation of its Alfvenic turbulence, indicating the importance of this process for cores in the verge of forming a star.

The maps presented in Paper I are here used to infer the variation of the column densities of HCO+, DCO+, N2H+, and N2D+ as a function of distance from the dust peak. These results are interpreted with the aid of a crude chemical model which predicts the abundances of these species as a function of radius in a spherically symmetric model with radial density distribution inferred from the observations of dust emission at millimeter wavelengths and dust absorption in the infrared. Our main observational finding is that the N(N2D+)/N(N2H+) column density ratio is of order 0.2 towards the L1544 dust peak as compared to N(DCO+)/N(HCO+) = 0.04. We conclude that this result as well as the general finding that N2H+ and N2D+ correlate well with the dust is caused by CO being depleted to a much higher degree than molecular nitrogen in the high density core of L1544. Depletion also favors deuterium enhancement and thus N2D+, which traces the dense and highly CO-depleted core nucleus, is much more enhanced than DCO+. Our models do not uniquely define the chemistry in the high density depleted nucleus of L1544 but they do suggest that the ionization degree is a few times 10^{-9} and that the ambipolar diffusion time scale is locally similar to the free fall time. It seems likely that the lower limit which one obtains to ionization degree by summing all observable molecular ions is not a great underestimate of the true ionization degree. We predict that atomic oxygen is abundant in the dense core and, if so, H3O+ may be the main ion in the central highly depleted region of the core.

We present the results of an extensive mapping survey of 53 `starless' cores in the optically thick line of CS 2-1 and the optically thin lines of N2H+ 1-0 and C18O 1-0. The purpose of this survey was to search for signatures of extended inward motions. This study finds 10 `strong' and 9 `probable' infall candidates, based on $\delta V_{CS}$ analysis and on the spectral shapes of CS lines. From our analysis of the blue-skewed CS spectra and the $\delta V_{CS}$ parameter, we find typical infall radii of 0.06-0.14 pc. Also, using a simple two layer radiative transfer model to fit the profiles, we derive one-dimensional infall speeds, half of whose values lie in the range of 0.05-0.09 km s$^{-1}$. These values are similar to those found in L1544 by Tafalla et al., and this result confirms that infall speeds in starless cores are generally faster than expected from ambipolar diffusion in a strongly sub-critical core. In addition, the observed infall regions are too extended to be consistent with the `inside-out' collapse model applied to a very low-mass star. In the largest cores, the spatial extent of the CS spectra with infall asymmetry is larger than the extent of the $\rm N_2H^+$ core by a factor of 2-3. All these results suggest that extended inward motions are a common feature in starless cores, and that they could represent a necessary stage in the condensation of a star-forming dense core.

New, high-sensitivity interferometric CO J=2-1 observations of the HH 7-11 outflow show that despite previous doubts, this system is powered by the Class I source SVS 13. The molecular outflow from SVS 13 is formed by a shell with a large opening angle at the base, which is typical of outflows from Class I sources, but it also contains an extremely-high-velocity jet composed of ``molecular bullets'', which is more typical of Class 0 outflows. This suggests that SVS 13 could be a very young Class I, which still keeps some features of the previous evolutionary stage. We briefly discuss the nature of some sources in the SVS 13 vicinity which are emitters of cm-wave continuum, but have no counterpart at mm wavelengths.

We present molecular-line observations toward the dense core in L483 and its bipolar outflow powered by the Class 0 object IRAS 18148-0440. CO (carbon monoxide) maps show that the outflow is well collimated and asymmetric, and that its gas is warmer than the surrounding cloud by at least a factor of 2. H2CO (formaldehyde) lines toward the outflow show prominent high-velocity wings and evidence for an H2CO abundance enhancement of a factor of 20. At ambient velocities, these lines show strong self-absorption and a brighter blue peak, a characteristic signature of inward motions. Finally, and in contrast with the outflows from other Class 0 objects, the CH3OH (methanol) lines in L483 do not show high velocity wings, and no evidence for abundance enhancement is found in this molecule or in SiO (silicon monoxide). Comparing the physical and chemical properties of the outflow in L483 with those of other outflows from Class 0 and Class I sources, we find that the L483 outflow is somewhat intermediate between these types. This suggests that the L483 central source has already started its transition between Class 0 and Class I, and that its mixed properties illustrate how this transition occurs.

In this paper, we present the first results of a survey of 220 starless cores selected primarily by their optical obscuration (Lee & Myers 1999) and observed in CS (2-1), N$_2$H$^+$ (1-0), and C$^{18}$O (1-0) using the NEROC Haystack 37-m telescope, to quantify the incidence of infall motions toward starless cores and find infall candidates for further study in the future.

We report observations of 47 candidate protostars in two optically thick lines (H2CO 212-111 and CS 2-1) and one optically thin line (N2H+ 1-0) using the IRAM 30-m, SEST 15-m, and Haystack 37-m radio telescopes. The sources were selected from the redness of their spectra (Tbol < 200 K) and their near distance (d < 400 pc). Most of the sources have asymmetric optically thick lines. The observed distribution of velocity differences, deltaV = (Vthick-Vthin)/FWHMthin, is skewed toward negative (blue-shifted) velocities for both the H2CO and CS samples. This excess is much more significant for Class 0 than for Class I sources, suggesting that we detect infall motions toward Class 0 and not toward Class I sources. This indicates a difference in the physical conditions in the circumstellar envelopes around Class I and Class 0 sources, but does not rule out the presence of infall onto Class I sources by e.g. lower opacity gas. Bipolar outflows alone, or rotation alone, cannot reproduce these statistics if the sample of sources has randomly oriented symmetry axes. We identify 15 spectroscopic infall candidates, of which 6 are new. Most of these infall candidates have primarily turbulent rather than thermal motions, and are associated with clusters rather than being isolated.