ALMA cycle 3 observations of $^{12}$CO $ J = 3\rightarrow2$, $^{13}$CO $J = 4\rightarrow 3$, SiO J = $8 \rightarrow 7$, and HCN J = $ 5 \rightarrow 4$ are presented. Significant extended emission is detected in $^{12}$CO J $ = 3\rightarrow2$ with a morphology that is indicative of m = 2 tidal features, suggesting gas inflow. In addition, outflow for both nuclei are found in the $^{12}$CO J $ = 3\rightarrow2$. Significant SiO absorption is detected in the western nucleus. HCN that is morphologically distinct from CO is detected in both nuclei. These observations are compared to non-LTE radiative transfer models created using the Line Modeling Engine (LIME) for simple gas dynamics to gain insight into how physical parameters, such as rotational velocity, turbulent velocity, gas temperature, dust temperature, and gas mass can reproduce the observed kinematic and spatial features. The eastern nucleus is found to be best modeled with an inclusion of a temperature asymmetry from one side of the disk to the other. It is also found that the western nucleus is optically thick even in the less abundant species of $^{13}$CO, absorbing significant amounts of continuum radiation.

Astrophysics spans an enormous range of questions on scales from individual planets to the entire cosmos. To address the richness of 21st century astrophysics requires a corresponding richness of telescopes spanning all bands and all messengers. Much scientific benefit comes from having the multi-wavelength capability available at the same time. Most of these bands,or measurement sensitivities, require space-based missions. Historically, NASA has addressed this need for breadth with a small number of flagship-class missions and a larger number of Explorer missions. While the Explorer program continues to flourish, there is a large gap between Explorers and strategic missions. A fortunate combination of new astrophysics technologies with new, high capacity, low dollar-per-kg to orbit launchers, and new satellite buses allow for cheaper missions with capabilities approaching strategic mission levels. NASA has recognized these developments by calling for Probe-class mission ideas for mission studies, spanning most of the electromagnetic spectrum from GeV gamma-rays to the far infrared, and the new messengers of neutrinos and ultra-high energy cosmic rays. The key insight from the Probes exercise is that order-of-magnitude advances in science performance metrics are possible across the board for initial total cost estimates in the range 500M-1B dollars.

SuperSpec is an on-chip filter-bank spectrometer designed for wideband moderate-resolution spectroscopy at millimeter wavelengths, employing TiN kinetic inductance detectors. SuperSpec technology will enable large-format spectroscopic integral field units suitable for high-redshift line intensity mapping and multi-object spectrographs. In previous results we have demonstrated noise performance in individual detectors suitable for photon noise limited ground-based observations at excellent mm-wave sites. In these proceedings we present the noise performance of a full $R\sim 275$ spectrometer measured using deployment-ready RF hardware and software. We report typical noise equivalent powers through the full device of $\sim 3 \times 10^{-16} \ \mathrm{W}/\sqrt{\mathrm{Hz}}$ at expected sky loadings, which are photon noise dominated. Based on these results, we plan to deploy a six-spectrometer demonstration instrument to the Large Millimeter Telescope in early 2020.

We investigate how the properties of Galactic giant molecular clouds (GMCs) and their denser substructures (clumps) correlate with the local star formation rate. We trace clouds using the $^{12}$CO(3-2) transition, as observed by the CO High Resolution Survey (COHRS). We identify their constituent clumps using thermal dust emission, as observed by the Herschel infrared GALactic plane survey (Hi-GAL). We estimate star formation rates in these clouds using 70 $\mu$m emission. In total, we match 3,674 clumps to 473 clouds in position-position-velocity space spanning the Galactic longitude range $10^\circ<\ell<56^\circ$. We find that more massive clouds produce more clumps and more massive clumps. These clumps have average number densities an order of magnitude greater than their host clouds. We find a mean clump mass fraction of $0.20^{+0.13}_{-0.10}$. This mass fraction weakly varies with mass and mass surface density of clouds, and shows no clear dependence on the virial parameter and line width of the clouds. The average clump mass fraction is only weakly dependent upon Galactocentric radius. Although the scatter in our measured properties is significant, the star formation rate for clouds is independent of clump mass fraction. However, there is a positive correlation between the depletion times for clouds and clump mass fraction. We find a star formation efficiency per free fall time of $\epsilon_{\mathrm{ff}}=0.15\%$ for GMCs but $\epsilon_{\mathrm{ff}}=0.37\%$ for clumps.

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.

We define the molecular cloud properties of the Milky Way first quadrant using data from the JCMT CO(3-2) High Resolution Survey. We apply the Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES) algorithm to extract objects from the full-resolution dataset, creating the first catalog of molecular clouds with a large dynamic range in spatial scale. We identify $>85\,000$ clouds with two clear sub-samples: $\sim35\,500$ well-resolved objects and $\sim540$ clouds with well-defined distance estimations. Only 35% of the cataloged clouds (as well as the total flux encompassed by them) appear enclosed within the Milky Way spiral arms. The scaling relationships between clouds with known distances are comparable to the characteristics of the clouds identified in previous surveys. However, these relations between integrated properties, especially from the full catalog, show a large intrinsic scatter ($\sim0.5$ dex), comparable to other cloud catalogs of the Milky Way and nearby galaxies. The mass distribution of molecular clouds follows a truncated-power law relationship over three orders of magnitude in mass with a form $dN/dM \propto M^{-1.7}$ with a clearly defined truncation at an upper mass of $M_0 \sim 3 \times 10^6~M_\odot$, consistent with theoretical models of cloud formation controlled by stellar feedback and shear. Similarly, the cloud population shows a power-law distribution of size with $dN/dR \propto R^{-2.8}$ with a truncation at $R_0 = 70$ pc.

In the quest to understand high-mass star formation, it is necessary to understand that from which the high-mass stars will form -- the dense molecular gas clouds and clumps. The Herschel infrared GALactic plane survey (Hi-GAL) is a comprehensive survey of thermal dust emission that can be used to characterize the properties and Galactic distribution of molecular gas. We have analyzed the survey maps within the Galactic longitude range $10^\circ < \ell < 56^\circ$ and have transformed these 2D maps into a 3D dense gas map of the Galactic plane using distance probability density functions. This range corresponds to the extent of the Galactic Ring Survey (GRS), which provided the majority of our line-of-sight velocities, and thus kinematic distances. In this section of Hi-GAL, we identify 19,886 clumps, out of which 10,124 have, at minimum, a kinematic distance, and 5,405 have their distances well constrained through Bayesian techniques. Of those with well-constrained distances, clump masses tend to decrease with Galactocentric radius, whereas clump radii are independent of Galactocentric radius.

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.

