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

We present the performance of the upGREAT heterodyne array receivers on the SOFIA telescope after several years of operations. This instrument is a multi-pixel high resolution (R > 10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receivers use 7-pixel subarrays configured in a hexagonal layout around a central pixel. The low frequency array receiver (LFA) has 2x7 pixels (dual polarization), and presently covers the 1.83-2.06 THz frequency range, which allows to observe the [CII] and [OI] lines at 158 um and 145 um wavelengths. The high frequency array (HFA) covers the [OI] line at 63 um and is equipped with one polarization at the moment (7 pixels, which can be upgraded in the near future with a second polarization array). The 4.7 THz array has successfully flown using two separate quantum-cascade laser local oscillators from two different groups. NASA completed the development, integration and testing of a dual-channel closed-cycle cryocooler system, with two independently operable He compressors, aboard SOFIA in early 2017 and since then, both arrays can be operated in parallel using a frequency separating dichroic mirror. This configuration is now the prime GREAT configuration and has been added to SOFIA's instrument suite since observing cycle 6.

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

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

We study the effects of an asymmetric radiation field on the properties of a molecular cloud envelope. We employ observations of carbon monoxide (12CO and 13CO), atomic carbon, ionized carbon, and atomic hydrogen to analyze the chemical and physical properties of the core and envelope of L1599B, a molecular cloud forming a portion of the ring at approximately 27 pc from the star Lambda Ori. The O III star provides an asymmetric radiation field that produces a moderate enhancement of the external radiation field. Observations of the [CII] fine structure line with the GREAT instrument on SOFIA indicate a significant enhanced emission on the side of the cloud facing the star, while the [Ci], 12CO and 13CO J = 1-0 and 2-1, and 12CO J = 3-2 data from the PMO and APEX telescopes suggest a relatively typical cloud interior. The atomic, ionic, and molecular line centroid velocities track each other very closely, and indicate that the cloud may be undergoing differential radial motion. The HI data from the Arecibo GALFA survey and the SOFIA/GREAT [CII] data do not suggest any systematic motion of the halo gas, relative to the dense central portion of the cloud traced by 12CO and 13CO.

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.

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

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.

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

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

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

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