In the Solar System, minor bodies and dust deliver various materials to planetary surfaces. Several exoplanetary systems are known to host inner and outer belts, analogues of the main asteroid belt and the Kuiper belt. We study the possibility that exominor bodies and exodust deliver volatiles and refractories to the exoplanets in the system HR8799 by performing N-body simulations. The model consists of the host star, four giant planets (HR8799 e, d, c, and b), 650000 test particles representing the inner belt, and 1450000 test particles representing the outer belt. Moreover we modelled dust populations that originate from both belts. Within a million years, the two belts evolve towards the expected dynamical structure (also derived in other works), where mean-motion resonances with the planets carve the analogues of Kirkwood gaps. We find that, after this point, the planets suffer impacts by objects from the inner and outer belt at rates that are essentially constant with time, while dust populations do not contribute significantly to the delivery process. We convert the impact rates to volatile and refractory delivery rates using our best estimates of the total mass contained in the belts and their volatile and refractory content. Over their lifetime, the four giant planets receive between $10^{-4}$ and $10^{-3}M_\bigoplus$ of material from both belts. The total amount of delivered volatiles and refractories, ${5\times10^{-3}\textrm{M}_\bigoplus}$, is small compared to the total mass of the planets, $11\times10^{3}\textrm{M}_\bigoplus$. However, if the planets were formed to be volatile-rich, their exogenous enrichment in refractory material may well be significant and observable, for example with JWST-MIRI. If terrestrial planets exist within the snow line of the system, volatile delivery would be an important astrobiological mechanism and may be observable as atmospheric trace gases.

M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution.

The Leiden Atomic and Molecular Database (LAMDA) collects spectroscopic information and collisional rate coefficients for molecules, atoms, and ions of astrophysical and astrochemical interest. We describe the developments of the database since its inception in 2005, and outline our plans for the near future. Such a database is constrained both by the nature of its uses and by the availability of accurate data: we suggest ways to improve the synergies among users and suppliers of data. We summarize some recent developments in computation of collisional cross sections and rate coefficients. We consider atomic and molecular data that are needed to support astrophysics and astrochemistry with upcoming instruments that operate in the mid- and far-infrared parts of the spectrum.

Assuming our Solar System as typical, exomoons may outnumber exoplanets. If their habitability fraction is similar, they would thus constitute the largest portion of habitable real estate in the Universe. Icy moons in our Solar System, such as Europa and Enceladus, have already been shown to possess liquid water, a prerequisite for life on Earth. We intend to investigate under what circumstances small, icy moons may sustain subsurface oceans and thus be "subsurface habitable". We pay specific attention to tidal heating. We made use of a phenomenological approach to tidal heating. We computed the orbit averaged flux from both stellar and planetary (both thermal and reflected stellar) illumination. We then calculated subsurface temperatures depending on illumination and thermal conduction to the surface through the ice shell and an insulating layer of regolith. We adopted a conduction only model, ignoring volcanism and ice shell convection as an outlet for internal heat. In doing so, we determined at which depth, if any, ice melts and a subsurface ocean forms. We find an analytical expression between the moon's physical and orbital characteristics and the melting depth. Since this expression directly relates icy moon observables to the melting depth, it allows us to swiftly put an upper limit on the melting depth for any given moon. We reproduce the existence of Enceladus' subsurface ocean; we also find that the two largest moons of Uranus (Titania & Oberon) could well sustain them. Our model predicts that Rhea does not have liquid water. Habitable exomoon environments may be found across an exoplanetary system, largely irrespective of the distance to the host star. Small, icy subsurface habitable moons may exist anywhere beyond the snow line. This may, in future observations, expand the search area for extraterrestrial habitable environments beyond the circumstellar habitable zone.

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

