We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $\nu_1$ and $\nu_3$ vibrational bands at 2.7 $\mu$m, from the $\nu_2$ vibrational band at 6.1 $\mu$m, and from pure rotational transitions near 10-13 $\mu$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission.

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.

We report the detection of emission in the v=1-0 P(1) (3.51629 micron) and P(2) (3.60776 micron) rovibrational lines of the helium hydride cation (HeH+) from the planetary nebula NGC 7027. These detections were obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. The confirm the discovery of HeH+ reported recently by Guesten et al. (2019), who used the GREAT instrument on the SOFIA airborne observatory to observe its pure rotational J=1-0 transition at 149.137 micron. The flux measured for the HeH+ v=1-0 P(1) line is in good agreement with our model for the formation, destruction and excitation of HeH+ in NGC 7027. The measured strength of the J=1-0 pure rotational line, however, exceeds the model prediction significantly, as does that of the v=1-0 P(2) line, by factors of 2.9 and 2.3 respectively. Possible causes of these discrepancies are discussed. Our observations of NGC 7027, covering the 3.26 - 3.93 micron spectral region, have led to the detection of more than sixty spectral lines including nine rovibrational emissions from CH+. The latter are detected for the first time in an astronomical source.

Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters.

We report the results from a new, highly sensitive ($\Delta T_{mb} \sim 3 $mK) survey for thermal OH emission at 1665 and 1667 MHz over a dense, 9 x 9-pixel grid covering a $1\deg$ x $1\deg$ patch of sky in the direction of $l = 105\deg, b = +2.50\deg$ towards the Perseus spiral arm of our Galaxy. We compare our Green Bank Telescope (GBT) 1667 MHz OH results with archival CO J=1-0 observations from the Five College Radio Astronomy Observatory (FCRAO) Outer Galaxy Survey within the velocity range of the Perseus Arm at these galactic coordinates. Out of the 81 statistically-independent pointings in our survey area, 86% show detectable OH emission at 1667 MHz, and 19% of them show detectable CO emission. We explore the possible physical conditions of the observed features using a set of diffuse molecular cloud models. In the context of these models, both OH and CO disappear at current sensitivity limits below an A$_{\rm v}$ of 0.2, but the CO emission does not appear until the volume density exceeds 100-200 cm$^{-3}$. These results demonstrate that a combination of low column density A$_{\rm v}$ and low volume density $n_{H}$ can explain the lack of CO emission along sight lines exhibiting OH emission. The 18-cm OH main lines, with their low critical density of $n^{*}$ $ \sim 1 $ cm$^{-3}$, are collisionally excited over a large fraction of the quiescent galactic environment and, for observations of sufficient sensitivity, provide an optically-thin radio tracer for diffuse H$_2$.

Spectrally-resolved observations of three pure rotational lines of H$_2$, conducted with the EXES instrument on SOFIA toward the classic bow shock HH7, reveal systematic velocity shifts between the S(5) line of ortho-H$_2$ and the two para-H$_2$ lines [S(4) and S(6)] lying immediately above and below it on the rotational ladder. These shifts, reported here for the first time, imply that we are witnessing the conversion of para-H$_2$ to ortho-H$_2$ within a shock wave driven by an outflow from a young stellar object. The observations are in good agreement with the predictions of models for non-dissociative, C-type molecular shocks. They provide a clear demonstration of the chemical changes wrought by interstellar shock waves, in this case the conversion of para-H$_2$ to ortho-H$_2$ in reactive collisions with atomic hydrogen, and provide among the most compelling evidence yet obtained for C-type shocks in which the flow velocity changes continuously.

During the dawn of chemistry when the temperature of the young Universe had fallen below $\sim$4000 K, the ions of the light elements produced in Big Bang nucleosynthesis recombined in reverse order of their ionization potential. With its higher ionization potentials, He$^{++}$ (54.5 eV) and He$^+$ (24.6 eV) combined first with free electrons to form the first neutral atom, prior to the recombination of hydrogen (13.6 eV). At that time, in this metal-free and low-density environment, neutral helium atoms formed the Universe's first molecular bond in the helium hydride ion HeH$^+$, by radiative association with protons (He + H$^+$ $\rightarrow$ HeH$^+$ + h$\nu$). As recombination progressed, the destruction of HeH$^+$ (HeH$^+$ + H $\rightarrow$ He + H$_2^+$) created a first path to the formation of molecular hydrogen, marking the beginning of the Molecular Age. Despite its unquestioned importance for the evolution of the early Universe, the HeH$^+$ molecule has so far escaped unequivocal detection in interstellar space. In the laboratory, the ion was discovered as long ago as 1925, but only in the late seventies was the possibility that HeH$^+$ might exist in local astrophysical plasmas discussed. In particular, the conditions in planetary nebulae were shown to be suitable for the production of potentially detectable HeH$^+$ column densities: the hard radiation field from the central hot white dwarf creates overlapping Str\"omgren spheres, where HeH$^+$ is predicted to form, primarily by radiative association of He$^+$ and H. With the GREAT spectrometer onboard SOFIA, the HeH$^+$ rotational ground-state transition at $\lambda$149.1 $\mu$m is now accessible. We report here its detection towards the planetary nebula NGC7027.

To constrain the properties of dark matter (DM) that interacts with nucleons, we have conducted an experimental search for any anomalous heating of ordinary baryonic matter at 77 K. Our tabletop experiment is motivated by the possibility (discussed in a previous paper) that DM particles with masses in the $\sim 1 - 2 m_{\rm p}$ range could be captured by and concentrated within the Earth. For suitable parameters, this phenomenon could lead to a substantial density $\sim 10^{14}\,\rm cm^{-3}$ of thermalized (300 K) DM particles at Earth's surface that would heat cooler baryonic matter. Our experiment involves precise differential measurements of the evaporation rate of liquid nitrogen in a storage dewar within which various materials are immersed. The results revealed no statistically-significant detections of heating in the 27 elements with molar fractions $> 10^{-5}$ in Earth's crust. For material with the average composition of Earth's crust, our measurements imply a $3 \sigma$ upper limit of $1.32 \times 10^{-27}\, n_{\rm 14}^{-1} (m_{\rm DM}/2\,m_{\rm p})^{-1/2}\,\rm cm^2$ on the mean cross-section for scattering with thermal HIDM at 300 K, where $10^{14}\,n_{14} \, \rm cm^{-3}$ is the particle density at Earth's surface. In combination with a lower limit on the scattering cross-section, obtained from a consideration of the heat flow through the Earth's crust (Neufeld, Farrar & McKee 2018), our experiment places an upper limit of $1.6 \times 10^{13}\,\rm cm^{-3}$ on the density of DM at the Earth's surface. This in turn, significantly constrains the properties for any DM candidate that interacts with baryons.

The transition between atomic and molecular hydrogen is associated with important changes in the structure of interstellar clouds, and marks the beginning of interstellar chemistry. Because of the relatively simple networks controlling their abundances, molecular ions are usually good probes of the underlying physical conditions including for instance the fraction of gas in molecular form or the fractional ionization. In this paper we focus on three possible probes of the molecular hydrogen column density, HCO+, HOC+, and CF+. We presented high sensitivity ALMA absorption data toward a sample of compact HII regions and bright QSOs with prominent foreground absorption, in the ground state transitions of the molecular ions HCO+, HOC+, and CF+ and the neutral species HCN and HNC, and from the excited state transitions of C3H+(4-3) and 13CS(2-1). These data are compared with Herschel absorption spectra of the ground state transition of HF and p-H2O. We show that the HCO+, HOC+, and CF+ column densities are well correlated with each other. HCO+ and HOC+ are tightly correlated with p-H2O, while they exhibit a different correlation pattern with HF depending on whether the absorbing matter is located in the Galactic disk or in the central molecular zone. We report new detections of C3H+ confirming that this ion is ubiquitous in the diffuse matter, with an abundance relative to H2 of ~7E-11. We confirm that the CF+ abundance is lower than predicted by simple chemical models and propose that the rate of the main formation reaction is lower by a factor of about 3 than usually assumed. In the absence of CH or HF data, we recommend to use the ground state transitions of HCO+, CCH, and HOC+ to trace diffuse molecular hydrogen, with mean abundances relative to H2 of 3E-9, 4E-8 and 4E-11.

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $\nu_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO$_2$ gas in the massive YSO MonR2 IRS3. This yields a SO$_2$/$\mathrm{H}$ abundance lower limit of $5.6\pm0.5\times10^{-7}$, or $>\!4\%$ of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of $b<3.20\;\mathrm{km\;s}^{-1}$. The small line widths and high temperature ($T_\mathrm{ex}=234\pm15\;\mathrm{K}$) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov--Gerasimenko), which is followed by SO$_2$ formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO$_2$ from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

We have performed a 5-8 $\mu$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $\mu$m using the iSHELL instrument on the Infrared Telescope Facility (IRTF), and the full N band window from 8-13 $\mu$m using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the IRTF. Here we present the first detection of ro-vibrational transitions of CS in this source. The absorption lines are centred on average around -10 kms$^{-1}$ and the line widths of CS compare well with the hot component of $^{13}$CO (around 10 kms$^{-1}$). Temperatures for CS, hot $^{13}$CO and $^{12}$CO v=1-2 agree well and are around 700 K. We derive a CS abundance of 8$\times$10$^{-3}$ and 2$\times$10$^{-6}$ with respect to CO and H$_2$ respectively. This enhanced CS abundance with respect to the surrounding cloud (1$\times$10$^{-8}$) may reflect sublimation of H$_2$S ice followed by gas-phase reactions to form CS. Transitions are in LTE and we derive a density of $>$10$^7$ cm$^{-3}$, which corresponds to an absorbing region of $<$0.04$''$. EXES observations of CS are likely to probe deeply into the hot core, to the base of the outflow. Submillimeter and infrared observations trace different components of the hot core as revealed by the difference in systemic velocities, line widths and temperatures, as well as the CS abundance.

Orion BN has a large proper motion and radial velocity with respect to the gas and other stars in the region where it is presumed to have formed. Multiple dynamical interaction scenarios have been proposed to explain this motion. In one case BN is thought to have interacted with stars in the Trapezium cluster, while in another it is thought to have interacted with source I while deeply embedded in molecular gas. If there is dense gas that has been retained in close proximity to BN, it may be evidence that the latter scenario is favored. We observed BN at high spectral resolution in three windows near 6 $\mu$m using SOFIA/EXES targeting the $\nu_2$ vibrational band of H$_2$O. Absorption from only three transitions of H$_2$O is detected, and through kinematic analysis is associated with cool, dense foreground gas, not BN itself. We find no evidence for H$_2$O absorption or emission at the systemic velocity of BN.

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.

For dark matter (DM) particles with masses in the 0.6 - 6 m_p range, we set stringent constraints on the interaction cross-sections for scattering with ordinary baryonic matter. These constraints follow from the recognition that such particles can be captured by - and thermalized within - the Earth, leading to a substantial accumulation and concentration of DM that interact with baryons. Here, we discuss the probability that DM intercepted by the Earth will be captured, the number of DM particles thereby accumulated over Earth's lifetime, the fraction of such particles retained in the face of evaporation, and the density distribution of such particles within the Earth. In the latter context, we note that a previous treatment of the density distribution of DM, presented by Gould and Raffelt and applied subsequently to DM in the Sun, is inconsistent with considerations of hydrostatic equilibrium. Our analysis provides an estimate of the DM particle density at Earth's surface, which may exceed 1.E+14 cm-3 for the mass range under consideration. Based upon our determination of the DM density at Earth's surface, we derive constraints on the scattering cross-sections. These constraints are placed by four considerations: (1) the lifetime of the relativistic proton beam at the Large Hadron collider (LHC); (2) the orbital decay of spacecraft in low Earth orbit (LEO); (3) the vaporization rate of cryogenic liquids in well-insulated storage dewars; and (4) the thermal conductivity of Earth's crust. As an example application of our results, we show that for the scattering cross-sections that were invoked recently in Barkana's original explanation for the anomalously deep 21 cm absorption reported by EDGES, DM particle masses in the 0.6 - 4 m_p range are ruled out.