As the precursors to stellar clusters, it is imperative that we understand the distribution and physical properties of dense molecular gas clouds and clumps. Such a study has been done with the ground-based Bolocam Galactic Plane Survey (BGPS). Now the Herschel infrared GALactic plane survey (Hi-GAL) allows us to do the same with higher quality data and complete coverage of the Galactic plane. We have made a pilot study comparing dense molecular gas clumps identified in the Hi-GAL and BGPS surveys, using six $2^\circ \times 2^\circ$ regions centered at Galactic longitudes of $\ell = 11^\circ$, 30$^\circ$, 41$^\circ$, 50$^\circ$, 202$^\circ$, and $217^\circ$. We adopted the BGPS methodology for identifying clumps and estimating distances, leading to 6198 clumps being identified in our substudy, with 995 of those having well-constrained distances. These objects were evenly distributed with Galactic longitude, a consequence of Hi-GAL being source confusion limited. These clumps range in mass from 10$^{-2}$ M$_\odot$ to 10$^{5}$ M$_\odot$, and have heliocentric distances of up to 16 kpc. When clumps found in both surveys are compared, we see that distances agree within 1 kpc and ratios of masses are of order unity. This serves as an external validation for BGPS and instills confidence as we move forward to cataloging the clumps from the entirety of Hi-GAL. In addition to the sources that were in common with BGPS, Hi-GAL found many additional sources, primarily due to the lack of atmospheric noise. We expect Hi-GAL to yield $2\times10^5$ clumps, with 20% having well-constrained distances, an order of magnitude above what was found in BGPS.

Selecting sources with rising flux densities towards longer wavelengths from Herschel/SPIRE maps is an efficient way to produce a catalogue rich in high-redshift (z > 4) dusty star-forming galaxies. The effectiveness of this approach has already been confirmed by spectroscopic follow-up observations, but the previously available catalogues made this way are limited by small survey areas. Here we apply a map-based search method to 274 deg$^2$ of the HerMES Large Mode Survey (HeLMS) and create a catalogue of 477 objects with SPIRE flux densities $S_{500} > S_{350} >S_{250}$ and a $5 \sigma$ cut-off $S_{500} > $ 52 mJy. From this catalogue we determine that the total number of these "red" sources is at least an order of magnitude higher than predicted by galaxy evolution models. These results are in agreement with previous findings in smaller HerMES fields; however, due to our significantly larger sample size we are also able to investigate the shape of the red source counts for the first time. We have obtained spectroscopic redshift measurements for two of our sources using the Atacama Large Millimeter/submillimeter Array (ALMA). The redshifts z = 5.1 and z = 3.8 confirm that with our selection method we can indeed find high-redshift dusty star-forming galaxies.

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 present Sunyaev-Zel'dovich (SZ) effect measurements from wide-field images towards the galaxy cluster RX J1347.5-1145 obtained from the Caltech Submillimeter Observatory with the Multiwavelength Submillimeter Inductance Camera (MUSIC) at 147, 213, 281, and 337 GHz and with Bolocam at 140 GHz. As part of our analysis, we have used higher frequency data from Herschel-SPIRE and previously published lower frequency radio data to subtract the signal from the brightest dusty star-forming galaxies behind RX J1347.5-1145 and from the AGN in RX J1347.5-1145's BCG. Using these five-band SZ effect images, combined with X-ray spectroscopic measurements of the temperature of the intra-cluster medium (ICM) from Chandra, we constrain the ICM optical depth to be $\tau_e = 7.33^{+0.96}_{-0.97} \times 10^{-3}$ and the ICM line of sight peculiar velocity to be $v_{pec} = -1040^{+870}_{-840}$ km s$^{-1}$. The errors for both quantities are limited by measurement noise rather than calibration uncertainties or astrophysical contamination, and significant improvements are possible with deeper observations. Our best-fit velocity is in good agreement with one previously published SZ effect analysis and in mild tension with the other, although some or all of that tension may be because that measurement samples a much smaller cluster volume. Furthermore, our best-fit optical depth implies a gas mass slightly larger than the Chandra-derived value, implying the cluster is elongated along the line of sight.

We sort $4683$ molecular clouds between $10^\circ< \ell <65^\circ$ from the Bolocam Galactic Plane Survey based on observational diagnostics of star formation activity: compact $70$ $\mu{\rm m}$ sources, mid-IR color-selected YSOs, ${\rm H_2O}$ and ${\rm CH_3OH}$ masers, and UCHII regions. We also present a combined ${\rm NH_3}$-derived gas kinetic temperature and ${\rm H_2O}$ maser catalog for $1788$ clumps from our own GBT 100m observations and from the literature. We identify a subsample of $2223$ ($47.5\%$) starless clump candidates, the largest and most robust sample identified from a blind survey to date. Distributions of flux density, flux concentration, solid angle, kinetic temperature, column density, radius, and mass show strong ($>1$ dex) progressions when sorted by star formation indicator. The median starless clump candidate is marginally sub-virial ($\alpha \sim 0.7$) with $>75\%$ of clumps with known distance being gravitationally bound ($\alpha < 2$). These samples show a statistically significant increase in the median clump mass of $\Delta M \sim 170-370$ M$_\odot$ from the starless candidates to clumps associated with protostars. This trend could be due to (i) mass growth of the clumps at $\dot{M}\sim200-440$ Msun Myr$^{-1}$ for an average free-fall $0.8$ Myr time-scale, (ii) a systematic factor of two increase in dust opacity from starless to protostellar phases, (iii) and/or a variation in the ratio of starless to protostellar clump lifetime that scales as $\sim M^{-0.4}$. By comparing to the observed number of ${\rm CH_3OH}$ maser containing clumps we estimate the phase-lifetime of massive ($M>10^3$ M$_\odot$) starless clumps to be $0.37 \pm 0.08 \ {\rm Myr} \ (M/10^3 \ {\rm M}_\odot)^{-1}$; the majority ($M<450$ M$_\odot$) have phase-lifetimes longer than their average free-fall time.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

SuperSpec is a novel on-chip spectrometer we are developing for multi-object, moderate resolution (R = 100 - 500), large bandwidth (~1.65:1) submillimeter and millimeter survey spectroscopy of high-redshift galaxies. The spectrometer employs a filter bank architecture, and consists of a series of half-wave resonators formed by lithographically-patterned superconducting transmission lines. The signal power admitted by each resonator is detected by a lumped element titanium nitride (TiN) kinetic inductance detector (KID) operating at 100 - 200 MHz. We have tested a new prototype device that achieves the targeted R = 100 resolving power, and has better detector sensitivity and optical efficiency than previous devices. We employ a new method for measuring photon noise using both coherent and thermal sources of radiation to cleanly separate the contributions of shot and wave noise. We report an upper limit to the detector NEP of $1.4\times10^{-17}$ W Hz$^{-1/2}$, within 10% of the photon noise limited NEP for a ground-based R=100 spectrometer.

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.

