Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.

Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences.

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time.

[abridged]CO rovibrational lines are efficient probes of warm molecular gas and can give unique insights into the inner 10 AU of proto-planetary disks. Recent studies have found a relation between the ratio of lines originating from the second and first vibrationally excited state, denoted as $v2/v1$, and the emitting radius of CO. In disks around Herbig Ae stars the vibrational excitation is low when CO lines come from close to the star, and high when lines only probe gas at large radii (more than 5 AU). We aim to find explanations for the observed trends between CO vibrational ratio, emitting radii, and NIR excess, and identify their implications in terms of the physical and chemical structure of inner disks around Herbig stars. Slab models and full disk thermo chemical models are calculated. Simulated observations from the models are directly compared to the data. Broad CO lines with low vibrational ratios are best explained by a warm (400-1300 K) inner disk surface with gas-to-dust ratios below 1000; no CO is detected within/at the inner dust rim, due to dissociation at high temperatures. In contrast, explaining the narrow lines with high vibrational ratios requires an inner cavity of a least 5 AU in both dust and gas, followed by a cool (100-300 K) molecular gas reservoir with gas-to-dust ratios greater than 10000 at the cavity wall. In all cases the CO gas must be close to thermalization with the dust. The high gas-to-dust ratios needed to explain high $v2/v1$ in narrow CO lines for a subset of group I disks can naturally be interpreted as due to the dust traps that have been proposed to explain millimeter dust cavities. The broad lines seen in most group II objects indicate a very flat disk in addition to inner disk substructures within 10 AU that can be related to the substructures recently observed with ALMA.

The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.

The extent of the gas in protoplanetary disks is observed to be universally larger than the extent of the dust. This is often attributed to radial drift and grain growth of the mm grains, but line optical depth produces a similar observational signature. We investigate in what parts of the disk structure parameter space dust evolution and line optical depth are the dominant drivers of the observed gas and dust size difference. Using the thermochemical model DALI with dust evolution included we ran a grid of models aimed at reproducing the observed gas and dust size dichotomy. The relation between R_dust and dust evolution is non-monotonic and depends on the disk structure. R_gas is directly related to the radius where the CO column density drops below 10^15 cm^-2 and CO becomes photodissociated. R_gas is not affected by dust evolution but scales with the total CO content of the disk. R_gas/R_dust > 4 is a clear sign for dust evolution and radial drift in disks, but these cases are rare in current observations. For disks with a smaller R_gas/R_dust, identifying dust evolution from R_gas/R_dust requires modelling the disk structure including the total CO content. To minimize the uncertainties due to observational factors requires FWHM_beam < 1x the characteristic radius and a peak SNR > 10 on the 12CO emission moment zero map. For the dust outer radius to enclose most of the disk mass, it should be defined using a high fraction (90-95%) of the total flux. For the gas, any radius enclosing > 60% of the 12CO flux will contain most of the disk mass. To distinguish radial drift and grain growth from line optical depth effects based on size ratios requires disks to be observed at high enough angular resolution and the disk structure should to be modelled to account for the total CO content of the disk.

The formation of planets occurs within protoplanetary disks surrounding young stars, resulting in perturbation of the gas and dust surface densities. Here, we report the first evidence of spatially resolved gas surface density ($\Sigma_{g}$) perturbation towards the AS~209 protoplanetary disk from the optically thin C$^{18}$O ($J=2-1$) emission. The observations were carried out at 1.3~mm with ALMA at a spatial resolution of about 0.3$\arcsec$ $\times$ 0.2$\arcsec$ (corresponding to $\sim$ 38 $\times$ 25 au). The C$^{18}$O emission shows a compact ($\le$60~au), centrally peaked emission and an outer ring peaking at 140~au, consistent with that observed in the continuum emission and, its azimuthally averaged radial intensity profile presents a deficit that is spatially coincident with the previously reported dust map. This deficit can only be reproduced with our physico-thermochemical disk model by lowering $\Sigma_{gas}$ by nearly an order of magnitude in the dust gaps. Another salient result is that contrary to C$^{18}$O, the DCO$^{+}$ ($J=3-2$) emission peaks between the two dust gaps. We infer that the best scenario to explain our observations (C$^{18}$O deficit and DCO$^{+}$ enhancement) is a gas perturbation due to forming-planet(s), that is commensurate with previous continuum observations of the source along with hydrodynamical simulations. Our findings confirm that the previously observed dust gaps are very likely due to perturbation of the gas surface density that is induced by a planet of at least 0.2~M$\rm_{Jupiter}$ in formation. Finally, our observations also show the potential of using CO isotopologues to probe the presence of saturn mass planet(s).

High resolution ALMA observations of protoplanetary disks have revealed that many, if not all primordial disks consist of ring-like dust structures. The origin of these dust rings remains unclear, but a common explanation is the presence of planetary companions that have cleared gaps along their orbit and trapped the dust at the gap edge. A signature of this scenario is a decrease of gas density inside these gaps. In recent work, Isella et al. 2016 derived drops in gas density consistent with Saturn-mass planets inside the gaps in the HD163296 disk through spatially resolved CO isotopologue observations. However, as CO abundance and temperature depends on a large range of factors, the interpretation of CO emission is non-trivial. We use the physical-chemical code DALI to show that the gas temperature increases inside dust density gaps, implying that any gaps in the gas, if present, would have to be much deeper, consistent with planet masses higher than a Jupiter mass. Furthermore, we show that a model with increased grain growth at certain radii, as expected at a snowline, can reproduce the dust rings in HD163296 equally well without the need for companions. This scenario can explain both younger and older disks with observed gaps, as gaps have been seen in systems as young <1 Myr. While the origin of the rings in HD163296 remains unclear, these modeling results demonstrate that care has to be taken when interpreting CO emission in protoplanetary disk observations.

It is key to constrain the gas surface density distribution, Sigma_gas, as function of disk radius in protoplanetary disks. In this work we investigate if spatially resolved observations of rarer CO isotopologues may be good tracers of Sigma_gas. Physical-chemical models with different input Sigma_gas(R) are run. The input disk surface density profiles are compared with the simulated 13CO intensity radial profiles to check if and where the two follow each other. There is always an intermediate region in the disk where the slope of the 13CO radial emission profile and Sigma_gas(R) coincide. At small radii the line radial profile underestimates Sigma_gas, as 13CO emission becomes optically thick. The same happens at large radii where the column densities become too low and 13CO is not able to efficiently self-shield. If the gas surface density profile is a simple power-law of the radius, the input power-law index can be retrieved within 20% uncertainty if one choses the proper radial range. If instead Sigma_gas(R) follows the self-similar solution for a viscously evolving disk, retrieving the input power-law index becomes challenging, in particular for small disks. Nevertheless, it is found that the power-law index can be in any case reliably fitted at a given line intensity contour around 6 K km/s, and this produces a practical method to constrain the slope of Sigma_gas(R). Application of such a method is shown in the case study of the TW Hya disk. Spatially resolved 13CO line radial profiles are promising to probe the disk surface density distribution, as they directly trace Sigma_gas(R)profile at radii well resolvable by ALMA. There, chemical processes like freeze-out and isotope selective photodissociation do not affect the emission, and, assuming that the volatile carbon does not change with radius, no chemical model is needed when interpreting the observations.

Aims: The two stable isotopes of nitrogen, 14N and 15N, exhibit a range of abundance ratios both inside and outside the solar system. The elemental ratio in the solar neighborhood is 440. Recent ALMA observations showed HCN/HC15N ratios from 83 to 156 in six T Tauri and Herbig disks and a CN/C15 N ratio of 323 +/- 30 in one T Tauri star. We aim to determine the dominant mechanism responsible for these enhancements of 15N: low-temperature exchange reactions or isotope-selective photodissociation of N2. Methods: Using the thermochemical code DALI, we model the nitrogen isotope chemistry in circumstellar disks with a 2D axisymmetric geometry. Our chemical network is the first to include both fractionation mechanisms for nitrogen. The model produces abundance profiles and isotope ratios for several key N-bearing species. We study how these isotope ratios depend on various disk parameters. Results: The formation of CN and HCN is closely coupled to the vibrational excitation of H2 in the UV-irradiated surface layers of the disk. Isotope fractionation is completely dominated by isotope-selective photodissociation of N2. The column density ratio of HCN over HC15N in the disk's inner 100 au does not depend strongly on the disk mass, the flaring angle or the stellar spectrum, but it is sensitive to the grain size distribution. For larger grains, self-shielding of N2 becomes more important relative to dust extinction, leading to stronger isotope fractionation. Between disk radii of ~50 and 200 au, the models predict HCN/HC15N and CN/C15N abundance ratios consistent with observations of disks and comets. The HCN/HC15N and CN/C15N column density ratios in the models are a factor of 2-3 higher than those inferred from the ALMA observations.