We report the discovery of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz 8(27)-7(34) transition of water toward three oxygen-rich evolved stars: W Hya, U Her, and VY CMa. An upper limit on the 1.296 THz line flux was obtained toward R Aql. Near-simultaneous observations of the 22.23508 GHz 6(16)-5(23) water maser transition were carried out towards all four sources using the Effelsberg 100m telescope. The measured line fluxes imply 22 GHz / 1.296 THz photon luminosity ratios of 0.012, 0.12, and 0.83 respectively for W Hya, U Her, and VY CMa, values that confirm the 22 GHz maser transition to be unsaturated in W Hya and U Her. We also detected the 1.884888 THz 8(45)-7(53) transition toward W Hya and VY CMa, and the 1.278266 THz 7(43)-6(52) transition toward VY CMa. Like the 22 GHz maser transition, all three of the THz emission lines detected here originate from the ortho-H2O spin isomer. Based upon a model for the circumstellar envelope of W Hya, we estimate that stimulated emission is responsible for ~ 85% of the observed 1.296 THz line emission, and thus that this transition may be properly described as a terahertz-frequency maser. In the case of the 1.885 THz transition, by contrast, our W Hya model indicates that the observed emission is dominated by spontaneous radiative decay, even though a population inversion exists.

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.

We have obtained estimates for the cosmic-ray ionization rate (CRIR) in the Galactic disk, using a detailed model for the physics and chemistry of diffuse interstellar gas clouds to interpret previously-published measurements of the abundance of four molecular ions: ArH$^+$, OH$^+$, H$_2$O$^+$ and H$_3^+$. For diffuse $atomic$ clouds at Galactocentric distances in the range $R_g \sim 4 - 9$ kpc, observations of ArH$^+$, OH$^+$, and H$_2$O$^+$ imply a mean primary CRIR of $(2.2 \pm 0.3) \exp [(R_0-R_g)/4.7\,\rm{kpc}] \times 10^{-16} \rm \, s^{-1}$ per hydrogen atom, where $R_0=8.5$ kpc. Within diffuse $molecular$ clouds observed toward stars in the solar neighborhood, measurements of H$_3^+$ and H$_2$ imply a primary CRIR of $(2.3 \pm 0.6) \times 10^{-16}\,\,\rm s^{-1}$ per H atom, corresponding to a total ionization rate per H$_2$ molecule of $(5.3 \pm 1.1) \times 10^{-16}\,\,\rm s^{-1},$ in good accord with previous estimates. These estimates are also in good agreement with a rederivation, presented here, of the CRIR implied by recent observations of carbon and hydrogen radio recombination lines along the sight-line to Cas A. Here, our best-fit estimate for the primary CRIR is $2.9 \times 10^{-16}\,\,\rm s^{-1}$ per H atom. Our results show marginal evidence that the CRIR in diffuse molecular clouds decreases with cloud extinction, $A_{\rm V}({\rm tot})$, with a best-fit dependence $\propto A_{\rm V}({\rm tot})^{-1}$ for $A_{\rm V}({\rm tot}) \ge 0.5$.

Recent submillimeter and far-infrared wavelength observations of absorption in the rotational ground-state lines of various simple molecules against distant Galactic continuum sources have opened the possibility of studying the chemistry of diffuse molecular clouds throughout the Milky Way. In order to calculate abundances, the column densities of molecular and atomic hydrogen, HI, must be known. We aim at determining the atomic hydrogen column densities for diffuse clouds located on the sight lines toward a sample of prominent high-mass star-forming regions that were intensely studied with the HIFI instrument onboard Herschel. Based on Jansky Very Large Array data, we employ the 21 cm HI absorption-line technique to construct profiles of the HI opacity versus radial velocity toward our target sources. These profiles are combined with lower resolution archival data of extended HI emission to calculate the HI column densities of the individual clouds along the sight lines. We employ Bayesian inference to estimate the uncertainties of the derived quantities. Our study delivers reliable estimates of the atomic hydrogen column density for a large number of diffuse molecular clouds at various Galactocentric distances. Together with column densities of molecular hydrogen derived from its surrogates observed with HIFI, the measurements can be used to characterize the clouds and investigate the dependence of their chemistry on the molecular fraction, for example.

We report the results of a sensitive search for the 443.952902 GHz $J=1-0$ transition of the LiH molecule toward two interstellar clouds in the Milky Way, W49N and Sgr B2 (Main), that has been carried out using the Atacama Pathfinder Experiment (APEX) telescope. The results obtained toward W49N place an upper limit of $1.9 \times 10^{-11}\, (3\sigma)$ on the LiH abundance, $N({\rm LiH})/N({\rm H}_2)$, in a foreground, diffuse molecular cloud along the sight-line to W49N, corresponding to 0.5% of the solar system lithium abundance. Those obtained toward Sgr B2 (Main) place an abundance limit $N({\rm LiH})/N({\rm H}_2) < 3.6 \times 10^{-13} \,(3\sigma)$ in the dense gas within the Sgr B2 cloud itself. These limits are considerably smaller that those implied by the tentative detection of LiH reported previously for the $z=0.685$ absorber toward B0218+357.

We present a general parameter study, in which the abundance of interstellar argonium (ArH$^+$) is predicted using a model for the physics and chemistry of diffuse interstellar gas clouds. Results have been obtained as a function of UV radiation field, cosmic-ray ionization rate, and cloud extinction. No single set of cloud parameters provides an acceptable fit to the typical ArH$^+$, OH$^+$ and $\rm H_2O^+$ abundances observed in diffuse clouds within the Galactic disk. Instead, the observed abundances suggest that ArH$^+$ resides primarily in a separate population of small clouds of total visual extinction of at most 0.02 mag per cloud, within which the column-averaged molecular fraction is in the range $10^{-5} - 10^{-2}$, while OH$^+$ and $\rm H_2O^+$ reside primarily in somewhat larger clouds with a column-averaged molecular fraction $\sim 0.2$. This analysis confirms our previous suggestion that the argonium molecular ion is a unique tracer of almost purely atomic gas.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level).

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. The chemistry of the diffuse interstellar medium rests upon three pillars: exothermic ion-neutral reactions (" cold chemistry "), endothermic neutral-neutral reactions with significant activation barriers (" warm chemistry "), and reactions on the surfaces of dust grains. While warm chemistry becomes important in the shocks associated with turbulent dissipation regions, the main path for the formation of interstellar OH and H2O is that of cold chemistry. Aims. The aim of this study is to observationally confirm the association of atomic oxygen with both atomic and molecular gas phases, and to understand the measured abundances of OH and OH + as a function of the available reservoir of H2. Methods. We obtained absorption spectra of the ground states of OH, OH+ and OI with high-velocity resolution, with GREAT on-board SOFIA, and with the THz receiver at the APEX. We analyzed them along with ancillary spectra of HF and CH from HIFI. To deconvolve them from the hyperfine structure and to separate the blend that is due to various velocity components on the sightline, we fit model spectra consisting of an appropriate number of Gaussian profiles using a method combining simulated annealing with downhill simplex minimization. Together with HF and/or CH as a surrogate for H2, and HI $\lambda$21 cm data, the molecular hydrogen fraction f^N\_H2 = N(H 2)/(N(H) + 2N(H 2)) can be determined. We then investigated abundance ratios as a function of f^N\_H2. Results. The column density of OI is correlated at a high significance with the amount of available molecular and atomic hydrogen, with an atomic oxygen abundance of $3 \times 10 ^{-4}$ relative to H nuclei. While the velocities of the absorption features of OH and OH+ are loosely correlated and reflect the spiral arm crossings on the sightline, upon closer inspection they display an anticorrespondence. The arm-to-interarm density contrast is found to be higher in OH than in OH+. While both species can coexist, with a higher abundance in OH than in OH+, the latter is found less frequently in absence of OH than the other way around, which is a direct consequence of the rapid destruction of OH+ by dissociative recombination when not enough H2 is available. This conjecture has been substantiated by a comparison between the OH/OH+ ratio with f^N\_H2, showing a clear correlation. The hydrogen abstraction reaction chain OH+ (H2,H) H2O+ (H2,H)H3O+ is confirmed as the pathway for the production of OH and H 2 O. Our estimate of the branching ratio of the dissociative recombination of H3O+ to OH and H2O is confined within the interval of 84 to 91%, which matches laboratory measurements (74 to 83%). -- A correlation between the linewidths and column densities of OH+ features is found to be significant with a false-alarm probability below 5%. Such a correlation is predicted by models of interstellar MHD turbulence. For OH the same correlation is found to be insignificant because there are more narrow absorption features. Conclusions. While it is difficult to assess the contributions of warm neutral-neutral chemistry to the observed abundances, it seems fair to conclude that the predictions of cold ion-neutral chemistry match the abundance patterns we observed.

Context. The recent detection of warm H$_2$O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

Interstellar hydrides -- that is, molecules containing a single heavy element atom with one or more hydrogen atoms -- were among the first molecules detected outside the solar system. They lie at the root of interstellar chemistry, being among the first species to form in initially-atomic gas, along with molecular hydrogen and its associated ions. Because the chemical pathways leading to the formation of interstellar hydrides are relatively simple, the analysis of the observed abundances is relatively straightforward and provides key information about the environments where hydrides are found. Recent years have seen rapid progress in our understanding of interstellar hydrides, thanks largely to far-IR and submillimeter observations performed with the Herschel Space Observatory. In this review, we will discuss observations of interstellar hydrides, along with the advanced modeling approaches that have been used to interpret them, and the unique information that has thereby been obtained.

We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass-flow rates, the different classes of young stellar objects lie in the sequence Class 0 -- Class I/flat-spectrum -- Class II, indicating that the trend is an evolutionary sequence in which mass outflow and accretion rates decrease together with increasing age, while maintaining rough proportionality. The survey results include two which are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b, and limits on the distribution of outflow speeds. Neither rule out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (0.31) of the sample have branching ratio sufficiently small that only disk winds, launched on scales as large as several AU, have been demonstrated to account for them.

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.

To study the usefulness of \CFP\ as a tracer of the regions where C\p\ and \HH\ coexist in the interstellar medium. We used the Plateau de Bure Interferometer to synthesize \CFP\ J=1-0 absorption at 102.6 GHz toward the core of the distant HII region W49N at l = 43.2\degr, b=0.0\degr, and we modeled the fluorine chemistry in diffuse/translucent molecular gas. We detected \CFP\ absorption over a broad range of velocity showing that \CFP\ is widespread in the \HH-bearing Galactic disk gas. Originally detected in dense gas in the Orion Bar and Horsehead PDR, \CFP\ was subsequently detected in absorption from diffuse and translucent clouds seen toward \bll\ and 3C111. Here we showed that \CFP\ is distributed throughout the diffuse and translucent molecular disk gas with N(\CFP)/N(\HH) $= 1.5-2.0\times10^{-10}$, increasing to N(\CFP)/N(\HH) $= 3.5\times10^{-10}$ in one cloud at 39 \kms\ having higher N(\HH) $\approx 3\times10^{21}\pcc$. Models of the fluorine chemistry reproduce the observed column densities and relative abundance of HF, from which \CFP\ forms, but generally overpredict the the column density of \CFP\ by factors of 1.4-4. We show that a free space photodissociation rate $\Gamma \ga 10^{-9}\ps$, comparable to that of CH, might account for much of the discrepancy but a recent calculation finds a value about ten times smaller. In the heavily blended and kinematically complex spectra seen toward W49, \CFP\ absorption primarily traces the peaks of the \HH\ distribution.

We combine Herschel observations of a total of 12 sources to construct the most uniform survey of HF and H2O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-Infrared (HIFI) allows us to compare the HF and H2O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H2O velocity distributions follow each other almost perfectly and establish that HF and H2O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H2 in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H2 in our sample, we find that the N(H2O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H2O might be used as a tracer of H2 -within a factor 2.5- in the diffuse interstellar medium. We show that the measured factor of ~2.5 variation around the median is driven by true local variations in the H2O abundance relative to H2 throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H2O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse interstellar medium n_H < 1000 cm^-3.

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

We present high spectral resolution (~3 km/s) observations of the nu_2 ro-vibrational band of H2O in the 6.086--6.135 micron range toward the massive protostar AFGL 2591 using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). Ten absorption features are detected in total, with seven caused by transitions in the nu_2 band of H2O, two by transitions in the first vibrationally excited nu_2 band of H2O, and one by a transition in the nu_2 band of H2{18}O. Among the detected transitions is the nu_2 1(1,1)--0(0,0) line which probes the lowest lying rotational level of para-H2O. The stronger transitions appear to be optically thick, but reach maximum absorption at a depth of about 25%, suggesting that the background source is only partially covered by the absorbing gas, or that the absorption arises within the 6 micron emitting photosphere. Assuming a covering fraction of 25%, the H2O column density and rotational temperature that best fit the observed absorption lines are N(H2O)=(1.3+-0.3)*10^{19} cm^{-2} and T=640+-80 K.

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2\Pi_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3.

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

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

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.

Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries.

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board \textit{Herschel}\thanks{\textit{Herschel} is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} and with the IRAM\thanks{This work was based on observations carried out with the IRAM 30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain)} 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.

