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.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

The fastest molecular component to the protostellar outflows -- extremely high-velocity (EHV) molecular jets -- are still puzzling since they are seen only rarely. The first aim is to analyze the interaction between the EHV jet and the slow outflow by comparing their outflow force content. The second aim is to analyze the chemical composition of the different outflow velocity components and to reveal the spatial location of molecules. ALMA 3 mm and 1.3 mm observations of five outflow sources at 130 -- 260 au resolution in the Serpens Main cloud are presented. Observations of CO, SiO, H$_2$CO and HCN reveal the kinematic and chemical structure of those flows. Three velocity components are distinguished: the slow and the fast wing, and the EHV jet. Out of five sources, three have the EHV component. Comparison of outflow forces reveals that only the EHV jet in the youngest source Ser-emb 8 (N) has enough momentum to power the slow outflow. The SiO abundance is generally enhanced with velocity, while HCN is present in the slow and the fast wing, but disappears in the EHV jet. For Ser-emb 8 (N), HCN and SiO show a bow-shock shaped structure surrounding one of the EHV peaks suggesting sideways ejection creating secondary shocks upon interaction with the surroundings. Also, the SiO abundance in the EHV gas decreases with distance from this protostar, whereas that in the fast wing increases. H$_2$CO is mostly associated with low-velocity gas but also appears surprisingly in one of the bullets in the Ser-emb~8~(N) EHV jet. The high detection rate suggests that the presence of the EHV jet may be more common than previously expected. The origin and temporal evolution of the abundances of SiO, HCN and H$_2$CO through high-temperature chemistry are discussed. The data are consistent with a low C/O ratio in the EHV gas versus high C/O ratio in the fast and slow wings.

The deuterium fractionation of water can serve as a tracer for the chemical and physical evolution of water during star formation and can constrain the origin of water in Solar System bodies. We determine the HDO/H$_2$O ratio in the inner warm gas toward three low-mass Class 0 protostars selected to be in isolated cores, i.e., not associated with any cloud complexes. Previous sources for which the HDO/H$_2$O ratio have been established were all part of larger star-forming complexes. Targeting these isolated protostars allows comparison of the water chemistry in isolated and clustered regions to determine the influence of local cloud environment. We present ALMA observations of the HDO $3_{1,2}$-$2_{2,1}$ and $2_{1,1}$-$2_{1,2}$ transitions at 225.897 GHz and 241.562 GHz along with the H$_2^{18}$O $3_{1,3}$-$2_{2,0}$ transition at 203.407 GHz. The high angular resolution (0\farcs3-1\farcs3) allow the study of the inner warm envelope gas. Model-independent estimates for the HDO/H$_2$O ratios are obtained and compared with previous determinations in the warm gas toward low-mass protostars. We detect the targeted water transitions toward the three sources with S/N > 5. We determine the HDO/H$_2$O ratio toward L483, B335 and BHR71-IRS1 to be ($2.2\pm0.4$)$\times 10^{-3}$, ($1.7\pm0.3$)$\times 10^{-3}$, and ($1.8\pm0.4$)$\times 10^{-3}$, respectively, assuming $T_\mathrm{ex} = 124$ K. The degree of water deuteration of these isolated protostars are a factor of 2-4 higher relative to Class 0 protostars that are members of known nearby clustered star-forming regions. The results indicate that the water deuterium fractionation is influenced by the local cloud environment. This effect can be explained by variations in either collapse timescales or temperatures, which depends on local cloud dynamics and could provide a new method to decipher the history of young stars.

Context. The relationship between outflow launching and formation of accretion disks around young stellar objects is still not entirely understood, which is why spectrally and spatially resolved observations are needed. Recently, the Atacama Large Millimetre/sub-millimetre Array (ALMA) has carried out long-baseline observations towards a handful of sources, revealing connections between outflows and the inner regions of disks. Aims. Here we aim to determine the small-scale kinematic and morphological properties of the outflow from the isolated protostar B335 for which no Keplerian disk has, so far, been observed on scales down to 10 au. Methods. We use ALMA in its longest-baseline configuration to observe emission from CO isotopologs, SiO, SO$_2$ and CH$_3$OH. The proximity of B335 provides a resolution of ~3 au (0.03''). We also combine our long-baseline data with archival data to produce a high-fidelity image covering scales up to 700 au (7''). Results. $^{12}$CO has a X-shaped morphology with arms ~50 au in width that we associate with the walls of an outflow cavity, similar to what is observed on larger scales. Long-baseline continuum emission is confined to <7 au of the protostar, while short-baseline continuum emission follows the $^{12}$CO outflow and cavity walls. Methanol is detected within ~30 au of the protostar. SiO is also detected in the vicinity of the protostar, but extended along the outflow. Conclusions. The $^{12}$CO outflow shows no clear signs of rotation at distances $\gtrsim$30 au from the protostar. SiO traces the protostellar jet on small scales, but without obvious rotation. CH$_3$OH and SO$_2$ trace a region <16 au in diameter, centred on the continuum peak, which is clearly rotating. Using episodic, high-velocity, $^{12}$CO features, we estimate the launching radius of the outflow to be <0.1 au and dynamical timescales on the order of a few years.

With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 {\mu}m, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in Serpens Emb8(N). We use 3 mm and 1 mm molecular tracers to investigate outflow and dense gas properties and their correlation with the polarization. These observations allow us to investigate the physical processes involved in the Radiative Alignment Torques (RATs) acting on dust grains along the outflow cavity walls, which experience irradiation from accretion processes and outflow shocks. The inner core of SMM1-a presents a polarization pattern with a poloidal magnetic field at the bases of the two lobes of the bipolar outflow. To the south of SMM1-a we see two polarized filaments, one of which seems to trace the redshifted outflow cavity wall. The other may be an accretion streamer of material infalling onto the central protostar. We propose that the polarized emission we see at millimeter wavelengths along the irradiated cavity walls can be reconciled with the expectations of RAT theory if the aligned grains present at '<' 500 au scales in Class 0 envelopes have grown larger than the 0.1 {\mu}m size of ISM dust grains. Our observations allow us to constrain the star-forming sources magnetic field morphologies within the central cores, along the outflow cavity walls, and in possible accretion streamers.

