The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems.

The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.

Earth's carbon deficit has been an outstanding problem in our understanding of the formation of our Solar System. A possible solution would be the sublimation of carbon grains at the so-called soot line (~300 K) early in the planet-formation process. Here, we argue that the most likely signatures of this process are an excess of hydrocarbons and nitriles inside the soot line, and a higher excitation temperature for these molecules compared to oxygen-bearing complex organics that desorb around the water snowline (~100 K). Such characteristics have been reported in the literature, for example, in Orion KL, although not uniformly, potentially due to differences in observational settings and analysis methods of different studies or related to the episodic nature of protostellar accretion. If this process is active, this would mean that there is a heretofore unknown component to the carbon chemistry during the protostellar phase that is acting from the top down - starting from the destruction of larger species - instead of from the bottom up from atoms. In the presence of such a top-down component, the origin of organic molecules needs to be re-explored.

The recent ALMA DSHARP survey provided illuminating results on the diversity of substructures in planet forming disks. These substructures trace pebble-sized grains accumulated at local pressure maxima, possibly due to planet-disk interactions or other planet formation processes. DSHARP sources are heavily biased to large and massive disks that only represent the high (dust flux) tail end of the disk population. Thus it is unclear whether similar substructures and corresponding physical processes also occur in the majority of disks which are fainter and more compact. Here we explore the presence and characteristics of features in a compact disk around GQ Lup A, the effective radius of which is 1.5 to 10 times smaller than those of DSHARP disks. We present our analysis of ALMA 1.3mm continuum observations of the GQ Lup system. By fitting visibility profiles of the continuum emission, we find substructures including a gap at ~ 10 au. The compact disk around GQ Lup exhibits similar substructures to those in the DSHARP sample, suggesting that mechanisms of trapping pebble-sized grains are at work in small disks as well. Characteristics of the feature at ~ 10 au, if due to a hidden planet, are evidence of planet formation at Saturnian distances. Our results hint at a rich world of substructures to be identified within the common population of compact disks, and subsequently a population of solar system analogs within these disks. Such study is critical to understanding the formation mechanisms and planet populations in the majority of protoplanetary disks.

Comets spend most of their lives at large distances from any star, during which time their interior compositions remain relatively unaltered. Cometary observations can therefore provide direct insight into the chemistry that occurred during their birth at the time of planet formation. To-date, there have been no confirmed observations of parent volatiles (gases released directly from the nucleus) of a comet from any planetary system other than our own. Here we present high-resolution, interferometric observations of 2I/Borisov, the first confirmed interstellar comet, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) on 15th-16th December 2019. Our observations reveal emission from hydrogen cyanide (HCN), and carbon monoxide (CO), coincident with the expected position of 2I/Borisov's nucleus, with production rates Q(HCN)=$(7.0\pm1.1)\times10^{23}$ s$^{-1}$ and Q(CO)=$(4.4\pm0.7)\times10^{26}$ s$^{-1}$. While the HCN abundance relative to water (0.06-0.16%) appears similar to that of typical, previously observed comets in our Solar System, the abundance of CO (35-105%) is among the highest observed in any comet within 2 au of the Sun. This shows that 2I/Borisov must have formed in a relatively CO-rich environment - probably beyond the CO ice-line in the very cold, outer regions of a distant protoplanetary accretion disk, as part of a population of small, icy bodies analogous to our Solar System's own proto-Kuiper Belt.

The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds.

Recent observations show that the CO gas abundance, relative to H$_2$, in many 1-10 Myr old protoplanetary disks may be heavily depleted, by a factor of 10-100 compared to the canonical interstellar medium value of 10$^{-4}$. When and how this depletion happens can significantly affect compositions of planetesimals and atmospheres of giant planets. It is therefore important to constrain if the depletion occurs already at the earliest protostellar disk stage. Here we present spatially resolved observations of C$^{18}$O, C$^{17}$O, and $^{13}$C$^{18}$O $J$=2-1 lines in three protostellar disks. We show that the C$^{18}$O line emits from both the disk and the inner envelope, while C$^{17}$O and $^{13}$C$^{18}$O lines are consistent with a disk origin. The line ratios indicate that both C$^{18}$O and C$^{17}$O lines are optically thick in the disk region, and only $^{13}$C$^{18}$O line is optically thin. The line profiles of the $^{13}$C$^{18}$O emissions are best reproduced by Keplerian gaseous disks at similar sizes as their mm-continuum emissions, suggesting small radial separations between the gas and mm-sized grains in these disks, in contrast to the large separation commonly seen in protoplanetary disks. Assuming a gas-to-dust ratio of 100, we find that the CO gas abundances in these protostellar disks are consistent with the ISM abundance within a factor of 2, nearly one order of magnitude higher than the average value of 1-10 Myr old disks. These results suggest that there is a fast, $\sim$1 Myr, evolution of the abundance of CO gas from the protostellar disk stage to the protoplanetary disk stage.

The atmospheric composition of giant planets carries the information of their formation history. Superstellar C/H ratios are seen in atmospheres of Jupiter, Saturn, and various giant exoplanets. Also, giant exoplanets show a wide range of C/O ratio. To explain these ratios, one hypothesis is that protoplanets accrete carbon-enriched gas when a large number of icy pebbles drift across the CO snowline. Here we report the first direct evidence of an elevated C/H ratio in disk gas. We use two thermo-chemical codes to model the $^{13}$C$^{18}$O, C$^{17}$O, and C$^{18}$O (2-1) line spectra of the HD 163296 disk. We show that the gas inside the CO snowline ($\sim$70 au) has a C/H ratio of 1-2 times higher than the stellar value. This ratio exceeds the expected value substantially, as only 25-60% of the carbon should be in gas at these radii. Although we cannot rule out the case of a normal C/H ratio inside 70 au, the most probable solution is an elevated C/H ratio of 2-8 times higher than the expectation. Our model also shows that the gas outside 70 au has a C/H ratio of 0.1$\times$ the stellar value. This picture of enriched C/H gas at the inner region and depleted gas at the outer region is consistent with numerical simulations of icy pebble growth and drift in protoplanetary disks. Our results demonstrate that the large-scale drift of icy pebble can occur in disks and may significantly change the disk gas composition for planet formation.

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.

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

Protoplanetary disks are known to posses a stunning variety of substructure in the distribution of their mm~sized grains, predominantly seen as rings and gaps (Andrews et al. 2018), which are frequently interpreted as due to the shepherding of large grains by either hidden, still-forming planets within the disk (Zhang et al. 2018) or (magneto-)hydrodynamic instabilities (Flock et al. 2015). The velocity structure of the gas offers a unique probe of both the underlying mechanisms driving the evolution of the disk, the presence of embedded planets and characterising the transportation of material within the disk, such as following planet-building material from volatile-rich regions to the chemically-inert midplane, or detailing the required removal of angular momentum. Here we present the radial profiles of the three velocity components of gas in upper disk layers in the disk of HD 163296 as traced by 12CO molecular emission. These velocities reveal significant flows from the disk surface towards the midplane of disk at the radial locations of gaps argued to be opened by embedded planets (Isella et al. 2016, 2018, Teague et al. 2018, Pinte et al. 2018), bearing striking resemblance to meridional flows, long predicted to occur during the early stages of planet formation (Szulagyi et al. 2014, Morbidelli et al. 2014, Fung & Chiang 2016, Dong et al. 2019). In addition, a persistent radial outflow is seen at the outer edge of the disk, potentially the base of a wind associated with previously detected extended emission (Klaassen et al. 2013).

We report the detection of spiral substructure in both the gas velocity and temperature structure of the disk around TW~Hya, suggestive of planet-disk interactions with an unseen planet. Perturbations from Keplerian rotation tracing out a spiral pattern are observed in the SE of the disk, while significant azimuthal perturbations in the gas temperature are seen in the outer disk, outside 90~au, extending the full azimuth of the disk. The deviation in velocity is either $\Delta v_{\phi} \, / \, v_{\rm kep} \sim 0.1$ or $\Delta v_{z} \, / \, v_{\rm kep} \sim 0.01$ depending on whether the perturbation is in the rotational or vertical direction, while radial perturbations can be ruled out. Deviations in the gas temperature are $\pm 4$ K about the azimuthally averaged profile, equivalent to deviations of $\Delta T_{\rm gas} \, / \, T_{\rm gas} \sim 0.05$. Assuming all three structures can be described by an Archimedean spiral, measurements of the pitch angles of both velocity and temperature spirals show a radially decreasing trend for all three, ranging from 9$^{\circ}$ at 70 au, dropping to 3$^{\circ}$ at 200 au. Such low pitch-angled spirals are not readily explained through the wake of an embedded planet in the location of previously reported at 94 au, but rather require a launching mechanism which results in much more tightly wound spirals. Molecular emission tracing distinct heights in the disk is required to accurately distinguish between spiral launching mechanisms.

Phosphorus is a key ingredient in terrestrial biochemistry, but is rarely observed in the molecular ISM and therefore little is known about how it is inherited during the star and planet formation sequence. We present observations of the phosphorus-bearing molecules PO and PN towards the Class I low-mass protostar B1-a using the IRAM 30m telescope, representing the second detection of phosphorus carriers in a Solar-type star forming region. The P/H abundance contained in PO and PN is ~10$^{-10}$-10$^{-9}$ depending on the assumed source size, accounting for just 0.05-0.5% of the solar phosphorus abundance and implying significant sequestration of phosphorus in refractory material. Based on a comparison of the PO and PN line profiles with the shock tracers SiO, SO$_2$, and CH$_3$OH, the phosphorus molecule emission seems to originate from shocked gas and is likely associated with a protostellar outflow. We find a PO/PN column density ratio of ~1-3, which is consistent with the values measured in the shocked outflow of the low-mass protostar L1157, the massive star-forming regions W51 and W3(OH), and the galactic center GMC G+0.693-0.03. This narrow range of PO/PN ratios across sources with a range of environmental conditions is surprising, and likely encodes information on how phosphorus carriers are stored in grain mantles.

Sulfur is one of the most abundant elements in the Universe, with important roles in astro-, geo-, and biochemistry. Its main reservoirs in planet-forming disks have previously eluded detection: gaseous molecules only account for $<1\,$\% of total elemental sulfur, with the rest likely in either ices or refractory minerals. Mechanisms such as giant planets can filter out dust from gas accreting onto disk-hosting stars. For stars above 1.4 solar masses, this leaves a chemical signature on the stellar photosphere that can be used to determine the fraction of each element that is locked in dust. Here, we present an application of this method to sulfur, zinc, and sodium. We analyse the accretion-contaminated photospheres of a sample of young stars and find $(89\pm8)\,$\% of elemental sulfur is in refractory form in their disks. The main carrier is much more refractory than water ice, consistent with sulfide minerals such as FeS.

CO is the most widely used gas tracer of protoplanetary disks. Its abundance is usually assumed to be an interstellar ratio throughout the warm molecular layer of the disk. But recent observations of low CO gas abundance in many protoplanetary disks challenge our understanding of physical and chemical evolutions in disks. Here we investigate the CO abundance structures in four well-studied disks and compare their structures with predictions of chemical processing of CO and transport of CO ice-coated dust grains in disks. We use spatially resolved CO isotopologue line observations and detailed thermo-chemical models to derive CO abundance structures. We find that the CO abundance varies with radius by an order of magnitude in these disks. We show that although chemical processes can efficiently reduce the total column of CO gas within 1 Myr under an ISM level of cosmic-ray ionization rate, the depletion mostly occurs at the deep region of a disk. Without sufficient vertical mixing, the surface layer is not depleted enough to reproduce weak CO emissions observed. The radial profiles of CO depletion in three disks are qualitatively consistent with predictions of pebble formation, settling, and drifting in disks. But the dust evolution alone cannot fully explain the high depletion observed in some disks. These results suggest that dust evolution may play a significant role in transporting volatile materials and a coupled chemical-dynamical study is necessary to understand what raw materials are available for planet formation at different distances from the central star.

Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of both molecules, because it is destroyed by CO and formed from N$_2$. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of N$_2$H$^+$ at 0.2$''$-0.4$''$ resolution in the disks around LkCa 15, GM Aur, DM Tau, V4046 Sgr, AS 209, and IM Lup. We find two distinctive emission morphologies: N$_2$H$^+$ is either present in a bright, narrow ring surrounded by extended tenuous emission, or in a broad ring. These emission patterns can be explained by two different kinds of vertical temperature structures. Bright, narrow N$_2$H$^+$ rings are expected in disks with a thick Vertically Isothermal Region above the Midplane (VIRaM) layer (LkCa 15, GM Aur, DM Tau) where the N$_2$H$^+$ emission peaks between the CO and N$_2$ snowlines. Broad N$_2$H$^+$ rings come from disks with a thin VIRaM layer (V4046 Sgr, AS 209, IM Lup). We use a simple model to extract the first sets of CO and N$_2$ snowline pairs and corresponding freeze-out temperatures towards the disks with a thick VIRaM layer. The results reveal a range of N$_2$ and CO snowline radii towards stars of similar spectral type, demonstrating the need for empirically determined snowlines in disks.

We emphasize that the completeness of the partition function, that is, the use of a converged partition function at the typical temperature range of the survey, is very important to decrease the uncertainty on this quantity and thus to derive reliable interstellar molecular densities. In that context, we show how the use of different approximations for the rovibrational partition function together with some interpolation and/or extrapolation procedures may affect the estimate of the interstellar molecular column density. For that purpose, we apply the partition function calculations to astronomical observations performed with the IRAM-30m telescope towards the NGC7538-IRS1 source of two N-bearing molecules: isocyanic acid (HNCO, a quasilinear molecule) and methyl cyanide (CH$_3$CN, a symmetric top molecule). The case of methyl formate (HCOOCH$_3$), which is an asymmetric top O-bearing molecule containing an internal rotor is also discussed. Our analysis shows that the use of different partition function approximations leads to relative differences in the resulting column densities in the range 9 to 43\%. Thus, we expect this work to be relevant for surveys of sources with temperatures higher than 300~K and to observations in the infrared.

We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations.

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

CO is thought to be the main reservoir of volatile carbon in protoplanetary disks, and thus the primary initial source of carbon in the atmospheres of forming giant planets. However, recent observations of protoplanetary disks point towards low volatile carbon abundances in many systems, including at radii interior to the CO snowline. One potential explanation is that gas phase carbon is chemically reprocessed into less volatile species, which are frozen on dust grain surfaces as ice. This mechanism has the potential to change the primordial C/H ratio in the gas. However, current observations primarily probe the upper layers of the disk. It is not clear if the low volatile carbon abundances extend to the midplane, where planets form. We have run a grid of 198 chemical models, exploring how the chemical reprocessing of CO depends on disk mass, dust grain size distribution, temperature, cosmic ray and X-ray ionization rate, and initial water abundance. Building on our previous work focusing on the warm molecular layer, here we analyze the results for our grid of models in the disk midplane at 12 au. We find that either an ISM level cosmic ray ionization rate or the presence of UV photons due to a low dust surface density are needed to chemically reduce the midplane CO gas abundance by at least an order of magnitude within 1 Myr. In the majority of our models CO does not undergo substantial reprocessing by in situ chemistry and there is little change in the gas phase C/H and C/O ratios over the lifetime of the typical disk. However, in the small sub-set of disks where the disk midplane is subject to a source of ionization or photolysis, the gas phase C/O ratio increases by up to nearly 9 orders of magnitude due to conversion of CO into volatile hydrocarbons.

Molecular lines observed towards protoplanetary disks carry information about physical and chemical processes associated with planet formation. We present ALMA Band 6 observations of C2H, HCN, and C18O in a sample of 14 disks spanning a range of ages, stellar luminosities, and stellar masses. Using C2H and HCN hyperfine structure fitting and HCN/H13CN isotopologue analysis, we extract optical depth, excitation temperature, and column density radial profiles for a subset of disks. C2H is marginally optically thick (\tau ~1-5) and HCN is quite optically thick (\tau ~ 5-10) in the inner 200 AU. The extracted temperatures of both molecules are low (10-30K), indicative of either sub-thermal emission from the warm disk atmosphere or substantial beam dilution due to chemical substructure. We explore the origins of C2H morphological diversity in our sample using a series of toy disk models, and find that disk-dependent overlap between regions with high UV fluxes and high atomic carbon abundances can explain a wide range of C2H emission features (e.g. compact vs. extended and ringed vs. ringless emission). We explore the chemical relationship between C2H, HCN, and C18O and find a positive correlation between C2H and HCN fluxes, but no relationship between C2H or HCN with C18O fluxes. We also see no evidence that C2H and HCN are enhanced with disk age. C2H and HCN seem to share a common driver, however more work remains to elucidate the chemical relationship between these molecules and the underlying evolution of C, N, and O chemistries in disks.

Line ratios for different transitions of the same molecule have long been used as a probe of gas temperature. Here we use ALMA observations of the N2H+ J~=~1-0 and J~=~4-3 lines in the protoplanetary disk around TW Hya to derive the temperature at which these lines emit. We find an averaged temperature of 39~K with a one sigma uncertainty of 2~K for the radial range 0.8-2'', significantly warmer than the expected midplane temperature beyond 0.5'' in this disk. We conclude that the N2H+ emission in TW Hya is not emitting from near the midplane, but rather from higher in the disk, in a region likely bounded by processes such as photodissociation or chemical reprocessing of CO and N2 rather than freeze out.

The chemistry of complex organic molecules in interstellar dark clouds is still highly uncertain in part because of the lack of constraining observations. Orion is the closest massive star-forming region, and observations making use of ALMA allow us to separate the emission regions of various complex organic molecules (COMs) in both velocity and space. Orion also benefits from an exceptional situation, in that it is the site of a powerful explosive event that occurred 550 years ago. We show that the closely surrounding Kleinmann-Low region has clearly been influenced by this explosion; some molecular species have been pushed away from the densest parts while others have remained in close proximity. This dynamical segregation reveals the time dependence of the chemistry and, therefore allows us to better constrain the formation sequence of COMs and other species, including deuterated molecules.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns.

Water is fundamental to our understanding of the evolution of planetary systems and the delivery of volatiles to the surfaces of potentially habitable planets. Yet, we currently have essentially no facilities capable of observing this key species comprehensively. With this white paper, we argue that we need a relatively large, cold space-based observatory equipped with a high-resolution spectrometer, in the mid- through far-infrared wavelength range (25-600~$\mu$m) in order to answer basic questions about planet formation, such as where the Earth got its water, how giant planets and planetesimals grow, and whether water is generally available to planets forming in the habitable zone of their host stars.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

One of the primary mechanisms for inferring the dynamical history of planets in our Solar System and in exoplanetary systems is through observation of elemental ratios (i.e. C/O). The ability to effectively use these observations relies critically on a robust understanding of the chemistry and evolutionary history of the observed abundances. Significant efforts have been devoted to this area from within astrochemistry circles, and these efforts should be supported going forward by the larger exoplanetary science community. In addition, the construction of a next-generation radio interferometer will be required to test many of these predictive models in situ, while simultaneously providing the resolution necessary to pinpoint the location of planets in formation.

The water snowline in protoplanetary disks is one of the most pivotal locations for planet formation: it sets a critical boundary of chemical composition in the disk and likely serves as a favorable site of planetesimal growth and planet formation. However, the water snowline at the mid-plane cannot be directly traced by water line observations as the dust emission around the water snowline is highly optically thick at these line frequencies. Alternatively, the water snowline maybe traced through the sharp changes in physical properties (surface density and dust size distribution) at this transition or through the snowlines of other molecules that sublimate co-spatially with water. In this chapter, we discuss the ngVLA's capability of tracing the mid-plane water snowline through mapping sharp transitions in the continuum spectral index caused by changes of dust size distributions and by using the NH3 snowline as a proxy. We found that with the current design of the ngVLA ability, a sharp transition in spectral index across the water snowline can be readily traced in relatively bright disks. A detection of the NH3 snowline is challenging under the current design, but may still be possible if the NH3-to-water abundance in the disk is close to that of the ISM ices. In summary, the proposed ngVLA is the only facility can provide sufficient sensitivity and spatial resolution to trace the evolution of water snowline at the disk mid-plane and to reveal its role in accelerating the planetesimal formation and in building the architecture of planetary systems.

We present an improved method to measure the rotation curves for disks with non-axisymmetric brightness profiles initially published in Teague et al. (2018a). Application of this method to the well studied AS$~$209 system shows substantial deviations from Keplerian rotation of up to $\pm 5\%$. These deviations are most likely due to perturbations in the gas pressure profile, including a perturbation located at $\approx 250~$au and spanning up to $\approx 50~$au which is only detected kinematically. Modelling the required temperature and density profiles required to recover the observed rotation curve we demonstrate that the rings observed in $\mu$m scattered light are coincident with the pressure maxima, and are radially offset from the rings observed in mm continuum emission. This suggests that if rings in the NIR are due to sub-$\mu$m grains trapped in pressure maxima that there is a vertical dependence on the radius of the pressure minima.

It is not known whether the original carriers of Earth's nitrogen were molecular ices or refractory dust. To investigate this question, we have used data and results of Herschel observations towards two protostellar sources: the high-mass hot core of Orion KL, and the low-mass protostar IRAS 16293-2422. Towards Orion KL, our analysis of the molecular inventory of Crockett et al. (2014) indicates that HCN is the organic molecule that contains by far the most nitrogen, carrying $74_{-9}^{+5}\%$ of nitrogen-in-organics. Following this evidence, we explore HCN towards IRAS 16293-2422, which we consider a solar analog. Towards IRAS 16293-2422, we have reduced and analyzed Herschel spectra of HCN, and fit these observations against "jump" abundance models of IRAS 16293-2422's protostellar envelope. We find an inner-envelope HCN abundance $X_{\textrm{in}} = 5.9\pm0.7 \times 10^{-8}$ and an outer-envelope HCN abundance $X_{\textrm{out}} = 1.3 \pm 0.1 \times 10^{-9}$. We also find the sublimation temperature of HCN to be $T_{\textrm{jump}} = 71 \pm 3$~K; this measured $T_{\textrm{jump}}$ enables us to predict an HCN binding energy $E_{\textrm{B}}/k = 3840 \pm 140$~K. Based on a comparison of the HCN/H2O ratio in these protostars to N/H2O ratios in comets, we find that HCN (and, by extension, other organics) in these protostars is incapable of providing the total bulk N/H2O in comets. We suggest that refractory dust, not molecular ices, was the bulk provider of nitrogen to comets. However, interstellar dust is not known to have 15N enrichment, while high 15N enrichment is seen in both nitrogen-bearing ices and in cometary nitrogen. This may indicate that these 15N-enriched ices were an important contributor to the nitrogen in planetesimals and likely to the Earth.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

In young circumstellar disks, accretion--the inspiral of disk material onto the central star--is important for both the buildup of stellar masses and the outcome of planet formation. Although the existence of accretion is well documented, understanding the angular momentum transport mechanism that enables disk accretion has proven to be an enduring challenge. The leading theory to date, the magnetorotational instability, which redistributes angular momentum within the disk, is increasingly questioned, and magnetothermal disk winds, which remove angular momentum from the disk, have emerged as an alternative theoretical solution. Here we investigate whether measurements of disk radii can provide useful insights into which, if either, of these mechanisms drive disk accretion, by searching for evidence of viscous spreading in gaseous disks, a potential signature of "in disk" angular momentum transport. We find that the large sizes of most Class II (T Tauri) gas disks compared to those of their earlier evolutionary counterparts, Class I gas disks, are consistent with expectations for viscous spreading in the Class II phase. There is, however, a large spread in the sizes of Class II gas disks at any age, including a population of very small Class II gas disks. Their small sizes may result from processes such as photoevaporation, disk winds, or truncation by orbiting low mass companions.

Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, $\zeta_{\rm H}$, should correlate with the star-formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star forming galaxies in the Universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH$^+$ and H$_2$O$^+$ are useful for inferring $\zeta_{\rm H}$ when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts $z\gtrsim2$ transitions out of the ground rotational states of OH$^+$ and H$_2$O$^+$ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135$-$0102 and SDP 17b. These detections enable an exploration of $\zeta_{\rm H}$ in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH$^+$ and H$_2$O$^+$ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of $\zeta_{\rm H}\sim10^{-16}$-$10^{-14}$ s$^{-1}$ in these extended gaseous halos. As our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher.

Current models of (exo)planet formation often rely on a large influx of so-called `pebbles' from the outer disk into the planet formation region. In this paper, we investigate how the formation of pebbles in the cold outer regions of protoplanetary disks and their subsequent migration to the inner disk can alter the gas-phase CO distribution both interior and exterior to the midplane CO snowline. By simulating the resulting CO abundances in the midplane as well as the warm surface layer, we identify observable signatures of large-scale pebble formation and migration that can be used as `smoking guns' for these important processes. Specifically, we find that after $1\mathrm{~Myr}$, the formation and settling of icy pebbles results in the removal of up to $80\%$ of the CO vapor in the warm ($T>22\mathrm{~K}$) disk layers outside the CO snowline, while the radial migration of pebbles results in the generation of a plume of CO vapor interior the snowline, increasing the CO abundance by a factor ${\sim}2{-}6$ depending on the strength of the turbulence and the sizes of the individual pebbles. The absence of this plume of CO vapor in young nearby disks could indicate efficient conversion of CO into a more refractory species, or a reduction in the radial mass flux of pebbles by, for example, disk inhomogeneities or early planetesimal formation.

Observations of multiple rotational transitions from a single molecule allow for unparalleled constraints on the physical conditions of the emitting region. We present an analysis of CS in TW~Hya using the $J=7-6$, $5-4$ and $3-2$ transitions imaged at $\sim 0.5^{\prime\prime}$ spatial resolution, resulting in a temperature and column density profile of the CS emission region extending out to 230~au, far beyond previous measurements. In addition, the 15~kHz resolution of the observations and the ability to directly estimate the temperature of the CS emitting gas, allow for one of the most sensitive searches for turbulent broadening in a disk to date. Limits of $v_{\rm turb} \lesssim 0.1 c_s$ can be placed across the entire radius of the disk. We are able to place strict limits of the local H$_2$ density due to the collisional excitations of the observed transitions. From these we find that a minimum disk mass of $3 \times 10^{-4}~M_{\rm sun}$ is required to be consistent with the CS excitation conditions and can uniquely constrain the gas surface density profile in the outer disk.

