{Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness.

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

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)

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

We propose a statistical tool to compare the scaling behaviour of turbulence in pairs of molecular cloud maps. Using artificial maps with well defined spatial properties, we calibrate the method and test its limitations to ultimately apply it to a set of observed maps. We develop the wavelet-based weighted cross-correlation (WWCC) method to study the relative contribution of structures of different sizes and their degree of correlation in two maps as a function of spatial scale, and the mutual displacement of structures in the molecular cloud maps. We test the WWCC for circular structures having a single prominent scale and fractal structures showing a self-similar behavior without prominent scales. Observational noise and a finite map size limit the scales where the cross-correlation coefficients and displacement vectors can be reliably measured. For fractal maps containing many structures on all scales, the limitation from the observational noise is negligible for signal-to-noise ratios >5. (abbrev). Application of the WWCC to the observed line maps of the giant molecular cloud G333 allows to add specific scale information to the results obtained earlier using the principle component analysis. It confirms the chemical and excitation similarity of $^{13}$CO and C$^{18}$O on all scales, but shows a deviation of HCN at scales of up to 7' (~7 pc). This can be interpreted as a chemical transition scale. The largest structures also show a systematic offset along the filament, probably due to a large-scale density gradient. The WWCC can compare correlated structures in different maps of molecular clouds identifying scales that represent structural changes such as chemical and phase transitions and prominent or enhanced dimensions.

We probe the chemical and energetic conditions in dense gas created by radiative feedback through observations of multiple CO, HCN and HCO$^+$ transitions toward the dense core of M17 SW. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain maps of the $J=16-15$, $J=12-11$, and $J=11-10$ transitions of $^{12}$CO. We compare these maps with corresponding APEX and IRAM 30m telescope data for low- and mid-$J$ CO, HCN and HCO$^+$ emission lines, including maps of the HCN $J=8-7$ and HCO$^+$ $J=9-8$ transitions. The excitation conditions of $^{12}$CO, HCO$^+$ and HCN are estimated with a two-phase non-LTE radiative transfer model of the line spectral energy distributions (LSEDs) at four selected positions. The energy balance at these positions is also studied. We obtained extensive LSEDs for the CO, HCN and HCO$^+$ molecules toward M17 SW. The LSED shape, particularly the high-$J$ tail of the CO lines observed with SOFIA/GREAT, is distinctive for the underlying excitation conditions. The critical magnetic field criterion implies that the cold cloudlets at two positions are partially controlled by processes that create and dissipate internal motions. Supersonic but sub-Alfv\'enic velocities in the cold component at most selected positions indicates that internal motions are likely MHD waves. Magnetic pressure dominates thermal pressure in both gas components at all selected positions, assuming random orientation of the magnetic field. The magnetic pressure of a constant magnetic field throughout all the gas phases can support the total internal pressure of the cold components, but it cannot support the internal pressure of the warm components. If the magnetic field scales as $B \propto n^{2/3}$, then the evolution of the cold cloudlets at two selected positions, and the warm cloudlets at all selected positions, will be determined by ambipolar diffusion.

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

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.

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.

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.

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

Column-density maps of molecular clouds are one of the most important observables in the context of molecular cloud- and star-formation (SF) studies. With the Herschel satellite it is now possible to determine the column density from dust emission. We use observations and simulations to demonstrate how LOS contamination affects the column density probability distribution function (PDF). We apply a first-order approximation (removing a constant level) to the molecular clouds of Auriga, Maddalena, Carina and NGC3603. In perfect agreement with the simulations, we find that the PDFs become broader, the peak shifts to lower column densities, and the power-law tail of the PDF flattens after correction. All PDFs have a lognormal part for low column densities with a peak at Av~2, a deviation point (DP) from the lognormal at Av(DP)~4-5, and a power-law tail for higher column densities. Assuming a density distribution rho~r^-alpha, the slopes of the power-law tails correspond to alpha(PDF)=1.8, 1.75, and 2.5 for Auriga, Carina, and NGC3603 (alpha~1.5-2 is consistent gravitational collapse). We find that low-mass and high-mass SF clouds display differences in the overall column density structure. Massive clouds assemble more gas in smaller cloud volumes than low-mass SF ones. However, for both cloud types, the transition of the PDF from lognormal shape into power-law tail is found at the same column density (at Av~4-5 mag). Low-mass and high-mass SF clouds then have the same low column density distribution, most likely dominated by supersonic turbulence. At higher column densities, collapse and external pressure can form the power-law tail. The relative importance of the two processes can vary between clouds and thus lead to the observed differences in PDF and column density structure.

