The ground-breaking image of a black hole's event horizon, which captured the public's attention and imagination in April 2019, was captured using the power of interferometry: many separate telescopes working together to observe the cosmos in incredible detail. Many recent astrophysical discoveries that have revolutionized the scientific community's understanding of the cosmos were made by interferometers such as LIGO, ALMA, and the Event Horizon Telescope. Astro 101 instructors who want their students to learn the science behind these discoveries must teach about interferometry. Decades of research show that using active learning strategies can significantly increase students' learning and reduces achievement gaps between different demographic groups over what is achieved from traditional lecture-based instruction. As part of an effort to create active learning materials on interferometry, we developed and tested a new Lecture-Tutorial to help Astro 101 students learn about key properties of astronomical interferometers. This paper describes this new Lecture-Tutorial and presents evidence for its effectiveness from a study conducted with 266 Astro 101 students at the University of North Carolina at Chapel Hill.

The current generation of radio and millimeter telescopes, particularly the Atacama Large Millimeter Array (ALMA), offers enormous advances in observing capabilities. While these advances represent an unprecedented opportunity to facilitate scientific understanding, the increased complexity in the spatial and spectral structure of these ALMA data cubes lead to challenges in their interpretation. In this paper, we perform a feasibility study for applying topological data analysis and visualization techniques never before tested by the ALMA community. Through techniques based on contour trees, we seek to improve upon existing analysis and visualization workflows of ALMA data cubes, in terms of accuracy and speed in feature extraction. We review our application development process in building effective analysis and visualization capabilities for the astrophysicists. We also summarize effective design practices by identifying domain-specific needs of simplicity, integrability, and reproducibility, in order to best target and service the large astrophysics community.

We present $\sim$0.1" resolution ($\sim$10 pc) ALMA observations of a molecular cloud identified in the merging Antennae galaxies with the potential to form a globular cluster, nicknamed the ``Firecracker.' Since star formation has not yet begun at an appreciable level, this cloud provides an example of what the birth environment of a globular cluster may have looked like before stars form and disrupt the natal conditions. Using emission from $^{12}$CO(2-1), $^{12}$CO(3-2), $^{13}$CO(2-1), HCN(4-3), and HCO$^+$(4-3), we are able to resolve the cloud's structure and find that it has a characteristic radius of 22 pc and a mass of 1--9$\times10^6 M_\odot$. We also constrain the abundance ratios of $^{12}$CO/$^{13}$CO and H$_2$/\twelveCO. Based on the calculated mass, we determine that the commonly used CO-to-H$_2$ conversion factor varies spatially, with average values in the range $X_{CO}=(0.12-1.1)\times10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We demonstrate that if the cloud is bound (as is circumstantially suggested by its bright, compact morphology), an external pressure of $P/k > 10^8$ K cm$^{-3}$ is required. This would be consistent with theoretical expectations that globular cluster formation requires high pressure environments. The position-velocity diagram of the cloud and its surrounding material suggests that this pressure may be produced by the collision of filaments. The radial profile of the column density can be fit with both a Gaussian and Bonnor-Ebert profile. The relative line strengths of HCN and HCO$^+$ in this region suggest that these molecular lines can be used as tracers for the evolutionary stage of a cluster.

Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.

Modeling of the spectral line energy distribution (SLED) of the CO molecule can reveal the physical conditions (temperature, density) of molecular gas in Galactic clouds and other galaxies. Recently, the Herschel Space Observatory and ALMA have offered, for the first time, a comprehensive view of the rotational J = 4-3 through J = 13-12 lines, which arise from a complex, diverse range of physical conditions that must be simplified to one, two, or three components when modeled. Here we investigate the recoverability of physical conditions from SLEDs produced by galaxy evolution simulations containing a large dynamical range in physical properties. These simulated SLEDs were generally fit well by one component of gas whose properties largely resemble or slightly underestimate the luminosity-weighted properties of the simulations when clumping due to non-thermal velocity dispersion is taken into account. If only modeling the first three rotational lines, the median values of the marginalized parameter distributions better represent the luminosity-weighted properties of the simulations, but the uncertainties in the fitted parameters are nearly an order of magnitude, compared to approximately 0.2 dex in the "best-case" scenario of a fully sampled SLED through J = 10-9. This study demonstrates that while common CO SLED modeling techniques cannot reveal the underlying complexities of the molecular gas, they can distinguish bulk luminosity-weighted properties that vary with star formation surface densities and galaxy evolution, if a sufficient number of lines are detected and modeled.