The oxygen-bearing molecular ions OH+, H2O+, and H3O+ are key species that probe the ionization rate of (partially) molecular gas that is ionized by X-rays and cosmic rays permeating the interstellar medium. We report Herschel far-infrared and submillimeter spectroscopic observations of OH+ in Mrk 231, showing both ground-state P-Cygni profiles, and excited line profiles with blueshifted absorption wings extending up to ~1000 km s^{-1}. In addition, OH+ probes an excited component peaking at central velocities, likely arising from the torus probed by the OH centimeter-wave megamaser. Four lines of H2O+ are also detected at systemic velocities, but H3O+ is undetected. Based on our earlier OH studies, we estimate an abundance ratio of OH/OH+~5-10 for the outflowing components and ~20 for the torus, and an OH+ abundance relative to H nuclei of ~>10^{-7}. We also find high OH+/H2O+ and OH+/H3O+ ratios, both are ~>4 in the torus and ~>10-20 in the outflowing gas components. Chemical models indicate that these high OH+ abundances relative to OH, H2O+, and H3O+ are characteristic of gas with a high ionization rate per unit density, \zeta/n_H~(1-5)x10^{-17} cm^3 s^{-1} and ~(1-2)x10^{-16} cm^3 s^{-1} for the above components, respectively, and an ionization rate of \zeta~(0.5-2)x10^{-12} s^{-1}. X-rays appear to be unable to explain the inferred ionization rate, and thus we suggest that low-energy (10-400 MeV) cosmic-rays are primarily responsible for the ionization with \dot{M}_{CR}~0.01 M_{sun} yr^{-1} and \dot{E}_{CR}~10^{44} erg s^{-1}, the latter corresponding to 1% of the AGN luminosity and similar to the energetics of the molecular outflow. We suggest that cosmic-rays accelerated in the forward shock associated with the molecular outflow are responsible for the ionization, as they diffuse through the outflowing molecular phase downstream.

This work aims to understand which midplane conditions are probed by the DCO$^+$ emission in the disk around the Herbig Ae star HD 169142. We explore the sensitivity of the DCO$^+$ formation pathways to the gas temperature and the CO abundance. The DCO$^+$ $J$=3-2 transition was observed with ALMA at a spatial resolution of 0.3". The HD 169142 DCO$^+$ radial intensity profile reveals a warm, inner component at radii <30 AU and a broad, ring-like structure from ~50-230 AU with a peak at 100 AU just beyond the millimeter grain edge. We modeled DCO$^+$ emission in HD 169142 with a physical disk structure adapted from the literature, and employed a simple deuterium chemical network to investigate the formation of DCO$^+$ through the cold deuterium fractionation pathway via H$_2$D$^+$. Contributions from the warm deuterium fractionation pathway via CH$_2$D$^+$ are approximated using a constant abundance in the intermediate disk layers. Parameterized models show that alterations to the midplane gas temperature and CO abundance of the literature model are both needed to recover the observed DCO$^+$ radial intensity profile. The best-fit model contains a shadowed, cold midplane in the region z/r < 0.1 with an 8 K decrease in gas temperature and a factor of five CO depletion just beyond the millimeter grain edge, and a 2 K decrease in gas temperature for r > 120 AU. The warm deuterium fractionation pathway is implemented as a constant DCO$^+$ abundance of 2.0$\times$10$^{-12}$ between 30-70 K. The DCO$^+$ emission probes a reservoir of cold material in the HD 169142 outer disk that is not revealed by the millimeter continuum, the SED, nor the emission from the 12CO, 13CO, or C18O $J$=2-1 lines.

Protoplanetary discs around brown dwarfs and very low mass stars offer some of the best prospects for forming Earth-sized planets in their habitable zones. To this end, we study the nature of the disc around the very low mass star V410 X-ray 1, whose SED is indicative of an optically thick and very truncated dust disc, with our modelling suggesting an outer radius of only 0.6 au. We investigate two scenarios that could lead to such a truncation, and find that the observed SED is compatible with both. The first scenario involves the truncation of both the dust and gas in the disc, perhaps due to a previous dynamical interaction or the presence of an undetected companion. The second scenario involves the fact that a radial location of 0.6 au is close to the expected location of the H$_2$O snowline in the disc. As such, a combination of efficient dust growth, radial migration, and subsequent fragmentation within the snowline leads to an optically thick inner dust disc and larger, optically thin outer dust disc. We find that a firm measurement of the CO $J=2$--1 line flux would enable us to distinguish between these two scenarios, by enabling a measurement of the radial extent of gas in the disc. Many models we consider contain at least several Earth-masses of dust interior to 0.6 au, suggesting that V410 X-ray 1 could be a precursor to a system with tightly-packed inner planets, such as TRAPPIST-1.

Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities hinting at recently formed, giant planets. However, many of these studies are biased towards the brightest disks in the nearby star forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (>20 AU radius) from a complete disk survey of the Lupus star forming region, using ALMA Band 7 observations at 0.3" (22-30 AU radius) resolution of the 345 GHz continuum, 13CO and C18O 3-2 observations and the Spectral Energy Distribution of each source. Gas and dust surface density profiles are derived using the physical-chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star forming region. The dust cavity sizes range from 20-90 AU radius and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is ~11%, which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities.

Bright ring-like structure emission of the CN molecule has been observed in protoplanetary disks. We investigate whether such structures are due to the morphology of the disk itself or if they are instead an intrinsic feature of CN emission. With the intention of using CN as a diagnostic, we also address to which physical and chemical parameters CN is most sensitive. A set of disk models were run for different stellar spectra, masses, and physical structures via the 2D thermochemical code DALI. An updated chemical network that accounts for the most relevant CN reactions was adopted. Ring-shaped emission is found to be a common feature of all adopted models; the highest abundance is found in the upper outer regions of the disk, and the column density peaks at 30-100 AU for T Tauri stars with standard accretion rates. Higher mass disks generally show brighter CN. Higher UV fields, such as those appropriate for T Tauri stars with high accretion rates or for Herbig Ae stars or for higher disk flaring, generally result in brighter and larger rings. These trends are due to the main formation paths of CN, which all start with vibrationally excited H2* molecules, that are produced through far ultraviolet (FUV) pumping of H2. The model results compare well with observed disk-integrated CN fluxes and the observed location of the CN ring for the TW Hya disk. CN rings are produced naturally in protoplanetary disks and do not require a specific underlying disk structure such as a dust cavity or gap. The strong link between FUV flux and CN emission can provide critical information regarding the vertical structure of the disk and the distribution of dust grains which affects the UV penetration, and could help to break some degeneracies in the SED fitting. In contrast with C2H or c-C3H2, the CN flux is not very sensitive to carbon and oxygen depletion.

The paper presents new high angular resolution ALMA 1.3 mm dust continuum observations of the protoplanetary system AS 209 in the Ophiuchus star forming region. The dust continuum emission is characterized by a main central core and two prominent rings at $r = 75\,$au and $r = 130\,$au intervaled by two gaps at at $r = 62\,$au and $r = 103\,$au. The two gaps have different widths and depths, with the inner one being narrower and shallower. We determined the surface density of the millimeter dust grains using the 3D radiative transfer disk code \textsc{dali}. According to our fiducial model the inner gap is partially filled with millimeter grains while the outer gap is largely devoid of dust. The inferred surface density is compared to 3D hydrodynamical simulations (FARGO-3D) of planet-disk interaction. The outer dust gap is consistent with the presence of a giant planet ($M_{\rm planet} \sim 0.8\,M_{\rm Staturn}$); the planet is responsible for the gap opening and for the pile-up of dust at the outer edge of the planet orbit. The simulations also show that the same planet can give origin to the inner gap at $r = 62\,$au. The relative position of the two dust gaps is close to the 2:1 resonance and we have investigated the possibility of a second planet inside the inner gap. The resulting surface density (including location, width and depth of the two dust gaps) are in agreement with the observations. The properties of the inner gap pose a strong constraint to the mass of the inner planet ($M_{\rm planet} < 0.1\,M_{\rm J}$). In both scenarios (single or pair of planets), the hydrodynamical simulations suggest a very low disk viscosity ($\alpha < 10^{-4}$). Given the young age of the system (0.5 - 1 Myr), this result implies that the formation of giant planets occurs on a timescale of $\lesssim$ 1\,Myr.

ALMA observations of protoplanetary disks confirm earlier indications that there is a clear difference between the dust and gas radial extents. The origin of this difference is still debated, with both radial drift of the dust and optical depth effects suggested in the literature. In this work, the feedback of realistic dust particle distributions onto the gas chemistry and molecular emissivity is investigated, with a particular focus on CO isotopologues. The radial dust grain size distribution is determined using dust evolution models that include growth, fragmentation and radial drift. A new version of the code DALI is used to take into account how dust surface area and density influence the disk thermal structure, molecular abundances and excitation. The difference of dust and gas radial sizes is largely due to differences in the optical depth of CO lines and millimeter continuum, without the need to invoke radial drift. The effect of radial drift is primarily visible in the sharp outer edge of the continuum intensity profile. The gas outer radius probed by $^{12}$CO emission can easily differ by a factor of $\sim 2$ between the models for a turbulent $\alpha$ ranging between typical values. Grain growth and settling concur in thermally decoupling the gas and dust components, due to the low collision rate with large grains. As a result, the gas can be much colder than the dust at intermediate heights, reducing the CO excitation and emission, especially for low turbulence values. Also, due to disk mid-plane shadowing, a second CO thermal desorption (rather than photodesorption) front can occur in the warmer outer mid-plane disk. The models are compared to ALMA observations of HD 163296 as a test case. In order to reproduce the observed CO snowline of the system, a binding energy for CO typical of ice mixtures needs to be used rather than the lower pure CO value.