Context. As a building block for amino acids, formamide (NH$_2$CHO) is an important molecule in astrobiology and astrochemistry, but its formation path in the interstellar medium is not understood well. Aims. We aim to find empirical evidence to support the chemical relationships of formamide to HNCO and H$_2$CO. Methods. We examine high angular resolution (~0.2") Atacama Large Millimeter/submillimeter Array (ALMA) maps of six sources in three high-mass star-forming regions and compare the spatial extent, integrated emission peak position, and velocity structure of HNCO and H$_2$CO line emission with that of NH$_2$CHO by using moment maps. Through spectral modeling, we compare the abundances of these three species. Results. In these sources, the emission peak separation and velocity dispersion of formamide emission is most often similar to HNCO emission, while the velocity structure is generally just as similar to H$_2$CO and HNCO (within errors). From the spectral modeling, we see that the abundances between all three of our focus species are correlated, and the relationship between NH$_2$CHO and HNCO reproduces the previously demonstrated abundance relationship. Conclusions. In this first interferometric study, which compares two potential parent species to NH$_2$CHO, we find that all moment maps for HNCO are more similar to NH$_2$CHO than H$_2$CO in one of our six sources (G24 A1). For the other five sources, the relationship between NH$_2$CHO, HNCO, and H$_2$CO is unclear as the different moment maps for each source are not consistently more similar to one species as opposed to the other.

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

The hydrogen fluoride (HF) molecule is seen in absorption in the interstellar medium (ISM) along many lines of sight. Surprisingly, it is observed in emission toward the Orion Bar, which is an interface between the ionized region around the Orion Trapezium stars and the Orion molecular cloud. We aim to understand the origin of HF emission in the Orion Bar by comparing its spatial distribution with other tracers. We examine three mechanisms to explain the HF emission: thermal excitation, radiative dust pumping, and chemical pumping. We used a Herschel/HIFI strip map of the HF J = 1-0 line, covering 0.5' by 1.5' that is oriented perpendicular to the Orion Bar. We used the RADEX non-local thermodynamic equilibrium (non-LTE) code to construct the HF column density map. We use the Meudon PDR code to explain the morphology of HF. The bulk of the HF emission at 10 km s^-1 emerges from the CO-dark molecular gas that separates the ionization front from the molecular gas that is deeper in the Orion Bar. The excitation of HF is caused mainly by collisions with H2 at a density of 10^5 cm^-3 together with a small contribution of electrons in the interclump gas of the Orion Bar. Infrared pumping and chemical pumping are not important. We conclude that the HF J = 1-0 line traces CO-dark molecular gas. Similarly, bright photodissociation regions associated with massive star formation may be responsible for the HF emission observed toward active galactic nuclei.

Phosphorus nitride (PN) is the P-bearing species with the highest number of detections in star-forming regions. Multi-line studies of the molecule have shown that the excitation temperature of PN is usually lower than the gas kinetic temperature, suggesting that PN is likely in conditions of sub-thermal excitation. We present an analysis of PN which takes the possible sub-thermal excitation conditions into account in a sample of 24 massive star-forming regions. We observed PN (2-1), (3-2), (4-3), and (6-5) with the IRAM-30m and APEX telescopes and detected PN lines in 15 of them. Together with 9 similar sources detected in PN in previous works, we have analysed the largest sample of star-forming regions to date, made of 33 sources with 24 detections in total (among which 13 are new detections). Hence, we have increased the number of star-forming regions detected in PN by more than a factor 2. Our analysis indicates that the PN lines are indeed sub-thermally excited, but well described by a single excitation temperature. We have compared line profiles and fractional abundances of PN and SiO, a typical shock tracer, and found that almost all objects detected in PN have high-velocity SiO wings. Moreover, the SiO and PN abundances with respect to H2 are correlated over several orders of magnitude, and uncorrelated with gas temperature. This clearly shows that the production of PN is strongly linked to the presence of shocked gas, and rules out alternative scenarios based on thermal evaporation from iced grain mantles.

Galaxy mergers and interactions are an integral part of our basic understanding of how galaxies grow and evolve over time. However, the effect that galaxy mergers have on star formation rates (SFR) is contested, with observations of galaxy mergers showing reduced, enhanced and highly enhanced star formation. We aim to determine the effect of galaxy mergers on the SFR of galaxies using statistically large samples of galaxies, totalling over 200\,000, over a large redshift range, 0.0 to 4.0. We train and use convolutional neural networks to create binary merger identifications (merger or non-merger) in the SDSS, KiDS and CANDELS imaging surveys. We then compare the galaxy main sequence subtracted SFR of the merging and non-merging galaxies to determine what effect, if any, a galaxy merger has on SFR. We find that the SFR of merging galaxies are not significantly different from the SFR of non-merging systems. The changes in the average SFR seen in the star forming population when a galaxy is merging are small, of the order of a factor of 1.2. However, the higher the SFR above the galaxy main sequence, the higher the fraction of galaxy mergers. Galaxy mergers have little effect on the SFR of the majority of merging galaxies compared to the non-merging galaxies. The typical change in SFR is less than 0.1~dex in either direction. Larger changes in SFR can be seen but are less common. The increase in merger fraction as the distance above the galaxy main sequence increases demonstrates that galaxy mergers can induce starbursts.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity and high resolution ($\sim 100$ au) observations show a clear hourglass magnetic field morphology centered on the protostellar system, but the central pattern is consistent with a toroidal field indicative of a circumstellar disk, although other interpretations are also possible, including field lines dragged by an equatorial accretion flow into a configuration parallel to the midplane. If a relatively large disk does exist, it would suggest that the magnetic braking catastrophe is averted in this system, not through a large misalignment between the magnetic and rotation axes, but rather through some other mechanisms, such as non-ideal magneto-hydrodynamic effects and/or turbulence. We have also found a relationship of decreasing polarization fractions with intensities and the various slopes of this relationship can be understood as multiple polarization mechanisms and/or depolarization from a changing field morphology. In addition, we found a prominent clumpy depolarization strip crossing the center perpendicular to the bipolar outflow. Moreover, a rough estimate of the magnetic field strength indicates that the field is strong enough to hinder formation of a rotationally supported disk, which is inconsistent with the feature of a central toroidal field.

We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $\mu$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr.

Mergers are an important aspect of galaxy formation and evolution. We aim to test whether deep learning techniques can be used to reproduce visual classification of observations, physical classification of simulations and highlight any differences between these two classifications. With one of the main difficulties of merger studies being the lack of a truth sample, we can use our method to test biases in visually identified merger catalogues. A convolutional neural network architecture was developed and trained in two ways: one with observations from SDSS and one with simulated galaxies from EAGLE, processed to mimic the SDSS observations. The SDSS images were also classified by the simulation trained network and the EAGLE images classified by the observation trained network. The observationally trained network achieves an accuracy of 91.5% while the simulation trained network achieves 65.2% on the visually classified SDSS and physically classified EAGLE images respectively. Classifying the SDSS images with the simulation trained network was less successful, only achieving an accuracy of 64.6%, while classifying the EAGLE images with the observation network was very poor, achieving an accuracy of only 53.0% with preferential assignment to the non-merger classification. This suggests that most of the simulated mergers do not have conspicuous merger features and visually identified merger catalogues from observations are incomplete and biased towards certain merger types. The networks trained and tested with the same data perform the best, with observations performing better than simulations, a result of the observational sample being biased towards conspicuous mergers. Classifying SDSS observations with the simulation trained network has proven to work, providing tantalizing prospects for using simulation trained networks for galaxy identification in large surveys.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

(Abridged) We use HIFI maps of the 987 GHz H2O 2(02)-1(11) emission to measure the sizes and shapes of 19 high-mass protostellar envelopes. To identify infall, we use HIFI spectra of the optically thin C18O 9-8 and H2O-18 1(11)-0(00) lines. The high-J C18O line traces the warm central material and redshifted H2O-18 1(11)-0(00) absorption indicates material falling onto the warm core. We probe small-scale chemical differentiation by comparing H2O 752 and 987 GHz spectra with those of H2O-18. Our measured radii of the central part of the H2O 2(02)-1(11) emission are 30-40% larger than the predictions from spherical envelope models, and axis ratios are <2, which we consider good agreement. For 11 of the 19 sources, we find a significant redshift of the H2O-18 1(11)-0(00) line relative to C18O 9-8. The inferred infall velocities are 0.6-3.2 km/s, and estimated mass inflow rates range from 7e-5 to 2e-2 M0/yr, with the highest mass inflow rates occurring toward the sources with the highest masses, and possibly the youngest ages. The other sources show either expanding motions or H2O-18 lines in emission. The H2O-18 1(11)-0(00) line profiles are remarkably similar to the "differences" between the H2O 2(02)-1(11) and 2(11)-2(02) profiles, suggesting that the H2O-18 line and the H2O 2(02)-1(11) absorption originate just inside the radius where water evaporates from grains, typically 1000-5000 au from the center. In some sources, the H2O-18 line is detectable in the outflow, where no C18O emission is seen. Together, the H2O-18 absorption and C18O emission profiles show that the water abundance around high-mass protostars has at least three levels: low in the cool outer envelope, high within the 100 K radius, and very high in the outflowing gas. Thus, despite the small regions, the combination of lines presented here reveals systematic inflows and chemical information about the outflows.

Starburst galaxies are often found to be the result of galaxy mergers. As a result, galaxy mergers are often believed to lie above the galaxy main sequence: the tight correlation between stellar mass and star formation rate. Here, we aim to test this claim. Deep learning techniques are applied to images from the Sloan Digital Sky Survey to provide visual-like classifications for over 340 000 objects between redshifts of 0.005 and 0.1. The aim of this classification is to split the galaxy population into merger and non-merger systems and we are currently achieving an accuracy of 91.5%. Stellar masses and star formation rates are also estimated using panchromatic data for the entire galaxy population. With these preliminary data, the mergers are placed onto the full galaxy main sequence, where we find that merging systems lie across the entire star formation rate - stellar mass plane.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended Herschel photometry, which reaches $\approx10$ times deeper than the 5$\sigma$ confusion limit in SPIRE. This provides more reliable star formation rates (SFRs), especially for the fainter galaxies, and hence a more reliable MS. We built a pipeline that uses the spectral energy distribution (SED) modelling and fitting tool CIGALE to generate flux density priors in the Herschel SPIRE bands. These priors were then fed into the de-blending tool XID+ to extract flux densities from the SPIRE maps. Multi-wavelength data were combined with the extracted SPIRE flux densities to constrain SEDs and provide stellar mass (M$_{\star}$) and SFRs. These M$_{\star}$ and SFRs were then used to populate the SFR-M$_{\star}$ plane over $0.2\leq z<6.0$. No significant evidence of a high-mass turn-over was found; the best fit is thus a simple two-parameter power law of the form log(SFR)$=\alpha[$log(M$_{\star})-10.5]+\beta$. The normalisation of the power law increases with redshift, rapidly at $z\lesssim1.8$, from $0.58\pm0.09$ at $z\approx0.37$ to $1.31\pm0.08$ at $z\approx1.8$. The slope is also found to increase with redshift, perhaps with an excess around $1.8\leq z<2.9$. The increasing slope indicates that galaxies become more self-similar as redshift increases, implying that the specific SFR of high-mass galaxies increases over $z=0.2$ to $z=6.0$, becoming closer to that of low-mass galaxies. The excess in the slope at $1.8\leq z<2.9$, if present, coincides with the peak of the cosmic star formation history.

In the high-mass star-forming region G35.20-0.74N, small scale (about 800 AU) chemical segregation has been observed in which complex organic molecules containing the CN group are located in a small location. We aim to determine the physical origin of the large abundance difference (about 4 orders of magnitude) in complex cyanides within G35.20-0.74 B, and we explore variations in age, gas and dust temperature, and gas density. We performed gas-grain astrochemical modeling experiments with exponentially increasing (coupled) gas and dust temperature rising from 10 to 500 K at constant H$_2$ densities of 10$^7$, 10$^8$, and 10$^9$ cm$^{-3}$. We tested the effect of varying the initial ice composition, cosmic-ray ionization rate, warm-up time (over 50, 200, and 1000 kyr), and initial (10, 15, and 25 K) and final temperatures (300 and 500 K). Varying the initial ice compositions within the observed and expected ranges does not noticeably affect the modeled abundances indicating that the chemical make-up of hot cores is determined in the warm-up stage. Complex cyanides vinyl and ethyl cyanide (CH$_2$CHCN and C$_2$H$_5$CN, respectively) cannot be produced in abundances (versus H$_2$) greater than 5x10$^{-10}$ for CH$_2$CHCN and 2x10$^{-10}$ for C$_2$H$_5$CN with a fast warm-up time (52 kyr), while the lower limit for the observed abundance of C$_2$H$_5$CN toward source B3 is 3.4x10$^{-10}$. Complex cyanide abundances are reduced at higher initial temperatures and increased at higher cosmic-ray ionization rates. Reproducing the observed abundances toward G35.20-0.74 Core B3 requires a fast warm-up at a high cosmic-ray ionization rate (1x10$^{-16}$ s$^{-1}$) at a high gas density (>10$^9$ cm$^{-3}$). G35.20-0.74 source B3 only needs to be about 2000 years older than B1/B2 for the observed chemical difference to be present. (This abstract has been shortened)

Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured Universe, but sensitivity has been limited because up to now it has not been possible to fly a telescope that is both large and cold. SPICA is a mission concept aimed at taking the next step in mid- and far-infrared observational capability by combining a large and cold telescope with instruments employing state-of-the-art ultra-sensitive detectors. The mission concept foresees a 2.5-meter diameter telescope cooled to below 8 K. With cooling provided by mechanical coolers instead of depending on a limited cryogen supply, the mission lifetime can extend significantly beyond the required three years. SPICA offers instrumentation with spectral resolving powers ranging from R ~50 through 11000 in the 17-230 $\mu$m domain as well as R~28.000 spectroscopy between 12 and 18 $\mu$m. Additionally SPICA will provide efficient 30-37 $\mu$m broad band mapping, and polarimetric imaging in the 100-350 $\mu$m range. SPICA will provide unprecedented spectroscopic sensitivity of ~5 x $10^{-20}$ W/m$^2$ (5$\sigma$/1hr) - at least two orders of magnitude improvement over what has been attained to date. With this exceptional leap in performance, new domains in infrared astronomy will become accessible, allowing us, for example, to unravel definitively galaxy evolution and metal production over cosmic time, to study dust formation and evolution from very early epochs onwards, and to trace the formation history of planetary systems.

Given rapid photodissociation and photodegradation, the recently discovered organics in the Martian subsurface and atmosphere were probably delivered in geologically recent times. Possible parent bodies are C-type asteroids, comets, and interplanetary dust particles (IDPs). The dust infall rate was estimated, using different methods, to be between $0.71$ and $2.96 \times 10^6$ kg/yr (Nesvorny et al., 2011, Borin et al., 2017, Crismani et al., 2017); assuming a carbon content of 10% (Flynn, 1996), this implies an IDP carbon flux of $0.07 - 0.3 \times 10^6$ kg/yr. We calculate for the first time the carbon flux from impacts of asteroids and comets. To this end, we perform dynamical simulations of impact rates on Mars. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We separately add comets and asteroids to the simulations as massless test particles, based on their current orbital elements, yielding Mars impact rates of $4.34\times10^{-3}$ comets/Myr and 3.3 asteroids/Myr. We estimate the global carbon flux on Mars from cometary impacts to be $\sim 0.013 \times 10^{6}$~kg/yr within an order of magnitude, while asteroids deliver $\sim 0.05 \times 10^6$~kg/yr. These values correspond to $\sim 4-19 \%$ and $\sim 17-71 \%$, respectively, of the IDP-borne carbon flux estimated by Nesvorny et al. 2011, Borin et al. 2017 and Crismani et al. 2017. Unlike the spatially homogeneous IDP infall, impact ejecta are distributed locally, concentrated around the impact site. We find organics from asteroids and comets to dominate over IDP-borne organics at distances up to 150~km from the crater center. Our results may be important for the interpretation of in situ detections of organics on Mars.

High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$\alpha$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$\alpha$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$\alpha$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region.

The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $\mu$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $\mu$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy.

The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA.

A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths.

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

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

In the study of high-mass star formation, hot cores are empirically defined stages where chemically rich emission is detected toward a massive YSO. It is unknown whether the physical origin of this emission is a disk, inner envelope, or outflow cavity wall and whether the hot core stage is common to all massive stars. We investigate the chemical make up of several hot molecular cores to determine physical and chemical structure. We use high spectral and spatial resolution Cycle 0 ALMA observations to determine how this stage fits into the formation sequence of a high mass star. We observed the G35.20-0.74N and G35.03+0.35 hot cores at 350 GHz. We analyzed spectra and maps from four continuum peaks (A, B1, B2 and B3) in G35.20, separated by 1000-2000 AU, and one continuum peak in G35.03. We made all possible line identifications across 8 GHz of spectral windows of molecular emission lines and determined column densities and temperatures for as many as 35 species assuming local thermodynamic equilibrium. In comparing the spectra of the four peaks, we find each has a distinct chemical composition expressed in over 400 different transitions. In G35.20, B1 and B2 contain oxygen- and sulfur-bearing organic and inorganic species but few nitrogen-bearing species whereas A and B3 are strong sources of O, S, and N-bearing species (especially those with the CN-bond). CH$_2$DCN is clearly detected in A and B3 with D/H ratios of 8 and 13$\%$, respectively, but is much weaker at B1 and undetected at B2. No deuterated species are detected in G35.03, but similar molecular abundances to G35.20 were found in other species. We also find co-spatial emission of HNCO and NH$_2$CHO in both sources indicating a strong chemical link between the two species. The chemical segregation between N-bearing organic species and others in G35.20 suggests the presence of multiple protostars, surrounded by a disk or torus.

It is not well known what drives the chemistry of a protostellar envelope, in particular the role of the stellar mass and the outflows on its chemical enrichment. We study the chemical structure of NGC 1333 IRAS 4A in order to (i) investigate the influence of the outflows on the chemistry, (ii) constrain the age of our object, (iii) compare it with a typical high-mass protostellar envelope. In our analysis we use JCMT line mapping and HIFI pointed spectra. To study the influence of the outflow on the degree of deuteration, we compare JCMT maps of HCO+ and DCO+ with non-LTE (RADEX) models in a region that spatially covers the outflow activity of IRAS 4A. To study the envelope chemistry, we derive empirical molecular abundance profiles for the observed species using the radiative transfer code (RATRAN) and adopting a 1D dust density/temperature profile from the literature. We compare our best-fit observed abundance profiles with the predictions from the time dependent gas grain chemical code (ALCHEMIC). The CO, HCN, HNC and CN abundance require an enhanced UV field which points towards an outflow cavity. The abundances (wrt H2) are 1 to 2 orders of magnitude lower than those observed in the high mass protostellar envelope (AFGL 2591), while they are found to be similar within factors of a few with respect to CO. Differences in UV radiation may be responsible for such chemical differentiation, but temperature differences seem a more plausible explanation. The CH3OH modeled abundance profile points towards an age of > 4x10^4 yrs for IRAS 4A. The spatial distribution of H2D+ differs from that of other deuterated species, indicating an origin from a foreground colder layer (<20 K). The observed abundances can be explained by passive heating towards the high mass protostellar envelope, while the presence of UV cavity channels become more important toward the low mass protostellar envelope.

Cosmological surveys in the far infrared are known to suffer from confusion. The Bayesian de-blending tool, XID+, currently provides one of the best ways to de-confuse deep Herschel SPIRE images, using a flat flux density prior. This work is to demonstrate that existing multi-wavelength data sets can be exploited to improve XID+ by providing an informed prior, resulting in more accurate and precise extracted flux densities. Photometric data for galaxies in the COSMOS field were used to constrain spectral energy distributions (SEDs) using the fitting tool CIGALE. These SEDs were used to create Gaussian prior estimates in the SPIRE bands for XID+. The multi-wavelength photometry and the extracted SPIRE flux densities were run through CIGALE again to allow us to compare the performance of the two priors. Inferred ALMA flux densities (F$^i$), at 870$\mu$m and 1250$\mu$m, from the best fitting SEDs from the second CIGALE run were compared with measured ALMA flux densities (F$^m$) as an independent performance validation. Similar validations were conducted with the SED modelling and fitting tool MAGPHYS and modified black body functions to test for model dependency. We demonstrate a clear improvement in agreement between the flux densities extracted with XID+ and existing data at other wavelengths when using the new informed Gaussian prior over the original uninformed prior. The residuals between F$^m$ and F$^i$ were calculated. For the Gaussian prior, these residuals, expressed as a multiple of the ALMA error ($\sigma$), have a smaller standard deviation, 7.95$\sigma$ for the Gaussian prior compared to 12.21$\sigma$ for the flat prior, reduced mean, 1.83$\sigma$ compared to 3.44$\sigma$, and have reduced skew to positive values, 7.97 compared to 11.50. These results were determined to not be significantly model dependent. This results in statistically more reliable SPIRE flux densities.

We aim at characterizing the large-scale distribution of H2O in G327.3-0.6, a massive star-forming region made of individual objects in different evolutionary phases. We investigate variations of H2O abundance as function of evolution. We present Herschel continuum maps at 89 and 179 $\mu$m of the whole region and an APEX map at 350 {\mu}m of the IRDC. New spectral HIFI maps toward the IRDC region covering low-energy H2O lines at 987 and 1113 GHz are also presented and combined with HIFI pointed observations of the G327 hot core. We infer the physical properties of the gas through optical depth analysis and radiative transfer modeling. The continuum emission at 89 and 179 {\mu}m follows the thermal continuum emission at longer wavelengths, with a peak at the position of the hot core, a secondary peak in the Hii region, and an arch-like layer of hot gas west of the Hii region. The same morphology is observed in the 1113 GHz line, in absorption toward all dust condensations. Optical depths of ~80 and 15 are estimated and correspond to column densities of 10^15 and 2 10^14 cm-2, for the hot core and IRDC position. These values indicate an H2O to H2 ratio of 3 10^-8 toward the hot core; the abundance of H2O does not change along the IRDC with values of some 10^-8. Infall (over ~ 20") is detected toward the hot core position with a rate of 1-1.3 10^-2 M_sun /yr, high enough to overcome the radiation pressure due to the stellar luminosity. The source structure of the hot core region is complex, with a cold outer gas envelope in expansion, situated between the outflow and the observer, extending over 0.32 pc. The outflow is seen face-on and centered away from the hot core. The distribution of H2O along the IRDC is roughly constant with an abundance peak in the more evolved object. These water abundances are in agreement with previous studies in other massive objects and chemical models.

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

[Abridged] Do some environments favor efficient conversion of molecular gas into stars? To answer this, we need to be able to estimate the H2 mass. Traditionally, this is done using CO and a few assumptions but the Herschel observations in the FIR make it possible to estimate the molecular gas mass independently of CO. Previous attempts to derive gas masses from dust emission suffered from biases. Generally, dust surface densities, HI column densities, and CO intensities are used to derive a gas-to-dust ratio (GDR) and the local CO intensity to H2 column density ratio (XCO), sometimes allowing for an additional CO-dark gas component (Kdark). We tested earlier methods, revealing degeneracies among the parameters, and then used a Bayesian formalism to derive the most likely values for each of the parameters mentioned above as a function of position in the nearby low metallicity spiral galaxy M33. The data are from the IRAM 30m CO(2-1) line, high-resolution HI and Herschel dust continuum observations. Solving for GDR, XCO, and Kdark in macro pixels 500 pc in size, we find that (i) allowing for CO-dark gas significantly improves fits; (ii) Kdark decreases with galactocentric distance; (iii) GDR is slightly higher than initially expected and increases with galactocentric distance; (iv) the total amount of dark gas closely follows the radially decreasing CO emission, as might be expected if the dark gas is H2 where CO is photodissociated. The total amount of H2, including dark gas, yields an average XCO of twice the galactic value of 2e20 cm^-2/(K km/s), 55% of this traced directly through CO. The rather constant fraction of dark gas suggests that there is no large population of diffuse H2 clouds (unrelated to GMCs) without CO emission. Unlike in large spirals, we detect no systematic radial trend in XCO, possibly linked to the absence of a radial decrease in CO line ratios.

Massive stars have a major influence on their environment yet their formation is difficult to study. W43 is a highly luminous galactic massive star forming region at a distance of 5.5 kpc and the MM1 part hosts a very massive dense core (1000 M$_{\odot}$ within 0.05 pc). We present new Herschel HIFI maps of the W43 MM1 region covering the main low-energy water lines at 557, 987, and 1113 GHz, their H$_2^{18}$O counterparts, and other lines such as $^{13}$CO(10-9) and C$^{18}$O(9-8) which trace warm gas. These water lines are, with the exception of line wings, observed in absorption. Herschel SPIRE and JCMT 450 $\mu$m data have been used to make a model of the continuum emission at these wavelengths. Analysis of the maps, and in particular the optical depth maps of each line and feature, shows that a velocity gradient, possibly due to rotation, is present in both the envelope and the protostellar core. Velocities increase in both components from SW to NE, following the general source orientation. While the H$_2$O lines trace essentially the cool envelope, we show that the envelope cannot account for the H$_2^{18}$O absorption, which traces motions close to the protostar. The core has rapid infall, 2.9 kms, as manifested by the H$_2^{18}$O absorption features which are systematically red-shifted with respect to the $^{13}$CO(10-9) emission line which also traces the inner material with the same angular resolution. Some H$_2^{18}$O absorption is detected outside the central core and thus outside the regions expected to be above 100 K - we attribute this to warm gas associated with the other massive dense cores in W43 MM1. Using the maps to identify absorption from cool gas on large scales, we subtract this component to model spectra for the inner envelope. Modeling the new spectra results in a lower water abundance, decreased from $8 10^{-8}$ to $8 10^{-9}$ , with no change in infall rate.

Protostellar evolution, following the formation of the protostar is becoming reasonably well characterized, but the evolution from a prestellar core to a protostar is not well known, although the first hydrostatic core (FHSC) must be a pivotal step. NGC 1333 IRAS 4C is a potentially very young object, that we directly compare with the nearby Class 0 IRAS 4A and IRAS 4B. Observational constraints are provided by spectral imaging from the JCMT Spectral Legacy Survey (330-373 GHz) and continuum and line observations from CARMA. We present integrated intensity and velocity maps of several species, including CO, H2CO and CH3OH. The velocity of an observed outflow, the degree of CO depletion, the deuteration of DCO+/HCO+ and gas kinetic temperatures are observational signatures that we present. We report differences between the three sources in four aspects: a) the kinetic temperature is much lower towards IRAS 4C, b) the line profiles of the detected species show strong outflow activity towards IRAS 4A and IRAS 4B but not towards IRAS 4C, c) the HCN/HNC is < 1 towards IRAS 4C, which confirms the cold nature of the source, d) the degree of CO depletion and the deuteration are the lowest towards the warmest of the sources, IRAS 4B. IRAS 4C seems to be in a different evolutionary state than the IRAS 4A and IRAS 4B sources. We can probably exclude the FHSC stage due to the relatively low Lsmm/Lbol (~6\%) and we investigate the earliest accretion phase of Class 0 stage and the transition between Class 0 to Class I. Our results do not show a consistent scenario for either case, with the major issue being the absence of outflow activity and the cold nature of IRAS 4C. The number of FHSC candidates in Perseus is ~10 times higher than current models predict, which suggests that the lifespan of these objects is greater than 1000 years, possibly due to an accretion rate lower than 4x10^-5 msun/yr.

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

