In recent years, a plethora of high spectral resolution observations of sub-mm and FIR transitions of methylidene (CH), have demonstrated this radical to be a valuable proxy for H2, that can be used for characterising molecular gas within the interstellar medium (ISM) on a Galactic scale, including the CO-dark component. Here we report the discovery of the 13CH isotopologue in the ISM using the upGREAT receiver on board SOFIA. We have detected the three hyperfine structure components of the 2THz frequency transition from its ground-state toward four high-mass star-forming regions and determine 13CH column densities. The ubiquity of molecules containing carbon in the ISM has turned the determination of the ratio between the abundances of carbon's two stable isotopes, 12C/13C, into a cornerstone for Galactic chemical evolution studies. Whilst displaying a rising gradient with Galactocentric distance, this ratio, when measured using observations of different molecules (CO, H2CO, and others) shows systematic variations depending on the tracer used. These observed inconsistencies may arise from optical depth effects, chemical fractionation or isotope-selective photo-dissociation. Formed from C+ either via UV-driven or turbulence-driven chemistry, CH reflects the fractionation of C+, and does not show any significant fractionation effects unlike other molecules previously used to determine the 12C/13C isotopic ratio which make it an ideal tracer for the 12C/13C ratio throughout the Galaxy. Therefore, by comparing the derived column densities of 13CH with previously obtained SOFIA data of the corresponding transitions of the main isotopologue 12CH, we derive 12C/13C isotopic ratios toward Sgr B2(M), G34.26+0.15, W49(N) and W51E. Adding our values derived from 12/13CH to previous calculations of the Galactic isotopic gradient we derive a revised value of 12C/13C = 5.85(0.50)R_GC + 15.03(3.40).

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

We have mapped the G287.84-0.82 cometary globule (with the Treasure Chest cluster embedded in it) in the South Pillars region of Carina (i) in [CII], 63micron [OI], and CO(11-10) using upGREAT on SOFIA and (ii) in J=2-1 transitions of CO, 13CO, C18O and J=3-2 transitions of H2CO using the APEX telescope in Chile. We probe the morphology, kinematics, and physical conditions of the molecular gas and the photon dominated regions (PDRs) in G287.84-0.82. The [CII] and [OI] emission suggest that the overall structure of the pillar (with red-shifted photo evaporating tails) is consistent with the effect of FUV radiation and winds from eta-Car and O stars in Trumpler 16. The gas in the head of the pillar is strongly influenced by the embedded cluster, whose brightest member is an O9.5V star, CPD-59 2661. The emission of the [CII] and [OI] lines peak at a position close to the embedded star, while all other tracers peak at another position lying to the north-east consistent with gas being compressed by the expanding PDR created by the embedded cluster. The molecular gas inside the globule is probed with the J=2-1 transitions of CO and isotopologues as well as H2CO, and analyzed using a non-LTE model (escape-probability approach), while we use PDR models to derive the physical conditions of the PDR. We identify at least two PDR gas components; the diffuse part (~10^4 cm^-3) is traced by [CII], while the dense (n~ 2-8x10^5 cm^-3) part is traced by [CII], [OI], CO(11-10). Using the F=2-1 transition of [13CII] detected at 50 positions in the region, we derive optical depths (0.9-5), excitation temperatures of [CII] (80-255 K), and N(C+) of 0.3-1x10^19 cm^-2. The total mass of the globule is ~1000 Msun, about half of which is traced by [CII]. The dense PDR gas has a thermal pressure of 10^7-10^8 K cm^-3, which is similar to the values observed in other regions.

We present Submillimeter Array (SMA) observations in the CO J=3-2, SiO J=5-4 and 8-7, and SO 9_8-8_7 lines, as well as Atacama Pathfinder EXperiment (APEX) observations in the CO J=6-5 line, of an extremely high-velocity and jet-like outflow in high-mass star-forming region HH 80--81. The outflow is known to contain two prominent molecular bullets, namely B1 and B2, discovered from our previous SMA CO J=2-1 observations. While B1 is detected in all the CO, SiO, and SO lines, B2 is only detected in CO lines. The CO 3-2/2-1 line ratio in B1 is clearly greater than that in B2. We perform a large velocity gradient analysis of the CO lines and derive a temperature of 70--210 K for B1 and 20--50 K for B2. Taking into account the differences in the velocity, distance from the central source, excitation conditions, and chemistry between the two bullets, we suggest that the bullets are better explained by direct ejections from the innermost vicinity of the central high-mass protostar, and that we are more likely observing the molecular component of a primary wind rather than entrained or swept-up material from the ambient cloud. These findings further support our previous suggestions that the molecular bullets indicate an episodic, disk-mediated accretion in the high-mass star formation process.

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.

High-mass stars are formed within massive molecular clumps, where a large number of stars form close together. The evolution of the clumps with different masses and luminosities is mainly regulated by its high-mass stellar content and the formation of such objects is still not well understood. In this work, we characterise the mid-J CO emission in a statistical sample of 99 clumps (Top100) selected from the ATLASGAL survey that are representative of the Galactic proto-cluster population. High-spatial resolution APEX-CHAMP+ maps of the CO(6-5) and CO(7-6) transitions were obtained and combined with additional single-pointing APEX-FLASH+ spectra of the CO(4-3) line. We study the correlations of the CO line luminosities and profiles for the three CO transitions with the clump properties and investigate if and how they change as a function of the evolution. All sources were detected above 3-$\sigma$ in all three CO transitions and most of the sources exhibit broad CO emission likely associated with molecular outflows. We found that the extension of the mid-J CO emission is correlated with the size of the dust emission traced by the Herschel-PACS 70 $\mu$m maps. The CO line luminosity is correlated with the luminosity and mass of the clumps. However, it does not correlate with the L/M ratio. The dependency of the CO luminosity with the properties of the clumps is steeper for higher-J transitions. Our data seem to exclude that this trend is biased by self-absorption features in the CO emission, but rather suggest that different J transitions arise from different regions of the inner envelope. Moreover, high-mass clumps show similar trends in CO luminosity as lower mass clumps, but are systematically offset towards larger values, suggesting that higher column density and/or temperature (of unresolved) CO emitters are found inside high-mass clumps.

Stretching a length reaching 10 pc projected in the plane of sky, the radio jet associated with Herbig-Haro objects 80 and 81 (HH 80-81) is known as the largest and best collimated protostellar jet in our Galaxy. The nature of the molecular outflow associated with this extraordinary jet remains an unsolved question which is of great interests to our understanding of the relationship between jets and outflows in high-mass star formation. Here we present Atacama Pathfinder EXperiment CO(6-5) and (7-6), James Clerk Maxwell Telescope CO(3-2), Caltech Submillimeter Observatory CO(2-1), and Submillimeter Array CO and $^{13}$CO(2-1) mapping observations of the outflow. We report on the detection of a two-component outflow consisting of a collimated component along the jet path and a wide-angle component with an opening angle of about $30^{\circ}$. The gas velocity structure suggests that each of the two components traces part of a primary wind. From LVG calculations of the CO lines, the outflowing gas has a temperature around 88 K, indicating that the gas is being heated by shocks. Based on the CO(6-5) data, the outflow mass is estimated to be a few $M_{\odot}$, which is dominated by the wide-angle component. A comparison between the HH 80-81 outflow and other well shaped massive outflows suggests that the opening angle of massive outflows continues to increase over time. Therefore, the mass loss process in the formation of early-B stars seems to be similar to that in low-mass star formation, except that a jet component would disappear as the central source evolves to an ultracompact HII region.

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.