In this article we review the methods used to determine the gas and dust masses of protoplanetary disks, with an emphasis on the lesser characterized total gas mass. Our review encompasses all the indirect tracers and the methodology that is be used to isolate the hidden H2 via dust, CO, and HD emission. We discuss the overall calibration of gaseous tracers which is based on decades of study of the dense phases of the interstellar medium. At present, disk gas masses determined via CO and HD are (in a few instances) different by orders of magnitude, hinting at either significant evolution in total disk mass or in the CO abundance . Either of these would represent a fundamental physical or chemical process that appears to dominate the system on ~million year timescales. Efforts to reconcile these differences using existing and future facilities are discussed.

In this chapter we outline some of the basic understanding of the chemistry that accompanies planet formation. We discuss the basic physical environment which dictates the dominant chemical kinetic pathways for molecule formation. We focus on three zones from both observational and theoretical perspectives: (1) the planet forming midplane and ice/vapor transition zones (snow-lines), (2) the warm disk surface that is shielded from radiation, which can be readily accessed by todays observational facilities, and (3) the surface photodissociation layers where stellar radiation dominates. We end with a discussion of how chemistry influences planet formation along with how to probe the link between formation and ultimate atmospheric composition for gas giants and terrestrial worlds.

The formation of asteroids, comets and planets occurs in the interior of protoplanetary disks during the early phase of star formation. Consequently, the chemical composition of the disk might shape the properties of the emerging planetary system. In this context, it is crucial to understand whether and what organic molecules are synthesized in the disk. In this Letter, we report the first detection of formic acid (HCOOH) towards the TW Hydrae protoplanetary disk. The observations of the trans-HCOOH 6$_{(1,6)-5(1,5)}$ transition were carried out at 129~GHz with ALMA. We measured a disk-averaged gas-phase t-HCOOH column density of $\sim$ (2-4)$\times$10$^{12}$~cm$^{-2}$, namely as large as that of methanol. HCOOH is the first organic molecules containing two oxygen atoms detected in a protoplanetary disk, a proof that organic chemistry is very active even though difficult to observe in these objects. Specifically, this simplest acid stands as the basis for synthesis of more complex carboxylic acids used by life on Earth.

We present the first kinematical detection of embedded protoplanets within a protoplanetary disk. Using archival ALMA observations of HD 163296, we demonstrate a new technique to measure the rotation curves of CO isotopologue emission to sub-percent precision relative to the Keplerian rotation. These rotation curves betray substantial deviations caused by local perturbations in the radial pressure gradient, likely driven by gaps carved in the gas surface density by Jupiter-mass planets. Comparison with hydrodynamic simulations shows excellent agreement with the gas rotation profile when the disk surface density is perturbed by two Jupiter mass planets at 83 au and 137 au. As the rotation of the gas is dependent on the pressure of the total gas component, this method provides a unique probe of the gas surface density profile without incurring significant uncertainties due to gas-to-dust ratios or local chemical abundances which plague other methods. Future analyses combining both methods promise to provide the most accurate and robust measures of embedded planetary mass. Furthermore, this method provides a unique opportunity to explore wide-separation planets beyond the mm continuum edge and to trace the gas pressure profile essential in modelling grain evolution in disks.

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

In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

At a time when ALMA produces spectacular high resolution images of gas and dust in circumstellar disks, the next observational frontier in our understanding of planet formation and the chemistry of planet-forming material may be found in the mid- to far-infrared wavelength range. A large, actively cooled far-infrared telescope in space will offer enormous spectroscopic sensitivity improvements of 3-4 orders of magnitude, making it possible to uniquely survey certain fundamental properties of planet formation. Specifically, the Origins Space Telescope (OST), a NASA flagship concept to be submitted to the 2020 decadal survey, will provide a platform that allows complete surveys of warm and cold water around young stars of all masses and across all evolutionary stages, and to measure their total planet-forming gas mass using the ground-state line of HD. While this white paper is formulated in the context of the NASA Origins Space Telescope concept, it can be applied in general to inform any future space-based, cold far-infrared observatory.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

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

CO is commonly used as a tracer of the total gas mass in both the interstellar medium and in protoplanetary disks. Recently there has been much debate about the utility of CO as a mass tracer in disks. Observations of CO in protoplanetary disks reveal a range of CO abundances, with measurements of low CO to dust mass ratios in numerous systems. One possibility is that carbon is removed from CO via chemistry. However, the full range of physical conditions conducive to this chemical reprocessing is not well understood. We perform a systematic survey of the time dependent chemistry in protoplanetary disks for 198 models with a range of physical conditions. We varying dust grain size distribution, temperature, comic ray and X-ray ionization rate, disk mass, and initial water abundance, detailing what physical conditions are necessary to activate the various CO depletion mechanisms in the warm molecular layer. We focus our analysis on the warm molecular layer in two regions: the outer disk (100 au) well outside the CO snowline and the inner disk (19 au) just inside the midplane CO snow line. After 1 Myr, we find that the majority of models have a CO abundance relative to H$_2$ less than $10^{-4}$ in the outer disk, while an abundance less than $10^{-5}$ requires the presence of cosmic rays. Inside the CO snow line, significant depletion of CO only occurs in models with a high cosmic ray rate. If cosmic rays are not present in young disks it is difficult to chemically remove carbon from CO. Additionally, removing water prior to CO depletion impedes the chemical processing of CO. Chemical processing alone cannot explain current observations of low CO abundances. Other mechanisms must also be involved.

Snowlines are key ingredients for planet formation. Providing observational constraints on the locations of the major snowlines is therefore crucial for fully connecting planet compositions to their formation mechanism. Unfortunately, the most important snowline, that of water, is very difficult to observe directly in protoplanetary disks due to its close proximity to the central star. Based on chemical considerations, HCO$^+$ is predicted to be a good chemical tracer of the water snowline, because it is particularly abundant in dense clouds when water is frozen out. This work maps the optically thin isotopologue H$^{13}$CO$^+$ ($J=3-2$) toward the envelope of the low-mass protostar NGC1333-IRAS2A (observed with NOEMA at ~0.9" resolution), where the snowline is at larger distance from the star than in disks. The H$^{13}$CO$^+$ emission peaks ~2" northeast of the continuum peak, whereas the previously observed H$_2^{18}$O shows compact emission on source. Quantitative modeling shows that a decrease in H$^{13}$CO$^+$ abundance by at least a factor of six is needed in the inner ~360 AU to reproduce the observed emission profile. Chemical modeling predicts indeed a steep increase in HCO$^+$ just outside the water snowline; the 50% decrease in gaseous H$_2$O at the snowline is not enough to allow HCO$^+$ to be abundant. This places the water snowline at 225 AU, further away from the star than expected based on the 1D envelope temperature structure for NGC1333-IRAS2A. In contrast, DCO$^+$ observations show that the CO snowline is at the expected location, making an outburst scenario unlikely. The spatial anticorrelation of the H$^{13}$CO$^+$ and H$_2^{18}$O emission provide a proof of concept that H$^{13}$CO$^+$ can be used as a tracer of the water snowline.

Starburst galaxies at the peak of cosmic star formation are among the most extreme starforming engines in the universe, producing stars over ~100 Myr. The star formation rates of these galaxies, which exceed 100 $M_\odot$ per year, require large reservoirs of cold molecular gas to be delivered to their cores, despite strong feedback from stars or active galactic nuclei. Starburst galaxies are therefore ideal targets to unravel the critical interplay between this feedback and the growth of a galaxy. The methylidyne cation, CH$^+$, is a most useful molecule for such studies because it cannot form in cold gas without supra-thermal energy input, so its presence highlights dissipation of mechanical energy or strong UV irradiation. Here, we report the detection of CH$^+$(J=1-0) emission and absorption lines in the spectra of six lensed starburst galaxies at redshifts z~2.5. This line has such a high critical density for excitation that it is emitted only in very dense ($>10^5$ cm$^{-3}$) gas, and is absorbed in low-density gas. We find that the CH$^+$ emission lines, which are broader than 1000 km s$^{-1}$, originate in dense shock waves powered by hot galactic winds. The CH$^+$ absorption lines reveal highly turbulent reservoirs of cool ($T\sim 100$K), low-density gas, extending far outside (>10 kpc) the starburst cores (radii <1 kpc). We show that the galactic winds sustain turbulence in the 10 kpc-scale environments of the starburst cores, processing these environments into multi-phase, gravitationally bound reservoirs. However, the mass outflow rates are found to be insufficient to balance the star formation rates. Another mass input is therefore required for these reservoirs, which could be provided by on-going mergers or cold stream accretion. Our results suggest that galactic feedback, coupled jointly to turbulence and gravity, extends the starburst phase instead of quenching it.

The Earth and other rocky bodies in the inner solar system contain significantly less carbon than the primordial materials that seeded their formation. These carbon-poor objects include the parent bodies of primitive meteorites, suggesting that at least one process responsible for solid-phase carbon depletion was active prior to the early stages of planet formation. Potential mechanisms include the erosion of carbonaceous materials by photons or atomic oxygen in the surface layers of the protoplanetary disk. Under photochemically generated favorable conditions, these reactions can deplete the near-surface abundance of carbon grains and polycyclic aromatic hydrocarbons by several orders of magnitude on short timescales relative to the lifetime of the disk out to radii of ~20-100+ au from the central star depending on the form of refractory carbon present. Due to the reliance of destruction mechanisms on a high influx of photons, the extent of refractory carbon depletion is quite sensitive to the disk's internal radiation field. Dust transport within the disk is required to affect the composition of the midplane. In our current model of a passive, constant-alpha disk, where alpha = 0.01, carbon grains can be turbulently lofted into the destructive surface layers and depleted out to radii of ~3-10 au for 0.1-1 um grains. Smaller grains can be cleared out of the planet-forming region completely. Destruction may be more effective in an actively accreting disk or when considering individual grain trajectories in non-idealized disks.

We report the first detection and high angular resolution (1.8" $\times$ 1.1") imaging of acetic acid (CH$_3$COOH) and gGg$^{\prime}$--ethylene glycol (gGg$^{\prime}$(CH$_2$OH)$_2$) towards the Orion Kleinmann--Low nebula. The observations were carried out at $\sim$1.3mm with ALMA during the Cycle~2. A notable result is that the spatial distribution of the acetic acid and ethylene glycol emission differs from that of the other O-bearing molecules within Orion-KL. Indeed, while the typical emission of O-bearing species harbors a morphology associated with a "V-shape" linking the Hot Core region to the Compact Ridge (with an extension towards the BN object), that of acetic acid and ethylene glycol mainly peaks at about 2" southwest from the hot core region (near sources I and n). We find that the measured CH$_3$COOH:aGg$^{\prime}$(CH$_2$OH)$_2$ and CH$_3$COOH:gGg$^{\prime}$(CH$_2$OH)$_2$ ratios differ from the ones measured towards the low-mass protostar IRAS 16293--2422 by more than one order of magnitude. Our best hypothesis to explain these findings is that CH$_3$COOH, aGg$^{\prime}$(CH$_2$OH)$_2$ and gGg$^{\prime}$(CH$_2$OH)$_2$ are formed on the icy-surface of grains and then released into the gas-phase, via co-desorption with water, due to a bullet of matter ejected during the explosive event that occurred in the heart of the Nebula about 500-700 years ago.

Understanding the formation of wide binary systems of very low mass stars (M $\le$ 0.1 Msun) is challenging. The most obvious route is via widely separated low-mass collapsing fragments produced through turbulent fragmentation of a molecular core. However, close binaries/multiples from disk fragmentation can also evolve to wide binaries over a few initial crossing times of the stellar cluster through tidal evolution. Finding an isolated low mass wide binary system in the earliest stage of formation, before tidal evolution could occur, would prove that turbulent fragmentation is a viable mechanism for (very) low mass wide binaries. Here we report high resolution ALMA observations of a known wide-separation protostellar binary, showing that each component has a circumstellar disk. The system is too young to have evolved from a close binary and the disk axes are misaligned, providing strong support for the turbulent fragmentation model. Masses of both stars are derived from the Keplerian rotation of the disks; both are very low mass stars.

We report the first detection of a substantial brightening event in an isotopologue of a key molecular ion, HCO$^+$, within a protoplanetary disk of a T Tauri star. The H$^{13}$CO$^+$ $J=3-2$ rotational transition was observed three times toward IM Lup between July 2014 and May 2015 with the Atacama Large Millimeter Array. The first two observations show similar spectrally integrated line and continuum fluxes, while the third observation shows a doubling in the disk integrated $J=3-2$ line flux compared to the continuum, which does not change between the three epochs. We explore models of an X-ray active star irradiating the disk via stellar flares, and find that the optically thin H$^{13}$CO$^+$ emission variation can potentially be explained via X-ray driven chemistry temporarily enhancing the HCO$^+$ abundance in the upper layers of the disk atmosphere during large or prolonged flaring events. If the HCO$^+$ enhancement is indeed caused by a X-ray flare, future observations should be able to spatially resolve these events and potentially enable us to watch the chemical aftermath of the high-energy stellar radiation propagating across the face of protoplanetary disks, providing a new pathway to explore ionization physics and chemistry, including electron density, in disks.

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

The initial mass distribution in the solar nebula is a critical input to planet formation models that seek to reproduce today's Solar System. Traditionally, constraints on the gas mass distribution are derived from observations of the dust emission from disks, but this approach suffers from large uncertainties in grain growth and gas-to-dust ratio. On the other hand, previous observations of gas tracers only probe surface layers above the bulk mass reservoir. Here we present the first partially spatially resolved observations of the $^{13}$C$^{18}$O J=3-2 line emission in the closest protoplanetary disk, TW Hya, a gas tracer that probes the bulk mass distribution. Combining it with the C$^{18}$O J=3-2 emission and the previously detected HD J=1-0 flux, we directly constrain the mid-plane temperature and optical depths of gas and dust emission. We report a gas mass distribution of 13$^{+8}_{-5}\times$(R/20.5AU)$^{-0.9^{+0.4}_{-0.3}}$ g cm$^{-2}$ in the expected formation zone of gas and ice giants (5-21AU). We find the total gas/millimeter-sized dust mass ratio is 140 in this region, suggesting that at least 2.4M_earth of dust aggregates have grown to >centimeter sizes (and perhaps much larger). The radial distribution of gas mass is consistent with a self-similar viscous disk profile but much flatter than the posterior extrapolation of mass distribution in our own and extrasolar planetary systems.

We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter and a Super-Earth. We study the dependency of the final orbital radius, mass, and atmospheric chemistry of planets forming in disk models with initial disk masses that vary by 0.02 $M_\odot$ above and below our fiducial model ($M_{disk,0} = 0.1 ~M_\odot$). We compute C/O and C/N for the atmospheres formed in our 3 models and find that C/O$_{\rm planet}\sim$ C/O$_{\rm disk}$, which does not vary strongly between different planets formed in our model. The nitrogen content of atmospheres can vary in planets that grow in different disk models. These differences are related to the formation history of the planet, the time and location that the planet accretes its atmosphere, and are encoded in the bulk abundance of NH$_3$. These results suggest that future observations of atmospheric NH$_3$ and an estimation of the planetary C/O and C/N can inform the formation history of particular planetary systems.

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 high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure.

We report the observation and physical characterization of the possible dwarf planet \UZ\ ("DeeDee"), a dynamically detached trans-Neptunian object discovered at 92 AU. This object is currently the second-most distant known trans-Neptunian object with reported orbital elements, surpassed in distance only by the dwarf planet Eris. The object was discovered with an $r$-band magnitude of 23.0 in data collected by the Dark Energy Survey between 2014 and 2016. Its 1140-year orbit has $(a,e,i) = (109~\mathrm{AU}, 0.65, 26.8^{\circ})$. It will reach its perihelion distance of 38 AU in the year 2142. Integrations of its orbit show it to be dynamically stable on Gyr timescales, with only weak interactions with Neptune. We have performed followup observations with ALMA, using 3 hours of on-source integration time to measure the object's thermal emission in the Rayleigh-Jeans tail. The signal is detected at 7$\sigma$ significance, from which we determine a $V$-band albedo of $13.1^{+3.3}_{-2.4}\mathrm{(stat)}^{+2.0}_{-1.4}\mathrm{(sys)}$ percent and a diameter of $635^{+57}_{-61}\mathrm{(stat)}^{+32}_{-39}\mathrm{(sys)}$~km, assuming a spherical body with uniform surface properties.

We present 2.5-square-degree C$_{2}$H N=1-0 and N$_2$H$^+$ J=1-0 maps of the $\rho$ Ophiuchi molecular cloud complex. These are the first large-scale maps of the $\rho$ Ophiuchi molecular cloud complex with these two tracers. The C$_{2}$H emission is spatially more extended than the N$_2$H$^+$ emission. One faint N$_2$H$^+$ clump Oph-M and one C$_{2}$H ring Oph-RingSW are identified for the first time. The observed C$_{2}$H to N$_{2}$H$^{+}$ abundance ratio ([C$_{2}$H]/[N$_{2}$H$^{+}$]) varies between 5 and 110. We modeled the C$_{2}$H and N$_2$H$^+$ abundances with 1-D chemical models which show a clear decline of [C$_2$H]/[N$_2$H$^+$] with chemical age. Such an evolutionary trend is little affected by temperatures when they are below 40 K. At high density (n$_H$ $>$ 10$^5$ cm$^{-3}$), however, the time it takes for the abundance ratio to drop at least one order of magnitude becomes less than the dynamical time (e.g., turbulence crossing time $\rm \sim$10$^5$ years). The observed [C$_2$H]/[N$_2$H$^+$] difference between L1688 and L1689 can be explained by L1688 having chemically younger gas in relatively less dense regions. The observed [C$_{2}$H]/[N$_{2}$H$^{+}$] values are the results of time evolution, accelerated at higher densities. For the relative low density regions in L1688 where only C$_2$H emission was detected, the gas should be chemically younger.

Many aspects of the Atacama Large Millimeter Array (ALMA) instrument are still unknown due to its young age. One such aspect is the true nature of the primary beam of each baseline, and how changes to the individual primary beams affect astronomical observations when said changes are ignored during imaging. This paper aims to create a more thorough understanding of the strengths and weaknesses of ALMA through realistic modeling of the primary beams and simulated observations, which in turn can inform the user of the necessity of implementing more computationally costly algorithms, such as A-Projection, and when simpler, quicker algorithms will suffice. We quantify our results by examining the dynamic range of each observation, along with the ability to reconstruct the Stokes I amplitude of the test sources. These tests conclude that for dynamic ranges of less than 1000, for point sources and sources much smaller than the main lobe of the primary beam, the accuracy of the primary beam model beyond the physical size of the aperture simply doesn't matter. In observations of large extended sources, deconvolution errors dominate the reconstructed images and the individual primary beam errors were indistinguishable from each other.

The relative populations in rotational transitions of CO can be useful for inferring gas conditions and excitation mechanisms at work in the interstellar medium. We present CO emission lines from rotational transitions observed with Herschel/HIFI in the star-forming cores Orion S, Orion KL, Sgr B2(M), and W49N. Integrated line fluxes from these observations are combined with those from Herschel/PACS observations of the same sources to construct CO spectral line energy distributions (SLEDs) from $5\leq J_u\leq48$. These CO SLEDs are compared to those reported in other galaxies, with the intention of empirically determining which mechanisms dominate excitation in such systems. We find that CO SLEDs in Galactic star-forming cores cannot be used to reproduce those observed in other galaxies, although the discrepancies arise primarily as a result of beam filling factors. The much larger regions sampled by the Herschel beams at distances of several Mpc contain significant amounts of cooler gas which dominate the extragalactic CO SLEDs, in contrast to observations of Galactic star-forming regions which are focused specifically on cores containing primarily hot molecular gas.

Photon Dominated Regions (PDRs) are interfaces between the mainly ionized and mainly molecular material around young massive stars. Analysis of the physical and chemical structure of such regions traces the impact of far-ultraviolet radiation of young massive stars on their environment. We present results on the physical and chemical structure of the prototypical high UV-illumination edge-on Orion Bar PDR from an unbiased spectral line survey with a wide spectral coverage. A spectral scan from 480-1250 GHz and 1410-1910 GHz at 1.1 MHz resolution was obtained by the HIFI instrument onboard the Herschel Space Observatory. For molecules with multiple transitions we used rotational diagrams to obtain excitation temperatures and column densities. For species with a single detected transition we used an optically thin LTE approximation. In case of species with available collisional rates, we also performed a non-LTE analysis to obtain kinetic temperatures, H2 volume densities, and column densities. About 120 lines corresponding to 29 molecules (including isotopologues) have been detected in the Herschel/HIFI line survey, including 11 transitions of CO, 7 transitions of 13CO, 6 transitions of C18O, 10 transitions of H2CO, and 6 transitions of H2O. Most species trace kinetic temperatures in the range between 100 and 150 K and H2 volume densities in the range between 10^5 and 10^6 cm^-3. The species with temperatures and / or densities outside of this range include the H2CO transitions tracing a very high temperature (315 K) and density (1.4x10^6 cm^-3) component and SO corresponding to the lowest temperature (56 K) measured as a part of this line survey. The observed lines/species reveal a range of physical conditions (gas density /temperature) involving structures at high density / high pressure, obsoleting the traditional 'clump / interclump' picture of the Orion Bar.

The compositions of nascent planets depend on the compositions of their birth disks. In particular, the elemental compositions of Gas Giant gaseous envelopes depend on the elemental composition of the disk gas from which the envelope is accreted. Previous models demonstrated that sequential freeze-out of O and C-bearing volatiles in disks will result in an supersolar C/O ratios and subsolar C/H ratios in the gas between water and CO snowlines. This result does not take into account, however, the expected grain growth and radial drift of pebbles in disks, and the accompanying re-distribution of volatiles from the outer to the inner disk. Using a toy model we demonstrate that when drift is considered, CO is enhanced between the water and CO snowline, resulting in both supersolar C/O and C/H ratios in the disk gas in the Gas Giant formation zone. This result appears robust to the details of the disk model as long as there is substantial pebble drift across the CO snowline, and the efficiency of CO vapor diffusion is limited. Gas Giants that accrete their gaseous envelopes exterior to the water snowline and do not experience substantial core-envelope mixing, may thus present both superstellar C/O and C/H ratios in their atmospheres. Pebble drift will also affect the nitrogen and noble gas abundances in the planet forming zones, which may explain some of Jupiter's peculiar abundance patterns.

The processes that govern the evolution of dust and water (in the form of vapor or ice) in protoplanetary disks are intimately connected. We have developed a model that simulates dust coagulation, dust dynamics (settling, turbulent mixing), vapor diffusion, and condensation/sublimation of volatiles onto grains in a vertical column of a protoplanetary disk. We employ the model to study how dust growth and dynamics influence the vertical distribution of water vapor and water ice in the region just outside the radial snowline. Our main finding is that coagulation (boosted by the enhanced stickiness of icy grains) and the ensuing vertical settling of solids results in water vapor being depleted, but not totally removed, from the region above the snowline on a timescale commensurate with the vertical turbulent mixing timescale. Depending on the strength of the turbulence and the temperature, the depletion can reach factors of up to ${\sim}50$ in the disk atmosphere. In our isothermal column, this vapor depletion results in the vertical snowline moving closer to the midplane (by up to 2 gas scale heights) and the gas-phase C/O ratio above the vertical snowline increasing. Our findings illustrate the importance of dynamical effects and the need for understanding coevolutionary dynamics of gas and solids in planet-forming environments.

We report observations of resolved C2H emission rings within the gas-rich protoplanetary disks of TW Hya and DM Tau using the Atacama Large Millimeter Array (ALMA). In each case the emission ring is found to arise at the edge of the observable disk of mm-sized grains (pebbles) traced by (sub)mm-wave continuum emission. In addition, we detect a C3H2 emission ring with an identical spatial distribution to C2H in the TW Hya disk. This suggests that these are hydrocarbon rings (i.e. not limited to C2H). Using a detailed thermo-chemical model we show that reproducing the emission from C2H requires a strong UV field and C/O > 1 in the upper disk atmosphere and outer disk, beyond the edge of the pebble disk. This naturally arises in a disk where the ice-coated dust mass is spatially stratified due to the combined effects of coagulation, gravitational settling and drift. This stratification causes the disk surface and outer disk to have a greater permeability to UV photons. Furthermore the concentration of ices that transport key volatile carriers of oxygen and carbon in the midplane, along with photochemical erosion of CO, leads to an elemental C/O ratio that exceeds unity in the UV-dominated disk. Thus the motions of the grains, and not the gas, lead to a rich hydrocarbon chemistry in disk surface layers and in the outer disk midplane.

How protostars accrete mass is one of the fundamental problems of star formation. High column densities and complex kinematical structures make direct observations challenging and they only provide a snapshot. Chemical tracers provide an interesting alternative to characterise the infall histories of protostars. Previous observations of H13CO+ towards the low-mass protostar IRAS15398-3359 showed a depression in the abundance. This is a sign of destruction of HCO+ by an enhanced presence of gaseous water in an extended region, possibly related to a recent burst in the accretion. Direct observations of water vapour can determine the exact extent of the emission and confirm the hypothesis that HCO+ is indeed a good tracer of the water snow-line. IRAS15398 was observed using ALMA at 0.5" resolution. Maps of HDO(101-000) and H218O(414-321) were taken simultaneously with observations of the CS(8-7) and N2H+(5-4) lines and continuum at 0.65 and 0.75 mm. The maps were interpreted using dust radiative transfer calculations of the protostellar infalling envelope with an outflow cavity. HDO is clearly detected and extended over the scales of the H13CO+ depression, although it is displaced by ~500 AU in the direction of the outflow. H218O is tentatively detected towards the red-shifted outflow lobe, but otherwise it is absent from the mapped region, which suggests that temperatures are low. Based on the temperature structure obtained from dust radiative transfer models, we conclude that the water was most likely released from the grains in an extended hour-glass configuration during a recent accretion burst. HDO is only detected in the region closest to the protostar, at distances of up to 500 AU. These signatures can only be explained if the luminosity has recently been increased by orders of magnitudes. Additionally, the densities in the outflow cones must be sufficiently low.

We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the HIFI instrument aboard the \textit{Herschel Space Observatory}, covering the frequency range 480 to 1900 GHz. We detect 685 spectral lines with S/N $>$ 3$\sigma$, originating from 52 different molecular and atomic species. We model each of the detected species assuming conditions of Local Thermodynamic Equilibrium. This analysis provides an estimate of the physical conditions of Orion South (column density, temperature, source size, \& V$_{LSR}$). We find evidence for three different cloud components: a cool (T$_{ex} \sim 20-40$ K), spatially extended ($> 60"$), and quiescent ($\Delta V_{FWHM} \sim 4$ km s $^{-1}$) component; a warmer (T$_{ex} \sim 80-100$ K), less spatially extended ($\sim 30"$), and dynamic ($\Delta V_{FWHM} \sim 8$ km s $^{-1}$) component, which is likely affected by embedded outflows; and a kinematically distinct region (T$_{ex}$ $>$ 100 K; V$_{LSR}$ $\sim$ 8 km s $^{-1}$), dominated by emission from species which trace ultraviolet irradiation, likely at the surface of the cloud. We find little evidence for the existence of a chemically distinct "hot core" component, likely due to the small filling factor of the hot core or hot cores within the \textit{Herschel} beam. We find that the chemical composition of the gas in the cooler, quiescent component of Orion South more closely resembles that of the quiescent ridge in Orion-KL. The gas in the warmer, dynamic component, however, more closely resembles that of the Compact Ridge and Plateau regions of Orion-KL, suggesting that higher temperatures and shocks also have an influence on the overall chemistry of Orion South.

This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality.

The CH$^+$ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigation CH$^+$ formation with the first velocity-resolved spectral mapping of the CH$^+$ $J=1-0, 2-1$ rotational transitions, three sets of CH $\Lambda$-doubled triplet lines, $^{12}$C$^+$ and $^{13}$C$^+$, and CH$_3$OH 835~GHz E-symmetry Q branch transitions, obtained with Herschel/HIFI over $\approx$12 arcmin$^2$ centered on the Orion BN/KL source. We present the spatial morphologies and kinematics, cloud boundary conditions, excitation temperatures, column densities, and $^{12}$C$^+$ optical depths. Emission from C$^+$, CH$^+$, and CH is indicated to arise in the diluted gas, outside of the explosive, dense BN/KL outflow. Our models show that UV-irradiation provides favorable conditions for steady-state production of CH$^+$ in this environment. Surprisingly, no spatial or kinematic correspondences of these species are found with H$_2$ S(1) emission tracing shocked gas in the outflow. We propose that C$^+$ is being consumed by rapid production of CO to explain the lack of C$^+$ and CH$^+$ in the outflow, and that fluorescence provides the reservoir of H$_2$ excited to higher ro-vibrational and rotational levels. Hence, in star-forming environments containing sources of shocks and strong UV radiation, a description of CH$^+$ formation and excitation conditions is incomplete without including the important --- possibly dominant --- role of UV irradiation.

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

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

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.

CO is widely used as a tracer of molecular gas. However, there is now mounting evidence that gas phase carbon is depleted in the disk around TW Hya. Previous efforts to quantify this depletion have been hampered by uncertainties regarding the radial thermal structure in the disk. Here we present resolved ALMA observations of 13CO 3-2, C18O 3-2, 13CO 6-5, and C18O 6-5 emission in TW Hya, which allow us to derive radial gas temperature and gas surface density profiles, as well as map the CO abundance as a function of radius. These observations provide a measurement of the surface CO snowline at ~30 AU and show evidence for an outer ring of CO emission centered at 53 AU, a feature previously seen only in less abundant species. Further, the derived CO gas temperature profile constrains the freeze-out temperature of CO in the warm molecular layer to < 21 K. Combined with the previous detection of HD 1-0, these data constrain the surface density of the warm H2 gas in the inner ~30 AU. We find that CO is depleted by two orders of magnitude from R=10-60 AU, with the small amount of CO returning to the gas phase inside the surface CO snowline insufficient to explain the overall depletion. Finally, this new data is used in conjunction with previous modeling of the TW Hya disk to constrain the midplane CO snowline to 17-23 AU.

The all-Galaxy CO survey of Dame, Hartmann, & Thaddeus (2001) is by far the most uniform, large-scale Galactic CO survey. Using a dendrogram-based decomposition of this survey, we present a catalog of 1064 massive molecular clouds throughout the Galactic plane. This catalog contains $2.5 \times 10^8$ solar masses, or $25^{+10.7}_{-5.8} \%$ of the Milky Way's estimated H$_2$ mass. We track clouds in some spiral arms through multiple quadrants. The power index of Larson's first law, the size-linewidth relation, is consistent with 0.5 in all regions - possibly due to an observational bias - but clouds in the inner Galaxy systematically have significantly (~ 30%) higher linewidths at a given size, indicating that their linewidths are set in part by Galactic environment. The mass functions of clouds in the inner Galaxy versus the outer Galaxy are both qualitatively and quantitatively distinct. The inner Galaxy mass spectrum is best described by a truncated power-law with a power index of $\gamma=-1.6\pm0.1$ and an upper truncation mass $M_0 = (1.0 \pm 0.2) \times 10^7 M_\odot$, while the outer Galaxy mass spectrum is better described by a non-truncating power law with $\gamma=-2.2\pm0.1$ and an upper mass $M_0 = (1.5 \pm 0.5) \times 10^6 M_\odot$, indicating that the inner Galaxy is able to form and host substantially more massive GMCs than the outer Galaxy. Additionally, we have simulated how the Milky Way would appear in CO from extragalactic perspectives, for comparison with CO maps of other galaxies.

We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed.

Deuterium-to-hydrogen (D/H) enrichments in molecular species provide clues about their original formation environment. The organic materials in primitive solar system bodies have generally higher D/H ratios and show greater D/H variation when compared to D/H in solar system water. We propose this difference arises at least in part due to 1) the availability of additional chemical fractionation pathways for organics beyond that for water, and 2) the higher volatility of key carbon reservoirs compared to oxygen. We test this hypothesis using detailed disk models, including a sophisticated, new disk ionization treatment with a low cosmic ray ionization rate, and find that disk chemistry leads to higher deuterium enrichment in organics compared to water, helped especially by fractionation via the precursors CH$_2$D$^+$/CH$_3^+$. We also find that the D/H ratio in individual species varies significantly depending on their particular formation pathways. For example, from $\sim20-40$ AU, CH$_4$ can reach $\rm{D/H\sim2\times10^{-3}}$, while D/H in CH$_3$OH remains locally unaltered. Finally, while the global organic D/H in our models can reproduce intermediately elevated D/H in the bulk hydrocarbon reservoir, our models are unable to reproduce the most deuterium-enriched organic materials in the solar system, and thus our model requires some inheritance from the cold interstellar medium from which the Sun formed.