(Abridged) We present Herschel/HIFI spectra of the H2O 1113 GHz and H2O+ 1115 GHz lines toward five nearby prototypical starburst/AGN systems, and OH+ 971 GHz spectra toward three of these. The beam size of 20" corresponds to resolutions between 0.35 and 7 kpc. The observed line profiles range from pure absorption (NGC 4945, M82) to P-Cygni indicating outflow (NGC 253, Arp 220) and inverse P-Cygni indicating infall (Cen A). The similarity of the H2O, OH+, and H2O+ profiles to each other and to HI indicates that diffuse and dense gas phases are well mixed. We estimate column densities assuming negligible excitation (for absorption features) and using a non-LTE model (for emission features), adopting calculated collision data for H2O and OH+, and rough estimates for H2O+. Column densities range from ~10^13 to ~10^15 cm^-2 for each species, and are similar between absorption and emission components, indicating that the nuclear region does not contribute much to the emission in these ground-state lines. The N(H2O)/N(H2O+) ratios of 1.4-5.6 indicate an origin of the lines in diffuse gas, and the N(OH+)/N(H2O+) ratios of 1.6-3.1 indicate a low H2 fraction (~11%) in the gas. Adopting recent Galactic values for the average gas density and the ionization efficiency, we find ionization rates for our sample galaxies of ~3x10^-16 s^-1 which are similar to the value for the Galactic disk, but ~10x below that of the Galactic Center and ~100x below estimates for AGN from excited-state H3O+ lines. We conclude that the ground-state lines of water and its associated ions probe primarily non-nuclear gas in the disks of these centrally active galaxies. Our data thus provide evidence for a decrease in ionization rate by a factor of ~10 from the nuclei to the disks of galaxies, as found before for the Milky Way.

We present Herschel, ALMA Compact Array (ACA), and Caltech Submillimeter Observatory (CSO) observations of the prestellar core in L1689N, which has been suggested to be interacting with a molecular outflow driven by the nearby solar type protostar IRAS 16293-2422. This source is characterized by some of the highest deuteration levels seen in the interstellar medium. The change in the NH2D line velocity and width across the core provides clear evidence of an interaction with the outflow, traced by the high-velocity water emission. Quiescent, cold gas, characterized by narrow line widths is seen in the NE part of the core, while broader, more disturbed line profiles are seen in the W/SW part. Strong N2D+ and ND3 emission is detected with the ACA, extending S/SW from the peak of the single-dish NH2D emission. The ACA data also reveal the presence a compact dust continuum source, with a mean size of ~1100 au, a central density of (1-2) 10^7 cm-3, and a mass of 0.2-0.4 Msun. The dust emission peak is displaced ~5" to the south with respect to the N2D+ and ND3 emission, as well as the single-dish dust continuum peak, suggesting that the northern, quiescent part of the core is characterized by spatially extended continuum emission, which is resolved out by the interferometer. We see no clear evidence of fragmentation in this quiescent part of the core, which could lead to a second generation of star formation, although a weak dust continuum source is detected in this region in the ACA data.

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

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

We aim to understand the contribution of the ionized, atomic and molecular phases of the ISM to the [CII] emission from clouds near the dynamical center and the BCLMP302 HII region in the north of the nearby galaxy M33 at a spatial resolution of 50pc. We combine high resolution [CII] spectra taken with the HIFI spectrometer onboard the Herschel satellite with [CII] Herschel-PACS maps and ground-based observations of CO(2-1) and HI. All data are at a common spatial resolution of 50pc. Typically, the [CII] lines have widths intermediate between the narrower CO(2-1) and broader HI line profiles. We decomposed the [CII] spectra in terms of contribution from molecular and atomic gas detected in CO(2-1) and HI, respectively. We find that the relative contribution of molecular and atomic gas traced by CO(2-1) and HI varies depends mostly on the local physical conditions and geometry. We estimate that 11-60% and 5-34% of the [CII] intensities in the center and in BCLMP302, respectively, arise at velocities showing no CO(2-1) or HI emission and could arise in CO-dark molecular gas. The deduced strong variation in the [CII] emission not associated with CO and HI cannot be explained in terms of differences in A_v, far-ultraviolet radiation field, and metallicity between the two studied regions. Hence the relative amounts of diffuse (CO-dark) and dense molecular gas possibly vary on spatial scales smaller than 50pc. Based on the emission measure observed at radio wavelengths we estimate the contribution of ionized gas at a few positions to lie between 10-25%. The correlations between the intensities of tracers corresponding to the same velocity range as [CII], differ from the correlation derived from PACS data. The results in this paper emphasize the need for velocity-resolved observations to discern the contribution of different components of the ISM to [CII] emission. (abridged)

This paper reviews recent observations of water in Galactic interstellar clouds and nearby galactic nuclei. Two results are highlighted: (1) Multi-line H$_2$O mapping of the Orion Bar shows that the water chemistry in PDRs is driven by photodissociation and -desorption, unlike in star-forming regions. (2) High-resolution spectra of H$_2$O and its ions toward 5 starburst / AGN systems reveal low ionization rates, unlike as found from higher-excitation lines. We conclude that the chemistry of water strongly depends on radiation environment, and that the ionization rates of interstellar clouds decrease by at least 10 between galactic nuclei and disks.

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

Context: Infrared dark clouds are the coldest and densest portions of giant molecular clouds. The most massive ones represent some of the most likely birthplaces for the next generation of massive stars in the Milky Way. Because a strong mid-IR background is needed to make them appear in absorption, they are usually assumed to be nearby. Aims: We use THz absorption spectroscopy to solve the distance ambiguity associated with kinematic distances for the IR-dark clouds in the TOP100 ATLASGAL sample, a flux-limited selection of massive clumps in different evolutionary phases of star formation. Methods: The para-H2O ground state transition at 1113.343 GHz, observed with Herschel/HIFI, was used to investigate the occurrence of foreground absorption along the line of sight directly towards infrared-dark clouds. Additional consistency checks were performed using MALT90 and HiGAL archival data and targeted Mopra and APEX spectroscopic observations. Results: We report the first discovery of five IRDCs in the TOP100 lying conclusively at the far kinematic distance, showing that the mere presence of low-contrast mid-IR absorption is not sufficient to unequivocally resolve the near/far ambiguity in favour of the former. All IRDCs are massive and actively forming high-mass stars; four of them also show infall signatures. Conclusions: We give a first estimate of the fraction of dark sources at the far distance (~11% in the TOP100) and describe their appearance and properties. The assumption that all dark clouds lie at the near distance may lead, in some cases, to underestimating masses, sizes, and luminosities, possibly causing clouds to be missed that will form very massive stars and clusters.

We try to understand the gas heating and cooling in the S 140 star forming region by spatially and spectrally resolving the distribution of the main cooling lines with GREAT/SOFIA. We mapped the fine structure lines of [OI] (63 {\mu}m) and [CII] (158 {\mu}m) and the rotational transitions of CO 13-12 and 16-15 with GREAT/SOFIA and analyzed the spatial and velocity structure to assign the emission to individual heating sources. We measure the optical depth of the [CII] line and perform radiative transfer computations for all observed transitions. By comparing the line intensities with the far-infrared continuum we can assess the total cooling budget and measure the gas heating efficiency. The main emission of fine structure lines in S 140 stems from a 8.3'' region close to the infrared source IRS 2 that is not prominent at any other wavelength. It can be explained by a photon-dominated region (PDR) structure around the embedded cluster if we assume that the [OI] line intensity is reduced by a factor seven due to self-absorption. The external cloud interface forms a second PDR at an inclination of 80-85 degrees illuminated by an UV field of 60 times the standard interstellar radiation field. The main radiation source in the cloud, IRS 1, is not prominent at all in the fine structure lines. We measure line-to-continuum cooling ratios below 10^(-4), i.e. values lower than in any other Galactic source, rather matching the far-IR line deficit seen in ULIRGs. In particular the low intensity of the [CII] line can only be modeled by an extreme excitation gradient in the gas around IRS 1. We found no explanation why IRS 1 shows no associated fine-structure line peak, while IRS 2 does. The inner part of S 140 mimics the far-IR line deficit in ULIRGs thereby providing a template that may lead to a future model.

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

In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $\mu$m and SCUBA at 450 $\mu$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.

The massive and luminous star-forming region W49A is a well known Galactic candidate to probe the physical conditions and chemistry similar to those expected in external starburst galaxies. We aim to probe the physical and chemical structure of W49A on a spatial scale of ~0.8 pc based on the JCMT Spectral Legacy Survey, which covers the frequency range between 330 and 373 GHz. The wide 2x2 arcminutes field and the high spectral resolution of the HARP instrument on JCMT provides information on the spatial structure and kinematics of the cloud. For species where multiple transitions are available, we estimate excitation temperatures and column densities. We detected 255 transitions corresponding to 60 species in the 330-373 GHz range at the center position of W49A. Excitation conditions can be probed for 16 molecules. The chemical composition suggests the importance of shock-, PDR-, and hot core chemistry. Many molecular lines show a significant spatial extent across the maps including high density tracers (e.g. HCN, HNC, CS, HCO+) and tracers of UV-irradiation (e.g. CN and C2H). Large variations are seen between the sub-regions with mostly blue-shifted emission toward the Eastern tail, mostly red-shifted emission toward the Northern clump, and emission peaking around the expected source velocity toward the South-west clump. A comparison of column density ratios of characteristic species observed toward W49A to Galactic PDRs suggests that while the chemistry toward the W49A center is driven by a combination of UV-irradiation and shocks, UV-irradiation dominates for the Northern Clump, Eastern tail, and South-west clump regions. A comparison to a starburst galaxy and an AGN suggests similar C2H, CN, and H2CO abundances (with respect to the dense gas tracer 34CS) between the ~0.8 pc scale probed for W49A and the >1 kpc regions in external galaxies with global star-formation.

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.

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

We aim to understand the rich chemical composition of AFGL 2591, a prototypical isolated high-mass star-forming region. Based on HIFI and JCMT data, the molecular abundances of species found in the protostellar envelope of AFGL 2591 were derived with the Monte Carlo radiative transfer code RATRAN, assuming either constant values or 1D stepwise radial profiles as abundance distributions. The reconstructed 1D abundances were compared with the results of time-dependent gas-grain chemical modeling, considering ages of 10,000 to 50,000 years, cosmic-ray ionization rates of 0.5 to 50 times 10^-16 s^-1, uniformly-sized 0.1-1 micron dust grains, a dust/gas ratio of 1%, and several sets of initial molecular abundances with C/O <1 and >1. Constant abundance models give good fits to the data for CO, CN, CS, HCO+, H2CO, N2H+, C2H, NO, OCS, OH, H2CS, O, C, C+, and CH. Models with an abundance jump at 100 K give good fits to the data for NH3, SO, SO2, H2S, H2O, HCl, and CH3OH. For HCN and HNC, the best models have an abundance jump at 230 K. The time-dependent chemical model can accurately explain abundance profiles of 15 out of these 24 species. The jump-like radial profiles for key species like HCO+, NH3, and H2O are consistent with the outcome of the time-dependent chemical modeling. The best-fit model has a chemical age of 10-50 kyr, a solar C/O ratio of 0.44, and a cosmic-ray ionization rate of 5 x 10^-17 s^-1; grain properties and external UV intensity do not affect the calculated chemical structure much. We thus demonstrate that simple constant or jump-like abundance profiles agree with time-dependent chemical modeling for most key C-, O-, N-, and S-bearing molecules. The main exceptions are species with very few observed transitions (C, O, C+, and CH), with a poorly established chemical network (HCl, H2S) or whose chemistry is strongly affected by surface processes (CH3OH).

Water is a sensitive tracer of physical conditions in star-forming regions because of its large abundance variations between hot and cold regions. We use spectrally resolved observations of rotational lines of H$_2$O and its isotopologs with Herschel/HIFI to constrain the physical conditions of the water emitting region toward the high-mass protostar AFGL2591. We use analytical estimates and rotation diagrams to estimate T$_{ex}$ and column densities of H$_2$O of the envelope, the outflow, and a foreground cloud. Furthermore, we use the non-LTE radiative transfer code to estimate the temperature and volume density of the H$_2$O emitting gas. Assuming LTE, we estimate an T$_{ex}$ of 42 K and a column density of 2$\times$10$^{14}$ cm$^{-2}$ for the envelope and 45 K and 4$\times$10$^{13}$ cm$^{-2}$ for the outflow, in beams of 4" and 30", respectively. Non-LTE models indicate a kinetic temperature of 60-230 K and a volume density of 7$\times$10$^6$-10$^8$ cm$^{-3}$ for the envelope, and a kinetic temperature of 70-90 K and a gas density of 10$^7$-10$^8$ cm$^{-3}$ for the outflow. The o/p ratio of the foreground absorption is 1.9$\pm$0.4, suggesting a low temperature. In contrast, the o/p ratio seen in absorption by the outflow is 3.5$\pm$1.0, as expected for warm gas. The water abundance in the envelope is 10$^{-9}$, similar to the low values found for other high- and low-mass protostars, suggesting that this abundance is constant during the embedded phase of high-mass star formation. The water abundance in the outflow is 10$^{-10}$, which is 10$\times$ lower than in the envelope and in the outflows of high- and low-mass protostars. Since beam size effects can only increase this estimate by a factor of 2, we suggest that the water in the outflow is affected by dissociating UV radiation due to the low extinction in the outflow lobe.

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.

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

Understanding water deuterium fractionation is important for constraining the mechanisms of water formation in interstellar clouds. Observations of HDO and H$_2^{18}$O transitions were carried out towards the high-mass star-forming region G34.26+0.15 with the HIFI instrument onboard the Herschel Space Observatory, as well as with ground-based single-dish telescopes. Ten HDO lines and three H$_2^{18}$O lines covering a broad range of upper energy levels (22-204 K) were detected. We used a non-LTE 1D analysis to determine the HDO/H$_2$O ratio as a function of radius in the envelope. Models with different water abundance distributions were considered in order to reproduce the observed line profiles. The HDO/H$_2$O ratio is found to be lower in the hot core ($\sim$3.5 $\times$ 10$^{-4}$ - 7.5 $\times$ 10$^{-4}$) than in the colder envelope ($\sim$1.0 $\times$ 10$^{-3}$ - 2.2 $\times$ 10$^{-3}$). This is the first time that a radial variation of the HDO/H$_2$O ratio has been found to occur in a high-mass source. The chemical evolution of this source was modeled as a function of its radius and the observations are relatively well reproduced. The comparison between the chemical model and the observations leads to an age of $\sim$10$^5$ years after the infrared dark cloud stage.