Protoplanetary disks around young stars are the sites of planet formation. While the dust mass can be estimated using standard methods, determining the gas mass - and thus the amount of material available to form giant planets - has proven to be very difficult. Hydrogen deuteride (HD) is a promising alternative to the commonly-used gas mass tracer, CO. We aim to examine the robustness of HD as tracer of the disk gas mass, specifically the effect of gas mass on the HD FIR emission and its sensitivity to the vertical structure. Deuterium chemistry reactions relevant for HD were implemented in the thermochemical code DALI and models were run for a range of disk masses and vertical structures. The HD J=1-0 line intensity depends directly on the gas mass through a sublinear power law relation with a slope of ~0.8. Assuming no prior knowledge about the vertical structure of a disk and using only the HD 1-0 flux, gas masses can be estimated to within a factor of 2 for low mass disks (M$_{\rm disk} < 10^{-3}$ M$_\odot$). For more massive disks, this uncertainty increases to more than an order of magnitude. Adding the HD 2-1 line or independent information about the vertical structure can reduce this uncertainty to a factor of ~3 for all disk masses. For TW Hya, using the radial and vertical structure from Kama et al. 2016b the observations constrain the gas mass to $6\cdot10^{-3}$ M$_\odot$ < M$_{\rm disk} < 9\cdot10^{-3}$ M$_\odot$. Future observations require a 5$\sigma$ sensitivity of $1.8\cdot10^{-20}$ W m$^{-2}$ ($2.5\cdot10^{-20}$ W m$^{-2}$) and a spectral resolving power R > 300 (1000) to detect HD 1-0 (HD 2-1) for all disk masses above $10^{-5}$ M$_\odot$ with a line-to-continuum ratio > 0.01. These results show that HD can be used as an independent gas mass tracer with a relatively low uncertainty and should be considered as an important science goal for future FIR missions.

We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are generally much lower than predictions of thermo-chemical models with canonical inputs. One ground-state transition is newly detected in the stacked spectrum of AA Tau, DM Tau, LkCa 15, and MWC 480. We run a grid of models to show that the abundance of gas-phase oxygen needs to be reduced by a factor of at least ~100 to be consistent with the observational upper limits (and positive detections) if a dust-to-gas mass ratio of 0.01 were to be assumed. As a continuation of previous ideas, we propose that the underlying reason for the depletion of oxygen (hence the low detection rate) is the freeze-out of volatiles such as water and CO onto dust grains followed by grain growth and settling/migration, which permanently removes these gas-phase molecules from the emissive upper layers of the outer disk. Such depletion of volatiles is likely ubiquitous among different disks, though not necessarily to the same degree. The volatiles might be returned back to the gas phase in the inner disk (within about 15 AU), which is consistent with current constraints. Comparison with studies on disk dispersal due to photoevaporation indicates that the timescale for volatile depletion is shorter than that of photoevaporation.

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. To study the interaction between the protoplanet and the disk itself observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3\,mm, $^{12}$CO, $^{13}$CO and C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 56 - 83\,$au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity ($R \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 180\,$au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution).

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and $^{12}$CO, $^{13}$CO, and C$^{18}$O J = 2-1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of ~0".25 (35 AU). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-$\mu$m-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius $\sim$70 AU that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 AU, or an edge that transitions smoothly over an annular width of 10 AU near 60 AU. In gas, the data are consistent with a cavity that is smaller, about 15 AU in radius, and whose surface density at 15 AU is $10^{3\pm1}$ times smaller than the surface density at 70 AU; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 AU or a double-drop model as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities, e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation.

[Abridged] The infrared ro-vibrational emission lines from organic molecules in the inner regions of protoplanetary disks are unique probes of the physical and chemical structure of planet forming regions and the processes that shape them. The non-LTE excitation effects of carbon dioxide (CO2) are studied in a full disk model to evaluate: (i) what the emitting regions of the different CO2 ro-vibrational bands are; (ii) how the CO2 abundance can be best traced using CO2 ro-vibrational lines using future JWST data and; (iii) what the excitation and abundances tell us about the inner disk physics and chemistry. CO2 is a major ice component and its abundance can potentially test models with migrating icy pebbles across the iceline. A full non-LTE CO2 excitation model has been built. The characteristics of the model are tested using non-LTE slab models. Subsequently the CO2 line formation has been modelled using a two-dimensional disk model representative of T-Tauri disks. The CO2 gas that emits in the 15 $\mu$m and 4.5 $\mu$m regions of the spectrum is not in LTE and arises in the upper layers of disks, pumped by infrared radiation. The v$_2$ 15 $\mu$m feature is dominated by optically thick emission for most of the models that fit the observations and increases linearly with source luminosity. Its narrowness compared with that of other molecules stems from a combination of the low rotational excitation temperature (~250 K) and the inherently narrower feature for CO2. The inferred CO2 abundances derived for observed disks are more than two orders of magnitude lower than those in interstellar ices (~10$^5$), similar to earlier LTE disk estimates. Line-to-continuum ratios are low, of order a few %, thus high signal-to-noise (S/N > 300) observations are needed for individual line detections. Prospects of accurate abundance retreival with JWST-MIRI and JWST-NIRSpec are discussed.

Despite intensive studies of protoplanetary disks, there is still no reliable way to determine their total mass and their surface density distribution, quantities that are crucial for describing both the structure and the evolution of disks up to the formation of planets. The goal of this work is to use less abundant CO isotopologues, whose detection is routine for ALMA, to infer the gas mass of disks. Isotope-selective effects need to be taken into account in the analysis, because they can significantly modify CO isotopologues line intensities. CO isotope-selective photodissociation has been implemented in the physical-chemical code DALI and 800 disk models have been run for a range of disk and stellar parameters. Dust and gas temperature structures have been computed self-consistently, together with a chemical calculation of the main species. Both disk structure and stellar parameters have been investigated. Total fluxes have been ray-traced for different CO isotopologues and for various transitions for different inclinations. A combination of 13CO and C18O total intensities allows inference of the total disk mass, although with non-negligible uncertainties. These can be overcome by employing spatially resolved observations, i.e. the disk's radial extent and inclination. Comparison with parametric models shows differences at the factor of a few level, especially for extremely low and high disk masses. Finally, total line intensities for different CO isotopologue and for various low-J transitions are provided and are fitted to simple formulae. The effects of a lower gas-phase carbon abundance and different gas/dust ratios are investigated as well, and comparison with other tracers is made. Disk masses can be determined within a factor of a few by comparing CO isotopologue lines observations with the simulated line fluxes, modulo the uncertainties in the volatile elemental abundances.

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $\leq0.1\,M_{\odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{\rm gas}=(0.1-1.5)\times 10^{-4}$. HD 100546 itself is a $\lambda\,$Bo\"{o}tis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{\rm gas}=(0.2-5.0)\times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{\rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems.

The gas temperature structure of protoplanetary disks is a key ingredient for interpreting various disk observations and for quantifying the subsequent evolution of these systems. The comparison of low- and mid-$J$ CO rotational lines is a powerful tool to assess the temperature gradient in the warm molecular layer of disks. Spectrally resolved high-$J$ ($J_{\rm u} > 14$) CO lines probe intermediate distances and heights from the star that are not sampled by (sub-)millimeter CO spectroscopy. This paper presents new {\it Herschel}/HIFI and archival PACS observations of $^{12}$CO, $^{13}$CO and \cii \ emission in 4 Herbig AeBe (HD 100546, HD 97048, IRS 48, HD 163296) and 3 T Tauri (AS 205, S CrA, TW Hya) disks. In the case of the T Tauri systems AS 205 and S CrA, the CO emission has a single-peaked profile, likely due to a slow wind. For all other systems, the {\it Herschel} CO spectra are consistent with pure disk emission and the spectrally-resolved lines (HIFI) and the CO rotational ladder (PACS) are analyzed simultaneously assuming power-law temperature and column density profiles, using the velocity profile to locate the emission in the disk. The temperature profile varies substantially from disk to disk. In particular, $T_{\rm gas}$ in the disk surface layers can differ by up to an order of magnitude among the 4 Herbig AeBe systems with HD 100546 being the hottest and HD 163296 the coldest disk of the sample. Clear evidence of a warm disk layer where $T_{\rm gas} > T_{\rm dust}$ is found in all the Herbig Ae disks. The observed CO fluxes and line profiles are compared to predictions of physical-chemical models. The primary parameters affecting the disk temperature structure are the flaring angle, the gas-to-dust mass ratio the scale height and the dust settling.

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.

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.

Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$\sigma$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $\lambda = 65.69\, \mu$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

Context: High-resolution studies of class 0 protostars represent the key to constraining protostar formation models. VLA16234-2417 represents the prototype of class 0 protostars, and it has been recently identified as a triple non-coeval system. Aim: We aim at deriving the physical properties of the jets in VLA16234-2417 using tracers of shocked gas. Methods: ALMA Cycle 0 Early Science observations of CO(2-1) in the extended configuration are presented in comparison with previous SMA CO(3-2) and Herschel-PACS [OI}] 63 micron observations. Gas morphology and kinematics were analysed to constrain the physical structure and origin of the protostellar outflows. Results: We reveal a collimated jet component associated with the [OI] 63 micron emission at about 8'' (about 960 AU) from source B. This newly detected jet component is inversely oriented with respect to the large-scale outflow driven by source A, and it is aligned with compact and fast jet emission very close to source B (about 0.3'') rather than with the direction perpendicular to the A disk. We also detect a cavity-like structure at low projected velocities, which surrounds the [OI] 63 micron emission and is possibly associated with the outflow driven by source A. Finally, no compact outflow emission is associated with source W. Conclusions: Our high-resolution ALMA observations seem to suggest there is a fast and collimated jet component associated with source B. This scenario would confirm that source B is younger than A, that it is in a very early stage of evolution, and that it drives a faster, more collimated, and more compact jet with respect to the large-scale slower outflow driven by A. However, a different scenario of a precessing jet driven by A cannot be firmly excluded from the present observations.