Observations with modern instruments such as Herschel reveal that stars form clustered inside filamentary arms of ~1 pc length embedded in Giant Molecular Clouds. On smaller scales (~1000 au), observations of, e.g., IRAS 16293--2422 show signs of filamentary `bridge' structures connecting young protostars in their birth environment. We investigate the formation of a deeply embedded triple protostellar system, and characterize its connection to the filamentary structure present on GMC scales. Using the magnetohydrodynamical code RAMSES, we carry out zoom-in simulations of low-mass star formation starting from Giant-Molecular-Cloud-scales. Colliding flows of gas in the filamentary arms induce the formation of two protostellar companions at distances of ~1000 au from the primary. After their birth, the stellar companions quickly (~10 kyr) approach and orbit the primary on eccentric orbits with separations of ~100 au. The colliding flows induce transient structures lasting for up to a few 10 kyr connecting two forming protostellar objects that are kinematically quiescent along the line-of-sight. Colliding flows compress the gas and trigger the formation of stellar companions via turbulent fragmentation. Our results suggest that protostellar companions form with a wide separation of ~1000 au. Smaller separations (~100 au or less) are a consequence of subsequent migration and capturing. Associated with the formation phase of the companion, the turbulent environment induces arc- and bridge-like structures. These bridges can become kinematically quiescent, when the velocity components of the colliding flows eliminate each other. However, the gas in bridges still contributes to stellar accretion later. Our results demonstrate: bridge-like structures are a transient phenomenon of stellar multiple formation.

(Abridged) The purpose of this paper is to explore the structure of a line-rich Class I protobinary source, Oph-IRS 67, and analyse the differences and similarities with Class 0 and Class II sources. We present a systematic molecular line study of IRS 67 with the Submillimeter Array (SMA) on 1 - 2" (150 - 300 AU) scales. The wide instantaneous band-width of the SMA observations (~30 GHz) provide detections of a range of molecular transitions that trace different physics, such as CO isotopologues, sulphur-bearing species, deuterated species, and carbon-chain molecules. We see significant differences between different groups of species. For example, the CO isotopologues and sulphur-bearing species show a rotational profile and are tracing the larger-scale circumbinary disc structure, while CN, DCN, and carbon-chain molecules peak at the southern edge of the disc at blue-shifted velocities. In addition, the cold gas tracer DCO+ is seen beyond the extent of the circumbinary disc. The detected molecular transitions can be grouped into three main components: cold regions far from the system, the circumbinary disc, and a UV-irradiated region likely associated with the surface layers of the disc that are reached by the UV radiation from the sources. The different components are consistent with the temperature structure derived from the ratio of two H2CO transitions, that is, warm temperatures are seen towards the outflow direction, lukewarm temperatures are associated with the UV-radiated region, and cold temperatures are related with the circumbinary disc structure. The chemistry towards IRS 67 shares similarities with both Class 0 and Class II sources, possibly due to the high gas column density and the strong UV radiation arising from the binary system. IRS 67 is, therefore, highlighting the intermediate chemistry between deeply embedded sources and T-Tauri discs.

(Abridged) The purpose of this paper is to explore and compare the physical and chemical structure of Class I low-mass protostellar sources on protoplanetary disc scales. We present a study of the dust and gas emission towards a representative sample of 12 Class I protostars from the Ophiuchus molecular cloud with the Atacama Large Millimeter/submillimeter Array (ALMA). The continuum at 0.87 mm and molecular transitions from C17O, C34S, H13CO+, CH3OH, SO2 , and C2H were observed at high angular resolution (0.4", ~60 au diameter) towards each source. Disc and stellar masses are estimated from the continuum flux and position-velocity diagrams, and six of the sources show disc-like structures. Towards the more luminous sources, compact emission and large line widths are seen for transitions of SO2 that probe warm gas (Eu ~200 K). In contrast, C17O emission is detected towards the least evolved and less luminous systems. No emission of CH3OH is detected towards any of the continuum peaks, indicating an absence of warm CH3OH gas towards these sources. A power-law relation is seen between the stellar mass and the bolometric luminosity, corresponding to a mass accretion rate of (2.4 +/- 0.6) x 10^-7 Msun/year for the Class I sources. This mass accretion rate is lower than the expected value if the accretion is constant in time and rather points to a scenario of accretion occurring in bursts. The differentiation between C17O and SO2 suggests that they trace different physical components: C17O traces the densest and colder regions of the disc-envelope system, while SO2 may be associated with regions of higher temperature, such as accretion shocks. The lack of warm CH3OH emission suggests that there is no hot-core-like region around any of the sources and that the CH3OH column density averaged over the disc is low.

[Abridged] The majority of stars form in binary or higher order systems. The Class 0 protostellar system IRAS16293-2422 contains two protostars, 'A' and 'B', separated by ~600 au and embedded in a single, 10^4 au scale envelope. Their relative evolutionary stages have been debated. We aim to study the relation and interplay between the two protostars A and B at spatial scales of 60 to ~1000 au. We selected molecular gas line transitions of CO, H2CO, HCN, CS, SiO, and CCH from the ALMA-PILS spectral imaging survey (329-363 GHz) and used them as tracers of kinematics, density, and temperature in the IRAS16293-2422 system. The angular resolution of the PILS data set allows us to study these quantities at a resolution of 0.5 arcsec (60 au [..]). Line-of-sight velocity maps of both optically thick and optically thin molecular lines reveal: (i) new manifestations of previously known outflows emanating from protostar A; (ii) a kinematically quiescent bridge of dust and gas spanning between the two protostars, with an inferred density between 4 10^4 and 3 10^7 cm^-3; and (iii) a separate, straight filament seemingly connected to protostar B seen only in CCH, with a flat kinematic signature. Signs of various outflows, all emanating from source A, are evidence of high-density and warmer gas; none of them coincide spatially and kinematically with the bridge. We hypothesize that the bridge arc is a remnant of filamentary substructure in the protostellar envelope material from which protostellar sources A and B have formed. One particular morphological structure appears to be due to outflowing gas impacting the quiescent bridge material. The continuing lack of clear outflow signatures unambiguously associated to protostar B and the vertically extended shape derived for its disk-like structure lead us to conclude that source B may be in an earlier evolutionary stage than source A.

Studies of deuterated isotopologues of complex organic molecules can provide important constraints on their origin in regions of star formation. In particular, the abundances of deuterated species are very sensitive to the physical conditions in the environment where they form. Due to the low temperatures in regions of star formation, these isotopologues are enhanced to significant levels, making detections of multiply-deuterated species possible. However, for complex organic species, only the multiply-deuterated variants of methanol and methyl cyanide have been reported so far. The aim of this paper is to initiate the characterisation of multiply-deuterated variants of complex organic species with the first detection of doubly-deuterated methyl formate, CHD2OCHO. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) of the protostellar binary IRAS 16293-2422, in the spectral range of 329.1 GHz to 362.9 GHz. We report the first detection of doubly-deuterated methyl formate CHD2OCHO in the ISM. The D/H ratio of CHD2OCHO is found to be 2-3 times higher than the D/H ratio of CH2DOCHO for both sources, similar to the results for formaldehyde from the same dataset. The observations are compared to a gas-grain chemical network coupled to a dynamical physical model, tracing the evolution of a molecular cloud until the end of the Class 0 protostellar stage. The overall D/H ratio enhancements found in the observations are of the same order of magnitude as the predictions from the model for the early stages of Class 0 protostars. However, the higher D/H ratio of CHD2OCHO compared to the D/H ratio of CH2DOCHO is still not predicted by the model. This suggests that a mechanism is enhancing the D/H ratio of singly- and doubly-deuterated methyl formate that is not in the model, e.g. mechanisms for H-D substitutions.

We present a characterization of the binary protostar system that is forming within a dense core in the isolated dark cloud BHR71. The pair of protostars, IRS1 and IRS2, are both in the Class 0 phase, determined from observations that resolve the sources from 1 um out to 250 um and from 1.3 mm to 1.3cm. The resolved observations enable the luminosities of IRS1 and IRS2 to be independently measured (14.7 and 1.7L_sun, respectively), in addition to the bolometric temperatures 68~K, and 38~K, respectively. The surrounding core was mapped in NH3 (1,1) with the Parkes radio telescope, and followed with higher-resolution observations from ATCA in NH3 (1,1) and 1.3cm continuum. The protostars were then further characterized with ALMA observations in the 1.3~mm continuum along with N2D+ (J=3-2), 12CO, 13CO, and C18O (J=2-1) molecular lines. The Parkes observations find evidence for a velocity gradient across the core surrounding the two protostars, while ATCA reveals more complex velocity structure toward the protostars within the large-scale gradient. The ALMA observations then reveal that the two protostars are at the same velocity in C18O, and N2H+ exhibits a similar velocity structure as NH3. However, the C18O kinematics reveal that the rotation on scales $<$1000~AU around IRS1 and IRS2 are in opposite directions. Taken with the lack of a systematic velocity difference between the pair, it is unlikely that their formation resulted from rotational fragmentation. We instead conclude that the binary system most likely formed via turbulent fragmentation of the core.

Hydroxylamine (NH$_{2}$OH) and methylamine (CH$_{3}$NH$_{2}$) have both been suggested as precursors to the formation of amino acids and are therefore of interest to prebiotic chemistry. Their presence in interstellar space and formation mechanisms, however, are not well established. We aim to detect both amines and their potential precursor molecules NO, N$_{2}$O and CH$_{2}$NH towards the low-mass protostellar binary IRAS 16293--2422, in order to investigate their presence and constrain their interstellar formation mechanisms around a young Sun-like protostar. ALMA observations from the unbiased, high angular resolution and sensitivity Protostellar Interferometric Line Survey (PILS) are used. Spectral transitions of the molecules under investigation are searched for with the CASSIS line analysis software. CH$_2$NH and N$_{2}$O are detected for the first time towards a low-mass source, the latter molecule through confirmation with the single-dish TIMASSS survey. NO is also detected. CH$_{3}$NH$_{2}$ and NH$_{2}$OH are not detected and stringent upper limit column densities are determined. The non-detection of CH$_{3}$NH$_{2}$ and NH$_{2}$OH limits the importance of formation routes to amino acids involving these species. The detection of CH$_{2}$NH makes amino acid formation routes starting from this molecule plausible. The low abundances of CH$_2$NH and CH$_{3}$NH$_{2}$ compared to Sgr B2 indicate that different physical conditions influence their formation in low- and high-mass sources.

(Abridged) Protostellar systems evolve from prestellar cores, through the deeply embedded stage and then disk-dominated stage, before they end up on the main sequence. Knowing how much time a system spends in each stage is crucial for understanding how stars and associated planetary systems form, because a key constraint is the time available to form such systems. Equally important is understanding what the spread in these time scales is. The most commonly used method for inferring protostellar ages is to assume the lifetime of one evolutionary stage, and then scale this to the relative number of protostars in the other stages, i.e., assuming steady state. This method does not account for the underlying age distribution and apparent stochasticity of star formation, nor that relative populations are not in steady state. To overcome this, we propose a new scheme where the lifetime of each protostellar stage follows a distribution based on the formalism of sequential nuclear decay. The main assumptions are: Class 0 sources follow a straight path to Class III sources, the age distribution follows a binomial distribution, and the star-formation rate is constant. The results are that the half-life of Class 0, Class I, and Flat sources are (2.4+/-0.2)%, (4.4+/-0.3)%, and (4.3+/-0.4)% of the Class II half-life, respectively, which translates to 47+/-4, 88+/-7, and 87+/-8 kyr, respectively, for a Class II half-life of 2 Myr for protostars in the Gould Belt clouds with more than 100 protostars. The mean age of these clouds is 1.2+/-0.1 Myr, and the star formation rate is (8.3+/-0.5)x10^-4 Msun/yr. The critical parameters in arriving at these numbers are the assumed half-life of the Class II stage, and the assumption that the star-formation rate and half-lives are constant. This method presents a first step in moving from steady-state to non-steady-state solutions of protostellar populations.

Methyl isocyanide (CH$_3$NC) is the isocyanide with the largest number of atoms confirmed in the interstellar medium (ISM), but it is not an abundant molecule, having only been detected towards a handful of objects. Conversely, its isomer, methyl cyanide (CH$_3$CN), is one of the most abundant complex organic molecules detected in the ISM, with detections in a variety of low- and high-mass sources. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) to search for methyl isocyanide and compare its abundance with that of its isomer methyl cyanide. We use a new line catalogue from the Cologne Database for Molecular Spectroscopy (CDMS) to identify methyl isocyanide lines. We also model the chemistry with an updated version of the three-phase chemical kinetics model {\em MAGICKAL}, presenting the first chemical modelling of methyl isocyanide to date. We detect methyl isocyanide for the first time in a solar-type protostar, IRAS 16293$-$2422 B, and present upper limits for its companion protostar, IRAS 16293$-$2422 A. Methyl isocyanide is found to be at least 20 times more abundant in source B compared to source A, with a CH$_3$CN/CH$_3$NC abundance ratio of 200 in IRAS 16293--2422 B and >5517 in IRAS 16293$-$2422 A. We also present the results of a chemical model of methyl isocyanide chemistry in both sources, and discuss the implications on methyl isocyanide formation mechanisms and the relative evolutionary stages of both sources.

The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary disks are more massive than their mature counterparts, this suggests that planet formation starts early, but evidence for grain growth that is spatially and temporally coincident with a massive reservoir in young disks remains scarce. Here, we report observations on a lack of emission of carbon monoxide isotopologues within the inner ~15 au of a very young (age ~100,000 years) disk around the Solar-type protostar TMC1A. By using the absence of spatially resolved molecular line emission to infer the gas and dust content of the disk, we conclude that shielding by millimeter-size grains is responsible for the lack of emission. This suggests that grain growth and millimeter-size dust grains can be spatially and temporally coincident with a mass reservoir sufficient for giant planet formation. Hence, planet formation starts during the earliest, embedded phases in the life of young stars.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

We present full spectral scans from 200-670$\mu$m of 26 Class 0+I protostellar sources, obtained with $Herschel$-SPIRE, as part of the "COPS-SPIRE" Open Time program, complementary to the DIGIT and WISH Key programs. Based on our nearly continuous, line-free spectra from 200-670 $\mu$m, the calculated bolometric luminosities ($L_{\rm bol}$) increase by 50% on average, and the bolometric temperatures ($T_{\rm bol}$) decrease by 10% on average, in comparison with the measurements without Herschel. Fifteen protostars have the same Class using $T_{\rm bol}$ and $L_{\rm bol}$/$L_{\rm submm}$. We identify rotational transitions of CO lines from J=4-3 to J=13-12, along with emission lines of $^{13}$CO, HCO$^+$, H$_{2}$O, and [CI]. The ratios of $^{12}$CO to $^{13}$CO indicate that $^{12}$CO emission remains optically thick for $J_{\rm up}$ < 13. We fit up to four components of temperature from the rotational diagram with flexible break points to separate the components. The distribution of rotational temperatures shows a primary population around 100 K with a secondary population at $\sim$350 K. We quantify the correlations of each line pair found in our dataset, and find the strength of correlation of CO lines decreases as the difference between $J$-level between two CO lines increases. The multiple origins of CO emission previously revealed by velocity-resolved profiles are consistent with this smooth distribution if each physical component contributes to a wide range of CO lines with significant overlap in the CO ladder. We investigate the spatial extent of CO emission and find that the morphology is more centrally peaked and less bipolar at high-$J$ lines. We find the CO emission observed with SPIRE related to outflows, which consists two components, the entrained gas and shocked gas, as revealed by our rotational diagram analysis as well as the studies with velocity-resolved CO emission.

Complex organic molecules are readily detected in the inner regions of the gaseous envelopes of forming protostars. In particular, molecules that contain nitrogen are interesting due to the role nitrogen plays in the development of life and the compact scales such molecules have been found to trace around forming protostars. The goal of this work is to determine the inventory of one family of nitrogen-bearing organic molecules, complex nitriles (molecules with a $-$CN functional group) towards two hot corino sources in the low-mass protostellar binary IRAS 16293$-$2422. This work explores the abundance differences between the two sources, the isotopic ratios, and the spatial extent derived from molecules containing the nitrile functional group. Using data from the Protostellar Interferometric Line Survey (PILS) obtained with ALMA we determine abundances and excitation temperatures for the detected nitriles. We also present a new method for determining the spatial structure of sources with high line density and large velocity gradients $-$ Velocity-corrected INtegrated emission (VINE) maps. We detect methyl cyanide (CH$_3$CN) as well as 5 of its isotopologues, including the detection of CHD$_2$CN which is the first detection in the ISM. We also detect ethyl cyanide (C$_2$H$_5$CN), vinyl cyanide (C$_2$H$_3$CN), and cyanoacetylene (HC$_3$N). We find that abundances are similar between IRAS 16293A and IRAS 16293B on small scales except for vinyl cyanide which is only detected towards the latter source. This suggests an important difference between the sources either in their evolutionary stage or warm-up timescales. We also detect a spatially double-peaked emission for the first time in molecular emission in the A source, suggesting that this source is showing structure related to a rotating toroid of material.

(Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C17O, H13CO+ , C34S, SO2, C2H and c-C3H2 molecular transitions with ALMA at 0.8 mm. The spectrally and spatially resolved maps reveal the kinematics of the circumbinary disk as well as its chemistry. Molecular abundances are estimated using RADEX. The continuum emission reveals the presence of a circumbinary disk around the two sources. This disk has a diameter of ~620 AU and is well traced by C17O and H13CO+ emission. C2H and c-C3H2 trace a higher-density region which is spatially offset from the sources (~430 AU). Finally, SO2 shows compact emission around one of the sources, Oph-IRS67 B. The molecular transitions which trace the circumbinary disk are consistent with a Keplerian profile on disk scales (< 200 AU) and an infalling profile for envelope scales (> 200 AU). The Keplerian fit leads to a mass of 2.2 Msun. Inferred CO abundances w.r.t. H2 are comparable to the canonical ISM value of 2.7e-4. This study proves the first detection of the circumbinary disk associated with Oph-IRS67. The disk is chemically differentiated from the nearby high-density region. The lack of methanol emission suggests the extended disk dominates the mass budget in the inner- most regions of the protostellar envelope, generating a flat density profile where less material is exposed to high temperatures. Thus, complex organic molecules would be associated with lower column densities. Finally, Oph-IRS67 is a promising candidate for the detection of both circumstellar disks with higher-angular-resolution observations.

Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars (Class 0/I). The Herschel-PACS spectral maps covering $\sim55-210$ $\mu$m with a field-of-view of $\sim$50'' are used to quantify the gas excitation conditions and spatial extent using rotational transitions of H$_{2}$O, high-$J$ CO, and OH, as well as [O I] and [C II]. We confirm that a warm ($\sim$300 K) CO reservoir is ubiquitous and that a hotter component ($760\pm170$ K) is frequently detected around protostars. The line emission is extended beyond $\sim$1000 AU spatial scales in 40/90 objects, typically in molecular tracers in Class 0 and atomic tracers in Class I objects. High-velocity emission ($\gtrsim90$ km s$^{-1}$) is detected in only 10 sources in the [O I] line, suggesting that the bulk of [O I] arises from gas that is moving slower than typical jets. Line flux ratios show an excellent agreement with models of $C$-shocks illuminated by UV photons for pre-shock densities of $\sim$$10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-IR molecular and atomic lines are a unique diagnostic of feedback from UV emission and shocks in envelopes of deeply embedded protostars.

Context. The Class 0 protostellar binary IRAS 16293-2422 is an interesting target for (sub)millimeter observations due to, both, the rich chemistry toward the two main components of the binary and its complex morphology. Its proximity to Earth allows the study of its physical and chemical structure on solar system scales using high angular resolution observations. Such data reveal a complex morphology that cannot be accounted for in traditional, spherical 1D models of the envelope. Aims. The purpose of this paper is to study the environment of the two components of the binary through 3D radiative transfer modeling and to compare with data from the Atacama Large Millimeter/submillimeter Array. Such comparisons can be used to constrain the protoplanetary disk structures, the luminosities of the two components of the binary and the chemistry of simple species. Methods. We present 13CO, C17O and C18O J=3-2 observations from the ALMA Protostellar Interferometric Line Survey (PILS), together with a qualitative study of the dust and gas density distribution of IRAS 16293-2422. A 3D dust and gas model including disks and a dust filament between the two protostars is constructed which qualitatively reproduces the dust continuum and gas line emission. Results and conclusions. Radiative transfer modeling of source A and B, with the density solution of an infalling, rotating collapse or a protoplanetary disk model, can match the constraints for the disk-like emission around source A and B from the observed dust continuum and CO isotopologue gas emission. If a protoplanetary disk model is used around source B, it has to have an unusually high scale-height in order to reach the dust continuum peak emission value, while fulfilling the other observational constraints. Our 3D model requires source A to be much more luminous than source B; LA ~ 18 $L_\odot$ and LB ~ 3 $L_\odot$.

Far-infrared spectroscopy reveals gas cooling and its underlying heating due to physical processes taking place in the surroundings of protostars. These processes are reflected in both the chemistry and excitation of abundant molecular species. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars from the WISH (van Dishoeck et al. 2011), DIGIT (Green et al. 2013), and WILL surveys (Mottram et al. 2017). The $5\times5$ spectra covering the $\sim50''\times50''$ field-of-view include rotational transitions of CO, H$_2$O, and OH lines, as well as fine-structure [O I] and [C II] in the $\sim$50-200 $\mu$m range. The CO rotational temperatures (for $J_\mathrm{u}\geq14)$ are typically $\sim$300 K, with some sources showing additional components with temperatures as high as $\sim$1000 K. The H$_2$O / CO and H$_2$O / OH flux ratios are low compared to stationary shock models, suggesting that UV photons may dissociate some H$_2$O and decrease its abundance. Comparison to C shock models illuminated by UV photons shows good agreement between the line emission and the models for pre-shock densities of $10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-infrared molecular and atomic lines are the unique diagnostic of shocks and UV fields in deeply-embedded sources.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (~0.1 pc) scales -- where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded -- and the intermediate and small scales probed by CARMA (~1000 AU resolution), the SMA (~350 AU resolution), and ALMA (~140 AU resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (~130 AU) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both CO(2-1) and SiO(5-4); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1.

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.

In molecular outflows from forming low-mass protostars, most oxygen is expected to be locked up in water. However, Herschel observations have shown that typically an order of magnitude or more of the oxygen is still unaccounted for. To test if the oxygen is instead in atomic form, SOFIA-GREAT observed the R1 position of the bright molecular outflow from NGC1333-IRAS4A. The [OI] 63 um line is detected and spectrally resolved. From an intensity peak at +15 km/s, the intensity decreases until +50 km/s. The profile is similar to that of high-velocity (HV) H2O and CO 16-15, the latter observed simultaneously with [OI]. A radiative transfer analysis suggests that ~15% of the oxygen is in atomic form toward this shock position. The CO abundance is inferred to be ~10^-4 by a similar analysis, suggesting that this is the dominant oxygen carrier in the HV component. These results demonstrate that a large portion of the observed [OI] emission is part of the outflow. Further observations are required to verify whether this is a general trend.

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar.

Context: Understanding how accretion proceeds is a key question of star formation, with important implications for both the physical and chemical evolution of young stellar objects. In particular, very little is known about the accretion variability in the earliest stages of star formation. Aims: To characterise protostellar accretion histories towards individual sources by utilising sublimation and freeze-out chemistry of CO. Methods: A sample of 24 embedded protostars are observed with the Submillimeter Array (SMA) in context of the large program "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). The size of the C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is measured towards each source and compared to the expected size of the region given the current luminosity. The SMA observations also include 1.3 mm continuum data, which are used to investigate whether a link can be established between accretion bursts and massive circumstellar disks. Results: Depending on the adopted sublimation temperature of the CO ice, between 20% and 50% of the sources in the sample show extended C$^{18}$O emission indicating that the gas was warm enough in the past that CO sublimated and is currently in the process of refreezing; something which we attribute to a recent accretion burst. Given the fraction of sources with extended C$^{18}$O emission, we estimate an average interval between bursts of 20000-50000 yr, which is consistent with previous estimates. No clear link can be established between the presence of circumstellar disks and accretion bursts, however the three closest known binaries in the sample (projected separations <20 AU) all show evidence of a past accretion burst, indicating that close binary interactions may also play a role in inducing accretion variability.

High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles.

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

We report ALMA observations of a one-sided, high-velocity ($\sim$80 km s$^{-1}$) CO($J = 2 \rightarrow 1$) jet powered by the intermediate-mass protostellar source Serpens SMM1-a. The highly collimated molecular jet is flanked at the base by a wide-angle cavity; the walls of the cavity can be seen in both 4 cm free-free emission detected by the VLA and 1.3 mm thermal dust emission detected by ALMA. This is the first time that ionization of an outflow cavity has been directly detected via free-free emission in a very young, embedded Class 0 protostellar source that is still powering a molecular jet. The cavity walls are ionized either by UV photons escaping from the accreting protostellar source, or by the precessing molecular jet impacting the walls. These observations suggest that ionized outflow cavities may be common in Class 0 protostellar sources, shedding further light on the radiation, outflow, and jet environments in the youngest, most embedded forming stars.

With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A simple model of a molecular jet precessing around an axis close to the plane of the sky can reproduce the observations. The velocity of the outflowing gas is difficult to constrain due to the velocity reversals but we estimate it to be between 240 and 850 km/s and the jet to precess with a period P=0.3-1.1 Myr. The CO emission is clumpy along the jet and the total molecular mass in the high-velocity (+-(60 to 150 km/s)) gas lies between 2e6 Msun (light jet) and 2e7 Msun (massive jet). There is also CO emission extending along the minor axis of NGC1377. It holds >40% of the flux in NGC1377 and may be a slower, wide-angle molecular outflow which is partially entrained by the molecular jet. We discuss the driving mechanism of the molecular jet and suggest that it is either powered by a very faint radio jet or by an accretion disk-wind similar to those found towards protostars. The nucleus of NGC1377 harbours intense embedded activity and we detect emission from vibrationally excited HCN J=4-3 v_2=1f which is consistent with hot gas and dust. We find large columns of H2 in the centre of NGC1377 which may be a sign of a high rate of recent gas infall. The dynamical age of the molecular jet is short (<1 Myr), which could imply that it is young and consistent with the notion that NGC1377 is caught in a transient phase of its evolution. However, further studies are required to determine the age of the molecular jet, its mass and the role it is playing in the growth of the nucleus of NGC1377.

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $<$60 au like the millimeter-sized grains) and two vertical distributions (near the midplane where water is expected to photodesorb off icy grains) to describe the (unknown) location of the molecules. These distributions capture the effects of radial drift and vertical settling of ice-covered grains. We use physical-chemical models to reproduce the fluxes with assuming that water and ammonia are co-spatial. We infer ammonia gas-phase masses of 0.7-11.0 $\times$10$^{21}$ g. For water, we infer gas-phase masses of 0.2-16.0 $\times$10$^{22}$ g. This corresponds to NH$_3$/H$_2$O abundance ratios of 7\%-84\%, assuming that water and ammonia are co-located. Only in the most compact and settled adopted configuration is the inferred NH$_3$/H$_2$O consistent with interstellar ices and solar system bodies of $\sim$ 5\%-10\%. Volatile release in the midplane may occur via collisions between icy bodies if the available surface for subsequent freeze-out is significantly reduced, e.g., through growth of small grains into pebbles or larger.

Using the Submillimeter Array (SMA) on Mauna Kea, the H2-16O 10_2,9-9_3,6 transition (E_up=1863K) at 321.2 GHz has been detected toward the embedded low-mass protostar HL Tau. The line centroid is blue-shifted by 15 km/s with respect to the source velocity, and it has a FWHM of 20 km/s. The emission is tentatively resolved and extends ~3-4" over the sky (~2 beams), or ~500 AU at the distance of Taurus. The velocity offset, and to a lesser degree the spatial extent of the emission, shows that the line originates in the protostellar jet or wind. This result suggests that at least some water emission observed toward embedded sources, and perhaps also disk sources, with Herschel and Spitzer contains a wind or jet component, which is crucial for interpreting these data. These pathfinder observations done with the SMA opens a new window to studying the origin of water emission with e.g. ALMA, thus providing new insights into where water is in protostellar systems.

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.

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

We present Herschel/HIFI observations (WISH KP) of 14 water lines in a small sample of galactic massive protostellar objects: NGC6334I(N), DR21(OH), IRAS16272-4837, and IRAS05358+3543. We analyze the gas dynamics from the line profiles. Through modeling of the observations using RATRAN, we estimate outflow, infall, turbulent velocities, molecular abundances, and investigate any correlation with the evolutionary status of each source. The molecular line profiles exhibit a broad component coming from the shocks along the cavity walls associated with the protostars, and an infalling (or expansion for IRAS05358+3543) and passively heated envelope component, with highly supersonic turbulence likely increasing with the distance from the center. Accretion rates between 6.3 10^{-5} and 5.6 10^{-4} \msun yr^{-1} are derived from the infall observed in three of our sources. The outer water abundance is estimated to be at the typical value of a few 10^{-8} while the inner abundance varies from 1.7 10^{-6} to 1.4 10^{-4} with respect to H2 depending on the source. We confirm that regions of massive star formation are highly turbulent and that the turbulence likely increases in the envelope with the distance to the star. The inner abundances are lower than the expected 10^{-4} perhaps because our observed lines do not probe deep enough into the inner envelope, or because photodissociation through protostellar UV photons is more efficient than expected. We show that the higher the infall/expansion velocity in the protostellar envelope, the higher is the inner abundance, maybe indicating that larger infall/expansion velocities generate shocks that will sputter water from the ice mantles of dust grains in the inner region. High-velocity water must be formed in the gas-phase from shocked material.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

We present millimetre (SMA) and sub-millimetre (SCUBA-2) continuum observations of the peculiar star KIC 8462852 which displayed several deep and aperiodic dips in brightness during the Kepler mission. Our observations are approximately confusion-limited at 850 $\mu$m and are the deepest millimetre and sub-millimetre photometry of the star that has yet been carried out. No significant emission is detected towards KIC 8462852. We determine upper limits for dust between a few 10$^{-6}$ M$_{\oplus}$ and 10$^{-3}$ M$_{\oplus}$ for regions identified as the most likely to host occluding dust clumps and a total overall dust budget of $<$7.7 M$_{\oplus}$ within a radius of 200 AU. Such low limits for the inner system make the catastrophic planetary disruption hypothesis unlikely. Integrating over the Kepler lightcurve we determine that at least 10$^{-9}$ M$_{\oplus}$ of dust is required to cause the observed Q16 dip. This is consistent with the currently most favoured cometary breakup hypothesis, but nevertheless implies the complete breakup of $\sim$ 30 Comet 1/P Halley type objects. Finally, in the wide SCUBA-2 field-of-view we identify another candidate debris disc system that is potentially the largest yet discovered.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

Water probes the dynamics in young stellar objects (YSOs) effectively, especially shocks in molecular outflows. It is a key molecule for exploring whether the physical properties of low-mass protostars can be extrapolated to massive YSOs. As part of the WISH key programme, we investigate the dynamics and the excitation conditions of shocks along the outflow cavity wall as function of source luminosity. Velocity-resolved Herschel-HIFI spectra of the H2O 988, 752, 1097 GHz and 12CO J=10-9, 16-15 lines were analysed for 52 YSOs with bolometric luminosities (L_bol) ranging from <1 to >10^5 L_sun. The profiles of the H2O lines are similar, indicating that they probe the same gas. We see two main Gaussian emission components in all YSOs: a broad component associated with non-dissociative shocks in the outflow cavity wall (cavity shocks) and a narrow component associated with quiescent envelope material. More than 60% of the total integrated intensity of the H2O lines (L_H2O) comes from the cavity shock component. The H2O line widths are similar for all YSOs, whereas those of 12CO 10-9 increase slightly with L_bol. The excitation analysis of the cavity shock component, performed with the non-LTE radiative transfer code RADEX, shows stronger 752 GHz emission for high-mass YSOs, likely due to pumping by an infrared radiation field. As previously found for CO, a strong correlation with slope unity is measured between log(L_H2O) and log(L_bol), which can be extrapolated to extragalactic sources. We conclude that the broad component of H2O and high-J CO lines originate in shocks in the outflow cavity walls for all YSOs, whereas lower-J CO transitions mostly trace entrained outflow gas. The higher UV field and turbulent motions in high-mass objects compared to their low-mass counterparts may explain the slightly different kinematical properties of 12CO 10-9 and H2O lines from low- to high-mass YSOs.

OH is a key molecule in H2O chemistry, a valuable tool for probing physical conditions, and an important contributor to the cooling of shock regions. OH participates in the re-distribution of energy from the protostar towards the surrounding ISM. Our aim is to assess the origin of the OH emission from the Cepheus A massive star-forming region and to constrain the physical conditions prevailing in the emitting gas. We thus want to probe the processes at work during the formation of massive stars. We present spectrally resolved observations of OH towards the outflows of Cepheus A with the GREAT spectrometer onboard the SOFIA telescope. Three triplets were observed at 1834.7 GHz, 1837.8 GHz, and 2514.3 GHz (163.4, 163.1, and 119.2 microns), at angular resolutions of 16.3", 16.3", and 11.9", respectively. We present the CO (16-15) spectrum at the same position. We compared the integrated intensities in the redshifted wings to shock models. The two triplets near 163 microns are detected in emission with blending hyperfine structure unresolved. Their profiles and that of CO can be fitted by a combination of 2 or 3 Gaussians. The observed 119.2 microns triplet is seen in absorption, since its blending hyperfine structure is unresolved, but with three line-of-sight components and a blueshifted emission wing consistent with that of the other lines. The OH line wings are similar to those of CO, suggesting that they emanate from the same shocked structure. Under this common origin assumption, the observations fall within the model predictions and within the range of use of our model only if we consider that four shock structures are caught in our beam. Our comparisons suggest that the observations might be consistently fitted by a J-type model with nH > 1e5 cm-3, v > 20 km/s, and with a filling factor of ~1. Such a high density is generally found in shocks associated to high-mass protostars.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

Most low-mass protostars form in clusters, in particular high-mass clusters; however, how low-mass stars form in high-mass clusters and what the mass distribution is, are still open questions both in our own Galaxy and elsewhere. To access the population of forming embedded low-mass protostars observationally, we propose to use molecular outflows as tracers. Because the outflow emission scales with mass, the effective contrast between low-mass protostars and their high-mass cousins is greatly lowered. In particular, maps of methanol emission at 338.4 GHz (J=7_0 - 6_0 A+) in low-mass clusters illustrate that this transition is an excellent probe of the low-mass population. We here present a model of a forming cluster where methanol emission is assigned to every embedded low-mass protostar. The resulting model image of methanol emission is compared to recent ALMA observations toward a high-mass cluster and the similarity is striking: the toy model reproduces observations to better than a factor of two and suggests that approximately 50\% of the total flux originates in low-mass outflows. Future fine-tuning of the model will eventually make it a tool for interpreting the embedded low-mass population of distant regions within our own Galaxy and ultimately higher-redshift starburst galaxies, not just for methanol emission but also water and high-J CO.

Intermediate-mass young stellar objects (YSOs) provide a link to understand how feedback from shocks and UV radiation scales from low to high-mass star forming regions. Aims: Our aim is to analyze excitation of CO and H$_2$O in deeply-embedded intermediate-mass YSOs and compare with low-mass and high-mass YSOs. Methods: Herschel/PACS spectral maps are analyzed for 6 YSOs with bolometric luminosities of $L_\mathrm{bol}\sim10^2 - 10^3$ $L_\odot$. The maps cover spatial scales of $\sim 10^4$ AU in several CO and H$_2$O lines located in the $\sim55-210$ $\mu$m range. Results: Rotational diagrams of CO show two temperature components at $T_\mathrm{rot}\sim320$ K and $T_\mathrm{rot}\sim700-800$ K, comparable to low- and high-mass protostars probed at similar spatial scales. The diagrams for H$_2$O show a single component at $T_\mathrm{rot}\sim130$ K, as seen in low-mass protostars, and about $100$ K lower than in high-mass protostars. Since the uncertainties in $T_\mathrm{rot}$ are of the same order as the difference between the intermediate and high-mass protostars, we cannot conclude whether the change in rotational temperature occurs at a specific luminosity, or whether the change is more gradual from low- to high-mass YSOs. Conclusions: Molecular excitation in intermediate-mass protostars is comparable to the central $10^{3}$ AU of low-mass protostars and consistent within the uncertainties with the high-mass protostars probed at $3\cdot10^{3}$ AU scales, suggesting similar shock conditions in all those sources.

During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow regulates star formation by e.g. inducing turbulence and heating the surrounding gas. Understanding the associated shocks is therefore essential to the study of star formation. We present comparisons of shock models with CO, H2, and SiO observations in a 'pure' shock position in the BHR71 bipolar outflow. These comparisons provide an insight into the shock and pre-shock characteristics, and allow us to understand the energetic and chemical feedback of star formation on Galactic scales. New CO (Jup = 16, 11, 7, 6, 4, 3) observations from the shocked regions with the SOFIA and APEX telescopes are presented and combined with earlier H2 and SiO data (from the Spitzer and APEX telescopes). The integrated intensities are compared to a grid of models that were obtained from a magneto-hydrodynamical shock code which calculates the dynamical and chemical structure of these regions combined with a radiative transfer module based on the 'large velocity gradient' approximation. The CO emission leads us to update the conclusions of our previous shock analysis: pre-shock densities of 1e4 cm-3 and shock velocities around 20-25 km s-1 are still constrained, but older ages are inferred ( 4000 years). We evaluate the contribution of shocks to the excitation of CO around forming stars. The SiO observations are compatible with a scenario where less than 4% of the pre-shock SiO belongs to the grain mantles. We infer outflow parameters: a mass of 1.8x1e-2 Msun was measured in our beam, in which a momentum of 0.4 Msun km s-1 is dissipated, for an energy of 4.2x1e43erg. We analyse the energetics of the outflow species by species. Comparing our results with previous studies highlights their dependence on the method: H2 observations only are not sufficient to evaluate the mass of outflows.

We present Herschel PACS mapping observations of the [OI]63 micron line towards protostellar outflows in the L1448, NGC1333-IRAS4, HH46, BHR71 and VLA1623 star forming regions. We detect emission spatially resolved along the outflow direction, which can be associated with a low excitation atomic jet. In the L1448-C, HH46 IRS and BHR71 IRS1 outflows this emission is kinematically resolved into blue- and red-shifted jet lobes, having radial velocities up to 200 km/s. In the L1448-C atomic jet the velocity increases with the distance from the protostar, similarly to what observed in the SiO jet associated with this source. This suggests that [OI] and molecular gas are kinematically connected and that this latter could represent the colder cocoon of a jet at higher excitation. Mass flux rates (\.M$_{jet}$(OI)) have been measured from the [OI]63micron luminosity adopting two independent methods. We find values in the range 1-4 10$^{-7}$ Mo/yr for all sources but HH46, for which an order of magnitude higher value is estimated. \.M$_{jet}$(OI) are compared with mass accretion rates (\.M$_{acc}$) onto the protostar and with \.M$_{jet}$ derived from ground-based CO observations. \.M$_{jet}$(OI)/\.M$_{acc}$ ratios are in the range 0.05-0.5, similar to the values for more evolved sources. \.M$_{jet}$(OI) in HH46 IRS and IRAS4A are comparable to \.M$_{jet}$(CO), while those of the remaining sources are significantly lower than the corresponding \.M$_{jet}$(CO). We speculate that for these three sources most of the mass flux is carried out by a molecular jet, while the warm atomic gas does not significantly contribute to the dynamics of the system.

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

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.

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.

Aims: Our aim is to determine the critical parameters in water chemistry and the contribution of water to the oxygen budget by observing and modelling water gas and ice for a sample of eleven low-mass protostars, for which both forms of water have been observed. Methods: A simplified chemistry network, which is benchmarked against more sophisticated chemical networks, is developed that includes the necessary ingredients to determine the water vapour and ice abundance profiles in the cold, outer envelope in which the temperature increases towards the protostar. Comparing the results from this chemical network to observations of water emission lines and previously published water ice column densities, allows us to probe the influence of various agents (e.g., FUV field, initial abundances, timescales, and kinematics). Results: The observed water ice abundances with respect to hydrogen nuclei in our sample are 30-80ppm, and therefore contain only 10-30% of the volatile oxygen budget of 320ppm. The keys to reproduce this result are a low initial water ice abundance after the pre-collapse phase together with the fact that atomic oxygen cannot freeze-out and form water ice in regions with T(dust)>15 K. This requires short prestellar core lifetimes of less than about 0.1Myr. The water vapour profile is shaped through the interplay of FUV photodesorption, photodissociation, and freeze-out. The water vapour line profiles are an invaluable tracer for the FUV photon flux and envelope kinematics. Conclusions: The finding that only a fraction of the oxygen budget is locked in water ice can be explained either by a short pre-collapse time of less than 0.1 Myr at densities of n(H)~1e4 cm-3, or by some other process that resets the initial water ice abundance for the post-collapse phase. A key for the understanding of the water ice abundance is the binding energy of atomic oxygen on ice.

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.

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.

Water in outflows from protostars originates either as a result of gas-phase synthesis from atomic oxygen at T > 200 K, or from sputtered ice mantles containing water ice. We aim to quantify the contribution of the two mechanisms that lead to water in outflows, by comparing observations of gas-phase water to methanol (a grain surface product) towards three low-mass protostars in NGC1333. In doing so, we also quantify the amount of methanol destroyed in outflows. To do this, we make use of JCMT and Herschel-HIFI data of H2O, CH3OH and CO emission lines and compare them to RADEX non-LTE excitation simulations. We find up to one order of magnitude decrease in the column density ratio of CH3OH over H2O as the velocity increases in the line wings up to ~15 km/s. An independent decrease in X(CH3OH) with respect to CO of up to one order of magnitude is also found in these objects. We conclude that gas-phase formation of H2O must be active at high velocities (above 10 km/s, relative to the source velocity) to re-form the water destroyed during sputtering. In addition, the transition from sputtered water at low velocities to formed water at high velocities must be gradual. We place an upper limit of two orders of magnitude on the destruction of methanol by sputtering effects.

We wish to test the origins proposed for the extended hot H2 at 2000K around the atomic jet from the T Tauri star DGTau, in order to constrain the wide-angle wind structure and the possible presence of an MHD disk wind. We present flux calibrated IFS observations in H2 1-0 S(1) obtained with SINFONI/VLT. Thanks to spatial deconvolution by the PSF and to accurate correction for uneven slit illumination, we performed a thorough analysis and modeled the morphology, kinematics, and surface brightness. We also compared our results with studies in [FeII], [OI], and FUV-pumped H2. The limb-brightened H2 emission in the blue lobe is strikingly similar to FUV-pumped H2 imaged 6yr later, confirming that they trace the same hot gas and setting an upper limit of 12km/s on any expansion proper motion. The wide-angle H2 rims are at lower blueshifts than probed by narrow long-slit spectra. We confirm that they extend to larger angle and to lower speed the onion-like velocity structure observed in optical atomic lines. The latter is shown to be steady over more/equal than 4yr but undetected in [FeII] by SINFONI, probably due to strong iron depletion. The H2 rim thickness less/equal than 14AU rules out excitation by C-shocks, and J-shock speeds are constrained to 10km/s. We find that explaining the H2 wide-angle emission with a shocked layer requires either a recent outburst (15yr) into a pre-existing ambient outflow or an excessive wind mass flux. A slow photoevaporative wind from the dense irradiated disk surface and an MHD disk wind heated by ambipolar diffusion seem to be more promising and need to be modeled in more detail.

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

Aims. The aim of this paper is to study deuterated water in the solar-type protostars NGC1333 IRAS4A and IRAS4B, to compare their HDO abundance distribution with other star-forming regions, and to constrain their HDO/H2O ratios. Methods. Using the Herschel/HIFI instrument as well as ground-based telescopes, we observed several HDO lines covering a large excitation range (Eup/k=22-168 K) towards these protostars and an outflow position. Non-LTE radiative transfer codes were then used to determine the HDO abundance profiles in these sources. Results. The HDO fundamental line profiles show a very broad component, tracing the molecular outflows, in addition to a narrower emission component and a narrow absorbing component. In the protostellar envelope of NGC1333 IRAS4A, the HDO inner (T>100 K) and outer (T<100 K) abundances with respect to H2 are estimated at 7.5x10^{-9} and 1.2x10^{-11}, respectively, whereas, in NGC1333 IRAS4B, they are 1.0x10^{-8} and 1.2x10^{-10}, respectively. Similarly to the low-mass protostar IRAS16293-2422, an absorbing outer layer with an enhanced abundance of deuterated water is required to reproduce the absorbing components seen in the fundamental lines at 465 and 894 GHz in both sources. This water-rich layer is probably extended enough to encompass the two sources as well as parts of the outflows. In the outflows emanating from NGC1333 IRAS4A, the HDO column density is estimated at about (2-4)x10^{13} cm^{-2}, leading to an abundance of about (0.7-1.9)x10^{-9}. An HDO/H2O ratio between 7x10^{-4} and 9x10^{-2} is derived in the outflows. In the warm inner regions of these two sources, we estimate the HDO/H2O ratios at about 1x10^{-4}-4x10^{-3}. This ratio seems higher (a few %) in the cold envelope of IRAS4A, whose possible origin is discussed in relation to formation processes of HDO and H2O.

Abridged abstract: For stars to form, material must fall inwards from core scales through the envelope towards the central protostar. The velocity profile around protostars is poorly constrained. 6 Class 0 protostars and one Class I protostars observed with HIFI on board Herschel as part of the "Water in star-forming regions with Herschel" (WISH) survey show infall signatures in water line observations. We use 1-D non-LTE RATRAN radiative transfer models of the observed water lines to constrain the infall velocity and chemistry in the protostellar envelopes of these sources. We assume a free-fall velocity profile and, having found the best fit, vary the radii over which infall takes place. In the well-studied Class 0 protostar NGC1333-IRAS4A we find that infall takes place over the whole envelope to which our observations are sensitive (r>~1000 AU). For 4 sources infall takes place on core to envelope scales (i.e. ~10000-3000 AU). In 2 sources the inverse P-Cygni profiles seen in the ground-state lines are more likely due to larger-scale motions or foreground clouds. Models including a simple consideration of the chemistry are consistent with the observations, while using step abundance profiles are not. The non-detection of excited water in the inner envelope in 6/7 protostars is further evidence that water must be heavily depleted from the gas-phase at these radii. Infall in four of the sources is supersonic and infall in all sources must take place at the outer edge of the envelope, which may be evidence that collapse is global or outside-in rather than inside-out. The mass infall rate in IRAS4A is large (>~10^-4 msol\yr), higher than the mass outflow rate and expected mass accretion rates onto the star, suggesting that any flattened disk-like structure on small scales will be gravitationally unstable, potentially leading to rotational fragmentation and/or episodic accretion.

We present Herschel/PACS observations of L1448-MM, a Class 0 protostar with a prominent outflow. Numerous emission lines are detected at 55 < lambda < 210 micrometer including CO, OH, H2O, and [O I]. We investigate the spatial distribution of each transition to find that lines from low energy levels tend to distribute along the outflow direction while lines from high energy levels peak at the central spatial pixel. Spatial maps reveal that OH emission lines are formed in a relatively small area, while [O I] emission is extended. According to the rotational diagram analysis, the CO emission can be fitted by two (warm and hot) temperature components. For H2O, the ortho-to-para ratio is close to 3. The non-LTE LVG calculations suggest that CO and H2O lines could instead be formed in a high kinetic temperature (T > 1000 K) environment, indicative of a shock origin. For OH, IR-pumping processes play an important role in the level population. The molecular emission in L1448-MM is better explained with a C-shock model, but the atomic emission of PACS [O I] and Spitzer/IRS [Si II] emission is not consistent with C-shocks, suggesting multiple shocks in this region. Water is the major line coolant of L1448-MM in the PACS wavelength range, and the best-fit LVG models predict that H2O and CO emit (50-80)% of their line luminosity in the PACS wavelength range.

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

We report on spectro-imaging observations employing Spitzer IRS and Herschel PACS, aiming to constrain the physical conditions around SMM3 and SMM4 in Serpens. The combined power of both instruments provides an almost complete wavelength coverage between 5 and 200 micron at an angular resolution of 10". We detect line emission from all major molecular (H2, CO, H2O and OH) and many atomic ([OI], [CII], [FeII], [SiII] and [SI]) coolants. Line emission tends to peak at distances of 10 - 20" from the protostellar sources, at positions of known outflow shocks. The only exception is [CII] which likely traces a PDR excited from the neighboring source SMM6. Excitation analysis indicates that H2 and CO originate from gas at two distinct rotational temperatures of 300 K and 1000 K, while H2O and OH emission corresponds to rotational temperatures of 100 - 200 K. The morphological and physical association between CO and H2 suggests a common excitation mechanism which allows direct comparisons between the two molecules. The CO/H2 abundance ratio varies from 10^-5 in the warm gas up to 10^-4 in the hotter regions. The occurrence of J-shocks is suggested by the strong atomic/ionic (except for [CII]) emission as well as a number of line ratio diagnostics. Both C- and J-shocks can account for the observed molecular emission, however J-shocks are strongly advocated by the atomic emission and provide simpler and more homogeneous solutions for CO and H2. C-shocks describe better the emission from H2O and OH. The variations in the CO/H2 abundance ratio for gas at different temperatures can be interpreted by their reformation rates in dissociative J-type shocks, or the simultaneous influence of both C and J shocks.

Herschel-HIFI spectra of H2O towards low-mass protostars show a distinct velocity component not seen in observations from the ground of CO or other species. The aim is to characterise this component in terms of excitation conditions and physical origin. A velocity component with an offset of ~10 km/s detected in spectra of the H2O 110-101 557 GHz transition towards six low-mass protostars in the 'Water in star-forming regions with Herschel' (WISH) programme is also seen in higher-excited H2O lines. The emission from this component is quantified and excitation conditions are inferred using 1D slab models. Data are compared to observations of hydrides (high-J CO, OH+, CH+, C+, OH) where the same component is uniquely detected. The velocity component is detected in all 6 targeted H2O transitions (Eup~50-250K), and in CO 16-15 towards one source, Ser SMM1. Inferred excitation conditions imply that the emission arises in dense (n~5x10^6-10^8 cm^-3) and hot (T~750K) gas. The H2O and CO column densities are ~10^16 and 10^18 cm^-2, respectively, implying a low H2O abundance of 10^-2 with respect to CO. The high column densities of ions such as OH+ and CH+ (both ~10^13 cm^-2) indicate an origin close to the protostar where the UV field is strong enough that these species are abundant. The estimated radius of the emitting region is 100AU. This component likely arises in dissociative shocks close to the protostar, an interpretation corroborated by a comparison with models of such shocks. Furthermore, one of the sources, IRAS4A, shows temporal variability in the offset component over a period of two years which is expected from shocks in dense media. High-J CO gas detected with Herschel-PACS with Trot~700K is identified as arising in the same component and traces the part of the shock where H2 reforms. Thus, H2O reveals new dynamical components, even on small spatial scales in low-mass protostars.

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.

The outflow force of molecular bipolar outflows is a key parameter in theories of young stellar feedback on their surroundings. The focus of many outflow studies is the correlation between the outflow force, bolometric luminosity and envelope mass. However, it is difficult to combine the results of different studies in large evolutionary plots over many orders of magnitude due to the range of data quality, analysis methods and corrections for observational effects such as opacity and inclination. We aim to determine the outflow force for a sample of low luminosity embedded sources. We will quantify the influence of the analysis method and the assumptions entering the calculation of the outflow force. We use the James Clerk Maxwell Telescope to map 12CO J=3-2 over 2'x2' regions around 16 Class I sources of a well-defined sample in Ophiuchus at 15" resolution. The outflow force is then calculated using seven different methods differing e.g. in the use of intensity-weighted emission and correction factors for inclination. The results from the analysis methods differ from each other by up to a factor of 6, whereas observational properties and choices in the analysis procedure affect the outflow force by up to a factor of 4. For the sample of Class I objects, bipolar outflows are detected around 13 sources including 5 new detections, where the three non-detections are confused by nearby outflows from other sources. When combining outflow forces from different studies, a scatter by up to a factor of 5 can be expected. Although the true outflow force remains unknown, the separation method (separate calculation of dynamical time and momentum) is least affected by the uncertain observational parameters. The correlations between outflow force, bolometric luminosity and envelope mass are further confirmed down to low luminosity sources.

We present 50-210 um spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel-PACS, and 0.5-1000 um SEDs, as part of the Dust, Ice, and Gas in Time (DIGIT) Key Program. Some sources exhibit up to 75 H2O lines ranging in excitation energy from 100-2000 K, 12 transitions of OH, and CO rotational lines ranging from J=14-13 up to J=40-39. [O I] is detected in all but one source in the entire sample; among the sources with detectable [O I] are two Very Low Luminosity Objects (VeLLOs). The mean 63/145 um [O I] flux ratio is 17.2 +/- 9.2. The [O I] 63 um line correlates with Lbol, but not with the time-averaged outflow rate derived from low-J CO maps. [C II] emission is in general not local to the source. The sample Lbol increased by 1.25 (1.06) and Tbol decreased to 0.96 (0.96) of mean (median) values with the inclusion of the Herschel data. Most CO rotational diagrams are characterized by two optically thin components (<N> = (0.70 +/- 1.12) x 10^49 total particles). N_CO correlates strongly with Lbol, but neither Trot nor N_CO(warm)/N_CO(hot) correlates with Lbol, suggesting that the total excited gas is related to the current source luminosity, but that the excitation is primarily determined by the physics of the interaction (e.g., UV- heating/shocks). Rotational temperatures for H2O (<Trot> = 194 +/- 85 K) and OH (<Trot> = 183 +/- 117 K) are generally lower than for CO, and much of the scatter in the observations about the best fit is attributed to differences in excitation conditions and optical depths amongst the detected lines.

(Abridged) Star and planet formation theories predict an evolution in the density, temperature, and velocity structure as the envelope collapses and forms an accretion disk. The aim of this work is to model the evolution of the molecular excitation, line profiles, and related observables during low-mass star formation. Specifically, the signatures of disks during the deeply embedded stage are investigated. Semi-analytic 2D axisymmetric models have been used to describe the evolution of the density, stellar mass, and luminosity from the pre-stellar to the T-Tauri phase. A full radiative transfer calculation is carried out to accurately determine the time-dependent dust temperatures and CO abundance structure. We present non-LTE near-IR, FIR, and submm lines of CO have been simulated at a number of time steps. In contrast to the dust temperature, the CO excitation temperature derived from submm/FIR lines does not vary during the protostellar evolution, consistent with C18O observations obtained with Herschel and from ground-based telescopes. The near-IR spectra provide complementary information to the submm lines by probing not only the cold outer envelope but also the warm inner region. The near-IR high-J (>8) absorption lines are particularly sensitive to the physical structure of the inner few AU, which does show evolution. High signal-to-noise ratio subarcsec resolution data with ALMA are needed to detect the presence of small rotationally supported disks during the Stage 0 phase and various diagnostics are discussed.

Evaporation of water ice above 100 K in the inner few 100 AU of low-mass embedded protostars (the so-called hot core) should produce quiescent water vapor abundances of ~10^-4 relative to H2. Observational evidence so far points at abundances of only a few 10^-6. However, these values are based on spherical models, which are known from interferometric studies to be inaccurate on the relevant spatial scales. Are hot cores really that much drier than expected, or are the low abundances an artifact of the inaccurate physical models? We present deep velocity-resolved Herschel-HIFI spectra of the 3(12)-3(03) lines of H2-16O and H2-18O (1097 GHz, Eup/k = 249 K) in the low-mass Class 0 protostar NGC1333 IRAS2A. A spherical radiative transfer model with a power-law density profile is unable to reproduce both the HIFI data and existing interferometric data on the H2-18O 3(13)-2(20) line (203 GHz, Eup/k = 204 K). Instead, the HIFI spectra likely show optically thick emission from a hot core with a radius of about 100 AU. The mass of the hot core is estimated from the C18O J=9-8 and 10-9 lines. We derive a lower limit to the hot water abundance of 2x10^-5, consistent with the theoretical predictions of ~10^-4. The revised HDO/H2O abundance ratio is 1x10^-3, an order of magnitude lower than previously estimated.

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

Context; Our understanding of the star formation process has traditionally been confined to certain mass or luminosity boundaries because most studies focus only on low-, intermediate- or high-mass star-forming regions. As part of the "Water In Star-forming regions with Herschel" (WISH) key program, water and other important molecules, such as CO and OH, have been observed in 51 embedded young stellar objects (YSOs). The studied sample covers a range of luminosities from <1 to >10^5 L_sol. Aims; We analyse the CO line emission towards a large sample of protostars in terms of both line intensities and profiles. Methods; Herschel-HIFI spectra of the 12CO 10-9, 13CO 10-9 and C18O 5-4, 9-8 and 10-9 lines are analysed for a sample of 51 YSOs. In addition, JCMT spectra of 12CO 3-2 and C18O 3-2 extend this analysis to cooler gas components. Results; All observed CO and isotopologue spectra show a strong linear correlation between the logarithms of the line and bolometric luminosities across six orders of magnitude on both axes. This suggests that the high-J CO lines primarily trace the amount of dense gas associated with YSOs. This relation can be extended to larger (extragalactic) scales. The majority of the detected 12CO line profiles can be decomposed into a broad and a narrow Gaussian component, while the C18O spectra are mainly fitted with a single Gaussian. A broadening of the line profile is also observed from pre-stellar cores to embedded protostars, which is due mostly to non-thermal motions (turbulence/infall). The widths of the broad 12CO 3-2 and 10-9 velocity components correlate with those of the narrow C18O 9-8 profiles, suggesting that the entrained outflowing gas and envelope motions are related independent of the mass of the protostar. These results indicate that physical processes in protostellar envelopes have similar characteristics across the studied luminosity range.

OH is a key species in the water chemistry of star-forming regions, because its presence is tightly related to the formation and destruction of water. This paper presents OH observations from 23 low- and intermediate-mass young stellar objects obtained with the PACS integral field spectrometer on-board Herschel in the context of the Water In Star-forming Regions with Herschel (WISH) key program. Most low-mass sources have compact OH emission (< 5000 AU scale), whereas the OH lines in most intermediate-mass sources are extended over the whole PACS detector field-of-view (> 20000 AU). The strength of the OH emission is correlated with various source properties such as the bolometric luminosity and the envelope mass, but also with the OI and H2O emission. Rotational diagrams for sources with many OH lines show that the level populations of OH can be approximated by a Boltzmann distribution with an excitation temperature at around 70 K. Radiative transfer models of spherically symmetric envelopes cannot reproduce the OH emission fluxes nor their broad line widths, strongly suggesting an outflow origin. Slab excitation models indicate that the observed excitation temperature can either be reached if the OH molecules are exposed to a strong far-infrared continuum radiation field or if the gas temperature and density are sufficiently high. Using realistic source parameters and radiation fields, it is shown for the case of Ser SMM1 that radiative pumping plays an important role in transitions arising from upper level energies higher than 300 K. The compact emission in the low-mass sources and the required presence of a strong radiation field and/or a high density to excite the OH molecules points towards an origin in shocks in the inner envelope close to the protostar.

Observations of higher-excited transitions of abundant molecules such as CO are important for determining where energy in the form of shocks is fed back into the parental envelope of forming stars. The nearby prototypical and protobinary low-mass hot core, IRAS16293-2422 (I16293) is ideal for such a study. The source was targeted with ALMA for science verification purposes in band 9, which includes CO J=6-5 (E_up/k_B ~ 116 K), at an unprecedented spatial resolution (~0.2", 25 AU). I16293 itself is composed of two sources, A and B, with a projected distance of 5". CO J=6-5 emission is detected throughout the region, particularly in small, arcsecond-sized hotspots, where the outflow interacts with the envelope. The observations only recover a fraction of the emission in the line wings when compared to data from single-dish telescopes, with a higher fraction of emission recovered at higher velocities. The very high angular resolution of these new data reveal that a bow shock from source A coincides, in the plane of the sky, with the position of source B. Source B, on the other hand, does not show current outflow activity. In this region, outflow entrainment takes place over large spatial scales, >~ 100 AU, and in small discrete knots. This unique dataset shows that the combination of a high-temperature tracer (e.g., CO J=6-5) and very high angular resolution observations is crucial for interpreting the structure of the warm inner environment of low-mass protostars.

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

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.

We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE radiative transfer model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s < 20 km s^-1) are needed to explain both the atomic lines and the hot CO and H2O lines respectively. Observations also show the signature of UV radiation and thus, most observed species likely arise in UV-irradiated shocks. Dissociative J-shocks produced by an atomic jet are the most probable origin of [OI] and OH emission and of a significant fraction of the warm CO emission. In addition, H2O photodissociation in UV-irradiated non-dissociative shocks can also contribute to the [OI] and OH emission.

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.

We present Herschel/HIFI observations of fourteen water lines in W43-MM1, a massive protostellar object in the luminous star cluster-forming region W43. We analyze the gas dynamics from the line profiles using Herschel-HIFI observations (WISH-KP) of fourteen far-IR water lines (H2O, H217O, H218O), CS(11-10), and C18O(9-8) lines, and using our modeling of the continuum spectral energy distribution. As for lower mass protostellar objects, the molecular line profiles are a mix of emission and absorption, and can be decomposed into 'medium', and 'broad' velocity components. The broad component is the outflow associated with protostars of all masses. Our modeling shows that the remainder of the water profiles can be well fitted by an infalling and passively heated envelope, with highly supersonic turbulence varying from 2.2 km/s in the inner region to 3.5 km/s in the outer envelope. Also, W43-MM1 has a high accretion rate, between 4.0 x 10^{-4} and 4.0 x 10^{-2} \msun /yr, derived from the fast (0.4-2.9 km/s) infall observed. We estimate a lower mass limit of gaseous water of 0.11 \msun and total water luminosity of 1.5 \lsun (in the 14 lines presented here). The central hot core is detected with a water abundance of 1.4 x 10^{-4} while the water abundance for the outer envelope is 8 x10^{-8}. The latter value is higher than in other sources, most likely related to the high turbulence and the micro-shocks created by its dissipation. Examining water lines of various energies, we find that the turbulent velocity increases with the distance to the center. While not in clear disagreement with the competitive accretion scenario, this behavior is predicted by the turbulent core model. Moreover, the estimated accretion rate is high enough to overcome the expected radiation pressure.

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

NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage 0 low-mass protostars which are driving prominent bipolar outflows. Most studies have so far concentrated on the colder parts (T<30K) of these regions. The aim is to characterize the warmer parts of the protostellar envelope in order to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation and outflow dispersal. Fully sampled large scale maps of the region were obtained; APEX-CHAMP+ was used for 12CO 6-5, 13CO 6-5 and [CI] 2-1, and JCMT-HARP-B for 12CO 3-2 emissions. Complementary Herschel-HIFI and ground-based lines of CO and its isotopologs, from 1-0 upto 10-9 (Eu/k 300K), are collected at the source positions. Radiative-transfer models of the dust and lines are used to determine temperatures and masses of the outflowing and UV-heated gas and infer the CO abundance structure. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~100K. At other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13CO 6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-heated and outflowing gas are found to be comparable from the 12CO and 13CO 6-5 maps, implying that UV photons can affect the gas as much as the outflows. Weak [CI] emission throughout the region indicates a lack of CO dissociating photons. Modeling of the C18O lines indicates the necessity of a "drop" abundance profile throughout the envelope where the CO freezes out and is reloaded back into the gas phase, thus providing quantitative evidence for the CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H_2=2.7x10^-4, indicating some processing of CO into other species on the grains.

NGC 7129 FIRS 2 is a young intermediate-mass (IM) protostar, which is associated with two energetic bipolar outflows and displays clear signs of the presence of a hot core. It has been extensively observed with ground based telescopes and within the WISH Guaranteed Time Herschel Key Program. We present new observations of the C18O 3-2 and the HDO 3_{12}-2_{21} lines towards NGC 7129 FIRS 2. Combining these observations with Herschel data and modeling their emissions, we constrain the C18O and HDO abundance profiles across the protostellar envelope. In particular, we derive the abundance of C18O and HDO in the hot core. The intensities of the C18O lines are well reproduced assuming that the C18O abundance decreases through the protostellar envelope from the outer edge towards the centre until the point where the gas and dust reach the CO evaporation temperature (~20-25 K) where the C18O is released back to the gas phase. Once the C18O is released to the gas phase, the modelled C18O abundance is found to be ~1.6x10^{-8}, which is a factor of 10 lower than the reference abundance. This result is supported by the non-detection of C18O 9-8, which proves that even in the hot core (T_k>100 K) the CO abundance must be 10 times lower than the reference value. Several scenarios are discussed to explain this C18O deficiency. One possible explanation is that during the pre-stellar and protostellar phase, the CO is removed from the grain mantles by reactions to form more complex molecules. Our HDO modeling shows that the emission of HDO 3_{12}-2_{21} line is maser and comes from the hot core (T_k>100 K). Assuming the physical structure derived by Crimier et al. (2010), we determine a HDO abundance of ~0.4 - 1x10^{-7} in the hot core of this IM protostar, similar to that found in the hot corinos NGC 1333 IRAS 2A and IRAS 16293-2422.

During the embedded phase of pre-main sequence stellar evolution, a disk forms from the dense envelope while an accretion-driven outflow carves out a cavity within the envelope. Highly excited H2O emission in spatially unresolved Spitzer/IRS spectra of a low-mass Class 0 object, NGC 1333 IRAS 4B, has previously been attributed to the envelope-disk accretion shock but could instead be produced in an outflow. As part of the survey of low-mass sources in the Water in Star Forming Regions with Herschel (WISH-LM) program, we used Herschel/PACS to obtain a far-IR spectrum and several Nyquist-sampled spectral images with to determine the origin of excited H2O emission from NGC 1333 IRAS 4B. The spectrum has high signal-to-noise in a rich forest of H2O, CO, and OH lines, providing a near-complete census of far-IR molecular emission from a Class 0 protostar. The excitation diagrams for the three molecules all require fits with two excitation temperatures, indicating the presence of two physical components. The highly excited component of H2O emission is characterized by subthermal excitation of 1500 K gas with a density of 10^5 - 10^7 cm-3, conditions that also reproduce the mid-IR H2O emission detected by Spitzer. On the other hand, a high density, low temperature gas can reproduce the H2O spectrum observed by Spitzer but underpredicts the H2O lines seen by Herschel. Nyquist-sampled spectral maps of several lines show two spatial components of H2O emission, one centered at 1200 AU south of the central source at the position of the blueshifted outflow lobe and a second centered on-source. Both spatial components of the far-IR H2O emission are consistent with emission from the outflow. The gas cooling from the IRAS 4B envelope cavity walls is dominated by far-IR H2O emission, in contrast to stronger [O I] and CO cooling from more evolved protostars. [one sentence truncated]

Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eup=4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, UV irradiation of the outflow cavity walls, and C-type shocks along the cavity walls. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases and physical characteristics. Results. The passively heated gas in the envelope accounts for 3-10% of the CO luminosity summed over all rotational lines up to J=40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.

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

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.

(Abridged) As the stellar X-ray and UV light penetration of a protoplanetary disk depends sensitively on the dust properties, trace molecular species like HCO+, HCN, and CN are expected to show marked differences from photoprocessing effects as the dust content in the disk evolves. We investigate the evolution of the UV irradiation of the molecular gas in a sample of classical T Tauri stars in Taurus that exhibit a wide range in grain growth and dust settling properties. We obtained HCO+ (J=3-2), HCN (J=3-2), and CN (J=2-1) observations of 13 sources with the JCMT. Our sample has 1.3mm fluxes in excess of 75mJy, indicating the presence of significant dust reservoirs; a range of dust settling as traced through their spectral slopes between 6, 13, and 25 microns; and varying degrees of grain growth as extrapolated from the strength of their 10-micron silicate emission features. We compare the emission line strengths with the sources' continuum flux and infrared features, and use detailed modeling based on two different model prescriptions to compare typical disk abundances for HCO+, HCN, and CN with the gas-line observations for our sample. We detected HCO+ (3-2) toward 6 disks, HCN (3-2) from 0 disks, and CN (2-1) toward 4 disks. For the complete sample, there is no correlation between the gas-line strengths or their ratios and either the sources' dust continuum flux or infrared slopes.

During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow phenomenon plays an important role in the regulation of star formation. Understanding the associated shocks and their energetic effects is therefore essential to the study of star formation. We present comparisons of shock models with observations of H$_2$ and SiO emission in the bipolar outflow BHR71, and predict water emission, under the basic assumption that the emission regions of the considered species coincide, at the resolution of currently available observations. New SiO observations from APEX are presented, and combined with \textit{Spitzer} and ground-based observations of H$_2$ to constrain shock models. The shock regions associated with targeted positions in the molecular outflow are studied by means of a 1D code that generates models of the propagation of stationary shock waves, and approximations to non-stationary ones. The SiO emission in the inner part of the outflow is concentrated near the apex of the corresponding bow-shock that is also seen in the pure rotational transitions of H$_2$. Simultaneous modelling is possible for H$_2$ and SiO and leads to constraints on the silicon pre-shock distribution on the grain mantles and/or cores. The best-fitting models are found to be of the non-stationary type, but the degeneracy of the solutions is still large. Pre-shock densities of 10$^4$ and 10$^5$ cm$^{-3}$ are investigated, and the associated best-model candidates have rather low velocity (respectively, 20-30 and 10-15 km s$^{-1}$) and are not older than 1000 years. We provide emission predictions for water, focusing on the brightest transitions, to be observed with the PACS and HIFI instruments of the \textit{Herschel} Telescope.

Gas-phase complex organic molecules are commonly detected toward high-mass protostellar hot cores. Detections toward low-mass protostars and outflows are comparatively rare, and a larger sample is key to investigate how the chemistry responds to its environment. Guided by the prediction that complex organic molecules form in CH3OH-rich ices and thermally or non-thermally evaporate with CH3OH, we have identified three sight-lines in the Serpens core - SMM1, SMM4 and SMM4-W - which are likely to be rich in complex organics. Using the IRAM 30m telescope, narrow lines (FWHM of 1-2 km s-1) of CH3CHO and CH3OCH3 are detected toward all sources, HCOOCH3 toward SMM1 and SMM4-W, and C2H5OH not at all. Beam-averaged abundances of individual complex organics range between 0.6 and 10% with respect to CH3OH when the CH3OH rotational temperature is applied. The summed complex organic abundances also vary by an order of magnitude, with the richest chemistry toward the most luminous protostar SMM1. The range of abundances compare well with other beam-averaged observations of low-mass sources. Complex organic abundances are of the same order of magnitude toward low-mass protostars and high-mass hot cores, but HCOOCH3 is relatively more important toward low-mass protostars. This is consistent with a sequential ice photochemistry, dominated by CHO-containing products at low temperatures and early times.

OH is an important molecule in the H2O chemistry and the cooling budget of star-forming regions. The goal of the Herschel key program `Water in Star-forming regions with Herschel' (WISH) is to study H2O and related species during protostellar evolution. Our aim in this letter is to assess the origin of the OH emission from star-forming regions and constrain the properties of the emitting gas. High-resolution observations of the OH 2Pi1/2 J = 3/2-1/2 triplet at 1837.8 GHz (163.1 micron) towards the high-mass star-forming region W3 IRS 5 with the Heterodyne Instrument for the Far-Infrared (HIFI) on Herschel reveal the first hyperfine velocity-resolved OH far-infrared spectrum of a star-forming region. The line profile of the OH emission shows two components: a narrow component (FWHM approx. 4-5 km/s) with partially resolved hyperfine structure resides on top of a broad (FWHM approx. 30 km/s) component. The narrow emission agrees well with results from radiative transfer calculations of a spherical envelope model for W3 IRS 5 with a constant OH abundance of approx. 8e-9. Comparison with H2O yields OH/H2O abundance ratios of around 1e-3 for T > 100 K and around unity for T < 100K, consistent with the current picture of the dense cloud chemistry with freeze-out and photodesorption. The broad component is attributed to outflow emission. An abundance ratio of OH/H2O > 0.028 in the outflow is derived from comparison with results of water line modeling. This ratio can be explained by a fast J-type shock or a slower UV-irradiated C-type shock.

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.