We present Herschel SPIRE-FTS intermediate-sampled mapping observations of the central ~8 kpc (~150") of M51, with a spatial resolution of 40". We detect 4 12CO transitions (J=4-3 to J=7-6) and the CI 3P2-3P1 and 3P1-3P0 transitions. We supplement these observations with ground based observations of 12CO J=1-0 to J=3-2 and perform a two-component non-LTE analysis. We find that the molecular gas in the nucleus and centre regions has a cool component (T_kin~10-20 K) with a moderate but poorly constrained density (n(H2)~10^3-10^6 cm^-3), as well as significant molecular gas in a warmer (T_kin~300-3000 K), lower density (n(H2)~10^1.6-10^2.5 cm^-3) component. We compare our CO line ratios and calculated densities along with ratios of CO to total infrared luminosity to a grid of photon dominated region (PDR) models and find that the cold molecular gas likely resides in PDRs with a field strength of G_o~10^2. The warm component likely requires an additional source of mechanical heating, from supernovae and stellar winds or possibly shocks produced in the strong spiral density wave. When compared to similar two-component models of other star-forming galaxies published as part of the Very Nearby Galaxies Survey (Arp 220, M82 and NGC 4038/39), M51 has the lowest density for the warm component, while having a warm gas mass fraction that is comparable to those of Arp 220 and M82, and significantly higher than that of NGC 4038/39.

We have conducted two-component, non-LTE modeling of the CO lines from J = 1-0 through J = 13-12 in 87 galaxies observed by the Herschel SPIRE Fourier Transform Spectrometer (FTS). We find the average pressure of the cold molecular gas, traced especially by CO J = 1-0, is $\sim 10^{5.0 \pm 0.5}$ K cm$^{-3}$. The mid- to high-J lines of CO trace higher-pressure gas at $10^{6.5 \pm 0.6}$ K cm$^{-3}$; this pressure is slightly correlated with LFIR. Two components are often necessary to accurately fit the Spectral Line Energy Distributions (SLEDs); a one-component fit often underestimates the flux of CO J = 1-0 and the mass. If low-J lines are not included, mass is underestimated by an order of magnitude. Even when modeling the low-J lines alone or using an $\alpha_{CO}$ conversion factor, the mass should be considered to be uncertain to a factor of at least 0.4 dex, and the vast majority of the CO luminosity will be missed (median, 65%). We find a very large spread in our derived values of $\alpha_{CO}$, though they do not have a discernible trend with LFIR; the best fit is a constant 0.7 M$_{\odot}$/ (K kms$^{-1}$ pc$^2$), with a standard deviation of 0.36 dex, and a range of 0.3-1.6 M$_{\odot}$/ (K kms$^{-1}$ pc$^2$). We find average molecular gas depletion times ($\tau_{dep}$) of $10^8$ yr that decrease with increasing SFR. Finally, we note that the J = 11-10/J = 1-0 line flux ratio is diagnostic of the warm component pressure, and discuss the implications of this comprehensive study of SPIRE FTS extragalactic spectra for future study post-Herschel.

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

We present a catalog of all CO (J=4-3 through J=13-12)), [CI], [NII] lines available from extragalactic spectra from the Herschel SPIRE Fourier Transform Spectrometer (FTS) archive combined with observations of the low-J CO lines from the literature and from the Arizona Radio Observatory. This work examines the relationships between LFIR, L'CO, and LCO/LCO(1-0). We also present a new method for estimating probability distribution functions (PDFs) from marginal signal-to-noise ratio Herschel} FTS spectra, which takes into account the instrumental "ringing" and the resulting highly correlated nature of the spectra. The slopes of log(LFIR) vs. log(L'CO) are linear for all mid- to high-J CO lines and slightly sublinear if restricted to (U)LIRGs. The mid- to high-J CO luminosity relative to CO J=1-0 increases with increasing LFIR, indicating higher excitement of the molecular gas, though these ratios do not exceed ~ 180. For a given bin in LFIR, the luminosities relative to CO J=1-0 remain relatively flat from J=6-5 through J=13-12, across three orders of magnitude of LFIR. A single component theoretical photon-dominated region (PDR) model cannot match these flat SLED shapes, though combinations of PDR models with mechanical heating added qualitatively match the shapes, indicating the need for further comprehensive modeling of the excitation processes of warm molecular gas in nearby galaxies.