An unsolved problem in step-wise core-accretion planet formation is that rapid radial drift in gas-rich protoplanetary disks should drive mm-/meter-sized particles inward to the central star before large bodies can form. One promising solution is to confine solids within small scale structures. Here we investigate dust structures in the (sub)mm continuum emission of four disks (TW Hya, HL Tau, HD 163296 and DM Tau), a sample of disks with the highest spatial resolution ALMA observations to date. We retrieve the surface brightness distributions using synthesized images and fitting visibilities with analytical functions. We find that the continuum emission of the four disks is ~axi-symmetric but rich in 10-30AU-sized radial structures, possibly due to physical gaps, surface density enhancements or localized dust opacity variations within the disks. These results suggest that small scale axi-symmetric dust structures are likely to be common, as a result of ubiquitous processes in disk evolution and planet formation. Compared with recent spatially resolved observations of CO snowlines in these same disks, all four systems show enhanced continuum emission from regions just beyond the CO condensation fronts, potentially suggesting a causal relationship between dust growth/trapping and snowlines.

We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.

The condensation fronts (snow lines) of H2O, CO and other abundant volatiles in the midplane of a protoplanetary disk affect several aspects of planet formation. Locating the CO snow line, where the CO gas column density is expected to drop substantially, based solely on CO emission profiles is challenging. This has prompted an exploration of chemical signatures of CO freeze-out. We present ALMA Cycle 1 observations of the N2H+ J=3-2 and DCO+ J=4-3 emission lines toward the disk around the Herbig Ae star HD~163296 at ~0.5" (60 AU) resolution, and evaluate their utility as tracers of the CO snow line location. The N2H+ emission is distributed in a ring with an inner radius at 90 AU, corresponding to a midplane temperature of 25 K. This result is consistent with a new analysis of optically thin C18O data, which implies a sharp drop in CO abundance at 90 AU. Thus N2H+ appears to be a robust tracer of the midplane CO snow line. The DCO+ emission also has a ring morphology, but neither the inner nor the outer radius coincides with the CO snow line location of 90 AU, indicative of a complex relationship between DCO+ emission and CO freeze-out in the disk midplane. Compared to TW Hya, CO freezes out at a higher temperature in the disk around HD 163296 (25 vs. 17 K in the TW Hya disk), perhaps due to different ice compositions. This highlights the importance of actually measuring the CO snow line location, rather than assuming a constant CO freeze-out temperature for all disks.

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

Argonium has recently been detected as a ubiquitous molecule in our Galaxy. Model calculations indicate that its abundance peaks at molecular fractions in the range of 1E-4 to 1E-3 and that the observed column densities require high values of the cosmic ray ionization rate. Therefore, this molecular cation may serve as an excellent tracer of the very diffuse interstellar medium (ISM), as well as an indicator of the cosmic ray ionization rate. We attempted to detect ArH+ in extragalactic sources to evaluate its diagnostic power as a tracer of the almost purely atomic ISM in distant galaxies. We obtained ALMA observations of a foreground galaxy at z = 0.89 in the direction of the lensed blazar PKS 1830-211. Two isotopologs of argonium, 36ArH+ and 38ArH+, were detected in absorption along two different lines of sight toward PKS 1830-211, known as the SW and NE images of the background blazar. The argonium absorption is clearly enhanced on the more diffuse line of sight (NE) compared to other molecular species. The isotopic ratio 36Ar/38Ar is 3.46 +- 0.16 toward the SW image, i.e., significantly lower than the solar value of 5.5. Our results demonstrate the suitability of argonium as a tracer of the almost purely atomic, diffuse ISM in high-redshift sources. The evolution of the isotopic ratio with redshift may help to constrain nucleosynthetic scenarios in the early Universe.

We use the C/N ratio as a monitor of the delivery of key ingredients of life to nascent terrestrial worlds. Total elemental C and N contents, and their ratio, are examined for the interstellar medium, comets, chondritic meteorites and terrestrial planets; we include an updated estimate for the Bulk Silicate Earth (C/N = 49.0 +/- 9.3). Using a kinetic model of disk chemistry, and the sublimation/condensation temperatures of primitive molecules, we suggest that organic ices and macro-molecular (refractory or carbonaceous dust) organic material are the likely initial C and N carriers. Chemical reactions in the disk can produce nebular C/N ratios of ~1-12, comparable to those of comets and the low end estimated for planetesimals. An increase of the C/N ratio is traced between volatile-rich pristine bodies and larger volatile-depleted objects subjected to thermal/accretional metamorphism. The C/N ratios of the dominant materials accreted to terrestrial planets should therefore be higher than those seen in carbonaceous chondrites or comets. During planetary formation, we explore scenarios leading to further volatile loss and associated C/N variations owing to core formation and atmospheric escape. Key processes include relative enrichment of nitrogen in the atmosphere and preferential sequestration of carbon by the core. The high C/N BSE ratio therefore is best satisfied by accretion of thermally processed objects followed by large-scale atmospheric loss. These two effects must be more profound if volatile sequestration in the core is effective. The stochastic nature of these processes hints that the surface/atmospheric abundances of biosphere-essential materials will likely be variable.

How complex organic - and potentially prebiotic - molecules are formed in regions of low- and high-mass star-formation remains a central question in astrochemistry. In particular, with just a few sources studied in detail, it is unclear what role environment plays in complex molecule formation. In this light, a comparison of relative abundances of related species between sources might be useful to explain observed differences. We seek to measure the relative abundance between three important complex organic molecules, ethylene glycol ((CH$_2$OH)$_2$), glycolaldehyde (CH$_2$OHCHO) and methyl formate (HCOOCH$_3$), toward high-mass protostars and thereby provide additional constraints on their formation pathways. We use IRAM 30-m single dish observations of the three species toward two high-mass star-forming regions - W51/e2 and G34.3+0.2 - and report a tentative detection of (CH2OH)2 toward both sources. Assuming that (CH$_2$OH)$_2$, CH$_2$OHCHO and HCOOCH$_3$ spatially coexist, relative abundance ratios, HCOOCH$_3$/(CH$_2$OH)$_2$, of 31 and 35 are derived for G34.3+0.2 and W51/e2, respectively. CH$_2$OHCHO is not detected, but the data provide lower limits to the HCOOCH$_3$/CH$_2$OHCHO abundance ratios of $\ge$193 for G34.3+0.2 and $\ge$550 for W51/e2. A comparison of these results to measurements from various sources in the literature indicates that the source luminosities may be correlated with the HCOOCH$_3$/(CH$_2$OH)$_2$ and HCOOCH$_3$/CH$_2$OHCHO ratios. This apparent correlation may be a consequence of the relative timescales each source spend at different temperatures-ranges in their evolution. Furthermore, we obtain lower limits to the ratio of (CH$_2$OH)$_2$/CH2OHCHO for G34.3+0.2 ($\ge$6) and W51/e2 ($\ge$16). This result confirms that a high (CH$_2$OH)$_2$/CH$_2$OHCHO abundance ratio is not a specific property of comets, as previously speculated.

We report the first high angular resolution imaging (3.4\arcsec $\times$ 3.0\arcsec) of deuterated formaldehyde (HDCO) toward Orion--KL, carried out with the Submillimeter Array (SMA). We find that the spatial distribution of the formaldehyde emission systematically differs from that of methanol: while methanol is found towards the inner part of the region, HDCO is found in colder gas that wraps around the methanol emission on four sides. The HDCO/H$_2$CO ratios are determined to be 0.003--0.009 within the region, up to an order of magnitude higher than the D/H measured for methanol. These findings strengthen the previously suggested hypothesis that there are differences in the chemical pathways leading to HDCO (via deuterated gas phase chemistry) and deuterated methanol (through conversion of formaldehyde into methanol on the surface of icy grain mantles).

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

We present chemical implications arising from spectral models fit to the Herschel/HIFI spectral survey toward the Orion Kleinmann-Low nebula (Orion KL). We focus our discussion on the eight complex organics detected within the HIFI survey utilizing a novel technique to identify those molecules emitting in the hottest gas. In particular, we find the complex nitrogen bearing species CH$_{3}$CN, C$_{2}$H$_{3}$CN, C$_{2}$H$_{5}$CN, and NH$_{2}$CHO systematically trace hotter gas than the oxygen bearing organics CH$_{3}$OH, C$_{2}$H$_{5}$OH, CH$_{3}$OCH$_{3}$, and CH$_{3}$OCHO, which do not contain nitrogen. If these complex species form predominantly on grain surfaces, this may indicate N-bearing organics are more difficult to remove from grain surfaces than O-bearing species. Another possibility is that hot (T$_{\rm kin}$$\sim$300 K) gas phase chemistry naturally produces higher complex cyanide abundances while suppressing the formation of O-bearing complex organics. We compare our derived rotation temperatures and molecular abundances to chemical models, which include gas-phase and grain surface pathways. Abundances for a majority of the detected complex organics can be reproduced over timescales $\gtrsim$ 10$^{5}$ years, with several species being under predicted by less than 3$\sigma$. Derived rotation temperatures for most organics, furthermore, agree reasonably well with the predicted temperatures at peak abundance. We also find that sulfur bearing molecules which also contain oxygen (i.e. SO, SO$_{2}$, and OCS) tend to probe the hottest gas toward Orion KL indicating the formation pathways for these species are most efficient at high temperatures.

An abundance decrease in carbon- and oxygen-bearing species relative to dust has been frequently found in planet-forming disks, which can be attributed to an overall reduction of gas mass. However, in the case of TW Hya, the only disk with gas mass measured directly with HD rotational lines, the inferred gas mass ($\lesssim$0.005 solar mass) is significantly below the directly measured value ($\gtrsim$0.05 solar mass). We show that this apparent conflict can be resolved if the elemental abundances of carbon and oxygen are reduced in the upper layers of the outer disk but are normal elsewhere (except for a possible enhancement of their abundances in the inner disk). The implication is that in the outer disk, the main reservoir of the volatiles (CO, water, ...) resides close to the midplane, locked up inside solid bodies that are too heavy to be transported back to the atmosphere by turbulence. An enhancement in the carbon and oxygen abundances in the inner disk can be caused by inward migration of these solid bodies. This is consistent with estimates based on previous models of dust grain dynamics. Indirect measurements of the disk gas mass and disk structure from species such as CO will thus be intertwined with the evolution of dust grains, and possibly also with the formation of planetesimals.

We examine changes in the molecular abundances resulting from increased heating due to a self-luminous planetary companion embedded within a narrow circumstellar disk gap. Using 3D models that include stellar and planetary irradiation, we find that luminous young planets locally heat up the parent circumstellar disk by many tens of Kelvin, resulting in efficient thermal desorption of molecular species that are otherwise locally frozen out. Furthermore, the heating is deposited over large regions of the disk, $\pm5$ AU radially and spanning $\lesssim60^\circ$ azimuthally. From the 3D chemical models, we compute rotational line emission models and full ALMA simulations, and find that the chemical signatures of the young planet are detectable as chemical asymmetries in $\sim10h$ observations. HCN and its isotopologues are particularly clear tracers of planetary heating for the models considered here, and emission from multiple transitions of the same species is detectable, which encodes temperature information in addition to possible velocity information from the spectra itself. We find submillimeter molecular emission will be a useful tool to study gas giant planet formation in situ, especially beyond $R\gtrsim10$ AU.

Many of the transformative processes in the Universe have taken place in regions obscured by dust, and are best studied with far-IR spectroscopy. We present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600 micron band. CALISTO will enable NASA and the world to study the rise of heavy elements in the Universe's first billion years, chart star formation and black hole growth in dust-obscured galaxies through cosmic time, and conduct a census of forming planetary systems in our region of the Galaxy. CALISTO will capitalize on rapid progress in both format and sensitivity of far-IR detectors. Arrays with a total count of a few 100,000 detector pixels will form the heart of a suite of imaging spectrometers in which each detector reaches the photon background limit. This document contains a large overview paper on CALISTO, as well as six 2-3 page scientific white papers, all prepared in response to NASA's Cosmic Origins Program Analysis Group (COPAG's) request for input on future mission concepts. The Far-IR Science Interest Group will meet from 3-5 June 2015 with the intention of reaching consensus on the architecture for the Far-IR Surveyor mission. This white paper describes one of the architectures to be considered by the community. One or more companion papers will describe alternative architectures.

Water and simple organic molecular ices dominate the mass of solid materials available for planetesimal and planet formation beyond the water snow line. Here we analyze ALMA long baseline 2.9, 1.3 and 0.87 mm continuum images of the young star HL Tau, and suggest that the emission dips observed are due to rapid pebble growth around the condensation fronts of abundant volatile species. Specifically, we show that the prominent innermost dip at 13 AU is spatially resolved in the 0.87 mm image, and its center radius is coincident with the expected mid-plane condensation front of water ice. In addition, two other prominent dips, at distances of 32 and 63 AU, cover the mid-plane condensation fronts of pure ammonia or ammonia hydrates and clathrate hydrates (especially with CO and N$_2$) formed from amorphous water ice. The spectral index map of HL Tau between 1.3 and 0.87 mm shows that the flux ratios inside the dips are statistically larger than those of nearby regions in the disk. This variation can be explained by a model with two dust populations, where most of solid mass resides in a component that has grown into decimeter size scales inside the dips. Such growth is in accord with recent numerical simulations of volatile condensation, dust coagulation and settling.

C$_2$H is one of the first radicals to be detected in the interstellar medium. Its higher rotational transitions have recently become available with the Herschel Space Observatory. We aim to constrain the physical parameters of the C$_2$H emitting gas toward the Orion Bar. We analyse the C$_2$H line intensities measured toward the Orion Bar CO$^+$ Peak and Herschel/HIFI maps of C$_2$H, CH, and HCO$^+$, and a NANTEN map of [CI]. We interpret the observed C$_2$H emission using radiative transfer and PDR models. Five rotational transitions of C$_2$H have been detected in the HIFI frequency range toward the CO$^+$ peak. A single component rotational diagram gives a rotation temperature of ~64 K and a beam-averaged C$_2$H column density of 4$\times$10$^{13}$ cm$^{-2}$. The measured transitions cannot be explained by any single parameter model. According to a non-LTE model, most of the C$_2$H column density produces the lower-$N$ C$_2$H transitions and traces a warm ($T_{\rm{kin}}$ ~ 100-150 K) and dense ($n$(H$_2$)~10$^5$-10$^6$ cm$^{-3}$) gas. A small fraction of the C$_2$H column density is required to reproduce the intensity of the highest-$N$ transitions ($N$=9-8 and N=10-9) originating from a high density ($n$(H$_2$)~5$\times$10$^6$ cm$^{-3}$) hot ($T_{\rm{kin}}$ ~ 400 K) gas. The total beam-averaged C$_2$H column density in the model is 10$^{14}$ cm$^{-2}$. Both the non-LTE radiative transfer model and a simple PDR model representing the Orion Bar with a plane-parallel slab of gas and dust suggest, that C$_2$H cannot be described by a single pressure component, unlike the reactive ion CH$^+$, which was previously analysed toward the Orion Bar CO$^+$ peak. The physical parameters traced by the higher rotational transitions ($N$=6-5,...,10-9) of C$_2$H may be consistent with the edges of dense clumps exposed to UV radiation near the ionization front of the Orion Bar.

The accretion histories of embedded protostars are an integral part of descriptions of their physical and chemical evolution. In particular, are the accretion rates smoothly declining from the earlier toward later stages or in fact characterized by variations such as intermittent bursts? We aim to characterize the impact of possible accretion variations in a sample of embedded protostars by measuring the size of the inner regions of their envelopes where CO is sublimated and relate those to their temperature profiles dictated by their current luminosities. Using observations from the Submillimeter Array we measure the extents of the emission from the C18O isotopologue toward 16 deeply embedded protostars. We compare these measurements to the predicted extent of the emission given the current luminosities of the sources through dust and line radiative transfer calculations. Eight out of sixteen sources show more extended C18O emission than predicted by the models. The modeling shows that the likely culprit for these signatures is sublimation due to increases in luminosities of the sources by about a factor five or more during the recent 10,000 years - the time it takes for CO to freeze-out again on dust grains. For four of those sources the increase would have had to have been a factor 10 or more. The compact emission seen toward the other half of the sample suggests that C18O only sublimates when the temperature exceeds 30 K - as one would expect if CO is mixed with H2O in the grain ice-mantles. The small-number statistics from this survey suggest that protostars undergo significant bursts about once every 20,000 years. This also illustrates the importance of taking the physical evolutionary histories into account for descriptions of the chemical structures of embedded protostars.

Aims: Accretion rates in low-mass protostars can be highly variable in time. Each accretion burst is accompanied by a temporary increase in luminosity, heating up the circumstellar envelope and altering the chemical composition of the gas and dust. This paper aims to study such chemical effects and discusses the feasibility of using molecular spectroscopy as a tracer of episodic accretion rates and timescales. Methods: We simulate a strong accretion burst in a diverse sample of 25 spherical envelope models by increasing the luminosity to 100 times the observed value. Using a comprehensive gas-grain network, we follow the chemical evolution during the burst and for up to 10^5 yr after the system returns to quiescence. The resulting abundance profiles are fed into a line radiative transfer code to simulate rotational spectra of C18O, HCO+, H13CO+, and N2H+ at a series of time steps. We compare these spectra to observations taken from the literature and to previously unpublished data of HCO+ and N2H+ 6-5 from the Herschel Space Observatory. Results: The bursts are strong enough to evaporate CO throughout the envelope, which in turn enhances the abundance of HCO+ and reduces that of N2H+. After the burst, it takes 10^3-10^4 yr for CO to refreeze and for HCO+ and N2H+ to return to normal. The chemical effects of the burst remain visible in the rotational spectra for as long as 10^5 yr after the burst has ended, highlighting the importance of considering luminosity variations when analyzing molecular line observations in protostars. The spherical models are currently not accurate enough to derive robust timescales from single-dish observations. As follow-up work, we suggest that the models be calibrated against spatially resolved observations in order to identify the best tracers to be used for statistically significant source samples.

In this Letter we model the chemistry of DCO$^{+}$ in protoplanetary disks. We find that the overall distribution of the DCO$^{+}$ abundance is qualitatively similar to that of CO but is dominated by thin layer located at the inner disk surface. To understand its distribution, we investigate the different key gas-phase deuteration pathways that can lead to the formation of DCO$^{+}$. Our analysis shows that the recent update in the exothermicity of the reaction involving CH$_2$D$^{+}$ as a parent molecule of DCO$^{+}$ favors deuterium fractionation in warmer conditions. As a result the formation of DCO$^{+}$ is enhanced in the inner warm surface layers of the disk where X-ray ionization occurs. Our analysis points out that DCO$^{+}$ is not a reliable tracer of the CO snow line as previously suggested. We thus predict that DCO$^{+}$ is a tracer of active deuterium and in particular X-ray ionization of the inner disk.

A fraction of the mid-$J$ ($J$= 14--13 to $J$= 24--23) CO emission detected by the \textit{Herschel}/PACS observations of embedded young stellar objects (YSOs) has been attributed to the UV-heated outflow cavity walls. We have applied our newly developed self-consistent models of Photon-Dominated Region (PDR) and non-local thermal equilibrium line Radiative transfer In general Grid (RIG) to the \textit{Herschel} FIR observations of 27 low mass YSOs and one intermediate mass YSO, NGC7129-FIRS2. When the contribution of the hot component (traced by transitions of $J> 24$) is removed, the rotational temperature of the warm component is nearly constant with $\sim250$ K. This can be reproduced by the outflow cavity wall ($n \geq 10^6\, \mathrm{cm}^{-3}$, $\log G_{0}/n \geq-4.5$, $\mathrm{log} G_0\ge 3$, $T_{\rm gas} \ge 300 $K, and X(CO)$ \ge 10^{-5}$) heated by a UV radiation field with a black body temperature of 15,000 K or 10,000 K. However, a shock model combined with an internal PDR will be required to determine the quantitative contribution of a PDR relative to a shock to the mid-$J$ CO emission.

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.

Hot molecular cores (HMCs) are intermediate stages of high-mass star formation and are also known for their rich emission line spectra at (sub-)mm wavebands. The observed spectral feature of HMCs such as total number of emission lines and associated line intensities are also found to vary with evolutionary stages. We developed various 3D models for HMCs guided by the evolutionary scenarios proposed by recent empirical and modeling studies. We then investigated the spatio-temporal variation of temperature and molecular abundances in HMCs by consistently coupling gas-grain chemical evolution with radiative transfer calculations. We explored the effects of varying physical conditions on molecular abundances including density distribution and luminosity evolution of the central protostar(s). The time-dependent temperature structure of the hot core models provides a realistic framework for investigating the spatial variation of ice mantle evaporation as a function of evolutionary timescales. With increasing protostellar luminosity, the water ice evaporation font ($\sim$100K) expands and the spatial distribution of gas phase abundances of these COMs also spreads out. We simulated the synthetic spectra for these models at different evolutionary timescales to compare with observations. A qualitative comparison of the simulated and observed spectra suggests that these self-consistent hot core models can reproduce the notable trends in hot core spectral variation within the typical hot core timescales of 10$^{5}$ year. These models predict that the spatial distribution of various emission line maps will also expand with evolutionary time. The model predictions can be compared with high resolution observation that can probe scales of a few thousand AU in high-mass star forming regions such as from ALMA.[Abridged]

The recent (apparent) passage of the Voyager 1 spacecraft into interstellar space provides us with front-row seats to the complex interplay between the solar wind and the protective surrounding bubble known as heliosphere. The heliosphere extends radially out to $\sim100$ AU from the sun, and within this sphere of influence, the solar wind modulates the incoming flux of galactic cosmic rays (CRs), especially those at low energies. Newly formed stars, which support both strong magnetic fields and winds, are expected to produce analogous regions of CR exclusion, perhaps at elevated levels. Such young stars are encircled by molecular gas-rich disks, and the net removal of CRs from the circumstellar environment significantly reduces the expected CR ionization rate in the disk gas, most likely by many orders-of-magnitude. The loss of ionization reduces disk turbulence, and thereby affects both planet-formation and active chemical processes in the disk. We present models of CR exclusion and explore the implications for turbulence and for predicted chemical abundances. We also discuss means by which ALMA can be used to search for extrasolar heliosphere-analogs around young stars.

We present an observational and theoretical study of the primary ionizing agents (cosmic rays and X-rays) in the TW Hya protoplanetary disk. We use a set of resolved and unresolved observations of molecular ions and other molecular species, encompassing eleven lines total, in concert with a grid of disk chemistry models. The molecular ion constraints comprise new data from the Submillimeter Array on HCO$^+$, acquired at unprecedented spatial resolution, and data from the literature, including ALMA observations of N$_2$H$^+$. We vary the model incident CR flux and stellar X-ray spectra and find that TW Hya's HCO$^+$ and N$_2$H$^+$ emission are best fit by a moderately hard X-ray spectra, as would be expected during the "flaring" state of the star, and a low CR ionization rate, $\zeta_{\rm CR}\lesssim10^{-19}$ s$^{-1}$. This low CR rate is the first indication of the presence of CR exclusion by winds and/or magnetic fields in an actively accreting T Tauri disk system. With this new constraint, our best fit ionization structure predicts a low turbulence "dead-zone" extending from the inner edge of the disk out to $50-65$ AU. This region coincides with an observed concentration of millimeter grains, and we propose that the inner region of TW Hya is a dust (and possibly planet) growth factory as predicted by previous theoretical work.

We present Herschel Space Observatory PACS spectra of T Tauri stars, in which we detect amorphous and crystalline water ice features. Using irradiated accretion disk models, we determine the disk structure and ice abundance in each of the systems. Combining a model-independent comparison of the ice feature strength and disk size with a detailed analysis of the model ice location, we estimate that the ice emitting region is at disk radii >30AU, consistent with a proto-Kuiper belt. Vertically, the ice emits most below the photodesorption zone, consistent with Herschel observations of cold water vapor. The presence of crystallized water ice at a disk location a) colder than its crystallization temperature and b) where it should have been re-amorphized in ~1 Myr suggests that localized generation is occurring; the most likely cause appears to be micrometeorite impact or planetesimal collisions. Based on simple tests with UV models and different ice distributions, we suggest that the SED shape from 20 to 50 micron may probe the location of the water ice snow line in the disk upper layers. This project represents one of the first extra-solar probes of the spatial structure of the cometary ice reservoir thought to deliver water to terrestrial planets.