The aim of this study is to investigate the distribution and excitation of the methanol gas around Cep A and to describe the physical conditions of the region. In addition the large scale abundance distribution is determined in order to understand the morphology and kinematics of star forming regions in which methanol masers occur. The spatial distribution of the methanol is studied by mapping line emission, as well as the column density and excitation temperature, which are estimated using rotation diagrams. For a limited number of positions the parameters are checked with non-LTE models. Furthermore, the distribution of the methanol abundance is derived in comparison with archival dust continuum maps. Methanol is detected over a 0.3x0.15 pc area centred on the Cep A HW2 source, showing an outflow signature. Most of the gas can be characterized by a moderately warm rotation temperature (30-60K). At the central position two velocity components are detected with different excitation characteristics, the first related to the large-scale outflow. The second component, uniquely detected at the central location, is probably associated with the maser emission on much smaller scales of 2". Detailed analysis reveals that the highest densities and temperatures occur for these inner components. In the inner region the dust and gas are shown to have different physical parameters. Abundances of methanol in the range 10E-9 - 10E-7 are inferred, with the abundance peaking at the maser position. The geometry of the large-scale methanol is in accordance with previous determinations of the Cep A geometry, in particular those from methanol masers. The dynamical and chemical time-scales are consistent with the methanol originating from a single driving source associated with the HW2 object and the masers in its equatorial region.

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 in NGC1333. The observations show several energetic components including shocked and quiescent gas. Radiative transfer models are used to quantify the C18O envelope abundance which require a jump in the abundance at an evaporation temperature, T_ev ~25 K, providing new direct evidence of a CO ice evaporation zone around protostars. The abundance in the outermost part of the envelope, X_0, is within the canonical value of 2x10^-4; however the inner abundance, X_in, is found around a factor of 3-5 lower than X_0.