Ionized carbon is the main gas-phase reservoir of carbon in the neutral diffuse interstellar medium and its 158 micron fine structure transition [CII] is the most important cooling line of the diffuse interstellar medium (ISM). We combine [CII] absorption and emission spectroscopy to gain an improved understanding of physical conditions in the different phases of the ISM. We present high resolution [CII] spectra obtained with the Herschel/HIFI instrument towards bright dust continuum sources regions in the Galactic plane, probing simultaneously the diffuse gas along the line of sight and the background high-mass star forming regions. These data are complemented by observations of the 492 and 809 GHz fine structure lines of atomic carbon and by medium spectral resolution spectral maps of the fine structure lines of atomic oxygen at 63 and 145 microns with Herschel/PACS. We show that the presence of foreground absorption may completely cancel the emission from the background source in medium spectral resolution data and that high spectral resolution spectra are needed to interpret the [CII] and [OI] emission and the [CII]/FIR ratio. This phenomenon may explain part of the [CII]/FIR deficit seen in external luminous infrared galaxies. The C+ and C excitation in the diffuse gas is consistent with a median pressure of 5900 Kcm-3 for a mean TK ~100 K. The knowledge of the gas density allows us to determine the filling factor of the absorbing gas along the selected lines of sight: the median value is 2.4 %, in good agreement with the CNM properties. The mean excitation temperature is used to derive the average cooling due to C+ in the Galactic plane : 9.5 x 10^{-26} erg/s/H. Along the observed lines of sight, the gas phase carbon abundance does not exhibit a strong gradient as a function of Galacto-centric radius and has a weighted average of C/H = 1.5 +/- 0.4 x 10^{-4}.

(abbreviated) We aim to study the inner wind of the well-known AGB star CW Leo. Different diagnostics probing different geometrical scales have pointed toward a non-homogeneous mass-loss process: dust clumps are observed at milli-arcsec scale, a bipolar structure is seen at arcsecond-scale and multi-concentric shells are detected beyond 1". We present the first ALMA Cycle 0 band 9 data around 650 GHz. The full-resolution data have a spatial resolution of 0".42x0".24, allowing us to study the morpho-kinematical structure within ~6". Results: We have detected 25 molecular lines. The emission of all but one line is spatially resolved. The dust and molecular lines are centered around the continuum peak position. The dust emission has an asymmetric distribution with a central peak flux density of ~2 Jy. The molecular emission lines trace different regions in the wind acceleration region and suggest that the wind velocity increases rapidly from about 5 R* almost reaching the terminal velocity at ~11 R*. The channel maps for the brighter lines show a complex structure; specifically for the 13CO J=6-5 line different arcs are detected within the first few arcseconds. The curved structure present in the PV map of the 13CO J=6-5 line can be explained by a spiral structure in the inner wind, probably induced by a binary companion. From modeling the ALMA data, we deduce that the potential orbital axis for the binary system lies at a position angle of ~10-20 deg to the North-East and that the spiral structure is seen almost edge-on. We infer an orbital period of 55 yr and a binary separation of 25 au (or ~8.2 R*). We tentatively estimate that the companion is an unevolved low-mass main-sequence star. The ALMA data hence provide us for the first time with the crucial kinematical link between the dust clumps seen at milli-arcsecond scale and the almost concentric arcs seen at arcsecond scale.

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

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

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.

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 the detection of hydrogen fluoride, HF, in two luminous nearby galaxies NGC 253 and NGC 4945 using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory. The HF line toward NGC 253 has a P-Cygni profile, while an asymmetric absorption profile is seen toward NGC 4945. The P-Cygni profile in NGC 253 suggests an outflow of molecular gas with a mass of M(H$_2$)$_{out}$ $\sim$ 1 $\times$ 10$^7$ M$_\odot$ and an outflow rate as large as \.{M} $\sim$ 6.4 M$_\odot$ yr$^{-1}$. In the case of NGC 4945, the axisymmetric velocity components in the HF line profile is compatible with the interpretation of a fast-rotating nuclear ring surrounding the nucleus and the presence of inflowing gas. The gas falls into the nucleus with an inflow rate of $\le$ 1.2 M$_\odot$ yr$^{-1}$, inside a inner radius of $\le$ 200 pc. The gas accretion rate to the central AGN is much smaller, suggesting that the inflow can be triggering a nuclear starburst. From these results, the HF $J = 1-0$ line is seen to provide an important probe of the kinematics of absorbing material along the sight-line to nearby galaxies with bright dust continuum and a promising new tracer of molecular gas in high-redshift galaxies.

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

Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts. This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditions. Three major routes to water formation are identified: low temperature ion-molecule chemistry, high-temperature neutral-neutral chemistry and gas-ice chemistry. The rate coefficients of several important processes entering the networks are discussed in detail; several of them have been determined only in the last decade through laboratory experiments and theoretical calculations. Astronomical examples of each of the different chemical routes are presented using data from powerful new telescopes, in particular the Herschel Space Observatory. Basic chemical physics studies remain critically important to analyze astronomical data.

We would like to understand the chemistry of dense clouds and their hot cores more quantitatively by obtaining more complete knowledge of the chemical species present in them. We have obtained high-resolution infrared absorption spectroscopy at 3-4 um toward the bright infrared source CRL 2136. The fundamental vibration-rotation band of HCl has been detected within a dense cloud for the first time. The HCl is probably located in the warm compact circumstellar envelope or disk of CRL 2136. The fractional abundance of HCl is (4.9-8.7)e-8, indicating that approximately 20 % of the elemental chlorine is in gaseous HCl. The kinetic temperature of the absorbing gas is 250 K, half the value determined from infrared spectroscopy of 13CO and water. The percentage of chlorine in HCl is approximately that expected for gas at this temperature. The reason for the difference in temperatures between the various molecular species is unknown.

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

We present ground-based observations of the nu_1 and nu_3 fundamental bands of H2O toward the massive protostar AFGL 2136 IRS 1, identifying absorption features due to 47 different ro-vibrational transitions between 2.468 micron and 2.561 micron. Analysis of these features indicates the absorption arises in warm (T=506+-25 K), very dense (n(H2)>5*10^9 cm^-3) gas, suggesting an origin close to the central protostar. The total column density of warm water is estimated to be N(H2O)=(1.02+-0.02)*10^19 cm^-2, giving a relative abundance of N(H2O)/N(H2)~10^-4. Our study represents the first extensive use of water vapor absorption lines in the near infrared, and demonstrates the utility of such observations in deriving physical parameters.

We present a 52-671um spectral scan toward SgrA* taken with the PACS and SPIRE spectrometers onboard Herschel. The achieved angular resolution allows us to separate, for the first time at far-IR wavelengths, the emission toward the central cavity (gas in the inner central parsec of the galaxy) from that of the surrounding circum-nuclear disk. The spectrum toward SgrA* is dominated by strong [OIII], [OI], [CII], [NIII], [NII], and [CI] fine structure lines (in decreasing order of luminosity) arising in gas irradiated by UV-photons from the central stellar cluster. In addition, rotationally excited lines of 12CO (from J=4-3 to 24-23), 13CO, H2O, OH, H3O+, HCO+ and HCN, as well as ground-state absorption lines of OH+, H2O+, H3O+, CH+, H2O, OH, HF, CH and NH are detected. The excitation of the 12CO ladder is consistent with a hot isothermal component at Tk ~ 10^{3.1} K and n(H2)< 10^4 cm^{-3}. It is also consistent with a distribution of temperature components at higher density with most CO at Tk<300 K. The detected molecular features suggest that, at present, neither very enhanced X-ray, nor cosmic-ray fluxes play a dominant role in the heating of the hot molecular gas. The hot CO component (either the bulk of the CO column or just a small fraction depending on the above scenario) results from a combination of UV- and shock-driven heating. If irradiated dense clumps/clouds do not exist, shocks likely dominate the heating of the hot molecular gas. This is consistent with the high-velocity gas detected toward SgrA*.

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have performed mapping observations of the 620.701 GHz 5(32)-4(41) transition of ortho-H2O within a roughly 1.5 x 1.5 arcmin region encompassing the Kleinmann-Low nebula in Orion, and pointed observations of that transition toward the Orion South condensation and the W49N region of high-mass star formation. Using the Effelsberg 100 m radio telescope, we obtained ancillary observations of the 22.23508 GHz 6(16)-5(23) water maser transition; in the case of Orion-KL, the 621 GHz and 22 GHz observations were carried out within 10 days of each other. The 621 GHz water line emission shows clear evidence for strong maser amplication in all three sources, exhibiting narrow (roughly 1 km/s FWHM) emission features that are coincident (kinematically and/or spatially) with observed 22 GHz features. Moreover, in the case of W49N - for which observations were available at three epochs spanning a two year period - the spectra exhibited variability. The observed 621 GHz/22 GHz line ratios are consistent with a maser pumping model in which the population inversions arise from the combined effects of collisional excitation and spontaneous radiative decay, and the inferred physical conditions can plausibly arise in gas heated by either dissociative or non-dissociative shocks. The collisional excitation model also predicts that the 22 GHz population inversion will be quenched at higher densities than that of the 621 GHz transition, providing a natural explanation for the observational fact that 22 GHz maser emission appears to be a necessary but insufficient condition for 621 GHz maser emission.

We report the detection of hydrogen chloride, HCl, in diffuse molecular clouds on the line of sight towards the star-forming region W31C (G10.6-0.4). The J = 1-0 lines of the two stable HCl isotopologues, H35Cl and H37Cl, are observed using the 1b receiver of the Heterodyne Instrument for the Far-Infrared (HIFI) aboard the Herschel Space Observatory. The HCl line is detected in absorption, over a wide range of velocities associated with diffuse clouds along the line of sight to W31C. The analysis of the absorption strength yields a total HCl column density of few 10^13 cm^-2, implying that HCl accounts for ~0.6 % of the total gasphase chlorine, which exceeds by a factor of ~6 the theoretical model predictions. This result is comparable to those obtained from the chemically-related species H2Cl+ and HCl+, for which large column densities have also been reported on the same line of sight. The source of discrepancy between models and observations is still unknown; however, the detection of these Cl-bearing molecules, provides key constraints for the chlorine chemistry in the diffuse gas.

We report the results of a sensitive search for the minor isotopologues of water, H2-17O and H2-18O, toward the carbon-rich AGB star IRC+10216 (a.k.a. CW Leonis) using the HIFI instrument on the Herschel Space Observatory. This search was motivated by the fact that any detection of isotopic enhancement in the H2-17O and H2-18O abundances would have strongly implicated CO photodissociation as the source of the atomic oxygen needed to produce water in a carbon-rich circumstellar envelope. Our observations place an upper limit of 1/470 on the H2-17O/H2-16O abundance ratio. Given the isotopic 17O/16O ratio of 1/840 inferred previously for the photosphere of IRC+10216, this result places an upper limit of a factor 1.8 on the extent of any isotope-selective enhancement of H2-17O in the circumstellar material, and provides an important constraint on any model that invokes CO photodissociation as the source of O for H2O production. In the context of the clumpy photodissociation model proposed previously for the origin of water in IRC+10216, our limit implies that 12C-16O (not 13C-16O or SiO) must be the dominant source of 16O for H2O production, and that the effects of self-shielding can only have reduced the 12C-16O photodissociation rate by at most a factor ~ 2.

We report the first detection of the v=1--0, R(0) ro-vibrational transition of HF at 2.499385 microns arising from interstellar gas. The line is seen in absorption toward 3 background sources---HD 154368, Elias 29, and AFGL 2136 IRS 1---all of which have reported H2 column densities determined from observations of H2. This allows for the first direct determination of the HF/H2 abundance ratio. We find values of N(HF)/N(H2)=1.15*10^-8 and 0.69*10^-8 for HD 154368 and Elias 29, respectively. The sight line toward AFGL 2136 IRS 1 also shows absorption from the v=1--0, R(1) transition of HF, indicating warm, dense (n_H > 10^9 cm^-3) gas, likely very close to the central protostar. Ascribing portions of the HF absorption to warm and cold gas, we find N(HF)/N(H2)=(1.7--2.9)*10^-8 and (0.33--0.58)*10^-8 for the two components, respectively. Except for the warm component toward AFGL 2136 IRS 1, all observed HF/H2 ratios are well below N(HF)/N(H2)=3.6*10^-8, the value predicted if nearly all gas-phase fluorine is in the form of HF. Models of fluorine chemistry that account for depletion onto grains are able to reproduce the results toward HD 154368, but not in the cold, dense gas toward AFGL 2136 IRS 1 and Elias 29. Most likely, some combination of simplifying assumptions made in the chemical models are responsible for these discrepancies.