Context. The earliest phases of high-mass star formation are poorly understood. Aims. Our goal is to determine the physical conditions and kinematic structure of massive star-forming cloud clumps. Methods. We analyze H$_2$O 557 GHz line profiles observed with HIFI toward four positions in two infrared-dark cloud clumps . By comparison with ground-based C$^{17}$O, N$_2$H$^+$, CH$_3$OH and NH$_3$ line observations, we constrain the volume density and kinetic temperature of the gas and estimate the column density and abundance of H$_2$O and N$_2$H$^+$. Results. The observed water lines are complex with emission and absorption components. The absorption is red shifted and consistent with a cold envelope, while the emission is interpreted as resulting from protostellar outflows. The gas density in the clumps is $\sim$ 10$^7$ cm$^{-3}$. The o-H2O outflow column density is 0.3 to 3.0 10$^{14}$ cm$^{-2}$, the o-H2O absorption column density is between 1.5 10$^{14}$ and 2.6 10$^{15}$ cm$^-2$ with cold o-H2O abundances between 1.5 10$^{-9}$ and 3.1 10$^{-8}$. Conclusions. All clumps have high gas densities ($\sim$ 10$^7$ cm$^{-3}$) and display infalling gas. Three of the four clumps have outflows. The clumps form an evolutionary sequence as probed by H$_2$O N$_2$H$^+$, NH$_3$ and CH$_3$OH. We find that G28-MM is the most evolved, followed by G11-MM then G28-NH3. The least evolved clump is G11-NH3 which shows no sign-posts of star-formation. G11-NH3 is a high-mass pre-stellar core.

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.

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

In order to study the ISM and the interplay between the atomic and molecular components in a low-metallicity environment, we present a complete high angular and spectral resolution map and data cube of the 12CO(2-1) emission from the Local Group galaxy M33. Its metallicity is roughly half-solar, such that we can compare its ISM with that of the Milky Way with the main changes being the metallicity and the gas mass fraction. The data have a 12" angular resolution (50pc) with a spectral resolution of 2.6 km/s and a mean noise level of 20 mK per channel in antenna temperature. A radial cut along the major axis was also observed in the 12CO(1-0) line. The CO data cube and integrated intensity map are optimal when using HI data to define the baseline window and the velocities over which the CO emission is integrated. Great care was taken when building these maps, testing different windowing and baseline options and investigating the effect of error beam pickup. The total CO(2-1) luminosity is 2.8e7 K km/s pc2, following the spiral arms in the inner disk. There is no clear variation in the CO(2-1/1-0) intensity ratio with radius and the average value is roughly 0.8. The total molecular gas mass is estimated, using a N(H2)/Ico(1-0)=4e20 cm-2/(K km/s) conversion factor, to be 3.1e8 Msol. The CO spectra in the cube were shifted to zero velocity by subtracting the velocity of the HI peak from the CO spectra. Hence, the velocity dispersion between the atomic and molecular components is extremely low, independently justifying the use of the HI line in building our maps. Stacking the spectra in concentric rings shows that the CO linewidth and possibly the CO-HI velocity dispersion decrease in the outer disk. Using the CO(2-1) emission to trace the molecular gas, the probability distribution function of the H2 column density shows an excess at high column density above a log normal distribution.

This paper presents the richness of submillimeter spectral features in the high-mass star forming region AFGL 2591. As part of the CHESS (Chemical Herschel Survey of Star Forming Regions) Key Programme, AFGL 2591 was observed by the Herschel/HIFI instrument. The spectral survey covered a frequency range from 480 up to 1240 GHz as well as single lines from 1267 to 1901 GHz (i.e. CO, HCl, NH3, OH and [CII]). Rotational and population diagram methods were used to calculate column densities, excitation temperatures and the emission extents of the observed molecules associated with AFGL 2591. The analysis was supplemented with several lines from ground-based JCMT spectra. From the HIFI spectral survey analysis a total of 32 species were identified (including isotopologues). In spite of the fact that lines are mostly quite week, 268 emission and 16 absorption lines were found (excluding blends). Molecular column densities range from 6e11 to 1e19 cm-2 and excitation temperatures range from 19 to 175 K. One can distinguish cold (e.g. HCN, H2S, NH3 with temperatures below 70 K) and warm species (e.g. CH3OH, SO2) in the protostellar envelope.

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

We investigate the physical and chemical processes at work during the formation of a massive protostar based on the observation of water in an outflow from a very young object previously detected in H2 and SiO in the IRAS 17233-3606 region. We estimated the abundance of water to understand its chemistry, and to constrain the mass of the emitting outflow. We present new observations of shocked water obtained with the HIFI receiver onboard Herschel. We detected water at high velocities in a range similar to SiO. We self-consistently fitted these observations along with previous SiO data through a state-of-the-art, one-dimensional, stationary C-shock model. We found that a single model can explain the SiO and H2O emission in the red and blue wings of the spectra. Remarkably, one common area, similar to that found for H2 emission, fits both the SiO and H2O emission regions. This shock model subsequently allowed us to assess the shocked water column density, N(H2O)=1.2x10^{18} cm^{-2}, mass, M(H2O)=12.5 M_earth, and its maximum fractional abundance with respect to the total density, x(H2O)=1.4x10^{-4}. The corresponding water abundance in fractional column density units ranges between 2.5x10^{-5} and 1.2x10^{-5}, in agreement with recent results obtained in outflows from low- and high-mass young stellar objects.

We present Submillimeter Array (SMA) $\lambda =$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

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

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

We report the first detection of a Galactic source of OH+ line emission: the Orion Bar, a bright nearby photon-dominated region. Line emission is detected over ~1' (0.12 pc), tracing the Bar itself as well as the Southern tip of the Orion Ridge. The line width of ~4 km/s suggests an origin of the OH+ emission close to the PDR surface, at a depth of A_V ~0.3-0.5 into the cloud where most hydrogen is in atomic form. Steady-state collisional and radiative excitation models require unrealistically high OH+ column densities to match the observed line intensity, indicating that the formation of OH+ in the Bar is rapid enough to influence its excitation. Our best-fit OH+ column density of ~1x10^14 cm^-2 is similar to that in previous absorption line studies, while our limits on the ratios of OH+/H2O+ (>~40) and OH+/H3O+ (>~15) are higher than seen before. The column density of OH+ is consistent with estimates from a thermo-chemical model for parameters applicable to the Orion Bar, given the current uncertainties in the local gas pressure and the spectral shape of the ionizing radiation field. The unusually high OH+/H2O+ and OH+/H3O+ ratios are probably due to the high UV radiation field and electron density in this object. In the Bar, photodissociation and electron recombination are more effective destroyers of OH+ than the reaction with H2, which limits the production of H2O+. The appearance of the OH+ lines in emission is the result of the high density of electrons and H atoms in the Orion Bar, since for these species, inelastic collisions with OH+ are faster than reactive ones. In addition, chemical pumping, far-infrared pumping by local dust, and near-UV pumping by Trapezium starlight contribute to the OH+ excitation. Similar conditions may apply to extragalactic nuclei where OH+ lines are seen in emission.

[Context] Two competing models describe the formation of massive stars in objects like the Orion Trapezium. In the turbulent core accretion model, the resulting stellar masses are directly related to the mass distribution of the cloud condensations. In the competitive accretion model, the gravitational potential of the protocluster captures gas from the surrounding cloud for which the individual cluster members compete. [Aims] With high resolution submillimeter observations of the structure, kinematics, and chemistry of the proto-Trapezium cluster W3 IRS5, we aim to determine which mode of star formation dominates. [Methods] We present 354 GHz Submillimeter Array observations at resolutions of 1"-3" (1800-5400 AU) of W3 IRS5. ...... [Results] The observations show five emission peaks (SMM1-5). SMM1 and SMM2 contain massive embedded stars (~20 Msun); SMM3-5 are starless or contain low-mass stars (<8 Msun). The inferred densities are high, >= 10^7 cm^-3, but the core masses are small, 0.2-0.6 Msun. The detected molecular emission reveals four different chemical zones. ...... [Conclusions] The proto-Trapezium cluster W3 IRS5 is an ideal test case to discriminate between models of massive star formation. Either the massive stars accrete locally from their local cores; in this case the small core masses imply that W3 IRS5 is at the very end stages (1000 yr) of infall and accretion, or the stars are accreting from the global collapse of a massive, cluster forming core. We find that the observed masses, densities and line widths observed toward W3 IRS 5 and the surrounding cluster forming core are consistent with the competitive accretion of gas at rates of Macc~10^-4 Msun yr^-1 by the massive young forming stars. ......

Based on astrochemical theory, the complex molecular composition around high-mass YSOs should evolve from the outer envelope in toward the central hot region as a sequence of temperature dependent chemical pathways are activated in ices and in the gas-phase. The resulting complex molecules have been broadly classified into three generations dependent on the temperature (<25, >25, and >100 K) required for formation. We combine IRAM 30m and Submillimeter Array observations to explore the spatial distribution of organic molecules around the high-mass young stellar object NGC 7538 IRS9, whose weak complex molecule emission previously escaped detection, quantifying the emission and abundance profiles of key organic molecules as a function of distance from the central protostar. We find that emission from N-bearing organics and saturated O-bearing organics present large increases in emission around 8000 AU and R<3000 AU, while O-bearing molecules and hydrocarbons do not. The increase in flux from some complex molecules in the envelope, around 8000 AU or 25 K, is consistent with recent model predictions of an onset of complex ice chemistry at 20-30 K. The emission increase for some molecules at R<3000 AU suggests the presence of a weak hot core, where thermal ice evaporation and hot gas-phase chemistry drives the chemistry. Complex organics thus form at all radii and temperatures around this protostar, but the composition changes dramatically as the temperature increases, which is used to constrain the chemical generation(s) to which different classes of molecule belong.

To understand the origin of water line emission and absorption during high-mass star formation, we decompose high-resolution Herschel-HIFI line spectra toward 19 high-mass star-forming regions into three distinct physical components. Protostellar envelopes are usually seen as narrow absorptions or emissions in the H2O 1113 and 1669 GHz ground-state lines, the H2O 987 GHz excited-state line, and the H2O-18 1102 GHz ground-state line. Broader features due to outflows are usually seen in absorption in the H2O 1113 and 1669 GHz lines, in 987 GHz emission, and not seen in H2O-18, indicating a low column density and a high excitation temperature. The H2O 1113 and 1669 GHz spectra show narrow absorptions by foreground clouds along the line of sight, which have a low column density and a low excitation temperature, although their H2O ortho/para ratios are close to 3. The intensities of the H2O 1113 and 1669 GHz lines do not show significant trends with luminosity, mass, or age. In contrast, the 987 GHz line flux increases with luminosity and the H2O-18 line flux decreases with mass. Furthermore, appearance of the envelope in absorption in the 987 GHz and H2O-18 lines seems to be a sign of an early evolutionary stage. We conclude that the ground state transitions of H2O trace the outer parts of the envelopes, so that the effects of star formation are mostly noticeable in the outflow wings. These lines are heavily affected by absorption, so that line ratios of H2O involving the ground states must be treated with caution. The average H2O abundance in high-mass protostellar envelopes does not change much with time. The 987 GHz line appears to be a good tracer of the mean weighted dust temperature of the source, which may explain why it is readily seen in distant galaxies.

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

We aim to reveal the gas energetics in the circumstellar environment of the prototypical high-mass protostellar object AFGL2591 using space-based far-infrared observations of linear rotor molecules. Rotational spectral line signatures of CO, HCO+, CS, HCN and HNC from a 490-1240 GHz survey with Herschel/HIFI, complemented by ground-based JCMT and IRAM 30m spectra, cover transitions with E(up)/k between 5 and ~300 K (750K for 12C16O, using selected frequency settings up to 1850 GHz). The resolved spectral line profiles are used to separate and study various kinematic components. The line profiles show two emission components, the widest and bluest of which is attributed to an approaching outflow and the other to the envelope. We find evidence for progressively more redshifted and wider line profiles from the envelope gas with increasing energy level, qualitatively explained by residual outflow contribution picked up in the systematically decreasing beam size. Integrated line intensities for each species decrease as E(up)/k increases from <50 to 700K. We constrain the following: n(H2)~10^5-10^6 cm^-3 and T~60-200K for the outflow gas; T=9-17K and N(H2)~3x10^21 cm^-2 for a known foreground absorption cloud; N(H2)<10^19 cm^-2 for a second foreground component. Our spherical envelope radiative transfer model systematically underproduces observed line emission at E(up)/k > 150 K for all species. This indicates that warm gas should be added to the model and that the model's geometry should provide low optical depth pathways for line emission from this warm gas to escape, for example in the form of UV heated outflow cavity walls viewed at a favorable inclination angle. Physical and chemical conditions derived for the outflow gas are similar to those in the protostellar envelope, possibly indicating that the modest velocity (<10 km/s) outflow component consists of recently swept-up gas.

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv)

We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios.

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

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

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

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.

Most studies of high-mass star formation focus on massive luminous clumps, but the physical properties of their larger scale environment are poorly known. In this work, we aim at characterising the effects of clustered star formation and feedback of massive stars on the surrounding medium by studying the distribution of warm gas through mid-J 12CO and 13CO data. We present APEX 12CO(6-5), (7-6), 13CO(6-5), (8-7) and HIFI 13CO(10-9) maps of the star forming region G327.36-0.6. We infer the physical properties of the emitting gas on large scales through a LTE analysis, while we apply a more sophisticated LVG approach on selected positions. Maps of all lines are dominated in intensity by the PDR around the Hii region G327.3-0.5. Mid-J 12CO emission is detected over the whole extent of the maps with excitation temperatures ranging from 20K up to 80K in the gas around the Hii region, and H2 column densities from few 10^21 cm-2 in the inter-clump gas to 3 10^22 cm-2 towards the hot core G327.3-0.6. The warm gas (traced by 12 and 13CO(6-5) emission) is only a small percentage (10%) of the total gas in the IRDC, while it reaches values up to 35% of the total gas in the ring surrounding the Hii region. The 12CO ladders are qualitatively compatible with PDR models for high density gas, but the much weaker than predicted 13CO emission suggests that it comes from a large number of clumps along the line of sight. All lines are detected in the inter-clump gas when averaged over a large region with an equivalent radius of 50" (~0.8pc), implying that the mid-J 12CO and 13CO inter-clump emission is due to high density components with low filling factor. Finally, the detection of the 13CO(10-9) line allows to disentangle the effects of gas temperature and gas density on the CO emission, which are degenerate in the APEX observations alone.

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

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

This work presents high spectral resolution observations of the \CII\ line at 158 \micron, one of the major cooling lines of the interstellar medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an \HII\ region far north (3.3 kpc) in the disk of M 33, the \CII\ and CO line profiles show similar velocities within $0.5 \kms$, while the \HI\ line velocities are systematically shifted towards lower rotation velocities by $\sim 5\kms$. Observed at the same $12"$ angular resolution, the \CII\ lines are broader than those of CO by about 50% but narrower than the \HI\ lines. The \CII\ line to far-infrared continuum ratio suggests a photoelectric heating efficiency of 1.1%. The data, together with published models indicate a UV field $G_0 \sim 100$ in units of the solar neighborhood value, a gas density $n_H \sim 1000 \cc$, and a gas temperature $T\sim 200$ K. Adopting these values, we estimate the C$^+$ column density to be $N_{C^+} \approx 1.3 \times 10^{17} \cmt$. The \CII\ emission comes predominantly from the warm neutral region between the \HII\ region and the cool molecular cloud behind it. From published abundances, the inferred C$^+$ column corresponds to a hydrogen column density of $N_H \sim 2 \times 10^{21} \cmt$. The CO observations suggest that $N_H = 2 N_{H_2} \sim 3.2 \times 10^{21} \cmt$ and 21cm measurements, also at $12"$ resolution, yield $N_\HI \approx 1.2 \times 10^{21} \cmt$ within the \CII\ velocity range. Thus, some H$_2$ not detected in CO must be present, in agreement with earlier findings based on the SPIRE 250 -- 500 $\mu$m emission.

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 study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).

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.