We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3mm spectral scans between 84-114GHz for 15 galaxies and targeted ALMA 1mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [CI] , [NII] , H_2O and NH_3. We further present APEX [CII] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new mm/submm line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z=3.9+/-0.4, and the highest-redshift source in our sample is at z=5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4mm-selected sources with a median redshift of z=3.1+/-0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.

We present Atacama Large Millimeter Array (ALMA) Cycle-0 observations of the CO J = 6-5 line in the advanced galaxy merger Arp 220. This line traces warm molecular gas, which dominates the total CO luminosity. The CO emission from the two nuclei is well resolved by the 0.39" x 0.22" beam and the exceptional sensitivity and spatial/spectral resolution reveal new complex features in the morphology and kinematics of the warm gas. The line profiles are asymmetric between the red and blue sides of the nuclear disks and the peak of the line emission is offset from the peak of the continuum emission in both nuclei by about 100 pc in the same direction. CO self-absorption is detected at the centers of both nuclei but it is much deeper in the eastern nucleus. We also clearly detect strong, highly redshifted CO absorption located near the southwest side of each nucleus. For the eastern nucleus, we reproduce the major line profile features with a simple kinematic model of a highly turbulent, rotating disk with a substantial line center optical depth and a large gradient in the excitation temperature. The red/blue asymmetries and line-to-continuum offset are likely produced by absorption of the blue (SW) sides of the two nuclei by blue-shifted, foreground molecular gas; the mass of the absorber is comparable to the nuclear warm gas mass (10^8 M_solar). We measure an unusually high L_CO/L_FIR ratio in the eastern nucleus, suggesting there is an additional energy source, such as mechanical energy from shocks, present in this nucleus.

With a substantial nuclear molecular gas reservoir and broad, high-velocity CO molecular line wings previously interpreted as an outflow, NGC 1266 is a rare SB$0$ galaxy. Previous analyses of interferometry, spectrally resolved low-$J$ CO emission lines, and unresolved high-$J$ emission lines have established basic properties of the molecular gas and the likely presence of an AGN. Here, new spectrally resolved CO $J = 5 - 4$ to $J = 8 - 7$ lines from {\it Herschel Space Observatory} HIFI observations are combined with ground-based observations and high-$J$ {\it Herschel} SPIRE observations to decompose the nuclear and putative outflow velocity components and to model the molecular gas to quantify its properties. Details of the modeling and results are described, with comparisons to previous results and exploration of the implications for the gas excitation mechanisms. Among the findings, like for other galaxies, the nuclear and putative outflow molecular gas are well represented by components that are cool ($T_{nuclear} = 6^{+10}_{-2}$ K and $T_{outflow} \sim 30$ K), comprising bulk of the mass (Log $M_{nuclear}/M_{\odot} = 8.3^{+0.5}_{-0.4}$ and Log $M_{outflow}/M_{\odot} = 7.6^{+0.3}_{-0.3}$), and the minority of the luminosity (Log $L_{nuclear}/L_{\odot} = 5.44^{+0.22}_{-0.18}$ and Log $L_{outflow}/L_{\odot} \sim 6.5$) and warm ($T_{nuclear} = 74^{+130}_{-26}$ K and $T_{outflow} > 100$ K), comprising a minority of the mass (Log $M_{nuclear}/M_{\odot} = 7.3^{+0.5}_{-0.5}$ and Log $M_{outflow}/M_{\odot} \sim 6.3$) but the majority of the luminosity (Log $L_{nuclear}/L_{\odot} = 6.90^{+0.16}_{-0.16}$ and Log $L_{outflow}/L_{\odot} \sim 7.2$). The outflow has an anomalously high $L_\mathrm{CO}/L_\mathrm{FIR}$ of $1.7 \times 10^{-3}$ and is almost certainly shock excited.

Since the launch of the Herschel Space Observatory, our understanding about the photo-dissociation regions (PDR) has taken a step forward. In the bandwidth of the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) on board Herschel, ten CO rotational transitions, including J=4-3 to J=13-12, and three fine structure lines, including [CI] 609, [CI] 370, and [NII] 250 micron, are covered. In this paper, we present our findings from the FTS observations at the nuclear region of M83, based on the spatially resolved physical parameters derived from the CO spectral line energy distribution (SLED) map and the comparisons with the dust properties and star-formation tracers. We discuss (1) the potential of using [NII] 250 and [CI] 370 micron as star-formation tracers; (2) the reliability of tracing molecular gas with CO; (3) the excitation mechanisms of warm CO; (4) the possibility of studying stellar feedback by tracing the thermal pressure of molecular gas in the nuclear region of M83.