Although deuterium enrichment of water may provide an essential piece of information in the understanding of the formation of comets and protoplanetary systems, only a few studies up to now have aimed at deriving the HDO/H2O ratio in low-mass star forming regions. Previous studies of the molecular deuteration toward the solar-type class 0 protostar, IRAS 16293-2422, have shown that the D/H ratio of water is significantly lower than other grain-surface-formed molecules. It is not clear if this property is general or particular to this source. In order to see if the results toward IRAS 16293-2422 are particular, we aimed at studying water deuterium fractionation in a second low-mass solar-type protostar, NGC1333-IRAS2A. Using the 1-D radiative transfer code RATRAN, we analyzed five HDO transitions observed with the IRAM 30m, JCMT, and APEX telescopes. We assumed that the abundance profile of HDO in the envelope is a step function, with two different values in the inner warm (T>100 K) and outer cold (T<100 K) regions of the protostellar envelope. The inner and outer abundance of HDO is found to be well constrained at the 3 sigma level. The obtained HDO inner and outer fractional abundances are x_in=6.6e-8 - 1e-7 and x_out=9e-11 - 1.8e-9 (3 sigma). These values are close to those in IRAS 16293-2422, which suggests that HDO may be formed by the same mechanisms in these two solar-type protostars. Taking into account the (rather poorly constrained) H2O abundance profile deduced from Herschel observations, the derived HDO/H2O in the inner envelope is larger than 1% and in the outer envelope it is 0.9%-18%. These values are more than one order of magnitude higher than what is measured in comets. If the same ratios apply to the protosolar nebula, this would imply that there is some efficient reprocessing of the material between the protostellar and cometary phases. The H2O inner fractional [...]

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

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.

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.

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.

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.