[Context] Recent detections of disks around young high-mass stars support the idea of massive star formation through accretion rather than coalescence, but the detailed kinematics in the equatorial region of the disk candidates is not well known, which limits our understanding of the accretion process. [Aims] This paper explores the kinematics of the gas around a young massive star with millimeter-wave interferometry to improve our understanding of the formation of massive stars though accretion. [Methods] We use Plateau de Bure interferometric images to probe the environment of the nearby (~1 kpc) and luminous (~20000 Lsun) high-mass (10-16 Msun) young star AFGL 2591-VLA3 in continuum and in lines of HDO, H218O and SO2 in the 115 and 230 GHz bands. Radiative transfer calculations are employed to investigate the kinematics of the source. [Results] At ~0.5" (500 AU) resolution, the line images clearly resolve the velocity field of the central compact source (diameter of ~ 800 AU) and show linear velocity gradients in the northeast-southwest direction. Judging from the disk-outflow geometry, the observed velocity gradient results from rotation and radial expansion in the equatorial region of VLA3. Radiative transfer calculations suggest that the velocity field is consistent with sub-Keplerian rotation plus Hubble-law like expansion. The line profiles of the observed molecules suggest a layered structure, with HDO emission arising from the disk mid-plane, H218O from the warm mid-layer, and SO2 from the upper disk. [Conclusions] We propose AFGL 2591-VLA3 as a new massive disk candidate, with peculiar kinematics. The rotation of this disk is sub-Keplerian, probably due to magnetic braking, while the stellar wind may be responsible for the expansion of the disk. The expansion motion [...]

This paper contains a summary of the results from the first years of observations with the HIFI instrument onboard ESA's Herschel space observatory. The paper starts by outlining the goals and possibilities of far-infrared and submillimeter astronomy, the limitations of the Earth's atmosphere, and the scientific scope of the Herschel-HIFI mission. The presentation of science results from the mission follows the life cycle of gas in galaxies as grouped into five themes: Structure of the interstellar medium, First steps in interstellar chemistry, Formation of stars and planets, Solar system results and Evolved stellar envelopes. The HIFI observations paint a picture where the interstellar medium in galaxies has a mixed, rather than a layered structure; the same conclusion may hold for protoplanetary disks. In addition, the HIFI data show that exchange of matter between comets and asteroids with planets and moons plays a large role. The paper concludes with an outlook to future instrumentation in the far-infrared and submillimeter wavelength ranges.

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

The star formation rates in starburst galaxies are orders of magnitude higher than in local star-forming regions, and the origin of this difference is not well understood. We use sub-mm spectral line maps to characterize the physical conditions of the molecular gas in the luminous Galactic star-forming region W49A and compare them with the conditions in starburst galaxies. We probe the temperature and density structure of W49A using H_2CO and HCN line ratios over a 2'x2' (6.6x6.6 pc) field with an angular resolution of 15" (~0.8 pc) provided by the JCMT Spectral Legacy Survey. We analyze the rotation diagrams of lines with multiple transitions with corrections for optical depth and beam dilution, and estimate excitation temperatures and column densities. Comparing the observed line intensity ratios with non-LTE radiative transfer models, our results reveal an extended region (about 1'x1', equivalent to ~3x3 pc at the distance of W49A) of warm (> 100 K) and dense (>10^5 cm^-3) molecular gas, with a mass of 2x10^4 - 2x10^5 M_Sun (by applying abundances derived for other regions of massive star-formation). These temperatures and densities in W49A are comparable to those found in clouds near the center of the Milky Way and in starburst galaxies. The highly excited gas is likely to be heated via shocks from the stellar winds of embedded, O-type stars or alternatively due to UV irradiation, or possibly a combination of these two processes. Cosmic rays, X-ray irradiation and gas-grain collisional heating are less likely to be the source of the heating in the case of W49A.

The H3+ ion plays a key role in the chemistry of dense interstellar gas clouds where stars and planets are forming. The low temperatures and high extinctions of such clouds make direct observations of H3+ impossible, but lead to large abundances of H2D+ and D2H+ which are very useful probes of the early stages of star and planet formation. Maps of H2D+ and D2H+ pure rotational line emission toward star-forming regions show that the strong deuteration of H3+ is the result of near-complete molecular depletion of CNO-bearing molecules onto grain surfaces, which quickly disappears as cores warm up after stars have formed. In the warmer parts of interstellar gas clouds, H3+ transfers its proton to other neutrals such as CO and N2, leading to a rich ionic chemistry. The abundances of such species are useful tracers of physical conditions such as the radiation field and the electron fraction. Recent observations of HF line emission toward the Orion Bar imply a high electron fraction, and we suggest that observations of OH+ and H2O+ emission may be used to probe the electron density in the nuclei of external galaxies.

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

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

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.

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

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

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

The estimation of molecular abundances in interstellar clouds from spectroscopic observations requires radiative transfer calculations, which depend on basic molecular input data. This paper reviews recent developments in the fields of molecular data and radiative transfer. The first part is an overview of radiative transfer techniques, along with a "road map" showing which technique should be used in which situation. The second part is a review of measurements and calculations of molecular spectroscopic and collisional data, with a summary of recent collisional calculations and suggested modeling strategies if collision data are unavailable. The paper concludes with an overview of future developments and needs in the areas of radiative transfer and molecular data.

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

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

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

The understanding of the formation process of massive stars (>8 Msun) is limited, due to theoretical complications and observational challenges. We investigate the physical structure of the large-scale (~10^4-10^5 AU) molecular envelope of the high-mass protostar AFGL2591 using spectral imaging in the 330-373 GHz regime from the JCMT Spectral Legacy Survey. Out of ~160 spectral features, this paper uses the 35 that are spatially resolved. The observed spatial distributions of a selection of six species are compared with radiative transfer models based on a static spherically symmetric structure, a dynamic spherical structure, and a static flattened structure. The maps of CO and its isotopic variations exhibit elongated geometries on scales of ~100", and smaller scale substructure is found in maps of N2H+, o-H2CO, CS, SO2, CCH, and methanol lines. A velocity gradient is apparent in maps of all molecular lines presented here, except SO, SO2, and H2CO. We find two emission peaks in warm (Eup~200K) methanol separated by 12", indicative of a secondary heating source in the envelope. The spherical models are able to explain the distribution of emission for the optically thin H13CO+ and C34S, but not for the optically thick HCN, HCO+, and CS, nor for the optically thin C17O. The introduction of velocity structure mitigates the optical depth effects, but does not fully explain the observations, especially in the spectral dimension. A static flattened envelope viewed at a small inclination angle does slightly better. We conclude that a geometry of the envelope other than an isotropic static sphere is needed to circumvent line optical depth effects. We propose that this could be achieved in envelope models with an outflow cavity and/or inhomogeneous structure at scales smaller than ~10^4 AU. The picture of inhomogeneity is supported by observed substructure in at least six species.

We use the H3O+ molecule to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium. Using the JCMT, we have observed the 3+_2 - 2-_2 364 GHz line of p-H3O+ towards the centers of seven active galaxies. We have detected p-H3O+ towards IC342, NGC253, NGC1068, NGC4418, and NGC6240. Upper limits were obtained for IRAS15250 and Arp299. We find large H3O+ abundances (N(H3O+)/N(H2)>10^{-8}) in all detected galaxies apart from in IC342 where it is about one order of magnitude lower. We note, however, that uncertainties in N(H3O+) may be significant due to lack of definite information on source size and excitation. We furthermore compare the derived N(H3O+) with N(HCO+) and find that the H3O+ to HCO+ column density ratio is large in NGC1068 (24), moderate in NGC4418 and NGC253 (4-5), slightly less than unity in NGC6240 (0.7) and lowest in IC342 (0.2-0.6). We compare our results with models of X-ray and photon dominated regions (XDRs and PDRs). For IC342 we find that a starburst PDR chemistry can explain the observed H3O+ abundance. For the other galaxies, the large H3O+ columns are generally consistent with XDR models. In particular for NGC1068 the elevated N(H3O+)/N(HCO+) ratio suggests a low column density XDR. For NGC4418 however, large HC3N abundances are inconsistent with the XDR interpretation. An alternative possibility is that H3O+ forms through H2O evaporating off dust grains and reacting with HCO+ in warm, dense gas. This scenario could also potentially fit the results for NGC253. Further studies of the excitation and distribution of H3O+ - as well as Herschel observations of water abundances - will help to further constrain the models.

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.

This paper reviews the first results of observations of H2O line emission with Herschel-HIFI towards high-mass star-forming regions, obtained within the WISH guaranteed time program. The data reveal three kinds of gas-phase H2O: `cloud water' in cold tenuous foreground clouds, `envelope water' in dense protostellar envelopes, and `outflow water' in protostellar outflows. The low H2O abundance (1e-10 -- 1e-9) in foreground clouds and protostellar envelopes is due to rapid photodissociation and freeze-out on dust grains, respectively. The outflows show higher H2O abundances (1e-7 -- 1e-6) due to grain mantle evaporation and (probably) neutral-neutral reactions.

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

W49A is a giant molecular cloud which harbors some of the most luminous embedded clusters in the Galaxy. However, the explanation for this starburst-like phenomenon is still under debate. Methods. We investigated large-scale Spitzer mid-infrared images together with a Galatic Ring Survey 13CO J = 1-0 image, complemented with higher resolution (~ 11 arcsec) 13CO J = 2-1 and C18O J = 2-1 images over a ~ 15 x 13 pc^2 field obtained with the IRAM 30m telescope. Two expanding shells have been identified in the mid-infrared images, and confirmed in the position-velocity diagrams made from the 13CO J = 2-1 and C18O J = 2-1 data. The mass of the averaged expanding shell, which has an inner radius of ~ 3.3 pc and a thickness of ~ 0.41 pc, is about 1.9 x 10^4 M*. The total kinetic energy of the expanding shells is estimated to be ~ 10^49 erg which is probably provided by a few massive stars, whose radiation pressure and/or strong stellar winds drive the shells. The expanding shells are likely to have a common origin close to the two ultracompact Hii regions (source O and source N), and their expansion speed is estimated to be ~ 5 km/s, resulting in an age of ~ 3-7 x 10^5 years. In addition, on larger (~ 35 x 50 pc^2) scales, remnants of two gas ejections have been identified in the 13CO J = 1 - 0 data. Both ejections seem to have the same center as the expanding shells with a total energy of a few times 10^50 erg. The main driving mechanism for the gas ejections is unclear, but likely related to the mechanism which triggers the starburst in W49A.

The astrochemistry of the HnO+ (n=1..3) ions is important as the main gas-phase formation route for water, and as tracer of the interstellar ionization rate by cosmic rays and other processes. While interstellar H3O+ has been known since the early 1990's, interstellar OH+ and H2O+ have only recently been detected using the Herschel space observatory and also from the ground. This paper reviews detections of HnO+ toward external galaxies and compares with ground-based work. The similarities and differences of the HnO+ chemistry within the Galaxy and beyond are discussed. Special attention is given to the low H2O/H3O+ ratio in M82 of only 3.3, suggesting rapid H2O photodissociation, and the high apparent OH+ and H2O+ abundances in Mrk 231, suggesting radiative excitation and/or formation pumping. Photodissociation rates for H3O+ and collisional cross-sections for OH+ and H2O+ with H, He and electrons are needed to test these interpretations.

The HCN, HCO+, and HNC molecules are commonly used as tracers of dense star-forming gas in external galaxies, but such observations are spatially unresolved. Reliably inferring the properties of galactic nuclei and disks requires detailed studies of sources whose structure is spatially resolved. We compare the spatial distributions and abundance ratios of HCN, HCO+, and HNC in W49A, the most massive and luminous star-forming region in the Galactic disk, based on maps of a 2' (6.6 pc) field at 14" (0.83 pc) resolution of the J=4-3 transitions of HCN, H13CN, HC15N, HCO+, H13CO+, HC18O+ and HNC. The kinematics of the molecular gas in W49A appears complex, with a mixture of infall and outflow motions. Both the line profiles and comparison of the main and rarer species show that the main species are optically thick. Two 'clumps' of infalling gas appear to be at ~40 K, compared to ~100 K at the source centre, and may be ~10x denser than the rest of the outer cloud. Chemical modelling suggests that the HCN/HNC ratio probes the current gas temperature, while the HCN/HCO+ ratio and the deuterium fractionation were set during an earlier, colder phase of evolution. The data suggest that W49A is an appropriate analogue of an extragalactic star forming region. Our data show that the use of HCN/HNC/HCO+ line ratios as proxies for the abundance ratios is incorrect for W49A, suggesting the same for galactic nuclei. Our observed isotopic line ratios such as H13CN/H13CO+ approach our modeled abundance ratios quite well in W49A. The 4-3 lines of HCN and HCO+ are much better tracers of the dense star-forming gas in W49A than the 1-0 lines. Our observed HCN/HNC and HCN/HCO+ ratios in W49A are inconsistent with homogeneous PDR or XDR models, indicating that irradiation hardly affects the gas chemistry in W49A. Overall, the W49A region appears to be a useful template for starburst galaxies.

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

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

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

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

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

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

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.

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

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

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.

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

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

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

We derive the dense core structure and the water abundance in four massive star-forming regions which may help understand the earliest stages of massive star formation. We present Herschel-HIFI observations of the para-H2O 1_11-0_00 and 2_02-1_11 and the para-H2-18O 1_11-0_00 transitions. The envelope contribution to the line profiles is separated from contributions by outflows and foreground clouds. The envelope contribution is modelled using Monte-Carlo radiative transfer codes for dust and molecular lines (MC3D and RATRAN), with the water abundance and the turbulent velocity width as free parameters. While the outflows are mostly seen in emission in high-J lines, envelopes are seen in absorption in ground-state lines, which are almost saturated. The derived water abundances range from 5E-10 to 4E-8 in the outer envelopes. We detect cold clouds surrounding the protostar envelope, thanks to the very high quality of the Herschel-HIFI data and the unique ability of water to probe them. Several foreground clouds are also detected along the line of sight. The low H2O abundances in massive dense cores are in accordance with the expectation that high densities and low temperatures lead to freeze-out of water on dust grains. The spread in abundance values is not clearly linked to physical properties of the sources.

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

Despite its major role in the evolution of the interstellar medium, the formation of high-mass stars (M > 10 Msol) is still poorly understood. Two types of massive star cluster precursors, the so-called Massive Dense Cores (MDCs), have been observed, which differ in their mid-infrared brightness. The origin of this difference is not established and could be the result of evolution, density, geometry differences, or a combination of these. We compare several molecular tracers of physical conditions (hot cores, shocks) observed in a sample of mid-IR weak emitting MDCs with previous results obtained in a sample of exclusively mid-IR bright MDCs. The aim is to understand the differences between these two types of object. We present single-dish observations of HDO, H2O-18, SO2 and CH3OH lines at lambda = 1.3 - 3.5 mm. We study line profiles and estimate abundances of these molecules, and use a partial correlation method to search for trends in the results. The detection rates of thermal emission lines are found to be very similar between mid-IR quiet and bright objects. The abundances of H2O, HDO (1E-13 to 1E-9 in the cold outer envelopes), SO2 and CH3OH differ from source to source but independently of their mid-IR flux. In contrast, the methanol class I maser emission, a tracer of outflow shocks, is found to be strongly anti-correlated with the 12 micron source brightnesses. The enhancement of the methanol maser emission in mid-IR quiet MDCs may indicate a more embedded nature. Since total masses are similar between the two samples, we suggest that the matter distribution is spherical around mid-IR quiet sources but flattened around mid-IR bright ones. In contrast, water emission is associated with objects containing a hot molecular core, irrespective of their mid-IR brightness. These results indicate that the mid-IR brightness of MDCs is an indicator of their evolutionary stage.

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

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

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

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

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.

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.

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

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

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

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

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

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

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

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