We present a comprehensive spectral and morphological analysis of the remnant of Supernova (SN) 1987A with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA). The non-thermal and thermal components of the radio emission are investigated in images from 94 to 672 GHz ($\lambda$ 3.2 mm to 450 $\mu$m), with the assistance of a high-resolution 44 GHz synchrotron template from the ATCA, and a dust template from ALMA observations at 672 GHz. An analysis of the emission distribution over the equatorial ring in images from 44 to 345 GHz highlights a gradual decrease of the east-to-west asymmetry ratio with frequency. We attribute this to the shorter synchrotron lifetime at high frequencies. Across the transition from radio to far infrared, both the synchrotron/dust-subtracted images and the spectral energy distribution (SED) suggest additional emission beside the main synchrotron component ($S_{\nu}\propto\nu^{-0.73}$) and the thermal component originating from dust grains at $T\sim22$ K. This excess could be due to free-free flux or emission from grains of colder dust. However, a second flat-spectrum synchrotron component appears to better fit the SED, implying that the emission could be attributed to a pulsar wind nebula (PWN). The residual emission is mainly localised west of the SN site, as the spectral analysis yields $-0.4\lesssim\alpha\lesssim-0.1$ across the western regions, with $\alpha\sim0$ around the central region. If there is a PWN in the remnant interior, these data suggest that the pulsar may be offset westward from the SN position.

The 12CO J=4-3 to J=13-12 lines of the interstellar medium from nearby galaxies, newly observable with the Herschel SPIRE Fourier Transform Spectrometer (FTS), offer an opportunity to study warmer, more luminous molecular gas than that traced by 12CO J=1-0. Here we present a survey of 17 nearby infrared-luminous galaxy systems (21 pointings). In addition to photometric modeling of dust, we modeled full 12CO spectral line energy distributions from J=1-0 to J=13-12 with two components of warm and cool CO gas, and included LTE analysis of [CI], [CII], [NII] and H2 lines. CO is emitted from a low-pressure/high-mass component traced by the low-J lines and a high-pressure/low-mass component which dominates the luminosity. We found that, on average, the ratios of the warm/cool pressure, mass, and 12CO luminosity are 60 +/- 30, 0.11 +/- 0.02, and 15.6 +/- 2.7. The gas-to-dust-mass ratios are < 120 throughout the sample. The 12CO luminosity is dominated by the high-J lines and is 4 $\times 10^{-4}$ LFIR on average. We discuss systematic effects of single-component and multi-component CO modeling (e.g., single-component J < 3 models overestimate gas pressure by ~ 0.5 dex), as well as compare to Galactic star-forming regions. With this comparison, we show the molecular interstellar medium of starburst galaxies is not simply an ensemble of Galactic-type GMCs. The warm gas emission is likely dominated by regions resembling the warm extended cloud of Sgr B2.

We present the first extragalactic detections of several CH rotational transitions in the far-infrared (FIR) in four nearby galaxies: NGC 1068, Arp 220, M 82 and NGC 253 using the \textit{Herschel Space Observatory}. The CH lines in all four galaxies are a factor of 2 - 4 brighter than the adjacent HCN and HCO+ J = 6-5 lines (also detected in the same spectra). In the star formation dominated galaxies, M 82, NGC 253 and Arp 220, the CH/CO abundance ratio is low (1e-5), implying that the CH is primarily arising in diffuse and translucent gas where the chemistry is driven by UV radiation as found in the Milky Way ISM. In NGC 1068, which has a luminous AGN, the CH/CO ratio is an order of magnitude higher suggesting that CH formation is driven by an X-ray dominated region. Our XDR models show that both the CH and CO abundances in NGC 1068 can be explained by an XDR-driven chemistry for gas densities and molecular hydrogen column densities that are well constrained by the CO observations. We conclude that the CH/CO ratio may a good indicator of the presence of AGN in galaxies. We also discuss the feasibility of detecting CH in intermediate- to high-z galaxies with ALMA.

We present Herschel SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205\mu m$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63\mu m$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $\alpha_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$ in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63\mu m$ to be $\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength $G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.

We present observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. We detect CO transitions from J_up = 4 to 12, as well as the two [CI] fine structure lines at 492 and 809 GHz and the [NII] 461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n_H2 ~ 10^3.2 - 10^3.9 cm^-3 and T_kin ~ 300 - 800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, however it is possibly related to the star-formation activity in these galaxies. Our modeling of the [CI] emission suggests that it is produced in cold (T_kin < 30 K) and dense (n_H2 > 10^3 cm^-3) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J=1-0 transition at 1232 GHz is detected in absorption in UGC05101 and in emission in NGC7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC7130.