We use the distance probability density function (DPDF) formalism of Ellsworth-Bowers et al. (2013, 2015) to derive physical properties for the collection of 1,710 Bolocam Galactic Plane Survey (BGPS) version 2 sources with well-constrained distance estimates. To account for Malmquist bias, we estimate that the present sample of BGPS sources is 90% complete above 400 $M_\odot$ and 50% complete above 70 $M_\odot$. The mass distributions for the entire sample and astrophysically motivated subsets are generally fitted well by a lognormal function, with approximately power-law distributions at high mass. Power-law behavior emerges more clearly when the sample population is narrowed in heliocentric distance (power-law index $\alpha = 2.0\pm0.1$ for sources nearer than 6.5 kpc and $\alpha = 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq \alpha \leq 2.05$ for various subsamples of sources, intermediate between that of giant molecular clouds and the stellar initial mass function. The fit to the entire sample yields a high-mass power-law $\hat{\alpha} = 1.94_{-0.10}^{+0.34}$. Physical properties of BGPS sources are consistent with large molecular cloud clumps or small molecular clouds, but the fractal nature of the dense interstellar medium makes difficult the mapping of observational categories to the dominant physical processes driving the observed structure. The face-on map of the Galactic disk's mass surface density based on BGPS dense molecular cloud structures reveals the high-mass star-forming regions W43, W49, and W51 as prominent mass concentrations in the first quadrant. Furthermore, we present a 0.25-kpc resolution map of the dense gas mass fraction across the Galactic disk that peaks around 5%.

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 an expanded distance catalog for 1,710 molecular cloud structures identified in the Bolocam Galactic Plane Survey (BGPS) version 2, representing a nearly threefold increase over the previous BGPS distance catalog. We additionally present a new method for incorporating extant data sets into our Bayesian distance probability density function (DPDF) methodology. To augment the dense-gas tracers (e.g., HCO+(3-2), NH3(1,1)) used to derive line-of-sight velocities for kinematic distances, we utilize the Galactic Ring Survey 13CO(1-0) data to morphologically extract velocities for BGPS sources. The outline of a BGPS source is used to select a region of the GRS 13CO data, along with a reference region to subtract enveloping diffuse emission, to produce a line profile of 13CO matched to the BGPS source. For objects with a HCO+(3-2) velocity, \approx 95% of the new 13CO(1-0) velocities agree with that of the dense gas. A new prior DPDF for kinematic distance ambiguity (KDA) resolution, based on a validated formalism for associating molecular cloud structures with known objects from the literature, is presented. We demonstrate this prior using catalogs of masers with trigonometric parallaxes and HII regions with robust KDA resolutions. The distance catalog presented here contains well-constrained distance estimates for 20% of BGPS V2 sources, with typical distance uncertainties \lesssim 0.5 kpc. Approximately 75% of the well-constrained sources lie within 6 kpc of the Sun, concentrated in the Scutum-Centarus arm. Galactocentric positions of objects additionally trace out portions of the Sagittarius, Perseus, and Outer arms in the first and second Galactic quadrants, and we also find evidence for significant regions of interarm dense gas.

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.

Photometric instruments operating at far infrared to millimetre wavelengths often have broad spectral passbands (central wavelength/bandwidth ~ 3 or less), especially those operating in space. A broad passband can result in significant variation of the beam profile and aperture efficiency across the passband, effects which thus far have not generally been taken into account in the flux calibration of such instruments. With absolute calibration uncertainties associated with the brightness of primary calibration standards now in the region of 5% or less, variation of the beam properties across the passband can be a significant contributor to the overall calibration accuracy for extended emission. We present a calibration framework which takes such variations into account for both antenna-coupled and absorber-coupled focal plane architectures. The scheme covers point source and extended source cases, and also the intermediate case of a semi-extended source profile. We apply the new method to the Herschel-SPIRE space-borne photometer.

We present new images of Arp 220 from the Atacama Large Millimeter/submillimeter Array with the highest combination of frequency (691 GHz) and resolution ($0.36 \times 0.20^{\prime\prime}$) ever obtained for this prototypical ultraluminous infrared galaxy. The western nucleus is revealed to contain warm (200 K) dust that is optically thick ($\tau_{434\mu m} = 5.3$), while the eastern nucleus is cooler (80 K) and somewhat less opaque ($\tau_{434\mu m} = 1.7$). We derive full-width half-maximum diameters of $ 76 \times \le 70$ pc and $123 \times 79$ pc for the western and eastern nucleus, respectively. The two nuclei combined account for ($83 ^{+65}_{-38}$ (calibration) $^{+0}_{-34}$ (systematic))% of the total infrared luminosity of Arp 220. The luminosity surface density of the western nucleus ($ \log(\sigma T^4) = 14.3\pm 0.2 ^{+0}_{-0.7}$ in units of L$_\odot$ kpc$^{-2}$) appears sufficiently high to require the presence of an AGN or a "hot starburst", although the exact value depends sensitively on the brightness distribution adopted for the source. Although the role of any central AGN remains open, the inferred mean gas column densities of $0.6-1.8 \times 10^{25}$ cm$^{-2}$ mean that any AGN in Arp 220 must be Compton-thick.

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 a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $\lambda_{\rm rest} <$ 3.5$ \mu$m; a hot dust component peaking at $\lambda_{rest} \sim$ 5$\,\mu$m; and cold dust component which becomes significant for $\lambda_{\rm rest} >$ 50$\,\mu$m. Six objects in the sample are detected at 250 and 350$ \mu$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ \mu$m detections corresponding to a rest-frame 4.5$ \mu$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

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 a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

The Bolocam Galactic Plane Survey (BGPS) is a 1.1 mm continuum survey of dense clumps of dust throughout the Galaxy covering 170 square degrees. We present spectroscopic observations using the Heinrich Hertz Submillimeter Telescope of the dense gas tracers, HCO+ and N2H+ 3-2, for all 6194 sources in the Bolocam Galactic Plane Survey v1.0.1 catalog between 7.5 <= l <= 194 degrees. This is the largest targeted spectroscopic survey of dense molecular gas in the Milky Way to date. We find unique velocities for 3126 (50.5%) of the BGPS v1.0.1 sources observed. Strong N2H+ 3-2 emission (T_{mb} > 0.5 K) without HCO+ 3-2 emission does not occur in this catalog. We characterize the properties of the dense molecular gas emission toward the entire sample. HCO+ is very sub-thermally populated and the 3-2 transitions are optically thick toward most BGPS clumps. The median observed line width is 3.3 km/s consistent with supersonic turbulence within BGPS clumps. We find strong correlations between dense molecular gas integrated intensities and 1.1 mm peak flux and the gas kinetic temperature derived from previously published NH3 observations. These intensity correlations are driven by the sensitivity of the 3-2 transitions to excitation conditions rather than by variations in molecular column density or abundance. We identify a subset of 113 sources with stronger N2H+ than HCO+ integrated intensity, but we find no correlations between the N2H+ / HCO+ ratio and 1.1 mm continuum flux density, gas kinetic temperature, or line width. Self-absorbed profiles are rare (1.3%).

The Sunyaev-Zel'dovich (SZ) effect spectrum crosses through a null where dT_CMB = 0 near nu_0 = 217 GHz. In a cluster of galaxies, nu_0 can be shifted from the canonical thermal SZ effect value by corrections to the SZ effect scattering due to the properties of the inter-cluster medium. We have measured the SZ effect in the hot galaxy cluster RX J1347.5-1145 with Z-Spec, an R ~ 300 grating spectrometer sensitive between 185 and 305 GHz. These data comprise a high spectral resolution measurement around the null of the SZ effect and clearly exhibit the transition from negative to positive dT_CMB over the Z-Spec band. The SZ null position is measured to be nu_0 = 225.8 \pm 2.5 (stat.) \pm 1.2 (sys.) GHz, which differs from the canonical null frequency by 3.0 sigma and is evidence for modifications to the canonical thermal SZ effect shape. Assuming the measured shift in nu_0 is due only to relativistic corrections to the SZ spectrum, we place the limit T_e = 17.1 \pm 5.3 keV from the zero-point measurement alone. By simulating the response of the instrument to the sky, we are able to generate likelihood functions in {y_0, T_e, v_pec} space. For v_pec = 0 km/s, we measure the best fitting SZ model to be y_0 = 4.6 (+0.6, -0.9) x 10^-4, T_e,0 = 15.2 (+12,-7.4) keV. When v_pec is allowed to vary, a most probable value of v_pec = +450 \pm 810 km/s is found.

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-\sigma RMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

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.

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB.

The Spectral and Photometric Imaging Receiver (SPIRE) on the European Space Agency's Herschel Space Observatory utilizes a pioneering design for its imaging spectrometer in the form of a Fourier Transform Spectrometer (FTS). The standard FTS data reduction and calibration schemes are aimed at objects with either a spatial extent much larger than the beam size or a source that can be approximated as a point source within the beam. However, when sources are of intermediate spatial extent, neither of these calibrations schemes is appropriate and both the spatial response of the instrument and the source's light profile must be taken into account and the coupling between them explicitly derived. To that end, we derive the necessary corrections using an observed spectrum of a fully extended source with the beam profile and the source's light profile taken into account. We apply the derived correction to several observations of planets and compare the corrected spectra with their spectral models to study the beam coupling efficiency of the instrument in the case of partially extended sources. We find that we can apply these correction factors for sources with angular sizes up to \theta_{D} ~ 17". We demonstrate how the angular size of an extended source can be estimated using the difference between the sub-spectra observed at the overlap bandwidth of the two frequency channels in the spectrometer, at 959<\nu<989 GHz. Using this technique on an observation of Saturn, we estimate a size of 17.2", which is 3% larger than its true size on the day of observation. Finally, we show the results of the correction applied on observations of a nearby galaxy, M82, and the compact core of a Galactic molecular cloud, Sgr B2.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

We present a new distance estimation method for dust-continuum-identified molecular cloud clumps. Recent (sub-)millimeter Galactic plane surveys have cataloged tens of thousands of these objects, but detailed study of their physical properties requires robust distance determinations. We derive Bayesian distance probability density functions (DPDFs) for 770 objects from the Bolocam Galactic Plane Survey in the longitude range 7.5 < l < 65. The DPDF formalism is based on kinematic distances, and uses external data sets to place prior distance probabilities to resolve the kinematic distance ambiguity (KDA) for objects in the inner Galaxy. We present priors related to the mid-infrared absorption of dust in dense molecular regions and the Galactic distribution of molecular gas. By assuming a numerical model of Galactic mid-infrared emission and simple radiative transfer, we match the morphology of millimeter thermal dust emission with mid-infrared absorption to compute a prior DPDF for distance discrimination. Selecting objects first from (sub-)millimeter source catalogs avoids a bias towards the darkest infrared dark clouds (IRDCs) and extends the range of heliocentric distance and contrast range probed by mid-infrared extinction. We derive well-constrained KDA resolutions for 618 molecular cloud clumps, with approximately 15% placed at or beyond the tangent distance. Objects sufficiently dark to be cataloged as IRDCs are generally placed at the near kinematic distance. Distance comparisons with Galactic Ring Survey KDA resolutions yield a 92% agreement, and a Galactic face-on view reveals sections of the Sagittarius and Scutum-Centaurus Arms. This KDA-resolution method for large catalogs of sources through the combination of (sub-)millimeter and mid-infrared observations of molecular cloud clumps is generally applicable to other dust-continuum Galactic plane surveys.

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.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

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.

Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star formation in the Universe. In order to understand the physical processes involved in this mode of star formation, it is essential to gain knowledge about the dark matter structures - most importantly their masses - that sub-millimetre galaxies live in. Here we use the magnification effect of gravitational lensing to determine the average mass and dust content of sub-millimetre galaxies with 250mu flux densities of S_250>15mJy selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions of background Lyman-break galaxies at z~3-5 selected using optical data from the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation signal at the 7-sigma level over a sky area of one square degree, with ~80% of this signal being due to magnification, whereas the remaining ~20% comes from dust extinction. Adopting some simple assumptions for the dark matter and dust profiles and the redshift distribution enables us to estimate the average mass of the halos hosting the sub-millimetre galaxies to be log(M_200/M_sun)=13.17+0.05-0.08(stat.) and their average dust mass fraction (at radii of >10kpc) to be M_dust/M_200~6x10^-5. This supports the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside in massive group-sized halos, and are a crucial phase in the assembly and evolution of structure in the Universe.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

We present the status of MUSIC, the MUltiwavelength Sub/millimeter Inductance Camera, a new instrument for the Caltech Submillimeter Observatory. MUSIC is designed to have a 14', diffraction-limited field-of-view instrumented with 2304 detectors in 576 spatial pixels and four spectral bands at 0.87, 1.04, 1.33, and 1.98 mm. MUSIC will be used to study dusty star-forming galaxies, galaxy clusters via the Sunyaev-Zeldovich effect, and star formation in our own and nearby galaxies. MUSIC uses broadband superconducting phased-array slot-dipole antennas to form beams, lumped-element on-chip bandpass filters to define spectral bands, and microwave kinetic inductance detectors to senseincoming light. The focal plane is fabricated in 8 tiles consisting of 72 spatial pixels each. It is coupled to the telescope via an ambient temperature ellipsoidal mirror and a cold reimaging lens. A cold Lyot stop sits at the image of the primary mirror formed by the ellipsoidal mirror. Dielectric and metal mesh filters are used to block thermal infrared and out-of-band radiation. The instrument uses a pulse tube cooler and 3He/3He/4He closed-cycle cooler to cool the focal plane to below 250 mK. A multilayer shield attenuates Earth's magnetic field. Each focal plane tile is read out by a single pair of coaxes and a HEMT amplifier. The readout system consists of 16 copies of custom-designed ADC/DAC and IF boards coupled to the CASPER ROACH platform. We focus on recent updates on the instrument design and results from the commissioning of the full camera in 2012.

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.

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $\mu$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Bo\"{o}tes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $\mu$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $\mu$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

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 describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

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.

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

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.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

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 report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

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.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

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.

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 {\mu}m images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 {\mu}m background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 {\mu}m, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

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.

We examine the physical properties and evolutionary stages of a sample of 17 clumps within 8 Infrared Dark Clouds (IRDCs) by combining existing infrared, millimeter, and radio data with new Bolocam Galactic Plane Survey (BGPS) 1.1 mm data, VLA radio continuum data, and HHT dense gas (HCO+ and N2H+) spectroscopic data. We combine literature studies of star formation tracers and dust temperatures within IRDCs with our search for ultra-compact (UC) HII regions to discuss a possible evolutionary sequence for IRDC clumps. In addition, we perform an analysis of mass tracers in IRDCs and find that 8 micron extinction masses and 1.1 mm Bolocam Galactic Plane Survey (BGPS) masses are complementary mass tracers in IRDCs except for the most active clumps (notably those containing UCHII regions), for which both mass tracers suffer biases. We find that the measured virial masses in IRDC clumps are uniformly higher than the measured dust continuum masses on the scale of ~1 pc. We use 13CO, HCO+, and N2H+ to study the molecular gas properties of IRDCs and do not see any evidence of chemical differentiation between hot and cold clumps on the scale of ~1 pc. However, both HCO+ and N2H+ are brighter in active clumps, due to an increase in temperature and/or density. We report the identification of four UCHII regions embedded within IRDC clumps and find that UCHII regions are associated with bright (>1 Jy) 24 micron point sources, and that the brightest UCHII regions are associated with "diffuse red clumps" (an extended enhancement at 8 micron). The broad stages of the discussed evolutionary sequence (from a quiescent clump to an embedded HII region) are supported by literature dust temperature estimates; however, no sequential nature can be inferred between the individual star formation tracers.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

We present an analysis of an 8 arcminute diameter map of the area around the galaxy cluster Abell 1835 from jiggle map observations at a wavelength of 1.1 mm using the Bolometric Camera (Bolocam) mounted on the Caltech Submillimeter Observatory (CSO). The data is well described by a model including an extended Sunyaev-Zel'dovich (SZ) signal from the cluster gas plus emission from two bright background submm galaxies magnified by the gravitational lensing of the cluster. The best-fit values for the central Compton value for the cluster and the fluxes of the two main point sources in the field: SMM J140104+0252, and SMM J14009+0252 are found to be $y_{0}=(4.34\pm0.52\pm0.69)\times10^{-4}$, 6.5$\pm{2.0}\pm0.7$ mJy and 11.3$\pm{1.9}\pm1.1$ mJy, where the first error represents the statistical measurement error and the second error represents the estimated systematic error in the result. This measurement assumes the presence of dust emission from the cluster's central cD galaxy of $1.8\pm0.5$ mJy, based on higher frequency observations of Abell 1835. The cluster image represents one of the highest-significance SZ detections of a cluster in the positive region of the thermal SZ spectrum to date. The inferred central intensity is compared to other SZ measurements of Abell 1835 and this collection of results is used to obtain values for $y_{0} = (3.60\pm0.24)\times10^{-4}$ and the cluster peculiar velocity $v_{z} = -226\pm275$ km/s.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 {\mu}m) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