We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for GMC-GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region (PDR) based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter, and compare isolated versus colliding clouds. We find factors of ~2-3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on 13CO(J=2-1), 13CO(J=3-2), and 12CO(J=8-7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find the ratio of J=8-7 to lower-J emission is a powerful diagnostic probe of GMC collisions.

(Abridged*) Models of the young solar nebula assume a hot initial disk with most volatiles are in the gas phase. The question remains whether an actively accreting disk is warm enough to have gas-phase water up to 50 AU radius. No detailed studies have yet been performed on the extent of snowlines in an embedded accreting disk (Stage 0). Quantify the location of gas-phase volatiles in embedded actively accreting disk system. Two-dimensional physical and radiative transfer models have been used to calculate the temperature structure of embedded protostellar systems. Gas and ice abundances of H$_2$O, CO$_2$, and CO are calculated using the density-dependent thermal desorption formulation. The midplane water snowline increases from 3 to 55 AU for accretion rates through the disk onto the star between $10^{-9}$-$10^{-4} \ M_{\odot} \ {\rm yr^{-1}}$. CO$_2$ can remain in the solid phase within the disk for $\dot{M} \leq 10^{-5} \ M_{\odot} \ {\rm yr^{-1}}$ down to $\sim 20$ AU. Most of the CO is in the gas phase within an actively accreting disk independent of disk properties and accretion rate. The predicted optically thin water isotopolog emission is consistent with the detected H$_2^{18}$O emission toward the Stage 0 embedded young stellar objects, originating from both the disk and the warm inner envelope (hot core). An accreting embedded disk can only account for water emission arising from $R < 50$ AU, however, and the extent rapidly decreases for low accretion rates. Thus, the radial extent of the emission can be measured with ALMA observations and compared to this limit. Volatiles sublimate out to 50 AU in young disks and can reset the chemical content inherited from the envelope in periods of high accretion rates. A hot young solar nebula out to 30 AU can only have occurred during the deeply embedded Stage 0, not during the T-Tauri phase of our early solar system.

Due to instrumental limitations and a lack of disk detections, the structure between the envelope and the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologs and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A which has a confirmed Keplerian disk, is constrained from ALMA observations of DCO+ 3-2 and C18O 2-1. The physical structure of VLA1623 is obtained from the large-scale SED and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12m array and Cycle 2 ACA observations of VLA1623. DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk, suggesting a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11-16K, lower than expected from models constrained by continuum and SED. Altering the density has little effect on the DCO+ position, but increased density is needed to reproduce the disk traced in C18O. The DCO+ emission around VLA1623A is the product of shadowing of the envelope by the disk. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200AU scales, encouraging deuterated molecule production. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer.

ALMA observations of a small sample of transitional disks with large dust cavities observed in Cycle 0 and 1 are summarized. The gas and dust surface density structures are inferred from the continuum and 12CO, 13CO and C18O line data using the DALI physical-chemical code. Thanks to its ability to self-shield, CO can survive inside dust cavities in spite of being exposed to intense UV radiation and can thus be used as a probe of the gas structure. Modeling of the existing data shows that gas is present inside the dust cavities in all cases, but at a reduced level compared with the gas surface density profile of the outer disk. The gas density decrease inside the dust cavity radius by factors of up to 10^4 suggests clearing by one or more planetary-mass companions. The accompanying pressure bumps naturally lead to trapping of the mm-sized dust grains observed in the ALMA images.

Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the origin of the dust cavity. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6--5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in 12CO 3--2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the drop delta-gas of the gas density inside the dust cavity. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR21, HD135344B, LkCa15, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve constrain the dust surface density. Subsequently, the same model is used to calculate the 12CO emission, which is compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the delta-gas parameter. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter Array (SMA) data. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping the millimeter-sized dust at the edge of the cavity.

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.

Abridged: Recent simulations have explored different ways to form accretion disks around low-mass stars. We aim to present observables to differentiate a rotationally supported disk from an infalling rotating envelope toward deeply embedded young stellar objects and infer their masses and sizes. Two 3D magnetohydrodynamics (MHD) formation simulations and 2D semi-analytical model are studied. The dust temperature structure is determined through continuum radiative transfer RADMC3D modelling. A simple temperature dependent CO abundance structure is adopted and synthetic spectrally resolved submm rotational molecular lines up to $J_{\rm u} = 10$ are simulated. All models predict similar compact components in continuum if observed at the spatial resolutions of 0.5-1$"$ (70-140 AU) typical of the observations to date. A spatial resolution of $\sim$14 AU and high dynamic range ($> 1000$) are required to differentiate between RSD and pseudo-disk in the continuum. The peak-position velocity diagrams indicate that the pseudo-disk shows a flatter velocity profile with radius than an RSD. On larger-scales, the CO isotopolog single-dish line profiles are similar and are narrower than the observed line widths of low-$J$ lines, indicating significant turbulence in the large-scale envelopes. However a forming RSD can provide the observed line widths of high-$J$ lines. Thus, either RSDs are common or a higher level of turbulence ($b \sim 0.8 \ {\rm km \ s^{-1}}$ ) is required in the inner envelope compared with the outer part. Multiple spatially and spectrally resolved molecular line observations are needed. The continuum data give a better estimate on disk masses whereas the disk sizes can be estimated from the spatially resolved molecular lines observations. The general observable trends are similar between the 2D semi-analytical models and 3D MHD RSD simulations.

We present multi-epoch high-spectral resolution observations with VLT/CRIRES of the OH doublet $^2\Pi_{3/2}$ P4.5 (1+,1-) (2.934 $\mu$m) towards the protoplanetary disk around HD 100546. The OH doublet is detected at all epochs and is spectrally resolved while nearby H$_2$O lines remains undetected. The OH line velocity profile is different in the three datasets: in the first epoch (April 2012, PA=26$^{\circ}$) the OH lines are symmetric and line broadening is consistent with the gas being in Keplerian rotation around the star. No OH emission is detected within a radius of 8-11 au from the star: the line emitting region is similar in size and extent to that of the CO ro-vibrational lines. In the other two epochs (March 2013 and April 2014, PA=90$^{\circ}$ and 10$^{\circ}$, respectively) the OH lines appear asymmetric and fainter compared to April 2012. We investigate the origin of these line asymmetries which were taken by previous authors as evidence for tidal interaction between an (unseen) massive planet and the disk. We show that the observed asymmetries can be fully explained by a misalignment of the slit of order 0\farcs04-0\farcs20 with respect to the stellar position. The disk is spatially resolved and the slit misalignment is likely caused by the extended dust emission which is brighter than the stellar photosphere at near-infrared wavelengths which is the wavelength used for the pointing. This can cause the photo-center of HD 100546 to be mis-aligned with the stellar position at near-infrared wavelengths.

We investigate new and archival 3-5 $\mu$m high resolution ($\sim3$ km s$^{-1}$) spectroscopy of molecular gas in the inner disk of the young solar-mass star EX Lupi, taken during and after the strong accretion outburst of 2008. The data were obtained using the CRIRES spectrometer at the ESO Very Large Telescope in 2008 and 2014. In 2008, emission lines from CO, H$_{2}$O, and OH were detected with broad profiles tracing gas near and within the corotation radius (0.02-0.3 AU). In 2014, the spectra display marked differences. The CO lines, while still detected, are much weaker, and the H$_{2}$O and OH lines have disappeared altogether. At 3 $\mu$m a veiled stellar photospheric spectrum is observed. Our analysis finds that the molecular gas mass in the inner disk has decreased by an order of magnitude since the outburst, matching a similar decrease in the accretion rate onto the star. We discuss these findings in the context of a rapid depletion of material accumulated beyond the disk corotation radius during quiescent periods, as proposed by models of episodic accretion in EXor type young stars.

(Abridged) Organic molecules are important constituents of protoplanetary disks. Their ro-vibrational lines observed in the near- and mid-infrared are commonly detected toward T Tauri disks. These lines are the only way to probe the chemistry in the inner few au where terrestrial planets form. To understand this chemistry, accurate molecular abundances have to be determined. This is complicated by excitation effects. Most analyses so far have made the assumption of local thermal equilibrium (LTE). Starting from estimates for the collisional rate coefficients of HCN, non-LTE slab models of the HCN emission were calculated to study the importance of different excitation mechanisms. Using a new radiative transfer model, the HCN emission from a full two-dimensional disk was then modeled to study the effect of the non-LTE excitation, together with the line formation. We ran models tailored to the T Tauri disk AS 205 (N) where HCN lines in both the 3 {\mu}m and 14 {\mu}m bands have been observed by VLT-CRIRES and the Spitzer Space Telescope. Reproducing the observed 3 {\mu}m / 14 {\mu}m flux ratios requires very high densities and kinetic temperatures ($n > 10^{14}$ cm$^{-3}$ and $T > 750$ K), if only collisional excitation is accounted for. Radiative pumping can, however, excite the lines easily out to considerable radii $\sim$ 10 au. Consequently, abundances derived from LTE and non-LTE models do not differ by more than a factor of about 3. Models with both a strongly enhanced abundance within $\sim$ 1 au (jump abundance) and constant abundance can reproduce the current observations, but future observations with the MIRI instrument on JWST and METIS on the E-ELT can easily distinguish between the scenarios and test chemical models. Depending on the scenario, ALMA can detect rotational lines within vibrationally excited levels.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

One of the methods for deriving disk masses relies on direct observations of the gas, whose bulk mass is in the outer cold ($T\lesssim30$K) regions. This zone can be well traced by rotational lines of less abundant CO isotopologues, that probe the gas down to the midplane. The total CO gas mass is then obtained with the isotopologue ratios taken to be constant at the elemental isotope values found in the local ISM. This approach is however imprecise, because isotope selective processes are ignored. The aim of this work is an isotopologue selective treatment of CO isotopologues, in order to obtain a more accurate determination of disk masses. The isotope-selective photodissociation, the main process controlling the abundances of CO isotopologues in the CO-emissive layer, is properly treated for the first time in a full disk model (DALI, Bruderer et al. 2012; Bruderer 2013). The chemistry, thermal balance, line and continuum radiative transfer are all considered together with a chemical network that treats $^{13}$CO, C$^{18}$O, C$^{17}$O, isotopes of all included atoms, and molecules, as independent species. Isotope selective processes lead to regions in the disk where the isotopologues abundance ratios are considerably different from the elemental ratio. The results of this work show that considering CO isotopologue ratios as constants can lead to an underestimate of disk masses by up to almost two orders of magnitude if grains have grown to larger sizes. This may explain observed discrepancies in mass determinations from different tracers. The dependence of the various isotopologues emission on stellar and disk parameters is investigated. Including CO isotope selective processes is crucial to determine the gas mass of the disk accurately (through ALMA observations) and thus to provide the amount of gas which may eventually form planets or change the dynamics of forming planetary systems.

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.

Simple molecules like H2CO and CH3OH in protoplanetary disks are the starting point for the production of more complex organic molecules. So far, the observed chemical complexity in disks has been limited due to freeze out of molecules onto grains in the bulk of the cold outer disk. Complex molecules can be studied more directly in transitional disks with large inner holes, as these have a higher potential of detection, through UV heating of the outer disk and the directly exposed midplane at the wall. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 9 (~680 GHz) line data of the transitional disk Oph IRS 48, previously shown to have a large dust trap, to search for complex molecules in regions where planetesimals are forming. We report the detection of the H2CO 9(1,8)-8(1,7) line at 674 GHz, which is spatially resolved as a semi-ring at ~60 AU radius centered south from the star. The inferred H2CO abundance is ~10^{-8} derived by combining a physical disk model of the source with a non-LTE excitation calculation. Upper limits for CH3OH lines in the same disk give an abundance ratio H2CO/CH3OH>0.3, which points to both ice formation and gas-phase routes playing a role in the H2CO production. Upper limits on the abundances of H13CO+, CN and several other molecules in the disk are also derived and found to be consistent with full chemical models. The detection of the H2CO line demonstrates the start of complex organic molecules in a planet-forming disk. Future ALMA observations should be able to push down the abundance detection limits of other molecules by 1-2 orders of magnitude and test chemical models of organic molecules in (transitional) disks.

(Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are <~ 100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii > 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.

(Abridged) Transition disks are recognized by the absence of emission of small dust grains inside a radius of up to several 10s of AUs. Due to the lack of angular resolution and sensitivity, the gas content of such dust holes has not yet been determined, but is of importance to constrain the mechanism leading to the dust holes. Transition disks are thought to currently undergo the process of dispersal, setting an end to the giant planet formation process. We present new high-resolution observations with the Atacama Large Millimeter/ submillimeter Array (ALMA) of gas lines towards the transition disk Oph IRS 48 previously shown to host a large dust trap. ALMA has detected the $J=6-5$ line of $^{12}$CO and C$^{17}$O around 690 GHz (434 $\mu$m) at a resolution of $\sim$0.25$''$ corresponding to $\sim$30 AU (FWHM). The observed gas lines are used to set constraints on the gas surface density profile. New models of the physical-chemical structure of gas and dust in Oph IRS 48 are developed to reproduce the CO line emission together with the spectral energy distribution (SED) and the VLT-VISIR 18.7 $\mu$m dust continuum images. Integrated intensity cuts and the total spectrum from models having different trial gas surface density profiles are compared to observations. Using the derived surface density profiles, predictions for other CO isotopologues are made, which can be tested by future ALMA observations of the object. The derived gas surface density profile points to the clearing of the cavity by one or more massive planet/companion rather than just photoevaporation or grain-growth.

(Abridged) We study the response of the gas to energetic processes associated with high-mass star formation and compare it with studies on low- and intermediate-mass young stellar objects (YSOs) using the same methods. The far-IR line emission and absorption of CO, H$_2$O, OH, and [OI] reveals the excitation and the relative contribution of different species to the gas cooling budget. Herschel-PACS spectra covering 55-190 um are analyzed for ten high-mass star forming regions of various luminosities and evolutionary stages at spatial scales of ~10^4 AU. Radiative transfer models are used to determine the contribution of the envelope to the far-IR CO emission. The close environments of high-mass YSOs show strong far-IR emission from molecules, atoms, and ions. Water is detected in all 10 objects even up to high excitation lines. CO lines from J=14-13 up to typically 29-28 show a single temperature component, Trot~300 K. Typical H$_2$O temperatures are Trot~250 K, while OH has Trot~80 K. Far-IR line cooling is dominated by CO (~75 %) and to a smaller extent by OI (~20 %), which increases for the most evolved sources. H$_2$O is less important as a coolant for high-mass sources because many lines are in absorption. Emission from the envelope is responsible for ~45-85 % of the total CO luminosity in high-mass sources compared with only ~10 % for low-mass YSOs. The highest-J lines originate most likely from shocks, based on the strong correlation of CO and H$_2$O with physical parameters of the sources from low- to high-masses. Excitation of warm CO is very similar for all mass regimes, whereas H$_2$O temperatures are ~100 K higher for high-mass sources than the low-mass YSOs. Molecular cooling is ~4 times more important than cooling by [OI]. The total far-IR line luminosity is about 10$^{-3}$ and 10$^{-5}$ times lower than the dust luminosity for the low- and high-mass YSOs.

There is growing observational evidence that disk evolution is stellar-mass dependent. Here, we show that these dependencies extend to the atomic and molecular content of disk atmospheres. We analyze a unique dataset of high-resolution Spitzer/IRS spectra from 8 very low-mass star and brown dwarf disks. We report the first detections of Ne+, H2, CO2, and tentative detections of H2O toward these faint and low-mass disks. Two of our [NeII] 12.81 micron emission lines likely trace the hot (>5,000 K) disk surface irradiated by X-ray photons from the central stellar/sub-stellar object. The H2 S(2) and S(1) fluxes are consistent with arising below the fully or partially ionized surface traced by the [NeII] emission, in gas at about 600 K. We confirm the higher C2H2/HCN flux and column density ratio in brown dwarf disks previously noted from low-resolution IRS spectra. Our high-resolution spectra also show that the HCN/H2O fluxes of brown dwarf disks are on average higher than those of T Tauri disks. Our LTE modeling hints that this difference extends to column density ratios if H2O lines trace warm > 600 K disk gas. These trends suggest that the inner regions of brown dwarf disks have a lower O/C ratio than those of T Tauri disks which may result from a more efficient formation of non-migrating icy planetesimals. A O/C=1, as inferred from our analysis, would have profound implications on the bulk composition of rocky planets that can form around very low-mass stars and brown dwarfs.

Context: Rotationally supported disks are critical in the star formation process. The questions of when do they form and what factors influence or hinder their formation have been studied but are largely unanswered. Observations of early stage YSOs are needed to probe disk formation. Aims: VLA1623 is a triple non-coeval protostellar system, with a weak magnetic field perpendicular to the outflow, whose Class 0 component, VLA1623A, shows a disk-like structure in continuum with signatures of rotation in line emission. We aim to determine whether this structure is in part or in whole a rotationally supported disk, i.e. a Keplerian disk, and what are its characteristics. Methods: ALMA Cycle 0 Early Science 1.3 mm continuum and C$^{18}$O (2-1) observations in the extended configuration are presented here and used to perform an analysis of the disk-like structure using PV diagrams and thin disk modelling with the addition of foreground absorption. Results: The PV diagrams of the C$^{18}$O line emission suggest the presence of a rotationally supported component with a radius of at least 50 AU. Kinematical modelling of the line emission shows that the disk out to 180 AU is actually rotationally supported, with the rotation being well described by Keplerian rotation out to at least 150 AU, and the central source mass to be $\sim$0.2 M$_{sun}$ for an inclination of 55$^{\circ}$. Pure infall and conserved angular momentum rotation models are excluded. Conclusions: VLA1623A, a very young Class 0 source, presents a disk with an outer radius $R_{\rm out}$ = 180 AU with a Keplerian velocity structure out to at least 150 AU. The weak magnetic fields and recent fragmentation in this region of rho Ophiuchus may have played a lead role in the formation of the disk.

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.

Herschel/HIFI spectroscopic observations of CO J=10-9, CO J=16-15 and [CII] towards HD 100546 are presented. The objective is to resolve the velocity profile of the lines to address the emitting region of the transitions and directly probe the distribution of warm gas in the disk. The spectra reveal double-peaked CO line profiles centered on the systemic velocity, consistent with a disk origin. The J=16-15 line profile is broader than that of the J=10-9 line, which in turn is broader than those of lower J transitions (6-5, 3-2, observed with APEX), thus showing a clear temperature gradient of the gas with radius. A power-law flat disk model is used to fit the CO line profiles and the CO rotational ladder simultaneously, yielding a temperature of T_0=1100 \pm 350 K (at r_0 = 13 AU) and an index of q=0.85 \pm 0.1 for the temperature radial gradient. This indicates that the gas has a steeper radial temperature gradient than the dust (mean q_{dust} ~ 0.5), providing further proof of the thermal decoupling of gas and dust at the disk heights where the CO lines form. The [CII] line profile shows a strong single-peaked profile red-shifted by 0.5 km s-1 compared to the systemic velocity. We conclude that the bulk of the [CII] emission has a non-disk origin (e.g., remnant envelope or diffuse cloud).

The cosmic abundance of hydrides depends critically on high-energy UV, X-ray, and particle irradiation. Here we study hydrides in star-forming regions where irradiation by the young stellar object can be substantial, and density and temperature can be much enhanced over interstellar values. Lines of OH, CH, NH, SH and their ions OH+, CH+, NH+, SH+, H2O+, and H3O+ were observed in star-forming regions by the HIFI spectrometer onboard the Herschel Space Observatory. Molecular column densities are derived from observed ground-state lines, models, or rotational diagrams. We report here on two prototypical high-mass regions, AFGL 2591 and W3 IRS5, and compare them to chemical calculations making assumptions on the high-energy irradiation. A model assuming no ionizing protostellar emission is compared with (i) a model assuming strong protostellar X-ray emission and (ii) a two-dimensional (2D) model including emission in the far UV (FUV, 6 -- 13.6 eV) irradiating the outflow walls that separate the outflowing gas and infalling envelope material. We confirm that the effect of FUV in two dimensional models with enlarged irradiated surfaces is clearly noticeable. A molecule that is very sensitive to FUV irradiation is CH+, enhanced in abundance by more than 5 orders of magnitude. The HIFI observations of CH+ lines agree with the two-dimensional FUV model by Bruderer et al. which computes abundances, non-LTE excitation and line radiative transfer.{Ref 20} It is concluded that CH+ is a good FUV tracer in star-forming regions. The effect of potential X-ray irradiation is not excluded, but cannot be demonstrated by the present data.

CO is an important component of a protoplanetary disc as it is one of the most abundant gas phase species. Furthermore, observations of the CO transitions can be used as a diagnostic of the gas, tracing conditions in both the inner and outer disc. We present Herschel/PACS spectroscopy of a sample of 22 Herbig Ae/Be (HAEBEs) and 8 T Tauri stars (TTS), covering the pure rotational CO transitions from J = 14-13 up to J = 49-48. CO is detected in only 5 HAEBEs, AB Aur, HD 36112, HD 97048, HD 100546 and IRS 48, and in 4 TTS, AS 205, S CrA, RU Lup, and DG Tau. The highest transition seen is J = 36-35, with Eup of 3669 K, detected in HD 100546 and DG Tau. We construct rotational diagrams for the discs with at least 3 CO detections to derive Trot, and find average temperatures of 270 K for the HAEBEs and 485 K for the TTS. HD 100546 requires an extra temperature component with Trot ~ 900-1000 K, suggesting a range of temperatures in its disc atmosphere, consistent with thermo-chemical disc models. In HAEBEs, the objects with CO detections all have flared discs in which the gas and dust are thermally decoupled. We use a small model grid to analyse our observations and find that an increased amount of flaring means higher line flux, as it increases the mass in warm gas. CO is not detected in our flat discs as the emission is below the detection limit. We find that HAEBE sources with CO detections have high Luv and strong PAH emission, again connected to the heating of the gas. In TTS, the objects with CO detections are all sources with evidence for a disc wind or outflow. For both groups of objects, sources with CO detections generally have high UV luminosity (either stellar in HAEBEs or due to accretion in TTS), but this is not a sufficient condition for the detection of the far-IR CO lines.

(Abridged) Planet formation is closely related to the structure and dispersal of protoplanetary disks. A certain class of disks, called transition disks, exhibit cavities in dust images at scales of up to a few 10s of AU. The formation mechanism of the cavities is still unclear. The gas content of such cavities can be spatially resolved for the first time using the Atacama Large Millimeter/submillimeter Array (ALMA). We have developed a new series of models to simulate the physical conditions and chemical abundances of the gas in cavities to address the question whether the gas is primarily atomic or molecular inside the dust free cavities exposed to intense UV radiation. Molecular/atomic line emission by carbon monoxide (CO), its isotopologues (13CO, C18O, C17O, and 13C18O) and related species ([CI], [CII], and [OI]) is predicted for comparison with ALMA and the Herschel Space Observatory. The gas can remain in molecular form down to very low amounts of gas in the cavity. Shielding of the stellar radiation by a dusty inner disk (pre-transition disk) allows CO to survive down to lower gas masses in the cavity. The main parameter for the CO emission from cavity is the gas mass. Other parameters such as the outer disk mass, bolometric luminosity, shape of the stellar spectrum or PAH abundance are less important. Since the CO pure rotational lines readily become optically thick, the CO isotopologues need to be observed in order to quantitatively determine the amount of gas in the cavity. A wide range of gas masses in the cavity of transition disks (~4 orders of magnitude) can be probed using combined observations of CO isotopologue lines with ALMA. Measuring the gas mass in the cavity will ultimately help to distinguish between different cavity formation theories.

[abridged] We present far-infrared spectroscopic observations of PMS stars taken with Herschel/PACS as part of the DIGIT key project. The sample includes 22 Herbig AeBe and 8 T Tauri sources. Multiple atomic fine structure and molecular lines are detected at the source position: [OI], [CII], CO, OH, H_2O, CH^+. The most common feature is the [OI] 63micron line detected in almost all of the sources followed by OH. In contrast with CO, OH is detected toward both Herbig AeBe groups (flared and non-flared sources). An isothermal LTE slab model fit to the OH lines indicates column densities of 10^13 < N_OH < 10^16 cm^-2, emitting radii 15 < r < 100 AU and excitation temperatures 100 < T_ex < 400 K. The OH emission thus comes from a warm layer in the disk at intermediate stellar distances. Warm H_2O emission is detected through multiple lines toward the T Tauri systems AS 205, DG Tau, S CrA and RNO 90 and three Herbig AeBe systems HD 104237, HD 142527, HD 163296 (through line stacking). Overall, Herbig AeBe sources have higher OH/H_2O abundance ratios across the disk than do T Tauri disks, from near- to far-infrared wavelengths. Far-infrared CH^+ emission is detected toward HD 100546 and HD 97048. The slab model suggests moderate excitation (T_ex ~ 100 K) and compact (r ~ 60 AU) emission in the case of HD 100546. The [CII] emission is spatially extended in all sources where the line is detected. This suggests that not all [CII] emission is associated with the disk and that there is a substantial contribution from diffuse material around the young stars. The flux ratios of the atomic fine structure lines are consistent with a disk origin for the oxygen lines for most of the sources.

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 statistics of discovered exoplanets suggest that planets form efficiently. However, there are fundamental unsolved problems, such as excessive inward drift of particles in protoplanetary disks during planet formation. Recent theories invoke dust traps to overcome this problem. We report the detection of a dust trap in the disk around the star Oph IRS 48 using observations from the Atacama Large Millimeter/submillimeter Array (ALMA). The 0.44-millimeter-wavelength continuum map shows high-contrast crescent-shaped emission on one side of the star originating from millimeter-sized grains, whereas both the mid-infrared image (micrometer-sized dust) and the gas traced by the carbon monoxide 6-5 rotational line suggest rings centered on the star. The difference in distribution of big grains versus small grains/gas can be modeled with a vortex-shaped dust trap triggered by a companion.

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

Stars may be assembled in large growth spurts, however the evidence for this hypothesis is circumstantial. Directly studying the accretion at the earliest phases of stellar growth is challenging because young stars are deeply embedded in optically thick envelopes, which have spectral energy distributions that peak in the far-IR, where observations are difficult. In this paper, we consider the feasibility of detecting accretion outbursts from these younger stars by investigating the timescales for how the protostellar envelope responds to changes in the emission properties of the central source. The envelope heats up in response to an outburst, brightening at all wavelengths and with the emission peak moving to shorter wavelengths. The timescale for this change depends on the time for dust grains to heat and re-emit photons and the time required for the energy to escape the inner, optically-thick portion of the envelope. We find that the dust response time is much shorter than the photon propagation time and thus the timescale over which the emission varies is set by time delays imposed by geometry. These times are hours to days near the peak of the spectral energy distribution and weeks to months in the sub-mm. The ideal location to quickly detect continuum variability is therefore in the mid- to far-IR, near the peak of the spectral energy distribution, where the change in emission amplitude is largest. Searching for variability in sub-mm continuum emission is also feasible, though with a longer time separation and a weaker relationship between the amount of detected emission amplitude and change in central source luminosity. Such observations would constrain accretion histories of protostars and would help to trace the disk/envelope instabilities that lead to stellar growth.

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

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.

We report on Herschel/PACS observations of absorption lines of OH+, H2O+ and H3O+ in NGC 4418 and Arp 220. Excited lines of OH+ and H2O+ with E_lower of at least 285 and \sim200 K, respectively, are detected in both sources, indicating radiative pumping and location in the high radiation density environment of the nuclear regions. Abundance ratios OH+/H2O+ of 1-2.5 are estimated in the nuclei of both sources. The inferred OH+ column and abundance relative to H nuclei are (0.5-1)x10^{16} cm-2 and \sim2x10^{-8}, respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have OH+/H2O+\sim5-10. H3O+ is detected in both sources with N(H3O+)\sim(0.5-2)x10^{16} cm-2, and in Arp 220 the pure inversion, metastable lines indicate a high rotational temperature of ~500 K, indicative of formation pumping and/or hot gas. Simple chemical models favor an ionization sequence dominated by H+ - O+ - OH+ - H2O+ - H3O+, and we also argue that the H+ production is most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular ions arise in a relatively low density (\gtrsim10^4 cm-3) interclump medium, in which case the ionization rate per H nucleus (including secondary ionizations) is zeta>10^{-13} s-1, a lower limit that is severalx10^2 times the highest rate estimates for Galactic regions. In Arp 220, our lower limit for zeta is compatible with estimates for the cosmic ray energy density inferred previously from the supernova rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that X-ray ionization due to an AGN is responsible for the molecular ion production.

We present observations of far-infrared (50-200 micron) OH and H2O emission of the disk around the Herbig Ae star HD 163296 obtained with Herschel/PACS in the context of the DIGIT key program. In addition to strong [OI] emission, a number of OH doublets and a few weak highly excited lines of H2O are detected. The presence of warm H2O in this Herbig disk is confirmed by a line stacking analysis, enabled by the full PACS spectral scan, and by lines seen in Spitzer data. The line fluxes are analyzed using an LTE slab model including line opacity. The water column density is 10^14 - 10^15 cm^-2, and the excitation temperature is 200-300 K implying warm gas with a density n > 10^5 cm^-3. For OH we find a column density of 10^14 - 2x10^15 cm^-2 and T_ex ~ 300-500 K. For both species we find an emitting region of r ~ 15-20 AU from the star. We argue that the molecular emission arises from the protoplanetary disk rather than from an outflow. This far-infrared detection of both H2O and OH contrasts with near- and mid-infrared observations, which have generally found a lack of water in the inner disk around Herbig AeBe stars due to strong photodissociation of water. Given the similarity in column density and emitting region, OH and H2O emission seems to arise from an upper layer of the disk atmosphere of HD 163296, probing a new reservoir of water. The slightly lower temperature of H2O compared to OH suggests a vertical stratification of the molecular gas with OH located higher and water deeper in the disk, consistent with thermo-chemical models.

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

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

(Abridged) Context. With the Herschel Space Observatory, lines of simple molecules (C+, O, and CO) have been observed in the atmosphere of protoplanetary disks. When combined with ground-based [CI], all principle forms of carbon can be studied. The absence of neutral carbon [CI], which is predicted by models to be strong, can then be interpreted together with [CII] and carbon monoxide. Aims. We study the gas temperature, excitation, and chemical abundance of the simple carbon-bearing species by the method of chemical-physical modeling. We explore the sensitivity of the lines to the entering parameters and constrain the region from which the line radiation emerges. Methods. Numerical models of the radiative transfer are used together with a chemical network simulation and a calculation of the gas energetics to obtain the gas temperature. We present our new model, which is based on our previous models but includes several improvements. Results. A model of the disk around the Herbig Be star HD 100546 is able to reproduce the CO ladder together with the atomic fine-structure lines of [OI] and either [CI] or [CII]. We find that the high-J lines of CO can only be reproduced by a warm atmosphere with Tgas>>Tdust. The high-J CO observable with PACS are dominated from regions within some tens of AU. Conclusions. Only a warm atmosphere with Tgas>>Tdust can reproduce the CO ladder. The CO ladder together with [O I] and the upper limit to [CI] can be reproduced by models with a high gas/dust ratio and a low abundance of volatile carbon. These models however produce too small amounts of [CII]. Models with a low gas/dust ratio and more volatile carbon also reproduce CO and [OI], are in closer agreement with observations of [CII], but overproduce [CI]. Due to the uncertain origin of the [CII] emission, we prefer the high gas/dust ratio models, indicating a low abundance of volatile carbon.

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]

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have observed para-chloronium (H2Cl+) toward six sources in the Galaxy. We detected interstellar chloronium absorption in foreground molecular clouds along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) and W31C. Both the para-H2-35Cl+ and para-H2-37Cl+ isotopologues were detected, through observations of their 1(11)-0(00) transitions at rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho-to-para ratio of 3, the observed optical depths imply that chloronium accounts for ~ 4 - 12% of chlorine nuclei in the gas phase. We detected interstellar chloronium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar photodissociation region and the Orion South condensation. For an assumed ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply total beam-averaged column densities of ~ 2.0E+13 cm-2 and ~ 1.2E+13 cm-2, respectively, for chloronium in these two sources. We obtained upper limits on the para-H2-35Cl+ line strengths toward H2 Peak 1 in the Orion Molecular cloud and toward the massive young star AFGL 2591. The chloronium abundances inferred in this study are typically at least a factor ~10 larger than the predictions of steady-state theoretical models for the chemistry of interstellar molecules containing chlorine. Several explanations for this discrepancy were investigated, but none has proven satisfactory, and thus the large observed abundances of chloronium remain puzzling.

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.

