The Orion Molecular Cloud is the nearest massive-star forming region. Massive stars have profound effects on their environment due to their strong radiation fields and stellar winds. Velocity-resolved observations of the [CII] $158\,\mu\mathrm{m}$ fine-structure line allow us to study the kinematics of UV-illuminated gas. Here, we present a square-degree-sized map of [CII] emission from the Orion Nebula complex obtained by the upGREAT instrument onboard SOFIA, covering the entire Orion Nebula (M42) plus M43 and the nebulae NGC 1973, 1975, and 1977. We compare the stellar characteristics of these three regions with the kinematics of the expanding bubbles surrounding them. The bubble blown by the O7V star $\theta^1$ Ori C in the Orion Nebula expands rapidly, at $13\,\mathrm{km\,s^{-1}}$. Simple analytical models reproduce the characteristics of the hot interior gas and the neutral shell of this wind-blown bubble and give us an estimate of the expansion time of $0.2\,\mathrm{Myr}$. M43 with the B0.5V star NU Ori also exhibits an expanding bubble structure, with an expansion velocity of $6\,\mathrm{km\,s^{-1}}$. Comparison with analytical models for the pressure-driven expansion of H\,{\sc ii} regions gives an age estimate of $0.02\,\mathrm{Myr}$. The bubble surrounding NGC 1973, 1975, and 1977 with the central B1V star 42 Orionis expands at $1.5\,\mathrm{km\,s^{-1}}$, likely due to the over-pressurized ionized gas as in the case of M43. We derive an age of $0.4\,\mathrm{Myr}$ for this structure. We conclude that the bubble of the Orion Nebula is driven by the mechanical energy input by the strong stellar wind from $\theta^1$ Ori C, while the bubbles associated with M43 and NGC 1977 are caused by the thermal expansion of the gas ionized by their central later-type massive stars.

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback.

Context. The [CII] 158 {\mu}m far-infrared (FIR) fine-structure line is one of the most important cooling lines of the star-forming interstellar medium (ISM). High spectral resolution observations have shown complex structures in the line profiles of the [CII] emission. Aims. Our aim is to determine whether the complex profiles observed in [^{12}CII] are due to individual velocity components along the line-of-sight or to self-absorption based on a comparison of the [^{12}CII] and isotopic [^{13}CII] line profiles. Methods. Deep integrations with the SOFIA/upGREAT 7-pixel array receiver in M43, Horsehead~PDR, Monoceros~R2, and M17~SW allow for the detection of optically thin [^{13}CII] emission lines, along with the [^{12}CII] emission lines, with a high signal-to-noise ratio. We first derived the [^{12}CII] optical depth and the [CII] column density from a single component model. However, the complex line profiles observed require a double layer model with an emitting background and an absorbing foreground. A multi-component velocity fit allows us to derive the physical conditions of the [CII] gas: column density and excitation temperature. Results. We find moderate to high [^{12}CII] optical depths in all four sources and self-absorption of [^{12}CII] in Mon R2 and M17 SW. The high column density of the warm background emission corresponds to an equivalent Av of up to 41 mag. The foreground absorption requires substantial column densities of cold and dense [CII] gas, with an equivalent Av ranging up to about 13 mag. Conclusions. The column density of the warm background material requires multiple photon-dominated region (PDR) surfaces stacked along the line of sight and in velocity. The substantial column density of dense and cold foreground [CII] gas detected in absorption cannot be explained with any known scenario and we can only speculate on its origins

We report the first detection of the isotopologues 13CCC and C13CC. We used the heterodyne receivers GREAT and upGREAT on board SOFIA to search for the ro-vibrational transitions Q(2) and Q(4) of 13CCC and C13CC at 1.9 THz along the line of sight towards SgrB2(M). For both species the ro-vibrational absorption lines Q(2) and Q(4) have been identified, primarily arising from the warm gas physically associated with the strong continuum source SgrB2(M). In addition, to determine the local excitation temperature we analyzed data from nine ro-vibrational transitions of the main isotopologue CCC in the frequency range between 1.6-1.9 THz which were taken from the Herschel Science Data Archive, and derived a gas excitation temperature of Tex = 44.4(+4.7/-3.9) K and a total column density of N(CCC)=3.88(+0.39/-0.35)x10^15 cm^-2. Assuming the excitation temperatures of C13CC and 13CCC to be the same as for CCC, we obtained column densities of the 13C-isotopologues of N(C13CC) = 2.1(+0.9/-0.6)X10^14 cm^-2 and N(13CCC)=2.4(+1.2/-0.8)x10^14 cm^-2. The derived 12C/13C abundance ratio in the C3 molecules is 20.5(4.2), which is in agreement with the elemental ratio of 20, typically observed in SgrB2(M). However, we find the N(13CCC) / N(C13CC) ratio to be 1.2(0.1), which is shifted from the statistically expected value of 2. We propose that the discrepant abundance ratio arises due to the lower zero-point energy of C13CC which makes position exchange reaction converting 13CCC to C13CC energetically favorable.

[13CII] observations in several Galactic sources show that the fine-structure [12CII] emission is often optically thick (the optical depths around 1 to a few). The aim of our study is to test whether this also affects the [12CII] emission from nearby galaxies like the Large Magellanic Cloud (LMC). We observed three star-forming regions in the LMC with upGREAT on board SOFIA at the frequency of the [CII] line. The 4GHz band width covers all three hyperfine lines of [13CII] simultaneously. For the analysis, we combined the [13CII] F=1-0 and F=1-1 hyperfine components, as they do not overlap with the [12CII] line in velocity. Three positions in N159 and N160 show an enhancement of [13CII] compared to the abundance-ratio-scaled [12CII] profile. This is likely due to the [12CII] line being optically thick, supported by the fact that the [13CII] line profile is narrower than [12CII], the enhancement varies with velocity, and the peak velocity of [13CII] matches the [OI] 63um self-absorption. The [12CII] line profile is broader than expected from a simple optical depth broadening of the [13CII] line, supporting the scenario of several PDR components in one beam having varying [12CII] optical depths. The derived [12CII] optical depth at three positions (beam size of 14arcsec, corresponding to 3.4pc) is 1--3, which is similar to values observed in several Galactic sources shown in previous studies. If this also applies to distant galaxies, the [CII] intensity will be underestimated by a factor of approximately 2.

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

The D/H ratio in cometary water has been shown to vary between 1 and 3 times the Earth's oceans value, in both Oort cloud comets and Jupiter-family comets originating from the Kuiper belt. We present new sensitive spectroscopic observations of water isotopologues in the Jupiter-family comet 46P/Wirtanen carried out using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The derived D/H ratio of $(1.61 \pm 0.65) \times 10^{-4}$ is the same as in the Earth's oceans. Although the statistics are limited, we show that interesting trends are already becoming apparent in the existing data. A clear anti-correlation is seen between the D/H ratio and the active fraction, defined as the ratio of the active surface area to the total nucleus surface. Comets with an active fraction above 0.5 typically have D/H ratios in water consistent with the terrestrial value. These hyperactive comets, such as 46P/Wirtanen, require an additional source of water vapor in their coma, explained by the presence of subliming icy grains expelled from the nucleus. The observed correlation may suggest that hyperactive comets belong to a population of ice-rich objects that formed just outside the snow line, or in the outermost regions of the solar nebula, from water thermally reprocessed in the inner disk that was transported outward during the early disk evolution. The observed anti-correlation between the active fraction and the nucleus size seems to argue against the first interpretation, as planetesimals near the snow line are expected to undergo rapid growth. Alternatively, isotopic properties of water outgassed from the nucleus and icy grains may be different due to fractionation effects at sublimation. In this case, all comets may share the same Earth-like D/H ratio in water, with profound implications for the early solar system and the origin of Earth's oceans.

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

Understanding the matter cycle in the interstellar medium of galaxies from the assembly of clouds to star formation and stellar feedback remains an important and exciting field in comtemporary astrophysics. Many open questions regarding cloud and structure formation, the role of turbulence, and the relative importance of the various feedback processes can only be addressed with observations of spectrally resolved lines. We here stress the importance of two specific sets of lines: the finestructure lines of atomic carbon as a tracer of the dark molecular gas and mid-J CO lines as tracers of the warm, active molecular gas in regions of turbulence dissipation and feedback. The observations must cover a wide range of environments (i.e., physical conditions), which will be achieved by large scale surveys of Galactic molecular clouds, the Galactic Center, the Magellanic clouds, and nearby galaxies. To date, such surveys are completely missing and thus constitute an important science opportunity for the next decade and beyond. For the successful interpretation of the observations, it will be essential to combine them with results from (chemical) modelling and simulations of the interstellar medium.