We present Herschel/HIFI observations of 30 transitions of water isotopologues toward the high-mass star forming region NGC 6334 I. The line profiles of H_2^{16}O, H_2^{17}O, H_2^{18}O, and HDO show a complex pattern of emission and absorption components associated with the embedded hot cores, a lower-density envelope, two outflow components, and several foreground clouds, some associated with the NGC 6334 complex, others seen in projection against the strong continuum background of the source. Our analysis reveals an H2O ortho/para ratio of 3 +/- 0.5 in the foreground clouds, as well as the outflow. The water abundance varies from ~10^{-8} in the foreground clouds and the outer envelope to ~10^{-6} in the hot core. The hot core abundance is two orders of magnitude below the chemical model predictions for dense, warm gas, but within the range of values found in other Herschel/HIFI studies of hot cores and hot corinos. This may be related to the relatively low gas and dust temperature (~100 K), or time dependent effects, resulting in a significant fraction of water molecules still locked up in dust grain mantles. The HDO/H_2O ratio in NGC 6334 I, ~2 10^{-4}, is also relatively low, but within the range found in other high-mass star forming regions.

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

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

We present Herschel/HIFI observations of nine transitions of \hho and \hheo towards six high-mass star-forming regions, obtained as part of the PRISMAS Key Program. Water vapor in translucent clouds is detected in absorption along every sightline. We derive the column density of \hho or \hheo for the lower energy level of each transition observed. The total water column density is about a few $10^{13} \rm{cm^{-2}}$. We find that the abundance of water relative to hydrogen nuclei is $1\times10^{-8}$ in agreement with models for oxygen chemistry with high cosmic ray ionization rates. Relative to \hh, the abundance of water is remarkably constant at $5\times10^{-8}$. The abundance of water in excited levels is at most 15%, implying that the excitation temperature $T_{ex}$ in the ground state transitions is below 10 K. The column densities derived from the two ortho ground state transitions indicates that $T_{ex}\simeq5$ K and that the density $n($\hh$)$ in the clouds is $\le10^4 \rm{cm^{-3}}$. For most clouds we derive a water ortho-to-para ratio consistent with the value of 3 expected in thermodynamic equilibrium in the high temperature limit. Two clouds with large column densities exhibit a ratio significantly below 3. This may argue that the history of water molecules includes a cold phase, either when the molecules were formed on cold grains, or when they later become at least partially thermalized with the cold gas ($\sim25$ K) in the shielded, low temperature regions of the clouds; evidently, they have not yet fully thermalized with the warmer ($\sim50$ K) translucent portions of the clouds.

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

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

Absorption lines from the molecules OH+, H2O+, and H3+ have been observed in a diffuse molecular cloud along a line of sight near W51 IRS2. We present the first chemical analysis that combines the information provided by all three of these species. Together, OH+ and H2O+ are used to determine the molecular hydrogen fraction in the outskirts of the observed cloud, as well as the cosmic-ray ionization rate of atomic hydrogen. H3+ is used to infer the cosmic-ray ionization rate of H2 in the molecular interior of the cloud, which we find to be zeta_2=(4.8+-3.4)x10^-16 per second. Combining the results from all three species we find an efficiency factor---defined as the ratio of the formation rate of OH+ to the cosmic-ray ionization rate of H---of epsilon=0.07+-0.04, much lower than predicted by chemical models. This is an important step in the future use of OH+ and H2O+ on their own as tracers of the cosmic-ray ionization rate.

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

Hydrogen fluoride has been established to be an excellent tracer of molecular hydrogen in diffuse clouds. In denser environments, however, the HF abundance has been shown to be approximately two orders of magnitude lower. We present Herschel/HIFI observations of HF J=1-0 toward two high-mass star formation sites, NGC6334 I and AFGL 2591. In NGC6334 I the HF line is seen in absorption in foreground clouds and the source itself, while in AFGL 2591 HF is partially in emission. We find an HF abundance with respect to H2 of 1.5e-8 in the diffuse foreground clouds, whereas in the denser parts of NGC6334 I, we derive a lower limit on the HF abundance of 5e-10. Lower HF abundances in dense clouds are most likely caused by freeze out of HF molecules onto dust grains in high-density gas. In AFGL 2591, the view of the hot core is obstructed by absorption in the massive outflow, in which HF is also very abundant 3.6e-8) due to the desorption by sputtering. These observations provide further evidence that the chemistry of interstellar fluorine is controlled by freeze out onto gas grains.

We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<~1 and one at A_V~3-8, the exact values depending on the ratio of incident ultraviolet flux to gas density. For relatively low values of the incident far ultraviolet flux on the cloud ({\chi}<~ 1000; {\chi}= 1= local interstellar value), the columns of OH+ and H2O+ scale roughly as the cosmic ray primary ionization rate {\zeta}(crp) divided by the hydrogen nucleus density n. The H3O+ column is dominated by the second peak, and we show that if PAHs are present, N(H3O+) ~ 4x10^{13} cm^{-2} independent of {\zeta}(crp) or n. If there are no PAHs or very small grains at the second peak, N(H3O+) can attain such columns only if low ionization potential metals are heavily depleted. We also model diffuse and translucent clouds in the interstellar medium, and show how observations of N(OH+)/N(H) and N(OH+)/N(H2O+) can be used to estimate {\zeta}(crp)/n, {\chi}/n and A_V in them. We compare our models to Herschel observations of these two ions, and estimate {\zeta}(crp) ~ 4-6 x 10^-16 (n/100 cm^-3) s^-1 and \chi/n = 0.03 cm^3 for diffuse foreground clouds towards W49N.

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

The chemical composition of the interstellar medium is determined by gas phase chemistry, assisted by grain surface reactions, and by shock chemistry. The aim of this study is to measure the abundance of the hydroxyl radical (OH) in diffuse spiral arm clouds as a contribution to our understanding of the underlying network of chemical reactions. Owing to their high critical density, the ground states of light hydrides provide a tool to directly estimate column densities by means of absorption spectroscopy against bright background sources. We observed onboard the SOFIA observatory the 2Pi3/2, J = 5/2 3/2 2.5 THz line of ground-state OH in the diffuse clouds of the Carina-Sagittarius spiral arm. OH column densities in the spiral arm clouds along the sightlines to W49N, W51 and G34.26+0.15 were found to be of the order of 10^14 cm^-2, which corresponds to a fractional abundance of 10^-7 to 10^-8, which is comparable to that of H_2O. The absorption spectra of both species have similar velocity components, and the ratio of the derived H_2O to OH column densities ranges from 0.3 to 1.0. In W49N we also detected the corresponding line of ^18OH.

We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the local standard of rest, we infer a total SH column density ~ 2.6 E+12 cm-2, corresponding to an abundance of ~ 7 E-9 relative to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic neutral-neutral reaction SH + H2 -> H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and SH+ in diffuse molecular clouds.

Motivated by recent observations with Herschel/PACS, and the availability of new rate coefficients for the collisional excitation of CO (Yang et al. 2010), the excitation of warm astrophysical CO is revisited with the use of numerical and analytic methods. For the case of an isothermal medium, results have been obtained for a wide range of gas temperatures (100 to 5000 K) and H2 densities (1E+3 to 1E+9 cm-3), and presented in the form of rotational diagrams, in which the logarithm of the column density per magnetic substate, log (N[J]/g[J]), is plotted for each state, as a function of its energy, E[J]. For rotational transitions in the wavelength range accessible to Herschel/PACS, such diagrams are nearly linear when n(H2) > 1E+8 cm-3. When log10(n[H2]) = 6.8 to 8, they exhibit significant negative curvature, whereas when log10(n[H2]) < 4.8 the curvature is uniformly positive throughout the PACS-accessible range. Thus, the observation of a positively-curved CO rotational diagram does not NECESSARILY require the presence of multiple temperature components. Indeed, for some sources observed with Herschel/PACS, the CO rotational diagrams show a modest positive curvature that can be explained by a single isothermal component. Typically, the required physical parameters are H2 densities in the 1E+4 to 1E+5 cm-3 range and temperatures, T, close to the maximum at which CO can survive. Other sources exhibit rotational diagrams with more curvature than can be accounted for by a single temperature component. For the case of a medium with a power-law distribution of gas temperatures, with dN/dT proportional to T to the power -b, results have been obtained for H2 densities 1E+3 to 1E+9 cm-3 and power-law indices, b, in the range 1 to 5; such a medium can account for a CO rotational diagram that is more positively curved than any resulting from an isothermal medium.

Aims. The HIFI instrument onboard Herschel has allowed high spectral resolution and sensitive observations of ground-state transi- tions of three molecular ions: the methylidyne cation CH+, its isotopologue 13CH+, and sulfanylium SH+. Because of their unique chemical properties, a comparative analysis of these cations provides essential clues to the link between the chemistry and dynamics of the diffuse interstellar medium. Methods. The CH+, 13CH+, and SH+ lines are observed in absorption towards the distant high-mass star-forming regions (SFRs) DR21(OH), G34.3+0.1, W31C, W33A, W49N, and W51, and towards two sources close to the Galactic centre, SgrB2(N) and SgrA*+50. All sight lines sample the diffuse interstellar matter along pathlengths of several kiloparsecs across the Galactic Plane. In order to compare the velocity structure of each species, the observed line profiles were deconvolved from the hyperfine structure of the SH+ transition and the CH+, 13CH+, and SH+ spectra were independently decomposed into Gaussian velocity components. To analyse the chemical composition of the foreground gas, all spectra were divided, in a second step, into velocity intervals over which the CH+, 13CH+, and SH+ column densities and abundances were derived. Results. SH+ is detected along all observed lines of sight, with a velocity structure close to that of CH+ and 13CH+. The linewidth distributions of the CH+, SH+, and 13CH+ Gaussian components are found to be similar. These distributions have the same mean (<\delta\u{psion}> ~ 4.2 km s-1) and standard deviation (\sigma(\delta\u{psion}) ~ 1.5 km s-1). This mean value is also close to that of the linewidth distribution of the CH+ visible transitions detected in the solar neighbourhood. We show that the lack of absorption components narrower than 2 km s-1 is not an artefact caused by noise: the CH+, 13CH+, and SH+ line profiles are therefore statistically broader than those of most species detected in absorption in diffuse interstellar gas (e. g. HCO+, CH, or CN). The SH+/CH+ column density ratio observed in the components located away from the Galactic centre spans two orders of magnitude and correlates with the CH+ abundance. Conversely, the ratio observed in the components close to the Galactic centre varies over less than one order of magnitude with no apparent correlation with the CH+ abundance. The observed dynamical and chemical properties of SH+ and CH+ are proposed to trace the ubiquitous process of turbulent dissipation, in shocks or shears, in the diffuse ISM and the specific environment of the Galactic centre regions.

We report the first detection of hydrogen fluoride (HF) toward a high redshift quasar. Using the Caltech Submillimeter Observatory (CSO) we detect the HF J = 1 - 0 transition in absorption toward the Cloverleaf, a broad absorption line (BAL) quasi-stellar object (QSO) at z=2.56. The detection is statistically significant at the ~ 6 sigma level. We estimate a lower limit of 4 \times 1014 cm-2 for the HF column density and using a previous estimate of the hydrogen column density, we obtain a lower limit of 1.7 \times 10-9 for the HF abundance. This value suggests that, assuming a Galactic N(HF)/NH ratio, HF accounts for at least ~10% of the fluorine in the gas phase along the line of sight to the Cloverleaf quasar. This observation corroborates the prediction that HF should be a good probe of the molecular gas at high redshift. Measurements of the HF abundance as a function of redshift are urgently needed to better constrain the fluorine nucleosynthesis mechanism(s).

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

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

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.