Aims: We employ archival Spitzer slit-scan observations of the HH211 outflow in order to investigate its warm gas content, assess the jet mass flux in the form of H2 and probe for the existence of an embedded atomic jet. Methods: Detected molecular and atomic lines are interpreted by means of emission line diagnostics and an existing grid of molecular shock models. The physical properties of the warm gas are compared against other molecular jet tracers and to the results of a similar study towards the L1448-C outflow. Results: We have detected and mapped the v=0-0 S(0) - S(7) H2 lines and fine-structure lines of S, Fe+, and Si+. H2 is detected down to 5" from the source and is characterized by a "cool" T~300K and a "warm" T~1000 K component, with an extinction Av ~ 8 mag. The amount of cool H2 towards the jet agrees with that estimated from CO assuming fully molecular gas. The warm component is well fitted by C-type shocks with a low beam filling factor ~ 0.01-0.04 and a mass-flux similar to the cool H2. The fine-structure line emission arises from dense gas with ionization fraction ~0.5 - 5 x 10e-3, suggestive of dissociative shocks. Line ratios to sulfur indicate that iron and silicon are depleted compared to solar abundances by a factor ~10-50. Conclusions: Spitzer spectral mapping observations reveal for the first time a cool H$_2$ component towards the CO jet of HH211 consistent with the CO material being fully molecular and warm at ~ 300 K. The maps also reveal for the first time the existence of an embedded atomic jet in the HH211 outflow that can be traced down to the central source position. Its significant iron and silicon depletion excludes an origin from within the dust sublimation zone around the protostar. The momentum-flux seems insufficient to entrain the CO jet, although current uncertainties on jet speed and shock conditions are too large for a definite conclusion.

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.

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

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

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