We present the modeling results of deuterium fractionation of water ice, H2, and the primary deuterium isotopologues of H3+ adopting physical conditions associated with the star and planet formation process. We calculated the deuterium chemistry for a range of gas temperatures (T_gas ~ 10 - 30 K), molecular hydrogen density (n(H2)~ 10^4 - 10^7), and ortho/para ratio (opr) of H2 based on state-to-state reaction rates and explore the resulting fractionation including the formation of a water ice mantle coating grain surfaces. We find that the deuterium fractionation exhibits the expected temperature dependence of large enrichments at low gas temperature. More significantly the inclusion of water ice formation leads to large D/H ratios in water ice (>= 10^-2 at 10 K) but also alters the overall deuterium chemistry. For T < 20 K the implantation of deuterium into ices lowers the overall abundance of HD which reduces the efficiency of deuterium fractionation at high density. In agreement with an earlier study, under these conditions HD may not be the primary deuterium reservoir in the cold dense interstellar medium and H3+ will be the main charge carrier in the dense centers of pre-stellar cores and the protoplanetary disk midplane.

We present ~2'x2' spectral-maps of Orion BN/KL outflows taken with Herschel at ~12'' resolution. For the first time in the far-IR domain, we spatially resolve the emission associated with the bright H2 shocked regions "Peak 1" and "Peak 2" from that of the Hot Core and ambient cloud. We analyze the ~54-310um spectra taken with the PACS and SPIRE spectrometers. More than 100 lines are detected, most of them rotationally excited lines of 12CO (up to J=48-47), H2O, OH, 13CO, and HCN. Peaks 1/2 are characterized by a very high L(CO)/L(FIR)~5x10^{-3} ratio and a plethora of far-IR H2O emission lines. The high-J CO and OH lines are a factor ~2 brighter toward Peak 1 whereas several excited H2O lines are ~50% brighter toward Peak 2. A simplified non-LTE model allowed us to constrain the dominant gas temperature components. Most of the CO column density arises from Tk~200-500 K gas that we associate with low-velocity shocks that fail to sputter grain ice mantles and show a maximum gas-phase H2O/CO~10^{-2} abundance ratio. In addition, the very excited CO (J>35) and H2O lines reveal a hotter gas component (Tk~2500 K) from faster (v_S>25 km/s) shocks that are able to sputter the frozen-out H2O and lead to high H2O/CO>~1 abundance ratios. The H2O and OH luminosities cannot be reproduced by shock models that assume high (undepleted) abundances of atomic oxygen in the preshock gas and/or neglect the presence of UV radiation in the postshock gas. Although massive outflows are a common feature in other massive star-forming cores, Orion BN/KL seems more peculiar because of its higher molecular luminosities and strong outflows caused by a recent explosive event.

The dominant form of nitrogen provided to most solar system bodies is currently unknown, though available measurements show that the detected nitrogen in solar system rocks and ices is depleted with respect to solar abundances and the interstellar medium. We use a detailed chemical/physical model of the chemical evolution of a protoplanetary disk to explore the evolution and abundance of nitrogen-bearing molecules. Based on this model we analyze how initial chemical abundances, provided as either gas or ice during the early stages of disk formation, influence which species become the dominant nitrogen bearers at later stages. We find that a disk with the majority of its initial nitrogen in either atomic or molecular nitrogen is later dominated by atomic and molecular nitrogen as well as NH$_{3}$ and HCN ices, where the dominant species varies with disk radius. When nitrogen is initially in gaseous ammonia, it later becomes trapped in ammonia ice except in the outer disk where atomic nitrogen dominates. For a disk with the initial nitrogen in the form of ammonia ice the nitrogen remains trapped in the ice as NH$_{3}$ at later stages. The model in which most of the initial nitrogen is placed in atomic N best matches the ammonia abundances observed in comets. Furthermore the initial state of nitrogen influences the abundance of N$_{2}$H$^{+}$, which has been detected in protoplanetary disks. Strong N$_{2}$H$^{+}$ emission is found to be indicative of an N$_{2}$ abundance greater than $n_{\mathrm{N_{2}}}/n_{\mathrm{H_{2}}}>10^{-6}$, in addition to tracing the CO snow line. Our models also indicate that NO is potentially detectable, with lower N gas abundances leading to higher NO abundances.

We have surveyed a sample of massive star-forming regions located over a range of distances from the Galactic centre for methyl formate, HCOOCH$_{3}$, and its isotopologues H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$. The observations were carried out with the APEX telescope in the frequency range 283.4--287.4~GHz. Based on the APEX observations, we report tentative detections of the $^{13}$C-methyl formate isotopologue HCOO$^{13}$CH$_{3}$ towards the following four massive star-forming regions: Sgr~B2(N-LMH), NGC~6334~IRS~1, W51 e2 and G19.61-0.23. In addition, we have used the 1~mm ALMA science verification observations of Orion-KL and confirm the detection of the $^{13}$C-methyl formate species in Orion-KL and image its spatial distribution. Our analysis shows that the $^{12}$C/$^{13}$C isotope ratio in methyl formate toward Orion-KL Compact Ridge and Hot Core-SW components (68.4$\pm$10.1 and 71.4$\pm$7.8, respectively) are, for both the $^{13}$C-methyl formate isotopologues, commensurate with the average $^{12}$C/$^{13}$C ratio of CO derived toward Orion-KL. Likewise, regarding the other sources, our results are consistent with the $^{12}$C/$^{13}$C in CO. We also report the spectroscopic characterization, which includes a complete partition function, of the complex H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$ species. New spectroscopic data for both isotopomers H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$, presented in this study, has made it possible to measure this fundamentally important isotope ratio in a large organic molecule for the first time.

Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk. Utilizing a comprehensive treatment of disk ionization, we find that ion-driven deuterium pathways are inefficient, curtailing the disk's deuterated water formation and its viability as the sole source for the solar system's water. This finding implies that if the solar system's formation was typical, abundant interstellar ices are available to all nascent planetary systems.

The ortho-to-para ratio (OPR) of H$_2$O is thought to be sensitive to the temperature of water formation. The OPR of H$_2$O is thus useful to study the formation mechanism of water. We investigate the OPR of water in the Orion PDR (Photon-dominated region), at the Orion Bar and Orion S positions, using data from {\it Herschel}/HIFI. We detect the ground-state lines of ortho- and para-H$_2$$^{18}$O in the Orion Bar and Orion S and we estimate the column densities using LTE and non-LTE methods. Based on our calculations, the ortho-to-para ratio (OPR) in the Orion Bar is 0.1 $-$ 0.5, which is unexpectedly low given the gas temperature of $\sim$ 85 K, and also lower than the values measured for other interstellar clouds and protoplanetary disks. Toward Orion S, our OPR estimate is below 2. This low OPR at 2 positions in the Orion PDR is inconsistent with gas phase formation and with thermal evaporation from dust grains, but it may be explained by photodesorption.

The chemical properties of protoplanetary disks are especially sensitive to their ionization environment. Sources of molecular gas ionization include cosmic rays, stellar X-rays and short-lived radionuclides, each of which varies with location in the disk. This behavior leads to a significant amount of chemical structure, especially in molecular ion abundances, which is imprinted in their submillimeter rotational line emission. Using an observationally motivated disk model, we make predictions for the dependence of chemical abundances on the assumed properties of the ionizing field. We calculate the emergent line intensity for abundant molecular ions and simulate sensitive observations with the Atacama Large Millimeter/Sub-millimeter Array (ALMA) for a disk at D=100 pc. The models readily distinguish between high ionization rates ($\zeta\gtrsim10^{-17}$ s$^{-1}$ per H$_2$) and below, but it becomes difficult to distinguish between low ionization models when $\zeta\lesssim10^{-19}$ s$^{-1}$. We find that \htdp\ emission is not detectable for sub-interstellar CR rates with ALMA (6h integration), and that \ntdp\ emission may be a more sensitive tracer of midplane ionization. HCO$^+$ traces X-rays and high CR rates ($\zeta_{\rm{CR}}\gtrsim10^{-17}$ s$^{-1}$), and provides a handle on the warm molecular ionization properties where CO is present in the gas. Furthermore, species like HCO$^+$, which emits from a wide radial region and samples a large gradient in temperature, can exhibit ring-like emission as a consequence of low-lying rotational level de-excitation near the star. This finding highlights a scenario where rings are not necessarily structural or chemical in nature, but simply a result of the underlying line excitation properties.

Water vapor has been detected in protoplanetary disks. In this work we model the distribution of water vapor in protoplanetary disks with a thermo-chemical code. For a set of parameterized disk models, we calculate the distribution of dust temperature and radiation field of the disk with a Monte Carlo method, and then solve the gas temperature distribution and chemical composition. The radiative transfer includes detailed treatment of scattering by atomic hydrogen and absorption by water of Lyman alpha photons, since the Lyman alpha line dominates the UV spectrum of accreting young stars. In a fiducial model, we find that warm water vapor with temperature around 300 K is mainly distributed in a small and well-confined region in the inner disk. The inner boundary of the warm water region is where the shielding of UV field due to dust and water itself become significant. The outer boundary is where the dust temperature drops below the water condensation temperature. A more luminous central star leads to a more extended distribution of warm water vapor, while dust growth and settling tends to reduce the amount of warm water vapor. Based on typical assumptions regarding the elemental oxygen abundance and the water chemistry, the column density of warm water vapor can be as high as 1E22 cm^{-2}. A small amount of hot water vapor with temperature higher than ~300 K exists in a more extended region in the upper atmosphere of the disk. Cold water vapor with temperature lower than 100 K is distributed over the entire disk, produced by photodesorption of the water ice.

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

A photon dominated region (PDR) is one of the leading candidate mechanisms for the origin of the warm CO gas with near universal $\sim$300~K rotational temperature inferred from the CO emission detected towards embedded protostars by Herschel/PACS. We have developed a PDR model in general coordinates, where we can use the most adequate coordinate system for an embedded protostar having outflow cavity walls, to solve chemistry and gas energetics self-consistently for given UV radiation fields with different spectral shapes. Simple 1D tests and applications show that FIR mid-$J$ ($14\leq J\leq 24$) CO lines are emitted from near the surface of a dense region exposed to high UV fluxes. We apply our model to HH46 and find the UV-heated outflow cavity wall can reproduce the mid-$J$, CO transitions observed by Herschel/PACS. A model with UV radiation corresponding to a blackbody of 10,000~K results in the rotational temperature lower than 300~K, while models with the Draine interstellar radiation field and the 15,000 K blackbody radiation field predict the rotational temperature similar to the observed one.

Given the central role of carbon in the chemistry of life, it is a fundamental question as to how carbon is supplied to the Earth, in what form and when. We provide an accounting of carbon found in solar system bodies, in particular a comparison between the organic content of meteorites and that in identified organics in the dense interstellar medium (ISM). Based on this accounting identified organics created by the chemistry of star formation could contain at most ~15% of the organic carbon content in primitive meteorites and significantly less for cometary organics, which represent the putative contributors to starting materials for the Earth. In the ISM ~30% of the elemental carbon is found in CO, either in the gas or ices, with a typical abundance of ~10^-4 (relative to H2). Recent observations of the TW Hya disk find that the gas phase abundance of CO is reduced by an order of magnitude compared to this value. We explore a solution where the volatile CO is destroyed via a gas phase processes, providing an additional source of carbon for organic material to be incorporated into planetesimals and cometesimals. This chemical processing mechanism requires warm grains (> 20 K), partially ionized gas, and sufficiently small <10 micron grains, i.e. a larger total grain surface area, such that freeze-out is efficient. Under these conditions static (non-turbulent) chemical models predict that a large fraction of the carbon nominally sequestered in CO can be the source of carbon for a wide variety of organics that are present as ice coatings on the surfaces of warm pre-planetesimal dust grains.

We present a comprehensive analysis of a broad band spectral line survey of the Orion Kleinmann-Low nebula (Orion KL), one of the most chemically rich regions in the Galaxy, using the HIFI instrument on board the Herschel Space Observatory. This survey spans a frequency range from 480 to 1907 GHz at a resolution of 1.1 MHz. These observations thus encompass the largest spectral coverage ever obtained toward this high-mass star-forming region in the sub-mm with high spectral resolution, and include frequencies $>$ 1 THz where the Earth's atmosphere prevents observations from the ground. In all, we detect emission from 39 molecules (79 isotopologues). Combining this dataset with ground based mm spectroscopy obtained with the IRAM 30 m telescope, we model the molecular emission from the mm to the far-IR using the XCLASS program which assumes local thermodynamic equilibrium (LTE). Several molecules are also modeled with the MADEX non-LTE code. Because of the wide frequency coverage, our models are constrained by transitions over an unprecedented range in excitation energy. A reduced $\chi^{2}$ analysis indicates that models for most species reproduce the observed emission well. In particular, most complex organics are well fit by LTE implying gas densities are high ($>$10$^6$ cm$^{-3}$) and excitation temperatures and column densities are well constrained. Molecular abundances are computed using H$_{2}$ column densities also derived from the HIFI survey. The distribution of rotation temperatures, $T_{\rm rot}$, for molecules detected toward the hot core is significantly wider than the compact ridge, plateau, and extended ridge $T_{\rm rot}$ distributions, indicating the hot core has the most complex thermal structure.

Context: The R CrA cloud hosts a handful of Class 0/I low-mass young stellar objects. The chemistry and physics at scales $>500$ AU in this cloud are dominated by the irradiation from the nearby Herbig Be star R CrA. The luminous large-scale emission makes it necessary to use high-resolution spectral imaging to study the chemistry and dynamics of the inner envelopes and discs of the protostars. Aims: We aim to better understand the structure of the inner regions of these protostars and, in particular, the interplay between the chemistry and the presence of discs. Methods: Using Atacama Large Millimeter/submillimeter Array (ALMA) high-resolution spectral imaging interferometry observations, we study the molecular line and dust continuum emission at submillimetre wavelengths. Results: We detect dust continuum emission from four circumstellar discs around Class 0/I objects within the R CrA cloud. Towards IRS7B we detect C$^{17}$O emission showing a rotation curve consistent with a Keplerian disc with a well-defined edge that gives a good estimate for the disc radius at 50 AU. We derive the central object mass to $2.3M_{\odot}$ and the disc mass to $0.024M_{\odot}$. The observations are also consistent with a model of material infalling under conservation of angular momentum; however, this model provides a worse fit to the data. We also report a likely detection of faint CH$_3$OH emission towards this point source, as well as more luminous CH$_3$OH emission in an outflow orthogonal to the major axis of the C$^{17}$O emission. Conclusions: The faint CH$_3$OH emission seen towards IRS7B can be explained by a flat density profile of the inner envelope caused by the disc with a radius $\lesssim50$ AU. We propose that the regions of the envelopes where complex organic molecules are present in Class 0/I young stellar objects can become quenched as the disc grows.

A sensitive broadband molecular line survey of the Sagittarius B2(N) star-forming region has been obtained with the HIFI instrument on the Herschel Space Observatory, offering the first high-spectral resolution look at this well-studied source in a wavelength region largely inaccessible from the ground (625-157 um). From the roughly 8,000 spectral features in the survey, a total of 72 isotopologues arising from 44 different molecules have been identified, ranging from light hydrides to complex organics, and arising from a variety of environments from cold and diffuse to hot and dense gas. We present an LTE model to the spectral signatures of each molecule, constraining the source sizes for hot core species with complementary SMA interferometric observations, and assuming that molecules with related functional group composition are cospatial. For each molecule, a single model is given to fit all of the emission and absorption features of that species across the entire 480-1910 GHz spectral range, accounting for multiple temperature and velocity components when needed to describe the spectrum. As with other HIFI surveys toward massive star forming regions, methanol is found to contribute more integrated line intensity to the spectrum than any other species. We discuss the molecular abundances derived for the hot core, where the local thermodynamic equilibrium approximation is generally found to describe the spectrum well, in comparison to abundances derived for the same molecules in the Orion KL region from a similar HIFI survey.

We describe the assignment of a previously unidentified interstellar absorption line to ArH$^+$ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H$_2$ fraction. The column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing interstellar medium (ISM) components with different H$_2$ abundances. Chemical models are constructed, taking UV radiation and cosmic ray ionization into account. Due to the detection of two isotopologues, $^{36}$ArH$^+$ and $^{38}$ArH$^+$, we are confident about the carrier assignment to ArH$^+$. NeH$^+$ is not detected with a limit of [NeH$^+$]/[ArH$^+$] $\le$ 0.1. The derived column densities agree well with the predictions of chemical models. ArH$^+$ is a unique tracer of gas with a fractional H$_2$ abundance of $10^{-4}- 10^{-3}$ and shows little correlation with H$_2$O$^+$, which traces gas with a fractional H$_2$ abundance of $\approx $0.1. A careful analysis of variations in the ArH$^+$, OH$^+$, H$_2$O$^+$ and HF column densities promises to be a faithful tracer of the distribution of the H$_2$ fractional abundance, providing unique information on a poorly known phase in the cycle of interstellar matter, its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD) simulations of the interstellar medium, and potentially could evolve into a tool to characterize the ISM. Paradoxically, the ArH$^+$ molecule is a better tracer of \new{almost} purely atomic hydrogen gas than H{\sc i} itself, since H{\sc i} can also be present in gas with a significant molecular content, but ArH$^+$ singles out gas that is $>99.9$\% atomic.

We present new observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ~ 500 K, rotational temperatures characterizing the population of the highly-excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone, but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ~ 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic ray fluxes, shocks and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly-excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation to the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process, not limited to the active environments associated with galactic nuclei.

The far-ultraviolet (FUV; 912--1700 \AA) radiation field from accreting central stars in Classical T Tauri systems influences the disk chemistry during the period of giant planet formation. Previous efforts to measure the true stellar+accretion-generated FUV luminosity (both hot gas emission lines and continua) have been complicated by a combination of low-sensitivity and/or low-spectral resolution and did not include the contribution from the bright Ly$\alpha$ emission line. In this work, we present a high-resolution spectroscopic study of the FUV radiation fields of 16 T Tauri stars whose dust disks display a range of evolutionary states. We include reconstructed Ly$\alpha$ line profiles and remove atomic and molecular disk emission (from H$_{2}$ and CO fluorescence) to provide robust measurements of both the FUV continuum and hot gas lines (e.g., Ly$\alpha$, \ion{N}{5}, \ion{C}{4}, \ion{He}{2}) for an appreciable sample of T Tauri stars for the first time. We find that the flux of the typical Classical T Tauri Star FUV radiation field at 1 AU from the central star is $\sim$ 10$^{7}$ times the average interstellar radiation field. The Ly$\alpha$ emission line contributes an average of 88% of the total FUV flux, with the FUV continuum accounting for an average of 8%. Both the FUV continuum and Ly$\alpha$ flux are strongly correlated with \ion{C}{4} flux, suggesting that accretion processes dominate the production of both of these components. On average, only $\sim$ 0.5% of the total FUV flux is emitted between the Lyman limit (912 \AA) and the H$_{2}$ (0 -- 0) absorption band at 1110 \AA. The total and component-level high-resolution radiation fields are made publicly available in machine-readable format.

Results from recent space missions, in particular Spitzer and Herschel, have lead to significant progress in our understanding of the formation and transport of water from clouds to disks, planetesimals, and planets. In this review, we provide the underpinnings for the basic molecular physics and chemistry of water and outline these advances in the context of water formation in space, its transport to a forming disk, its evolution in the disk, and finally the delivery to forming terrestrial worlds and accretion by gas giants. Throughout, we pay close attention to the disposition of water as vapor or solid and whether it might be subject to processing at any stage. The context of the water in the solar system and the isotopic ratios (D/H) in various bodies are discussed as grounding data point for this evolution. Additional advances include growing knowledge of the composition of atmospheres of extra-solar gas giants, which may be influenced by the variable phases of water in the protoplanetary disk. Further, the architecture of extra-solar systems leaves strong hints of dynamical interactions, which are important for the delivery of water and subsequent evolution of planetary systems. We conclude with an exploration of water on Earth and note that all of the processes and key parameters identified here should also hold for exoplanetary systems.