Water is a key molecule in the star formation process, but its spatial distribution in star-forming regions is not well known. We study the distribution of dust continuum and H2O and 13CO line emission in DR21, a luminous star-forming region with a powerful outflow and a compact HII region. Herschel-HIFI spectra near 1100 GHz show narrow 13CO 10-9 emission and H2O 1(11)-0(00) absorption from the dense core and broad emission from the outflow in both lines. The H2O line also shows absorption by a foreground cloud known from ground-based observations of low-J CO lines. The dust continuum emission is extended over 36" FWHM, while the 13CO and H2O lines are confined to ~24" or less. The foreground absorption appears to peak further North than the other components. Radiative transfer models indicate very low abundances of ~2e-10 for H2O and ~8e-7 for 13CO in the dense core, and higher H2O abundances of ~4e-9 in the foreground cloud and ~7e-7 in the outflow. The high H2O abundance in the warm outflow is probably due to the evaporation of water-rich icy grain mantles, while the H2O abundance is kept down by freeze-out in the dense core and by photodissociation in the foreground cloud.

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.

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

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.

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

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

Photon-dominated regions (PDRs) are expected to show a layered structure in molecular abundances and emerging line emission, which is sensitive to the physical structure of the region as well as the UV radiation illuminating it. We aim to study this layering in the Orion Bar, a prototypical nearby PDR with a favorable edge-on geometry. We present new maps of 2 by 2 arcminute fields at 14-23 arcsecond resolution toward the Orion Bar in the SO 8_8-9_9, H2CO 5_(1,5)-4_(1,4), 13CO 3-2, C2H 4_(9/2)-3_(7/2) and 4_(7/2)-3_(5/2), C18O 2-1 and HCN 3-2 transitions. The data reveal a clear chemical stratification pattern. The C2H emission peaks close to the ionization front, followed by H2CO and SO, while C18O, HCN and 13CO peak deeper into the cloud. A simple PDR model reproduces the observed stratification, although the SO emission is predicted to peak much deeper into the cloud than observed while H2CO is predicted to peak closer to the ionization front than observed. In addition, the predicted SO abundance is higher than observed while the H2CO abundance is lower than observed. The discrepancies between the models and observations indicate that more sophisticated models, including production of H2CO through grain surface chemistry, are needed to quantitatively match the observations of this region.

The far-infrared spectra of circumstellar envelopes around various oxygen-rich stars were observed using the ISO Long Wavelength Spectrometer (LWS). We have examined high signal-to-noise ISO LWS observations of the luminous supergiant star, VY CMa, with the aim of identifying all of the spectral lines. By paying particular attention to water lines, we aim to separate the lines due to other species, in particular, to prepare for forthcoming observations that will cover the same spectral range using Herschel PACS and at higher spectral resolution using Herschel HIFI and SOFIA. We have developed a fitting method to account for blended water lines using a simple weighting scheme to distribute the flux. We have applied this approach to several other stars which we compare with VY CMa. We present line fluxes for the unblended H2O and CO lines, and present detections of several possible nu_2=1 vibrationally excited water lines. We also identify blended lines of OH, one unblended and several blended lines of NH3, and one possible detection of H3O+.

We have observed the rotational ground-state (J = 1-0) transitions of DCO+, HN13C and DNC with the IRAM 30m telescope toward the dark cloud LDN 1512 which has exceptionally narrow lines permitting hyperfine splitting to be resolved in part. The measured splittings of 50-300 kHz are used to derive nuclear quadrupole and spin-rotation parameters for these species. The measurements are supplemented by high-level quantum-chemical calculations using coupled-cluster techniques and large atomic-orbital basis sets. We find eQq = +151.12 (400) kHz and C_I = -1.12 (43) kHz for DCO+, eQq = 272.5 (51) kHz for HN13C, and eQq(D) = 265.9 (83) kHz and eQq(N) = 288.2 (71) kHz for DNC. The numbers for DNC are consistent with previous laboratory data, while our constants for DCO+ are somewhat smaller than previous results based on astronomical data. For both DCO+ and DNC, our results are more accurate than previous determinations. Our results are in good agreement with the corresponding best theoretical estimates. We also derive updated rotational constants for HN13C: B = 43545.6000 (47) MHz and D = 93.7 (20) kHz. The hyperfine splittings of the DCO+, DNC and HN13C J = 1-0 lines range over 0.47-1.28 km/s, which is comparable to typical line widths in pre-stellar cores and to systematic gas motions on ~1000 AU scales in protostellar cores. We present tabular information to allow inclusion of the hyperfine splitting in astronomical data interpretation. The large differences in the 14N quadrupole parameters of DNC and HN13C have been traced to differences in the vibrational corrections caused by significant non-rigidity of these molecules, particularly along the bending coordinate.

With the incorporation of high-J molecular lines, we aim to constrain the physical conditions of the dense gas in the central region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. We used the James Clerk Maxwell Telescope to observe the J=4-3 transition of HCN, HNC, and HCO+, as well as the CN N_J=2_{5/2}-1_{3/2} and N_J=3_{5/2}-2_{5/2}, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line intensity ratios and radiative transfer models. We find that the bulk emission of HCN, HNC, CN, and the high-J HCO+ emerge from dense gas n(H_2)>=10^5 cm^-3). However, the low-J HCO+ lines (dominating the HCO+ column density) trace less dense (n(H_2)<10^5 cm^-3) and colder (T_K<=20 K) gas, whereas the high-J HCO+ emerges from warmer (>30 K) gas than the other molecules. The HCO+ J=4-3 line intensity, compared with the lower transition lines and with the HCN J=4-3 line, support the influence of a local XDR environment. The estimated N(CN)/N(HCN)~1-4 column density ratios are indicative of an XDR/AGN environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks.

We present a survey of the ortho-H2D+(1_{1,0}-1_{1,1}) line toward a sample of 10 starless cores and 6 protostellar cores, carried out at the Caltech Submillimeter Observatory. The high diagnostic power of this line is revealed for the study of the chemistry, and the evolutionary and dynamical status of low-mass dense cores. The line is detected in 7 starless cores and in 4 protostellar cores. N(ortho-H2D+) ranges between 2 and 40x10^{12} cm^{-2} in starless cores and between 2 and 9x10^{12} cm^{-2} in protostellar cores. The brightest lines are detected toward the densest and most centrally concentrated starless cores, where the CO depletion factor and the deuterium fractionation are also largest. The large scatter observed in plots of N(ortho-H2D+) vs. the observed deuterium fractionation and vs. the CO depletion factor is likely to be due to variations in the ortho-to-para (o/p) ratio of H2D+ from >0.5 for T_{kin} < 10 K gas in pre-stellar cores to ~0.03 (consistent with T_{kin} ~15 K for protostellar cores). The two Ophiuchus cores in our sample also require a relatively low o/p ratio (~0.3). Other parameters, including the cosmic-ray ionization rate, the CO depletion factor (or, more in general, the depletion factor of neutral species), the volume density, the fraction of dust grains and PAHs also largely affect the ortho-H2D+ abundance. The most deuterated and H2D+-rich objects (L429, L1544, L694-2 and L183) are reproduced by chemical models of centrally concentrated (central densties ~10^{6} cm^{-3}) cores with chemical ages between 10^4 and 10^6 yr. Upper limits of the para-H3O+ (1_1- -2_1+) and para-D2H+ (1_{1,0}-1_{0,1}) lines are also given. (Abridged)

Although most young massive stars appear to be part of multiple systems, it is poorly understood how this multiplicity influences the formation of massive stars. The high-mass star-forming region W3 IRS5 is a prime example of a young massive cluster where the cluster center is resolved into multiple subsources at cm and infrared wavelengths, a potential proto-Trapezium system. The region W3 IRS5 was mapped with the PdBI at 1.4mm and 3.4mm in the AB configurations, observing shock-tracing SiO and SO_2 emission. In the continuum we detect five sources, one of them for the first time, while counterparts were detected in the NIR, MIR or at radio wavelengths for the remaining four sources. Three of the detected sources are within the inner 2100AU, where the protostellar number density exceeds 10^6 protostars pc^-3 assuming spherical symmetry. Lower limits for the circumstellar masses of the detected sources were calculated, although they were strongly affected by the spatial filtering of the interferometer. However, the projected separations of the sources ranging between ~750 and ~4700AU indicate a multiple, Trapezium-like system. We detected five molecular outflows in SiO, two of them nearly in the line of sight direction, which allowed us to see the collapsing protostars in the NIR through the cavities carved by the outflows. The SO_2 velocity structure indicates a rotating, bound system, and we find tentative signatures of converging flows as predicted by the gravoturbulent star formation and converging flow theories. The obtained data strongly indicate that the clustered environment has a major influence on the formation of high-mass stars; however, our data do not clearly allow us to distinguish whether the ongoing star-forming process follows a monolithic collapse or a competitive accretion mechanism.

The study of physical and chemical properties of massive protostars is critical to better understand the evolutionary sequence which leads to the formation of high-mass stars. IRAS 18151-1208 is a nearby massive region (d = 3kpc, L ~ 20000 Lsun) which splits into three cores: MM1, MM2 and MM3 (separated by 1'-2'). We aim at (1) studying the physical and chemical properties of the individual MM1, MM2 and MM3 cores; (2) deriving their evolutionary stages; (3) using these results to improve our view of the evolutionary sequence of massive cores. The region was observed in the CS, C34S, H2CO, HCO+, H13CO+, and N2H+ lines at mm wavelengths with the IRAM 30m and Mopra telescopes. We use 1D and 2D modeling of the dust continuum to derive the density and temperature distributions, which are then used in the RATRAN code to model the lines and constrain the abundances of the observed species. All the lines were detected in MM1 and MM2. MM3 shows weaker emission, or even is undetected in HCO+ and all isotopic species. MM2 is driving a newly discovered CO outflow and hosts a mid-IR-quiet massive protostar. The abundance of CS is significantly larger in MM1 than in MM2, but smaller than in a reference massive protostar such as AFGL2591. In contrast the N2H+ abundance decreases from MM2 to MM1, and is larger than in AFGL2591. Both MM1 and MM2 host an early phase massive protostar, but MM2 (and mid-IR-quiet sources in general) is younger and more dominated by the host protostar than MM1 (mid-IR-bright). The MM3 core is probably in a pre-stellar phase. We find that the N2H+/C34S ratio varies systematically with age in the massive protostars for which the data are available. It can be used to identify young massive protostars.

We discuss some implications of our recent detection of extragalactic H3O+: the location of the gas in M82, the origin of energetic radiation in M82, and the possible feedback effects of star formation on the cosmic ray flux in galaxies.

This review discusses recent results on the astrochemistry of (mostly high-mass) star-forming regions. After an introduction on the use of chemistry in astrophysics and some basic concepts of astrochemistry, specific results are presented. Highlighted areas are the use of chemistry in the search for massive circumstellar disks, the interaction of molecular clouds with cosmic rays, and the feedback effects of protostellar irradiation on the parent molecular cloud. The review concludes with a discussion of future observational opportunities.

The H3O+ molecule probes the oxygen chemistry and the ionization rate of dense circumnuclear gas in galaxies. In particular, recent H3O+ observations show variations in the cosmic-ray ionization rate by factors of $>$10 within our Galaxy. Using the JCMT, we have observed the 364 GHz line of p-H3O+ in the centers of M82 and Arp 220. In Arp 220, the line profile suggests that the emission originates in the Western nucleus. In M82, both the eastern molecular peak and the circumnuclear region contribute to the emission. The derived column densities, abundances, and H3O+ / H2O ratios indicate ionization rates similar to or even exceeding that in the Galactic Center. Model calculations of the chemistry of irradiated molecular gas indicate a likely origin of this high ionization rate in the extended, evolved starburst of M82. In contrast, irradiation by X-rays from the AGN disk is the most likely model for Arp 220.

Context: observations of water lines are a sensitive probe of the geometry, dynamics and chemical structure of dense molecular gas. The launch of Herschel with on board HIFI and PACS allow to probe the behaviour of multiple water lines with unprecedented sensitivity and resolution. Aims: we investigate the diagnostic value of specific water transitions in high-mass star-forming regions. As a test case, we apply our models to the AFGL2591 region. Results: in general, for models with a constant water abundance, the ground state lines, i.e., 1_(10)-1_(01), 1_(11)-0_(00), and 2_(12)-1_(01), are predicted in absorption, all the others in emission. This behaviour changes for models with a water abundance jump profile in that the line profiles for jumps by a factor of ~10-100 are similar to the line shapes in the constant abundance models, whereas larger jumps lead to emission profiles. Asymmetric line profiles are found for models with a cavity outflow and depend on the inclination angle. Models with an outflow cavity are favoured to reproduce the SWAS observations of the 1_(10)-1_(01) ground-state transition. PACS spectra will tell us about the geometry of these regions, both through the continuum and through the lines. Conclusions: it is found that the low-lying transitions of water are sensitive to outflow features, and represent the excitation conditions in the outer regions. High-lying transitions are more sensitive to the adopted density and temperature distribution which probe the inner excitation conditions. The Herschel mission will thus be very helpful to constrain the physical and chemical structure of high-mass star-forming regions such as AFGL2591.

The large quantity and high quality of modern radio and infrared line observations require efficient modeling techniques to infer physical and chemical parameters such as temperature, density, and molecular abundances. We present a computer program to calculate the intensities of atomic and molecular lines produced in a uniform medium, based on statistical equilibrium calculations involving collisional and radiative processes and including radiation from background sources. Optical depth effects are treated with an escape probability method. The program is available on the World Wide Web at http://www.sron.rug.nl/~vdtak/radex/index.shtml . The program makes use of molecular data files maintained in the Leiden Atomic and Molecular Database (LAMDA), which will continue to be improved and expanded. The performance of the program is compared with more approximate and with more sophisticated methods. An Appendix provides diagnostic plots to estimate physical parameters from line intensity ratios of commonly observed molecules. This program should form an important tool in analyzing observations from current and future radio and infrared telescopes.

This paper presents an overview of several modeling tools for analyzing molecular line observations at submillimeter wavelengths. These tools are already proving to be very valuable for the interpretation of data from current telescopes, and will be indispensable for data obtained with ALMA. The tools are: (1) the Leiden Atomic and Molecular DAtabase (LAMDA), a collection of spectroscopic data and collisional excitation rates; (2) RADEX, an on-line and off-line program to calculate non-LTE excitation and emission from a homogeneous medium, based on the escape probability approximation; (3) RATRAN, an accelerated Monte Carlo program to solve molecular excitation and radiative transfer in spherical and cylindrical symmetry. The paper presents examples of how to use these tools in conjunction with existing data reduction packages to quantitatively interpret submillimeter single-dish and interferometric observations. The described tools are publically available at http://www.strw.leidenuniv.nl/~moldata . The paper concludes with a discussion of future needs in the fields of molecular data and radiative transfer.