We study for the first time the complete sub-millimeter spectra (450 GHz to 1550 GHz) of a sample of nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. The CO ladder (from Jup = 4 to 12) is the most prominent spectral feature in this range. These CO lines probe warm molecular gas that can be heated by ultraviolet photons, shocks, or X-rays originated in the active galactic nucleus or in young star-forming regions. In these proceedings we investigate the physical origin of the CO emission using the averaged CO spectral line energy distribution (SLED) of six Seyfert galaxies. We use a radiative transfer model assuming an isothermal homogeneous medium to estimate the molecular gas conditions. We also compare this CO SLED with the predictions of photon and X-ray dominated region (PDR and XDR) models.

The first complete submillimetre spectrum (190-670um) of the Seyfert 2 galaxy NGC1068 has been observed with the SPIRE Fourier Transform Spectrometer onboard the {\it Herschel} Space Observatory. The sequence of CO lines (Jup=4-13), lines from water, the fundamental rotational transition of HF, two o-H_2O+ lines and one line each from CH+ and OH+ have been detected, together with the two [CI] lines and the [NII]205um line. The observations in both single pointing mode with sparse image sampling and in mapping mode with full image sampling allow us to disentangle two molecular emission components, one due to the compact circum-nuclear disk (CND) and one from the extended region encompassing the star forming ring (SF-ring). Radiative transfer models show that the two CO components are characterized by density of n(H_2)=10^4.5 and 10^2.9 cm^-3 and temperature of T=100K and 127K, respectively. The comparison of the CO line intensities with photodissociation region (PDR) and X-ray dominated region (XDR) models, together with other observational constraints, such as the observed CO surface brightness and the radiation field, indicate that the best explanation for the CO excitation of the CND is an XDR with density of n(H_2) 10^4 cm^-3 and X-ray flux of 9 erg s^-1 cm^-2, consistent with illumination by the active galactic nucleus, while the CO lines in the SF-ring are better modeled by a PDR. The detected water transitions, together with those observed with the \her \sim PACS Spectrometer, can be modeled by an LVG model with low temperature (T_kin \sim 40K) and high density (n(H_2) in the range 10^6.7-10^7.9 cm^-3).

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

We present new Herschel-SPIRE imaging spectroscopy (194-671 microns) of the bright starburst galaxy M82. Covering the CO ladder from J=4-3 to J=13-12, spectra were obtained at multiple positions for a fully sampled ~ 3 x 3 arcminute map, including a longer exposure at the central position. We present measurements of 12CO, 13CO, [CI], [NII], HCN, and HCO+ in emission, along with OH+, H2O+ and HF in absorption and H2O in both emission and absorption, with discussion. We use a radiative transfer code and Bayesian likelihood analysis to model the temperature, density, column density, and filling factor of multiple components of molecular gas traced by 12CO and 13CO, adding further evidence to the high-J lines tracing a much warmer (~ 500 K), less massive component than the low-J lines. The addition of 13CO (and [CI]) is new and indicates that [CI] may be tracing different gas than 12CO. No temperature/density gradients can be inferred from the map, indicating that the single-pointing spectrum is descriptive of the bulk properties of the galaxy. At such a high temperature, cooling is dominated by molecular hydrogen. Photon-dominated region (PDR) models require higher densities than those indicated by our Bayesian likelihood analysis in order to explain the high-J CO line ratios, though cosmic-ray enhanced PDR models can do a better job reproducing the emission at lower densities. Shocks and turbulent heating are likely required to explain the bright high-J emission.