We probe the column densities and masses traced by the ionized and neutral atomic carbon with spectrally resolved maps, and compare them to the diffuse and dense molecular gas traced by [C I] and low-$J$ CO lines toward the star-forming region M17SW. We mapped a 4.1pc x 4.7pc region in the [C I] 609 m$\mu$ line using the APEX telescope, as well as the CO isotopologues with the IRAM 30m telescope. We analyze the data based on velocity channel maps that are 1 km/s wide. We correlate their spatial distribution with that of the [C II] map obtained with SOFIA/GREAT. Optically thin approximations were used to estimate the column densities of [C I] and [C II] in each velocity channel. The spatial distribution of the [C I] and all CO isotopologues emission was found to be associated with that of [C II] in about 20%-80% of the mapped region, with the high correlation found in the central (15-23 km/s ) velocity channels. The excitation temperature of [C I] ranges between 40 K and 100 K in the inner molecular region of M17 SW. Column densities in 1 km/s channels between ~10$^{15}$ and ~10$^{17}$ cm$^{-2}$ were found for [C I]. Just ~20% of the velocity range (~40 km/s) that the [C II] line spans is associated with the star-forming material traced by [C I] and CO. The total gas mass estimated from the [C II] emission gives a lower limit of ~4.4x10$^3$ $M_{\odot}$. At least 64% of this mass is not associated with the star-forming material in M17SW. We also found that about 36%, 17%, and 47% of the [C II] emission is associated with the HII, HI, and H_2 regimes, respectively. Comparisons with the H41$\alpha$ line shows an ionization region mixed with the neutral and part of the molecular gas, in agreement with the clumped structure and dynamical processes at play in M17SW. These results are also relevant to extra-galactic studies in which [C II] is often used as a tracer of star-forming material.

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.

Context: How do molecular clouds form out of the atomic phase? And what are the relative fractions of carbon in the ionized, atomic and molecular phase? These are questions at the heart of cloud and star formation. Methods: Using multiple observatories from Herschel and SOFIA to APEX and the IRAM 30m telescope, we mapped the ionized, atomic and molecular carbon ([CII]@1900GHz, [CI]@492GHz and C18O(2-1)@220GHz) at high spatial resolution (12"-25") in four young massive infrared dark clouds (IRDCs). Results: The three carbon phases were successfully mapped in all four regions, only in one source the [CII] line remained a non-detection. Both the molecular and atomic phases trace the dense structures well, with [CI] also tracing material at lower column densities. [CII] exhibits diverse morphologies in our sample, from compact to diffuse structures probing the cloud environment. In at least two out of the four regions, we find kinematic signatures strongly indicating that the dense gas filaments have formed out of a dynamically active and turbulent atomic/molecular cloud, potentially from converging gas flows. The atomic-to-molecular carbon gas mass ratios are low between 7% and 12% with the lowest values found toward the most quiescent region. In the three regions where [CII] is detected, its mass is always higher by a factor of a few than that of the atomic carbon. The ionized carbon emission depends as well on the radiation field, however, we also find strong [CII] emission in a region without significant external sources, indicating that other processes, e.g., energetic gas flows can contribute to the [CII] excitation as well.

CCAT will be a 25m diameter sub-millimeter telescope capable of operating in the 0.2 to 2.1mm wavelength range. It will be located at an altitude of 5600m on Cerro Chajnantor in northern Chile near the ALMA site. The anticipated first generation instruments include large format (60,000 pixel) kinetic inductance detector (KID) cameras, a large format heterodyne array and a direct detection multi-object spectrometer. The paper describes the architecture of the CCAT software and the development strategy.

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars.

The Heterodyne Instrument for the Far Infrared (HIFI) on board the Herschel Space Observatory is composed of a set of fourteen double sideband mixers. We discuss the general problem of the sideband ratio (SBR) determination and the impact of an imbalanced sideband ratio on the line calibration in double sideband heterodyne receivers. The HIFI SBR is determined from a combination of data taken during pre-launch gas cell tests and in-flight. The results and some of the calibration artefacts discovered in the gas cell test data are presented here along with some examples of how these effects appear in science data taken in orbit.

We study how the estimation of the sonic Mach number ($M_s$) from $^{13}$CO linewidths relates to the actual 3D sonic Mach number. For this purpose we analyze MHD simulations which include post-processing to take radiative transfer effects into account. As expected, we find very good agreement between the linewidth estimated sonic Mach number and the actual sonic Mach number of the simulations for optically thin tracers. However, we find that opacity broadening causes $M_s$ to be overestimated by a factor of ~ 1.16 - 1.3 when calculated from optically thick $^{13}$CO lines. We also find that there is a dependency on the magnetic field: super-Alfv\'enic turbulence shows increased line broadening as compared with sub-Alfv\'enic turbulence for all values of optical depth for supersonic turbulence. Our results have implications for the observationally derived sonic Mach number--density standard deviation ($\sigma_{\rho/<\rho>}$) relationship, $\sigma^2_{\rho/<\rho>}=b^2M_s^2$, and the related column density standard deviation ($\sigma_{N/<N>}$) sonic Mach number relationship. In particular, we find that the parameter b, as an indicator of solenoidal vs. compressive driving, will be underestimated as a result of opacity broadening. We compare the $\sigma_{N/<N>}-M_s$ relation derived from synthetic dust extinction maps and $^{13}$CO linewidths with recent observational studies and find that solenoidally driven MHD turbulence simulations have values of $\sigma_{N/<N>}$ which are lower than real molecular clouds. This may be due to the influence of self-gravity which should be included in simulations of molecular cloud dynamics.

CCAT will be a 25-meter telescope for sub millimeter astronomy located at 5600m altitude on Cerro Chajnantor in northern Chile. CCAT will combine high sensitivity, a wide field of view, and a broad wavelength range (0.35 to 2.1mm) to provide an unprecedented capability for deep, large-area multicolor submillimeter surveys. It is planned to have a suite of instruments including large format KID cameras, a large heterodyne array and a KID-based direct detection multi-object spectrometer. The remote location drives a desire for fully autonomous observing coupled with data reduction pipelines and fast feedback to principal investigators.

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

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. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics.