We report the detection of emission in the J=1-0 rotational transition of hydrogen fluoride (HF), together with observations of the J=1-0 to J=3-2 rotational lines of H35Cl and H37Cl, towards the envelope of the carbon star IRC +10216. High-sensitivity, high-spectral resolution observations have been carried out with the HIFI instrument on board Herschel, allowing us to resolve the line profiles and providing insights into the spatial distribution of the emission. Our interpretation of the observations, with the use of radiative transfer calculations, indicates that both HF and HCl are formed in the inner regions of the envelope close to the AGB star. Thermochemical equilibrium calculations predict HF and HCl to be the major reservoirs of fluorine and chlorine in the atmospheres of AGB stars. The abundances relative to H2 derived for HF and HCl, 8e-9 and 1e-7 respectively, are substantially lower than those predicted by thermochemical equilibrium, indicating that F and Cl are likely affected by significant depletion onto dust grains, although some chlorine may be in the form of atomic Cl. The H35Cl/H37Cl abundance ratio is 3.3 +/- 0.3. The low abundance derived for HF in IRC +10216 makes it likely that the fluorine abundance is not enhanced over the solar value by nucleosynthesis in the AGB star, although this conclusion may not be robust because the HF abundance we derive is a lower limit to the elemental abundance of F. These observations suggest that both HF and HCl should be detectable through low J rotational transitions in other evolved stars.

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 report a detection of the excited 220-211 rotational transition of para-H2O in APM 08279+5255 using the IRAM Plateau de Bure interferometer. At z = 3.91, this is the highest-redshift detection of interstellar water to date. From LVG modeling, we conclude that this transition is predominantly radiatively pumped and on its own does not provide a good estimate of the water abundance. However, additional water transitions are predicted to be detectable in this source, which would lead to an improved excitation model. We also present a sensitive upper limit for the HF J = 1 - 0 absorption toward APM 08279+5255. While the face-on geometry of this source is not favorable for absorption studies, the lack of HF absorption is still puzzling and may be indicative of a lower fluorine abundance at z = 3.91 compared with the Galactic ISM.

Spitzer-IRS maps of H_2 pure rotational lines from S(0) to S(7) in 3 outflows from Class 0 sources - L1448, BHR71, and NGC2071 - are presented. These lines are used, in conjunction with available rovibrational, near-IR H_2 lines, to probe the physical conditions of the warm gas between hundreds and thousands of Kelvin. We have constructed maps of the molecular hydrogen column density, ortho-to-para ratio and volume density, together with the index beta of the power law describing the distribution of gas temperature. In all three outflows, the present ortho-to-para ratio significantly deviates from the high temperature equilibrium of 3, being on average between 2.0 and 2.3. These low values, that reflect the young age of these flows, are found also in regions of relatively high temperature (~ 1000 K), likely indicating that shocks are occurring in a time shorter than that needed for a complete para to ortho conversion. Density maps indicate upper limits close to LTE conditions, i.e. between 10^6-10^7 cm^-3; moreover we demonstrate, based on the detections of HD emission spots (R(3) and R(4) lines), that a density stratification does exist, with the low density components (10^4-10^5 cm^-3) associated with the coldest gas. The beta index is found in all flows to be above 3.8: this value is consistent with predictions of multiple C-type bow shocks with a range of velocities, some of which insufficient to achieve the temperature at which H_2 partially dissociates. The contribution of H_2 to the total cooling is quantitatively similar to that of other abundant molecules emitting in the far-infrared, such as water and CO; moreover, the luminosity radiated away is comparable with the estimated kinetic energy of the swept-out outflow: this supports the scenario in which the shock from which H_2 is emitted is also capable of accelerating the molecular outflow driven at the shock working surface.

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.

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

We report the results of observations of ten rotational transitions of water vapor toward the carbon-rich AGB (asymptotic giant branch) star IRC+10216 (CW Leonis), carried out with Herschel's HIFI instrument. Each transition was securely detected by means of observations using the dual beam switch mode of HIFI. The measured line ratios imply that water vapor is present in the inner outflow at small distances (few x 1.E+14 cm) from the star, confirming recent results reported by Decin et al. from observations with Herschel's PACS and SPIRE instruments. This finding definitively rules out the hypothesis that the observed water results from the vaporization of small icy objects in circular orbits. The origin of water within the dense C-rich envelope of IRC+10216 remains poorly understood. We derive upper limits on the H2-17O/H2-16O and H2-18O/H2-16O isotopic abundance ratios of ~ 5.E-3 (3 sigma), providing additional constraints on models for the origin of the water vapor in IRC+10216.

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 have studied the shock-excited molecular regions associated with four supernova remnants (SNRs) - IC443C, W28, W44 and 3C391 - and two Herbig-Haro objects, HH7 and HH54, using Spitzer's Infrared Spectrograph (IRS). The physical conditions within the observed areas are inferred from spectroscopic data obtained from IRS and from SWS and LWS onboard ISO, together with photometric data from Spitzer's Infrared Array Camera (IRAC). Adopting a power-law distribution for the gas temperature in the observed region, the H2 S(0) to S(7) spectral line maps obtained with IRS were used to constrain the gas density, yielding estimated n(H2) in the range 2-4*10^3 cm^-3. The excitation of H2 S(9) to S(12) and high-J CO pure rotational lines, however, require environments several times denser. The inconsistency among the best-fit densities estimated from different species can be explained by density fluctuations within the observed regions. The best-fit power-law index b is smaller than the value 3.8 predicted for a paraboloidal C-type bow shock, suggesting that the shock front has a "flatter" shape than that of a paraboloid. The best-fit parameters for SNRs and Herbig-Haro objects do not differ significantly between the two classes of sources, except that for the SNRs the ortho-to-para ratio (OPR) of hot gas (T> 1000 K) is close to the LTE value 3, while for HH7 and HH54 even the hottest gas exhibits an OPR smaller than 3. Finally, we mapped the physical parameters within the regions observed with IRS and found that the mid-lying H2 emissions - S(3) to S(5) - tend to trace the hot component of the gas, while the intensities of S(6) and S(7) are more sensitive to the density of the gas compared to S(3) to S(5).

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.

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

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

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.

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.

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

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.

Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predictions from conventional chemical models by many orders of magnitude. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen-rich stars using the NH3 ortho J_K = 1_0 - 0_0 ground-state line. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 J_K = 1_0 - 0_0 transition near 572.5 GHz, simultaneously with the ortho-H2O J_Ka,Kc = 1_1,0 -1_0,1 transition, toward VY CMa, OH 26.5+0.6, IRC+10420, and IK Tau. We conducted non-LTE radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the latter two stars, Very Large Array imaging of NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those obtained for the ground state H2O line. The NH3 abundances relative to H2 are very high and range from 2 x 10-7 to 3 x 10-6 for the objects we have studied. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield satisfactory fits to the observed line profiles, it leads to abundance estimates that confirm the very high values found in earlier studies. New ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array.

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.

The HIFI instrument on board the Herschel Space Observatory has been used to observe interstellar nitrogen hydrides along the sight-line towards G10.6-0.4 in order to improve our understanding of the interstellar chemistry of nitrogen. We report observations of absorption in NH N=1-0, J=2-1 and ortho-NH2 1_1,1-0_0,0. We also observed ortho-NH3 1_0-0_0, and 2_0-1_0, para-NH3 2_1-1_1, and searched unsuccessfully for NH+. All detections show emission and absorption associated directly with the hot-core source itself as well as absorption by foreground material over a wide range of velocities. All spectra show similar, non-saturated, absorption features, which we attribute to diffuse molecular gas. Total column densities over the velocity range 11-54 km/s are estimated. The similar profiles suggest fairly uniform abundances relative to hydrogen, approximately 6*10^-9, 3*10^-9, and 3*10^-9 for NH, NH2, and NH3, respectively. These abundances are discussed with reference to models of gas-phase and surface chemistry.

We discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

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 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 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 performed Herschel/HIFI observations of several CO lines in the far-infrared/sub-mm in the protoplanetary nebula CRL618. The high spectral resolution provided by HIFI allows measurement of the line profiles. Since the dynamics and structure of the nebula is well known from mm-wave interferometric maps, it is possible to identify the contributions of the different nebular components (fast bipolar outflows, double shells, compact slow shell) to the line profiles. The observation of these relatively high-energy transitions allows an accurate study of the excitation conditions in these components, particularly in the warm ones, which cannot be properly studied from the low-energy lines. The 12CO J=16-15, 10-9, and 6-5 lines are easily detected in this source. 13CO J=10-9 and 6-5 are also detected. Wide profiles showing spectacular line wings have been found, particularly in 12CO 16-15. Other lines observed simultaneously with CO are also shown. Our analysis of the CO high-J transitions, when compared with the existing models, confirms the very low expansion velocity of the central, dense component, which probably indicates that the shells ejected during the last AGB phases were driven by radiation pressure under a regime of maximum transfer of momentum. No contribution of the diffuse halo found from mm-wave data is identified in our spectra, because of its low temperature. We find that the fast bipolar outflow is quite hot, much hotter than previously estimated; for instance, gas flowing at 100 km/s must have a temperature higher than ~ 200 K. Probably, this very fast outflow, with a kinematic age < 100 yr, has been accelerated by a shock and has not yet cooled down. The double empty shell found from mm-wave mapping must also be relatively hot, in agreement with the previous estimate.

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

A set of new, sensitive, and spectrally resolved, sub-millimeter line observations are used to probe the warm circumstellar gas around the S-type AGB star chi Cyg. The observed lines involve high rotational quantum numbers, which, combined with previously obtained lower-frequency data, make it possible to study in detail the chemical and physical properties of, essentially, the entire circumstellar envelope of chi Cyg. The data were obtained using the HIFI instrument aboard Herschel, whose high spectral resolution provides valuable information about the line profiles. Detailed, non-LTE, radiative transfer modelling, including dust radiative transfer coupled with a dynamical model, has been performed to derive the temperature, density, and velocity structure of the circumstellar envelope. We report the first detection of circumstellar H2O rotational emission lines in an S-star. Using the high-J CO lines to derive the parameters for the circumstellar envelope, we modelled both the ortho- and para-H2O lines. Our modelling results are consistent with the velocity structure expected for a dust-driven wind. The derived total H2O abundance (relative to H2) is (1.1 +/- 0.2)E-5, much lower than that in O-rich stars. The derived ortho-to-para ratio of 2.1 +/- 0.6 is close to the high-temperature equilibrium limit, consistent with H2O being formed in the photosphere.

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

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

We report the first detection of the ground-state rotational transition of the methylidyne cation CH+ towards the massive star-forming region DR21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR21 molecular ridge and foreground gas. These observations allow us to determine a CH+ J=1-0 line frequency of 835137 +/- 3 MHz, in good agreement with a recent experimental determination. We estimate the CH+ column density to be a few 1e13 cm^-2 in the gas seen in emission, and > 1e14 cm^-2 in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH+]/[H] > 1.2x10^-8. We show that the CH+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

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

We report the discovery of water vapour toward the carbon star V Cygni. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1(11) - 0(00) para-water transition at 1113.3430 GHz in the upper sideband of the Band 4b receiver. The observed spectral line profile is nearly parabolic, but with a slight asymmetry associated with blueshifted absorption, and the integrated antenna temperature is 1.69 \pm 0.17 K km/s. This detection of thermal water vapour emission, carried out as part of a small survey of water in carbon-rich stars, is only the second such detection toward a carbon-rich AGB star, the first having been obtained by the Submillimeter Wave Astronomy Satellite toward IRC+10216. For an assumed ortho-to-para ratio of 3 for water, the observed line intensity implies a water outflow rate ~ (3 - 6) E-5 Earth masses per year and a water abundance relative to H2 of ~ (2-5) E-6. This value is a factor of at least 1E+4 larger than the expected photospheric abundance in a carbon-rich environment, and - as in IRC+10216 - raises the intriguing possibility that the observed water is produced by the vapourisation of orbiting comets or dwarf planets. However, observations of the single line observed to date do not permit us to place strong constraints upon the spatial distribution or origin of the observed water, but future observations of additional transitions will allow us to determine the inner radius of the H2O-emitting zone, and the H2O ortho-to-para ratio, and thereby to place important constraints upon the origin of the observed water emission.

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

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

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

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 report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

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.

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.