We present a comparison of archival Spitzer spectra of the strongly variable T Tauri EX Lupi, observed before and during its 2008 outburst. We analyze the mid-infrared emission from gas-phase molecules thought to originate in a circumstellar disk. In quiescence the emission shows a forest of H2O lines, highly excited OH lines, and the Q branches of the organics C2H2, HCN, and CO2, similar to the emission observed toward several T Tauri systems. The outburst emission shows instead remarkable changes: H2O and OH line fluxes increase, new OH, H2, and HI transitions are detected, and organics are no longer seen. We adopt a simple model of a single-temperature slab of gas in local thermal equilibrium, a common approach for molecular analyses of Spitzer spectra, and derive the excitation temperature, column density, and emitting area of H2O and OH. We show how model results strongly depend on the selection of emission lines fitted, and that spectrally-resolved observations are essential for a correct interpretation of the molecular emission from disks, particularly in the case of water. Using H2O lines that can be approximated as thermalized to a single temperature, our results are consistent with a column density decrease in outburst while the emitting area of warm gas increases. A rotation diagram analysis suggests that the OH emission can be explained with two temperature components, which remarkably increase in column density in outburst. The relative change of H2O and OH emission suggests a key role for UV radiation in the disk surface chemistry.

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

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.

Molecules containing one or a few hydrogen atoms and a heavier atom (hydrides) have been predicted to trace FUV radiation. In some chemical models, FUV emission by the central object or protostar of a star forming region greatly enhances some of the hydride abundances. Two massive regions, W3 IRS5 and AFGL 2591, have been observed in hydride lines by HIFI onboard the {\it Herschel Space Observatory}. We use published results as well as new observations of CH$^+$ towards W3 IRS5. Molecular column densities are derived from ground state absorption lines, radiative transfer modeling or rotational diagrams. Models assuming no internal FUV are compared with two-dimensional models including FUV irradiation of outflow walls. We confirm that the effect of FUV is clearly noticeable and greatly improves the fit. The most sensitive molecules to FUV irradiation are CH$^+$ and OH$^+$, enhanced in abundance by many orders of magnitude. Modeling in addition also full line radiative transfer, Bruderer et al (2010b) achieve good agreement of a two-dimensional FUV model with observations of CH$^+$ in AFGL 2591. It is concluded that CH$^+$ and OH$^+$ are good FUV tracers in star-forming regions.

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.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

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.

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

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.