We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas.

The Fourier power spectrum is one of the most widely used statistical tools to analyze the nature of magnetohydrodynamic turbulence in the interstellar medium. Lazarian & Pogosyan (2004) predicted that the spectral slope should saturate to -3 for an optically thick medium and many observations exist in support of their prediction. However, there have not been any numerical studies to-date testing these results. We analyze the spatial power spectrum of MHD simulations with a wide range of sonic and Alfv\'enic Mach numbers, which include radiative transfer effects of the $^{13}$CO transition. We confirm numerically the predictions of Lazarian & Pogosyan (2004) that the spectral slope of line intensity maps of an optically thick medium saturates to -3. Furthermore, for very optically thin supersonic CO gas, where the density or CO abundance values are too low to excite emission in all but the densest shock compressed gas, we find that the spectral slope is shallower than expected from the column density. Finally, we find that mixed optically thin/thick CO gas, which has average optical depths on order of unity, shows mixed behavior: for super-Alfv\'enic turbulence, the integrated intensity power spectral slopes generally follow the same trend with sonic Mach number as the true column density power spectrum slopes. However, for sub-Alfv\'enic turbulence the spectral slopes are steeper with values near -3 which are similar to the very optically thick regime.

We study the effects of radiative transfer on the Probability Distribution Functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied $^{13}$CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have $\tau \approx 1$ best matching the PDF of the column density. We fit a 2D variance-sonic Mach number relationship to our logarithmic PDFs of the form $\sigma_{ln(\Sigma/\Sigma_0)}^2=A\times ln(1+b^2{\cal M}_s^2)$ and find that, for parameter $b=1/3$, parameter $A$ depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the $\tau \approx 1$ case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always sub-thermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed and consequently, the observable PDFs, are less determined by turbulence and more-often determined by radiative transfer effects.

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$\mu$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

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)

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 revised the treatment of interstellar dust in the KOSMA-\tau\ PDR model code to achieve a consistent description of the dust-related physics in the code. The detailed knowledge of the dust properties is then used to compute the dust continuum emission together with the line emission of chemical species. We coupled the KOSMA-\tau\ PDR code with the MCDRT (multi component dust radiative transfer) code to solve the frequency-dependent radiative transfer equations and the thermal balance equation in a dusty clump under the assumption of spherical symmetry, assuming thermal equilibrium in calculating the dust temperatures, neglecting non-equilibrium effects. We updated the calculation of the photoelectric heating and extended the parametrization range for the photoelectric heating toward high densities and UV fields. We revised the computation of the H2 formation on grain surfaces to include the Eley-Rideal effect, thus allowing for high-temperature H2 formation. We demonstrate how the different optical properties, temperatures, and heating and cooling capabilities of the grains influence the physical and chemical structure of a model cloud. The most influential modification is the treatment of H2 formation on grain surfaces that allows for chemisorption. This increases the total H2 formation significantly and the connected H2 formation heating provides a profound heating contribution in the outer layers of the model clumps. The contribution of PAH surfaces to the photoelectric heating and H2 formation provides a boost to the temperature of outer cloud layers, which is clearly traced by high-J CO lines. Increasing the fraction of small grains in the dust size distribution results in hotter gas in the outer cloud layers caused by more efficient heating and cooler cloud centers, which is in turn caused by the more efficient FUV extinction.

We upgraded the chemical network from the UMIST Database for Astrochemistry 2006 to include isotopes such as ^{13}C and ^{18}O. This includes all corresponding isotopologues, their chemical reactions and the properly scaled reaction rate coefficients. We study the fractionation behavior of astrochemically relevant species over a wide range of model parameters, relevant for modelling of photo-dissociation regions (PDRs). We separately analyze the fractionation of the local abundances, fractionation of the total column densities, and fractionation visible in the emission line ratios. We find that strong C^+ fractionation is possible in cool C^+ gas. Optical thickness as well as excitation effects produce intensity ratios between 40 and 400. The fractionation of CO in PDRs is significantly different from the diffuse interstellar medium. PDR model results never show a fractionation ratio of the CO column density larger than the elemental ratio. Isotope-selective photo-dissociation is always dominated by the isotope-selective chemistry in dense PDR gas. The fractionation of C, CH, CH^+, and HCO^+ is studied in detail, showing that the fractionation of C, CH and CH^+ is dominated by the fractionation of their parental species. The light hydrides chemically derive from C^+, and, consequently, their fractionation state is coupled to that of C^+. The fractionation of C is a mixed case depending on whether formation from CO or HCO^+ dominates. Ratios of the emission lines of [C II], [C I], ^{13}CO, and H^{13}CO^+ provide individual diagnostics to the fractionation status of C^+, C, and CO.

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

The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data.

We present a cluster analysis of the bright main-sequence and faint pre--main-sequence stellar populations of a field ~ 90 x 90 pc centered on the HII region NGC 346/N66 in the Small Magellanic Cloud, from imaging with HST/ACS. We extend our earlier analysis on the stellar cluster population in the region to characterize the structuring behavior of young stars in the region as a whole with the use of stellar density maps interpreted through techniques designed for the study of the ISM structuring. In particular, we demonstrate with Cartwrigth & Whitworth's Q parameter, dendrograms, and the Delta-variance wavelet transform technique that the young stellar populations in the region NGC 346/N66 are hierarchically clustered, in agreement with other regions in the Magellanic Clouds observed with HST. The origin of this hierarchy is currently under investigation.

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud.