An artificial neural network (ANN) is investigated as a tool for estimating rate coefficients for the collisional excitation of molecules. The performance of such a tool can be evaluated by testing it on a dataset of collisionally-induced transitions for which rate coefficients are already known: the network is trained on a subset of that dataset and tested on the remainder. Results obtained by this method are typically accurate to within a factor ~ 2.1 (median value) for transitions with low excitation rates and ~ 1.7 for those with medium or high excitation rates, although 4% of the ANN outputs are discrepant by a factor of 10 more. The results suggest that ANNs will be valuable in extrapolating a dataset of collisional rate coefficients to include high-lying transitions that have not yet been calculated. For the asymmetric top molecules considered in this paper, the favored architecture is a cascade-correlation network that creates 16 hidden neurons during the course of training, with 3 input neurons to characterize the nature of the transition and one output neuron to provide the logarithm of the rate coefficient.

Although they are only minor constituents of the interstellar medium, halogen-containing molecules are of special interest because of their unique thermochemistry. Here, we present a theoretical study of the chemistry of interstellar molecules containing the halogen elements chlorine and fluorine. We have modeled both diffuse and dense molecular clouds, making use of updated estimates for the rates of several key chemical processes. We present predictions for the abundances of the three halogen molecules that have been detected to date in the interstellar medium: HF, CF+ and HCl. As in our previous study of fluorine-bearing interstellar molecules, we predict HF to be the dominant gas-phase reservoir of fluorine within both diffuse and dense molecular clouds; we expect the Herschel Space Observatory to detect widespread absorption in the HF J=1-0 transition. Our updated model now overpredicts the CF+ abundance by a factor ~10 relative to observations of the Orion Bar; this discrepancy has widened because we now adopt a laboratory measurement of the CF+ dissociative recombination rate that is smaller than the estimate we adopted previously. This disagreement suggests that the reaction of C+ with HF proceeds more slowly than the capture rate assumed in our model; a laboratory measurement of this reaction rate would be very desirable. Our model predicts diffuse cloud HCl abundances that are similar to those predicted previously and detected tentatively toward zeta Oph. Two additional species are potentially detectable from photodissociation regions: the H2Cl+ and HCl+ molecular ions. Ortho-H2Cl+ has its lowest-lying transition in the millimeter spectral region observable from the ground, and the lowest rotational transition of HCl+ is observable with Herschel's HIFI instrument.

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

This white paper is submitted to the Astronomy and Astrophysics 2010 Decadal Survey (Astro2010)1 Committee on the State of the Profession to emphasize the potential of the Stratospheric Observatory for Infrared Astronomy (SOFIA) to contribute to the training of instrumentalists and observers, and to related technology developments. This potential goes beyond the primary mission of SOFIA, which is to carry out unique, high priority astronomical research. SOFIA is a Boeing 747SP aircraft with a 2.5 meter telescope. It will enable astronomical observations anywhere, any time, and at most wavelengths between 0.3 microns and 1.6 mm not accessible from ground-based observatories. These attributes, accruing from the mobility and flight altitude of SOFIA, guarantee a wealth of scientific return. Its instrument teams (nine in the first generation) and guest investigators will do suborbital astronomy in a shirt-sleeve environment. The project will invest $10M per year in science instrument development over a lifetime of 20 years. This, frequent flight opportunities, and operation that enables rapid changes of science instruments and hands-on in-flight access to the instruments, assure a unique and extensive potential - both for training young instrumentalists and for encouraging and deploying nascent technologies. Novel instruments covering optical, infrared, and submillimeter bands can be developed for and tested on SOFIA by their developers (including apprentices) for their own observations and for those of guest observers, to validate technologies and maximize observational effectiveness.

We explore the design of a space mission, Project Lyman, which has the goal of quantifying the ionization history of the universe from the present epoch to a redshift of z ~ 3. Observations from WMAP and SDSS show that before a redshift of z >~ 6 the first collapsed objects, possibly dwarf galaxies, emitted Lyman continuum (LyC) radiation shortward of 912 A, reionizing most of the universe. How LyC escapes from galactic environments, whether it induces positive or negative feedback on the local and global collapse of structures, and the role played by clumping, molecules, metallicity and dust are major unanswered theoretical questions, requiring observational constraint. Numerous intervening Lyman limit systems, which frustrate the detection of LyC from high z objects, thin below z ~ 3 where there are a few objects with apparently very high fesc. At low z there are only controversial detections and a handful of upper limits. A wide-field multi-object spectroscopic survey with moderate spectral and spatial resolution can quantify fesc within diverse spatially resolved galactic environments over redshifts with significant evolution in galaxy assemblage and quasar activity. It can also calibrate LyC escape against Ly-alpha escape, providing an essential tool to JWST for probing the beginnings of reionization. We present calculations showing the evolution of the characteristic apparent magnitude of star-forming galaxy luminosity functions at 900 A, as a function of redshift and assumed escape fraction to determine the required aperture for detecting LyC. We review our efforts to build a pathfinding dual order multi-object spectro/telescope with a (0.5deg)^2 field-of-view, using a GSFC microshutter array, and crossed delay-line micro-channel plate detector.

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.

Photometric maps, obtained with Spitzer's Infrared Array Camera (IRAC), can provide a valuable probe of warm molecular hydrogen within the interstellar medium. IRAC maps of the supernova remnant IC443, extracted from the Spitzer archive, are strikingly similar to spectral line maps of the H2 pure rotational transitions that we obtained with the Infrared Spectrograph (IRS) instrument on Spitzer. IRS spectroscopy indicates that IRAC Bands 3 and 4 are indeed dominated by the H2 v=0-0 S(5) and S(7) transitions, respectively. Modeling of the H2 excitation suggests that Bands 1 and 2 are dominated by H2 v=1-0 O(5) and v=0-0 S(9). Large maps of the H2 emission in IC433, obtained with IRAC, show band ratios that are inconsistent with the presence of gas at a single temperature. The relative strengths of IRAC Bands 2, 3, and 4 are consistent with pure H2 emission from shocked material with a power-law distribution of gas temperatures. CO vibrational emissions do not contribute significantly to the observed Band 2 intensity. Assuming that the column density of H2 at temperatures T to T+dT is proportional to T raised to the power -b for temperatures up to 4000 K, we obtained a typical estimate of 4.5 for b. The power-law index, b, shows variations over the range 3 to 6 within the set of different sight-lines probed by the maps, with the majority of sight-lines showing b in the range 4 to 5. The observed power-law index is consistent with the predictions of simple models for paraboloidal bow shocks.

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.

The first map of interstellar acetylene (C2H2) has been obtained with the infrared spectrograph onboard the Spitzer Space Telescope. A spectral line map of the $\nu_5$ vibration-rotation band at 13.7 microns carried out toward the star-forming region Cepheus A East, shows that the C2H2 emission peaks in a few localized clumps where gas-phase CO2 emission was previously detected with Spitzer. The distribution of excitation temperatures derived from fits to the C2H2 line profiles ranges from 50 to 200 K, a range consistent with that derived for gaseous CO2 suggesting that both molecules probe the same warm gas component. The C2H2 molecules are excited via radiative pumping by 13.7 microns continuum photons emanating from the HW2 protostellar region. We derive column densities ranging from a few x 10^13 to ~ 7 x 10^14 cm^-2, corresponding to C2H2 abundances of 1 x 10^-9 to 4 x 10^-8 with respect to H2. The spatial distribution of the C2H2 emission along with a roughly constant N(C2H2)/N(CO2) strongly suggest an association with shock activity, most likely the result of the sputtering of acetylene in icy grain mantles.

We have modeled the emission of H2O rotational lines from the extreme C-rich star IRC+10216. Our treatment of the excitation of H2O emissions takes into account the excitation of H2O both through collisions, and through the pumping of the nu2 and nu3 vibrational states by dust emission and subsequent decay to the ground state. Regardless of the spatial distribution of the water molecules, the H2O 1_{10}-1_{01} line at 557 GHz observed by the Submillimeter Wave Astronomy Satellite (SWAS) is found to be pumped primarily through the absorption of dust-emitted photons at 6 $\mu$m in the nu2 band. As noted by previous authors, the inclusion of radiative pumping lowers the ortho-H2O abundance required to account for the 557 GHz emission, which is found to be (0.5-1)x10^{-7} if the presence of H2O is a consequence of vaporization of orbiting comets or Fischer-Tropsch catalysis. Predictions for other submillimeter H2O lines that can be observed by the Herschel Space Observatory (HSO) are reported. Multitransition HSO observations promise to reveal the spatial distribution of the circumstellar water vapor, discriminating among the several hypotheses that have been proposed for the origin of the H2O vapor in the envelope of IRC+10216. We also show that, for observations with HSO, the H2O 1_{10}-1_{01} 557 GHz line affords the greatest sensitivity in searching for H2O in other C-rich AGB stars.

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.

By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the planets in our solar system. The technical challenges associated with interferometry in the far-IR are greatly relaxed relative to those encountered at shorter wavelengths or when starlight nulling is required. A structurally connected far-IR interferometer with a maximum baseline length of 36 m can resolve the interesting spatial structures in nascent and developed exoplanetary systems and measure exozodiacal emission at a sensitivity level critical to TPF-I mission planning. The Space Infrared Interferometric Telescope was recommended in the Community Plan for Far-IR/Submillimeter Space Astronomy, studied as a Probe-class mission, and estimated to cost 800M dollars. The scientific communities in Europe, Japan, and Canada have also demonstrated a keen interest in far-IR interferometry through mission planning workshops and technology research, suggesting the possibility of an international collaborative effort.

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.

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

On 4 July 2005 at 5:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the 1(10)-1(01) ortho-water ground-state rotational transition in comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS and we derive an upper limit of 1.8e7 kg on the water ice evaporated by the impact. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before and after the impact. Episodes of increased activity with Q(H2O)~1e28 molecule/s alternated with periods with low outgassing (Q(H2O)<~5e27 molecule/s). We estimate that 9P/Tempel 1 vaporized a total of N~4.5e34 water molecules (~1.3e9 kg) during June-September 2005. Our observations indicate that only a small fraction of the nucleus of Tempel 1 appears to be covered with active areas. Water vapor is expected to emanate predominantly from topographic features periodically facing the Sun as the comet rotates. We calculate that appreciable asymmetries of these features could lead to a spin-down or spin-up of the nucleus at observable rates.

We discuss the first astronomical detection of the CF+ (fluoromethylidynium) ion, obtained by observations of the J=1-0 (102.6 GHz), J=2-1 (205.2 GHz) and J=3-2 (307.7 GHz) rotational transitions toward the Orion Bar region. Our search for CF+, carried out using the IRAM 30m and APEX 12m telescopes, was motivated by recent theoretical models that predict CF+ abundances of a few times 1.E-10 in UV-irradiated molecular regions where C+ is present. The CF+ ion is produced by exothermic reactions of C+ with HF. Because fluorine atoms can react exothermically with H2, HF is predicted to be the dominant reservoir of fluorine, not only in well-shielded regions but also in the surface layers of molecular clouds where the C+ abundance is large. The observed CF+ line intensities imply the presence of CF+ column densities of at least 1.E+12 cm-2 over a region of size at least ~ 1 arcmin, in good agreement with theoretical predictions. They provide support for our current theories of interstellar fluorine chemistry, which suggest that hydrogen fluoride should be ubiquitous in interstellar gas clouds and widely detectable in absorption by future satellite and airborne observatories.

On 4 July 2005 at 1:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. The 370 kg impactor collided with the sunlit side of the nucleus with a relative velocity of 10.2 km/s. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the 1(10)-1(01) ortho-water ground-state rotational transition in comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before the impact.