We present a detailed analysis of the far-IR properties of the bright, lensed, z = 2.3, SMG, SMM J2135-0102, using new observations with Herschel, SCUBA-2 and the VLA. These data allow us to constrain the galaxy's SED and show that it has an intrinsic rest-frame 8-1000um luminosity, L(bol), of (2.3 +/- 0.2) x 10^12 L(sun) and a likely SFR of ~400 M(sun)/yr. The galaxy sits on the far-IR/radio correlation for far-IR-selected galaxies. At ~>70um, the SED can be described adequately by dust components with T(d) ~ 30 and 60K. Using SPIRE's Fourier Transform Spectrometer we report a detection of the [CII] 158um cooling line. If the [CII], CO and far-IR continuum arise in photo-dissociation regions, we derive a characteristic gas density, n ~ 10^3 cm^-3, and a far-UV radiation field, G_0, 10^3x stronger than the Milky Way. L([CII])/L(bol) is significantly higher than in local ULIRGs but similar to the values found in local star-forming galaxies and starburst nuclei. This is consistent with SMM J2135-0102 being powered by starburst clumps distributed across ~2 kpc, evidence that SMGs are not simply scaled-up ULIRGs. Our results show that SPIRE's FTS has the ability to measure the redshifts of distant, obscured galaxies via the blind detection of atomic cooling lines, but it will not be competitive with ground-based CO-line searches. It will, however, allow detailed study of the integrated properties of high-redshift galaxies, as well as the chemistry of their ISM, once more suitably bright candidates have been found.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 micron. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K-band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, (ii) the connection between the dust content and composition and the other phases of the interstellar medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multi-frequency data to carry out an analysis of the statistical properties of the sample.