We present Atacama Pathfinder EXperiment (APEX) observations toward the massive star-forming region G240.31+0.07 in the CO J = 3--2, 6--5, and 7--6 lines. We detect a parsec-sized, bipolar, and high velocity outflow in all the lines, which allow us, in combination with the existing CO J = 2--1 data, to perform a multi-line analysis of physical conditions of the outflowing gas. The CO 7--6/6--5, 6--5/3--2, and 6--5/2--1 ratios are found to be nearly constant over a velocity range of $\sim$5--25 km s$^{-1}$ for both blueshifted and redshifted lobes. We carry out rotation diagram and large velocity gradient (LVG) calculations of the four lines, and find that the outflow is approximately isothermal with a gas temperature of $\sim$50 K, and that the the CO column density clearly decreases with the outflow velocity. If the CO abundance and the velocity gradient do not vary much, the decreasing CO column density indicates a decline in the outflow gas density with velocity. By comparing with theoretical models of outflow driving mechanisms, our observations and calculations suggest that the massive outflow in G240.31+0.07 is being driven by a wide-angle wind and further support a disk mediated accretion at play for the formation of the central high-mass star.

We analyze a [C II] 158 micron map obtained with the L2 GREAT receiver on SOFIA of the emission/reflection nebula illuminated by the early B star S1 in the rho-OphA cloud core. This data set has been complemented with maps of CO(3-2), 13CO(3-2) and C18O(3-2), observed as a part of the JCMT Gould Belt Survey, with archival HCO^+(4-3) JCMT data, as well as with [O I] 63 and 145 micron imaging with Herschel/PACS. The [C II] emission is completely dominated by the strong PDR emission from the nebula surrounding S1 expanding into the dense Oph A molecular cloud west and south of S1. The [C II] emission is significantly blue shifted relative to the CO spectra and also relative to the systemic velocity, particularly in the northwestern part of the nebula. The [C II] lines are broader towards the center of the S1 nebula and narrower towards the PDR shell. The [C II] lines are strongly self-absorbed over an extended region in the S1 PDR. Based on the strength of the [13C II] F = 2-1 hyperfine component, [C II] is significantly optically thick over most of the nebula. CO and 13CO(3-2) spectra are strongly self-absorbed, while C18O(3-2) is single peaked and centered in the middle of the self-absorption. We have used a simple two-layer LTE model to characterize the background and foreground cloud contributing to the [C II] emission. From this analysis we estimate the extinction due to the foreground cloud to be ~9.9 mag, which is slightly less than the reddening estimated towards S1. Since some of the hot gas in the PDR is not traced by low J CO emission, this result appears quite plausible. Using a plane parallel PDR model with the observed [OI(145)]/[C II] brightness ratio and an estimated FUV intensity of 3100-5000 G0 suggests that the density of the [C II] emitting gas is ~3-4x10^3 cm^-3.

We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC and its immediate surroundings. This strongest emission component has a velocity similar to that of the $^{13}$CO(2-1) emission of the IRDC, but offset by $\sim3\:{\rm km\:s}^{-1}$ and with a larger velocity width of $\sim9\:{\rm km\:s}^{-1}$. The spatial distribution of the [CII] emission of this component is also offset predominantly to one side of the dense filamentary structure of the IRDC. The CII column density is estimated to be of the order of $\sim10^{17}-10^{18}\,{\rm cm}^{-2}$. We compare these results to the [CII] emission from numerical simulations of magnetized, dense gas filaments formed from giant molecular cloud (GMC) collisions, finding similar spatial and kinematic offsets. These observations and modeling of [CII] add further to the evidence that IRDC G035.39-00.33 has been formed by a process of GMC-GMC collision, which may thus be an important mechanism for initiating star cluster formation.

