The new field of gravitational wave astrophysics requires a growing pool of students and researchers with unique, interdisciplinary skill sets. It also offers an opportunity to build a diverse, inclusive astronomy community from the ground up. We describe the efforts used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) NSF Physics Frontiers Center to foster such growth by involving students at all levels in low-frequency gravitational wave astrophysics with pulsar timing arrays (PTAs) and establishing collaboration policies that ensure broad participation by diverse groups. We describe and illustrate the impact of these techniques on our collaboration as a case study for other distributed collaborations.

We have fabricated planar amorphous Indium Oxide superconducting resonators ($T_c\sim2.8$ K) that are sensitive to frequency-selective radiation in the range of 7 to 10 GHz. Those values lay far below twice the superconducting gap that worths about 200 GHz. The photons detection consists in a shift of the fundamental resonance frequency. We show that the detected frequency can be adjusted by modulating the total length of the superconducting resonator. We attribute those observations to the excitation of higher-order resonance modes. The coupling between the fundamental lumped and the higher order distributed resonance is due to the kinetic inductance non-linearity with current. These devices, that we have called Sub-gap Kinetic Inductance Detectors (SKIDs), are to be distinguished from the standard Kinetic Inductance Detectors (KIDs) in which quasi-particles are generated when incident light breaks down Cooper pairs.

Context: We describe the new SEPIA (Swedish-ESO PI Instrument for APEX) receiver, which was designed and built by the Group for Advanced Receiver Development (GARD), at Onsala Space Observatory (OSO) in collaboration with ESO. It was installed and commissioned at the APEX telescope during 2015 with an ALMA Band 5 receiver channel and updated with a new frequency channel (ALMA Band 9) in February 2016. Aims: This manuscript aims to provide, for observers who use the SEPIA receiver, a reference in terms of the hardware description, optics and performance as well as the commissioning results. Methods: Out of three available receiver cartridge positions in SEPIA, the two current frequency channels, corresponding to ALMA Band 5, the RF band 158--211 GHz, and Band 9, the RF band 600--722 GHz, provide state-of-the-art dual polarization receivers. The Band 5 frequency channel uses 2SB SIS mixers with an average SSB noise temperature around 45K with IF (intermediate frequency) band 4--8 GHz for each sideband providing total 4x4 GHz IF band. The Band 9 frequency channel uses DSB SIS mixers with a noise temperature of 75--125K with IF band 4--12 GHz for each polarization. Results: Both current SEPIA receiver channels are available to all APEX observers.

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

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

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

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

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

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

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

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

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

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

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

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

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