We report on the detection of an atomic jet associated with the protostellar source SMM1 (FIRS1) in Serpens. The jet is revealed in [FeII] and [NeII] line maps observed with Spitzer/IRS, and further confirmed in HiRes IRAC and MIPS images. It is traced very close to SMM1 and peaks at ~5 arcsec" from the source at a position angle of $\sim 125 degrees. In contrast, molecular hydrogen emission becomes prominent at distances > 5" from the protostar and extends at a position angle of 160 degrees. The morphological differences suggest that the atomic emission arises from a companion source, lying in the foreground of the envelope surrounding the embedded protostar SMM1. In addition the molecular and atomic Spitzer maps disentangle the large scale CO (3-2) emission observed in the region into two distinct bipolar outflows, giving further support to a proto-binary source setup. Analysis at the peaks of the [FeII] jet show that emission arises from warm and dense gas (T ~1000 K, n(electron) 10^5 - 10^6 cm^-3). The mass flux of the jet derived independently for the [FeII] and [NeII] lines is 10^7 M(sun)/yr, pointing to a more evolved Class~I/II protostar as the driving source. All existing evidence converge to the conclusion that SMM1 is a non-coeval proto-binary source.

Volatiles are compounds with low sublimation temperatures, and they make up most of the condensible mass in typical planet-forming environments. They consist of relatively small, often hydrogenated, molecules based on the abundant elements carbon, nitrogen and oxygen. Volatiles are central to the process of planet formation, forming the backbone of a rich chemistry that sets the initial conditions for the formation of planetary atmospheres, and act as a solid mass reservoir catalyzing the formation of planets and planetesimals. This growth has been driven by rapid advances in observations and models of protoplanetary disks, and by a deepening understanding of the cosmochemistry of the solar system. Indeed, it is only in the past few years that representative samples of molecules have been discovered in great abundance throughout protoplanetary disks - enough to begin building a complete budget for the most abundant elements after hydrogen and helium. The spatial distributions of key volatiles are being mapped, snow lines are directly seen and quantified, and distinct chemical regions within protoplanetary disks are being identified, characterized and modeled. Theoretical processes invoked to explain the solar system record are now being observationally constrained in protoplanetary disks, including transport of icy bodies and concentration of bulk condensibles. The balance between chemical reset - processing of inner disk material strong enough to destroy its memory of past chemistry, and inheritance - the chemically gentle accretion of pristine material from the interstellar medium in the outer disk, ultimately determines the final composition of pre-planetary matter. This chapter focuses on making the first steps toward understanding whether the planet formation processes that led to our solar system are universal.

We present Herschel/HIFI observations of the light hydride H$_{2}$S obtained from the full spectral scan of the Orion Kleinmann-Low nebula (Orion KL) taken as part of the HEXOS GT key program. In total, we observe 52, 24, and 8 unblended or slightly blended features from H$_{2}$$^{32}$S, H$_{2}$$^{34}$S, and H$_{2}$$^{33}$S, respectively. We only analyze emission from the so called hot core, but emission from the plateau, extended ridge, and/or compact ridge are also detected. Rotation diagrams for ortho and para H$_{2}$S follow straight lines given the uncertainties and yield T$_{\rm rot}$=141$\pm$12 K. This indicates H$_{2}$S is in LTE and is well characterized by a single kinetic temperature or an intense far-IR radiation field is redistributing the population to produce the observed trend. We argue the latter scenario is more probable and find that the most highly excited states (E$_{\rm up}$>1000 K) are likely populated primarily by radiation pumping. We derive an H$_{2}$$^{32}$S column density, N$_{\rm tot}$(H$_{2}$$^{32}$S)=9.5$\pm$1.9$\times$10$^{17}$ cm$^{-2}$, gas kinetic temperature, T$_{\rm kin}$=120$\pm^{13}_{10}$ K, and constrain the H$_{2}$ volume density, n$_{H2}$>9$\times$10$^{7}$ cm$^{-3}$, for the H$_{2}$S emitting gas. These results point to an H$_{2}$S origin in markedly dense, heavily embedded gas, possibly in close proximity to a hidden self-luminous source (or sources), which are conceivably responsible for Orion KL's high luminosity. We also derive an H$_{2}$S ortho/para ratio of 1.7$\pm$0.8 and set an upper limit for HDS/H2S of <4.9$\times$10$^{-3}$.

Low-mass protostars have been suggested to show highly variable accretion rates through-out their evolution. Such changes in accretion, and related heating of their ambient envelopes, may trigger significant chemical variations on different spatial scales and from source-to-source. We present images of emission from C17O, H13CO+, CH3OH, C34S and C2H toward the low-mass protostar IRAS 15398-3359 on 0.5" (75 AU diameter) scales with the Atacama Large Millimeter/submillimeter Array (ALMA) at 340 GHz. The resolved images show that the emission from H13CO+ is only present in a ring-like structure with a radius of about 1-1.5" (150-200 AU) whereas the CO and other high dipole moment molecules are centrally condensed toward the location of the central protostar. We propose that HCO+ is destroyed by water vapor present on small scales. The origin of this water vapor is likely an accretion burst during the last 100-1000 years increasing the luminosity of IRAS 15398-3359 by a factor of 100 above its current luminosity. Such a burst in luminosity can also explain the centrally condensed CH3OH and extended warm carbon-chain chemistry observed in this source and furthermore be reflected in the relative faintness of its compact continuum emission compared to other protostars.

We present simple analytic solutions for the ionization rate $\zeta_{\rm{SLR}}$ arising from the decay of short-lived radionuclides (SLRs) within protoplanetary disks. We solve the radiative transfer problem for the decay products within the disk, and thereby allow for the loss of radiation at low disk surface densities; energy loss becomes important outside $R\gtrsim30$ for typical disk masses $M_g=0.04$ M$_\odot$. Previous studies of chemistry/physics in these disks have neglected the impact of ionization by SLRs, and often consider only cosmic rays (CRs), because of the high CR-rate present in the ISM. However, recent work suggests that the flux of CRs present in the circumstellar environment could be substantially reduced by relatively modest stellar winds, resulting in severely modulated CR ionization rates, $\zeta_{\rm{CR}}$, equal to or substantially below that of SLRs ($\zeta_{\rm{SLR}}\lesssim10^{-18}$ s$^{-1}$). We compute the net ionizing particle fluxes and corresponding ionization rates as a function of position within the disk for a variety of disk models. The resulting expressions are especially simple for the case of vertically gaussian disks (frequently assumed in the literature). Finally, we provide a power-law fit to the ionization rate in the midplane as a function of gas disk surface density and time. Depending on location in the disk, the ionization rates by SLRs are typically in the range $\zeta_{\rm{SLR}}\sim(1-10)\times10^{-19}$ s$^{-1}$.

In this Letter we report the CO abundance relative to H2 derived toward the circumstellar disk of the T-Tauri star TW Hya from the HD (1-0) and C18O (2-1) emission lines. The HD (1-0) line was observed by the Herschel Space Observatory Photodetector Array Camera and Spectrometer whereas C18O (2-1) observations were carried out with the Submillimeter Array at a spatial resolution of 2.8" x 1.9" (corresponding to 142 x 97 AU). In the disk's warm molecular layer (T>20 K) we measure a disk-averaged gas-phase CO abundance relative to H2 of $\chi{\rm(CO)}=(0.1-3)x10^{-5}$, substantially lower than the canonical value of $\chi{\rm(CO)}=10^{-4}$. We infer that the best explanation of this low $\chi$(CO) is the chemical destruction of CO followed by rapid formation of carbon chains, or perhaps CO2, that can subsequently freeze-out, resulting in the bulk mass of carbon locked up in ice grain mantles and oxygen in water. As a consequence of this likely time-dependent carbon sink mechanism, CO may be an unreliable tracer of H2 gas mass.

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

In this review we explore aspects of the field of astrobiology from an astronomical viewpoint. We therefore focus on the origin of life in the context of planetary formation, with additional emphasis on tracing the most abundant volatile elements, C, H, O, and N that are used by life on Earth. We first explore the history of life on our planet and outline the current state of our knowledge regarding the delivery of the C, H, O, N elements to the Earth. We then discuss how astronomers track the gaseous and solid molecular carriers of these volatiles throughout the process of star and planet formation. It is now clear that the early stages of star formation fosters the creation of water and simple organic molecules with enrichments of heavy isotopes. These molecules are found as ice coatings on the solid materials that represent microscopic beginnings of terrestrial worlds. Based on the meteoritic and cometary record, the process of planet formation, and the local environment, lead to additional increases in organic complexity. The astronomical connections towards this stage are only now being directly made. Although the exact details are uncertain, it is likely that the birth process of star and planets likely leads to terrestrial worlds being born with abundant water and organics on the surface.

We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.

We present a comprehensive study of the deuterated molecules detected in the fullband HIFI survey of the Orion KL region. Ammonia, formaldehyde, and methanol and their singly deuterated isotopologues are each detected through numerous transitions in this survey with a wide range in optical depths and excitation conditions. In conjunction with a recent study of the abundance of HDO and H$_2$O in Orion KL, this study yields the best constraints on deuterium fractionation in an interstellar molecular cloud to date. As previous studies have found, both the Hot Core and Compact Ridge regions within Orion KL contain significant abundances of deuterated molecules, suggesting an origin in cold grain mantles. In the Hot Core, we find that ammonia is roughly a factor of 2 more fractionated than water. In the Compact Ridge, meanwhile, we find similar deuterium fractionation in water, formaldehyde, and methanol, with D/H ratios of (2---8) $\times$ $10^{-3}$. The [CH$_2$DOH]/[CH$_3$OD] ratio in the Compact Ridge is found to be $1.2 \pm 0.3$. The Hot Core generally has lower deuterium fractionation than the Compact Ridge, suggesting a slightly warmer origin, or a greater contribution from warm gas phase chemistry.

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

Planets form in the disks around young stars. Their formation efficiency and composition are intimately linked to the protoplanetary disk locations of "snow lines" of abundant volatiles. We present chemical imaging of the CO snow line in the disk around TW Hya, an analog of the solar nebula, using high spatial and spectral resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of N2H+, a reactive ion present in large abundance only where CO is frozen out. The N2H+ emission is distributed in a large ring, with an inner radius that matches CO snow line model predictions. The extracted CO snow line radius of ~ 30 AU helps to assess models of the formation dynamics of the Solar System, when combined with measurements of the bulk composition of planets and comets.

We present Herschel observations of water isotopologues in the atmosphere of the Jupiter-family comet 45P/Honda-Mrkos-Pajdusakova. No HDO emission is detected, with a 3 sigma upper limit of 2.0 10-4 for the D/H ratio. This value is consistent with the earlier Herschel measurement in the Jupiter-family comet 103P/Hartley 2. The canonical value of 3 10-4 measured pre-Herschel in a sample of Oort-cloud comets can be excluded at a 4.5 sigma level. The observations presented here further confirm that a diversity of D/H ratios exists in the comet population and emphasize the need for additional measurements with future ground-based facilities, such as CCAT, in the post-Herschel era.

(Abridged) Cosmic rays (CRs) are thought to provide an important source of ionization in the outermost and densest regions of protoplanetary disks; however, it is unknown to what degree they are physically present. As is observed in the Solar System, stellar winds can inhibit the propagation of cosmic rays within the circumstellar environment and subsequently into the disk. In this work, we explore the hitherto neglected effects of cosmic ray modulation by both stellar winds and magnetic field structures and study how these processes act to reduce disk ionization rates. We construct a two-dimensional protoplanetary disk model of a T-Tauri star system, focusing on ionization from stellar and interstellar FUV, stellar X-ray photons, and cosmic rays. We show that stellar winds can power a Heliosphere-like analogue, i.e., a "T-Tauriosphere," diminishing cosmic ray ionization rates by several orders of magnitude at low to moderate CR energies (E_CR<1 GeV). We explore models of both the observed solar wind cosmic ray modulation and a highly simplified estimate for "elevated" cosmic ray modulation as would be expected from a young T-Tauri star. In the former (solar analogue) case, we estimate the ionization rate from galactic cosmic rays to be zeta_CR~(0.23-1.4) x 1e-18 s^-1. This range of values, which we consider to be the maximum CR ionization rate for the disk, is more than an order of magnitude lower than what is generally assumed in current models for disk chemistry and physics. In the latter case, i.e., for a "T-Tauriosphere," the ionization rate by cosmic rays is zeta_CR<1e-20 s^-1, which is 1000 times smaller than the interstellar value. Indeed, if winds are as efficient at cosmic ray modulation as predicted here, short-lived radionuclides (now extinct) would have provided the major source of ionization (zeta_RN~7.3e-19 s^-1) in the planet-forming zone of the young Solar Nebula.

Using a broadband, high spectral resolution survey toward Orion KL acquired with Herschel/HIFI as part of the HEXOS key program, we derive the abundances of H$_2$O and HDO in the different spatial/velocity components associated with this massive star-forming region: the Hot Core, Compact Ridge, and Plateau. A total of 20 transitions of H$_2$$^{18}$O, 14 of H$_2$$^{17}$O, 37 of HD$^{16}$O, 6 of HD$^{18}$O, and 6 of D$_2$O are used in the analysis, spanning from ground state transitions to over 1200 K in upper-state energy. Low-excitation lines are detected in multiple components, but the highest-excitation lines ($E_u >$ 500 K) are well modeled as emitting from a small ($\sim 2"$) clump with a high abundance of H$_2$O ($\chi = 6.5 \times 10^{-4}$ relative to H$_2$) and a HDO/H$_2$O ratio of 0.003. Using high spatial resolution ($1.5" \times 1.1"$) images of two transitions of HDO measured by ALMA as part of its science verification phase, we identify this component as located near, but not directly coincident with, known continuum sources in the Hot Core region. Significant HDO/H$_2$O fractionation is also seen in the Compact Ridge and Plateau components. The outflowing gas, observed with both emission and absorption components, has a lower HDO/H$_2$O ratio than the compact components in Orion KL, which we propose could be due to modification by gas-phase shock chemistry.

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

We explore a set of single-point chemical models to study the fundamental chemical aspects of episodic accretion in low-mass embedded protostars. Our goal is twofold: (1) to understand how the repeated heating and cooling of the envelope affects the abundances of CO and related species; and (2) to identify chemical tracers that can be used as a novel probe of the timescales and other physical aspects of episodic accretion. We develop a set of single-point models that serve as a general prescription for how the chemical composition of a protostellar envelope is altered by episodic accretion. The main effect of each accretion burst is to drive CO ice off the grains in part of the envelope. The duration of the subsequent quiescent stage (before the next burst hits) is similar to or shorter than the freeze-out timescale of CO, allowing the chemical effects of a burst to linger long after the burst has ended. We predict that the resulting excess of gas-phase CO can be observed with single-dish or interferometer facilities as evidence of an accretion burst in the past 10^3 - 10^4 yr.

For Classical T Tauri Stars (CTTSs), the resonance lines of N V, Si IV, and C IV, as well as the He II 1640 A line, act as diagnostics of the accretion process. Here we assemble a large high-resolution dataset of these lines in CTTSs and Weak T Tauri Stars (WTTSs). We present data for 35 stars: one Herbig Ae star, 28 CTTSs, and 6 WTTSs. We decompose the C IV and He II lines into broad and narrow Gaussian components (BC & NC). The most common (50 %) C IV line morphology in CTTSs is that of a low-velocity NC together with a redshifted BC. The velocity centroids of the BCs and NCs are such that V_BC > 4 * V_NC, consistent with the predictions of the accretion shock model, in at most 12 out of 22 CTTSs. We do not find evidence of the post-shock becoming buried in the stellar photosphere due to the pressure of the accretion flow. The He II CTTSs lines are generally symmetric and narrow, less redshifted than the CTTSs C IV lines, by ~10 km/sec. The flux in the BC of the He II line is small compared to that of the C IV line, consistent with models of the pre-shock column emission. The observations are consistent with the presence of multiple accretion columns with different densities or with accretion models that predict a slow-moving, low-density region in the periphery of the accretion column. For HN Tau A and RW Aur A, most of the C IV line is blueshifted suggesting that the C IV emission is produced by shocks within outflow jets. In our sample, the Herbig Ae star DX Cha is the only object for which we find a P-Cygni profile in the C IV line, which argues for the presence of a hot (10^5 K) wind. For the overall sample, the Si IV and N V line luminosities are correlated with the C IV line luminosities, although the relationship between Si IV and C IV shows large scatter about a linear relationship and suggests that TW Hya, V4046 Sgr, AA Tau, DF Tau, GM Aur, and V1190 Sco are silicon-poor.

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

From the masses of planets orbiting our Sun, and relative elemental abundances, it is estimated that at birth our Solar System required a minimum disk mass of ~0.01 solar masses within ~100 AU of the star. The main constituent, gaseous molecular hydrogen, does not emit from the disk mass reservoir, so the most common measure of the disk mass is dust thermal emission and lines of gaseous carbon monoxide. Carbon monoxide emission generally probes the disk surface, while the conversion from dust emission to gas mass requires knowledge of the grain properties and gas-to-dust mass ratio, which likely differ from their interstellar values. Thus, mass estimates vary by orders of magnitude, as exemplified by the relatively old (3--10 Myr) star TW Hya, with estimates ranging from 0.0005 to 0.06 solar masses. Here we report the detection the fundamental rotational transition of hydrogen deuteride, HD, toward TW Hya. HD is a good tracer of disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The HD detection, combined with existing observations and detailed models, implies a disk mass >0.05 solar masses, enough to form a planetary system like our own.

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

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

The abundances of interstellar CH+ and SH+ are not well understood as their most likely formation channels are highly endothermic. Using data from Herschel, we study the formation of CH+ and SH+ in a typical high UV-illumination photon-dominated region (PDR), the Orion Bar. Herschel/HIFI provides velocity-resolved data of CH+ 1-0 and 2-1 and three hyperfine transitions of SH+. Herschel/PACS provides information on the excitation and spatial distribution of CH+ (up to J=6-5). The widths of the CH+ 2-1 and 1-0 transitions are of ~5 km/s, significantly broader than the typical width of dense gas tracers in the Orion Bar (2-3 km/s) and are comparable to the width of tracers of the interclump medium such as C+ and HF. The detected SH+ transitions are narrower compared to CH+ and have line widths of 3 km/s, indicating that SH+ emission mainly originates in denser condensations. Non-LTE radiative transfer models show that electron collisions affect the excitation of CH+ and SH+, and that reactive collisions need to be taken into account to calculate the excitation of CH+. Comparison to PDR models shows that CH+ and SH+ are tracers of the warm surface region (AV<1.5) of the PDR with temperatures between 500-1000 K. We have also detected the 5-4 transition of CF+ (FWHM=1.9 km/s) with an intensity that is consistent with previous observations of lower-J CF+ transitions toward the Orion Bar. A comparison to PDR models indicate that the internal vibrational energy of H2 can explain the formation of CH+ for typical physical conditions in the Orion Bar near the ionization front. H2 vibrational excitation is the most likely explanation of SH+ formation as well. The abundance ratios of CH+ and SH+ trace the destruction paths of these ions, and through that, indirectly, the ratios of H, H2 and electron abundances as a function of depth into the cloud.

Recently, we introduced detailed isotopic chemistry into the KOSMA-tau model for photon-dominated regions (PDRs) to give theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative transfer computations for specific geometries, we estimated the possible intensity ratio of the [CII]/[13CII] lines. Here, we compare these predictions with new observations. We performed Herschel/HIFI observations of the [CII] 158micron line in a number of PDRs. In all sources we observed at least two hyperfine components of the [13CII] transition allowing to determine the [CII]/[13CII] intensity ratio, after some revision of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-tau model, we identify cases dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. An observable enhancement of the [CII]/[13CII] intensity ratio due to chemical fractionation depends mostly on geometry and velocity structure, and less on the gas density and radiation field. In our observations the [CII]/[13CII] ratio for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact HII region Mon R2, the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC7023. The mapping of the [13CII] lines in the Orion Bar allows to derive a C+ column density map confirming the temperature stratification of the C+ layer, in agreement with the chemical stratification of the Bar. The C+ column densities for all sources show that at the position of the [CII] peak emission, a dominant fraction of the gas-phase carbon is in the form of C+.

High-energy radiation from T Tauri stars (TTS) influences the amount and longevity of gas in disks, thereby playing a crucial role in the creation of gas giant planets. Here we probe the high-energy ionizing radiation from TTS using high-resolution mid-infrared (MIR) Spitzer IRS Neon forbidden line detections in a sample of disks from IC 348, NGC 2068, and Chamaeleon. We report three new detections of [Ne III] from CS Cha, SZ Cha, and T 54, doubling the known number of [Ne III] detections from TTS. Using [Ne III]-to-[Ne II] ratios in conjunction with X-ray emission measurements, we probe high-energy radiation from TTS. The majority of previously inferred [Ne III]/[Ne II] ratios based on [Ne III] line upper limits are significantly less than 1, pointing to the dominance of either X-ray radiation or soft Extreme-Ultraviolet (EUV) radiation in producing these lines. Here we report the first observational evidence for hard EUV dominated Ne forbidden line production in a T Tauri disk: [Ne III]/[Ne II]~1 in SZ Cha. Our results provide a unique insight into the EUV emission from TTS, by suggesting that EUV radiation may dominate the creation of Ne forbidden lines, albeit in a minority of cases.

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

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

Far-ultraviolet (FUV) radiation plays an important role in determining chemical abundances in protoplanetary disks. HI Lyman alpha is suspected to be the dominant component of the FUV emission from Classical T Tauri Stars (CTTSs), but is difficult to measure directly due to circumstellar and interstellar HI absorption. To better characterize the intrinsic Lyman alpha radiation, we present FUV spectra of 14 CTTSs taken with the Hubble Space Telescope COS and STIS instruments. H2 fluorescence, commonly seen in the spectra of CTTSs, is excited by Lyman alpha photons, providing an indirect measure of the Lyman alpha flux incident upon the warm disk surface. We use observed H2 progression fluxes to reconstruct the CTTS Lyman alpha profiles. The Lyman alpha flux correlates with total measured FUV flux, in agreement with an accretion-related source of FUV emission. With a geometry-independent analysis, we confirm that in accreting T Tauri systems Lyman alpha radiation dominates the FUV flux (~1150 - 1700 Angstroms). In the systems surveyed this one line comprises 70 - 90 % of the total FUV flux.

Few scientific discoveries have captured the public imagination like the explosion of exoplanetary science during the past two decades. This work has fundamentally changed our picture of Earth's place in the Universe and led NASA to make significant investments towards understanding the demographics of exoplanetary systems and the conditions that lead to their formation. The story of the formation and evolution of exoplanetary systems is essentially the story of the circumstellar gas and dust that are initially present in the protostellar environment; in order to understand the variety of planetary systems observed, we need to understand the life cycle of circumstellar gas from its initial conditions in protoplanetary disks to its endpoint as planets and their atmospheres. In this white paper response to NASA's Request for Information "Science Objectives and Requirements for the Next NASA UV/Visible Astrophysics Mission Concepts (NNH12ZDA008L)", we describe scientific programs that would use the unique capabilities of a future NASA ultraviolet (UV)/visible space observatory to make order-of-magnitude advances in our understanding of the life cycle of circumstellar gas.

The Heterodyne Instrument for the Far Infrared (HIFI) aboard the Herschel Space Observatory has acquired high-resolution broadband molecular spectra of star-forming regions in a wavelength range that is mostly inaccessible from ground-based astronomical observatories. These spectral surveys provide new insight into the chemical composition and physical properties of molecular clouds. In this manuscript, we present initial results from the HIFI spectral survey of the Sagittarius B2(N) molecular cloud, which contains spectral features assigned to at least 40 different molecules in a range of physical environments. While extensive line blending is observed due to the chemical complexity of this region, reliable molecular line identifications can be made, down to the noise floor, due to the large number of transitions detected for each species in the 1.2 THz survey bandwidth. This allows for the extraction of new weakly emitting species from the line forest. These HIFI surveys will be an invaluable archival resource for future investigations into interstellar chemistry.

The formation timescale and final architecture of exoplanetary systems are closely related to the properties of the molecular disks from which they form. Observations of the spatial distribution and lifetime of the molecular gas at planet-forming radii (r < 10 AU) are important for understanding the formation and evolution of exoplanetary systems. Towards this end, we present the largest spectrally resolved survey of H2 emission around low-mass pre-main sequence stars compiled to date. We use a combination of new and archival far-ultraviolet spectra from the COS and STIS instruments on the Hubble Space Telescope to sample 34 T Tauri stars (27 actively accreting CTTSs and 7 non-accreting WTTSs) with ages ranging from roughly 1-10 Myr. We observe fluorescent H2 emission, excited by LyA photons, in 100 of the accreting sources, including all of the transitional disks in our sample (CS Cha, DM Tau, GM Aur, UX Tau A, LkCa15, HD 135344B and TW Hya). The spatial distribution of the emitting gas is inferred from spectrally resolved H2 line profiles. Some of the emitting gas is produced in outflowing material, but the majority of H2 emission appears to originate in a rotating disk. For the disk-dominated targets, the H2 emission originates predominately at r < 3 AU. The emission line-widths and inner molecular radii are found to be roughly consistent with those measured from mid-IR CO spectra.

We present a far-ultraviolet (FUV) spectral atlas consisting of spectra of 91 pre-main sequence stars. Most stars in this sample were observed with the Space Telescope Imaging Spectrograph (STIS) and Advanced Camera for Surveys (ACS) on the \emph{Hubble Space Telescope} (\emph{HST}). We find strong correlations among the \ion{O}{1} $\lambda$1304 triplet, %\ion{C}{2} $\lambda$1335, the \ion{Si}{4} $\lambda\lambda$1394/1403 doublet, the \ion{C}{4} $\lambda$1549 doublet, and the \ion{He}{2} $\lambda$1640 line luminosities. For classical T Tauri stars (CTTSs), we also find strong correlations between these lines and the accretion luminosity, suggesting that these lines form in processes related to accretion. These FUV line fluxes and X-ray luminosity correlate loosely with large scatters. The FUV emission also correlates well with H$\alpha$, H$\beta$, and \ion{Ca}{2} K line luminosities. These correlations between FUV and optical diagostics can be used to obtain rough estimates of FUV line fluxes from optical observations. Molecular hydrogen (H$_{2}$) emission is generally present in the spectra of actively accreting CTTSs but not the weak-lined T Tauri stars (WTTSs) that are not accreting. The presence of H$_2$ emission in the spectrum of HD 98800 N suggests that the disk should be classified as actively accreting rather than a debris disk. The spectra in the atlas are available at http://archive.stsci.edu/prepds/ttauriatlas.

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

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]

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

Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and here, we investigate the kinematic structure and energy content. We find that IRDCs overall display organized velocity fields, with only localized disruptions due to embedded star formation. The local effects seen in NH3 emission are not high velocity outflows but rather moderate (few km/s) increases in the line width that exhibit maxima near or coincident with the mid-infrared emission tracing protostars. These line width enhancements could be the result of infall or (hidden in NH3 emission) outflow. Not only is the kinetic energy content insufficient to support the IRDCs against collapse, but also the spatial energy distribution is inconsistent with a scenario of turbulent cloud support. We conclude that the velocity signatures of the IRDCs in our sample are due to active collapse and fragmentation, in some cases augmented by local feedback from stars.