Context: Hydrides of the most abundant heavier elements are fundamental molecules in cosmic chemistry. Some of them trace gas irradiated by UV or X-rays. Aims: We explore the abundances of major hydrides in W3 IRS5, a prototypical region of high-mass star formation. Methods: W3 IRS5 was observed by HIFI on the Herschel Space Observatory with deep integration (about 2500 s) in 8 spectral regions. Results: The target lines including CH, NH, H3O+, and the new molecules SH+, H2O+, and OH+ are detected. The H2O+ and OH+ J=1-0 lines are found mostly in absorption, but also appear to exhibit weak emission (P-Cyg-like). Emission requires high density, thus originates most likely near the protostar. This is corroborated by the absence of line shifts relative to the young stellar object (YSO). In addition, H2O+ and OH+ also contain strong absorption components at a velocity shifted relative to W3 IRS5, which are attributed to foreground clouds. Conclusions: The molecular column densities derived from observations correlate well with the predictions of a model that assumes the main emission region is in outflow walls, heated and irradiated by protostellar UV radiation.

The Herschel Space Observatory opens the sky for observations in the far infrared at high spectral and spatial resolution. A particular class of molecules will be directly observable; light diatomic hydrides and their ions (CH, OH, SH, NH, CH+, OH+, SH+, NH+). These simple constituents are important both for the chemical evolution of the region and as tracers of high-energy radiation. If outflows of a forming star erode cavities in the envelope, protostellar far UV (FUV; 6 < E_gamma < 13.6 eV) radiation may escape through such low-density regions. Depending on the shape of the cavity, the FUV radiation then irradiates the quiescent envelope in the walls along the outflow. The chemical composition in these outflow walls is altered by photoreactions and heating via FUV photons in a manner similar to photo dominated regions (PDRs). In this work, we study the effect of cavity shapes, outflow density, and of a disk with the two-dimensional chemical model of a high-mass young stellar object introduced in the second paper in this series. We find that the shape of the cavity is particularly important in the innermost part of the envelope, where the dust temperatures are high enough (> 100 K) for water ice to evaporate. If the cavity shape allows FUV radiation to penetrate this hot-core region, the abundance of FUV destroyed species (e.g. water) is decreased. In particular, diatomic hydrides and their ions CH$+, OH+ and NH+ are enhanced by many orders of magnitude in the outflow walls due to the combination of high gas temperatures and rapid photodissociation of more saturated species. The enhancement of these diatomic hydrides is sufficient for a detection using the HIFI and PACS instruments onboard Herschel. The effect of X-ray ionization on the chemistry is found to be small, due to the much larger luminosity in FUV bands compared to X-rays.

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

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

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