Submillimeter lines of CN, NO, CO+ and SO+, and upper limits on SH+ and N2O are observed with the James Clerk Maxwell Telescope in two high-mass and up to nine low-mass young stellar objects and compared with chemical models. Constant fractional abundances derived from radiative transfer modeling of the line strengths are x(CN) ~ a few x 10^{-11}-10^{-8}, x(NO) ~ 10^{-9}-10^{-8} and x(CO+) ~ 10^{-12}-10^{-10}. SO+ has abundances of a few x 10^{-11} in the high-mass objects and upper limits of ~ 10^{-12}-10^{-11} in the low-mass sources. All abundances are up to 1-2 orders of magnitude higher if the molecular emission is assumed to originate mainly from the inner region (< 1000 AU) of the envelope. For high-mass sources, the CN, SO+ and CO+ abundances and abundance ratios are best explained by an enhanced far-ultraviolet (FUV) field impacting gas at temperatures of a few hundred K. The observed column densities require that this region of enhanced FUV has scales comparable to the observing beam, such as in a geometry in which the enhanced FUV irradiates outflow walls. For low-mass sources, the required temperatures within the FUV models of T > 300 K are much higher than found in models, so that an X-ray enhanced region close to the protostar (r < 500 AU) is more plausible. The observed abundances imply X-ray fluxes for the Class 0 objects of L_X ~ 10^{29}-10^{31} erg s^{-1}, comparable to those observed from low-mass Class I protostars. Spatially resolved data are needed to clearly distinguish the effects of FUV and X-rays for individual species.

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

The cosmic-ray ionization rate (zeta) of dense molecular clouds is a key parameter for their dynamics and chemistry. Variations of zeta are well established, but it is unclear if these are related to source column density or to Galactic location. Using the APEX telescope, we have mapped the 364 GHz line of H3O+ in the Sgr B2 region and observed the 307 GHz line at selected positions. With the IRAM 30-m telescope we have observed the 203 GHz line of H2O-18 at the same positions. Strong H3O+ emission is detected over a ~3x2 pc region, indicating H3O+ column densities of 10^15 - 10^16 cm^-2 in an 18" beam. The H3O+ abundance of ~3 x 10^-9 and H3O+/H2O ratio of ~1/50 in the Sgr B2 envelope are consistent with models with zeta ~4 x 10^-16 s^-1, 3x lower than derived from H3+ observations toward Sgr A, but 10x that of local dense clouds. The ionization rates of interstellar clouds thus seem to be to first order determined by the ambient cosmic-ray flux, while propagation effects cause a factor of ~3 decrease from diffuse to dense clouds.

The presence of H2D+ in dense cloud cores underlies ion-molecule reactions that strongly enhance the deuterium fractionation of many molecular species. We determine the H2D+ abundance in one starless core, Barnard 68, that has a particularly well established physical, chemical, and dynamical structure. We observed the ortho-H2D+ ground-state line 1_10-1_11, the N2H+ J=4-3 line, and the H13CO+ 4-3 line with the APEX telescope. We report the probable detection of the o-H2D+ line at an intensity Tmb=0.22 +- 0.08 K and exclusively thermal line width, and find only upper limits to the N2H+ 4-3 and H13CO+ 4-3 intensities. Within the uncertainties in the chemical reaction rates and the collisional excitation rates, chemical model calculations and excitation simulations reproduce the observed intensities and that of o-H2D+ in particular.

Aims: Our goal is to set constraints on the evolutionary state of the dense core Cha-MMS1 in the Chamaeleon I molecular cloud. Methods: We analyze molecular line observations carried out with the new submillimeter telescope APEX. We look for outflow signatures around the dense core and probe its chemical structure, which we compare to predictions of models of gas-phase chemistry. We also use the public database of the Spitzer Space Telescope (SST) to compare Cha-MMS1 with the two Class 0 protostars IRAM 04191 and L1521F, which are at the same distance. Results: We measure a large deuterium fractionation for N2H+ (11 +/- 3 %), intermediate between the prestellar core L1544 and the very young Class 0 protostar L1521F. It is larger than for HCO+ (2.5 +/- 0.9 %), which is probably the result of depletion removing HCO+ from the high-density inner region. Our CO(3-2) map reveals the presence of a bipolar outflow driven by the Class I protostar Ced 110 IRS 4 but we do not find evidence for an outflow powered by Cha-MMS1. We also report the detection of Cha-MMS1 at 24, 70 and 160 microns by the instrument MIPS of the SST, at a level nearly an order of magnitude lower than IRAM 04191 and L1521F. Conclusions: Cha-MMS1 appears to have already formed a compact object, either the first hydrostatic core at the very end of the prestellar phase, or an extremely young protostar that has not yet powered any outflow, at the very beginning of the Class 0 accretion phase.

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

The H3+ and H2D+ ions are important probes of the physical and chemical conditions in regions of the interstellar medium where new stars are forming. This paper reviews how observations of these species and of heavier ions such as HCO+ and H3O+ can be used to derive chemical and kinematic properties of nearby pre-stellar cores, and the cosmic-ray ionisation rate toward more distant regions of high-mass star formation. Future prospects in the field are outlined at the end.

Single-dish spectra and interferometric maps of (sub)mm lines of H2O-18 and HDO are used to study the chemistry of water in eight regions of high-mass star formation. The spectra indicate HDO excitation temperatures of ~110 K and column densities in an 11'' beam of ~2x10^14 cm^-2 for HDO and ~2x10^17 cm^-2 for H2O, with the N(HDO)/N(H2O) ratio increasing with decreasing temperature. Simultaneous observations of CH3OH and SO2 indicate that 20-50% of the single-dish line flux arises in the molecular outflows of these objects. The outflow contribution to the H2O-18 and HDO emission is estimated to be 10-20%. Radiative transfer models indicate that the water abundance is low (~10^-6) outside a critical radius corresponding to a temperature in the protostellar envelope of ~100 K, and `jumps' to H2O/H2 ~ 10^-4 inside this radius. This value corresponds to the observed abundance of solid water and together with the derived HDO/H2O abundance ratios of ~1/1000 suggests that the origin of the observed water is evaporation of grain mantles. This idea is confirmed in the case of AFGL 2591 by interferometer observations of HDO, H2O-18 and SO2 lines, which reveal compact (~800 AU) emission with a systematic velocity gradient. This size is similar to that of the 1.3 mm continuum towards AFGL 2591, from which we estimate a mass of ~0.8 M0, or ~5% of the mass of the central star. We speculate that we may be observing a circumstellar disk in an almost face-on orientation.

This paper reviews the chemistry of star-forming regions, with an emphasis on the formation of high-mass stars. We first outline the basic molecular processes in dense clouds, their implementation in chemical models, and techniques to measure molecular abundances. Then, recent observational, theoretical and laboratory developments are reviewed on the subjects of hot molecular cores, cosmic-ray ionization, depletion and deuteration, and oxygen chemistry. The paper concludes with a summary of outstanding problems and future opportunities.

Radio observations in the ammonia inversion lines of four comets, C/2001 A2 (LINEAR), 153P/Ikeya-Zhang, C/2001 Q4 (NEAT) and C/2002 T7 (LINEAR), were performed at the Effelsberg 100-m Radio Telescope during their respective close approaches to Earth. None of the four lowest energy metastable lines (J,K=J),J=1--4, could be detected in these comets. We derive the following 3 sigma upper bounds on the NH_3 production rate, and c omparing to the corresponding water production rates, percentage NH_3 abundances relative to H_2O: Q(NH_3) < 1.9x10^26 s^-1 (0.63%) for C/2001 A2 (LINEAR), Q(NH_3) < 2.7x10^26 s^-1 (0.13%) for C/2001 Q4 (NEAT), Q(NH_3) < 2.3x10^27 s^-1 (0.74%) for C/2002 T7 (LINEAR) and Q(NH_3) <= 6.3x10^26 s^-1 (0.63%) for Comet 153P/Ikeya-Zhang. At 0.74% or less, the ammonia-to-water ratios are factors of 2 below the value for C/1995 O1 (Hale-Bopp) and 1P/Halley, suggesting chemical diversity between comets. The 18-cm lines of OH were clearly detected in the two comets observed during the 2004 campaign, thereby validating the cometary ephemerides.

Velocity profiles of ground state lines of H2D+, HC18O+ and N2H+, observed previously with the CSO and IRAM 30m telescopes, are modeled with a Monte Carlo radiative transfer program to study the temperature, density and velocity structure of the pre-stellar core LDN 1544. The H2D+ line is double-peaked like that of the other ions, but previous models that fit the HC18O+ and N2H+ profiles are found not to fit the H2D+ data. Matching the H2D+ observations requires at least three modifications to the model at small radii: (1) the density profile must continue to rise inward and not flatten off toward the center; (2) the gas temperature must be nearly constant and not drop inwards significantly; (3) the infall velocity must increase inward, in a fashion intermediate between `quasi-static' (ambipolar diffusion) and `fully dynamic' (Larson-Penston) collapse. The C18O emission indicates a chemical age of <~0.1 Myr. The effects of a flattened structure and rotation on the line profiles are shown to be unimportant, at least on the scales probed by single-dish telescopes. Alternatively, the H2D+ profile is affected by absorption in the outer layers of the core, if gas motions in these layers are sufficiently small.

We use spectra and maps of NH2D, ND2H, and ND3, obtained with the CSO, IRAM 30m and Arecibo telescopes, to study deuteration processes in dense cores. The data include the first detection of the hyperfine structure of ND2H. The emission of ND2H and ND3 does not seem to peak at the positions of the embedded protostars, but instead at offset positions, where outflow interactions may occur. A constant ammonia fractionation ratio in star-forming regions is generally assumed to be consistent with an origin on dust grains. However, in the pre-stellar cores studied here, the fractionation varies significantly when going from NH3 to ND3. We present a steady state model of the gas-phase chemistry for these sources, which includes passive depletion onto dust grains and multiply saturated deuterated species up to five deuterium atoms (e.g. CD5+). The observed column density ratios of all four ammonia isotopologues are reproduced within a factor of 3 for a gas temperature of 10 K. We also predict that deuterium fractionation remains significant at temperatures up to 20 K. ND and NHD, which have rotational transitions in the submillimeter domain are predicted to be abundant.

This paper reports 43 GHz imaging of the high-mass protostars W33A, AFGL 2591 and NGC 7538 IRS9 at 0.04'' and 0.6'' resolution. In each case, weak (~mJy), compact (~100 AU) emission is detected, which has an elongated shape (axis ratio ~3). For AFGL 2591 and NGC 7538 IRS9, the emission is single-peaked, while for the highest luminosity source, W33A, a `mini-cluster' of three sources is detected. The derived sizes, flux densities, and broad-band radio spectra of the sources support recent models where the initial expansion of HII regions around very young O-type stars is prevented by stellar gravity. In these models, accretion flows onto high-mass stars originate in large-scale molecular envelopes and become ionized close to the star. These models reproduce our observations of ionized gas as well as the structure of the molecular envelopes of these sources on 10^3--10^4 AU scales derived previously from single-dish submillimeter continuum and line mapping. For AFGL 2591, the 43 GHz flux density is also consistent with dust emission from a disk seen in near-infrared `speckle' images. However, the alignment of the 43 GHz emission with the large-scale molecular outflow argues against an origin in a disk for AFGL 2591 and NGC 7538 IRS9. In contrast, the outflow from W33A does not appear to be collimated. Together with previously presented case studies of W3 IRS5 and AFGL 2136, our results indicate that the formation of stars and stellar clusters with luminosities up to ~10^5 L0 proceeds through accretion and produces collimated outflows as in the solar-type case, with the `additional feature' that the accretion flow becomes ionized close to the star. Above ~10^5 L0, clusters of HII regions appear, and outflows are no longer collimated, possibly as the result of mergers of protostars or pre-stellar cores.

Observations at mid-infrared (4.8-17.65 um) and radio (0.7-1.3 cm) wavelengths are used to constrain the structure of the high-mass star-forming region W3 IRS5 on 0.1'' (200 AU) scales. Two bright mid-infrared sources are detected, as well as diffuse emission. The bright sources have associated compact radio emission and probably are young high-mass stars. The measured sizes and estimated temperatures indicate that these sources together can supply the observed far-infrared luminosity. However, an optically thick radio source with a possible mid-infrared counterpart may also contribute significant luminosity; if so, it must be extremely deeply embedded. The infrared colour temperatures of 350-390 K and low radio brightness suggest gravitational confinement of the H II regions and ongoing accretion at a rate of a few 10^-8 M0/yr or more. Variations in the accretion rate would explain the observed radio variability. The low estimated foreground extinction suggests the existence of a cavity around the central stars, perhaps blown by stellar winds. At least three radio sources without mid-infrared counterparts appear to show proper motions of ~100 km/s, and may be deeply embedded young runaway OB stars, but more likely are clumps in the ambient material which are shock-ionized by the OB star winds.

Atomic and molecular data for the transitions of a number of astrophysically interesting species are summarized, including energy levels, statistical weights, Einstein A-coefficients and collisional rate coefficients. Available collisional data from quantum chemical calculations and experiments are extrapolated to higher energies. These data, which are made publically available through the WWW at http://www.strw.leidenuniv.nl/~moldata, are essential input for non-LTE line radiative transfer programs. An online version of a computer program for performing statistical equilibrium calculations is also made available as part of the database. Comparisons of calculated emission lines using different sets of collisional rate coefficients are presented. This database should form an important tool in analyzing observations from current and future (sub)millimetre and infrared telescopes.

We present deep images of dust continuum emission at 450, 800, and 850 micron of the dark cloud LDN 1689N which harbors the low-mass young stellar objects (YSOs) IRAS 16293-2422A and B (I16293A and I16293B) and the cold prestellar object I16293E. Toward the positions of I16293A and E we also obtained spectra of CO-isotopomers and deep submillimeter observations of chemically related molecules with high critical densities. To I16293A we report the detection of the HDO 1_01 - 0_00 and H2O 1_10 - 1_01 ground-state transitions as broad self-reversed emission profiles with narrow absorption, and a tentative detection of H2D+ 1_10 - 1_11. To I16293E we detect weak emission of subthermally excited HDO 1_01 - 0_00. Based on this set of submillimeter continuum and line data we model the envelopes around I16293A and E. The density and velocity structure of I16293A is fit by an inside-out collapse model, yielding a sound speed of a=0.7 km/s, an age of t=(0.6--2.5)e4 yr, and a mass of 6.1 Msun. The density in the envelope of I16293E is fit by a radial power law with index -1.0+/-0.2, a mass of 4.4 Msun, and a constant temperature of 16K. These respective models are used to study the chemistry of the envelopes of these pre- and protostellar objects. The [HDO]/[H2O] abundance ratio in the warm inner envelope of I16293A of a few times 1e-4 is comparable to that measured in comets. This supports the idea that the [HDO]/[H2O] ratio is determined in the cold prestellar core phase and conserved throughout the formation process of low-mass stars and planets.

We report 5-43 GHz radio observations of the CRL 2136 region at 0.6-6 arcsec resolution. We detect weak (mJy intensity) radio emission from the deeply embedded high-mass protostar IRS 1, which has an optically thick spectrum up to frequencies of 22 GHz, flattening at higher frequencies, which might be explained by emission from a jet. Water maser mapping shows that the strong emission observed redshifted relative to the systemic velocity is spatially coincident with the optically thick continuum emission. The H2O maser emission from this object (and others we know of) seems to have a different origin than most of these masers, which are frequently tracing bipolar high-velocity outflows. Instead, the CRL 2136 H2O emission arises in the close circumstellar environment of the protostar (within 1000 AU). We speculate that most of it is excited in the hot, dense infalling gas after the accretion shock, although this cannot explain all the H2O emission. An accretion shock nature for the continuum emission seems unlikely.