We report the detection of an extremely bright ($\sim$37 mJy at 1100 $\mu$m and $\sim$91 mJy at 880 $\mu$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $\mu$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $\mu$m and SMA 880 $\mu$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $\sigma$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs.

We present Herschel SPIRE-FTS observations of Arp~220, a nearby ULIRG. The FTS continuously covers 190 -- 670 microns, providing a good measurement of the continuum and detection of several molecular and atomic species. We detect luminous CO (J = 4-3 to 13-12) and water ladders with comparable total luminosity; very high-J HCN absorption; OH+, H2O+, and HF in absorption; and CI and NII. Modeling of the continuum yields warm dust, with T = 66 K, and an unusually large optical depth of ~5 at 100 microns. Non-LTE modeling of the CO shows two temperature components: cold molecular gas at T ~ 50 K and warm molecular gas at T ~1350 K. The mass of the warm gas is 10% of the cold gas, but dominates the luminosity of the CO ladder. The temperature of the warm gas is in excellent agreement with H2 rotational lines. At 1350 K, H2 dominates the cooling (~20 L_sun/M_sun) in the ISM compared to CO (~0.4 L_sun/M_sun). We found that only a non-ionizing source such as the mechanical energy from supernovae and stellar winds can excite the warm gas and satisfy the energy budget of ~20 L_sun/M_sun. We detect a massive molecular outflow in Arp 220 from the analysis of strong P-Cygni line profiles observed in OH+, H2O+, and H2O. The outflow has a mass > 10^{7} M_sun and is bound to the nuclei with velocity < 250 km/s. The large column densities observed for these molecular ions strongly favor the existence of an X-ray luminous AGN (10^{44} ergs/s) in Arp 220.

We present the rest-frame 200--320 \mm\ spectrum of the z=3.91 quasar \apm, obtained with Z-Spec at the Caltech Submillimeter Observatory. In addition to the \jeight\ to \jthirteen\ CO rotational transitions which dominate the CO cooling, we find six transitions of water originating at energy levels ranging up to 643 K. Most are first detections at high redshift, and we have confirmed one transition with CARMA. The CO cooling is well-described by our XDR model, assuming L$_{\rm 1-100\,keV}\sim1\times10^{46}\rm\,erg\,s^{-1}$, and that the gas is distributed over a 550-pc sizescale, per the now-favored $\mu$=4 lensing model. The total observed cooling in water corresponds to 6.5$\times10^{9}$ \ls, comparable to that of CO. We compare the water spectrum with that of Mrk 231, finding that the intensity ratios among the high-lying lines are similar, but with a total luminosity scaled up by a factor of $\sim$50. Using this scaling, we estimate an average water abundance relative to \hh\ of 1.4$\times10^{-7}$, a good match to the prediction of the chemical network in the XDR model. As with Mrk 231, the high-lying water transitions are excited radiatively via absorption in the rest-frame far-infrared, and we show that the powerful dust continuum in \apm\ is more than sufficient to pump this massive reservoir of warm water vapor.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We present a 190-307 GHz broadband spectrum obtained with Z-Spec of NGC 1068 with new measurements of molecular rotational transitions. After combining our measurements with those previously published and considering the specific geometry of this Seyfert 2 galaxy, we conduct a multi-species Bayesian likelihood analysis of the density, temperature, and relative molecular abundances of HCN, HNC, CS, and HCO+. We find that these molecules trace warm (T > 100 K) gas of H2 number densities 10^4.2 - 10^4.9 cm^-3. Our models also place strong constraints on the column densities and relative abundances of these molecules, as well as on the total mass in the circumnuclear disk. Using the uniform calibration afforded by the broad Z-Spec bandpass, we compare our line ratios to X-ray dominated region (XDR) and photon-dominated region models. The majority of our line ratios are consistent with the XDR models at the densities indicated by the likelihood analysis, lending substantial support to the emerging interpretation that the energetics in the circumnuclear disk of NGC 1068 are dominated by accretion onto an active galactic nucleus.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

We present a three-pointing study of the molecular gas in the starburst nucleus of M82 based on 190 - 307 GHz spectra obtained with Z-Spec at the Caltech Submillimeter Observatory. We present intensity measurements, detections and upper limits, for 20 transitions, including several new detections of CS, HNC, C2H, H2CO, and CH3CCH lines. We combine our measurements with previously-published measurements at other frequencies for HCN, HNC, CS, C34S, and HCO+ in a multi-species likelihood analysis constraining gas mass, density and temperature, and the species' relative abundances. We find some 1.7 - 2.7 x 10^8 M_sun of gas with n_H2 between 1 - 6 x 10^4 cm^-3 and T > 50 K. While the mass and temperature are comparable to values inferred from mid-J CO transitions, the thermal pressure is a factor of 10 - 20 greater. The molecular interstellar medium is largely fragmented and is subject to ultraviolet irradiation from the star clusters. It is also likely subject to cosmic rays and mechanical energy input from the supernovae, and is warmer on average than the molecular gas in the massive star formation regions in the Milky Way. The typical conditions in the dense gas in M82's central kpc appear unfavorable for further star formation; if any appreciable stellar populations are currently forming, they are likely biased against low mass stars, producing a top-heavy initial mass function.