The clumpy density structure of photon-dominated regions is well established, but the physical properties of the clumps and of the surrounding interclump medium are only approximately known. The aim of this paper is to constrain the physical and chemical conditions in the Orion Bar, a prototypical nearby photon-dominated region. We present observations of the HF J=1-0 line, which appears in emission toward the Orion Bar, and compare the brightness of the line to non-LTE radiative transfer calculations. The large width of the HF line suggests an origin of the emission in the interclump gas, but collisional excitation by H2 in the interclump gas underpredicts the observed line intensity by factors of 3-5. In contrast, an origin of the line in the dense clumps requires a density of ~10^9 cm^-3, 10-100 times higher than previous estimates, which is unlikely. However, electron impact excitation reproduces our observations for T = 100 K and n(e) = 10 cm^-3, as expected for the interclump gas. We conclude that HF emission is a signpost of molecular gas with a high electron density. Similar conditions may apply to active galactic nuclei where HF also appears in emission.

We explore the chemical structure of a disk that contains a large central gap of R ~ 45 AU, as is commonly seen in transitional disk systems. In our chemical model of a disk with a cleared inner void, the midplane becomes revealed to the central star so that it is directly irradiated. The midplane material at the truncation radius is permissive to reprocessed optical heating radiation, but opaque to the photo-dissociating ultraviolet, creating an environment abundant in gas-phase molecules. Thus the disk midplane, which would otherwise for a full disk be dominated by near complete heavy element freeze-out, should become observable in molecular emission. If this prediction is correct this has exciting prospects for observations with the Atacama Large Millimeter/Submillimeter Array (ALMA), as the inner transition region should thus be readily detected and resolved, especially using high-J rotational transitions excited in the high density midplane gas. Therefore such observations will potentially provide us with a direct probe of the physics and chemistry at this actively evolving interface.

Young stars surrounded by disks with very low mass accretion rates are likely in the final stages of inner disk evolution and therefore particularly interesting to study. We present ultraviolet (UV) observations of the ~5-9 Myr old stars RECX-1 and RECX-11, obtained with the Cosmic Origins Spectrograph (COS) and Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST), as well as optical and near infrared spectroscopic observations. The two stars have similar levels of near UV emission, although spectroscopic evidence indicates that RECX-11 is accreting and RECX-1 is not. The line profiles of Halpha and He I 10830 in RECX-11 show both broad and narrow redshifted absorption components that vary with time, revealing the complexity of the accretion flows. We show that accretion indicators commonly used to measure mass accretion rates, e.g. U band excess luminosity or the Ca II triplet line luminosity are unreliable for low accretors, at least in the middle K spectral range. Using RECX-1 as a template for the intrinsic level of photospheric and chromospheric emission, we determine an upper limit of 3 x 10^-10 solar masses/ year for RECX-11. At this low accretion rate, recent photoevaporation models predict that an inner hole should have developed in the disk. However, the spectral energy distribution of RECX-11 shows fluxes comparable to the median of Taurus in the near infrared, indicating that substantial dust remains. Fluorescent H_2 emission lines formed in the innermost disk are observed in RECX-11, showing that gas is present in the inner disk, along with the dust.

The C/O ratio is predicted to regulate the atmospheric chemistry in hot Jupiters. Recent observations suggest that some exo-planets, e.g. Wasp 12- b, have atmospheric C/O ratios substantially different from the solar value of 0.54. In this paper we present a mechanism that can produce such atmospheric deviations from the stellar C/O ratio. In protoplanetary disks, different snowlines of oxygen- and carbon-rich ices, especially water and carbon monoxide, will result in systematic variations in the C/O ratio both in the gas and in the condensed phase. In particular, between the H2O and CO snowlines most oxygen is present in icy grains - the building blocks of planetary cores in the core accretion model - while most carbon remains in the gas-phase. This region is coincidental with the giant-planet forming zone for a range of observed protoplanetary disks. Based on standard core accretion models of planet formation, gas giants that sweep up most of their atmospheres from disk gas outside of the water snowline will have C/O?1, while atmospheres significantly contaminated by evaporating planetesimals will have stellar or sub-stellar C/O when formed at the same disk radius. The overall metallicity will also depend on the atmosphere formation mechanism, and exoplanetary atmospheric compositions may therefore provide constraints on where and how a specific planet formed.

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

In this paper we discuss the astronomical search for water vapor in order to understand the disposition of water in all its phases throughout the process of star and planet formation. Our ability to detect and study water vapor has recently received a tremendous boost with the successful launch and operations of the Herschel Space Observatory. Herschel spectroscopic detections of numerous transitions in a variety of astronomical objects, along with previous work by other space-based observatories, will be threaded throughout this contribution. In particular, we present observations of water tracing the earliest stage of star birth where it is predominantly frozen as ice. When a star is born the local energy release by radiation liberates ices in its surrounding envelope and powers energetic outflows that appear to be water factories. In these regions water plays an important role in the gas physics. Finally, we end with an exploration of water in planet forming disks surrounding young stars. The availability of accurate molecular data (frequencies, collisional rate coefficients, and chemical reaction rates) are crucial to analyze the observations at each of these steps.

Herschel/HIFI observations have revealed the presence of widespread absorption by hydrogen fluoride (HF) J = 1-0 rotational transition, toward a number of Galactic sources. We present observations of HF J = 1-0 toward the high-mass star-forming region Sagittarius B2(M). The spectrum obtained shows a complex pattern of absorption, with numerous features covering a wide range of local standard of rest velocities (-130 to 100 km -1). An analysis of this absorption yields HF abundances relative to H2 of ~1.3 {\times}10-8, in most velocity intervals. This result is in good agreement with estimates from chemical models, which predict that HF should be the main reservoir of gas-phase fluorine under a wide variety of interstellar conditions. Interestingly, we also find velocity intervals in which the HF spectrum shows strong absorption features that are not present, or are very weak, in spectra of other molecules, such as 13CO (1-0) and CS (2-1). HF absorption reveals components of diffuse clouds with small extinction that can be studied for the first time. Another interesting observation is that water is significantly more abundant than hydrogen fluoride over a wide range of velocities toward Sagittarius B2(M), in contrast to the remarkably constant H2O/HF abundance ratio with average value close to unity measured toward other Galactic sources.

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.

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

We present Very Large Array observations of NH3 (1,1) and (2,2) in a sample of six infrared-dark clouds (IRDCs) with distances from 2 to 5 kpc. We find that ammonia serves as an excellent tracer of dense gas in IRDCs, showing no evidence of depletion, and the average abundance in these clouds is 8.1 x 10^-7 relative to H2. Our sample consists of four IRDCs with 24 micron embedded protostars and two that appear starless. We calculate the kinetic temperature of the gas in IRDCs and find no significant difference between starless and star-forming IRDCs. We find that the bulk of the gas is between 8 and 13 K, indicating that any embedded or nearby stars or clusters do not affect the gas temperature dramatically. Though IRDCs have temperatures and volume densities on par with local star formation regions of lower mass, they consist of much more mass which induces very high internal pressures. In order for IRDCs to survive as coherent structures, the internal pressure must be balanced by a confining pressure provided by the high concentration of molecular clouds in the spiral arm in which they reside. The high molecular concentration and pressure is roughly consistent with gas dynamics of a bar galaxy.

We provide simple polynomial fits to the X-ray photoelectric cross-sections (0.03 < E < 10keV) for mixtures of gas and dust found in protoplanetary disks. Using the solar elemental abundances of Asplund et al. (2009) we treat the gas and dust components separately, facilitating the further exploration evolutionary processes such as grain settling and gain growth. We find that blanketing due to advanced grain-growth (a_max > 1 micron) can reduce the X-ray opacity of dust appreciably at E_X ~ 1keV, coincident with the peak of typical T Tauri X-ray spectra. However, the reduction of dust opacity by dust settling, which is known to occur in protoplanetary disks, is probably a more significant effect. The absorption of 1-10keV X-rays is dominated by gas opacity once the dust abundance has been reduced to about 1% of its diffuse interstellar value. The gas disk establishes a floor to the opacity at which point X-ray transport becomes insensitive to further dust evolution. Our choice of fitting function follows that of Morrison & McCammon (1983), providing a degree of backward-compatibility.

We study the role resonant scattering plays in the transport of Lyman-alpha photons in accreting protoplanetary disk systems subject to varying degrees of dust settling. While the intrinsic stellar FUV spectrum of accreting T Tauri systems may already be dominated by a strong, broad Lyman-alpha line (~80% of the FUV luminosity), we find that resonant scattering further enhances the Lyman-alpha density in the deep molecular layers of the disk. Lyman-alpha is scattered downwards efficiently by the photodissociated atomic hydrogen layer that exists above the molecular disk. In contrast, FUV-continuum photons pass unimpeded through the photodissociation layer, and (forward-)scatter inefficiently off dust grains. Using detailed, adaptive grid Monte Carlo radiative transfer simulations we show that the resulting Lyman-alpha/FUV-continuum photon density ratio is strongly stratified; FUV-continuum dominated in the photodissociation layer, and Lyman-alpha dominated field in the molecular disk. The enhancement is greatest in the interior of the disk (r~1AU) but is also observed in the outer disk (r~100AU). The majority of the total disk mass is shown to be increasingly Lyman-alpha dominated as dust settles towards the midplane.

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.

As part of a far-infrared (FIR) spectral scan with Herschel/PACS, we present the first detection of the hydroxyl radical (OH) towards the Orion Bar photodissociation region (PDR). Five OH rotational Lambda-doublets involving energy levels out to E_u/k~511 K have been detected (at ~65, ~79, ~84, ~119 and ~163um). The total intensity of the OH lines is I(OH)~5x10^-4 erg s^-1 cm^-2 sr^-1. The observed emission of rotationally excited OH lines is extended and correlates well with the high-J CO and CH^+ J=3-2 line emission (but apparently not with water vapour), pointing towards a common origin. Nonlocal, non-LTE radiative transfer models including excitation by the ambient FIR radiation field suggest that OH arises in a small filling factor component of warm (Tk~160-220 K) and dense (n_H~10^{6-7} cm^-3) gas with source-averaged OH column densities of ~10^15 cm^-2. High density and temperature photochemical models predict such enhanced OH columns at low depths (A_V<1) and small spatial scales (~10^15 cm), where OH formation is driven by gas-phase endothermic reactions of atomic oxygen with molecular hydrogen. We interpret the extended OH emission as coming from unresolved structures exposed to far-ultraviolet (FUV) radiation near the Bar edge (photoevaporating clumps or filaments) and not from the lower density "interclump" medium. Photodissociation leads to OH/H2O abundance ratios (>1) much higher than those expected in equally warm regions without enhanced FUV radiation fields.

This is the second in a series of papers based on data from DISCS, a Submillimeter Array observing program aimed at spatially and spectrally resolving the chemical composition of 12 protoplanetary disks. We present data on six Southern sky sources - IM Lup, SAO 206462 (HD 135344b), HD 142527, AS 209, AS 205 and V4046 Sgr - which complement the six sources in the Taurus star forming region reported previously. CO 2-1 and HCO+ 3-2 emission are detected and resolved in all disks and show velocity patterns consistent with Keplerian rotation. Where detected, the emission from DCO+ 3-2, N2H+ 3-2, H2CO 3-2 and 4-3,HCN 3-2 and CN 2-1 are also generally spatially resolved. The detection rates are highest toward the M and K stars, while the F star SAO 206462 has only weak CN and HCN emission, and H2CO alone is detected toward HD 142527. These findings together with the statistics from the previous Taurus disks, support the hypothesis that high detection rates of many small molecules depend on the presence of a cold and protected disk midplane, which is less common around F and A stars compared to M and K stars. Disk-averaged variations in the proposed radiation tracer CN/HCN are found to be small, despite two orders of magnitude range of spectral types and accretion rates. In contrast, the resolved images suggest that the CN/HCN emission ratio varies with disk radius in at least two of the systems. There are no clear observational differences in the disk chemistry between the classical/full T Tauri disks and transitional disks. Furthermore, the observed line emission does not depend on measured accretion luminosities or the number of infrared lines detected, which suggests that the chemistry outside of 100 AU is not coupled to the physical processes that drive the chemistry in the innermost few AU.

We exploit the high sensitivity and moderate spectral resolution of the $HST$-Cosmic Origins Spectrograph to detect far-ultraviolet spectral features of carbon monoxide (CO) present in the inner regions of protoplanetary disks for the first time. We present spectra of the classical T Tauri stars HN Tau, RECX-11, and V4046 Sgr, representative of a range of CO radiative processes. HN Tau shows CO bands in absorption against the accretion continuum. We measure a CO column density and rotational excitation temperature of N(CO) = 2 +/- 1 $\times$ 10$^{17}$ cm$^{-2}$ and T_rot(CO) 500 +/- 200 K for the absorbing gas. We also detect CO A-X band emission in RECX-11 and V4046 Sgr, excited by ultraviolet line photons, predominantly HI LyA. All three objects show emission from CO bands at $\lambda$ $>$ 1560 \AA, which may be excited by a combination of UV photons and collisions with non-thermal electrons. In previous observations these emission processes were not accounted for due to blending with emission from the accretion shock, collisionally excited H$_{2}$, and photo-excited H2; all of which appeared as a "continuum" whose components could not be separated. The CO emission spectrum is strongly dependent upon the shape of the incident stellar LyA emission profile. We find CO parameters in the range: N(CO) 10$^{18-19}$ cm$^{-2}$, T_{rot}(CO) > 300 K for the LyA-pumped emission. We combine these results with recent work on photo- and collisionally-excited H$_{2}$ emission, concluding that the observations of ultraviolet-emitting CO and H2 are consistent with a common spatial origin. We suggest that the CO/H2 ratio in the inner disk is ~1, a transition between the much lower interstellar value and the higher value observed in solar system comets today, a result that will require future observational and theoretical study to confirm.

We present new continuum and line observations, along with modelling, of the faint (6-8) Myr old T Tauri star ET Cha belonging to the eta Chamaeleontis cluster. We have acquired HERSCHEL/PACS photometric fluxes at 70 mic and 160 mic, as well as a detection of the [OI] 63 mic fine-structure line in emission, and derived upper limits for some other far-IR OI, CII, CO and o-H2O lines. The HERSCHEL data is complemented by new ANDICAM B-K photometry, new HST/COS and HST/STIS UV-observations, a non-detection of CO J=3-2 with APEX, re-analysis of a UCLES high-resolution optical spectrum showing forbidden emission lines like [OI] 6300A, [SII] 6731A and 6716A, and [NII] 6583A, and a compilation of existing broad-band photometric data. We used the thermo-chemical disk code ProDiMo and the Monte-Carlo radiative transfer code MCFOST to model the protoplanetary disk around ET Cha. Based on these models we can determine the disk dust mass Mdust = (2.E-8 - 5.E-8) Msun, whereas the total disk gas mass is found to be only little constrained, Mgas = (5.E-5 - 3.E-3) Msun. In the models, the disk extends from 0.022 AU (just outside of the co-rotation radius) to only about 10 AU. Larger disks are found to be inconsistent with the CO J=3-2 non-detection. The low velocity component of the [OI] 6300A emission line is consistent with being emitted from the inner disk. The model can also reproduce the line flux of H2 v=1-0 S(1) at 2.122 mic. An additional high-velocity component of the [OI] 6300A emission line, however, points to the existence of an additional jet/outflow of low velocity (40 - 65) km/s with mass loss rate ~1.E-9 Msun/yr. In relation to our low estimations of the disk mass, such a mass loss rate suggests a disk lifetime of only ~(0.05 - 3) Myr, substantially shorter than the cluster age. The evolutionary state of this unusual protoplanetary disk is discussed.

We present new X-ray and Far Ultraviolet (FUV) observations of T Tauri stars covering the age range 1 to 10 Myr. Our goals are to observationally constrain the intensity of radiation fields responsible for evaporating gas from the circumstellar disk and to assess the feasibility of current photoevaporation models, focusing on X-ray and UV radiation. We greatly increase the number of 7-10 Myr old T Tauri stars observed in the X-rays by including observations of the well populated 25 Ori aggregate in the Orion OB1a subassociation. With these new 7-10 Myr objects, we confirm that X-ray emission remains constant from 1-10 Myr. We also show, for the first time, observational evidence for the evolution of FUV radiation fields with a sample of 56 accreting and non-accreting young stars spanning 1 Myr to 1 Gyr. We find that the FUV emission decreases on timescales consistent with the decline of accretion in classical T Tauri stars until reaching the chromospheric level in weak T Tauri stars and debris disks. Overall, we find that the observed strength of high energy radiation is consistent with that required by photoevaporation models to dissipate the disks in timescales of approximately 10 Myr. Finally, we find that the high energy fields that affect gas evolution are not similarly affecting dust evolution; in particular, we find that disks with inner clearings, the transitional disks, have similar levels of FUV emission as full disks.

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

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

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

We report the results of a large-area study of water vapor along the Orion Molecular Cloud ridge, the purpose of which was to determine the depth-dependent distribution of gas-phase water in dense molecular clouds. We find that the water vapor measured toward 77 spatial positions along the face-on Orion ridge, excluding positions surrounding the outflow associated with BN/KL and IRc2, display integrated intensities that correlate strongly with known cloud surface tracers such as CN, C2H, 13CO J =5-4, and HCN, and less well with the volume tracer N2H+. Moreover, at total column densities corresponding to Av < 15 mag., the ratio of H2O to C18O integrated intensities shows a clear rise approaching the cloud surface. We show that this behavior cannot be accounted for by either optical depth or excitation effects, but suggests that gas-phase water abundances fall at large Av. These results are important as they affect measures of the true water-vapor abundance in molecular clouds by highlighting the limitations of comparing measured water vapor column densities with such traditional cloud tracers as 13CO or C18O. These results also support cloud models that incorporate freeze-out of molecules as a critical component in determining the depth-dependent abundance of water vapor.

We present results from a model of the chemical evolution of protoplanetary disks. In our models we directly calculate the changing propagation and penetration of a high energy radiation field with Lyman alpha radiation included. We also explore the effect on our models of including dust grain settling. We find that, in agreement with earlier studies, the evolution of dust grains plays a large role in determining how deep the UV radiation penetrates into the disk. Significant grain settling at the midplane leads to much smaller freeze-out regions and a correspondingly larger molecular layer, which leads to an increase in column density for molecular species such as CO, CN and SO. The inclusion of Lyman alpha radiation impacts the disk chemistry through specific species that have large photodissociation cross sections at 1216 A. These include HCN, NH3 and CH4, for which the column densities are decreased by an order of magnitude or more due to the presence of Lyman alpha radiation in the UV spectrum. A few species, such as CO2 and SO, are enhanced by the presence of Lyman alpha radiation, but rarely by more than a factor of a few.

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.

We present an analysis of Spitzer-IRS spectroscopic maps of the L1157 protostellar outflow in the H2 pure-rotational lines from S(0) to S(7). The aim of this work is to derive the physical conditions pertaining to the warm molecular gas and study their variations within the flow. The mid-IR H2 emission follows the morphology of the precessing flow, with peaks correlated with individual CO clumps and H2 2.12{\mu}m ro-vibrational emission. More diffuse emission delineating the CO cavities is detected only in the low-laying transitions, with J(lower) less or equal to 2. The H2 line images have been used to construct 2D maps of N(H2), H2 ortho-to-para ratio and temperature spectral index beta, in the assumption of a gas temperature stratification where the H2 column density varies as T^(beta). Variations of these parameters are observed along the flow. In particular, the ortho-to-para ratio ranges from 0.6 to 2.8, highlighting the presence of regions subject to recent shocks where the ortho-to-para ratio has not had time yet to reach the equilibrium value. Near-IR spectroscopic data on ro-vibrational H2 emission have been combined with the mid-IR data and used to derive additional shock parameters in the brightest blue- and red-shifted emission knots. A high abundance of atomic hydrogen (H/H2 about 0.1-0.3) is implied by the observed H2 column densities, assuming n(H2) values as derived by independent SiO observations. The presence of a high fraction of atomic hydrogen, indicates that a partially-dissociative shock component should be considered for the H2 excitation in these localized regions. However, planar shock models, either of C- or J-type, are not able to consistently reproduce all the physical parameters derived from our analysis of the H2 emission. Globally, H2 emission contributes to about 50% of the total shock radiated energy in the L1157 outflow.

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

Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow.

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.

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

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.

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

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

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

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77{\deg}). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75{\deg} but overestimate the volume mixing ratio compared to our measurements.

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

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

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

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

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.

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

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

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

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

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

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

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

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

Chemistry plays an important role in the structure and evolution of protoplanetary disks, with implications for the composition of comets and planets. This is the first of a series of papers based on data from DISCS, a Submillimeter Array survey of the chemical composition of protoplanetary disks. The six Taurus sources in the program (DM Tau, AA Tau, LkCa 15, GM Aur, CQ Tau and MWC 480) range in stellar spectral type from M1 to A4 and offer an opportunity to test the effects of stellar luminosity on the disk chemistry. The disks were observed in 10 different lines at ~3" resolution and an rms of ~100 mJy beam-1 at ~0.5 km s-1. The four brightest lines are CO 2-1, HCO+ 3-2, CN 2_3-1_2 and HCN 3-2 and these are detected toward all sources (except for HCN toward CQ Tau). The weaker lines of CN 2_2-1_1, DCO+ 3-2, N2H+ 3-2, H2CO 3_03-2_02 and 4_14-3_13 are detected toward two to three disks each, and DCN 3-2 only toward LkCa 15. CH3OH 4_21-3_12 and c-C3H2 are not detected. There is no obvious difference between the T Tauri and Herbig Ae sources with regard to CN and HCN intensities. In contrast, DCO+, DCN, N2H+ and H2CO are detected only toward the T Tauri stars, suggesting that the disks around Herbig Ae stars lack cold regions for long enough timescales to allow for efficient deuterium chemistry, CO freeze-out, and grain chemistry.

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

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.

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

The high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

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

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

We report on the initial analysis of a Herschel/PACS full range spectrum of Neptune, covering the 51-220 micrometer range with a mean resolving power of ~ 3000, and complemented by a dedicated observation of CH4 at 120 micrometers. Numerous spectral features due to HD (R(0) and R(1)), H2O, CH4, and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~ 3 K than inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5+/-1) X 10**-5, confirming the enrichment of Neptune in deuterium over the protosolar value (~ 2.1 X 10**-5); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+/-0.2) X 10**-3, and CH4 appears to decrease in the lower stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south polar region; (iv) the H2O stratospheric column is (2.1+/-0.5) X 10**14 cm-2 but its vertical distribution is still to be determined, so the H2O external flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming the dual origin of CO suspected from ground-based millimeter/submillimeter observations.

We have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPIRE Fourier Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum contains the characteristic transitions of 12CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow the comparison to global atmospheric models allowing the mean mixing ratios of water and 12CO to be investigated. We find that it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10-4 and a 12CO mixing ratio of 9 x 10-4

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

High-resolution far-infrared and sub-millimetre spectroscopy of water lines is an important tool to understand the physical and chemical properties of cometary atmospheres. We present observations of several rotational ortho- and para-water transitions in comet C/2008 Q3 (Garradd) performed with HIFI on Herschel. These observations have provided the first detection of the 2_{12}-1_{01} (1669 GHz) ortho and 1_{11}-0_{00} (1113 GHz) para transitions of water in a cometary spectrum. In addition, the ground-state transition 1_{10}-1_{01} at 557 GHz is detected and mapped. By detecting several water lines quasi-simultaneously and mapping their emission we can constrain the excitation parameters in the coma. Synthetic line profiles are computed using excitation models which include excitation by collisions, solar infrared radiation, and radiation trapping. We obtain the gas kinetic temperature, constrain the electron density profile, and estimate the coma expansion velocity by analyzing the map and line shapes. We derive water production rates of 1.7-2.8 x 10^{28} s^{-1} over the range r_h = 1.83-1.85 AU.

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

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

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

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

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

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel. The PACS instrument acquired images of the dust coma in 70- and 160-micrometers filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired with SPIRE. The comet emission continuum from 70 to 672 micrometers was measured, but no lines were detected. The spectral energy distribution indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and Nancay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1.