We report the first astronomical detection of the CF+ (fluoromethylidynium) ion obtained by recent observations of its J = 1 - 0 (102.6 GHz), J = 2 - 1 (205.2 GHz), and J = 3 - 2 (307.7 GHz) pure rotational emissions toward the Orion Bar. Our search for CF+, carried out using the IRAM 30m and APEX 12m telescopes, was motivated by recent theoretical models that predict CF+ abundances of a few x E-10 in UV-irradiated molecular regions where C+ is present. The measurements confirm the predictions. They provide support for our current theories of interstellar fluorine chemistry, which suggest that hydrogen fluoride should be ubiquitous in interstellar gas clouds.

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 a theoretical investigation of the chemistry of fluorine-bearing molecules in diffuse and dense interstellar gas clouds. The chemistry of interstellar fluorine is qualitatively different from that of any other element, because - unlike the neutral atoms of any other element found in diffuse or dense molecular clouds - atomic fluorine undergoes an exothermic reaction with molecular hydrogen. Over a wide range of conditions attained within interstellar gas clouds, the product of that reaction - hydrogen fluoride - is predicted to be the dominant gas-phase reservoir of interstellar fluorine nuclei. Our model predicts HF column densities ~ 1.E+13 cm-2 in dark clouds and column densities as large as 1.E-11 cm-2 in diffuse interstellar gas clouds with total visual extinctions as small as 0.1 mag. Such diffuse clouds will be detectable by means of absorption line spectroscopy of the J = 1 - 0 transition at 243.2 micron using the Stratospheric Observatory for Infrared Astronomy (SOFIA) and the Herschel Space Observatory (HSO). The CF+ ion is predicted to be the second most abundant fluorine-bearing molecule, with typical column densities a factor ~ 100 below those of HF; with its lowest two rotational transitions in the millimeter-wave spectral region, CF+ may be detectable from ground-based observatories. HF absorption in quasar spectra is a potential probe of molecular gas at high redshift, providing a possible bridge between the UV/optical observations capable of probing H2 in low column density systems and the radio/millimeter-wavelength observations that probe intervening molecular clouds of high extinction and large molecular fraction; at redshifts beyond ~ 0.3, HF is potentially detectable from ground-based submillimeter observatories in several atmospheric transmission windows.

We report the detection of H2CO (formaldehyde) around the carbon-rich AGB star, IRC+10216. We find a fractional abundance with respect to molecular hydrogen of x(H2CO)= (1.3 {+1.5}{-0.8}) x 10^{-8}. This corresponds to a formaldehyde abundance with respect to water vapor of x(H2CO)/x(H2O)=(1.1 +/- 0.2) x 10^{-2}, in line with the formaldehyde abundances found in Solar System comets, and indicates that the putative extrasolar cometary system around IRC+10216 may have a similar chemical composition to Solar System comets. However, we also failed to detect CH3OH (methanol) around IRC+10216 and our upper limit of x(CH3OH)/x(H2O) < 7.7 x 10^{-4}, (3 sigma), indicates that methanol is substantially underabundant in IRC+10216, compared to Solar System comets. We also conclude, based on offset observations, that formaldehyde has an extended source in the envelope of IRC+10216 and may be produced by the photodissociation of a parent molecule, similar to the production mechanism for formaldehyde in Solar System comet comae. Preliminary mapping observations also indicate a possible asymmetry in the spatial distribution of formaldehyde around IRC+10216, but higher signal-to-noise observations are required to confirm this finding. This study is based on observations carried out with the IRAM 30m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). (abridged)

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.

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.

We report the detection of the 1665 and 1667 MHz main lines of OH (hydroxyl) and upper limits on the 1612 MHz satellite line of OH towards the carbon-rich AGB star, IRC+10216. We find a beam averaged fractional abundance x(OH)~4x10^{-8}. This detection supports the identification by Melnick et al. (2001) of the 1_10-1_01 transition of water vapor with a 556.936 GHz emission feature detected towards IRC+10216, since OH is the expected photodissociation product of water vapor. The shape of the OH lines, however, differs significantly from the shape expected based on the observations of Melnick et al. (2001). Possible explanations for the anomalous shapes of the 1665 and 1667 MHz lines are discussed. The most likely explanations for the unexpected OH line shapes are either masing or an asymmetric distribution of OH molecules around IRC+10216.

ESPI has been proposed for direct imaging and spectral analysis of giant planets orbiting solar-type stars. ESPI extends the concept suggested by Nisenson and Papaliolios (2001) for a square aperture apodized telescope that has sufficient dynamic range to directly detect exo-planets. With a 1.5 M square mirror, ESPI can deliver high dynamic range imagery as close as 0.3 arcseconds to bright sources, permitting a sensitive search for exoplanets around nearby stars and a study of their characteristics in reflected light.

Observations of W51 with the Submillimeter Wave Astronomy Satellite (SWAS) have yielded the first detection of water vapor in a diffuse molecular cloud. The water vapor lies in a foreground cloud that gives rise to an absorption feature at an LSR velocity of 6 km/s. The inferred H2O column density is 2.5E+13 cm-2. Observations with the Arecibo radio telescope of hydroxyl molecules at ten positions in W51 imply an OH column density of 8E+13 cm-2 in the same diffuse cloud. The observed H2O/OH ratio of ~ 0.3 is significantly larger than an upper limit derived previously from ultraviolet observations of the similar diffuse molecular cloud lying in front of HD 154368. The observed variation in H2O/OH likely points to the presence in one or both of these clouds of a warm (T > 400) gas component in which neutral-neutral reactions are important sources of OH and/or H2O.

We argue that the presence of water vapor in the circumstellar outflow of a carbon-rich AGB star is potentially a distinctive signature of extra-solar cometary systems. Detailed models show that at suitable distances from the star, water ice can survive well into the carbon-rich AGB phase; water vapor abundances as large as 10^-6 could result from the vaporization of a collection of orbiting icy bodies with a total mass comparable to what might have been originally present in the solar system's Kuiper Belt. In particular, the recently-reported detection by the Submillimeter Wave Astronomy Satellite of water vapor in the circumstellar outflow of IRC+10216 can be explained if ~10 Earth masses of ice is present at a distance ~300 AU from that carbon-rich star. Future observations with the Herschel Space Observatory (HSO, formerly known as FIRST) will facilitate sensitive multi-transition observations of water, yielding line ratios that can establish the radial distribution of water vapor in IRC+10216. The greater sensitivity of HSO will also allow searches for water vapor to be carried out in a much larger sample of carbon-rich AGB stars.