With new THz maps that cover an area of ~3.3x2.1 pc^2 we probe the spatial distribution and association of the ionized, neutral and molecular gas components in the M17 SW nebula. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain a 5'.7x3'.7 map of the 12CO J=13-12 transition and the [C II] 158 um fine-structure line in M17 SW and compare the spectroscopically resolved maps with corresponding ground-based data for low- and mid-J CO and [C I] emission. For the first time SOFIA/GREAT allow us to compare velocity-resolved [C II] emission maps with molecular tracers. We see a large part of the [C II] emission, both spatially and in velocity, that is completely non-associated with the other tracers of photon-dominated regions (PDR). Only particular narrow channel maps of the velocity-resolved [C II] spectra show a correlation between the different gas components, which is not seen at all in the integrated intensity maps. These show different morphology in all lines but give hardly any information on the origin of the emission. The [C II] 158 um emission extends for more than 2 pc into the M17 SW molecular cloud and its line profile covers a broader velocity range than the 12CO J=13-12 and [C I] emissions, which we interpret as several clumps and layers of ionized carbon gas within the telescope beam. The high-J CO emission emerges from a dense region between the ionized and neutral carbon emissions, indicating the presence of high-density clumps that allow the fast formation of hot CO in the irradiated complex structure of M17 SW. The [C II] observations in the southern PDR cannot be explained with stratified nor clumpy PDR models.

By observing radiation-affected gas in the Cepheus B molecular cloud we probe whether the sequential star formation in this source is triggered by the radiation from newly formed stars. We used the dual band receiver GREAT onboard SOFIA to map [C II] and CO 13--12 and 11--10 in Cep B and compared the spatial distribution and the spectral profiles with complementary ground-based data of low-$J$ transitions of CO isotopes, atomic carbon, and the radio continuum. The interaction of the radiation from the neighboring OB association creates a large photon-dominated region (PDR) at the surface of the molecular cloud traced through the photoevaporation of C^+. Bright internal PDRs of hot gas are created around the embedded young stars, where we detect evidence of the compression of material and local velocity changes; however, on the global scale we find no indications that the dense molecular material is dynamically affected.

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

We aim to better understand how the spatial structure of molecular clouds is governed by turbulence. For that, we study the large-scale spatial distribution of low density molecular gas and search for characteristic length scales. We employ a 35 square degrees 13CO 1-0 molecular line survey of Cygnus X and visual extinction (A_V) maps of 17 Galactic clouds to analyse the spatial structure using the Delta-variance method. This sample contains a large variety of different molecular cloud types with different star forming activity. The Delta-variance spectra obtained from the A_V maps show differences between low-mass star-forming (SF) clouds and massive giant molecular clouds (GMC) in terms of shape of the spectrum and its power-law exponent beta. Low-mass SF clouds have a double-peak structure with characteristic size scales around 1 pc (though with a large scatter around this value) and 4 pc. GMCs show no characteristic scale in the A_V-maps, which can partly be ascribed to a distance effect due to a larger line-of-sight (LOS) confusion. The Delta-variance for Cygnus, determined from the 13CO survey, shows characteristic scales at 4 pc and 40 pc, either reflecting the filament structure and large-scale turbulence forcing or - for the 4 pc scale - the scale below which the 13CO 1-0 line becomes optically thick. Though there are different processes that can introduce characteristic scales, i.e. geometry, decaying turbulence the transition scale from supersonic to subsonic turbulence (the sonic scale), line-of-sight effects and energy injection due to expanding supernova shells, outflows, HII-regions, and the relative contribution of these effects strongly varies from cloud to cloud, it is remarkable that the resulting turbulent structure of molecular clouds shows similar characteristics.

Aims: We aim at deriving the excitation conditions of the interstellar gas as well as the local FUV intensities in the molecular cloud surrounding NGC 3603 to get a coherent picture of how the gas is energized by the central stars. Methods: The NANTEN2-4m submillimeter antenna is used to map the [CI] 1-0, 2-1 and CO 4-3, 7-6 lines in a 2' x 2' region around the young OB cluster NGC 3603 YC. These data are combined with C18O 2-1 data, HIRES-processed IRAS 60 and 100 micron maps of the FIR continuum, and Spitzer/IRAC maps. Results: The NANTEN2 observations show the presence of two molecular clumps located south-east and south-west of the cluster and confirm the overall structure already found by previous CS and C18O observations. We find a slight position offset of the peak intensity of CO and [CI], and the atomic carbon appears to be further extended compared to the molecular material. We used the HIRES far-infrared dust data to derive a map of the FUV field heating the dust. We constrain the FUV field to values of \chi = 3 - 6 \times 10^3 in units of the Draine field across the clouds. Approximately 0.2 to 0.3 % of the total FUV energy is re-emitted in the [CII] 158 {\mu}m cooling line observed by ISO. Applying LTE and escape probability calculations, we derive temperatures (TMM1 = 43 K, TMM2 = 47 K), column densities (N(MM1) = 0.9 \times 10^22 cm^-2, N(MM2) = 2.5 \times 10^22 cm^-2) and densities (n(MM1) = 3 \times 10^3 cm^-3, n(MM2) = 10^3 -10^4 cm^-3) for the two observed molecular clumps MM1 and MM2. Conclusions: The cluster is strongly interacting with the ambient molecular cloud, governing its structure and physical conditions. A stability analysis shows the existence of gravitationally collapsing gas clumps which should lead to star formation. Embedded IR sources have already been observed in the outskirts of the molecular cloud and seem to support our conclusions.