Context: Chemical models of dense cloud cores often utilize the so-called pseudo-time-dependent approximation, in which the physical conditions are held fixed and uniform as the chemistry occurs. In this approximation, the initial abundances chosen, which are totally atomic in nature except for molecular hydrogen, are artificial. A more detailed approach to the chemistry of dense cold cores should include the physical evolution during their early stages of formation. Aims: Our major goal is to investigate the initial synthesis of molecular ices and gas-phase molecules as cold molecular gas begins to form behind a shock in the diffuse interstellar medium. The abundances calculated as the conditions evolve can then be utilized as reasonable initial conditions for a theory of the chemistry of dense cores. Methods: Hydrodynamic shock-wave simulations of the early stages of cold core formation are used to determine the time-dependent physical conditions for a gas-grain chemical network. We follow the cold post-shock molecular evolution of ices and gas-phase molecules for a range of visual extinction up to AV ~ 3, which increases with time. At higher extinction, self-gravity becomes important. Results: As the newly condensed gas enters its cool post-shock phase, a large amount of CO is produced in the gas. As the CO forms, water ice is produced on grains, while accretion of CO produces CO ice. The production of CO2 ice from CO occurs via several surface mechanisms, while the production of CH4 ice is slowed by gas-phase conversion of C into CO.

Despite a surface dominated by carbon-based life, the bulk composition of the Earth is dramatically carbon poor when compared to the material available at formation. Bulk carbon deficiency extends into the asteroid belt representing a fossil record of the conditions under which planets are born. The initial steps of planet formation involve the growth of primitive sub-micron silicate and carbon grains in the Solar Nebula. We present a solution wherein primordial carbon grains are preferentially destroyed by oxygen atoms ignited by heating due to stellar accretion at radii < 5 AU. This solution can account for the bulk carbon deficiency in the Earth and meteorites, the compositional gradient within the asteroid belt, and for growing evidence for similar carbon deficiency in rocks surrounding other stars.

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

We analyze the far-ultraviolet (FUV) spectra of 33 classical T Tauri stars (CTTS), including 20 new spectra obtained with the Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) on the Hubble Space Telescope. Of the sources, 28 are in the ~1 Myr old Taurus-Auriga complex or Orion Molecular Cloud, 4 in the 8-10 Myr old Orion OB1a complex and one, TW Hya, in the 10 Myr old TW Hydrae Association. We also obtained FUV ACS/SBC spectra of 10 non-accreting sources surrounded by debris disks with ages between 10 and 125 Myr. We use a feature in the FUV spectra due mostly to electron impact excitation of \h2 to study the evolution of the gas in the inner disk. We find that the \h2 feature is absent in non-accreting sources, but is detected in the spectra of CTTS and correlates with accretion luminosity. Since all young stars have active chromospheres which produce strong X-ray and UV emission capable of exciting \h2 in the disk, the fact that the non-accreting sources show no \h2 emission implies that the \h2 gas in the inner disk has dissipated in the non-accreting sources, although dust (and possibly gas) remains at larger radii. Using the flux at 1600 {\AA}, we estimate that the column density of \h2 left in the inner regions of the debris disks in our sample is less than ~ 3x10^-6 g cm^-2, nine orders of magnitude below the surface density of the minimum mass solar nebula at 1 AU.

The origins of planets, and perhaps life itself, is intrinsically linked to the chemistry of planet formation. In this chapter we will attempt to explore the chemistry of planet-forming disks from the perspective of knowledge gained from decades of solar system study, combined with our rapidly growing knowledge of extra-solar protoplanetary disks. Our aim is not to provide a comprehensive survey of the literature, but rather to present the chemical and physical tools that are used to explore chemical evolution during the earliest stages of planet formation. We include a discussion of observational constraints from our solar system and extra-solar disks, thermodynamic equilibrium and chemical kinetics, and the physical prescriptions for various important aspects (X-ray/UV irradiation, thermal structure, molecular adsorption/desorption, ..) that are used in current models. The chapter is aimed for young researchers.

Using a source selection biased towards high mass star forming regions, we used a Large Velocity Gradient (LVG) code to calculate the H2 densities and CS column densities for a sample of Midcourse Space Experiment (MSX) 8 micron infrared dark cores. Our average H2 density and CS column density were 1.14 x 10e6 cm-3 and 1.21 x 10e13 cm-2 respectively. In addition, we have calculated the Jeans mass and Virial mass for each core to get a better understanding of their gravitational stability. We found that core masses calculated from observations of N2H+ J = 1-0 and C18O J = 1-0 by Ragan et al. 2006 (Paper 1) were sufficient for collapse, though most regions are likely to form protoclusters. We have explored the star-forming properties of the molecular gas within our sample and find some diversity which extends the range of infrared dark clouds from very the massive clouds that will create large clusters, to clouds that are similar to some of our local counterparts (e.g. Serpens, Ophiuchus).

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

We have observed the Oph A-N6 prestellar core in the following transitions: N2D+ J=3 to 2, DCO+ J=3 to 2 and J=5 to 4, HCO+ J=3 to 2, CS J=5 to 4 and J=7 to 6, and H13CO+ J=3 to 2 and J=4 to 3, using the James Clerk Maxwell Telescope. We also observed the NH3 (1,1) and (2,2) inversion transitions towards the Oph A-N6 peak with the Green Bank Telescope. We have found that the N6 core is composed of shells of different chemical composition due to the freezing out of chemical species at different densities. The undepleted species N2D+ appears to trace the high-density interior of the core, DCO+ and H13CO+ trace an intermediate region, and CS traces the outermost edges of the core. A distinct blue-red spectral asymmetry, indicative of infall motion, is clearly detected in the HCO+ spectra, suggesting that N6 is undergoing gravitational collapse. This collapse was possibly initiated by a decrease in turbulent support suggested by the fact that the non-thermal contribution to the line widths is smaller for the molecular species closer to the center of the core. We also present a temperature profile for the core. These observations support the claim that the Oph A-N6 core is an extremely young prestellar core, which may have been recently cut off from MHD support and begun to collapse.

We have conducted a survey of a sample of infrared-dark clouds (IRDCs) with the Spitzer Space Telescope in order to explore their mass distribution. We present a method for tracing mass using dust absorption against the bright Galactic background at 8 microns. The IRDCs in this sample are comprised of tens of clumps, ranging in sizes from 0.02 to 0.3 pc in diameter and masses from 0.5 to a few 10 Msun, the broadest dynamic range in any clump mass spectrum study to date. Structure with this range in scales confirms that IRDCs are the the precursors to stellar clusters in an early phase of fragmentation. Young stars are distributed in the vicinity of the IRDCs, but the clumps are typically not associated with stars and appear pre-stellar in nature. We find an IRDC clump mass spectrum with a slope of 1.76 +/- 0.05 for masses from 30 to 3000 Msun. This slope is consistent with numerous studies, culled from a variety of observational techniques, of massive star formation regions and is close to the mass function of Galactic stellar clusters and star clusters in other galaxies. We assert that the shape of the mass function is an intrinsic and universal feature of massive star formation regions, that are the birth sites of stellar clusters. As these clouds evolve and their constituent clumps fragment, the mass spectrum will steepen and eventually assume the form of the core mass function that is observed locally.

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

We model the temperature and chemical structure of molecular clouds as a function of depth into the cloud, assuming a cloud of constant density n illuminated by an external FUV (6 eV < E < 13.6 eV) flux G_0 (scaling factor in multiples of the local interstellar field). Extending previous photodissociation region models, we include the freezing of species, simple grain surface chemistry, and desorption (including FUV photodesorption) of ices. We also treat the opaque cloud interior with time-dependent chemistry. Here, under certain conditions, gas phase elemental oxygen freezes out as water ice and the elemental C/O abundance ratio can exceed unity, leading to complex carbon chemistry. Gas phase H2O and O2 peak in abundance at intermediate depth into the cloud, roughly A_V~3-8 from the surface, the depth proportional to ln(G_0/n). Closer to the surface, molecules are photodissociated. Deeper into the cloud, molecules freeze to grain surfaces. At intermediate depths photodissociation rates are attenuated by dust extinction, but photodesorption prevents total freezeout. For G_0 < 500, abundances of H2O and O2 peak at values ~10^(-7), producing columns ~10^(15) per cm^2, independent of G_0 and n. The peak abundances depend primarily on the product of the photodesorption yield of water ice and the grain surface area per H nucleus. At higher values of G_0, thermal desorption of O atoms from grains enhances the gas phase H2O peak abundance and column slightly, whereas the gas phase O2 peak abundance rises to ~10^(-5) and the column to ~2x10^(16) per cm^2. We present simple analytic equations for the abundances as a function of depth which clarify the dependence on parameters. The models are applied to observations of H2O, O2, and water ice in a number of sources, including B68, NGC 2024, and Rho Oph.

Massive star formation exhibits an extremely rich chemistry. However, not much evolutionary details are known yet, especially at high spatial resolution. Therefore, we synthesize previously published Submillimeter Array high-spatial-resolution spectral line observations toward four regions of high-mass star formation that are in various evolutionary stages with a range of luminosities. Estimating column densities and comparing the spatially resolved molecular emission allows us to characterize the chemical evolution in more detail. Furthermore, we model the chemical evolution of massive warm molecular cores to be directly compared with the data. The four regions reveal many different characteristics. While some of them, e.g., the detection rate of CH3OH, can be explained by variations of the average gas temperatures, other features are attributed to chemical effects. For example, C34S is observed mainly at the core-edges and not toward their centers because of temperature-selective desorption and successive gas-phase chemistry reactions. Most nitrogen-bearing molecules are only found toward the hot molecular cores and not the earlier evolutionary stages, indicating that the formation and excitation of such complex nitrogen-bearing molecules needs significant heating and time to be fully developed. Furthermore, we discuss the observational difficulties to study massive accretion disks in the young deeply embedded phase of massive star formation. The general potential and limitations of such kind of dataset are discussed, and future directions are outlined. The analysis and modeling of this source sample reveals many interesting features toward a chemical evolutionary sequence. However, it is only an early step, and many observational and theoretical challenges in that field lie ahead.

We report the results of spectroscopic mapping observations carried out toward a ~1 min x 1 min region within the northern lobe of the outflow from NGC 2071 using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations covered the 5.2-37 um spectral region and have led to the detection of a number of ionic, atomic, and molecular lines, including fine-structure emission of Si+, Fe+, S++, S, the S(0)-S(7) pure rotational lines of H2, the R(3) and R(4) transitions of HD, and at least 11 transitions of H2O. In addition, the 6.2, 7.4, 7.6, 7.9, 8.6 and 11.3 um PAH emission bands were also observed and several transitions of OH were tentatively detected. Most of the detected line transitions were strong enough to map including, for the first time, three transitions of hot H2O. We find that: (1) the water emission is extended; (2) the extended emission is aligned with the outflow; and, (3) the spatial distribution of the water emission generally follows that observed for H2. Based on the measured line intensities, we derive an HD abundance relative to H2 of 1.1-1.8 10^-5 and an H2O number density of 12-2 cm^3. The H2 density in the water-emitting region is not well constrained by our observations, but is likely between 3 10^4 and 10^6 cm^3, yielding an H2O abundance relative to H2 of between 2 10^-5 and 6 10^-4. Future observations planned for the Herschel Space Observatory should greatly improve the density estimate, and thus our knowledge of the H2O abundance, for the water-emitting regions reported here. Finally, we note a possible departure from the H2O ortho-to-para ratio of 3:1 expected for water formed in hot post-shocked gas, suggesting that a significant fraction of the water vapor we detect may arise from H2O sputtered from cold dust grains.

We present results from a model of oxygen isotopic anomaly production through selective photodissociation of CO within the collapsing proto-Solar cloud. Our model produces a proto-Sun with a wide range of Delta_17O values depending on the intensity of the ultraviolet radiation field. Dramatically different results from two recent Solar wind oxygen isotope measurements indicate that a variety of compositions remain possible for the solar oxygen isotope composition. However, constrained by other measurements from comets and meteorites, our models imply the birth of the Sun in a stellar cluster with an enhanced radiation field and are therefore consistent with a supernova source for 60Fe in meteorites.

We present SWAS detections of the ground-state 1(10)-1(01) transition of o-H2O at 557 GHz in 18 molecular outflows. These results are combined with ground-based observations of the J=1-0 transitions of 12CO and 13CO obtained at the FCRAO and, for a subset of the outflows, data from ISO. Assuming the SWAS water line emission originates from the same gas traced by CO emission, we find that the outflowing gas in most outflows has an o-H2O abundance relative to H2 of between 10(-7) and 10(-6). Analysis of the water abundance as a function of outflow velocity reveals a strong dependence. The water abundance increases with velocity, and at the highest outflow velocities some outflows have relative o-H2O abundances of order 10(-4). However the mass of gas with such elevated water abundances represents less that 1% of the total outflow gas mass. The ISO LWS observations of high-J rotational lines of CO and the 179.5 micron transition of o-H2O provide evidence for a warmer outflow component than required to produce either the SWAS or FCRAO lines. The mass associated with the ISO emission is similar to that responsible for the highest velocity water emission detected by SWAS. The bulk of the outflowing gas has an abundance of o-H2O well below what would be expected if the gas has passed through a C-shock with shock velocities greater than 10 km/s. Gas-phase water can be depleted in the post-shock gas due to freeze-out onto grain mantles, however the rate of freeze-out is too slow to explain our results. Therefore we believe that only a small fraction of the outflowing molecular gas has passed through shocks strong enough to fully convert the gas-phase oxygen to water. This result has implications for the acceleration mechanism of the molecular gas in these outflows.

We report the first fully sampled maps of the distribution of interstellar CO2 ices, H2O ices and total hydrogen nuclei, as inferred from the 9.7 micron silicate feature, toward the star-forming region Cepheus A East with the IRS instrument onboard the Spitzer Space Telescope. We find that the column density distributions for these solid state features all peak at, and are distributed around, the location of HW2, the protostar believed to power one of the outflows observed in this star-forming region. A correlation between the column density distributions of CO2 and water ice with that of total hydrogen indicates that the solid state features we mapped mostly arise from the same molecular clumps along the probed sight lines. We therefore derive average CO2 ice and water ice abundances with respect to the total hydrogen column density of X(CO2)_ice~1.9x10^-5 and X(H2O)_ice~7.5x10^-5. Within errors, the abundances for both ices are relatively constant over the mapped region exhibiting both ice absorptions. The fraction of CO2 ice with respect to H2O ice is also relatively constant at a value of 22% over that mapped region. A clear triple-peaked structure is seen in the CO2 ice profiles. Fits to those profiles using current laboratory ice analogs suggest the presence of both a low-temperature polar ice mixture and a high-temperature methanol-rich ice mixture along the probed sightlines. Our results further indicate that thermal processing of these ices occurred throughout the sampled region.

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

While most sources above 10^5Lsun have already formed an Ultracompact HII region (UCHII), this is not necessarily the case for sources of lower luminosity. Characterizing sources in the transition phase, i.e., very luminous objects without any detectable free-free emission, is important for a general understanding of massive star formation. Therefore, we observed the luminous High-Mass Protostellar Object IRAS23151+5912 with the Submillimeter Array at 875mum in the submm continuum and spectral line emission at sub-arcsecond resolution. The 875mum submm continuum emission has been resolved into at least two condensations. The previously believed driving source of one of the outflows, the infrared source IRS1, is ~0.9'' offset from the main submm peak. Over the entire 4GHz bandwidth we detect an intermediate dense spectral line forest with 27 lines from 8 different species, isotopologues or vibrationally-torsionally excited states. Temperature estimates based on the CH3OH line series result in values of T(Peak1)~150+-50K and T(Peak2)~80~30K for the two submm peak positions, respectively. The SiO(8-7) red- and blue-shifted line maps indicate the presence of two molecular outflows. In contrast, the vibrationally-torsionally excited CH3OH line exhibits a velocity gradient approximately perpendicular to one of the outflows. With a size of approximately 5000AU and no Keplerian rotation signature, this structure does not resemble a genuine accretion disk but rather a larger-scale rotating toroid that may harbor a more common accretion disk at its so far unresolved center.

Here we present the Spitzer IRS spectrum of CS Cha, a member of the ~2 Myr old Chamaeleon star-forming region, which reveals an optically thick circumstellar disk truncated at ~43 AU, the largest hole modeled in a transitional disk to date. Within this inner hole, ~5x10^-5 lunar masses of dust are located in a small optically thin inner region which extends from 0.1 to 1 AU. In addition, the disk of CS Cha has bigger grain sizes and more settling than the previously modeled transitional disks DM Tau, GM Aur, and CoKu Tau/4, suggesting that CS Cha is in a more advanced state of dust evolution. The Spitzer IRS spectrum also shows [Ne II] 12.81 micron fine-structure emission with a luminosity of 1.3x10^29 ergs s^-1, indicating that optically thin gas is present in this ~43 AU hole, in agreement with H_alpha measurements and a UV excess which indicate that CS Cha is still accreting 1.2x10^-8 M_sun yr^-1. We do not find a correlation of the [Ne II] flux with L_X, however, there is a possible correlation with mass accretion rate, which if confirmed would suggest that EUV fluxes due to accretion are the main agent for formation of the [Ne II] line.

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

High resolution molecular line observations of CS, HCO+, C18O and N2H+ were obtained toward the starless globule FeSt 1-457 in order to investigate its kinematics and chemistry. The HCO+ and CS spectra show clear self-reversed and asymmetric profiles across the face of the globule. The sense of the observed asymmetry is indicative of the global presence of expansion motions in the outer layers of the globule. These motions appear to be subsonic and significantly below the escape velocity of the globule. Comparison of our observations with near-infrared extinction data indicate that the globule is gravitationally bound. Taken together these considerations lead us to suggest that the observed expansion has its origin in an oscillatory motion of the outer layers of the globule which itself is likely in a quasi-stable state near hydrostatic equilibrium. Analysis of the observed linewidths of CO and N2H+ confirm that thermal pressure is the dominant component of the cloud's internal support. A simple calculation suggests that the dominant mode of pulsation would be an l = 2 mode with a period of 0.3 Myr. Deformation of the globule due to the large amplitude l = 2 oscillation may be responsible for the double-peaked structure of the core detected in high resolution extinction maps. Detailed comparison of the molecular-line observations and extinction data provides evidence for significant depletion of C18O and perhaps HCO+ while N2H+ may be undepleted to a cloud depth of about 40 magnitudes of visual extinction.

High resolution molecular line observations of CS, HCO+, C18O and N2H+ were obtained toward the starless globule FeSt 1-457 in order to investigate its kinematics and chemistry. The HCO+ and CS spectra show clear self-reversed and asymmetric profiles across the face of the globule. The sense of the observed asymmetry is indicative of the global presence of expansion motions in the outer layers of the globule. These motions appear to be subsonic and significantly below the escape velocity of the globule. Comparison of our observations with near-infrared extinction data indicate that the globule is gravitationally bound. Taken together these considerations lead us to suggest that the observed expansion has its origin in an oscillatory motion of the outer layers of the globule which itself is likely in a quasi-stable state near hydrostatic equilibrium. Analysis of the observed linewidths of CO and N2H+ confirm that thermal pressure is the dominant component of the cloud's internal support. A simple calculation suggests that the dominant mode of pulsation would be an l = 2 mode with a period of 0.3 Myr. Deformation of the globule due to the large amplitude l = 2 oscillation may be responsible for the double-peaked structure of the core detected in high resolution extinction maps. Detailed comparison of the molecular-line observations and extinction data provides evidence for significant depletion of C18O and perhaps HCO+ while N2H+ may be undepleted to a cloud depth of about 40 magnitudes of visual extinction.

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

We report the results of spectroscopic mapping observations carried out toward small (1 x 1 arcmin) regions within the supernova remnants W44, W28, IC443, and 3C391 using the Infrared Spectrograph of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, have led to the detection of a total of 15 fine structure transitions of Ne+, Ne++, Si+, P+, S, S++, Cl+, Fe+, and Fe++; the S(0) - S(7) pure rotational lines of molecular hydrogen; and the R(3) and R(4) transitions of hydrogen deuteride. In addition to these 25 spectral lines, the 6.2, 7.7, 8.6, 11.3 and 12.6 micron PAH emission bands were also observed. Most of the detected line transitions have proven strong enough to map in several sources, providing a comprehensive picture of the relative distribution of the various line emissions observable in the Spitzer/IRS bandpass. A principal component analysis of the spectral line maps reveals that the observed emission lines fall into five distinct groups, each of which may exhibit a distinct spatial distribution: (1) lines of S and H2 (J > 2); (2) the H2 S(0) line; (3) lines of ions with appearance potentials less than 13.6 eV; (4) lines of ions with appearance potentials greater than 13.6 eV, not including S++; (5) lines of S++. Lines of group (1) likely originate in molecular material subject to a slow, nondissociative shock that is driven by the overpressure within the supernova remnant, and lines in groups (3) - (5) are associated primarily with dissociative shock fronts with a range of (larger) shock velocities. The H2 S(0) line shows a low-density diffuse emission component, and - in some sources - a shock-excited component.

Aiming at a better understand of the physical and chemical processes in the hot molecular core stage of high-mass star formation, we observed the prototypical hot core G29.96-0.02 in the 862mu band with the Submillimeter Array (SMA) at sub-arcsecond spatial resolution. The observations resolved the hot molecular core into six submm continuum sources with the finest spatial resolution of 0.36''x0.25'' (~1800AU) achieved so far. Four of them located within 7800(AU)^2 comprise a proto-Trapezium system with estimated protostellar densities of 1.4x0^5 protostars/pc^3. The plethora of ~80 spectral lines allows us to study the molecular outflow(s), the core kinematics, the temperature structure of the region as well as chemical effects. The derived hot core temperatures are of the order 300K. We find interesting chemical spatial differentiations, e.g., C34S is deficient toward the hot core and is enhanced at the UCHII/hot core interface, which may be explained by temperature sensitive desorption from grains and following gas phase chemistry. The SiO(8-7) emission outlines likely two molecular outflows emanating from this hot core region. Emission from most other molecules peaks centrally on the hot core and is not dominated by any individual submm peak. Potential reasons for that are discussed. A few spectral lines that are associated with the main submm continuum source, show a velocity gradient perpendicular to the large-scale outflow. Since this velocity structure comprises three of the central protostellar sources, this is not a Keplerian disk. While the data are consistent with a gas core that may rotate and/or collapse, we cannot exclude the outflow(s) and/or nearby expanding UCHII region as possible alternative causes of this velocity pattern.

We investigate the structure of the core surrounding the recently identified deeply embedded young stellar object Barnard 1c which has an unusual polarization pattern as traced in submillimeter dust emission. Barnard 1c lies within the Perseus molecular cloud at a distance of 250 pc. It is a deeply embedded core of 2.4 solar masses (Kirk et al.) and a luminosity of 4 +/- 2 solar luminosities. Observations of CO, 13CO, C18O, HCO+ and N2H+ were obtained with the BIMA array, together with the continuum at 3.3 mm and 2.7 mm. Single-dish measurements of N2H+ and HCO+ with FCRAO reveal the larger scale emission in these lines, The CO and HCO+ emission traces the outflow, which coincides in detail with the S-shaped jet recently found in Spitzer IRAC imaging. The N2H+ emission, which anticorrelates spatially with the C18O emission, originates from a rotating envelope with effective radius ~ 2400 AU and mass 2.1 - 2.9 solar masses. N2H+ emission is absent from a 600 AU diameter region around the young star. The remaining N2H+ emission may lie in a coherent torus of dense material. With its outflow and rotating envelope, B1c closely resembles the previously studied object L483-mm, and we conclude that it is a protostar in an early stage of evolution. We hypothesize that heating by the outflow and star has desorbed CO from grains which has destroyed N2H+ in the inner region and surmise that the presence of grains without ice mantles in this warm inner region can explain the unusual polarization signature from B1c.

We report the first detection of gas-phase CO2 emission in the star-forming region Cepheus A East, obtained by spectral line mapping of the v2 bending mode at 14.98 micron with the Infrared Spectrograph (IRS) instrument onboard the Spitzer Space Telescope. The gaseous CO2 emission covers a region about 35'' x 25'' in extent, and results from radiative pumping by 15 micron continuum photons emanating predominantly from the HW2 protostellar region. The gaseous CO2 exhibits a temperature distribution ranging from 50 K to 200 K. A correlation between the gas-phase CO2 distribution and that of H2 S(2), a tracer of shock activity, indicates that the CO2 molecules originate in a cool post-shock gas component associated with the outflow powered by HW2. The presence of CO2 ice absorption features at 15.20 micron toward this region and the lack of correlation between the IR continuum emission and the CO2 gas emission distribution further suggest that the gaseous CO2 molecules are mainly sputtered off grain mantles -- by the passage of slow non-dissociative shocks with velocities of 15-30 km/s -- rather than sublimated through grain heating.

We report the detection of interstellar hydrogen deuteride (HD) toward the supernova remnant IC443, and the tentative detection of HD toward the Herbig Haro objects HH54 and HH7 and the star forming region GGD37 (Cepheus A West). Our detections are based upon spectral line mapping observations of the R(3) and R(4) rotational lines of HD, at rest wavelengths of 28.502 and 23.034 micron respectively, obtained using the Infrared Spectrograph onboard the Spitzer Space Telescope. The HD R(4)/R(3) line intensity ratio promises to be a valuable probe of the gas pressure in regions where it can be observed. The derived HD/H2 abundance ratios are 1.19(+0.35/-0.24)E-5, 1.80(+0.54/-0.32)E-5, and 1.41(+0.46/-0.33)E-5 respectively (68.3% confidence limits, based upon statistical errors alone) for IC443 (clump C), HH54, and HH7. If HD is the only significant reservoir of gas-phase deuterium in these sources, the inferred HD/H2 ratios are all consistent with a gas-phase elemental abundance [n(D)/n(H)](gas) ~ 7.5E-6, a factor 2 - 3 below the values obtained previously from observations of atomic deuterium in the local bubble and the Galactic halo. However, similarly low gas-phase deuterium abundances have been inferred previously for molecular gas clouds in the Orion region, and in atomic clouds along sight-lines within the Galactic disk to stars more distant than 500 pc from the Sun.