Mon R2, at a distance of 830 pc, is the only ultracompact HII region (UC HII) where the photon-dominated region (PDR) between the ionized gas and the molecular cloud can be resolved with Herschel. HIFI observations of the abundant compounds 13CO, C18O, o-H2-18O, HCO+, CS, CH, and NH have been used to derive the physical and chemical conditions in the PDR, in particular the water abundance. The 13CO, C18O, o-H2-18O, HCO+ and CS observations are well described assuming that the emission is coming from a dense (n=5E6 cm-3, N(H2)>1E22 cm-2) layer of molecular gas around the UC HII. Based on our o-H2-18O observations, we estimate an o-H2O abundance of ~2E-8. This is the average ortho-water abundance in the PDR. Additional H2-18O and/or water lines are required to derive the water abundance profile. A lower density envelope (n~1E5 cm-3, N(H2)=2-5E22 cm-2) is responsible for the absorption in the NH 1_1-0_2 line. The emission of the CH ground state triplet is coming from both regions with a complex and self-absorbed profile in the main component. The radiative transfer modeling shows that the 13CO and HCO+ line profiles are consistent with an expansion of the molecular gas with a velocity law, v_e =0.5 x (r/Rout)^{-1} km/s, although the expansion velocity is poorly constrained by the observations presented here.

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

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

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

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