Since 1995, astronomers have discovered planets with masses comparable to that of Jupiter (318 times Earth's mass) in orbit around approximately 60 stars. Although unseen directly, the presence of these planets is inferred by the small reflex motions that they gravitationally induce on the star they orbit; these result in small periodic wavelength shifts in the stellar spectrum. Since this method favors the detection of massive objects orbiting in close proximity to the star, the question of whether these systems also contain analogs of the smaller constituents of our Solar System has remained unanswered. Using an alternative approach, we report here observations of an aging carbon-star, IRC+10216, that reveal the presence of circumstellar water vapor, a molecule not expected in measurable abundances around such a star and thus a distinctive signature of an orbiting cometary system. The only plausible explanation for this water vapor is that the recent evolution of IRC+10216 - which is accompanied by a prodigious increase in its luminosity - is now causing the vaporization of a collection of orbiting icy bodies, a process first considered in a previous theoretical study.

We have used the Submillimeter Wave Astronomy Satellite to map the ground-state 1_{10}-1_{01} transition of ortho-water at 557 GHz in the Orion molecular cloud. Water emission was detected in Orion over an angular extent of about 20 arcmin, or nearly 3 pc. The water emission is relatively weak, with line widths (3-6 km s^{-1}) and V_{LSR} velocities (9-11 km s^{-1}) consistent with an origin in the cold gas of the molecular ridge. We find that the ortho-water abundance relative to H_2 in the extended gas in Orion varies between 1 and 8x10^{-8}, with an average of 3x10^{-8}. The absence of detectable narrow-line ortho-H_2^{18}O emission is used to set a 3-sigma upper limit on the relative ortho-water abundance of 7x10^{-8}.

We mapped a 13 by 22 pc region in emission from 492 GHz [C I] and, for the first time, 551 GHz ^{13}CO(5-4) in the giant molecular cloud M17SW, using the Submillimeter Wave Astronomy Satellite. The morphologies of the [C I] and ^{13}CO emission are strikingly similar. The extent and intensity of the [C I] and ^{13}CO(5-4) emission is explained as arising from photodissociation regions on the surfaces of embedded molecular clumps. Modeling of the ^{13}CO(5-4) emission in comparison to ^{13}CO(1-0) indicates a temperature gradient across the cloud, peaking to at least 63 K near the M17 ionization front and decreasing to at least 20 K at the western edge of the cloud. We see no correlation between gas density and column density. The beam-averaged column density of C I in the core is 1x10^{18} cm^-2, and the mean column density ratio N(C I)/N(CO) is about 0.4. The variations of N(C I)/N(CO) with position in M17SW indicate a similar clump size distribution throughout the cloud.

We present Submillimeter Wave Astronomy Satellite (SWAS) observations of the 1_{10}-1_{01} transition of ortho-water at 557 GHz toward 12 molecular cloud cores. The water emission was detected in NGC 7538, Rho Oph A, NGC 2024, CRL 2591, W3, W3(OH), Mon R2, and W33, and was not detected in TMC-1, L134N, and B335. We also present a small map of the water emission in S140. Observations of the H_2^{18}O line were obtained toward S140 and NGC 7538, but no emission was detected. The abundance of ortho-water relative to H_2 in the giant molecular cloud cores was found to vary between 6x10^{-10} and 1x10^{-8}. Five of the cloud cores in our sample have previous water detections; however, in all cases the emission is thought to arise from hot cores with small angular extents. The water abundance estimated for the hot core gas is at least 100 times larger than in the gas probed by SWAS. The most stringent upper limit on the ortho-water abundance in dark clouds is provided in TMC-1, where the 3-sigma upper limit on the ortho-water fractional abundance is 7x10^{-8}.

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.

Circumnuclear water masers can respond in two ways to changes in the luminosity of an active galactic nucleus. First, an increase in the X-ray luminosity can lead to an increase in the maser emissivity; and second, an increase in the intrinsic bolometric luminosity may result in a temporary decrease in the difference between the gas and dust temperature and a consequent decrease in the maser output. Whilst the latter effect can occur over a period shorter than the thermal timescale, the former effect cannot. Quantitative estimates of the response of the water maser emissivity to changes in either the X-ray or bolometric luminosity are presented, together with estimates of the relevant timescales. Either mechanism could account for recent observations by Gallimore et al. which suggest that the water maser variability in two widely separated regions of the circumnuclear gas in NGC 1068 have been coordinated by a signal from the active nucleus. For either mechanism, a minimum H2 density ~ 1.E+8 cm-3 is needed to explain the observed variability timescale.

SAFIRE is a versatile imaging Fabry-Perot spectrograph covering 145 to 655 microns, with spectral resolving powers ranging over 5-10,000. Selected as a "PI" instrument for the airborne Stratospheric Observatory for Infrared Astronomy (SOFIA), SAFIRE will apply two-dimensional pop-up bolometer arrays to provide background-limited imaging spectrometry. Superconducting transition edge bolometers and SQUID multiplexers are being developed for these detectors. SAFIRE is expected to be a "First Light" instrument, useable during the initial SOFIA operations. Although a PI instrument rather than a "Facility Class" science instrument, it will be highly integrated with the standard SOFIA planning, observation, and data analysis tools.

We report the detection of numerous far-infrared emission lines of water vapor toward the supergiant star VY Canis Majoris. A 29.5 - 45 micron grating scan of VY CMa, obtained using the Short Wavelength Spectrometer (SWS) of the Infrared Space Observatory (ISO) at a spectral resolving power of approximately 2000, reveals at least 41 spectral features due to water vapor that together radiate a total luminosity ~ 25 solar luminosities. In addition to pure rotational transitions within the ground vibrational state, these features include rotational transitions within the (010) excited vibrational state. The spectrum also shows the doublet Pi 1/2 (J=5/2) <-- doublet Pi 3/2 (J=3/2) OH feature near 34.6 micron in absorption. Additional SWS observations of VY CMa were carried out in the instrument's Fabry-Perot mode for three water transitions: the 7(25)-6(16) line at 29.8367 micron, the 4(41)-3(12) line 31.7721 micron, and the 4(32)-3(03) line at 40.6909 micron. The higher spectral resolving power of approximately 30,000 thereby obtained permits the line profiles to be resolved spectrally for the first time and reveals the "P Cygni" profiles that are characteristic of emission from an outflowing envelope.

Observations at far-infrared and submillimeter wavelengths promise to revolutionize the study of high redshift galaxies and AGN by providing a unique probe of the conditions within heavily extinguished regions of star formation and nuclear activity. Observational capabilities in this spectral region will expand greatly in the next decade as new observatories are developed both in space and on the ground. These facilities include the Space Infrared Telescope Facility (SIRTF), the far-infrared and submillimeter telescope (FIRST), and the millimeter array (MMA). In the longer term, the requirements of high angular resolution (comparable to that of HST), full wavelength coverage, and high sensitivity (approaching the fundamental limit imposed by photon counting statistics) will motivate the development of far-IR and submillimeter space interferometry using cold telescopes and incoherent detector arrays.

A major goal of modern astrophysics is to understand the processes by which the universe evolved from its initial simplicity, as seen in measurements of the Cosmic Microwave Background, to the universe we see today, with complexity on all scales. While the diffuse background measurements of COBE reveal the importance of the far infrared and submillimeter in early galaxy and star formation, the understanding of the development of complex structure requires high resolution imaging and spectroscopy. We present a concept for a space mission called SPECS, the Submillimeter Probe of the Evolution of Cosmic Structure, which provides these capabilities. SPECS is a cold, spatial and spectral Michelson interferometer with adjustable baselines ranging up to 1 km. It has Hubble sensitivity and angular resolution in the far IR and submillimeter, spectral resolution up to 10,000, and a 14' field of view. SPECS will be able to image thermal dust continuum and infrared cooling and diagnostic line emission over a wide range of redshifts, providing extinction-free astrophysical probes of young galaxies and early cosmic structures, and measures of the luminosity and heavy element formation history of the universe. SPECS would also have the potential to improve vastly our knowledge of protostars, protoplanetary systems, Active Galactic Nuclei and other objects in the local universe. We recommend that a concerted effort be made during the next decade to develop and test certain critical technologies (photon-counting far IR detectors, formation flying spacecraft, cold, lightweight mirrors, and active coolers), so that SPECS can be deployed in about 2015.

We discuss the ratio of ortho- to para-H2 in photodissociation regions (PDRs). We draw attention to an apparent confusion in the literature between the ortho-to-para ratio of molecules in FUV-pumped vibrationally excited states, and the H2 ortho-to-para abundance ratio. These ratios are not the same because the process of FUV-pumping of fluorescent H2 emission in PDRs occurs via optically thick absorption lines. Thus, gas with an equilibrium ratio of ortho- to para-H2 equal to 3 will yield FUV-pumped vibrationally excited ortho-to-para ratios smaller than 3, because the ortho-H2 pumping rates are preferentially reduced by optical depth effects. Indeed, if the ortho and para pumping lines are on the ``square root'' part of the curve-of-growth, then the expected ratio of ortho and para vibrational line strengths is the square root of 3, ~ 1.7, close to the typically observed value. Thus, contrary to what has sometimes been stated in the literature, most previous measurements of the ratio of ortho- to para-H2 in vibrationally excited states are entirely consistent with a total ortho-to-para ratio of 3, the equilibrium value for temperatures greater than 200 K. We present an analysis and several detailed models which illustrate the relationship between the total ratios of ortho- to para-H2 and the vibrationally excited ortho-to-para ratios in PDRs. Recent Infrared Space Observatory (ISO) measurements of pure rotational and vibrational H2 emissions from the PDR in the star-forming region S140 provide strong observational support for our conclusions.

We have used a coupled dynamical and chemical model to examine the chemical changes induced by the passage of an interstellar shock in well shielded regions. Using this model we demonstrate that the formation of water in a shock will be followed in the post--shock phase by depletion of the water molecules onto the grain surfaces. To attempt to discriminate between the creation of ices behind shocks and their production by means of grain surface chemistry, we examine the deuterium chemistry of water before, during, and after a shock. We show that chemical evolution in the post--shock gas can account for both the deuterium fractionation and the abundance of carbon dioxide relative to water observed in interstellar and cometary ices. Given the pervasiveness of shocks and turbulent motions within molecular clouds, the model presented here offers an alternate theory to grain surface chemistry for the creation of ices in the interstellar medium, ices that may ultimately be incorporated into comets.

We have detected the S(1), S(2), S(3), S(4) and S(5) pure rotational lines of molecular hydrogen toward the outflow source HH 54, using the Short Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO). The observed H2 line ratios indicate the presence of warm molecular gas with an H2 density of at least 100,000 per cubic centimeter and a temperature ~ 650 K in which the ortho- to para-H2 ratio is only 1.2 +/- 0.4, significantly smaller than the equilibrium ratio of 3 expected in gas at that temperature. These observations imply that the measured ortho- to para-H2 ratio is the legacy of an earlier stage in the thermal history of the gas when the gas had reached equilibrium at a temperature < 90 K. Based upon the expected timescale for equilibration, we argue that the non-equilibrium ortho- to para-H2 ratio observed in HH 54 serves as a chronometer that places a conservative upper limit ~ 5000 yr on the period for which the emitting gas has been warm. The S(2)/S(1) and S(3)/S(1) H2 line ratios measured toward HH 54 are consistent with recent theoretical models of Timmermann for the conversion of para- to ortho-H2 behind slow, `C'-type shocks, but only if the preshock ortho- to para-H2 ratio was < 0.2.

We have used a coupled time-dependent chemical and dynamical model to investigate the lifetime of the chemical legacy left in the wake of C-type shocks. We concentrate this study on the chemistry of H2O and O2, two molecules which are predicted to have abundances that are significantly affected in shock-heated gas. Two models are presented: (1) a three-stage model of pre-shock, shocked, and post-shock gas; and (2) a Monte-Carlo cloud simulation where we explore the effects of stochastic shock activity on molecular gas over a cloud lifetime. In agreement with previous studies, we find that shock velocities in excess of 10 km s^-1 are required to convert all of the oxygen not locked in CO into H2O before the gas has an opportunity to cool. For pure gas-phase models the lifetime of the high water abundances, or ``H2O legacy'', in the post-shock gas is 4 - 7 x 10^5 years. Through the Monte Carlo cloud simulation we demonstrate that the time-average abundance of H2O is a sensitive function of the frequency of shocks. Thus we predict that the abundance of H2O and other known outflow tracers can be used to trace the history of shock activity in molecular gas. For gas-grain models we find that the abundance of water-ice on grain surfaces can be quite large and is comparable to that observed in molecular clouds. This offers a possible alternative method to create water mantles without resorting to grain surface chemistry: gas heating and chemical modification due to a C-type shock and subsequent depletion of the gas-phase species onto grain mantles.

We report an upper limit of 9 x 10^{12} cm-2 on the column density of water in the translucent cloud along the line of sight toward HD 154368. This result is based upon a search for the C-X band of water near 1240 \AA carried out using the Goddard High Resolution Spectrograph of the Hubble Space Telescope. Our observational limit on the water abundance together with detailed chemical models of translucent clouds and previous measurements of OH along the line of sight constrain the branching ratio in the dissociative recombination of H_3O+ to form water. We find at the $3\sigma$ level that no more than 30% of dissociative recombinations of H_3O+ can lead to H_2O. The observed spectrum also yielded high-resolution observations of the Mg II doublet at 1239.9 \AA and 1240.4 \AA, allowing the velocity structure of the dominant ionization state of magnesium to be studied along the line of sight. The Mg II spectrum is consistent with GHRS observations at lower spectral resolution that were obtained previously but allow an additional velocity component to be identified.

Using the Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO), we have observed thermal water vapor emission from a roughly circular field of view approximately 75 arc seconds in diameter centered on the Orion BN-KL region. The Fabry-Perot line strengths, line widths, and spectral line shifts observed in eight transitions between 71 and 125 microns show good agreement with models of thermal emission arising from a molecular cloud subjected to a magnetohydrodynamic C-type shock. Both the breadth and the relative strengths of the observed lines argue for emission from a shock rather than from warm quiescent gas in the Orion core. Though one of the eight transitions appears anomalously strong, and may be subject to the effects of radiative pumping, the other seven indicate an H2O/H2 abundance ratio of order 5E-4, and a corresponding gas-phase oxygen-to-hydrogen abundance ratio of order 4E-4. Given current estimates of the interstellar, gas-phase, oxygen and carbon abundances in the solar vicinity, this value is consistent with theoretical shock models that predict the conversion into water of all the gas-phase oxygen that is not bound as CO. The overall cooling provided by rotational transitions of H2O in this region appears to be comparable to the cooling through rotational lines of CO, but is an order of magnitude lower than cooling through H2 emission. However, the model that best fits our observations shows cooling by H2O and CO dominant in that portion of the post-shock region where temperatures are below ~ 800 K and neither vibrational nor rotational radiative cooling by H2 is appreciable.

We report the first detection of interstellar hydrogen fluoride. Using the Long Wavelength Spectrometer (LWS) of the Infrared Space Observatory (ISO), we have detected the 121.6973 micron J = 2 - 1 line of HF in absorption toward the far-infrared continuum source Sagittarius B2. The detection is statistically significant at the 13 sigma level. On the basis of our model for the excitation of HF in Sgr B2, the observed line equivalent width of 1.0 nm implies a hydrogen fluoride abundance of 3E-10 relative to H2. If the elemental abundance of fluorine in Sgr B2 is the same as that in the solar system, then HF accounts for ~ 2% of the total number of fluorine nuclei. We expect hydrogen fluoride to be the dominant reservoir of gas-phase fluorine in Sgr B2, because it is formed rapidly in exothermic reactions of atomic fluorine with either water or molecular hydrogen; thus the measured HF abundance suggests a substantial depletion of fluorine onto dust grains. Similar conclusions regarding depletion have previously been reached for the case of chlorine in dense interstellar clouds. We also find evidence at a lower level of statistical significance (~ 5 sigma) for an emission feature at the expected position of the 4(3,2)-4(2,3) 121.7219 micron line of water. The emission line equivalent width of 0.5 nm for the water feature is consistent with the water abundance of 5E-6 relative to H2 that has been inferred previously from observations of the hot core of Sgr B2.

We present a detailed model for the thermal balance, chemistry, and radiative transfer within quiescent dense molecular cloud cores that contain a central protostar. Large variations in the gas temperature are expected to affect the gas-phase chemistry dramatically; with the predicted H2O abundance varying by more than a factor of 1000 within cloud cores. Based on our predicitions for the thermal and chemical structure of the cores, we have constructed self-consistent radiative transfer models to compute line strengths and profiles for transitions of various isotopomers of CO, H2O, and OI. We predict the high lying transitions of water to be in absorption, and low gain maser emission at 183 GHz. We predict the 63 micron line of OI to be in absorption against the continuum for many sources. Finally, our model can also account successfully for recent ISO observations of absorption in rovibrational transitions of water toward the source AFGL 2591.

We have modeled the gas temperature structure in unstable C-type shocks and obtained predictions for the resultant CO and H2 rotational line emissions, using numerical simulations of the Wardle instability that were carried out by Stone (1997) and that have been described in a companion paper. We obtained results for the gas temperature distribution in -- and H2 and CO line emission from -- shocks of neutral Alfvenic Mach number 10 and velocity 20 or 40 km/sec in which the Wardle instability has saturated. Although the Wardle instability profoundly affects the density structure behind C-type shocks, most of the shock-excited molecular line emission is generated upstream of the region where the strongest effects of the instability are felt. Thus the Wardle instability has a relatively small effect upon the overall gas temperature distribution in -- and the emission line spectrum from -- C-type shocks, at least for the cases that we have considered.

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen within diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photo- dissociation rate of H$_2$ within such clouds was computed using a 26-ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H$_2$ line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 um continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 umu continuum intensity to HI column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be largest. Using a standard chi^2 analysis to assess the goodness of fit of our models, we find (at the 60sigma level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be very much smaller than its extent in the plane-of-the-sky. If D is the distance to G236+39, and given standard assumptions about the rate of grain-catalysed H_2 formation, we find that the cloud has an extent along the line of sight that is 0.9+-0.1 times its mean extent projected onto the plane of the sky; a gas density of (53+-8)(100 pc/D) cm^-3; and is illuminated by a radiation field of (1.1+-0.2) (100 pc/D) times the mean interstellar radiation field estimated by Draine (1978). The derived 100 um emissivity per nucleon is (1.13+-0.06)x10^-20 MJy sr^-1 cm^2.

We have modeled the far-infrared water line emission expected from circumstellar outflows from oxygen-rich late-type stars, as a function of the mass-loss rate and the terminal outflow velocity. For each mass-loss rate and terminal outflow velocity considered, we computed self-consistently the gas density, temperature, outflow velocity, and water abundance as a function of distance from the star. We then used an escape probability method to solve for the equilibrium level populations of 80 rotational states of water and thereby obtained predictions for the luminosity of a large number of far-infrared rotational transitions of water. In common with previous models, our model predicts that water will be copiously produced in the warm circumstellar gas, and that water rotational emission will dominate the radiative cooling. However, our use of a realistic radiative cooling function for water leads to a lower gas temperature than that predicted in previous models. Our predictions for the far-infrared water line luminosities are consequently significantly smaller than those obtained in previous studies. Observations to be carried out by the Infrared Space Observatory will provide a crucial test of the models presented here.