We report the first spectroscopic mapping of an atomic carbon line in an infrared dark cloud (IRDC). By observing the spatial distribution of the [CI] emission in an IRDC, comparing it with the 13CO emission and the known distribution of internal heating sources, we can quantify the role of internal and external UV irradiation in the production of atomic carbon. We used the 2x4 pixel SMART receiver of the KOSMA observatory on Gornergrat to map the [CI] 3P1 - 3P0 line in the IRDC G48.65-0.29 and compared the resulting spectra with data from the BU-FCRAO 13CO 1-0 Galactic Ring Survey. The [CI]/13CO effective beam temperature ratio falls at about 0.3 with local deviations by less than a factor two. All velocity components seen in 13CO are also detected in[CI]. We find, however, significant differences in the morphology of the brightest regions seen in the two tracers. While 13CO basically follows the column density distribution derived from the near-infrared (NIR) extinction and the submm continuum, the [CI] emission peaks at the locations of the two known NIR point sources. We find CI/CO abundance ratios between 0.07 and 0.13, matching the lower end of the range previously measured in star-forming regions. Evaluating the relative importance of the irradiation by embedded sources and by the Galactic interstellar UV field, we find that in G48.65-0.29 most [CI] emission can be attributed to externally illuminated surfaces. Embedded sources have a significant impact on the overall abundance distribution of atomic carbon as soon as they are found in an evolved state with noticeable NIR flux.

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.

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.

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.

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

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

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.

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

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

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

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

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.

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

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

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

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

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

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

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

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

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

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.

We present multi-molecular line maps obtained with the Mopra Telescope towards the southern giant molecular cloud (GMC) complex G333, associated with the HII region RCW 106. We have characterised the GMC by decomposing the 3D data cubes with GAUSSCLUMPS, and investigated spatial correlations among different molecules with principal component analysis (PCA). We find no correlation between clump size and line width, but a strong correlation between emission luminosity and line width. PCA classifies molecules into high and low density tracers, and reveals that HCO+ and N2H+ are anti-correlated.

Astronomic line mapping with single-pixel instruments is usually performed in an on-the-fly (OTF) or a raster-mapping mode depending on the capabilities of the telescope and the instrument. The observing efficiency can be increased by combining several source-point integrations with a common reference measurement. This is implemented at many telescopes, but a thorough investigation of the optimum calibration of the modes and the best way of performing these observations is still lacking. We use knowledge of the instrumental stability obtained by an Allan variance measurement to derive a mathematical formalism for optimizing the setup of mapping observations. Special attention has to be paid to minimizing of the impact of correlated noise introduced by the common OFF integrations and to the correction of instrumental drifts. Both aspects can be covered using a calibration scheme that interpolates between two OFF measurements and an appropriate OFF integration time. The total uncertainty of the calibrated data consisting of radiometric noise and drift noise can be minimized by adjusting the source integration time and the number of data points observed between two OFF measurements. It turns out that OTF observations are very robust. They provide a low relative noise, even if their setup deviates considerably from the optimum. Fast data readouts are often essential to minimize the drift contributions. In particular, continuum measurements may be easily spoiled by instrumental drifts. The main drawback of the described mapping modes is the limited use of the measured data at different spatial or spectroscopic resolutions obtained by additional rebinning.

In this paper we discuss the physical conditions of clumpy nature in the IC 348 molecular cloud. We combine new observations of fully sampled maps in [C I] at 492 GHz and 12CO 4--3, taken with the KOSMA 3 m telescope at about 1' resolution, with FCRAO data of 12CO 1--0, 13CO 1--0 and far-infrared continuum data observed by HIRES/IRAS. To derive the physical parameters of the region we analyze the three different line ratios. A first rough estimate of abundance is obtained from an LTE analysis. To understand the [C I] and CO emission from the PDRs in IC 348, we use a clumpy PDR model. With an ensemble of identical clumps, we constrain the total mass from the observed absolute intensities. Then we apply a more realistic clump distribution model with a power law index of 1.8 for clump-mass spectrum and a power law index of 2.3 for mass-size relation. We provide detailed fits to observations at seven representative positions in the cloud, revealing clump densities between 4 10$^{4}$ cm$^{-3}$ and 4 10$^{5}$ cm$^{-3}$ and C/CO column density ratios between 0.02 and 0.26. The derived FUV flux from the model fit is consistent with the field calculated from FIR continuum data, varying between 2 and 100 Draine units across the cloud. We find that both an ensemble of identical clumps and an ensemble with a power law clump mass distribution produce line intensities which are in good agreement (within a factor ~ 2) with the observed intensities. The models confirm the anti-correlation between the C/CO abundance ratio and the hydrogen column density found in many regions.