We report measurements of the fluctuations in atmospheric emission (atmospheric noise) above Mauna Kea recorded with Bolocam at 143 and 268 GHz from the Caltech Submillimeter Observatory (CSO). The 143 GHz data were collected during a 40 night observing run in late 2003, and the 268 GHz observations were made in early 2004 and early 2005 over a total of 60 nights. Below 0.5 Hz, the data time-streams are dominated by atmospheric noise in all observing conditions. The atmospheric noise data are consistent with a Kolmogorov-Taylor (K-T) turbulence model for a thin wind-driven screen, and the median amplitude of the fluctuations is 280 mK^2 rad^(-5/3) at 143 GHz and 4000 mK^2 rad^(-5/3) at 268 GHz. Comparing our results with previous ACBAR data, we find that the normalization of the power spectrum of the atmospheric noise fluctuations is a factor of 80 larger above Mauna Kea than above the South Pole at millimeter wavelengths. Most of this difference is due to the fact that the atmosphere above the South Pole is much drier than the atmosphere above Mauna Kea. However, the atmosphere above the South Pole is slightly more stable as well: the fractional fluctuations in the column depth of precipitable water vapor are a factor of sqrt(2) smaller at the South Pole compared to Mauna Kea. Based on our atmospheric modeling, we developed several algorithms to remove the atmospheric noise, and the best results were achieved when we described the fluctuations using a low-order polynomial in detector position over the 8 arcmin field of view (FOV). However, even with these algorithms, we were not able to reach photon-background-limited instrument photometer (BLIP) performance at frequencies below 0.5 Hz in any observing conditions.

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