Observations of CO+ suggest column densities on the order 10^12 cm^-2 that can not be reproduced by many chemical models. CO+ is more likely to be destroyed than excited in collisions with hydrogen. An anomalous excitation mechanism may thus have to be considered when interpreting CO^+ observations. Chemical models are used to perform a parameter study of CO^+ abundances. Line fluxes are calculated for N(CO+)=10^12 cm^-2 and different gas densities and temperatures using a non-LTE escape probability method. The chemical formation and destruction rates are considered explicitly in the detailed balance equations of the radiative transfer. In addition, the rotational levels of CO+ are assumed to be excited upon chemical formation according to a formation temperature. It is found, that chemical models are generally able to produce high fractional CO+ abundances (x(CO+) =10^-10). In a far-ultraviolet (FUV) dominated environment, however, high abundances of CO+ are only produced in regions with a Habing field G0 > 100 and T(kin) > 600 K, posing a strong constraint on the gas temperature. For gas densities >10^6 cm^-3 and temperatures > 600 K, the combination of chemical and radiative transfer analysis shows little effect on intensities of CO+ lines with upper levels N_up <= 3. Significantly different line fluxes are calculated with an anomalous excitation mechanism, however, for transitions with higher upper levels and densities >10^6 cm ^ -3. The Herschel Space Observatory is able to reveal such effects in the terahertz wavelength regime. Ideal objects to observe are protoplanetary disks with densities 10^6 cm^-3. It is finally suggested that the CO+ chemistry may be well understood and that the abundances observed so far can be explained with a high enough gas temperature and a proper geometry.

Observations of the high-mass star forming region AFGL 2591 reveal a large abundance of CO+, a molecule known to be enhanced by far UV (FUV) and X-ray irradiation. In chemical models assuming a spherically symmetric envelope, the volume of gas irradiated by protostellar FUV radiation is very small due to the high extinction by dust. The abundance of CO+ is thus underpredicted by orders of magnitude. In a more realistic model, FUV photons can escape through an outflow region and irradiate gas at the border to the envelope. Thus, we introduce the first 2D axi-symmetric chemical model of the envelope of a high-mass star forming region to explain the CO+ observations as a prototypical FUV tracer. The model assumes an axi-symmetric power-law density structure with a cavity due to the outflow. The local FUV flux is calculated by a Monte Carlo radiative transfer code taking scattering on dust into account. A grid of precalculated chemical abundances, introduced in the first part of this series of papers, is used to quickly interpolate chemical abundances. This approach allows to calculate the temperature structure of the FUV heated outflow walls self-consistently with the chemistry. Synthetic maps of the line flux are calculated using a raytracer code. Single-dish and interferometric observations are simulated and the model results are compared to published and new JCMT and SMA observations. The two-dimensional model of AFGL 2591 is able to reproduce the JCMT single-dish observations and also explains the non-detection by the SMA. We conclude that the observed CO+ line flux and its narrow width can be interpreted by emission from the warm and dense outflow walls irradiated by protostellar FUV radiation.

Upcoming facilities such as the Herschel Space Observatory or ALMA will deliver a wealth of molecular line observations of young stellar objects (YSOs). Based on line fluxes, chemical abundances can then be estimated by radiative transfer calculations. To derive physical properties from abundances, the chemical network needs to be modeled and fitted to the observations. This modeling process is however computationally exceedingly demanding, particularly if in addition to density and temperature, far UV (FUV) irradiation, X-rays, and multi-dimensional geometry have to be considered. We develop a fast tool, suitable for various applications of chemical modeling in YSOs. A grid of the chemical composition of the gas having a density, temperature, FUV irradiation and X-ray flux is pre-calculated as a function of time. A specific interpolation approach is developed to reduce the database to a feasible size. Published models of AFGL 2591 are used to verify the accuracy of the method. A second benchmark test is carried out for FUV sensitive molecules. The novel method for chemical modeling is more than 250,000 times faster than direct modeling and agrees within a mean factor of 1.35. The tool is distributed for public use. In the course of devloping the method, the chemical evolution is explored: We find that X-ray chemistry in envelopes of YSOs can be reproduced by means of an enhanced cosmic-ray ionization rate. We further find that the abundance of CH+ in low-density gas with high ionization can be enhanced by the recombination of doubly ionized carbon (C++) and suggest a new value for the initial abundance of the main sulphur carrier in the hot-core.