The fractal structure of the interstellar medium suggests that the interaction of UV radiation with the ISM as described in the context of photon-dominated regions (PDR) dominates most of the physical and chemical conditions, and hence the far-infrared and submm emission from the ISM in the Milky Way. We investigate to what extent the Galactic FIR line emission of the important species CO, C, C+, and O, as observed by the Cosmic Background Explorer (COBE) satellite can be modeled in the framework of a clumpy, UV-penetrated cloud scenario. The far-infrared line emission of the Milky Way is modeled as the emission from an ensemble of clumps with a power law clump mass spectrum and mass-size relation with power-law indices consistent with the observed ISM structure. The individual clump line intensities are calculated using the KOSMA-tau PDR-model for spherical clumps. The model parameters for the cylindrically symmetric Galactic distribution of the mass density and volume filling factor are determined by the observed radial distributions. A constant FUV intensity, in which the clumps are embedded, is assumed. We show that this scenario can explain, without any further assumptions and within a factor of about 2, the absolute FIR-line intensities and their distribution with Galactic longitude as observed by COBE.

The Delta-variance analysis is an efficient tool for measuring the structural scaling behaviour of interstellar turbulence in astronomical maps. In paper I we proposed essential improvements to the Delta-variance analysis. In this paper we apply the improved Delta-variance analysis to i) a hydrodynamic turbulence simulation with prominent density and velocity structures, ii) an observed intensity map of rho Oph with irregular boundaries and variable uncertainties of the different data points, and iii) a map of the turbulent velocity structure in the Polaris Flare affected by the intensity dependence on the centroid velocity determination. The tests confirm the extended capabilities of the improved Delta-variance analysis. Prominent spatial scales were accurately identified and artifacts from a variable reliability of the data were removed. The analysis of the hydrodynamic simulations showed that the injection of a turbulent velocity structure creates the most prominent density structures are produced on a scale somewhat below the injection scale. The new analysis of a rho Oph continuum map reveals an intermediate stage in the molecular cloud evolution showing both signatures of the typical molecular cloud scaling behaviour and the formation of condensed cores. When analysing the velocity structure of the Polaris Flare we show that a universal power law connects scales from 0.03 pc to 3 pc. However, a plateau in the Delta-variance spectrum around 5 pc indicates that the visible large-scale velocity gradient is not converted directly into a turbulent cascade.

The Delta-variance analysis, has proven to be an efficient and accurate method of characterising the power spectrum of interstellar turbulence. The implementation presently in use, however, has several shortcomings. We propose and test an improved Delta-variance algorithm for two-dimensional data sets, which is applicable to maps with variable error bars and which can be quickly computed in Fourier space. We calibrate the spatial resolution of the Delta-variance spectra. The new Delta-variance algorithm is based on an appropriate filtering of the data in Fourier space. It allows us to distinguish the influence of variable noise from the actual small-scale structure in the maps and it helps for dealing with the boundary problem in non-periodic and/or irregularly bounded maps. We try several wavelets and test their spatial sensitivity using artificial maps with well known structure sizes. It turns out that different wavelets show different strengths with respect to detecting characteristic structures and spectral indices, i.e. different aspects of map structures. As a reasonable universal compromise for the optimum Delta-variance filter, we propose the Mexican-hat filter with a ratio between the diameters of the core and the annulus of 1.5.

Star formation at earlier cosmological times takes place in an interstellar medium with low metallicity. The Large Magellanic Cloud (LMC) is ideally suited to study star formation in such an environment. The physical and chemical state of the ISM in a star forming environment can be constrained by observations of submm and FIR spectral lines of the main carbon carrying species, CO, CI and CII, which originate in the surface layers of molecular clouds illuminated by the UV radiation of the newly formed, young stars. We present high-angular resolution sub-millimeter observations in the N159W region in the LMC obtained with the NANTEN2 telescope of the 12CO J = 4-3, J = 7-6, and 13CO J = 4-3 rotational and [CI] 3P1-3P0 and 3P2-3P1 fine-structure transitions. The 13CO J =4-3 and [CI] 3P2-3P1 transitions are detected for the first time in the LMC. We derive the physical and chemical properties of the low-metallicity molecular gas using an escape probability code and a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. The separate excitation analysis of the submm CO lines and the carbon fine structure lines shows that the emitting gas in the N159W region has temperatures of about 80 K and densities of about 10^4 cm^-3. The estimated C to CO abundance ratio close to unity is substantially higher than in dense massive star-forming regions in the Milky Way. The analysis of all observed lines together, including the [CII] line intensity reported in the literature, in the context of a clumpy cloud PDR model constrains the UV intensity to about \chi ~220 and an average density of the clump ensemble of about 10^5 cm^-3, thus confirming the presence of high density material in the LMC N159W region.

The Allan variance is a standard technique to characterise the stability of spectroscopic instruments used in astronomical observations. The period for switching between source and reference measurement is often derived from the Allan minimum time. We propose a new approach for the computation of the Allan variance of spectrometer data combining the advantages of the two existing methods into a unified scheme. Using the Allan variance spectrum we derive the optimum strategy for symmetric observing schemes minimising the total uncertainty of the data resulting from radiometric and drift noise. The unified Allan variance computation scheme is designed to trace total-power and spectroscopic fluctuations within the same framework. The method includes an explicit error estimate both for the individual Allan variance spectra and for the derived stability time. A new definition of the instrument stability time allows to characterise the instrument even in the case of a fluctuation spectrum shallower than 1/f, as measured for the total power fluctuations in high-electron-mobility transistors. We find a non-linear impact of the binning of spectrometer channels on the resulting noise and the Allan time deviating from the description in existing theoretical treatments.

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

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.