H2 formation is an important process in post-shock regions, since H2 is an active participant in the cooling and shielding of the environment. The onset of H2 formation therefore has a strong effect on the temperature and chemical evolution in the post shock regions. We recently developed a model for H2 formation on a graphite surface in warm conditions. The graphite surface acts as a model system for grains containing large areas of polycyclic aromatic hydrocarbon structures. Here this model is used to obtain a new description of the H2 formation rate as a function of gas temperature that can be implemented in molecular shock models. The H2 formation rate is substantially higher at high gas temperatures as compared to the original implementation of this rate in shock models, because of the introduction of H atoms which are chemically bonded to the grain (chemisorption). Since H2 plays such a key role in the cooling, the increased rate is found to have a substantial effect on the predicted line fluxes of an important coolant in dissociative shocks [O I] at 63.2 and 145.5 micron. With the new model a better agreement between model and observations is obtained. Since one of the goals of Herschel/PACS will be to observe these lines with higher spatial resolution and sensitivity than the former observations by ISO-LWS, this more accurate model is very timely to help with the interpretation of these future results.

Observations of Serpens have been performed at the JCMT using Harp-B. Maps over a 4.5'x5.4' region were made in a frequency window around 338 GHz, covering the 7-6 transitions of methanol. Emission is extended over each source, following the column density of H2 but showing up also particularly strongly around outflows. The rotational temperature is low, 15-20 K, and does not vary with position within each source. The abundance is typically 10^-9 - 10^-8 with respect to H2 in the outer envelope, whereas "jumps" by factors of up to 10^2 -10^3 inside the region where the dust temperature exceeds 100 K are not excluded. A factor of up to ~ 10^3 enhancement is seen in outflow gas. In one object, SMM4, the ice abundance has been measured to be ~ 3x10^-5 with respect to H2 in the outer envelope, i.e., a factor of 10^3 larger than the gas-phase abundance. Comparison with C18O J=3-2 emission shows that strong CO depletion leads to a high gas-phase abundance of CH3OH not just for the Serpens sources, but for a larger sample of protostars. The observations illustrate the large-scale, low-level desorption of CH3OH from dust grains, extending out to and beyond 7500 AU from each source, a scenario which is consistent with non-thermal (photo-)desorption from the ice. The observations also illustrate the usefulness of CH3OH as a tracer of energetic input in the form of outflows, where methanol is sputtered from the grain surfaces. Finally, the observations provide further evidence of CH3OH formation through CO hydrogenation proceeding on grain surfaces in low-mass envelopes.