We present the first broadband lambda = 1 mm spectrum toward the z=2.56 Cloverleaf Quasar, obtained with Z-Spec, a 1-mm grating spectrograph on the 10.4-meter Caltech Submillimeter Observatory. The 190-305 GHz observation band corresponds to rest-frame 272 to 444 microns, and we measure the dust continuum as well as all four transitions of carbon monoxide (CO) lying in this range. The power-law dust emission, F_nu = 14 mJy (nu/240GHz)^3.9 is consistent with the published continuum measurements. The CO J=6->5, J=8->7, and J=9->8 measurements are the first, and now provide the highest-J CO information in this source. Our measured CO intensities are very close to the previously-published interferometric measurements of J=7->6, and we use all available transitions and our 13CO upper limits to constrain the physical conditions in the Cloverleaf molecular gas disk. We find a large mass (2-50x10^9 Msun) of highly-excited gas with thermal pressure nT > 10^6 Kcm^-3. The ratio of the total CO cooling to the far-IR dust emission exceeds that in the local dusty galaxies, and we investigate the potential heating sources for this bulk of warm molecular gas. We conclude that both UV photons and X-rays likely contribute, and discuss implications for a top-heavy stellar initial mass function arising in the X-ray-irradiated starburst. Finally we present tentative identifications of other species in the spectrum, including a possible detection of the H20 2_0,2->1_1,1 transition at lambda_rest = 303 microns.

When and how did galaxies form and their metals accumulate? Over the last decade, this has moved from an archeological question to a live investigation: there is now a broad picture of the evolution of galaxies in dark matter halos: their masses, stars, metals and supermassive blackholes. Galaxies have been found and studied in which these formation processes are taking place most vigorously, all the way back in cosmic time to when the intergalactic medium (IGM) was still largely neutral. However, the details of how and why the interstellar medium (ISM) in distant galaxies cools, is processed, recycled and enriched in metals by stars, and fuels active galactic nuclei (AGNs) remain uncertain. In particular, the cooling of gas to fuel star formation, and the chemistry and physics of the most intensely active regions is hidden from view at optical wavelengths, but can be seen and diagnosed at mid- & far-infrared (IR) wavelengths. Rest-frame IR observations are important first to identify the most luminous, interesting and important galaxies, secondly to quantify accurately their total luminosity, and finally to use spectroscopy to trace the conditions in the molecular and atomic gas out of which stars form. In order to map out these processes over the full range of environments and large-scale structures found in the universe - from the densest clusters of galaxies to the emptiest voids - we require tools for deep, large area surveys, of millions of galaxies out to z~5, and for detailed follow-up spectroscopy. The necessary tools can be realized technically. Here, we outline the requirements for gathering the crucial information to build, validate and challenge models of galaxy evolution.

We have surveyed two science fields totaling one square degree with Bolocam at 2.1 mm to search for secondary CMB anisotropies caused by the Sunyaev- Zel'dovich effect (SZE). The fields are in the Lynx and Subaru/XMM SDS1 fields. Our survey is sensitive to angular scales with an effective angular multipole of l_eff = 5700 with FWHM_l = 2800 and has an angular resolution of 60 arcseconds FWHM. Our data provide no evidence for anisotropy. We are able to constrain the level of total astronomical anisotropy, modeled as a flat bandpower in C_l, with frequentist 68%, 90%, and 95% CL upper limits of 590, 760, and 830 uKCMB^2. We statistically subtract the known contribution from primary CMB anisotropy, including cosmic variance, to obtain constraints on the SZE anisotropy contribution. Now including flux calibration uncertainty, our frequentist 68%, 90% and 95% CL upper limits on a flat bandpower in C_l are 690, 960, and 1000 uKCMB^2. When we instead employ the analytic spectrum suggested by Komatsu and Seljak (2002), and account for the non-Gaussianity of the SZE anisotropy signal, we obtain upper limits on the average amplitude of their spectrum weighted by our transfer function of 790, 1060, and 1080 uKCMB^2. We obtain a 90% CL upper limit on sigma8, which normalizes the power spectrum of density fluctuations, of 1.57. These are the first constraints on anisotropy and sigma8 from survey data at these angular scales at frequencies near 150 GHz.

We present an unbiased census of deeply embedded protostars in Perseus, Serpens, and Ophiuchus, assembled by combining large-scale 1.1 mm Bolocam continuum and Spitzer Legacy surveys. We identify protostellar candidates based on their mid-infrared properties, correlate their positions with 1.1 mm core positions, and construct well-sampled SEDs using our extensive wavelength coverage (lam=1.25-1100 micron). Source classification based on the bolometric temperature yields a total of 39 Class 0 and 89 Class I sources in the three cloud sample. We compare to protostellar evolutionary models using the bolometric temperature-luminosity diagram, finding a population of low luminosity Class I sources that are inconsistent with constant or monotonically decreasing mass accretion rates. This result argues strongly for episodic accretion during the Class I phase, with more than 50% of sources in a ``sub-Shu'' (dM/dt < 1e-6 Msun/yr) accretion state. Average spectra are compared to protostellar radiative transfer models, which match the observed spectra fairly well in Stage 0, but predict too much near-IR and too little mid-IR flux in Stage I. Finally, the relative number of Class 0 and Class I sources are used to estimate the lifetime of the Class 0 phase; the three cloud average yields a Class 0 lifetime of 1.7e5 yr, ruling out an extremely rapid early accretion phase. Correcting photometry for extinction results in a somewhat shorter lifetime (1.1e5 yr). In Ophiuchus, however, we find very few Class 0 sources (N(Class0)/N(ClassI)=0.1-0.2), similar to previous studies of that cloud. The observations suggest a consistent picture of nearly constant average accretion rate through the entire embedded phase, with accretion becoming episodic by at least the Class I stage, and possibly earlier.