We present an overview of a high-mass star formation region through the major (sub-)mm, and far-infrared cooling lines to gain insight into the physical conditions and the energy budget of the molecular cloud. We used the KOSMA 3m telescope to map the core ($10'\times 14'$) of the Galactic star forming region DR 21/DR 21 (OH) in the Cygnus X region in the two fine structure lines of atomic carbon CI and four mid-$J$ transitions of CO and $^{13}$CO, and CS $J=7\TO6$. These observations have been combined with FCRAO $J=1\TO0$ observations of $^{13}$CO and C$^{18}$O. Five positions, including DR21, DR21 (OH), and DR21 FIR1, were observed with the ISO/LWS grating spectrometer in the \OI 63 and 145 $\mu$m lines, the \CII 158 $\mu$m line, and four high-$J$ CO lines. We discuss the intensities and line ratios at these positions and apply Local Thermal Equilibrium (LTE) and non-LTE analysis methods in order to derive physical parameters such as masses, densities and temperatures. The CO line emission has been modeled up to J=20. From non-LTE modeling of the low- to high-$J$ CO lines we identify two gas components, a cold one at temperatures of T$_\RM{kin}\sim 30-40$ K, and one with T$_\RM{kin}\sim 80-150$ K at a local clump density of about n(H$_2$)$\sim 10^4-10^6$ cm$^{-3}$. While the cold quiescent component is massive containing typically more than 94 % of the mass, the warm, dense, and turbulent gas is dominated by mid- and high-$J$ CO line emission and its large line widths. The medium must be clumpy with a volume-filling of a few percent. The CO lines are found to be important for the cooling of the cold molecular gas, e.g. at DR21 (OH). Near the outflow of the UV-heated source DR21, the gas cooling is dominated by line emission of atomic oxygen and of CO.

The investigation of the statistical properties of maps of line centroids has been used for almost 50 years, but there is still no general agreement on their interpretation. We try to quantify which properties of underlying turbulent velocity fields can be derived from centroid velocity maps, and we test conditions under which the scaling behaviour of the centroid velocities matches the scaling of the three-dimensional velocity field. Using fractal cloud models we study systematically the relation between three-dimensional density and velocity fields and the statistical properties of the produced line centroid maps. We put special attention to cases with large density fluctuations resembling supersonic interstellar turbulence. Starting from the Delta-variance analysis we derive a new tool to compute the scaling behaviour of the three-dimensional velocity field from observed intensity and centroid velocity maps. We provide two criteria to decide whether the information from the centroid velocities directly reflects the properties of the underlying velocity field. Applying these criteria allows to understand the different results found so far in the literature on the interpretation of the statistics of velocity centroids. A new iteration scheme can be used to derive the three-dimensional velocity scaling from centroid velocity maps for arbitrary density and velocity fields, but it requires an accurate knowledge of the average density of the considered interstellar cloud.

We study the effects of a metallicity variation on the thermal balance and [CII] fine-structure line strengths in interstellar photon dominated regions (PDRs). We find that a reduction in the dust-to-gas ratio and the abundance of heavy elements in the gas phase changes the heat balance of the gas in PDRs. The surface temperature of PDRs decreases as the metallicity decreases except for high density ($n>10^6$ cm$^{-3}$) clouds exposed to weak ($\chi< 100$) FUV fields where vibrational H$_2$-deexcitation heating dominates over photoelectric heating of the gas. We incorporate the metallicity dependence in our KOSMA-$\tau$ PDR model to study the metallicity dependence of [CII]/CO line ratios in low metallicity galaxies. We find that the main trend in the variation of the observed CII/CO ratio with metallicity is well reproduced by a single spherical clump, and does not necessarily require an ensemble of clumps as in the semi-analytical model presented by Bolatto et al. (1999).

Characterizing the spatial and velocity structure of molecular clouds is a first step towards a better understanding of interstellar turbulence and its link to star formation. We present observations and structure analysis results for a large-scale (~ 7.10 deg^2) 13CO J = 2--1 and 12 CO J = 3--2 survey towards the nearby Perseus molecular cloud observed with the KOSMA 3m telescope. We study the spatial structure of line-integrated and velocity channel maps, measuring the $\Delta$-variance as a function of size scale. We determine the spectral index $\beta$ of the corresponding power spectrum and study its variation across the cloud and across the lines. We find that the spectra of all CO line-integrated maps of the whole complex show the same index, $\beta$ ~ 3.1, for scales between about 0.2 and 3pc, independent of isotopomer and rotational transition. A complementary 2MASS map of optical extinction shows a noticeably smaller index of 2.6. In contrast to the overall region, the CO maps of individual subregions show a significant variation of $\beta$. The 12CO 3--2 data provide e.g. a spread of indices between 2.9 in L1455 and 3.5 in NGC1333. In general, active star forming regions show a larger power-law exponent. We find that the $\Delta$-variance spectra of individual velocity channel maps are very sensitive to optical depth effects clearly indicating self-absorption in the densest regions. When studying the dependence of the channel-map spectra as a function of the velocity channel width, the expected systematic increase of the spectral index with channel width is only detected in the blue line wings.This could be explained by a filamentary, pillar-like structure which is left at low velocities while the overall molecular gas is swept up by a supernova shock wave.