This paper explores the sensitivity of line profiles of H2D+, HCO+ and N2H+, observed towards the center of L1544, to various kinematic and chemical parameters. The total width of the H2D+ line can be matched by a static model and by models invoking ambipolar diffusion and gravitational collapse. The derived turbulent line width is b=0.15 km/s for the static case and <~ 0.05 km/s for the collapse case. However, line profiles of HC18O+ and N2H+ rule out the static solution. The double-peaked H2D+ line shape requires either infall speeds in the center that are much higher than predicted by ambipolar diffusion models, or a shell-type distribution of H2D+, as is the case for HCO+ and N2H+. At an offset of ~20 arcsec from the dust peak, the H2D+ abundance drops by a factor of ~5.

The sulphur chemistry in nine regions in the earliest stages of high-mass star formation is studied through single-dish submillimeter spectroscopy. The line profiles indicate that 10-50% of the SO and SO2 emission arises in high-velocity gas, either infalling or outflowing. For the low-velocity gas, excitation temperatures are 25 K for H2S, 50 K for SO, H2CS, NS and HCS+, and 100 K for OCS and SO2, indicating that most observed emission traces the outer parts (T<100 K) of the molecular envelopes, except high-excitation OCS and SO2 lines. Abundances in the outer envelopes, calculated with a Monte Carlo program, using the physical structures of the sources derived from previous submillimeter continuum and CS line data, are ~10^-8 for OCS, ~10^-9 for H2S, H2CS, SO and SO2, and ~10^-10 for HCS+ and NS. In the inner envelopes (T>100 K) of six sources, the SO2 abundance is enhanced by factors of ~100-1000. This region of hot, abundant SO2 has been seen before in infrared absorption, and must be small, <~ 0.2 arcsec (180 AU radius). The derived abundance profiles are consistent with models of envelope chemistry which invoke ice evaporation at T~100$ K. Shock chemistry is unlikely to contribute. A major sulphur carrier in the ices is probably OCS, not H2S as most models assume. The source-to-source abundance variations of most molecules by factors of ~10 do not correlate with previous systematic tracers of envelope heating. Without observations of H2S and SO lines probing warm (>~ 100 K) gas, sulphur-bearing molecules cannot be used as evolutionary tracers during star formation.

This review covers hot cores in the context of high-mass star formation. After giving an overview of chemical processes and diversity during high-mass star formation, it reviews the `warm envelope' phase which probably precedes the formation of hot cores. Some recent determinations of the cosmic-ray ionization rate are discussed, as well as recent evidence for hot cores around low-mass stars. Routes for future hot core research are outlined.

We have detected the 372 GHz line of ortho-H2D+ towards the pre-stellar core L1544. The strongest emission (T_mb ~ 1 K) occurs at the peak of the millimeter continuum emission, while measurements at offset positions indicate that H2D+ is confined within ~ 20 arcsec, where CO is highly depleted. The derived H2D+ abundance of ~ 10^{-9} is comparable with previous estimates of the electron abundance in the core, which suggests that H2D+ is the main molecular ion in the central 20 arcsec (2800 AU) of L1544. This confirms the expectations that H2D+ is dramatically enhanced in gas depleted of molecules other than H2. The measured abundance even exceeds the present model predictions by about a factor ten. One possibility is that all CNO-bearing neutral species, including atomic oxygen, are almost completely (> 98%) frozen within a radius of ~2800 AU.

We present submillimeter observations of SO, SO2, H2S, H2CS, OCS, NS and HCS+ toward nine massive young stars. The outflow contributes ~50% to the SO and SO2 emission in 15-20'' beams, more than for CS, where it is 10%. The SO2 abundance increases from dark cloud levels in the outer envelope (T<100 K) to levels seen in hot cores and shocks in the inner envelope (T>100 K). Molecular abundances are consistent with a model of ice evaporation in an envelope with gradients in temperature and density for a chemical age of ~30000 yr. The high observed abundance of OCS, the fact that T_ex (OCS) >> T_ex (H2S), and the data on solid OCS and H2S all suggest that the major sulphur carrier in grain mantles is OCS rather than H2S. For most other sulphur-bearing molecules, the source-to-source abundance variations by factors of up to 10 do not correlate with previously established evolutionary trends in temperature tracers. These species probe the chemically inactive outer envelope. Our data set does not constrain the abundances of H2S and SO in the inner envelope, which, together with SO2, are required to use sulphur as a clock.

Research on interstellar ND3 is reviewed and updated from the discovery papers. Results from observations of a dozen sources at centimeter and submillimeter wavelengths are presented. The two data sets are consistent, but do not constrain the excitation conditions in the ND3-emitting gas. The column density ratios of NH3:NH2D, NH2D:ND2H and ND2H:ND3, observed in similar sized beams, are ~10, and present a challenge for both the gas-phase and the grain-surface chemistry scenarios of deuterium fractionation. The role of shocks and of CO depletion in deuterium chemistry is discussed.

Molecular line data for the rotational transitions of a number of astrophysically interesting species have been collected and are made publically available through the www. These data, including energy levels, statistical weights, Einstein A-coefficients and collisional rate coefficients, are essential input for non-LTE molecular radiative transfer programs. A computer program for performing statistical equilibrium calculations is made publically available for use as well. This database should form an important tool in analyzing observations from current and future infrared and (sub)millimetre telescopes.

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

We present detailed thermal and gas-phase chemical models for the envelope of the massive star-forming region AFGL 2591. Time- and space-dependent chemistry are used to study the physical structure proposed by van der Tak et al. (1999; 2000), and the chemical evolution of this region. The model is compared with observations for 29 species covering a wide range of conditions within the source. Taking appropriate care when comparing models with both emission and absorption measurements, we find that the majority of the chemical structure can be well-explained. In particular, we find that the nitrogen and hydrocarbon chemistry can be significantly affected by temperature, with the possibility of high-temperature pathways to HCN. While we cannot determine the sulphur reservoir, the observations can be explained by models with the majority of the sulphur in CS in the cold gas, SO$_{2}$ in the warm gas, and atomic sulphur in the warmest gas. The observed abundances of ions such as HCO$^+$ and N$_2$H$^+$ and the cold gas-phase production of HCN constrain the cosmic-ray ionization rate to $\sim 5.6 \times 10^{-17}$ s$^{-1}$, to within a factor of three. Finally, we find that the model and observations can simultaneously agree at a reasonable level and often to within a factor of three for $7 \times 10^{3} \leq t(\mathrm{yrs}) \leq 5 \times 10^{4}$, with a strong preference for $t \sim 3 \times 10^{4}$ yrs since the collapse and formation of the central luminosity source.

The Caltech Submillimeter Observatory has detected triply deuterated ammonia, ND3, through its 10a-00s transition near 310 GHz. Emission is found in the NGC 1333 region, both towards IRAS 4A and a position to the South-East where DCO+ peaks. In both cases, the hyperfine ratio indicates that the emission is optically thin. Column densities of ND3 are 3--6 x 10^11 cm^-2 for T_ex=10 K and twice as high for T_ex=5 K. Using a Monte Carlo radiative transfer code and a model of the structure of the IRAS source with temperature and density gradients, the estimated ND3 abundance is 3.2 x 10^-12 if ND3/H2 is constant throughout the envelope. In the more likely case that ND3/H2D+ is constant, ND3/H2 peaks in the cold outer parts of the source at a value of 1.0 x 10^-11. To reproduce the observed NH3/ND3 abundance ratio of ~1000, grain surface chemistry requires an atomic D/H ratio of ~0.15 in the gas phase, >10 times higher than in recent chemical models. More likely, the deuteration of NH3 occurs by ion-molecule reactions in the gas phase, in which case the data indicate that deuteron transfer reactions are much faster than proton transfers.

This paper is a review and an update on recent work on the physical and chemical structure of the envelopes of newly born massive stars, at the stages preceding ultracompact H II regions. It discusses methods and results to determine total mass, temperature and density structure, ionization rate, and depth-dependent chemical composition.

We present observations of the v2=0 and vibrationally excited v2=1 J=9-8 rotational lines of HCN at 797 GHz toward the deeply embedded massive young stellar object GL 2591, which provide the missing link between the extended envelope traced by lower-J line emission and the small region of hot (T_ex >= 300 K), abundant HCN seen in 14 micron absorption with the Infrared Space Observatory (ISO). The line ratio yields T_ex=720^+135_-100 K and the line profiles reveal that the hot gas seen with ISO is at the velocity of the protostar, arguing against a location in the outflow or in shocks. Radiative transfer calculations using a depth-dependent density and temperature structure show that the data rule out a constant abundance throughout the envelope, but that a model with a jump of the abundance in the inner part by two orders of magnitude matches the observations. Such a jump is consistent with the sharp increase in HCN abundance at temperatures >~230 K predicted by recent chemical models in which atomic oxygen is driven into water at these temperatures. Together with the evidence for ice evaporation in this source, this result suggests that we may be witnessing the birth of a hot core. Thus, GL 2591 may represent a rare class of objects at an evolutionary stage just preceding the `hot core' stage of massive star formation.

We present a numerical method and computer code to calculate the radiative transfer and excitation of molecular lines. Formulating the Monte Carlo method from the viewpoint of cells rather than photons allows us to separate local and external contributions to the radiation field. This separation is critical to accurate and fast performance at high optical depths (tau>100). The random nature of the Monte Carlo method serves to verify the independence of the solution to the angular, spatial, and frequency sampling of the radiation field. These features allow use of our method in a wide variety of astrophysical problems without specific adaptations: in any axially symmetric source model and for all atoms or molecules for which collisional rate coefficients are available. Continuum emission and absorption by dust is explicitly taken into account but scattering is neglected. We illustrate these features in calculations of (i) the HCO+ J=1-0 and 3-2 emission from a flattened protostellar envelope with infall and rotation, (ii) the CO, HCO+, CN and HCN emission from a protoplanetary disk and (iii) HCN emission from a high-mass young stellar object, where infrared pumping is important. The program can be used for optical depths up to 1000-10,000, depending on source model. We expect this program to be an important tool in analysing data from present and future infrared and (sub) millimetre telescopes.

The chemical evolution of thirteen regions of massive star formation is studied through single-dish and interferometer line observations at submillimeter wavelengths. The wealth of CH3OH and H2CO lines covering a large range of excitation conditions allows us to calculate radial abundance profiles, using the physical structures of the sources derived earlier from submillimeter continuum and CS lines. The data indicate three types of abundance profiles: flat profiles at CH3OH/H2 ~10^-9 for the coldest sources, profiles with a jump in its abundance from ~10^-9 to ~10^-7 for the warmer sources, and flat profiles at CH3OH/H2 ~few 10^-8 for the hot cores. The location of the jump at T~100 K suggests that it is due to evaporation of grain mantles, followed by destruction in gas-phase reactions in the hot core stage. In contrast, the H2CO data can be well fit with a constant abundance of a few x 10^-9 throughout the envelope, providing limits on its grain surface formation. These results indicate that T_rot(CH3OH) can be used as evolutionary indicator during the embedded phase of massive star formation, independent of source optical depth or orientation.

Recent models of the envelopes of seven massive protostars are used to analyze observations of H3+ infrared absorption and H13CO+ submillimeter emission lines toward these stars, and to constrain the cosmic-ray ionization rate zeta. The H13CO+ gives best-fit values of zeta=(2.6+/- 1.8) x 10^-17 s^-1, in good agreement with diffuse cloud models and with recent Voyager/Pioneer data but factors of up to 7 lower than found from the H3+ data. No relation of zeta with luminosity or total column density is found, so that local (X-ray) ionization and shielding against cosmic rays appear unimportant for these sources. The difference between the H3+ and H13CO+ results and the correlation of N(H3+) with heliocentric distance suggest that intervening clouds contribute significantly to the H3+ absorptions in the more distant regions. The most likely absorbers are low-density (<~10^4 cm^-3) clouds with most carbon in neutral form or in CO.

Recent observational studies of intermediate- and high-mass star-forming regions at submillimeter and infrared wavelengths are reviewed, and chemical diagnostics of the different physical components associated with young stellar objects are summarized. Procedures for determining the temperature, density and abundance profiles in the envelopes are outlined. A detailed study of a set of infrared-bright massive young stars reveals systematic increases in the gas/solid ratios, the abundances of evaporated molecules, and the fraction of heated ices with increasing temperature. Since these diverse phenomena involve a range of temperatures from < 100 K to 1000 K, the enhanced temperatures must be communicated to both the inner and outer parts of the envelopes. This `global heating' plausibly results from the gradual dispersion of the envelopes with time. Similarities and differences with low-mass YSOs are discussed. The availability of accurate physical models will allow chemical models of ice evaporation followed by `hot core' chemistry to be tested in detail.

The physical structure of the envelopes around a sample of fourteen massive (1000-100,000 solar L) young stars is investigated on 100- 100,000 AU scales using maps and spectra in submillimeter continuum and lines of C17O, CS and H2CO. The total column densities and the temperature profiles are obtained by fitting self-consistent dust models to submillimeter photometry. Both the molecular line and dust emission data indicate density gradients ~r^{-alpha}, with alpha=1.0-1.5, significantly flatter than the alpha=2.0 generally found for low-mass objects. This flattening may indicate that in massive young stellar objects, nonthermal pressure is more important for the support against gravitational collapse, while thermal pressure dominates for low-mass sources. We find alpha=2 for two hot core-type sources, but regard this as an upper limit since in these objects, the CS abundance may be enhanced in the warm gas close to the star.

We report the detection of the 1_{10}-1_{11} ground state transition of ortho-H_2D^+ at 372.421 GHz in emission from the young stellar object NGC 1333 IRAS 4A. Detailed excitation models with a power-law temperature and density structure yield a beam-averaged H_2D^+ abundance of 3 x 10^{-12} with an uncertainty of a factor of two. The line was not detected toward W 33A, GL 2591, and NGC 2264 IRS, in the latter source at a level which is 3-8 times lower than previous observations. The H_2D^+ data provide direct evidence in support of low-temperature chemical models in which H_2D^+ is enhanced by the reaction of H_3^+ and HD. The H_2D^+ enhancement toward NGC 1333 IRAS 4A is also reflected in the high DCO^+/HCO^+ abundance ratio. Simultaneous observations of the N_2H^+ 4-3 line show that its abundance is about 50-100 times lower in NGC 1333 IRAS 4A than in the other sources, suggesting significant depletion of N_2. The N_2H^+ data provide independent lower limits on the H_3^+ abundance which are consistent with the abundances derived from H_2D^+. The corresponding limits on the H_3^+$ column density agree with recent near-infrared absorption measurements of H_3^+ toward W 33A and GL 2591.

The temperature, density and kinematics of the gas and dust surrounding the luminous young stellar object GL~2591 are investigated on scales as small as 100 AU, probed by 4.7 micron absorption spectroscopy, to over 60,000 AU, probed by single-dish submillimeter spectroscopy. These two scales are connected by interferometric 86-226 GHz images of size 60,000 AU and resolution 2000 AU in continuum and molecular lines. The data are used to constrain the physical structure of the envelope and investigate the influence of the young star on its immediate surroundings.