The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.

We present an unbiased census of starless cores in Perseus, Serpens, and Ophiuchus, assembled by comparing large-scale Bolocam 1.1 mm continuum emission maps with Spitzer c2d surveys. We use the c2d catalogs to separate 108 starless from 92 protostellar cores in the 1.1 mm core samples from Enoch et al. (2006), Young et al. (2006), and Enoch et al. (2007). A comparison of these populations reveals the initial conditions of the starless cores. Starless cores in Perseus have similar masses but larger sizes and lower densities on average than protostellar cores, with sizes that suggest density profiles substantially flatter than r^-2. By contrast, starless cores in Serpens are compact and have lower masses than protostellar cores; future star formation will likely result in lower mass objects than the currently forming protostars. Comparison to dynamical masses estimated from the NH_3 survey of Perseus cores by Rosolowsky et al. (2007) suggests that most of the starless cores are likely to be gravitationally bound, and thus prestellar. The combined prestellar core mass distribution includes 108 cores and has a slope of -2.3+/-0.4 for M>0.8 Msun. This slope is consistent with recent measurements of the stellar initial mass function, providing further evidence that stellar masses are directly linked to the core formation process. We place a lower limit on the core-to-star efficiency of 25%. There are approximately equal numbers of prestellar and protostellar cores in each cloud, thus the dense prestellar core lifetime must be similar to the lifetime of embedded protostars, or 4.5x10^5 years, with a total uncertainty of a factor of two. Such a short lifetime suggests a dynamic, rather than quasi-static, core evolution scenario, at least at the relatively high mean densities (n>2x10^4 cm^-3) to which we are sensitive.

AzTEC is a mm-wavelength bolometric camera utilizing 144 silicon nitride micromesh detectors. Herein we describe the AzTEC instrument architecture and its use as an astronomical instrument. We report on several performance metrics measured during a three month observing campaign at the James Clerk Maxwell Telescope, and conclude with our plans for AzTEC as a facility instrument on the Large Millimeter Telescope.

We have undertaken an unprecedentedly large 1.1 millimeter continuum survey of three nearby star forming clouds using Bolocam at the Caltech Submillimeter Observatory. We mapped the largest areas in each cloud at millimeter or submillimeter wavelengths to date: 7.5 sq. deg in Perseus (Paper I), 10.8 sq. deg in Ophiuchus (Paper II), and 1.5 sq. deg in Serpens with a resolution of 31", detecting 122, 44, and 35 cores, respectively. Here we report on results of the Serpens survey and compare the three clouds. Average measured angular core sizes and their dependence on resolution suggest that many of the observed sources are consistent with power-law density profiles. Tests of the effects of cloud distance reveal that linear resolution strongly affects measured source sizes and densities, but not the shape of the mass distribution. Core mass distribution slopes in Perseus and Ophiuchus (alpha=2.1+/-0.1 and alpha=2.1+/-0.3) are consistent with recent measurements of the stellar IMF, whereas the Serpens distribution is flatter (alpha=1.6+/-0.2). We also compare the relative mass distribution shapes to predictions from turbulent fragmentation simulations. Dense cores constitute less than 10% of the total cloud mass in all three clouds, consistent with other measurements of low star-formation efficiencies. Furthermore, most cores are found at high column densities; more than 75% of 1.1 mm cores are associated with Av>8 mag in Perseus, 15 mag in Serpens, and 20-23 mag in Ophiuchus.

SPIRE, the Spectral and Photometric Imaging Receiver, is one of three instruments to fly on ESA's Herschel Space Observatory. It contains a three-band imaging photometer operating at 250, 360 and 520 microns, and an imaging Fourier Transform Spectrometer (FTS) covering 200-670 microns. It will be used for many extragalactic science programmes, a number of which will be implemented as Herschel Key Projects. The SPIRE consortium's Guaranteed Time (GT) programme will devote more than 1000 hours to Key Projects covering the high-z universe and local galaxies. It is also expected that substantial amounts of Herschel Open Time will be used for further extragalactic investigations. The high-z part of the SPIRE GT programme will focus on blank-field surveys with a range of depths and areas optimised to sample the luminosity-redshift plane and characterise the bolometric luminosity density of the universe at high z. Fields will be selected that are well covered at other wavelengths to facilitate source identifications and enable detailed studies of the redshifts, spectral energy distributions, and infrared properties of detected galaxies. The local galaxies programme will include a detailed spectral and photometric study of a sample of well resolved nearby galaxies, a survey of more than 300 local galaxies designed to provide a statistical survey of dust in the nearby universe, and a study of the ISM in low-metallicity environments, bridging the gap between the local universe and primordial galaxies.

We present 350 micron photometry of all 17 galaxy candidates in the Lockman Hole detected in a 1.1 mm Bolocam survey. Several of the galaxies were previously detected at 850 microns, at 1.2 mm, in the infrared by Spitzer, and in the radio. Nine of the Bolocam galaxy candidates were detected at 350 microns and two new candidates were serendipitously detected at 350 microns (bringing the total in the literature detected in this way to three). Five of the galaxies have published spectroscopic redshifts, enabling investigation of the implied temperature ranges and a comparison of photometric redshift techniques. Lambda = 350 microns lies near the spectral energy distribution peak for z = 2.5 thermally emitting galaxies. Thus, luminosities can be measured without extrapolating to the peak from detection wavelengths of lambda > 850 microns. Characteristically, the galaxy luminosities lie in the range 1.0 - 1.2 x 10^13 L_solar, with dust temperatures in the range of 40 K to 70 K, depending on the choice of spectral index and wavelength of unit optical depth. The implied dust masses are 3 - 5 x 10^8 M_solar. We find that the far-infrared to radio relation for star-forming ULIRGs systematically overpredicts the radio luminosities and overestimates redshifts on the order of Delta z ~ 1, whereas redshifts based on either on submillimeter data alone or the 1.6 micron stellar bump and PAH features are more accurate.