This paper presents results from modelling 450 micron and 850 micron continuum and HCO+ line observations of three distinct cores of an infrared dark cloud (IRDC) directed toward the W51 GMC. In the sub-mm continuum these cores appear as bright, isolated emission features. One of them coincides with the peak of 8.3 micron extinction as measured by the Midcourse Space Experiment satellite. Detailed radiative transfer codes are applied to constrain the cores' physical conditions to address the key question: Do these IRDC-cores harbour luminous sources? The results of the continuum model, expressed in the $\chi^2$ quality-of-fit parameter, are also constrained by the absence of 100 micron emission from IRAS. For the sub-mm emission peaks this shows that sources of 300 solar luminosities are embedded within the cores. For the extinction peak, the combination of continuum and HCO+ line modelling indicates that a heating source is present as well. Furthermore, the line model provides constraints on the clumpiness of the medium. All three cores have similar masses of about 70-150 solar masses and similar density structures. The extinction peak differs from the other two cores by hosting a much weaker heating source, and the sub-mm emission core at the edge of the IRDC deviates from the other cores by a higher internal clumpiness.

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 compare velocity structure in the Polaris Flare molecular cloud at scales ranging from 0.015 pc to 20 pc to simulations of supersonic hydrodynamic and MHD turbulence computed with the ZEUS MHD code. We use several different statistical methods to compare models and observations. The Delta-variance wavelet transform is most sensitive to characteristic scales and scaling laws, but is limited by a lack of intensity weighting. The scanning-beam size-linewidth relation is more robust with respect to noisy data. Obtaining the global velocity scaling behaviour requires that large-scale trends in the maps not be removed but treated as part of the turbulent cascade. We compare the true velocity PDF in our models to velocity centroids and average line profiles in optically thin lines, and find that the line profiles reflect the true PDF better unless the map size is comparable to the total line-of-sight thickness of the cloud. Comparison of line profiles to velocity centroid PDFs can thus be used to measure the line-of-sight depth of a cloud. The observed density and velocity structure is consistent with supersonic turbulence with a driving scale at or above the size of the molecular cloud and dissipative processes below 0.05 pc. Ambipolar diffusion could explain the dissipation. The velocity PDFs exclude small-scale driving such as that from stellar outflows as a dominant process in the observed region. In the models, large-scale driving is the only process that produces deviations from a Gaussian PDF shape consistent with observations. Strong magnetic fields impose a clear anisotropy on the velocity field, reducing the velocity variance in directions perpendicular to the field. (abridged)

To study the interaction of star-formation and turbulent molecular cloud structuring, we analyse numerical models and observations of self-gravitating clouds using the Delta-variance as statistical measure for structural characteristics. In the models we resolve the transition from purely hydrodynamic turbulence to gravitational collapse associated with the formation and mass growth of protostellar cores. We compare models of driven and freely decaying turbulence with and without magnetic fields. Self-gravitating supersonic turbulence always produces a density structure that contains most power on the smallest scales provided by collapsed cores as soon as local collapse sets in. This is in contrast to non-self-gravitating hydrodynamic turbulence where the Delta-variance is dominated by large scale structures. To detect this effect in star-forming regions observations have to resolve the high density contrast of protostellar cores with respect to their ambient molecular cloud. Using the 3mm continuum map of a star-forming cluster in Serpens we show that the dust emission traces the full density evolution. On the contrary, the density range accessible by molecular line observations is insufficient for this analysis. Only dust emission and dust extinction observations are able to to determine the structural parameters of star-forming clouds following the density evolution during the gravitational collapse.

We present observations of 15 massive cores in three different CS transitions from the FCRAO 14m and the KOSMA 3m telescope. We derive physical parameters of these cores using different approaches to the line radiative transfer problem. The local radiative transfer approximations fail to provide reliable values except for the column densities. A self-consistent explanation of the observed line profiles is only possible when taking density gradients and an internal turbulent structure of the cores into account. The observational data can be fitted by a spherically symmetric radiative transfer model including such gradients and a turbulent clumping. We find that the observed massive cores are approximately virialised with a clumpy density profile that decays with a radial exponent of about -1.6 down to a relatively sharp outer boundary. We show that a careful analysis of spatially unresolved multi-line observations using a physical radiative transfer model can provide values for physical parameters that could be obtained otherwise only by direct observations with much higher spatial resolution. This applies to all quantities directly affecting the line excitation, like the mass and size of dense cores. Information on the exact location or number of clumps, of course, always has to rely on high-resolution observations e.g. from interferometers.

Modeling the structure of molecular clouds depends on good methods to statistically compare simulations with observations in order to constrain the models. Here we characterize a suite of hydrodynamical and magnetohydrodynamical (MHD) simulations of supersonic turbulence using an averaged wavelet transform, the Delta-variance, that has been successfully used to characterize observations. We find that, independent of numerical resolution and dissipation, the only models that produce scale-free, power-law Delta-variance spectra are those with hypersonic Mach numbers above M ~ 4, while slower supersonic turbulence tends to show characteristic scales and produce non-power-law spectra. Magnetic fields have only a minor influence on this tendency, though they tend to reduce the scale-free nature of the turbulence, and increase the transfer of energy from large to small scales. The evolution of the characteristic length scale seen in supersonic turbulence follows exactly the t^{1/2} power-law predicted from recent studies of the kinetic energy decay rate.