We have identified 41 infrared dark clouds from the 8 micron maps of the Midcourse Space Experiment (MSX), selected to be found within one square degree areas centered on known ultracompact HII regions. We have mapped these infrared dark clouds in N2H+(1-0), CS(2-1) and C18O(1-0) emission using the Five College Radio Astronomy Observatory. The maps of the different species often show striking differences in morphologies, indicating differences in evolutionary state and/or the presence of undetected, deeply embedded protostars. We derive an average mass for these clouds using N2H+ column densities of ~2500 solar masses, a value comparable to that found in previous studies of high mass star forming cores using other mass tracers. The linewidths of these clouds are typically ~2.0 - 2.9 km/s. Based on the fact that they are dark at 8 micron, compact, massive, and have large velocity dispersions, we suggest that these clouds may be the precursor sites of intermediate and high mass star formation.

We report the results of spectroscopic mapping observations carried out toward the Herbig-Haro objects HH7-11 and HH54 over the 5.2 - 37 micron region using the Infrared Spectrograph of the Spitzer Space Telescope. These observations have led to the detection and mapping of the S(0) - S(7) pure rotational lines of molecular hydrogen, together with emissions in fine structure transitions of Ne+, Si+, S, and Fe+. The H2 rotational emissions indicate the presence of warm gas with a mixture of temperatures in the range 400 - 1200 K, consistent with the expected temperature behind nondissociative shocks of velocity ~ 10 - 20 km/s, while the fine structure emissions originate in faster shocks of velocity 35 - 90 km/s that are dissociative and ionizing. Maps of the H2 line ratios reveal little spatial variation in the typical admixture of gas temperatures in the mapped regions, but show that the H2 ortho-to-para ratio is quite variable, typically falling substantially below the equilibrium value of 3 attained at the measured gas temperatures. The non-equilibrium ortho-to-para ratios are characteristic of temperatures as low as ~ 50 K, and are a remnant of an earlier epoch, before the gas temperature was elevated by the passage of a shock. Correlations between the gas temperature and H2 ortho-to-para ratio show that ortho-to-para ratios < 0.8 are attained only at gas temperatures below ~ 900 K; this behavior is consistent with theoretical models in which the conversion of para- to ortho-H2 behind the shock is driven by reactive collisions with atomic hydrogen, a process which possesses a substantial activation energy barrier (E_A/k ~ 4000 K) and is therefore very inefficient at low temperature.

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

In this review we re-evaluate our observational and theoretical understanding of the chemical evolution of protoplanetary disks. We discuss how improved observational capabilities have enabled the detection of numerous molecules exposing an active disk chemistry that appears to be in disequilibrium. We outline the primary facets of static and dynamical theoretical chemical models. Such models have demonstrated that the observed disk chemistry arises from warm surface layers that are irradiated by X-ray and FUV emission from the central accreting star. Key emphasis is placed on reviewing areas where disk chemistry and physics are linked: including the deuterium chemistry, gas temperature structure, disk viscous evolution (mixing), ionization fraction, and the beginnings of planet formation.

We present the first images of the 691.473 GHz CO J=6-5 line in a protoplanetary disk, obtained along with the 690 GHz dust continuum, toward the classical T Tauri star TW Hya using the Submillimeter Array. Imaging in the CO J=6-5 line reveals a rotating disk, consistent with previous observations of CO J=3-2 and 2-1 lines. Using an irradiated accretion disk model and 2D Monte Carlo radiative transfer, we find that additional surface heating is needed to fit simultaneously the absolute and relative intensities of the CO J=6-5, 3-2 and 2-1 lines. In particular, the vertical gas temperature gradient in the disk must be steeper than that of the dust, mostly likely because the CO emission lines probe nearer to the surface of the disk. We have used an idealized X-ray heating model to fit the line profiles of CO J=2-1 and 3-2 with Chi-square analysis, and the prediction of this model yields CO J=6-5 emission consistent with the observations.

We compare an evolutionary chemical model with simple empirical models of the abundance and with static chemical models. We focus on the prediction of molecular line profiles that are commonly observed in low mass star forming cores. We show that empirical models can be used to constrain evaporation radii and infall radii using lines of some species. Species with more complex abundance profiles are not well represented by the empirical models. Static chemical models produce abundance profiles different from those obtained from an evolutionary calculation because static models do not account for the flow of matter inward from the outer regions. The resulting profiles of lines used to probe infall may differ substantially.

We present the first Spitzer Infrared Spectrograph observations of the 15.2 micron bending mode of CO2 ice towards field stars behind a quiescent dark cloud. CO2 ice is detected towards 2 field stars (Elias 16, Elias 3) and a single protostar (HL Tau) with anabundance of ~15-20% relative to water ice. CO2 ice is not detected towards the source with lowest extinction in our sample, Tamura 17 (A_V = 3.9m). A comparison of the Elias 16 spectrum with laboratory data demonstrates that the majority of CO2 ice is embedded in a polar H2O-rich ice component, with ~15% of CO2 residing in an apolar H2O-poor mantle. This is the first detection of apolar CO2 towards a field star. We find that the CO2 extinction threshold is A_V = 4m +/- 1m, comparable to the threshold for water ice, but significantly less than the threshold for CO ice, the likely precursor of CO2. Our results confirm CO2 ice forms in tandem with H2O ice along quiescent lines of sight. This argues for CO2 ice formation via a mechanism similar to that responsible for H2O ice formation, viz. simple catalytic reactions on grain surfaces.

We present new observations of the FUV (1100-2200 Angstrom) radiation field and the near- to mid-IR (3--13.5 micron) spectral energy distribution (SED) of a sample of T Tauri stars selected on the basis of bright molecular disks (GM Aur, DM Tau, LkCa15). In each source we find evidence for Ly alpha induced H2 fluorescence and an additional source of FUV continuum emission below 1700 Angstroms. Comparison of the FUV spectra to a model of H2 excitation suggests that the strong continuum emission is due to electron impact excitation of H2. The ultimate source of this excitation is likely X-ray irradiation which creates hot photo-electrons mixed in the molecular layer. Analysis of the SED of each object finds the presence of inner disk gaps with sizes of a few AU in each of these young (~1 Myr) stellar systems. We propose that the presence of strong H2 continuum emission and inner disk clearing are related by the increased penetration power of high energy photons in gas rich regions with low grain opacity.

Understanding the chemical evolution in star-forming cores is a necessary pre-condition to correctly assess physical conditions when using molecular emission. We follow the evolution of chemistry and molecular line profiles through the entire star formation process, including a self-consistent treatment of dynamics, dust continuum radiative transfer, gas energetics, chemistry, molecular excitation, and line radiative transfer. In particular, the chemical code follows a gas parcel as it falls toward the center, passing through regimes of density, dust temperature, and gas temperature that are changing both because of the motion of the parcel and the evolving luminosity of the central source. We combine a sequence of Bonnor-Ebert spheres and the inside-out collapse model to describe dynamics from the pre-protostellar stage to later stages. The overall structures of abundance profiles show complex behavior that can be understood as interactions between freeze-out and evaporation of molecules. We find that the presence or absence of gas-phase CO has a tremendous effect on the less abundant species. In addition, the ambient radiation field and the grain properties have important effects on the chemical evolution, and the variations in abundance have strong effects on the predicted emission line profiles. Multi-transition and multi-position ob servations are necessary to constrain the parameters and interpret observations correctly in terms of physical conditions. Good spatial and spectral resolution is also important in distinguishing evolutionary stages.

We develop a method of analyzing radio frequency spectral line observations to derive data on the temperature, density, velocity, and molecular abundance of the emitting gas. The method incorporates a radiative transfer code with a new technique for handling overlapping hyperfine emission lines within the accelerated lambda iteration algorithm and a heuristic search algorithm based on simulated annnealing. We apply this method to new observations of N_2H^+ in three Lynds clouds thought to be starless cores in the first stages of star formation and determine their density structure. A comparison of the gas densities derived from the molecular line emission and the millimeter dust emission suggests that the required dust mass opacity is about kappa_{1.3mm}=0.04 cm^2/g, consistent with models of dust grains that have opacities enhanced by ice mantles and fluffy aggregrates.

We examine the formation of molecular gas behind shocks in atomic gas using a chemical/dynamical model, particular emphasis is given to constraints the chemistry places on the dynamical evolution. The most important result of this study is to stress the importance of shielding the molecular gas from the destructive effects of UV radiation. For shock ram pressures comparable to or exceeding typical local ISM pressures, self-shielding controls the formation time of H2 but CO formation requires shielding of the interstellar radiation field by dust grains. We find that the molecular hydrogen fractional abundance can become significant well before CO forms. The timescale for (CO) molecular cloud formation is not set by H2 formation, but rather by the timescale for accumulating a sufficient column density or extinction, A_V > 0.7. The local ratio of atomic to molecular gas (4:1), coupled with short estimates for cloud lifetimes (3-5 Myr), suggests that the timescales for accumulating molecular clouds from atomic material typically must be no longer than about 12-20 Myr. Based on the shielding requirement, this implies that the typical product of pre-shock density and velocity must be n*v > 20 cm^-3 km s^-1. Based on these results we find that flow-driven formation of molecular clouds in the local interstellar medium can occur sufficiently rapidly to account for observations. We also provide detailed predictions of atomic and molecular emission and absorption that track molecular cloud formation, with a view toward helping to verify cloud formation by shock waves. Finally, we provide an analytic solution for time-dependent H2 formation which may be of use in numerical hydrodynamic calculations.

We have observed 6 clouds along the line-of-sight toward W49A using the Submillimeter Wave Astronomy Satellite (SWAS) and several ground-based observatories. The ortho-H2O 1-0 and OH (1665 and 1667 MHz) transitions are observed in absorption, whereas the low-J CO, 13CO, and C18O lines, as well as the [CI] 1-0 transition, are seen in emission. By using both the o-H218O and o-H2O absorption lines, we are able to constrain the column-averaged o-H_2O abundances in each line-of-sight cloud to within about an order of magnitude. Assuming the standard N(H2)/N(CO) ratio of 10^4, we find N(o-H2O)/N(H2) = 8.1 x 10^-8 - 4 x 10^-7 for three clouds with optically thin water lines. In three additional clouds, the H$_2$O lines are saturated so we have used observations of the H218O ground-state transition to find upper limits to the water abundance of 8.2x 10^-8 - 1.5x10^-6. We measure the OH abundance from the average of the 1665 and 1667 MHz observations and find N(OH)/N(H2) = 2.3x10^-7 - 1.1x10^-6. The o-H2O and OH abundances are similar to those determined for line-of-sight water absorption features towards W51 and Sgr B2 but are higher than those seen from water emission lines in molecular clouds. However, the clouds towards W49 have lower ratios of OH relative to H2O column densities than are predicted by simple models which assume that dissociative recombination is the primary formation pathway for OH and H2O. Building on the work of Neufeld et al. (2002), we present photo-chemistry models including additional chemical effects, which can also explain the observed OH and H2O column densities as well as the observed H2O/CO abundance ratios.

A major limitation and a continuing source of confusion in the interpretation of molecular line observations has been the large degree of chemical complexity that is observed in star-forming molecular cores. The past decade has seen dramatic improvements in our ability to extract chemical information from molecular line observations. These advances have been driven by new observational platforms and techniques that use combined studies of molecular emission and dust emission/absorption to derive chemical abundances. I will outline the various methods and uncertainties in computing chemical information from observations and discuss these new techniques that allow for the estimation of line of sight structure in molecular abundance in both quiescent and active star forming regions. I will show how improved chemical knowledge can be used to place stronger constraints on theoretical astrochemical models, but also open new avenues in the use of molecular emission to examine the physics of molecular clouds and star formation. Thus, we can begin the process of unraveling the chemical complexities and more fully realize the long-standing promise of molecular line observations.

We show in this Letter that the spectral details of the FUV radiation fields have a large impact on the chemistry of protoplanetary disks surrounding T Tauri stars. We show that the strength of a realistic stellar FUV field is significantly lower than typically assumed in chemical calculations and that the radiation field is dominated by strong line emission, most notably Lyman alpha radiation. The effects of the strong Lyman alpha emission on the chemical equilibrium in protoplanetary disks has previously been unrecognized. We discuss the impact of this radiation on molecular observations in the context of a radiative transfer model that includes both direct attenuation and scattering. In particular, Lyman alpha radiation will directly dissociate water vapor and may contribute to the observed enhancements of CN/HCN in disks.

Observations of the 1(10)-1(01) 556.936 GHz transition of ortho-water with the Submillimeter Wave Astronomy Satellite (SWAS) have revealed the presence of widespread emission and absorption by water vapor around the strong submillimeter continuum source Sagittarius B2. An incompletely-sampled spectral line map of a region of size 26 x 19 arcmin around Sgr B2 reveals three noteworthy features. First, absorption by foreground water vapor is detectable at local standard-of-rest (LSR) velocities in the range -100 to 0 km/s at almost every observed position. Second, spatially-extended emission by water is detectable at LSR velocities in the range 80 to 120 km/s at almost every observed position. This emission is attributable to the 180-pc molecular ring identified from previous observations of CO. The typical peak antenna temperature of 0.075 K for this component implies a typical water abundance of 1.2E-6 to 8E-6 relative to H2. Third, strong absorption by water is observed within 5 arcmin of Sgr B2 at LSR velocities in the range 60 to 82 km/s. An analysis of this absorption yields a H2O abundance ~ 2E-7 to 4E-7 relative to H2 if the absorbing water vapor is located within the core of Sgr B2 itself; or, alternatively, a water column density ~ 2.5E+16 to 4E+16 per cm2 if the water absorption originates in the warm, foreground layer of gas proposed previously as the origin of ammonia absorption observed toward Sgr B2.

The past decade has led to significant improvements in our understanding of the physical structure of the molecular cores of cold dark clouds. Observational efforts, in combination with improved knowledge of cloud structure, now provide clear evidence that the chemistry of dark clouds is dominated by the depletion of gaseous species onto grain surfaces. We outline the basis of these observational efforts and show how the abundance determinations have moved beyond single point analyses to the derivation of abundance profiles. We discuss the basic physics of the interaction between molecules and grain surfaces and show that when physics is coupled into a chemical model there is excellent agreement, for a limited set of species, between theory and observations. We discuss our improved understanding of cloud chemistry can be used as a new tool for studies of the formation of stars and planetary systems.

Abundances for C18O, CS, NH3, H2CO, C3H2, and N2H+ and upper limits for the abundances of 13CO and HCO+ are derived for gas within the Bok globule B68 using our own BIMA array data, single-dish data from the literature, and Monte Carlo radiative transfer models. B68 has had its density structure well determined, removing a major uncertainty from abundance determinations. All abundances for B68 are lower than those derived for translucent and cold dense clouds, but perhaps only significantly for N2H+, NH3, and C3H2. CS depletion toward the extinction peak of B68 is hinted at by the large offset between the extinction peak and the position of maximum CS line brightness. C18O and N2H+ abundances derived here are consistent with recently determined values at positions observed in common.

We report sensitive, high resolution molecular-line observations of the dark cloud Barnard 68 obtained with the IRAM 30-m telescope. We analyze spectral-line observations of C18O, CS(2--1), C34S(2--1), and N2H+(1--0) in order to investigate the kinematics and dynamical state of the cloud. We find extremely narrow linewidths in the central regions of the cloud. These narrow lines are consistent with thermally broadened profiles for the measured gas temperature of 10.5 K. We determine the thermal pressure to be a factor 4 -- 5 times greater than the non-thermal (turbulent) pressure in the central regions of the cloud, indicating that thermal pressure is the primary source of support against gravity in this cloud. This confirms the inference of a thermally supported cloud drawn previously from deep infrared extinction measurements. The rotational kinetic energy is found to be only a few percent of the gravitational potential energy, indicating that the contribution of rotation to the overall stability of the cloud is insignificant. Finally, our observations show that CS line is optically thick and self-reversed across nearly the entire projected surface of the cloud. The shapes of the self-reversed profiles are asymmetric and are found to vary across the cloud in such a manner that the presence of both inward and outward motions are observed within the cloud. Moreover, these motions appear to be globally organized in a clear and systematic alternating spatial pattern which is suggestive of a small amplitude, non-radial oscillation or pulsation of the outer layers of the cloud about an equilibrium configuration.

We present SWAS observations of water vapor in two cold star-less clouds, B68 and Core D in rho Ophiuchus. Sensitive non-detections of the 1(10)-1(01) transition of o-H2O are reported for each source. Both molecular cores have been previously examined by detailed observations that have characterized the physical structure. Using these rather well defined physical properties and a Monte-Carlo radiation transfer model we have removed one of the largest uncertainties from the abundance calculation and set the lowest water abundance limit to date in cold low-mass molecular cores. These limits are < 3 x 10^{-8} (relative to H2) and < 8 x 10^{-9} in B68 and rho Oph D, respectively. Such low abundances confirm the general lack of ortho-water vapor in cold (T < 20 K) cores. Provided that the ortho/para ratio of water is not near zero, these limits are well below theoretical predictions and appear to support the suggestion that most of the water in dense low-mass cores is frozen onto the surfaces of cold dust grains.

Using NASA's Submillimeter Wave Astronomy Satellite (SWAS) we have examined the production of water in quiescent and shocked molecular gas through a survey of the 556.936 GHz transition of ortho-H2O in the NGC 1333 molecular core. These observations reveal broad emission lines associated with the IRAS~2, IRAS~4, IRAS~7, and HH7-11 outflows. Towards 3 positions we detect narrow (~2-3 km/s) emission lines clearly associated with the ambient gas. The SWAS observations, with a resolution of 4', are supplemented with observations from the Infrared Space Observatory (ISO) and by an unbiased survey of a 17' x 15' area, with 50" resolution, in the low-J transitions of CO, 13CO, C18O, N2H+, CH3OH, and SiO. Using these combined data sets, with consistent assumptions, we find beam-averaged o-H2O abundances of > 10^{-6} relative to H2O for all four outflows. A comparison of SWAS and ISO water data is consistent with non-dissociative shock models, provided the majority of the 557 GHz emission arises from cool post-shock material with enhanced abundances. In the ambient gas the o-H2O abundance is found to lie between 0.1-1 x 10^{-7} relative to H2 and is enhanced when compared to cold pre-stellar molecular cores. A comparison of the water emission with tracers of dense condensations and shock chemistry finds no clear correlation. However, the water emission appears to be associated with the presence of luminous external heating sources which power the reflection nebula and the photodissociation (PDR) region. Simple PDR models are capable of reproducing the water and high-J 13CO emission, suggesting that a PDR may account for the excitation of water in low density undepleted gas as suggested by Spaans & van Dishoeck (2001).

Abundances for 5 molecules (C18O, CS, NH3, H2CO, and C3H2) and 1 molecular ion (N2H+) and upper limits for the abundances of 1 molecule (13CO) and 1 molecular ion (HCO+) are derived for gas within the Bok globule Barnard 68 (B68). The abundances were determined using our own BIMA millimeter interferometer data and single-dish data gathered from the literature, in conjunction with a Monte Carlo radiative transfer model. Since B68 is the only starless core to have its density structure strongly constrained via extinction mapping, a major uncertainty has been removed from these determinations. All abundances for B68 are lower than those derived for translucent and cold dense clouds, but perhaps only significantly for N2H+, NH3, and C3H2. Depletion of CS toward the extinction peak of B68 is hinted at by the large offset between the extinction peak and the position of maximum CS line brightness. Abundances derived here for C18O and N2H+ are consistent with other, recently determined values at positions observed in common.

In order to obtain a census of the pre-stellar and star-forming molecular cores, we have begun an unbiased survey in CS and N2H+ of the L1630 and L1641 molecular clouds. The use of these two molecular species enables us to quantify and disentangle the effects of depletion often seen in CS observations of dense cores. The spectral line data will provide essential kinematical information not present in similar studies of the sub-millimeter dust-continuum, enabling us to examine the overall core to core velocity dispersion and study the effects of infall and outflows around known sub-mm and infra-red sources. Here we present our initial observations of part of L1630, taken during the commissioning phase of the FCRAO On-the-Fly Mapping system in January 2002.

We present sensitive, high angular resolution molecular-line observations of C18O and N2H+ toward the dark globule B68. We directly compare these data with the near-infrared extinction measurements of Alves, Lada, & Lada (2001) to derive the first evidence for the depletion of N2H+, and by inference N2, in a pre-stellar dark cloud. We also find widespread C18O depletion throughout the centrally condensed core of the B68 cloud. Specifically, we find the N2H+ emission to peak in a shell partially surrounding the peak of dust extinction. Moreover, the N2H+ peaks inside the much larger C18O depletion hole and has a smaller depletion zone, confirming theoretical predictions. These data are analyzed through a direct coupling of time dependent chemical models to a radiation transfer code. This analysis highlights the importance of photodissociation at cloud edges and suggests that the CO abundance declines by two orders of magnitude from edge to center. In contrast N2H+ declines in abundance, at minimum, by at least a factor of two. Indeed it is entirely possible that both N2H+ and N2 are completely absent from the central regions of the B68 core. The depletion of N2H+, and its parent molecule N2, opens the possibility that the centers of dense cores, prior to the formation of a star, may evade detection by conventional methods of probing cores using molecular emission. Under these conditions H2D+ may be the sole viable molecular probe of the innermost regions of star forming cores.

We report the first detection of thermal water vapor emission in the 557 GHz, $1_{10} - 1_{01}$ ground state transition of ortho-H$_2$O toward VY Canis Majoris. In observations obtained with the Submillimeter Wave Astronomy Satellite (SWAS), we measured a flux of $\sim 450 $Jy, in a spectrally resolved line centered on a velocity $v_{LSR} = 25 $km s$^{-1}$ with a full width half maximum of $\sim 35$ km s$^{-1}$, somewhat dependent on the assumed line shape. We analyze the line shape in the context of three different radial outflow models for which we provide analytical expressions. We also detected a weaker 557 GHz emission line from W Hydrae. We find that these and other H$_2$O emission line strengths scale as suggested by Zubko and Elitzur (2000).

We examine the ability of stellar X-rays to desorb molecules from grains in outer protoplanetary disks. In particular, we consider the possibility of spot heating by X-rays and examine its effectiveness, compared to whole grain heating by X-rays, in sustaining a gas phase abundance of CO. As found in previous studies, whole grain heating is effective only for small grains (< 500 A). As a result, large grains are a permanent sink of CO and whole grain heating cannot sustain an equilibrium gas phase abundance of CO. Spot heating, in which the incident X-ray deposits energy in only a restricted region of the grain volume and which may occur as a consequence of the aggregate nature of grains, proves to be more promising. Assuming that grains are comprised of many thermal subunits that are poorly connected thermally, we find that spot heating is efficient at all grain sizes if large grains are effectively ``coated'' with small grains < 100 A in size. We discuss the implications of X-ray spot heating for the interpretation of millimeter emission line studies of outer protoplanetary disks.

Observations of both star-forming regions and young, gas-free stellar associations indicate that most nearby molecular clouds form stars only over a short time span before dispersal; large-scale flows in the diffuse interstellar medium have the potential for forming clouds sufficiently rapidly, and for producing stellar populations with ages much less than the lateral crossing times of their host molecular clouds. We identify four important factors for understanding rapid star formation and short cloud lifetimes. First, much of the accumulation and dispersal of clouds near the solar circle might occur in the atomic phase; only the high-density portion of a cloud's lifecycle is spent in the molecular phase, thus helping to limit molecular cloud ``lifetimes''. Second, once a cloud achieves a high enough column density to form $\h2$ and CO, gravitational forces become larger than typical interstellar pressure forces; thus star formation can follow rapidly upon molecular gas formation and turbulent dissipation in limited areas of each cloud complex. Third, typical magnetic fields are not strong enough to prevent rapid cloud formation and gravitational collapse. Fourth, rapid dispersal of gas by newly-formed stars, and reduction of shielding by a small expansion of the cloud after the first events of star formation, might limit the length of the star formation epoch and the lifetime of a cloud in its molecular state. This picture emphasizes the importance of large-scale boundary conditions for understanding molecular cloud formation, and implies that star formation is a highly dynamic, rather than quasi-static, process.

We present a combined near-infrared and molecular-line study of 25' x 8' area in the Northern streamer of the IC 5146 cloud. Using the technique pioneered by Lada et al 1994, we construct a Gaussian smoothed map of the infrared extinction with the same resolution as the molecular line observations in order to examine correlations of integrated intensities and molecular abundances with extinction for C17O, C34S, and N2H+. We find that over a visual extinction range of 0 to 40 magnitudes, there is good evidence for the presence of differential gas-phase depletions in the densest portions of IC 5146. Both CO and CS exhibit a statistically significant (factor of ~3) abundance reduction near Av ~ 12 magnitudes while, in direct contrast, at the highest extinctions, Av > 10 magnitudes, N2H+ appears relatively undepleted. Moreover, for Av < 4 magnitudes there exists little or no N2H+. This pattern of depletions is consistent with the predictions of chemical theory.

We present a 17-point map of the M17SW cloud core in the 1_{10}-1_{01} transition of ortho-water at 557 GHz obtained with the Submillimeter Wave Astronomy Satellite. Water emission was detected in 11 of the 17 observed positions. The line widths of the water emission vary between 4 and 9 km s^{-1}, and are similar to other emission lines that arise in the M17SW core. A direct comparison is made between the spatial extent of the water emission and the ^{13}CO J = 5\to4 emission; the good agreement suggests that the water emission arises in the same warm, dense gas as the ^{13}CO emission. A spectrum of the H_2^{18}O line was also obtained at the center position of the cloud core, but no emission was detected. We estimate that the average abundance of ortho-water relative to H_2 within the M17 dense core is approximately 1x10^{-9}, 30 times smaller than the average for the Orion core. Toward the H II region/molecular cloud interface in M17SW the ortho-water abundance may be about 5 times larger than in the dense core.