We present a large-scale millimeter continuum map of the Ophiuchus molecular cloud. Nearly 11 square degrees, including all of the area in the cloud with visual extinction more than 3 magnitudes, was mapped at 1.1 mm with Bolocam on the Caltech Submillimeter Observatory (CSO). By design, the map also covers the region mapped in the infrared with the Spitzer Space Telescope. We detect 44 definite sources, and a few likely sources are also seen along a filament in the eastern streamer. The map indicates that dense cores in Ophiuchus are very clustered and often found in filaments within the cloud. Most sources are round, as measured at the half power point, but elongated when measured at lower contour levels, suggesting spherical sources lying within filaments. The masses, for an assumed dust temperature of 10 K, range from 0.24 to 3.9 solar masses, with a mean value of 0.96 solar masses. The total mass in distinct cores is 42 solar masses, 0.5 to 2% of the total cloud mass, and the total mass above 4 sigma is about 80 solar masses. The mean densities in the cores are quite high, with an average of 1.6 x 10^6 per cc, suggesting short free-fall times. The core mass distribution can be fitted with a power law with slope of 2.1 plus or minus 0.3 for M>0.5 solar masses, similar to that found in other regions, but slightly shallower than that of some determinations of the local IMF. In agreement with previous studies, our survey shows that dense cores account for a very small fraction of the cloud volume and total mass. They are nearly all confined to regions with visual extinction at least 9 mag, a lower threshold than found previously.

We have completed a 1.1 mm continuum survey of 7.5 sq deg of the Perseus Molecular Cloud using Bolocam at the Caltech Submillimeter Observatory. This represents the largest millimeter or submillimeter continuum map of Perseus to date. Our map covers more than 30,000 31" (FWHM) resolution elements to a 1 sigma RMS of 15 mJy/beam. We detect a total of 122 cores above a 5 sigma point source mass detection limit of 0.18 M_sun, assuming a dust temperature of 10 K, 60 of which are new millimeter or submillimeter detections. The 1.1 mm mass function is consistent with a broken power law of slope -1.3 (0.5 M_sun<M<2.5 M_sun) and -2.6 (M>2.5 M_sun), similar to the local initial mass function slope. No more than 5% of the total cloud mass is contained in discrete 1.1 mm cores, which account for a total mass of 285 M_sun. We suggest an extinction threshold for millimeter cores of Av~5 mag, based on our calculation of the probability of finding a 1.1 mm core as a function of Av. Much of the cloud is devoid of compact millimeter emission; despite the significantly greater area covered compared to previous surveys, only 5-10 of the newly identified sources lie outside previously observed areas. The two-point correlation function confirms that dense cores in the cloud are highly structured, with significant clustering on scales as large as 2e5 AU. These 1.1 mm results, especially when combined with recently acquired c2d Spitzer Legacy data, will provide a census of dense cores and protostars in Perseus and improve our understanding of the earliest stages of star formation in molecular clouds.

We perform a fluctuation analysis of the 1.1mm Bolocam Lockman Hole Survey, which covers 324 square arcmin to a very uniform point source-filtered RMS noise level of 1.4 mJy/beam. The fluctuation analysis has the significant advantage of utilizing all of the available data. We constrain the number counts in the 1-10 mJy range, and derive significantly tighter constraints than in previous work: the power-law index is 2.7 (+0.18, -0.15), while the amplitude is equal to 1595 (+85,-238) sources per mJy per square degree, or N(>1 mJy) = 940 (+50,-140) sources/square degree (95% confidence). Our results agree extremely well with those derived from the extracted source number counts by Laurent et al (2005). Our derived normalization is about 2.5 times smaller than determined by MAMBO at 1.2mm by Greve et al (2004). However, the uncertainty in the normalization for both data sets is dominated by the systematic (i.e., absolute flux calibration) rather than statistical errors; within these uncertainties, our results are in agreement. We estimate that about 7% of the 1.1mm background has been resolved at 1 mJy.

We present results of a new deep 1.1 mm survey using Bolocam, a millimeter-wavelength bolometer array camera designed for mapping large fields at fast scan rates, without chopping. A map, galaxy candidate list, and derived number counts are presented. This survey encompasses 324 arcmin^2 to an rms noise level (filtered for point sources) of 1.4 mJy/beam and includes the entire regions surveyed by the published 8 mJy 850 micron JCMT SCUBA and 1.2 mm IRAM MAMBO surveys. We reduced the data using a custom software pipeline to remove correlated sky and instrument noise via a principal component analysis. Extensive simulations and jackknife tests were performed to confirm the robustness of our source candidates and estimate the effects of false detections, bias, and completeness. In total, 17 source candidates were detected at a significance > 3.0 sigma, with six expected false detections. Nine candidates are new detections, while eight candidates have coincident SCUBA 850 micron and/or MAMBO 1.2 mm detections. From our observed number counts, we estimate the underlying differential number count distribution of submillimeter galaxies and find it to be in general agreement with previous surveys. Modeling the spectral energy distributions of these submillimeter galaxies after observations of dusty nearby galaxies suggests extreme luminosities of L = 1.0-1.6 x 10^13 L_solar and, if powered by star formation, star formation rates of 500-800 M_solar/yr.

We have monitored the SiO (v = 1, J = 2 -> 1) maser polarization in 17 variable stars (Miras, OH-IR stars, and supergiants) to investigate the long-term persistence of masers. The 8 epochs of observations span 2.5 years, thereby sampling multiple cycles for these stars with typical periods of ~1 year. The average polarization was 23% with a typical dispersion of 7%, although the variability differed substantially from star to star. In the Stokes q and u spectra of individual stars, a few strong maser features tended to dominate the polarization, with the maser features persisting for less than one stellar cycle for some stars, and for multiple cycles for a few stars. Because individual masers are not resolved in our beam-averaged total intensity spectra, we correlated the polarization spectra between epochs to measure the characteristic lifetimes of the features rather than attempting to trace the evolution of separate line components. We found that individual maser feature lifetimes ranged from a few months or less to >2 years. These data indicate that for the sample of stars as a whole, the masers are not reset at a particular stellar phase.

We describe the development of a bolometric receiver designed to measure the arcminute-scale polarization of the cosmic microwave background (CMB). The Polatron will be mounted at the Cassegrain focus of the 5.5 m telescope at the Owens Valley Radio Observatory (OVRO). The receiver will measure both the Q and U Stokes parameters over a 20% pass-band centered near 100 GHz, with the input polarization signal modulated at ~0.6 Hz by a rotating, birefringent, quartz half-wave plate. In six months of observation we plan to observe ~400 2.5 arcminute pixels in a ring about the North Celestial Pole to a precision of ~6 \mu K/pixel in each of Q and U, adequate to unambiguously detect CMB polarization at levels predicted by current models.