We present near-infrared imaging and spectroscopic high spatial resolution observations of the SMC region N88 containing the bright, excited, extincted and compact H II region N88A of size ~ 1 pc. To investigate its stellar content and reddening, N88 was observed using spectroscopy and imagery in the JHKs- and L'-band at a spatial resolution of ~ 0.1-0.3", using the VLT UT4 equipped with the NAOS adaptive optics system. In order to attempt to establish if the origin of the infra-red (IR) excess is due to bright nebulosity, circumstellar material and/or local dust, we used Ks vs J-K colour-magnitude (CM) and JHK colour-colour (CC) diagrams, as well as L' imagery.Our IR-data reveal in the N88 area an IR-excess fraction of geq 30 per cent of the detected stars,as well as an unprecedently detailed morphology of N88A. It consists of an embedded cluster of ~3.5" (~ 1 pc) in diameter, of at least thirteen resolved stars superposed with an unusual bright continuum centered on a very bright star. The four brightest stars in this cluster lie red-ward of H-K geq 0.45 mag, and could be classified as young stellar object (YSO) candidates. Four other probable YSO candidates are also detected in N88 along a south-north bow-shaped thin H2 filament at ~ 7" east of the young central bright star. At 0.2" east of this star, a heavily embedded core is detected that could be a massive class I protostar candidate. The 2.12 mu H2 image of N88A resembles a shell of diameter ~ 3" ~ 0.9 pc) centered on the bright star. The line ratios of H2 2-1 S(1) and 1-0 S(0) relative to 1-0 S(1), as well as the presence of high v lines, are indicative of photodissociation regions, rather than shocks.