Aims: We study the spatial distribution and kinematics of water emission in a ~64 pc$^2$ region of the Galactic Center (GC) around Sgr A*. We also analyze the water excitation to derive the physical conditions and water abundances in the CND and the `quiescent clouds'. Methods: We presented the integrated intensity maps of the ortho 1$_{10}-1_{01}$, and para 2$_{02}-1_{11}$ and 1$_{11}-0_{00}$ water transitions observed with the HIFI instrument on board Herschel. To study the water excitation we used ground state ortho and para H$_2^{18}$O transitions. In our study, we also used SPIRE continuum measurements of the CND. Using a non-LTE radiative transfer code, the water line profiles and dust continuum were modeled. We also used a rotating ring model to reproduce the CND kinematics represented by the PV diagram. Results: We identify the water emission arising from the CND, the Western Streamer, and the 20 and 50 km s$^{-1}$ clouds. The ortho water maps show absorption structures in the range of [-220,10] km s$^{-1}$. The PV diagram shows that the 2$_{02}-1_{11}$ H$_2$O emission traces the CND. We derive high X$_{H_2O}$ of $\sim$(0.1-1.3)$\times$10$^{-5}$, V$_t$ of 14-23 km s$^{-1}$ and T$_d$ of 15-45 K for the CND, and the lower X$_{\rm H_2O}$ of 4$\times$10$^{-8}$ and V$_t$ of 9 km s$^{-1}$ for the 20 km s$^{-1}$ cloud. Collisional excitation and dust effects are responsible for the water excitation in the southwest lobe of the CND and the 20 km s$^{-1}$ cloud, whereas only collisions can account for the water excitation in the northeast lobe of the CND. We propose that the water vapor in the CND is caused by grain sputtering by shocks of 10-20 km s$^{-1}$, with some contribution of high temperature and cosmic-ray chemistries plus a PDR chemistry. The low X$_{\rm H_2O}$ derived for the 20 km s$^{-1}$ cloud could be partially a consequence of the water freeze-out on grains.

Massive stars inject mechanical and radiative energy into the surrounding environment, which stirs it up, heats the gas, produces cloud and intercloud phases in the interstellar medium, and disrupts molecular clouds (the birth sites of new stars). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds. Theoretical studies predict that momentum injection by radiation should dominate that by stellar winds, but this has been difficult to assess observationally. Velocity-resolved large-scale images in the fine-structure line of ionized carbon ([C II]) provide an observational diagnostic for the radiative energy input and the dynamics of the interstellar medium around massive stars. Here we report observations of a one-square-degree region (about 7 parsecs in diameter) of Orion molecular core -- the region nearest to Earth that exhibits massive-star formation -- at a resolution of 16 arcseconds (0.03 parsecs) in the [C II] line at 1.9 terahertz (158 micrometres). The results reveal that the stellar wind originating from the massive star ${\theta}^{1}$ Orionis C has swept up the surrounding material to create a bubble roughly four parsecs in diameter with a 2,600-solar-mass shell, which is expanding at 13 kilometres per second. This finding demonstrates that the mechanical energy from the stellar wind is converted very efficiently into kinetic energy of the shell and causes more disruption of the Orion molecular core 1 than do photo-ionization and evaporation or future supernova explosions.

We present a [CII] 158um map of the entire M51 (including M51b) grand--design spiral galaxy observed with the FIFI-LS instrument on SOFIA. We compare the [CII] emission with the total far--infrared (TIR) intensity and star formation rate(SFR) surface density maps (derived using H_alpha and 24um emission) to study the relationship between [CII] and the star formation activity in a variety of environments within M51 on scales of 16" corresponding to ~660 pc. We find that [CII] and the SFR surface density are well correlated in the central, spiral arm, and inter-arm regions. The correlation is in good agreement with that found for a larger sample of nearby galaxies at kpc scales. We find that the SFR, and [CII] and TIR luminosities in M51 are dominated by the extended emission in M51's disk. The companion galaxy M51b, however, shows a deficit of [CII] emission compared with the TIR emission and SFR surface density, with [CII] emission detected only in the S-W part of this galaxy. The [CII] deficit is associated with an enhanced dust temperature in this galaxy. We interpret the faint [CII] emission in M51b to be a result of suppressed star formation in this galaxy, while the bright mid- and far-infrared emission, which drive the TIR and SFR values, are powered by other mechanisms. A similar but less pronounced effect is seen at the location of the black hole in M51's center. The observed [CII] deficit in M51b suggests that this galaxy is a valuable laboratory to study the origin of the apparent [CII] deficit observed in ultra-luminous galaxies.

The aim of our study is to investigate the physical properties of the star-forming interstellar medium (ISM) in the Large Magellanic Cloud (LMC) by separating the origin of the emission lines spatially and spectrally. Following Okada et al. (2015, Paper I), we investigate different phases of the ISM traced by carbon-bearing species in four star-forming regions in the LMC, and model the physical properties using the KOSMA-tau PDR model. We mapped 3--13 arcmin$^2$ areas in 30 Dor, N158, N160 and N159 along the molecular ridge of the LMC in [CII]158um with GREAT on board SOFIA, and in CO(2-1) to (6-5), $^{13}$CO(2-1) and (3-2), [CI]3P1-3P0 and 3P2-3P1 with APEX. In all four star-forming regions, the line profiles of CO, $^{13}$CO, and [CI] emission are similar, whereas [CII] typically shows wider line profiles or an additional velocity component. For selected positions in N159 and 30 Dor, we observed the velocity-resolved [OI] 145um and 63um lines for the first time with upGREAT. At some positions, the [OI] line profiles match those of CO, at other positions they are more similar to the [CII] profiles. We interpret the different line profiles of CO, [CII] and [OI] as contributions from spatially separated clouds and/or clouds in different physical phases, which give different line ratios depending on their physical properties. We model the emission from the CO, [CI], [CII], and [OI] lines and the far-infrared continuum emission using the latest KOSMA-tau PDR model, which treats the dust-related physics consistently and computes the dust continuum SED together with the line emission of the chemical species. We find that the line and continuum emissions are not well-reproduced by a single clump ensemble. Toward the CO peak at N159~W, we propose a scenario that the CO, [CII], and [OI] 63um emission are weaker than expected because of mutual shielding among clumps. (abridged for arXiv)

The CCAT-prime telescope is a 6-meter aperture, crossed-Dragone telescope, designed for millimeter and sub-millimeter wavelength observations. It will be located at an altitude of 5600 meters, just below the summit of Cerro Chajnantor in the high Atacama region of Chile. The telescope's unobscured optics deliver a field of view of almost 8 degrees over a large, flat focal plane, enabling it to accommodate current and future instrumentation fielding >100k diffraction-limited beams for wavelengths less than a millimeter. The mount is a novel design with the aluminum-tiled mirrors nested inside the telescope structure. The elevation housing has an integrated shutter that can enclose the mirrors, protecting them from inclement weather. The telescope is designed to co-host multiple instruments over its nominal 15 year lifetime. It will be operated remotely, requiring minimum maintenance and on-site activities due to the harsh working conditions on the mountain. The design utilizes nickel-iron alloy (Invar) and carbon-fiber-reinforced polymer (CFRP) materials in the mirror support structure, achieving a relatively temperature-insensitive mount. We discuss requirements, specifications, critical design elements, and the expected performance of the CCAT-prime telescope. The telescope is being built by CCAT Observatory, Inc., a corporation formed by an international partnership of universities. More information about CCAT and the CCAT-prime telescope can be found at www.ccatobservatory.org.

We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.

The central area (40"x40") of the bipolar nebula S106 was mapped in the OI line at 63.2 micron with high angular (6") and spectral resolution, using GREAT on board SOFIA. The OI emission distribution is compared to the CO 16-15, CII 158 micron, and CO 11-10 lines, mm-molecular lines, and continuum. It is composed of several velocity components in the range from -30 km/s to 25 km/s. The high-velocity blue- and redshifted emission can be explained as arising from accelerated photodissociated (PDR) gas associated with a dark lane close to the massive binary system S106 IR, and from shocks caused by the stellar wind and/or a disk--envelope interaction. At velocities from -9 to -4 km/s and 0.5 to 8 km/s line wings are observed that we attribute to cooling in PDRs created by the ionizing radiation impinging on the cavity walls. The bulk velocity range is dominated by PDR emission from the clumpy molecular cloud. Modelling the emission in the different velocity ranges with the KOSMA-tau code constrains a radiation field chi of a few times 10^4 and densities n of a few times 10^4 cm^-3. Considering self-absorption of the OI line results in higher densities (up to 10^6 cm^-3) only for the gas component seen at high blue- and red velocities. The dark lane has a mass of 275 Msun and shows a velocity difference of 1.4 km/s along its projected length of 1 pc, determined from H13CO+ 1-0 mapping. It can be interpreted as a massive accretion flow, or the remains of it, linked to S106 IR/FIR. The most likely explanation is that the binary system is at a stage of its evolution where gas accretion is counteracted by the stellar winds and radiation, leading to the very complex observed spatial and kinematic emission distribution of the various tracers.

The methylidyne radical CH is commonly used as a proxy for H$_2$ in the cold, neutral phase of the interstellar medium. The optical spectroscopy of CH is limited by interstellar extinction, whereas far-infrared observations provide an integral view through the Galaxy. While the HF ground state absorption, another H$_2$ proxy in diffuse gas, frequently suffers from saturation, CH remains transparent both in spiral-arm crossings and high-mass star forming regions, turning this light hydride into a universal surrogate for H$_2$. However, in slow shocks and in regions dissipating turbulence its abundance is expected to be enhanced by an endothermic production path, and the idea of a "canonical" CH abundance needs to be addressed. The $N=2\leftarrow 1$ ground state transition of CH at $\lambda 149\,\mu$m has become accessible to high-resolution spectroscopy thanks to GREAT aboard SOFIA. Its unsaturated absorption and the absence of emission makes it an ideal candidate for the determination of column densities with a minimum of assumptions. Here we present an analysis of four sightlines towards distant, far-infrared bright Galactic star forming regions. If combined with the sub-millimeter line of CH at $\lambda 560\,\mu$m, environments forming massive stars can be analyzed. For this we present a case study on the "proto-Trapezium" cluster W3 IRS5, and demonstrate that the sub-millimeter/far-infrared lines of CH reliably trace not only diffuse but also dense, molecular gas. While we confirm the global correlation between the column densities of HF and those of CH, clear signposts of an over-abundance of CH are observed towards lower densities. A quiescent ion-neutral chemistry alone cannot account for this over-abundance. Vortices forming in turbulent, diffuse gas may be the setting for an enhanced production path.

We present a systematic survey of multiple velocity-resolved H$_2$O spectra using Herschel/HIFI towards nine nearby actively star forming galaxies. The ground-state and low-excitation lines (E$_{\rm up}\,\le 130\,{\rm K}$) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines ($130\,{\rm K}\, \le\, E_{\rm up}\,\le\,350\,{\rm K}$) typically display line shapes similar to CO. We analyze the HIFI observation together with archival SPIRE/PACS H$_2$O data using a state-of-the-art 3D radiative transfer code which includes the interaction between continuum and line emission. The water excitation models are combined with information on the dust- and CO spectral line energy distribution to determine the physical structure of the interstellar medium (ISM). We identify two ISM components that are common to all galaxies: A warm ($T_{\rm dust}\,\sim\,40-70\,{\rm K}$), dense ($n({\rm H})\,\sim\,10^5-10^6\,{\rm cm^{-3}}$) phase which dominates the emission of medium-excitation H$_2$O lines. This gas phase also dominates the FIR emission and the CO intensities for $J_{\rm up} > 8$. In addition a cold ($T_{\rm dust}\,\sim\,20-30\,{\rm K}$), dense ($n({\rm H})\sim\,10^4- 10^5\,{\rm cm^{-3}}$) more extended phase is present. It outputs the emission in the low-excitation H$_2$O lines and typically also produces the prominent line absorption features. For the two ULIRGs in our sample (Arp 220 and Mrk 231) an even hotter and more compact (R$_s\,\le\,100$ pc) region is present which is possibly linked to AGN activity. We find that collisions dominate the water excitation in the cold gas and for lines with $E_{\rm up}\le300\,{\rm K}$ and $E_{\rm up}\le800\,{\rm K}$ in the warm and hot component, respectively. Higher energy levels are mainly excited by IR pumping.

We report on the detection of the ground-state rotational line of ortho-D2H+ at 1.477 THz (203 micron) using the German REceiver for Astronomy at Terahertz frequencies (GREAT) onboard the Stratospheric Observatory For Infrared Astronomy (SOFIA). The line is seen in absorption against far-infrared continuum from the protostellar binary IRAS 16293-2422 in Ophiuchus. The para-D2H+ line at 691.7 GHz was not detected with the APEX telescope toward this position. These D2H+ observations complement our previous detections of para-H2D+ and ortho-H2D+ using SOFIA and APEX. By modeling chemistry and radiative transfer in the dense core surrounding the protostars, we find that the ortho-D2H+ and para-H2D+ absorption features mainly originate in the cool (T<18 K) outer envelope of the core. In contrast, the ortho-H2D+ emission from the core is significantly absorbed by the ambient molecular cloud. Analyses of the combined D2H+ and H2D+ data result in an age estimate of ~500 000 yr for the core, with an uncertainty of ~200 000 yr. The core material has probably been pre-processed for another 500 000 years in conditions corresponding to those in the ambient molecular cloud. The inferred time scale is more than ten times the age of the embedded protobinary. The D2H+ and H2D+ ions have large and nearly equal total (ortho+para) fractional abundances of ~$10^{-9}$ in the outer envelope. This confirms the central role of H3+ in the deuterium chemistry in cool, dense gas, and adds support to the prediction of chemistry models that also D3+ should be abundant in these conditions.

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

The physical state of the gas in the central 500 pc of NGC~5128 (the radio galaxy Centaurus A - Cen A), was investigated using the far-infrared fine-structure lines of carbon, oxygen, and nitrogen, as well as the CO(4-3) molecular line. The circumnuclear disk (CND) is traced by emission from dust and the neutral gas ([CI] and CO). A gas outflow with a line-of-sight velocity of 60 km/s is evident in both species. The center of the CND is bright in [OI], [OIII], and [CII]; [OI]63mu emission dominates that of [CII] even though it is absorbed with optical depths of 1.0-1.5. The outflow is well-traced by the [NII] and [NIII] lines and also seen in the [CII] and [OIII] lines that peak in the center. Ionized gas densities are moderate in the CND and low everywhere else. Neutral gas densities range from 4000 per cm3 (outflow, extended thin disk ETD) to 20 000 per cm3 (CND). The CND radiation field is weak compared to the ETD starburst field. The outflow has a much stronger radiation field. The total mass of all the CND gas is 9 x 10^(7) M(o) and the mass of the outflowing gas is only 15%-30% of that. The outflow most likely originates from the shock-dominated CND cavity surrounding the central black hole. With a factor of three uncertainty, the mass outflow rate is about 2 M(o)/yr, a thousand times higher than the accretion rate of the black hole. Without replenishment, the CND will be depleted in 15-120 million years. However, the outflow velocity is well below the escape velocity.

We present a new multi-pixel high resolution (R >10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receiver uses 2 x 7-pixel subarrays in orthogonal polarization, each in an hexagonal array around a central pixel. We present the first results for this new instrument after commissioning campaigns in May and December 2015 and after science observations performed in May 2016 . The receiver is designed to ultimately cover the full 1.8-2.5 THz frequency range but in its first implementation, the observing range was limited to observations of the [CII] line at 1.9 THz in 2015 and extended to 1.83-2.07 THz in 2016. The instrument sensitivities are state-of-the-art and the first scientific observations performed shortly after the commissioning confirm that the time efficiency for large scale imaging is improved by more than an order of magnitude as compared to single pixel receivers. An example of large scale mapping around the Horsehead Nebula is presented here illustrating this improvement. The array has been added to SOFIA's instrument suite already for ongoing observing cycle 4.

We present large-scale submm observations towards the Sgr A Complex covering ~300 arcmin2. These data were obtained in the frame of the HEXGAL Program with the Herschel-HIFI satellite and are complemented with submm observations obtained with the NANTEN2/SMART telescope as part of the NANTEN2/SMART Central Nuclear Zone Survey. The observed species are CO(4-3) observed with the NANTEN2/SMART telescope, and [CI](1-0), [CI](2-1), [NII](1-0), and [CII](3/2-1/2) observed with the Herschel-HIFI satellite. The observations are presented in a 1 km/s spectral resolution and a spatial resolution ranging from 46" to 28". The spectral coverage of the three lower frequency lines is +-200 km/s, while in the two high frequency lines, the upper LSR velocity limit is +94 km/s and +145 km/s for the [NII] and [CII] lines, respectively. The spatial distribution of the emission in all lines is very widespread. The bulk of the CO emission is found towards Galactic latitudes below the Galactic plane, and all the known molecular clouds are identified. Both neutral atomic carbon lines have their brightest emission associated with the +50 km/s cloud. Their spatial distribution at this LSR velocity describes a crescent-shape structure, which is probably the result of interaction with the energetic event (one or several supernovae explosions) that gave origin to the non-thermal Sgr A-East source. The [CII] and [NII] emissions have most of their flux associated with the thermal Arched-Filaments and the H Region and bright spots in [CII] emission towards the Central Nuclear Disk (CND) are detected. Warm Gas at very high (|Vlsr| > 100 km/s) LSR velocities is also detected towards the line of sight to the Sgr A Complex, and it is most probably located outside the region, in the X1 orbits.

Aims: We investigate how much of the [CII] emission in the nucleus of the nearby spiral galaxy IC 342 is contributed by PDRs and by the ionized gas. We examine the spatial variations of starburst/PDR activity and study the correlation of the [CII] line with the [NII] 205{\textmu}m emission line coming exclusively from the HII regions. Methods: We present small maps of [CII] and [NII] lines recently observed with the GREAT receiver on board SOFIA. In particular we present a super-resolution method to derive how unresolved, kinematically correlated structures in the beam contribute to the observed line shapes. Results: We find that the emission coming from the ionized gas shows a kinematic component in addition to the general Doppler signature of the molecular gas. We interpret this as the signature of two bi-polar lobes of ionized gas expanding out of the galactic plane. We then show how this requires an adaptation of our understanding of the geometrical structure of the nucleus of IC~342. Examining the starburst activity we find ratios $I([CII])/I(^{12}\mathrm{CO} (1-0))$ between 400 and 1800 in energy units. Applying predictions from numerical models of HII and PDR regions to derive the contribution from the ionized phase to the total [CII] emission we find that 35-90% of the observed [CII] intensity stems from the ionized gas if both phases contribute. Averaged over the central few hundred parsec we find for the [CII] contribution a HII-to-PDR ratio of 70:30. Conclusions: The ionized gas in the center of IC 342 contributes more strongly to the overall [CII] emission than is commonly observed on larger scales and than is predicted. Kinematic analysis shows that the majority of the [CII] emission is related to the strong but embedded star formation in the nuclear molecular ring and only marginally emitted from the expanding bi-polar lobes of ionized gas.

The CDMS was founded 1998 to provide in its catalog section line lists of molecular species which may be observed in various astronomical sources using radio astronomy. The line lists contain transition frequencies with qualified accuracies, intensities, quantum numbers, as well as further auxilary information. They have been generated from critically evaluated experimental line lists, mostly from laboratory experiments, employing established Hamiltonian models. Seperate entries exist for different isotopic species and usually also for different vibrational states. As of December 2015, the number of entries is 792. They are available online as ascii tables with additional files documenting information on the entries. The Virtual Atomic and Molecular Data Centre was founded more than 5 years ago as a common platform for atomic and molecular data. This platform facilitates exchange not only between spectroscopic databases related to astrophysics or astrochemistry, but also with collisional and kinetic databases. A dedicated infrastructure was developed to provide a common data format in the various databases enabling queries to a large variety of databases on atomic and molecular data at once. For CDMS, the incorporation in VAMDC was combined with several modifications on the generation of CDMS catalog entries. Here we introduce related changes to the data structure and the data content in the CDMS. The new data scheme allows us to incorporate all previous data entries but in addition allows us also to include entries based on new theoretical descriptions. Moreover, the CDMS entries have been transferred into a mySQL database format. These developments within the VAMDC framework have been driven by the needs of the astronomical community to be able to deal efficiently with large data sets obtained with the Herschel Space Telescope or, more recently, with the Atacama Large Millimeter Array.

This paper presents [CII], [CI] and CO emission line maps of the star-forming regions N66, N25+N26, and N88 in the metal-poor Local Group dwarf galaxy SMC. The spatial and velocity structure of the large HII region N66 reveals an expanding ring of shocked molecular gas centered on the exciting star cluster NGC346, whereas a more distant dense molecular cloud is being eroded by UV radiation from the same cluster. In the N25+N26 and N88 maps, diffuse [CII] emission at a relatively low surface brightness extends well beyond the compact boundaries of the bright emission associated with the [CII] regions. In all regions, the distribution of this bright [CII] emission and the less prominent [CI] emission closely follows the outline of the CO complexes, but the intensity of the [CII] and [CI] emission is generally anticorrelated, which can be understood by the action of photodissociation and photoionization processes. Notwithstanding the overall similarity of CO and [CII] maps, the intensity ratio of these lines varies significantly, mostly due to changes in CO brightness. [CII] emission line profiles are up to 50% wider in velocity than corresponding CO profiles. A radiative transfer analysis shows that the [CII] line is the dominant tracer of (CO-dark) molecular hydrogen in the SMC. CO emission traces only a minor fraction of the total amount of gas. The similarity of the spatial distribution and line profile shape, and the dominance of molecular gas associated with [[CII] rather than CO emission imply that in the low-metallicity environment of the SMC the small amount of dense molecular gas traced by CO is embedded in the much more extended molecular gas traced only by [CII] emission. The contribution from neutral atomic and ionized hydrogen zones is negligible in the star-forming regions observed.

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

We present the upGREAT THz heterodyne arrays for far-infrared astronomy. The Low Frequency Array (LFA) is designed to cover the 1.9-2.5 THz range using 2x7-pixel waveguide-based HEB mixer arrays in a dual polarization configuration. The High Frequency Array (HFA) will perform observations of the [OI] line at ~4.745 THz using a 7-pixel waveguide-based HEB mixer array. This paper describes the common design for both arrays, cooled to 4.5 K using closed- cycle pulse tube technology. We then show the laboratory and telescope characterization of the first array with its 14 pixels (LFA), which culminated in the successful commissioning in May 2015 aboard the SOFIA airborne observatory observing the [CII] fine structure transition at 1.905 THz. This is the first successful demonstration of astronomical observations with a heterodyne focal plane array above 1 THz and is also the first time high- power closed-cycle coolers for temperatures below 4.5 K are operated on an airborne platform.

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

The [CII]158um line is one of the dominant cooling lines in star-forming active regions. The commonly assumed clumpy UV-penetrated cloud models predict a [CII] line profile similar to that of CO. However, recent spectral-resolved observations show that they are often very different, indicating a more complex origin of the line emission including the dynamics of the source region. The aim of our study is to investigate the physical properties of the star-forming ISM in the Large Magellanic Cloud (LMC) by separating the origin of the emission lines spatially and spectrally. In this paper, we focus on the spectral characteristics and the origin of the emission lines, and the phases of carbon-bearing species in the N159 star-forming region in the LMC. We mapped a 4'x(3-4)' region in N159 in [CII]158um and [NII]205um with the GREAT on board SOFIA, and in CO(3-2), (4-3), (6-5), 13CO(3-2), and [CI]3P1-3P0 and 3P2-3P1 with APEX. The emission of all transitions observed shows a large variation in the line profiles across the map and between the different species. At most positions the [CII] emission line profile is substantially wider than that of CO and [CI]. We estimated the fraction of the [CII] integrated line emission that cannot be fitted by the CO line profile to be 20%-50%. We derived the relative contribution from C+, C, and CO to the column density in each velocity bin. The contribution from C+ dominates the velocity range far from the velocities traced by the dense molecular gas, and the region located between the CO cores of N159 W and E. We estimate the contribution of the ionized gas to the [CII] emission using the ratio to the [NII] emission to be < 19% to the [CII] emission at its peak position, and <15% over the whole observed region. Using the integrated line intensities, we present the spatial distribution of I([CII])/I(FIR). (abridged for arXiv)

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.

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.

In the framework of the Herschel/PRISMAS Guaranteed Time Key Program, the line of sight to the distant ultracompact HII region W51e2 has been observed using several selected molecular species. Most of the detected absorption features are not associated with the background high-mass star-forming region and probe the diffuse matter along the line of sight. We present here the detection of an additional narrow absorption feature at ~70 km/s in the observed spectra of HDO, NH3 and C3. The 70 km/s feature is not uniquely identifiable with the dynamic components (the main cloud and the large-scale foreground filament) so-far identified toward this region. The narrow absorption feature is similar to the one found toward low-mass protostars, which is characteristic of the presence of a cold external envelope. The far-infrared spectroscopic data were combined with existing ground-based observations of 12CO, 13CO, CCH, CN, and C3H2 to characterize the 70 km/s component. Using a non-LTE analysis of multiple transitions of NH3 and CN, we estimated the density (n(H2) (1-5)x10^5 cm^-3) and temperature (10-30 K) for this narrow feature. We used a gas-grain warm-up based chemical model with physical parameters derived from the NH3 data to explain the observed abundances of the different chemical species. We propose that the 70 km/s narrow feature arises in a dense and cold clump that probably is undergoing collapse to form a low-mass protostar, formed on the trailing side of the high-velocity filament, which is thought to be interacting with the W51 main cloud. While the fortuitous coincidence of the dense clump along the line of sight with the continuum-bright W51e2 compact HII region has contributed to its non-detection in the continuum images, this same attribute makes it an appropriate source for absorption studies and in particular for ice studies of star-forming regions.

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.

We present new CO and C^o line measurements of the compact circumnuclear disk in the center of NGC 128 (Centaurus~A) obtained with the Herschel Space Observatory, as well as SEST, JCMT, and APEX. The Cen A center CO ladder is quite different from those of either star-burst galaxies or AGNs. In addition, the relative intensity of the central Cen A [CI] emission lines is much greate than that in any other galaxy. The CO surface brightness of the compact circumnuclear disk (CND) is significantly higher than that of the much more extended thin disk (ETD) in the same line of sight. Our LVG and PDR/XDR models suggest that much of the CND gas is relatively cool (25 - 80 K) and not very dense (~ 300 cm^{-3}) if the heating is by UV photons, although there is some gas in both the CND and the ETD with a much higher density of ~30 000 cm^{-3}. Finally, there is also high-excitation, high-density phase in the CND (but not in the ETD), either in the form of an extreme PDR but more likely in the form of an XDR. The total gas mass of the Cen A CND is 8.4 x 10^{7} M(sun), uncertain by a factor of two. The CO-H2 conversion factor is 4 x 10^{20} K km/s, also within a factor of two.

We present the first results from a new carbon monoxide (CO) survey of the southern Galactic plane being conducted with the Mopra radio telescope in Australia. The 12CO, 13CO and C18O J=1-0 lines are being mapped over the l = 305-345 deg, b = +/- 0.5 deg portion of the 4th quadrant of the Galaxy, at 35" spatial and 0.1 km/s spectral resolution. The survey is being undertaken with two principal science objectives: (i) to determine where and how molecular clouds are forming in the Galaxy and (ii) to probe the connection between molecular clouds and the "missing" gas inferred from gamma-ray observations. We describe the motivation for the survey, the instrumentation and observing techniques being applied, and the data reduction and analysis methodology. In this paper we present the data from the first degree surveyed, l = 323-324 deg, b = +/- 0.5 deg. We compare the data to the previous CO survey of this region and present metrics quantifying the performance being achieved; the rms sensitivity per 0.1 km/s velocity channel is ~1.5K for 12CO and ~0.7K for the other lines. We also present some results from the region surveyed, including line fluxes, column densities, molecular masses, 12CO/13CO line ratios and 12CO optical depths. We also examine how these quantities vary as a function of distance from the Sun when averaged over the 1 square degree survey area. Approximately 2 x 10E6 MSun of molecular gas is found along the G323 sightline, with an average H2 number density of nH2 ~ 1 cm-3 within the Solar circle. The CO data cubes will be made publicly available as they are published.

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 present the results of CO emission line observations toward the semicircular Galactic supernova remnant (SNR) 3C434.1 (G94.0+1.0). We mapped an area covering the whole SNR in the 12CO J=1-0 emission line using the Seoul Radio Astronomy Observatory (SRAO) 6-m telescope, and found a large molecular cloud superposed on the faint western part of the SNR. The cloud was elongated along the north-south direction and showed a very good spatial correlation with the radio features of the SNR. We carried out 12CO J=2-1 line observations of this cloud using the K\"{o}lner Observatorium f\"{u}r Sub-Millimeter Astronomie (KOSMA) 3-m telescope and found a region in which the 12CO J=2-1 to J=1-0 ratio was high (~ 1.6). This higher excitation, together with the morphological relation, strongly suggested that the molecular cloud was interacting with the SNR. The systemic velocity of the molecular cloud (-13 km/s) gave a kinematic distance of 3.0 kpc to the SNR-molecular cloud system. We derived the physical parameters of the SNR based on this new distance. We examined the variation of the radio spectral index over the remnant and found that it was flatter in the western part, wherein the SNR was interacting with the molecular cloud. We therefore propose that 3C434.1 is the remnant of a supernova explosion that occurred just outside the boundary of a relatively thin, sheetlike molecular cloud. We present a hydrodynamic model showing that its asymmetric radio morphology can result from its interaction with this blocking molecular cloud.

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.

The catalog section of the Cologne Database for Molecular Spectroscopy, CDMS, contains mostly rotational transition frequencies, with auxiliary information, of molecules observable in space. The frequency lists are generated mostly from critically evaluated laboratory data employing established Hamiltonian models. The CDMS has been online publicly for more than 12 years, e.g., via the short-cut http://www.cdms.de. Initially constructed as ascii tables, its inclusion into a database environment within the Virtual Atomic and Molecular Data Centre (VAMDC, http://www.vamdc.eu) has begun in June 2008. A test version of the new CDMS is about to be released. The CDMS activities have been part of the extensive laboratory spectroscopic investigations in Cologne. Moreover, these activities have also benefit from collaborations with other laboratory spectroscopy groups as well as with astronomers. We will provide some basic information on the CDMS and its participation in the VAMDC project. In addition, some recent detections of molecules as well as spectroscopic studies will be discussed to evaluate the spectroscopic data needs of Herschel, SOFIA, and ALMA in particular in terms of light hydrides, complex molecules, and metal containing species

(Abridged) We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654--1897 GHz using HIFI on board Herschel, in DR21(OH), a high mass star forming region. Several transitions of CCH and c-C3H2 have also been observed with HIFI and the IRAM 30m telescope. A gas and grain warm-up model was used to identify the primary C3 forming reactions in DR21(OH). We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is \sim6x10^{-10} and \sim3x10^{-9} respectively. For CCH and c-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3, and c-C3H2 abundances are most consistent with a chemical model with n(H2)\sim5x10^{6} cm^-3 post-warm-up dust temperature, T_max =30 K and a time of \sim0.7-3 Myr. Post warm-up gas phase chemistry of CH4 released from the grain at t\sim 0.2 Myr and lasting for 1 Myr can explain the observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. The observed lower C3 abundance in MM1 as compared to MM2 and the envelope could be indicative of destruction of C3 in the more evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived for DR21(OH) in other studies.

(Abridged) The 30 Dor region in the Large Magellanic Cloud (LMC) is the most vigorous star-forming region in the Local Group. Star formation in this region is taking place in low-metallicity molecular gas that is exposed to an extreme far--ultraviolet (FUV) radiation field powered by the massive compact star cluster R136. We used the NANTEN2 telescope to obtain high-angular resolution observations of the 12CO 4-3, 7-6, and 13CO 4-3 rotational lines and [CI] 3P1-3P0 and 3P2-3P1 fine-structure submillimeter transitions in 30Dor-10, the brightest CO and FIR-emitting cloud at the center of the 30Dor region. We derived the properties of the low-metallicity molecular gas using an excitation/radiative transfer code and found a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. We compared the derived properties with those in the N159W region, which is exposed to a more moderate far-ultraviolet radiation field compared with 30Dor-10, but has similar metallicity. We also combined our CO detections with previously observed low-J CO transitions to derive the CO spectral-line energy distribution in 30Dor-10 and N159W. The separate excitation analysis of the submm CO lines and the neutral carbon fine structure lines shows that the mid-J CO and [CI]-emitting gas in the 30Dor-10 region has a temperature of about 160 K and a H2 density of about 10^4 cm^-3. We find that the molecular gas in 30Dor-10 is warmer and has a lower beam filling factor compared to that of N159W, which might be a result of the effect of a strong FUV radiation field heating and disrupting the low--metallicity molecular gas. We use a clumpy PDR model (including the [CII] line intensity reported in the literature) to constrain the FUV intensity to about chi_0 ~ 3100 and an average total H density of the clump ensemble of about 10^5 cm^-3 in 30Dor-10.

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

Supplementing the publications based on the first-light observations with the German Receiver for Astronomy at Terahertz frequencies (GREAT) on SOFIA, we present background information on the underlying heterodyne detector technology. We describe the superconducting hot electron bolometer (HEB) detectors that are used as frequency mixers in the L1 (1400 GHz), L2 (1900 GHz), and M (2500 GHz) channels of GREAT. Measured performance of the detectors is presented and background information on their operation in GREAT is given. Our mixer units are waveguide-based and couple to free-space radiation via a feedhorn antenna. The HEB mixers are designed, fabricated, characterized, and flight-qualified in-house. We are able to use the full intermediate frequency bandwidth of the mixers using silicon-germanium multi-octave cryogenic low-noise amplifiers with very low input return loss. Superconducting HEB mixers have proven to be practical and sensitive detectors for high-resolution THz frequency spectroscopy on SOFIA. We show that our niobium-titanium-nitride (NbTiN) material HEBs on silicon nitride (SiN) membrane substrates have an intermediate frequency (IF) noise roll-off frequency above 2.8 GHz, which does not limit the current receiver IF bandwidth. Our mixer technology development efforts culminate in the first successful operation of a waveguide-based HEB mixer at 2.5 THz and deployment for radioastronomy. A significant contribution to the success of GREAT is made by technological development, thorough characterization and performance optimization of the mixer and its IF interface for receiver operation on SOFIA. In particular, the development of an optimized mixer IF interface contributes to the low passband ripple and excellent stability, which GREAT demonstrated during its initial successful astronomical observation runs.

Methods: We used the dual-band receiver GREAT on board the SOFIA airborne telescope to perform observations of the [C II] 158 {\mu}m fine-structure line at the postitions of two giant molecular clouds (GMC) in the center of IC 342 (GMCs C and E) and compared the spectra with corresponding ground-based data for low- and mid-J CO and [C I]. We performed model calculations assuming a clumpy photo-dissociation region (PDR) environment using the KOSMA-tau PDR model code to derive physical parameters of the local medium. Results: The [C II] 158 {\mu}m emission resembles the spectral signature of ground-based atomic and molecular lines, which indicates a common origin. The emission from GMC E can be decomposed into a cool, molecular component with weak far-ultraviolet (FUV) fields and low, mean densities of 103 cm^-3 and a strongly excited starburst/PDR region with higher densities of 104 cm^-3 and FUV intensities of 250-300 Draine fields. The emission from GMC C is consistent with gas densities of 5000 cm^-3, FUV intensities of a few Draine fields and total gas masses of 20\times10^6 M$_\odot$. Conclusions: The high spectral resolution of the GREAT receiver allowed us to decompose the [C II] emission of the GMC E into a strongly excited gas component resembling a PDR/starburst environment and a quieter, less excited gas component and to analyze the different components within a single beam individually.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

The BD+40\degree4124 region was observed with high angular and spectral resolution with the German heterodyne instrument GREAT in CO J = 13 \rightarrow 12 and [CII] on SOFIA. These observations show that the [CII] emission is very strong in the reflection nebula surrounding the young Herbig Ae/Be star BD+40\degree4124. A strip map over the nebula shows that the [CII] emission approximately coincides with the optical nebulosity. The strongest [CII] emission is centered on the B2 star and a deep spectrum shows that it has faint wings, which suggests that the ionized gas is expanding. We also see faint CO J = 13 \rightarrow 12 at the position of BD+40\degree4124, which suggests that the star may still be surrounded by an accretion disk.We also detected [CII] emission and strong CO J = 13 \rightarrow 12 toward V1318 Cyg. Here the [CII] emission is fainter than in BD+40\degree4124 and appears to come from the outflow, since it shows red and blue wings with very little emission at the systemic velocity, where the CO emission is quite strong. It therefore appears that in the broad ISO beam the [CII] emission was dominated by the reflection nebula surrounding BD+40\degree4124, while the high J CO lines originated from the adjacent younger and more deeply embedded binary system V1318 Cyg.

We report SOFIA/GREAT, Herschel/HIFI, and ground-based velocity-resolved spectroscopy of carbon monoxide (CO) rotational transitions from J=2-1 to J=16-15 toward two positions in the circum-nuclear disk (CND) in our Galactic center. Radiative transfer models were used to derive information on the physical state of the gas traced by CO. The excitation of the CO gas cannot be explained by a single physical component, but is clearly the superposition of various warm gas phases. In a two-component approach, our large velocity gradient (LVG) analysis suggests high temperatures of ~200 K with moderate gas densities of only ~10^4.5 cm^-3 for the bulk of the material. A higher excited phase, carrying ~20-30% of the column densities, is warmer (~300-500 K) but only slightly denser ~10^5.3 cm^-3. These densities are too low to self-stabilize the clumps against their high internal turbulence and fall below the Roche density (>10^7 cm^-3) at 1.5 pc galactocentric distance. We conclude that the bulk of the material in the CND is not organized by self-gravity nor stable against tidal disruption, and must be transient.

Context. The HII region/PDR/molecular cloud complex S106 is excited by a single O-star. The full extent of the warm and dense gas close to the star has not been mapped in spectrally resolved high-J CO or [CII] lines, so the kinematics of the warm, partially ionized gas, are unknown. Whether the prominent dark lane bisecting the hourglass-shaped nebula is due solely to the shadow cast by a small disk around the exciting star or also to extinction in high column foreground gas was an open question until now. Aims. To disentangle the morphology and kinematics of warm neutral and ionized gas close to the star, study their relation to the bulk of the molecular gas, and to investigate the nature of the dark lane. Methods. We use the heterodyne receiver GREAT on board SOFIA to observe velocity resolved spectral lines of [C II] and CO 11-10 in comparison with so far unpublished submm continuum data at 350 micron (SHARC-II) and complementary molecular line data. Results. The high angular and spectral resolution observations show a very complex morphology and kinematics of the inner S106 region, with many different components at different excitation conditions contributing to the observed emission. The [C II] lines are found to be bright and very broad, tracing high velocity gas close to the interface of molecular cloud and HII region. CO 11-10 emission is more confined, both spatially and in velocity, to the immediate surroundings of S106 IR showing the presence of warm, high density (clumpy) gas. Our high angular resolution submm continuum observations rule out the scenario where the dark lane separating the two lobes is due solely to the shadow cast by a small disk close to the star. The lane is clearly seen also as warm, high column density gas at the boundary of the molecular cloud and HII region.

We describe the design and construction of GREAT, the German REceiver for Astronomy at Terahertz frequencies operated on the Stratospheric Observatory for Infrared Astronomy (SOFIA). GREAT is a modular dual-color heterodyne instrument for highresolution far-infrared (FIR) spectroscopy. Selected for SOFIA's Early Science demonstration, the instrument has successfully performed three Short and more than a dozen Basic Science flights since first light was recorded on its April 1, 2011 commissioning flight. We report on the in-flight performance and operation of the receiver that - in various flight configurations, with three different detector channels - observed in several science-defined frequency windows between 1.25 and 2.5 THz. The receiver optics was verified to be diffraction-limited as designed, with nominal efficiencies; receiver sensitivities are state-of-the-art, with excellent system stability. The modular design allows for the continuous integration of latest technologies; we briefly discuss additional channels under development and ongoing improvements for Cycle 1 observations. GREAT is a principal investigator instrument, developed by a consortium of four German research institutes, available to the SOFIA users on a collaborative basis.

Context. Observing supernova remnants (SNRs) and modelling the shocks they are associated with is the best way to quantify the energy SNRs re-distribute back into the Interstellar Medium (ISM). Aims. We present comparisons of shock models with CO observations in the F knot of the W28 supernova remnant. These comparisons constitute a valuable tool to constrain both the shock characteristics and pre-shock conditions. Methods. New CO observations from the shocked regions with the APEX and SOFIA telescopes are presented and combined. The integrated intensities are compared to the outputs of a grid of models, which were combined from an MHD shock code that calculates the dynamical and chemical structure of these regions, and a radiative transfer module based on the 'large velocity gradient' (LVG) approximation. Results. We base our modelling method on the higher J CO transitions, which unambiguously trace the passage of a shock wave. We provide fits for the blue- and red-lobe components of the observed shocks. We find that only stationary, C-type shock models can reproduce the observed levels of CO emission. Our best models are found for a pre-shock density of 104 cm-3, with the magnetic field strength varying between 45 and 100 {\mu}G, and a higher shock velocity for the so-called blue shock (\sim25 km s-1) than for the red one (\sim20 km s-1). Our models also satisfactorily account for the pure rotational H2 emission that is observed with Spitzer.

Context: We analyze the NGC 2024 HII region and molecular cloud interface using [12CII] and [13CII] observations. Aims: We attempt to gain insight into the physical structure of the interface layer between the molecular cloud and the HII region. Methods. Observations of [12CII] and [13CII] emission at 158 {\mu}m with high spatial and spectral resolution allow us to study the detailed structure of the ionization front and estimate the column densities and temperatures of the ionized carbon layer in the PDR. Results: The [12CII] emission closely follows the distribution of the 8 mum continuum. Across most of the source, the spectral lines have two velocity peaks similar to lines of rare CO isotopes. The [13CII] emission is detected near the edge-on ionization front. It has only a single velocity component, which implies that the [12CII] line shape is caused by self-absorption. An anomalous hyperfine line-intensity ratio observed in [13CII] cannot yet be explained. Conclusions: Our analysis of the two isotopes results in a total column density of N(H)~1.6\times10^23 cm^-2 in the gas emitting the [CII] line. A large fraction of this gas has to be at a temperature of several hundred K. The self-absorption is caused by a cooler (T<=100 K) foreground component containing a column density of N(H)~10^22 cm^-2.

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.

We investigate the gas dynamics and the physical properties of photodissociation regions (PDRs) in IC1396A, which is an illuminated bright-rimmed globule with internal structures created by young stellar objects. Our mapping observations of the [CII] emission in IC1396A with GREAT onboard SOFIA revealed the detailed velocity structure of this region. We combined them with observations of the [CI] 3P_1 - 3P_0 and CO(4-3) emissions to study the dynamics of the different tracers and physical properties of the PDRs. The [CII] emission generally matches the IRAC 8 micron, which traces the polycyclic aromatic hydrocarbon (PAH) emissions. The CO(4-3) emission peaks inside the globule, and the [CI] emission is strong in outer regions, following the 8 micron emission to some degree, but its peak is different from that of [CII]. The [CII] emitting gas shows a clear velocity gradient within the globule, which is not significant in the [CI] and CO(4-3) emission. Some clumps that are prominent in [CII] emission appear to be blown away from the rim of the globule. The observed ratios of [CII]/[CI] and [CII]/CO(4-3) are compared to the KOSMA-tau PDR model, which indicates a density of 10^4-10^5 cm-3.

The GREAT observations need frequency-selective calibration across the passband for the residual atmospheric opacity at flight altitude. At these altitudes the atmospheric opacity has both narrow and broad spectral features. To determine the atmospheric transmission at high spectral resolution, GREAT compares the observed atmospheric emission with atmospheric model predictions, and therefore depends on the validity of the atmospheric models. We discusse the problems identified in this comparison with respect to the observed data and the models, and describe the strategy used to calibrate the science data from GREAT/SOFIA during the first observing periods.

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.

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.

RX J1713.7-3946 is one of the TeV {\gamma}-ray supernova remnants (SNRs) emitting synchrotron X rays. The SNR is associated with molecular gas located at ~1 kpc. We made new molecular observations toward the dense cloud cores, peaks A, C and D, in the SNR in the 12CO(J=2-1) and 13CO(J=2-1) transitions at angular resolution of 90". The most intense core in 13CO, peak C, was also mapped in the 12CO(J=4-3) transition at angular resolution of 38". Peak C shows strong signs of active star formation including bipolar outflow and a far-infrared protostellar source and has a steep gradient with a r^{-2.2$\pm$0.4} variation in the average density within radius r. Peak C and the other dense cloud cores are rim-brightened in synchrotron X rays, suggesting that the dense cloud cores are embedded within or on the outer boundary of the SNR shell. This confirms the earlier suggestion that the X rays are physically associated with the molecular gas (Fukui et al. 2003). We present a scenario where the densest molecular core, peak C, survived against the blast wave and is now embedded within the SNR. Numerical simulations of the shock-cloud interaction indicate that a dense clump can indeed survive shock erosion, since shock propagation speed is stalled in the dense clump. Additionally, the shock-cloud interaction induces turbulence and magnetic field amplification around the dense clump that may facilitate particle acceleration in the lower-density inter-clump space leading to the enhanced synchrotron X rays around dense cores.

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

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.

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

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

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

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

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

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

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

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

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 present observations of the rotational ortho-water ground transition, the two lowest para-water transitions, and the ground transition of ionised ortho-water in the archetypal starburst galaxy M82, performed with the HIFI instrument on the Herschel Space Observatory. These observations are the first detections of the para-H2O(111-000) (1113\,GHz) and ortho-H2O+(111-000) (1115\,GHz) lines in an extragalactic source. All three water lines show different spectral line profiles, underlining the need for high spectral resolution in interpreting line formation processes. Using the line shape of the para-H2O(111-000) and ortho-H2O+(111-000) absorption profile in conjunction with high spatial resolution CO observations, we show that the (ionised) water absorption arises from a ~2000 pc^2 region within the HIFI beam located about ~50 pc east of the dynamical centre of the galaxy. This region does not coincide with any of the known line emission peaks that have been identified in other molecular tracers, with the exception of HCO. Our data suggest that water and ionised water within this region have high (up to 75%) area-covering factors of the underlying continuum. This indicates that water is not associated with small, dense cores within the ISM of M82 but arises from a more widespread diffuse gas component.

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

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

We present spectrally resolved observations of triatomic carbon (C3) in several ro-vibrational transitions between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654-1897 GHz along the sight-lines to the submillimeter continuum sources W31C and W49N, using Herschel's HIFI instrument. We detect C3 in absorption arising from the warm envelope surrounding the hot core, as indicated by the velocity peak position and shape of the line profile. The sensitivity does not allow to detect C3 absorption due to diffuse foreground clouds. From the column densities of the rotational levels in the vibrational ground state probed by the absorption we derive a rotation temperature (T_rot) of ~50--70 K, which is a good measure of the kinetic temperature of the absorbing gas, as radiative transitions within the vibrational ground state are forbidden. It is also in good agreement with the dust temperatures for W31C and W49N. Applying the partition function correction based on the derived T_rot, we get column densities N(C3) ~7-9x10^{14} cm^{-2} and abundance x(C3)~10^{-8} with respect to H2. For W31C, using a radiative transfer model including far-infrared pumping by the dust continuum and a temperature gradient within the source along the line of sight we find that a model with x(C3)=10^{-8}, T_kin=30-50 K, N(C3)=1.5 10^{15} cm^{-2} fits the observations reasonably well and provides parameters in very good agreement with the simple excitation analysis.

We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.

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

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

Fukui et al. (2006) discovered two molecular loops in the Galactic center and argued that the foot points of the molecular loops, two bright spots at both loops ends, represent the gas accumulated by the falling motion along the loops, subsequent to magnetic flotation by the Parker instability. We have carried out sensitive CO observations of the foot points toward l=356 deg at a few pc resolution in the six rotational transitions of CO; 12CO(J=1-0, 3-2, 4-3, 7-6), 13CO(J=1-0) and C18O(J=1-0). The high resolution image of 12CO (J=3-2) has revealed the detailed distribution of the high excitation gas including U shapes, the outer boundary of which shows sharp intensity jumps accompanying strong velocity gradients. An analysis of the multi-J CO transitions shows that the temperature is in a range from 30-100 K and density is around 10^3-10^4 cm^-3, confirming that the foot points have high temperature and density although there is no prominent radiative heating source such as high mass stars in or around the loops. We argue that the high temperature is likely due to the shock heating under C-shock condition caused by the magnetic flotation. We made a comparison of the gas distribution with theoretical numerical simulations and note that the U shape is consistent with numerical simulations. We also find that the region of highest temperature of ~100 K or higher inside the U shape corresponds to the spur having an upward flow, additionally heated up either by magnetic reconnection or bouncing in the interaction with the narrow neck at the bottom of the U shape. We note these new findings further reinforce the magnetic floatation interpretation.

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

We report on the detection with the HIFI instrument on board the Herschel satellite of the two hydrogen chloride isotopologues, H35Cl and H37Cl, towards the massive star-forming region W3A. The J=1-0 line of both species was observed with receiver 1b of the HIFI instrument at 625.9 and 624.9 GHz. The different hyperfine components were resolved. The observations were modeled with a non-local, non-LTE radiative transfer model that includes hyperfine line overlap and radiative pumping by dust. Both effects are found to play an important role in the emerging intensity from the different hyperfine components. The inferred H35Cl column density (a few times 1e14 cm^-2), and fractional abundance relative to H nuclei (~7.5e^-10), supports an upper limit to the gas phase chlorine depletion of ~200. Our best-fit model estimate of the H35Cl/H37Cl abundance ratio is ~2.1+/-0.5, slightly lower, but still compatible with the solar isotopic abundance ratio (~3.1). Since both species were observed simultaneously, this is the first accurate estimation of the [35Cl]/[37Cl] isotopic ratio in molecular clouds. Our models indicate that even for large line opacities and possible hyperfine intensity anomalies, the H35Cl and H37Cl J=1-0 integrated line-intensity ratio provides a good estimate of the 35Cl/37Cl isotopic abundance ratio.

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

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

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

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

NGC 2024, a sites of massive star formation, have complex internal structures caused by cal heating by young stars, outflows, and stellar winds. These complex cloud structures lead to intricate emission line shapes. The goal of this paper is to show that the complex line shapes of 12 CO lines in NGC 2024 can be explained consistently with a model, whose temperature and velocity structure are based on the well-established scenario of a PDR and the Blister model. We present velocity-resolved spectra of seven CO lines ranging from J=3 to J=13, and we combined these data with CO high-frequency data from the ISO satellite. We find that the bulk of the molecular cloud associated with NGC 2024 consists of warm (75 K) and dense (9e5 cm-3) gas. An additional hot (~ 300 K) component, located at the interface of the HII region and the molecular cloud, is needed to explain the emission of the high-J CO lines. Deep absorption notches indicate that very cold material (20 K) exists in front of the warm material, too. A temperature and column density structure consistent with those predicted by PDR models, combined with the velocity structure of a Blister model, appropriately describes the observed emission line profiles of this massive star forming region. This case study of NGC 2024 shows that, with physical insights into these complex regions and careful modeling, multi-line observations of CO can be used to derive detailed physical conditions in massive star forming regions.

To date the onset of large-scale star formation in galaxies and its link to gravitational stability of the galactic disk have not been fully understood. The nearby face-on spiral galaxy M51 is an ideal target for studying this subject. This paper combines CO, dust, HI, and stellar maps of M51 and its companion galaxy to study the H2/HI transition, the gas-to-dust ratios, and the stability of the disk against gravitational collapse. We combine maps of the molecular gas using 12CO 2--1 map HERA/IRAM-30m data and HI VLA data to study the total gas surface density and the phase transition of atomic to molecular gas. The total gas surface density is compared to the dust surface density from 850 micron SCUBA data. Taking into account the velocity dispersions of the molecular and atomic gas, and the stellar surface densities derived from the 2MASS K-band survey, we derive the total Toomre Q parameter of the disk. The gas surface density in the spiral arms is approximately 2-3 higher compared to that of the interarm regions. The ratio of molecular to atomic surface density shows a nearly power-law dependence on the hydrostatic pressure P_hydro. The gas surface density distribution in M51 shows an underlying exponential distribution with a scale length of h_gas=7.6 kpc representing 55% of the total gas mass, comparable to the properties of the exponential dust disk. In contrast to the velocity widths observed in HI, the CO velocity dispersion shows enhanced line widths in the spiral arms compared to the interarm regions. The contribution of the stellar component in the Toomre Q-parameter analysis is significant and lowers the combined Q-parameter Q_tot by up to 70% towards the threshold for gravitational instability. The value of Q_tot varies from 1.5-3 in radial averages. A map of Q_tot shows values around 1 on the spiral arms.

Context: The earliest phases of massive star formation are currently much debated. Aims. In an effort to make progress, we took a census of Class0-like protostellar dense cores in the NGC 3576 region, one of the nearest and most luminous embedded sites of high-mass star formation in the Galaxy. Methods: We used the P-ArTeMiS bolometer camera on the APEX telescope to produce the first 450-micron dust continuum map of the filamentary dense clump associated with NGC 3576. Results: Combining our 450-micron observations with existing data at other wavelengths, we have identified seven massive protostellar sources along the NGC 3576 filament and placed them in the M_env - L_bol evolutionary diagram for protostars. Conclusions: Comparison with theoretical evolutionary tracks suggests that these seven protostellar sources will evolve into massive stars with masses M* ~ 15-50 Msun. Four sources are classified as candidate high-mass Class 0 objects, two sources as massive Class I objects, and one source appears to be at an intermediate stage.

Deuterated ions are abundant in cold (T=10 K), dense (n=10^5 cm^-3) regions, in which CO is frozen out onto dust grains. In such environments, the deuterium fractionation of such ions can exceed the elemental abundance ratio of D/H by a factor of 10^4. In this paper we use the deuterium fractionation to investigate the evolutionary state of Class 0 protostars. In a sample of 20 protostellar objects, we found a clear correlation between the N2D+/N2H+ ratio and evolutionary tracers. As expected, the coolest, i.e. the youngest, objects show the largest deuterium fractionation. Furthermore, we find that sources with a high N2D+/N2H+ ratio show clear indication for infall.

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.

Studying molecular gas in the central regions of the star burst galaxies NGC4945 and Circinus enables us to characterize the physical conditions and compare them to previous local and high-z studies. We estimate temperature, molecular density and column densities of CO and atomic carbon. Using model predictions we give a range of estimated CO/C abundance ratios. Using the new NANTEN2 4m sub-millimeter telescope in Pampa La Bola, Chile, we observed for the first time CO 4-3 and [CI] 3P1-3 P0 at the centers of both galaxies at linear scale of 682 pc and 732 pc respectively. We compute the cooling curves of 12CO and 13CO using radiative transfer models and estimate the physical conditions of CO and [CI]. The centers of NGC4945 and Circinus are very [CI] bright objects, exhibiting [CI] 3P1 - 3 P0 luminosities of 91 and 67Kkms-1kpc2, respectively. The [CI] 3P1-3 P0/CO 4-3 ratio of integrated intensities are large at 1.2 in NGC4945 and 2.8 in Circinus. Combining previous CO J= 1-0, 2-1 and 3-2 and 13CO J= 1-0, 2-1 studies with our new observations, the radiative transfer calculations give a range of densities, n(H2) = 10^3-3*104^cm-3, and a wide range of kinetic temperatures, Tkin = 20 - 100K, depending on the density. Future CO J= 7-6 and [CI] 2-1 observations will be important to resolve the ambiguity in the physical conditions and confirm the model predictions.

The Carina region is an excellent astrophysical laboratory for studying the feedback mechanisms of newly born, very massive stars within their natal giant molecular clouds (GMCs) at only 2.35 kpc distance. We use a clumpy PDR model to analyse the observed intensities of atomic carbon and CO and to derive the excitation conditions of the gas. The NANTEN2-4m submillimeter telescope was used to map the [CI] 3P1-3P0, 3P2-3P1 and CO 4-3, 7-6 lines in two 4'x4' regions of Carina where molecular material interfaces with radiation from the massive star clusters. One region is the northern molecular cloud near the compact OB cluster Tr14, and the second region is in the molecular cloud south of etaCar and Tr16. These data were combined with 13CO SEST spectra, HIRES/IRAS 60um and 100um maps of the FIR continuum, and maps of 8um IRAC/Spitzer and MSX emission. We used the HIRES far-infrared dust data to create a map of the FUV field heating the gas. The northern region shows an FUV field of a few 1000 in Draine units while the field of the southern region is about a factor 10 weaker. We constructed models consisting of an ensemble of small spherically symmetric PDR clumps within the 38" beam (0.43pc), which follow canonical power-law mass and mass-size distributions. We find that an average local clump density of 2x10**5 cm-3 is needed to reproduce the observed line emission at two selected interface positions. Stationary, clumpy PDR models reproduce the observed cooling lines of atomic carbon and CO at two positions in the Carina Nebula.

Star formation laws, like i.e. the Schmidt law relating star formation rate and total gas density, have been studied in several spiral galaxies but the underlying physics are not yet well understood. M51, as a nearby face-on, grand design spiral galaxy studied in many line transitions, is an ideal target to study the connection between physical conditions of the gas and star formation activity. In this contribution we combine molecular, atomic, total gas and stellar surface densities and study the gravitational stability of the gas (Schuster et al.2007, Hitschfeld et al. in prep.). From our IRAM-30m 12 CO2-1 map and complementary HI-, Radio Continuum- and ACS-HST B-band-data we derive maps of the total gas density and the stellar surface density to study the gravitational stability of the gas via the Toomre Q parameter.

We aim to study dust properties of massive star forming regions in the outer Galaxy, in a direction opposite to the Galactic center. We present observations of six outer Galaxy point sources IRAS 01045+6505, 01420+6401, 05271+3059, 05345+3556, 20222+3541 and 20406+4555, taken with the Submillimeter Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) at 450 and 850 micron. Single temperature greybody models are fitted to the Spectral Energy Distribution of the detected sub-mm cores to derive dust temperature, dust emissivity index and optical depth at 250 micron. The observed radial intensity profiles of the sub-mm cores were fitted with power laws to derive the indices describing the density distribution. At a resolution of 15" all six IRAS point sources show multiple emission peaks. Only four out of fourteen detected sub-mm cores show associated mid-infrared emission. For the sub-mm cores we derive dust temperatures of 32+-5 K and dust emissivity indices between 0.9 and 2.5. The density profiles of the sub-mm cores can be fitted by a single power law distribution with indices -1.5+-0.3, with most cores showing an index of -1.5. This is consistent with most observations of massive star forming regions and supports predictions of models of star formation which consider non-thermal support against gravitational collapse.

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.

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 THz atmospheric windows centered at roughly 1.3 and 1.5~THz, contain numerous spectral lines of astronomical importance, including three high-J CO lines, the N+ line at 205 microns, and the ground transition of para-H2D+. The CO lines are tracers of hot (several 100K), dense gas; N+ is a cooling line of diffuse, ionized gas; the H2D+ line is a non-depleting tracer of cold (~20K), dense gas. As the THz lines benefit the study of diverse phenomena (from high-mass star-forming regions to the WIM to cold prestellar cores), we have built the CO N+ Deuterium Observations Receiver (CONDOR) to further explore the THz windows by ground-based observations. CONDOR was designed to be used at the Atacama Pathfinder EXperiment (APEX) and Stratospheric Observatory For Infrared Astronomy (SOFIA). CONDOR was installed at the APEX telescope and test observations were made to characterize the instrument. The combination of CONDOR on APEX successfully detected THz radiation from astronomical sources. CONDOR operated with typical Trec=1600K and spectral Allan variance times of 30s. CONDOR's first light observations of CO 13-12 emission from the hot core Orion FIR4 (= OMC1 South) revealed a narrow line with T(MB) = 210K and delta(V)=5.4km/s. A search for N+ emission from the ionization front of the Orion Bar resulted in a non-detection. The successful deployment of CONDOR at APEX demonstrates the potential for making observations at THz frequencies from ground-based facilities.

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.

We present CI 3P1-3P0 spectra at four spiral arm positions and the nuclei of the nearby galaxies M83 and M51 obtained at the JCMT. This data is complemented with maps of CO 1-0, 2-1, and 3-2, and ISO/LWS far-infrared data of CII (158 micron), OI (63 micron), and NII (122 micron) allowing for the investigation of a complete set of all major gas cooling lines. From the intensity of the NII line, we estimate that between 15% and 30% of the observed CII emission originate from the dense ionized phase of the ISM. The analysis indicates that emission from the diffuse ionized medium is negligible. In combination with the FIR dust continuum, we find gas heating efficiencies below ~0.21% in the nuclei, and between 0.25 and 0.36% at the outer positions. Comparison with models of photon-dominated regions (PDRs) of Kaufman et al. (1999) with the standard ratios OI(63)/CII_PDR and (OI(63)+CII_PDR) vs. TIR, the total infrared intensity, yields two solutions. The physically most plausible solution exhibits slightly lower densities and higher FUV fields than found when using a full set of line ratios, CII_PDR/CI(1-0), CI(1-0)/CO(3-2), CO(3-2)/CO(1-0), CII/CO(3-2), and, OI(63)/CII_PDR. The best fits to the latter ratios yield densities of 10^4 cm^-3 and FUV fields of ~G_0=20-30 times the average interstellar field without much variation. At the outer positions, the observed total infrared intensities are in perfect agreement with the derived best fitting FUV intensities. The ratio of the two intensities lies at 4-5 at the nuclei, indicating the presence of other mechanisms heating the dust.

We used the KOSMA 3m telescope to map the core 7'x5' of the Galactic massive star forming region W3Main in the two fine structure lines of atomic carbon and four mid-J transitions of CO and 13CO. The maps are centered on the luminous infrared source IRS5 for which we obtained ISO/LWS data comprising four high-J CO transitions, CII, and OI at 63 and 145 micron. In combination with a KAO map of integrated line intensities of CII (Howe et al. 1991), this data set allows to study the physical structure of the molecular cloud interface regions where the occurence of carbon is believed to change from C+ to C0, and to CO. The molecular gas in W3Main is warmed by the far ultraviolet (FUV) field created by more than a dozen OB stars. Detailed modelling shows that most of the observed line intensity ratios and absolute intensities are consistent with a clumpy photon dominated region (PDR) of a few hundred unresolved clumps per 0.84pc beam, filling between 3 and 9% of the volume, with a typical clump radius of 0.025pc (2.2"), and typical mass of 0.44Msun. The high-excitation lines of CO stem from a 100-200K layer, as also the CI lines. The bulk of the gas mass is however at lower temperatures.

We have used the KOSMA 3m telescope to map the core 7'x5' of the Galactic massive star forming region W3Main in the two fine structure lines of atomic carbon and four mid-J transitions of CO and 13CO. In combination with a map of singly ionized carbon (Howe et al. 1991), and FIR fine structure line data observed by ISO/LWS at the center position, these data sets allow to study in detail the physical structure of the photon dominated cloud interface regions (PDRs) where the occurance of carbon changes from CII to CI, and to CO.

Observations of Carbon bearing species are among the most important diagnostic probes of ongoing star formation. CO is a surrogate for H$_2$ and is found in the vicinity of star formation sites. There, [CI] emission is thought to outline the dense molecular cores and extend into the lower density regions, where the impinging interstellar UV radiation field plays a critical role for the dissociation and ionization processes. Emission of ionized carbon ([CII]) is found to be even more extended than [CI] and is linking up with the ionized medium. These different tracers emphasize the importance of multi-wavelength studies to draw a coherent picture of the processes driving and driven by high mass star formation. Until now, large scale surveys were only done with low resolution, such as the COBE full sky survey, or were biased to a few selected bright sources (e.g. Yamamoto et al. 2001, Schneider et al. 2003). A broader basis of unbiased, high-resolution observations of [CI], CO, and [CII] may play a key role to probe the material processed by UV radiation.

We have used the Antarctic Submillimeter Telescope and Remote Observatory (AST/RO), a 1.7m diameter single-dish submillimeter-wave telescope at the geographic South Pole, to determine the physical state of gas in the Galactic Center region and assess its stability. We present an analysis based on data obtained as part of an ongoing AST/RO key project: the large-scale mapping of the dominant cooling lines of the molecular interstellar medium in the Milky Way. These data are released for general use.

Several galaxies, such as dwarfs and Irregulars as well as outer galactic clouds have low metallicity. At low metallicities a reduction in the amount of dust and heavy elements plays a significant role on the chemistry as well as the heating and cooling of the gas in the molecular regions, called as Photon Dominated Regions (PDRs). We present here the effects of reduced metallicity in PDRs and study the important PDR cooling lines ([CII], [CI] and CO). Moreover many observational evidences suggest that molecular clouds are clumpy. We model the molecular emission from galaxies incorporating a mass spectrum of clumps. We also compare our results with the semi-analytical results obtained by Bolatto et al. (1999).

We present 12CO(J=2->1) and CO(J=3->2) observations of quiescent Virgo Cluster spiral galaxies with the KOSMA 3m submm telescope. The beam sizes of 80" at 345 GHz and 120" at 230 GHz are well suited for the investigation of global properties of Virgo Cluster galaxies. The observed sample was selected based on previous 12CO(J=1->0) detections by Stark et al. (1986), performed with the AT&T Bell Laboratory 7m telescope (beam size 100"). We were able to detect 18 spiral galaxies in 12CO(2->1) and 16 in 12CO(3->2). Beam matched observations of the lowest three 12CO transitions allow us to compare our results with previous high spatial resolution studies of (moderate) starburst galaxies and galactic core regions. We discuss the global excitation conditions of the ISM in these quiescent spiral galaxies. The resulting CO (3--2)/(1--0) integrated line ratios vary over a relatively narrow range of values from 0.35 to 0.14 (on a K km/s-scale) with increasing CO (2--1)/(1--0) ratio (from 0.5 to 1.1). The line ratios between the three lowest rotational transitions of CO cannot be fitted by any radiative transfer model with a single source component. A two-component model, assuming a warm, dense nuclear and a cold, less dense disc component allows us to fit the observed line ratios for most of the galaxies individually by selecting suitable parameters. The two-component model, however, fails to explain the observed correlation of the line ratios. This is due to a variation of the relative filling factor of the warm gas alone, assuming a typical set of parameters for the two components common for all galaxies.

We report measurements of the 12C/13C abundance ratio in the three galactic regions G 333.0-0.4, NGC 6334 A and G 351.6-1.3 from observations of the 12CI 3P2-3P1 transition and the hyperfine components of the corresponding 13CI transition near 809 GHz. These transitions were observed simultaneously with the CO 7-6 line emission at 806 GHz with the AST/RO telescope located at the South Pole. From a simultaneous fit to the 12CI 3P2-3P1 transition and the HF components of the corresponding 13CI transition and an independent estimate of an upper limit to the optical depth of the 12CI emission we determine intrinsic 12CI/13CI column density ratios of 23+-1 for G 333.0-0.4, 56+-14 for NGC 6334 A and 69+-12 for G 351.6-1.3. As the regions observed are photon dominated, we argue that the apparent enhancement in the abundance of 13C towards G 333.0-0.4 may be due to strong isotope-selective photodissociation of 13CO, outweighing the effects of chemical isotopic fractionation as suggested by models of PDRs. Towards NGC 6334 A and G 351.6-1.3 these effects appear to be balanced, similar to the situation for the Orion Bar region observed by Keene et al. (1998).

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.

We report the detection of the (J,K) = (1,1), (2,2), and (3,3) inversion lines of ammonia (NH3) towards the south--western molecular lobe in M82. The relative intensities of the ammonia lines are characterized by a rotational temperature of T_rot=29+/-5 K which implies an average kinetic temperature of T_kin~60 K. A Gaussian decomposition of the observed spectra indicates increasing kinetic temperatures towards the nucleus of M82, consistent with recent findings based on CO observations. The observations imply a very low NH3 abundance relative to H2, X(NH3)~5x10^(-10). We present evidence for a decreasing NH3 abundance towards the central active regions in M82 and interpret this abundance gradient in terms of photodissociation of NH3 in PDRs. The low temperature derived here from NH3 also explains the apparent underabundance of complex molecules like CH_3OH and HNCO, which has previously been reported.

AST/RO, a 1.7 m diameter telescope for astronomy and aeronomy studies at wavelengths between 200 and 2000 microns, was installed at the South Pole during the 1994-1995 Austral summer. The telescope operates continuously through the Austral winter, and is being used primarily for spectroscopic studies of neutral atomic carbon and carbon monoxide in the interstellar medium of the Milky Way and the Magellanic Clouds. The South Pole environment is unique among observatory sites for unusually low wind speeds, low absolute humidity, and the consistent clarity of the submillimeter sky. Four heterodyne receivers, an array receiver, three acousto-optical spectrometers, and an array spectrometer are installed. A Fabry-Perot spectrometer using a bolometric array and a Terahertz receiver are in development. Telescope pointing, focus, and calibration methods as well as the unique working environment and logistical requirements of the South Pole are described.

This paper is understood as a supplement to the paper by [Stutzki et al, 1998], where we have shown the usefulness of the Allan-variance and its higher dimensional generalization, the Delta-variance, for the characterization of molecular cloud structures. In this study we present the connection between the Allan- and Delta-variance and a more popular structure analysis tool: the wavelet transform. We show that the Allan- and Delta-variances are the variances of wavelet transform coefficients.

We report the first extragalactic detection of the neutral carbon [CI] 3P2-3P1 fine structure line at 809 GHz. The line was observed towards M82 simultaneously with the 3P1-3P0 line at 492 GHz, providing a precise measurement of the J=2-1/J=1-0 integrated line ratio of 0.96 (on a [K km s^-1] -scale). This ratio constrains the [CI] emitting gas to have a temperature of at least 50 K and a density of at least 10^4 cm^-3. Already at this minimum temperature and density, the beam averaged CI-column density is large, 2.1 10^18 cm^-2, confirming the high CI/CO abundance ratio of approximately 0.5 estimated earlier from the 492 GHz line alone. We argue that the [CI] emission from M82 most likely arises in clouds of linear size around a few pc with a density of about 10^4 cm^-3 or slightly higher and temperatures of 50 K up to about 100 K.