We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $\mu$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we derive the gas kinetic temperature and spatial density using the measured p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$, 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$, and 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ ratios. The gas kinetic temperatures derived from the p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$ and 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$ line ratios are high, ranging from 43 to $>$300 K with an unweighted average of 91 $\pm$ 4 K. Deduced $T_{\rm kin}$ values from the $J$ = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the p-H$_2$CO 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ line ratios yield 0.6-8.3 $\times$ 10$^6$ cm$^{-3}$ with an unweighted average of 1.5 ($\pm$0.1) $\times$ 10$^6$ cm$^{-3}$. A comparison of kinetic temperatures derived from p-H$_2$CO, NH$_3$, and the dust emission indicates that p-H$_2$CO traces a distinctly higher temperature than the NH$_3$ (2,2)/(1,1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H$_2$CO linewidths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H$_2$CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H$_2$CO increase with time through the evolution of the clumps.

Aims: To reveal the morphology, chemical composition, kinematics and to establish the main processes prevalent in the gas at the foot points of the giant molecular loops (GMLs) in the Galactic center region Methods: Using the 22-m Mopra telescope, we mapped the M$-3.8+0.9$ molecular cloud, placed at the foot points of a giant molecular loop, in 3-mm range molecular lines. To derive the molecular hydrogen column density, we also observed the $^{13}$CO $(2-1)$ line at 1 mm using the 12-m APEX telescope. From the 3 mm observations 12 molecular species were detected, namely HCO$^+$, HCN, H$^{13}$CN, HNC, SiO, CS, CH$_3$OH, N$_2$H$^+$, SO, HNCO, OCS, and HC$_3$N. Results: Maps revealing the morphology and kinematics of the M$-3.8+0.9$ molecular cloud in different molecules are presented. We identified six main molecular complexes. We derive fractional abundances in 11 selected positions of the different molecules assuming local thermodynamical equilibrium. Conclusions: Most of the fractional abundances derived for the M$-3.8+0.9$ molecular cloud are very similar over the whole cloud. However, the fractional abundances of some molecules show significant difference with respect to those measured in the central molecular zone (CMZ). The abundances of the shock tracer SiO are very similar between the GMLs and the CMZ. The methanol emission is the most abundant specie in the GMLs. This indicates that the gas is likely affected by moderate $\sim $ 30 km s$^{-1}$ or even high velocity (50 km s$^{-1}$) shocks, consistent with the line profile observed toward one of the studied position. The origin of the shocks is likely related to the flow of the gas throughout the GMLs towards the foot points.

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar.

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 report a discovery of shocked gas from the supernova remnant (SNR) G357.7+0.3. Our millimeter and submillimeter observations reveal broad molecular lines of CO(2-1), CO(3-2), CO(4-3), 13CO (2-1) and 13CO (3-2), HCO^+ and HCN using HHSMT, Arizona 12-Meter Telescope, APEX and MOPRA Telescope. The widths of the broad lines are 15-30 kms, and the detection of such broad lines is unambiguous, dynamic evidence showing that the SNR G357.7+0.3 is interacting with molecular clouds. The broad lines appear in extended regions (>4.5'x5'). We also present detection of shocked H2 emission in mid-infrared but lacking ionic lines using the Spitzer IRS observations to map a few arcmin area. The H2 excitation diagram shows a best-fit with a two-temperature LTE model with the temperatures of ~200 and 660 K. We observed [C II] at 158um and high-J CO(11-10) with the GREAT on SOFIA. The GREAT spectrum of [C II], a 3 sigma detection, shows a broad line profile with a width of 15.7 km/s that is similar to those of broad CO molecular lines. The line width of [C~II] implies that ionic lines can come from a low-velocity C-shock. Comparison of H2 emission with shock models shows that a combination of two C-shock models is favored over a combination of C- and J-shocks or a single shock. We estimate the CO density, column density, and temperature using a RADEX model. The best-fit model with n(H2) = 1.7x10^{4} cm^{-3}, N(CO) = 5.6x10^{16} cm^{-2}, and T = 75 K can reproduce the observed millimeter CO brightnesses.

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.

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

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.

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.

The fine-structure line of [OI] at 63micron is an important diagnostic tool in different fields of astrophysics. However, our knowledge of this line relies on observations with low spectral resolution, and the real contribution of each component (PDR, jet) in complex environment of star-forming regions (SFRs) is poorly understood. We investigate the contribution of jet and PDR emission, and of absorption to the [OI]63micron line towards the ultra-compact H{\sc ii} region G5.89--0.39 and study its far-IR line luminosity in different velocity regimes through [OI], [CII], CO, OH, and H2O. We mapped G5.89--0.39 in [OI] and in CO(16--15) with the GREAT receiver onboard SOFIA. We observed the central position of the source in the OH^2\Pi_{3/2}, J=5/2\toJ=3/2 and ^2\Pi_{1/2}, J=3/2\toJ=1/2 lines. These data were complemented with APEX CO(6-5) and CO(7-6) and HIFI maps and single-pointing observations in [CII], H2O, and HF. The [OI] spectra in G5.89--0.39 are severely contaminated by absorptions from the envelope and from different clouds along the line of sight. Emission is detected only at HV, clearly associated with the compact north-south outflows traced by extremely HV low-J CO. The mass-loss rate and energetics of derived from [OI] agree well with estimates from CO, suggesting that the molecular outflows in G5.89--0.39 are driven by the jet system seen in [OI]. The far-IR line luminosity of G5.89--0.39 is dominated by [OI] at HV; the second coolant in this velocity regime is CO, while [CII], OH and H2O are minor contributors to the cooling in the outflow. Our study shows the importance of spectroscopically resolved data of [OI]63micron for using this line as diagnostic of SFRs. While this was not possible until now, the GREAT receiver onboard SOFIA has recently opened the possibility of detailed studies of this line to investigate its potential for probing different environments.

OH is a key molecule in H2O chemistry, a valuable tool for probing physical conditions, and an important contributor to the cooling of shock regions. OH participates in the re-distribution of energy from the protostar towards the surrounding ISM. Our aim is to assess the origin of the OH emission from the Cepheus A massive star-forming region and to constrain the physical conditions prevailing in the emitting gas. We thus want to probe the processes at work during the formation of massive stars. We present spectrally resolved observations of OH towards the outflows of Cepheus A with the GREAT spectrometer onboard the SOFIA telescope. Three triplets were observed at 1834.7 GHz, 1837.8 GHz, and 2514.3 GHz (163.4, 163.1, and 119.2 microns), at angular resolutions of 16.3", 16.3", and 11.9", respectively. We present the CO (16-15) spectrum at the same position. We compared the integrated intensities in the redshifted wings to shock models. The two triplets near 163 microns are detected in emission with blending hyperfine structure unresolved. Their profiles and that of CO can be fitted by a combination of 2 or 3 Gaussians. The observed 119.2 microns triplet is seen in absorption, since its blending hyperfine structure is unresolved, but with three line-of-sight components and a blueshifted emission wing consistent with that of the other lines. The OH line wings are similar to those of CO, suggesting that they emanate from the same shocked structure. Under this common origin assumption, the observations fall within the model predictions and within the range of use of our model only if we consider that four shock structures are caught in our beam. Our comparisons suggest that the observations might be consistently fitted by a J-type model with nH > 1e5 cm-3, v > 20 km/s, and with a filling factor of ~1. Such a high density is generally found in shocks associated to high-mass protostars.

The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

We have mapped the NGC 2023 reflection nebula in [CII] and CO(11--10) with the heterodyne receiver GREAT on SOFIA and obtained slightly smaller maps in 13CO(3--2), CO(3--2), CO(4--3), CO(6--5), and CO(7--6) with APEX in Chile. We use these data to probe the morphology, kinematics, and physical conditions of the C II region, which is ionized by FUV radiation from the B2 star HD37903. The [CII] emission traces an ellipsoidal shell-like region at a position angle of ~ -50 deg, and is surrounded by a hot molecular shell. In the southeast, where the C II region expands into a dense, clumpy molecular cloud ridge, we see narrow and strong line emission from high-J CO lines, which comes from a thin, hot molecular shell surrounding the [CII] emission. The [CII] lines are broader and show photo evaporating gas flowing into the C II region. Based on the strength of the [13CII] F=2--1 line, the [CII] line appears to be somewhat optically thick over most of the nebula with an optical depth of a few. We model the physical conditions of the surrounding molecular cloud and the PDR emission using both RADEX and simple PDR models. The temperature of the CO emitting PDR shell is ~ 90 -- 120 K, with densities of 10^5 -- 10^6 cm^-3, as deduced from RADEX modeling. Our PDR modeling indicates that the PDR layer where [CII] emission dominates has somewhat lower densities, 10^4 to a few times 10^5 cm^-3

The formation of massive stars play a crucial role in galaxies from numerous points of view. The protostar generates a strong ultraviolet radiation field that ionizes its surroundings, and it drives powerful shock waves in the neighbouring medium in the form of jets and bipolar outflows, whose structure can be partially organized by local, strong magnetic field. Such an ejection activity locally modifies the interstellar chemistry, contributing to the cycle of matter. It also significantly participates to the energetic balance of galaxies. In the latter stages of massive star formation, the protostar is surrounded by an ultra-compact HII region, and irradiates its bipolar outflows, where an intrinsically strong magnetic field structure is associated to the generally high densities. In the HII region, or in the bipolar outflows, the question of in situ cosmic rays acceleration can then be raised by the simultaneous presence of strong magnetic fields, significant ionization of the matter, and mechanical energy available in large quantities. In this contribution, we will only summarise the results of our study of potential in situ cosmic rays acceleration the irradiated shocks in the W28 A2 massive star forming region, based on observations from the APEX, IRAM 30m, and Herschel telescopes, and presented in Gusdorf et al., recently submitted to A&A.

During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow regulates star formation by e.g. inducing turbulence and heating the surrounding gas. Understanding the associated shocks is therefore essential to the study of star formation. We present comparisons of shock models with CO, H2, and SiO observations in a 'pure' shock position in the BHR71 bipolar outflow. These comparisons provide an insight into the shock and pre-shock characteristics, and allow us to understand the energetic and chemical feedback of star formation on Galactic scales. New CO (Jup = 16, 11, 7, 6, 4, 3) observations from the shocked regions with the SOFIA and APEX telescopes are presented and combined with earlier H2 and SiO data (from the Spitzer and APEX telescopes). The integrated intensities are compared to a grid of models that were obtained from a magneto-hydrodynamical shock code which calculates the dynamical and chemical structure of these regions combined with a radiative transfer module based on the 'large velocity gradient' approximation. The CO emission leads us to update the conclusions of our previous shock analysis: pre-shock densities of 1e4 cm-3 and shock velocities around 20-25 km s-1 are still constrained, but older ages are inferred ( 4000 years). We evaluate the contribution of shocks to the excitation of CO around forming stars. The SiO observations are compatible with a scenario where less than 4% of the pre-shock SiO belongs to the grain mantles. We infer outflow parameters: a mass of 1.8x1e-2 Msun was measured in our beam, in which a momentum of 0.4 Msun km s-1 is dissipated, for an energy of 4.2x1e43erg. We analyse the energetics of the outflow species by species. Comparing our results with previous studies highlights their dependence on the method: H2 observations only are not sufficient to evaluate the mass of outflows.

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.

We report mapping observations of the $^{12}$CO $J=3-2$, $4-3$, $6-5$, and $7-6$ transitions and the CI 492 GHz transition toward the central 40$''\times$40$''$ region of the Circinus galaxy, using the Atacama Pathfinder EXperiment (APEX) telescope. We also detected $^{13}$CO $J=3-2$ at the central position of Circinus. These observations are to date the highest CO transitions reported in Circinus. With large velocity gradient (LVG) modeling and likelihood analysis we try to obtain density, temperature, and column density of the molecular gas in three regions: the nuclear region ($D< 18''\sim$ 360 pc), the entire central 45$''$ ($D<45''\sim$ 900pc) region, and the star-forming (S-F) ring (18$''<D<45''$). In the nuclear region, we can fit the CO excitation with a single excitation component, yielding an average condition of $n_{\rm H_2} \sim 10^{3.2}$cm$^{-3}$, $T_{\rm kin}$ $\sim$ 200 K, and d$v$/d$r \sim$ 3 km s$^{-1}$pc$^{-1}$. In the entire 45$''$ region, two excitation components are needed with $n_{\rm H_2}$ $\sim$ 10$^{4.2}$ and 10$^{3.0}$ cm$^{-3}$, $T_{\rm kin}\sim$ 60 K and 30 K, and $M_{\rm H_2}\sim$ $2.3\times 10^7$ M$_\odot$ and $6.6 \times 10^7$ M$_\odot$, respectively. The gas excitation in the S-F ring can also be fitted with two LVG components, after subtracting the CO fluxes in the 18$''$ region. The S-F ring region contributes 80\% of the molecular mass in the 45$''$ region. For the 45$''$ region, we find a conversion factor of $N({\rm H_2})/I_{\rm CO\ 1-0}$ = $0.37\times 10^{20}$ cm$^{-2} ({\rm K\ km\ s}^{-1})^{-1}$, about 1/5 of the Galactic disk value. The luminosity ratios of CI and $^{12}$CO $J=3-2$ ($R_{\rm CI/CO\ J=3-2}$) in Circinus basically follow a linear trend. The average $R_{\rm CI/CO\ J=3-2}$ in Circinus is $\sim$ 0.2, lying at an intermediate value between non-AGN nuclear region and high-redshift galaxies.

Supernovae constitute a critical source of energy input to the interstellar medium (ISM). In this short review, we focus on their latest phase of evolution, the supernova remnants (SNRs). We present observations of three old SNRs that have reached the phase where they interact with the ambient ISM: W28, IC443, and 3C391. We show that such objects make up clean laboratories to constrain the physical and chemical processes at work in molecular shock environments. Our studies subsequently allow us to quantify the impact of SNRs on their environment in terms of mass, momentum, and energy dissipation. In turn, their contribution to the energy balance of galaxies can be assessed. Their potential to trigger a further generation of star formation can also be investigated. Finally, our studies provide strong support for the interpretation of gamma-ray emission in SNRs, a crucial step to answer questions related to cosmic rays population and acceleration.

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 report HCN J=4-3, HCO+ J=4-3, and CS J=7-6 observations in 20 nearby star-forming galaxies with the Acatama Pathfinder EXperiment 12-m telescope. Combined with 4 HCN, 3 HCO+, and 4 CS detections in literature, we probe the empirical link between the luminosity of molecular gas (L_gas) and that of infrared emission (L_IR), up to the highest gas densities (10^6 - 10^8 cm-3) that have been probed so far. For nearby galaxies with large radii, we measure the IR luminosity within the submm beam-size (14"-18") to match the molecular emission. We find linear slopes for L_CS76-L_IR and L_HCN43-L_IR, and a slightly super-linear slope for L_HCO+43-L_IR. The correlation of L_CS76-L_IR even extends over eight orders of luminosity magnitude down to Galactic dense cores, with a fit of log(L_IR)=1.00(\pm 0.01) \times log(L_CS76) + 4.03(\pm 0.04). Such linear correlations appear to hold for all densities >10^4 cm-3, and indicate that star formation rate is not related to free-fall time scale for dense molecular gas.

Within a few parsecs around the central Black Hole Sgr A*, chemistry in the dense molecular cloud material of the circumnuclear disk (CND) can be affected by many energetic phenomena such as high UV-flux from the massive central star cluster, X-rays from Sgr A*, shock waves, and an enhanced cosmic-ray flux. Recently, spectroscopic surveys with the IRAM 30 meter and the APEX 12 meter telescopes of substantial parts of the 80--500 GHz frequency range were made toward selected positions in and near the CND. These datasets contain lines from the molecules HCN, HCO$^+$, HNC, CS, SO, SiO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$, H$_3$O$^+$ and others. We conduct Large Velocity Gradient analyses to obtain column densities and total hydrogen densities, $n$, for each species in molecular clouds located in the southwest lobe of CND. The data for the above mentioned molecules indicate 10$^5\,$cm$^{-3} \lesssim n <10^6\,$cm$^{-3}$, which shows that the CND is tidally unstable. The derived chemical composition is compared with a chemical model calculated using the UCL\_CHEM code that includes gas and grain reactions, and the effects of shock waves. Models are run for varying shock velocities, cosmic-ray ionization rates, and number densities. The resulting chemical composition is fitted best to an extremely high value of cosmic-ray ionization rate $\zeta \sim 10^{-14}\,$s$^{-1}$, 3 orders of magnitude higher than the value in regular Galactic molecular clouds, if the pre-shock density is $n=10^5\,$cm$^{-3}$.

Observations of galaxies at sub-millimeter wavelengths, where the emission is mainly due to cold dust, are required to constrain the dust physical properties and provide important insight on the gas content of galaxies. We mapped NGC1365 at 870um with LABOCA, the Large APEX Bolometer Camera, allowing us to probe the central mass concentration as well as the rate at which the gas flows to the center. We obtained the dust physical properties both globally and locally for different locations in the galaxy. A 20 K modified black body represents about 98% of the total dust content of the galaxy, the rest can be represented by a warmer dust component of 40 K. The bar exhibits an east-west asymmetry in the dust distribution: The eastern bar is heavier than the western bar by more than a factor of 4. Integrating the dust SED, we derive a total infrared (IR) luminosity of 9.8 x 10^{10} L_{\odot} leading to a dust-enshrouded star formation rate of ~16.7 M_{\odot} per year in NGC1365. We derive the gas mass from the measurements of the dust emission resulting in a CO-to-H2 conversion factor of X_{CO}~1.2x10^{20} mol cm^{-2} (K km s^{-1})^{-1} in the central disk including the bar. Taking into account the metallicity variation, the central gas mass concentration is only ~20% at R<40" (3.6 kpc). On the other hand, the time-scale with which the gas flows into the center, ~300 Myr, is rather short. This indicates that the current central mass in NGC1365 is evolving fast due to the strong bar.

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

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

We report the detection of strong and compact molecular line emission (in the CO J=3-2, 4-3, 6-5, 7-6, 13CO J=3-2, HCN and HCO+ J=4-3 transitions) from a cometary-shaped object (Carina-frEGG1) in the Carina star-forming region (SFR) previously classified as a photoevaporating protoplanetary disk (proplyd). We derive a molecular mass of 0.35 Msun for Carina-frEGG1, which shows that it is not a proplyd, but belongs to a class of free-floating evaporating gas globules (frEGGs) recently found in the Cygnus SFR by Sahai, Morris & Claussen (2012). Archival Adaptive Optics near-IR (Ks) images show a central hourglass-shaped nebula. The derived source luminosity (about 8-18 Lsun), the hourglass morphology, and the presence of collimated jets seen in HST images, imply the presence of a jet-driving, young, low-mass star deeply embedded in the dust inside Carina-frEGG1. Our results suggest that the true nature of many or most such cometary-shaped objects seen in massive SFRs and previously labelled as proplyds has been misunderstood, and that these are really frEGGs.

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.

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.

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 report the detection toward \eta\ Carinae of six new molecules, CO, CN, HCO+, HCN, HNC, and N2H+, and of two of their less abundant isotopic counterparts, 13CO and H13CN. The line profiles are moderately broad (about 100 km /s) indicating that the emission originates in the dense, possibly clumpy, central arcsecond of the Homunculus Nebula. Contrary to previous claims, CO and HCO+ do not appear to be under-abundant in \eta\ Carinae. On the other hand, molecules containing nitrogen or the 13C isotope of carbon are overabundant by about one order of magnitude. This demonstrates that, together with the dust responsible for the dimming of eta Carinae following the Great Eruption, the molecules detected here must have formed in situ out of CNO-processed stellar material.

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.

Atacama Pathfinder Experiment (APEX) mapping observations in CO (9-8) and (4-3) toward a high-mass star-forming region, NGC 6334 I, are presented. The CO (9-8) map has a 6.4 arcsec resolution, revealing a ~0.5 pc, jet-like, and bipolar outflow. This is the first map of a molecular outflow in a THz line. The CO (9-8) and (4-3) lines arising from the outflow lobes both show extremely high-velocity line wings, and their ratios indicate a gas temperature greater than 100 K and a density higher than 10^4 cm^{-3}. The spatial-velocity structure of the CO (9-8) data is typical of a bow-shock-driven flow, which is consistent with the association between the bipolar outflow and the infrared bow-shaped tips. In short, the observations unveil a highly-excited and collimated component in a bipolar outflow that is powered by a high-mass protostar, and provide insights into the driving mechanism of the outflow. Meanwhile, the observations demonstrate that high-quality mapping observations can be performed with the new THz receiver on APEX.

We have used the Atacama Pathfinder Experiment 12 m telescope (APEX) to carry out an absorption study of submillimeter wavelength rotational ground-state lines of H35Cl, H37Cl, 13CH+, and, for the first time, of the SH+ radical (sulfoniumylidene or sulfanylium). We detected the quartet of ground-state hyperfine structure lines of SH+ near 683 GHz with the CHAMP+ array receiver against the strong continuum source Sagittarius B2, which is located close to the center of our Galaxy. In addition to absorption from various kinematic components of Galactic center gas, we also see absorption at the radial velocities belonging to intervening spiral arms. This demonstrates that SH+ is a ubiquitous component of the diffuse interstellar medium. We do not find clear evidence for other SH+ lines we searched for, which is partially due to blending with lines from other molecules. In addition to SH+, we observed absorption from H35Cl, H37Cl, and 13CH+. The observed submillimeter absorption is compared in detail with absorption in 3 mm transitions of H13CO+ and c-C3H2 and the CO J = 1-0 and 3-2 transitions.

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.

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

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.

Current generation millimeter wavelength detectors suffer from scaling limits imposed by complex cryogenic readout electronics. To circumvent this it is imperative to investigate technologies that intrinsically incorporate strong multiplexing. One possible solution is the kinetic inductance detector (KID). In order to assess the potential of this nascent technology, a prototype instrument optimized for the 2 mm atmospheric window was constructed. Known as the N\'eel IRAM KIDs Array (NIKA), it was recently tested at the Institute for Millimetric Radio Astronomy (IRAM) 30-meter telescope at Pico Veleta, Spain. The measurement resulted in the imaging of a number of sources, including planets, quasars, and galaxies. The images for Mars, radio star MWC349, quasar 3C345, and galaxy M87 are presented. From these results, the optical NEP was calculated to be around $1 \times 10^{-15}$ W$ / $Hz$^{1/2}$. A factor of 10 improvement is expected to be readily feasible by improvements in the detector materials and reduction of performance-degrading spurious radiation.

The Atacama Pathfinder Experiment (APEX) 12m telescope was used to observe the N=1-0, J=0-1 ground state transitions of OH+ at 909.1588 GHz with the CHAMP+ heterodyne array receiver. Two blended hyperfine structure transitions were detected in absorption against the strong continuum source Sagittarius B2(M) and in several pixels offset by 18". Both, absorption from Galactic center gas as well as absorption from diffuse clouds in intervening spiral arms in a velocity range from -116 to 38.5 km/s is observed. The total OH+ column density of absorbing gas is 2.4 \times 10^15 cm-2. A column density local to Sgr B2(M) of 2.6 \times 10^14 cm-2 is found. On the intervening line-of-sight the column density per unit velocity interval are in the range from 1 to 40 \times 10^12 cm-2/(km/s). OH+ is found to be on average more abundant than other hydrides such as SH+ and CH+. Abundance ratios of OH and atomic oxygen to OH+ are found to be in the range of 10^1-2 and 10^3-4, respectively. The detected absorption of a continuous velocity range on the line-of-sight shows OH+ to be an abundant component of diffuse clouds.

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

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

Observations of dense molecular gas lie at the basis of our understanding of the density and temperature structure of protostellar envelopes and molecular outflows. We aim to characterize the properties of the protostellar envelope, molecular outflow and surrounding cloud, through observations of high excitation molecular lines within a sample of 16 southern sources presumed to be embedded YSOs. Observations of submillimeter lines of CO, HCO+ and their isotopologues, both single spectra and small maps were taken with the FLASH and APEX-2a instruments mounted on APEX to trace the gas around the sources. The HARP-B instrument on the JCMT was used to map IRAS 15398-3359 in these lines. HCO+ mapping probes the presence of dense centrally condensed gas, a characteristic of protostellar envelopes. The rare isotopologues C18O and H13CO+ are also included to determine the optical depth, column density, and source velocity. The combination of multiple CO transitions, such as 3-2, 4-3 and 7-6, allows to constrain outflow properties, in particular the temperature. Archival submillimeter continuum data are used to determine envelope masses. Eleven of the sixteen sources have associated warm and/or dense quiescent as characteristic of protostellar envelopes, or an associated outflow. Using the strength and degree of concentration of the HCO+ 4-3 and CO 4-3 lines as a diagnostic, five sources classified as Class I based on their spectral energy distributions are found not to be embedded YSOs. The C18O 3-2 lines show that for none of the sources, foreground cloud layers are present. Strong molecular outflows are found around six sources, .. (continued in paper)

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

We present the first maps of high-excitation CO J=6-5 and 7-6 and isotopologue lines over 2'-5' regions at 10" resolution toward low-mass protostars to probe the origin of the warm gas in their surroundings. The data were obtained using the CHAMP+ 650/850 GHz heterodyne array receiver on APEX. Surprisingly strong quiescent extended narrow-line high-J 12CO 6-5 and 7-6 emission is seen toward all protostars, suggesting that heating by UV photons along the outflow cavity dominates the emission. At the source position itself, passive heating of the collapsing inner envelope by the luminosity of the source also contributes. The UV photons are generally not energetic enough to dissociate CO since the [C I] 2-1 emission, also probed by our data, is weak except at the bow shock at the tip of the outflow. The extended UV radiation is produced by the star-disk boundary layer as well as the jet- and bow-shocks, and will also affect the chemistry of species such as H2O and HCN around the outflow axis. Shock-heated warm gas characterized by broad CO line profiles is seen only toward the more massive Class 0 outflows. Outflow temperatures, estimated from the CO 6-5/3-2 line wing ratios, are ~100 K. These data illustrate the importance of getting spatial information to characterize the physical processes in YSO surroundings, important to interpret future Herschel and ALMA data.

Molecular outflows from high-mass young stellar objects provide an excellent way to study the star formation process, and investigate if they are scaled-up versions of their low-mass counterparts. We selected the nearby massive star forming region IRAS 17233-3606 in order to study the kinematics and physics along the molecular outflow(s) originating from this source. We observed IRAS 17233-3606 in CO, a typical tracer of gas associated with molecular outflow, with the Submillimeter Array in the (2-1) transition, and with the APEX telescope in the higher excitation (6-5) line. Additional infrared H2 observations were performed with the UKIRT telescope. The CO data were analysed using a LVG approach. Our data resolve the previously detected molecular outflow in at least three different components, one of them with a high collimation factor ~4, and characterised by emission at extremely high velocities (|v-v_{LSR}|>120 km s^{-1}). The estimate of the kinematical outflow parameters are typical of massive YSOs, and in agreement with the measured bolometric luminosity of the source. The kinematic ages of the flows are in the range 10^2-10^3 yr, and therefore point to young objects that still did not reach the main sequence.

The origin and heating mechanisms of warm (50<T<200 K) molecular gas in low-mass young stellar objects (YSOs) are strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by shocks and by UV associated with the outflows or accretion have been proposed. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, bipolar outflows and the quiescent cloud. High-J CO maps (12CO J=6--5 and 7--6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of the 13CO J=6--5 transition and [C I] 3P_2-3P_1 line were taken. Strong quiescent narrow-line 12CO 6--5 and 7--6 emission is seen toward all protostars. In the case of HH~46 and Ced 110 IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. Shock-heated warm gas is only detected for Class 0 flows and the more massive Class I sources such as HH~46. Outflow temperatures, estimated from the CO 6--5 and 3--2 line wings, are ~100 K, close to model predictions, with the exception of the L~1551 IRS 5 and IRAS 12496-7650, for which temperatures <50 K are found. APEX-CHAMP+ is uniquely suited to directly probe a protostar's feedback on its accreting envelope gas in terms of heating, photodissociation, and outflow dispersal by mapping 1'x1' regions in high-J CO and [C I] lines.

We present observations of the Sunyaev-Zel'dovich effect (SZE) in the Bullet cluster (1E 0657--56) using the APEX-SZ instrument at 150 GHz with a resolution of 1 arcmin. The main results are maps of the SZE in this massive, merging galaxy cluster. The cluster is detected with 23 sigma significance within the central 1 arcmin radius of the source position. The SZE map has a broadly similar morphology to that in existing X-ray maps of this system, and we find no evidence for significant contamination of the SZE emission by radio or IR sources. In order to make simple quantitative comparisons with cluster gas models derived from X-ray observations, we fit our data to an isothermal elliptical beta model, despite the inadequacy of such a model for this complex merging system. With an X-ray derived prior on the power-law index, beta = 1.04 +0.16 -0.10, we find a core radius r_c =142 +/- 18 arcsec, an axial ratio of 0.889 +/- 0.072, and a central temperature decrement of -771 +/- 71 micro-K_CMB, including a +/-5.5% flux calibration uncertainty. Combining the APEX-SZ map with a map of projected electron surface density from Chandra X-ray observations, we determine the mass-weighted temperature of the cluster gas to be T_mg=10.8 +/- 0.9 keV, significantly lower than some previously reported X-ray spectroscopic temperatures. Under the assumption of an isothermal cluster gas distribution in hydrostatic equilibrium, we compute the gas mass fraction for prolate and oblate spheroidal geometries and find it to be consistent with previous results from X-ray and weak lensing observations. This work is the first result from the APEX-SZ experiment, and represents the first reported scientific result from observations with a large array of multiplexed superconducting transition-edge sensor bolometers.

AIMS Our aim is to characterize the size, mass, density and temperature profiles of the protostellar envelope of HH~46 IRS 1 and its surrounding cloud material as well as the effect the outflow has on its environment.METHODS The CHAMP+ and LABOCA arrays on the APEX telescope, combined with lower frequency line receivers, are used to obtain a large continuum map and smaller heterodyne maps in various isotopologues of CO and HCO+. The high-J lines of CO (6--5 and 7--6) and its isotopologues together with [C I] 2--1, observed with CHAMP+, are used to probe the warm molecular gas in the inner few hundred AU and in the outflowing gas. The data are interpreted with continuum and line radiative transfer models. RESULTS Broad outflow wings are seen in CO low- and high-J lines at several positions, constraining the gas temperatures to a constant value of ~100 K along the red outflow axis and to ~60 K for the blue outflow. The derived outflow mass is of order 0.4--0.8 M_sol, significantly higher than previously found. The bulk of the strong high-J CO line emission has a surprisingly narrow width, however, even at outflow positions. These lines cannot be fit by a passively heated model of the HH 46 IRS envelope. We propose that it originates from photon heating of the outflow cavity walls by ultraviolet photons originating in outflow shocks and the accretion disk boundary layers. At the position of the bow shock itself, the UV photons are energetic enough to dissociate CO. The envelope mass of ~5 M_sol is strongly concentrated towards HH 46 IRS with a density power law of -1.8.

We present CO(2-1), $^{13}$CO(2-1), CO(6-5), CO(7-6), and SO(6$_5-5_4$) line observations made with the {\it IRAM 30 m} and {\it APEX} radiotelescopes and the {\it Submillimeter Array} toward the highly collimated and extended southwest lobe of the bipolar outflow {\it Ori-S6} located in the Orion South region. We report, for all these lines, the detection of velocity asymmetries about the flow axis, with velocity differences roughly on the order of 1 km s$^{-1}$ over distances of about 5000 AU, 4 km s$^{-1}$ over distances of about 2000 AU, and close to the source of between 7 and 11 km s$^{-1}$ over smaller scales of about 1000 AU. We interpret these velocity differences as a signature of rotation but also discuss some alternatives which we recognize as unlikely in view of the asymmetries' large downstream continuation. This rotation across the {\it Ori-S6} outflow is observed out to (projected) distances beyond 2.5 $\times$ 10$^4$ AU from the flow's presumed origin. Comparison of our large-scale and small-scale observations suggests the rotational velocity to decline not faster than 1/R with distance R from the axis; in the innermost few arcsecs an increase of rotational velocity with R is even indicated. The magnetic field lines threading the inner rotating CO shell may well be anchored in a disk of radius $\sim$ 50 AU; the field lines further out need a more extended rotating base.

The Large APEX Bolometer Camera, LABOCA, has been commissioned for operation as a new facility instrument t the Atacama Pathfinder Experiment 12m submillimeter telescope. This new 295-bolometer total power camera, operating in the 870 micron atmospheric window, combined with the high efficiency of APEX and the excellent atmospheric transmission at the site, offers unprecedented capability in mapping submillimeter continuum emission for a wide range of astronomical purposes.

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.

Using the Atacama Pathfinder Experiment (APEX) telescope we have detected the rotational ground-state transitions of ortho-ammonia and ortho-water toward the redshift ~0.89 absorbing galaxy in the PKS 1830-211 gravitational lens system. We discuss our observations in the context of recent space-borne data obtained for these lines with the SWAS and Odin satellites toward Galactic sources. We find commonalities, but also significant differences between the interstellar media in a galaxy at intermediate redshift and in the Milky Way. Future high-quality observations of the ground-state ammonia transition in PKS 1830-211, together with inversion line data, will lead to strong constraints on the variation of the proton to electron mass ratio over the past 7.2 Gyr.

We present large scale 870 micron maps of the nearby starburst galaxies NGC253, NGC4945 and the nearest giant elliptical radio galaxy Centaurus A (NGC 5128) obtained with the newly commissioned Large Apex Bolometer Camera (LABOCA) operated at the APEX telescope. Our continuum images reveal for the first time the distribution of cold dust at a angular resolution of 20" across the entire optical disks of NGC253 and NGC4945 out to a radial distance of 10' (7.5 kpc). In NGC5128 our LABOCA image also shows, for the first time at submillimeter wavelengths, the synchrotron emission associated with the radio jet and the inner radio lobes. From an analysis of the 870 micron emission in conjunction with ISO-LWS, IRAS and long wavelengths radio data we find temperatures for the cold dust in the disks of all three galaxies of 17-20 K, comparable to the dust temperatures in the disk of the Milky Way. The total gas mass in the three galaxies is determined to be 2.1, 4.2 and 2.8 x 10^9 solar masses for NGC253, NGC4945 and NGC5128, respectively. A detailed comparison between the gas masses derived from the dust continuum and the integrated CO(1-0) intensity in NGC253 suggests that changes of the CO luminosity to molecular mass conversion factor are mainly driven by a metallicity gradient and only to a lesser degree by variations of the CO excitation. An analysis of the synchrotron spectrum in the northern radio lobe of NGC5128 shows that the synchrotron emission from radio to the ultraviolet wavelengths is well described by a broken power law and that the break frequency is a function of the distance from the radio core as expected for aging electrons. We derive an outflow speed of ~0.5c at a distance of 2.6kpc from the center, consistent with the speed derived in the vicinity of the nucleus.

We present molecular line observations of the massive star forming region IRAS 17233-3606 aimed at studying the molecular core associated with the source. The observations were made using the Atacama Pathfinder Experiment telescope in the CO (3-2) and HCO^+ (4-3) transitions, and in the CH_3OH (6_K-5_K), (7_K-6_K) and CH_3CN (16_K-15_K) bands. For the CO(3-2) and HCO^+ (4-3) transitions, we obtained maps with a size of 70''\times 70''. The typical angular resolution of the data is ~18''. Our observations reveal an exceptionally rich molecular spectrum, a signpost of hot core activity. Comparisons with two other prominent southern hot cores were made through observations in the same frequency setups. We also detected a bipolar outflow in CO (3-2) and HCO^+ (4-3) lines. Modelling reveals a hot core of size ~3'' and a temperature of 150 K in the IRAS17233-3606 region. The parameters of the molecular outflow are derived through the analysis of the CO (3-2) emission, and are typical of outflows driven by high-mass young stellar objects.

We observed the high velocity outflow originating from NGC6334I in several CO transitions with the APEX telescope, with the goal of deriving the physical parameters of the gas. Using an LVG analysis, we studied line ratios between the CO(3-2), CO(4-3), and CO(7-6) data as a function of the density and of the kinetic temperature of the gas. An upper limit on the CO column density is derived by comparison with 13CO data. We constrained the temperature to be higher than 50 K and the H_2 density to values higher than n~10^4 cm^{-3} towards the peak position in the red lobe, while T>15 K and n>10^3 cm^{-3} are derived towards the peak position in the blue lobe. The physical parameters of the outflow, its mass and its energetics, have been computed using the temperatures derived from this analysis. We conclude that high kinetic temperatures are present in the outflow and traced by high excitation CO lines. Observations of high-J CO lines are thus needed to infer reliable values of the kinetic temperatures and of the other physical parameters in outflows.

Context: The M16 nebula is a relatively nearby Hii region, powered by O stars from the open cluster NGC 6611, which borders to a Giant Molecular Cloud. Radiation from these hot stars has sculpted columns of dense obscuring material on a few arcmin scales. The interface between these pillars and the hot ionised medium provides a textbook example of a Photodissociation Region (PDR). Aims: To constrain the physical conditions of the atomic and molecular material with submillimeter spectroscopic observations. Methods: We used the APEX submillimeter telescope to map a ~3'x3' region in the CO J=3-2, 4-3 and 7-6 rotational lines, and a subregion in atomic carbon lines. We also observed C18O(3-2) and CO(7-6) with longer integrations on five peaks found in the CO(3-2) map. The large scale structure of the pillars is derived from the molecular lines' emission distribution. We estimate the magnitude of the velocity gradient at the tips of the pillars and use LVG modelling to constrain their densities and temperatures. Excitation temperatures and carbon column densities are derived from the atomic carbon lines. Results: The atomic carbon lines are optically thin and excitation temperatures are of order 60 K to 100 K, well consistent with observations of other Hii region-molecular cloud interfaces. We derive somewhat lower temperatures from the CO line ratios, of order 40 K. The Ci/CO ratio is around 0.1 at the fingers tips.

Aims. The purpose of this study is to examine the prediction that the deuterated H3+ ion, H2D+, can be found exclusively in the coldest regions of molecular cloud cores. This is also a feasibility study for the detection of the ground-state line of ortho-H2D+ at 372 GHz with APEX. Methods. The 1(10)-1(11) transition of H2D+ at 372 GHz was searched towards selected positions in the massive star forming cloud OriB9, in the dark cloud L183, and in the low- to intermediate mass star-forming cloud R CrA. Results. The line was detected in cold, prestellar cores in the regions of OriB9 and L183, but only upper limits were obtained towards other locations which either have elevated temperatures or contain a newly born star. The H2D+ detection towards OriB9 is the first one in a massive star-forming region. The fractional ortho-H2D+ abundances (relative to H2) are estimated to be about 1.0E-10 in two cold cores in OriB9, and 3.0E-10 in the cold core of L183. Conclusions. The H2D+ detection in OriB9 shows that also massive star forming regions contain very cold prestellar cores which probably have reached matured chemical composition characterized, e.g., by a high degree of deuterium fractionation. Besides as a tracer of the interior parts of prestellar cores, H2D+ may therefore be used to put contraints on the timescales related to massive star formation.

Vibrationally excited emission from the SiO and H2O molecules probes the innermost circumstellar envelopes of oxygen-rich red giant and supergiant stars. VY CMa is the most prolific known stellar emission source in these molecules. Observations were made to search for rotational lines in the lowest vibrationally excited state of H2O. The APEX telescope was used for observations of H2O lines at frequencies around 300 GHz. Two vibrationally excited H2O lines were detected, a third one could not be found. In one of the lines we find evidence for weak maser action, similar to known (sub)millimeter H2O lines. We find that the other line's intensity is consistent with thermal excitation by the circumstellar infrared radiation field. Several SiO lines were detected together with the H2O lines.

We present a submm study of the massive hot core G327.3-0.6 that constrains its physical parameters and environment. The APEX telescope was used to image CO and N2H+ emission, to observe lines from other molecules toward a hot and a cold molecular core, and to measure the continuum flux density of the hot core. In the C18O J=3-2 line, two clumps were found, one associated with the HII region G327.3-0.5 and the other associated with the hot core. An additional cold clump is found 30 arcsec (0.4 pc) northeast of the hot core in bright N2H+ emission. From the the continuum data, we calculate a mass of 420 Msol and a size of 0.1 pc for the hot core. A new, more accurate position of the hot core is reported, which allows the association of the core with a bright mid-infrared source. The luminosity of the hot core is estimated to be between 5 and 15 10^4 Lsol. This study revealed several different evolutionary stages of massive star formation in the G327.3-0.6 region.

A study of 12 ultracompact HII regions was conducted to probe the physical conditions and kinematics in the inner envelopes of the molecular clumps harboring them. The APEX telescope was used to observe the sources in the CO (4-3) and 13CO (8-7) lines. Line intensities were modeled with the RATRAN radiative transfer code using power laws for the density and temperature to describe the physical structure of the clumps. All sources were detected in both lines. The optically thick CO (4-3) line shows predominantly blue skewed profiles reminiscent of infall. Line intensities can be reproduced well using the physical structure of the clumps taken from the literature. The optically thick line profiles show that CO is a sensitive tracer of ongoing infall in the outer envelopes of clumps harboring ultracompact HII regions and hot molecular cores.

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

Aims:To obtain insight into the physical structure of the warm gas in the inner envelope of protostars and the interaction with the outflow. Methods: Sub-millimeter observations of 12CO, 13CO and/or C18O in J=3-2, J=4-3 and J=7-6 were obtained with the APEX Telescope towards IRAS 12496-7650, an intermediate mass young stellar object. The data are compared to ISO-LWS observations of CO J=14-13 up to J=19-18 lines to test the different proposed origins of the CO lines. Results: The outflow is prominently detected in the 3-2 and 4-3 lines, but not seen at similar velocities in the 7-6 line, constraining the temperature in the high-velocity (> 5 km/s from line center) gas to less than 50 K, much lower than inferred from the analysis of the ISO-LWS data. In addition, no isothermal gas model can reproduce the emission in both the 7-6 and the higher-J ISO-LWS lines.The 7-6 line probably originates in the inner (<250 AU) region of the envelope at ~150 K. Detailed radiative transfer calculations suggest that the ISO-LWS lines are excited by a different mechanism, possibly related to the larger-scale outflow. All possible mechanisms on scales smaller than 8 " are excluded. High-resolution continuum as well as high-J 12CO and isotopic line mapping are needed to better constrain the structure of the warm gas in the inner envelope and the interaction with the outflow.

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

Context : Highly-collimated outflows are believed to be the earliest stage in outflow evolution, so their study is essential for understanding the processes driving outflows. The BHR71 Bok globule is known to harbour such a highly-collimated outflow, which is powered by a protostar belonging to a protobinary system. Aims : We aimed at investigating the interaction of collimated outflows with the ambient molecular cloud by using molecular tracers. Methods : We mapped the BHR71 highly-collimated outflow in CO(3-2) with the APEX telescope, and observed several bright points of the outflow in the molecular transitions CO(4-3), 13CO(3-2), C18O(3-2), and CH3OH(7-6). We use an LVG code to characterise the temperature enhancements in these regions. Results : In our CO(3-2) map, the second outflow driven by IRS2, which is the second source of the binary system, is completely revealed and shown to be bipolar. We also measure temperature enhancements in the lobes. The CO and methanol LVG modelling points to temperatures between 30 and 50K in the IRS1 outflow, while the IRS2 outflow seems to be warmer (up to 300K).

At a distance of 8 kpc only, the center of the Galaxy provides a unique testbed for studies about the physics of the interstellar medium and star formation in the nuclei of galaxies. Since it shares many properties of other, more spectacular nuclei, a detailed investigation of the Galactic Center is a necessary step to gain better understanding of the physical processes governing galactic nuclei in general. We review the characteristics of the central molecular cloud layer, the physical properties of these clouds, and their interplay with the magnetic field.