The origin and heating mechanisms of warm (50<T<200 K) molecular gas in low-mass young stellar objects (YSOs) are strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by shocks and by UV associated with the outflows or accretion have been proposed. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, bipolar outflows and the quiescent cloud. High-J CO maps (12CO J=6--5 and 7--6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of the 13CO J=6--5 transition and [C I] 3P_2-3P_1 line were taken. Strong quiescent narrow-line 12CO 6--5 and 7--6 emission is seen toward all protostars. In the case of HH~46 and Ced 110 IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. Shock-heated warm gas is only detected for Class 0 flows and the more massive Class I sources such as HH~46. Outflow temperatures, estimated from the CO 6--5 and 3--2 line wings, are ~100 K, close to model predictions, with the exception of the L~1551 IRS 5 and IRAS 12496-7650, for which temperatures <50 K are found. APEX-CHAMP+ is uniquely suited to directly probe a protostar's feedback on its accreting envelope gas in terms of heating, photodissociation, and outflow dispersal by mapping 1'x1' regions in high-J CO and [C I] lines.

AIMS Our aim is to characterize the size, mass, density and temperature profiles of the protostellar envelope of HH~46 IRS 1 and its surrounding cloud material as well as the effect the outflow has on its environment.METHODS The CHAMP+ and LABOCA arrays on the APEX telescope, combined with lower frequency line receivers, are used to obtain a large continuum map and smaller heterodyne maps in various isotopologues of CO and HCO+. The high-J lines of CO (6--5 and 7--6) and its isotopologues together with [C I] 2--1, observed with CHAMP+, are used to probe the warm molecular gas in the inner few hundred AU and in the outflowing gas. The data are interpreted with continuum and line radiative transfer models. RESULTS Broad outflow wings are seen in CO low- and high-J lines at several positions, constraining the gas temperatures to a constant value of ~100 K along the red outflow axis and to ~60 K for the blue outflow. The derived outflow mass is of order 0.4--0.8 M_sol, significantly higher than previously found. The bulk of the strong high-J CO line emission has a surprisingly narrow width, however, even at outflow positions. These lines cannot be fit by a passively heated model of the HH 46 IRS envelope. We propose that it originates from photon heating of the outflow cavity walls by ultraviolet photons originating in outflow shocks and the accretion disk boundary layers. At the position of the bow shock itself, the UV photons are energetic enough to dissociate CO. The envelope mass of ~5 M_sol is strongly concentrated towards HH 46 IRS with a density power law of -1.8.

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