{Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

The Heterodyne Instrument for the Far Infrared (HIFI) on board the Herschel Space Observatory is composed of a set of fourteen double sideband mixers. We discuss the general problem of the sideband ratio (SBR) determination and the impact of an imbalanced sideband ratio on the line calibration in double sideband heterodyne receivers. The HIFI SBR is determined from a combination of data taken during pre-launch gas cell tests and in-flight. The results and some of the calibration artefacts discovered in the gas cell test data are presented here along with some examples of how these effects appear in science data taken in orbit.

We reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$\sigma$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$\sigma$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics.

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

We report a new technique to select 1.6<z<4.6 dusty Lyman-alpha emitters (LAEs), over a third of which are `blobs' (LABs) with emission extended on scales of 30-100kpc. Combining data from the NASA Wide-field Infrared Survey Explorer (WISE) mission with optical spectroscopy from the W.M. Keck telescope, we present a color criteria that yields a 78% success rate in identifying rare, dusty LAEs of which at least 37% are LABs. The objects have a surface density of only ~0.1 per square degree, making them rare enough that they have been largely missed in narrow surveys. We measured spectroscopic redshifts for 92 of these WISE-selected, typically radio-quiet galaxies and find that the LAEs (LABs) have a median redshift of 2.3 (2.5). The WISE photometry coupled with data from Herschel reveals that these galaxies have extreme far-infrared luminosities (L_IR>10^{13-14}L_sun) and warm colors, typically larger than submillimeter-selected galaxies (SMGs) and dust-obscured galaxies (DOGs). These traits are commonly associated with the dust being energized by intense AGN activity. We hypothesize that the combination of spatially extended Lyman-alpha, large amounts of warm IR-luminous dust, and rarity (implying a short-lived phase) can be explained if the galaxies are undergoing strong `feedback' transforming them from an extreme dusty starburst to a QSO.

We present optical to far-infrared photometry of 31 reddened QSOs that show evidence for radiatively driven outflows originating from AGN in their rest-frame UV spectra. We use these data to study the relationships between the AGN-driven outflows, and the AGN and starburst infrared luminosities. We find that FeLoBAL QSOs are invariably IR-luminous, with IR luminosities exceeding 10^{12} Solar luminosities in all cases. The AGN supplies 76% of the total IR emission, on average, but with a range from 20% to 100%. We find no evidence that the absolute luminosity of obscured star formation is affected by the AGN-driven outflows. Conversely, we find an anticorrelation between the strength of AGN-driven outflows, as measured from the range of outflow velocities over which absorption exceeds a minimal threshold, and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no convincing evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN can have a dramatic, negative impact on luminous star formation in their host galaxies. We find that such outflows act to curtail star formation such that star formation contributes less than ~25% of the total IR luminosity. We also propose that the degree to which termination of star formation takes place is not deducible from the IR luminosity of the AGN.

We use SDSS spectra and optical to far-infrared photometry for a sample of 31 FeLoBAL QSOs to study the relationship between the AGN-driven outflows, and obscured star formation in their host galaxies. We find that FeLoBAL QSOs invariably have IR luminosities exceeding 10^12 Solar luminosities. The AGN supplies 75% of the total IR emission, on average, but with a range from 20% to 100%. We find a clear anticorrelation between the strength of the AGN-driven outflows and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN act to curtail obscured star formation in the host galaxies of reddened QSOs to less than ~25% of the total IR luminosity. This is the most direct evidence yet obtained for `quasar mode' AGN feedback.

In order to constrain the bolometric luminosities, dust properties and molecular gas content of giant Lyman alpha nebulae, the so-called Lyman alpha blobs, we have carried out a study of dust continuum and CO line emission in two well-studied representatives of this population at z ~ 3: a Lya blob discovered by its strong Spitzer MIPS 24um detection (LABd05; Dey et al. 2005) and the Steidel blob 1 (SSA22-LAB01; Steidel et al. 2000). We find that the spectral energy distribution of LABd05 is well described by an AGN-starburst composite template with L(FIR) = (4.0 +/- 0.5) x 10^12 Lsun, comparable to high-z sub-millimeter galaxies and ultraluminous infrared galaxies. New APEX/LABOCA 870um measurements rule out the reported SCUBA detection of the SSA22-LAB01 (S[850um] = 16.8 mJy) at the > 4sigma level. Consistent with this, ultra-deep Plateau de Bure Interferometer (PdBI) observations with ~2arcsec spatial resolution also fail to detect any 1.2mm continuum source down to ~0.45mJy per beam (3sigma). Combined with the existing (sub)mm observations in the literature, we conclude that the FIR luminosity of SSA22-LAB01 remains uncertain. No CO line is detected in either case down to integrated flux limits of (Snu dV) < 0.25--1.0 Jy km/s, indicating a modest molecular gas reservoir, M(H_2) < 1--3 x 10^10 Msun. The non-detections exclude, with high significance (12 sigma), the previous tentative detection of a CO(4-3) line in the SSA22-LAB01. The increased sensitivity afforded by ALMA will be critical in studying molecular gas and dust in these interesting systems.

We have imaged an 11.5 sq. deg. region of sky towards the South Ecliptic Pole (RA = 04h43m, Dec = -53d40m, J2000) at 24 and 70 microns with MIPS, the Multiband Imaging Photometer for Spitzer. This region is coincident with a field mapped at longer wavelengths by AKARI and the Balloon-borne Large Aperture Submillimeter Telescope. We discuss our data reduction and source extraction procedures. The median depths of the maps are 47 microJy/beam at 24 micron and 4.3 mJy/beam at 70 micron. At 24 micron, we identify 93098 point sources with signal-to-noise ratio (SNR) >5, and an additional 63 resolved galaxies; at 70 micron, we identify 891 point sources with SNR >6. From simulations, we determine a false detection rate of 1.8% (1.1%) for the 24 micron (70 micron) catalog. The 24 and 70 micron point-source catalogs are 80% complete at 230 microJy and 11 mJy, respectively. These mosaic images and source catalogs will be available to the public through the NASA/IPAC Infrared Science Archive.

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 outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

We present SHARC-2 350 micron data on 20 luminous z ~ 2 starbursts with S(1.2mm) > 2 mJy from the Spitzer-selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S(350um) > 25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass--temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6um--2mm measurements of local starbursts. We find characteristic single-component temperatures T1 ~ 35.5+-2.2 K and integrated infrared (IR) luminosities around 10^(12.9+-0.1) Lsun for the SWIRE-selected sources. Molecular gas masses are estimated at 4 x 10^(10) Msun, assuming kappa(850um)=0.15 m^2/kg and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply >~2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 um flux densities, provides an effective means for the study of SMGs at z ~ 1.5--2.5.

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We present a quantitative morphological analysis using HST NICMOS H160- and ACS I775- band imaging of 25 spectroscopically confirmed submillimetre galaxies (SMGs) which have redshifts between z=0.7-3.4. Our analysis also employs a comparison sample of more typical star-forming galaxies at similar redshifts (such as LBGs) which have lower far-infrared luminosities. This is the first large-scale study of the morphologies of SMGs in the near-infrared at ~0.1" resolution (<1kpc). We find that the half light radii of the SMGs (r_h=2.3+/-0.3 and 2.8+/-0.4kpc in the observed I- and H-bands respectively) and asymmetries are not statistically distinct from the comparison sample of star-forming galaxies. However, we demonstrate that the SMG morphologies differ more between the rest-frame UV and optical-bands than typical star-forming galaxies and interpret this as evidence for structured dust obscuration. We show that the composite observed H-band light profile of SMGs is better fit with a Sersic index with n~2, suggesting the stellar structure of SMGs is best described by a spheroid/elliptical galaxy light distribution. We also compare the sizes and stellar masses of SMGs to local and high-redshift populations, and find that the SMGs have stellar densities which are comparable to local early-type galaxies, as well as luminous, red and dense galaxies at z~1.5 which have been proposed as direct SMG descendants, although the SMG stellar masses and sizes are systematically larger. Overall, our results suggest that the physical processes occuring within the galaxies are too complex to be simply characterised by the rest-frame UV/optical morphologies which appear to be essentially decoupled from all other observables, such as bolometric luminosity, stellar or dynamical mass.

We present SHARC-II 350um imaging of twelve 24um-bright (F_24um > 0.8 mJy) Dust-Obscured Galaxies (DOGs) and CARMA 1mm imaging of a subset of 2 DOGs, all selected from the Bootes field of the NOAO Deep Wide-Field Survey. Detections of 4 DOGs at 350um imply IR luminosities which are consistent within a factor of 2 of expectations based on a warm dust spectral energy distribution (SED) scaled to the observed 24um flux density. The 350um upper limits for the 8 non-detected DOGs are consistent with both Mrk231 and M82 (warm dust SEDs), but exclude cold dust (Arp220) SEDs. The two DOGs targeted at 1mm were not detected in our CARMA observations, placing strong constraints on the dust temperature: T_dust > 35-60 K. Assuming these dust properties apply to the entire sample, we find dust masses of ~3x10^8 M_sun. In comparison to other dusty z ~ 2 galaxy populations such as sub-millimeter galaxies (SMGs) and other Spitzer-selected high-redshift sources, this sample of DOGs has higher IR luminosities (2x10^13 L_sun vs. 6x10^12 L_sun for the other galaxy populations), warmer dust temperatures (>35-60 K vs. ~30 K), and lower inferred dust masses (3x10^8 M_sun vs. 3x10^9 M_sun). Herschel and SCUBA-2 surveys should be able to detect hundreds of these power-law dominated DOGs. We use HST and Spitzer/IRAC data to estimate stellar masses of these sources and find that the stellar to gas mass ratio may be higher in our 24um-bright sample of DOGs than in SMGs and other Spitzer-selected sources. Although larger sample sizes are needed to provide a definitive conclusion, the data are consistent with an evolutionary trend in which the formation of massive galaxies at z~2 involves a sub-millimeter bright, cold-dust and star-formation dominated phase followed by a 24um-bright, warm-dust and AGN-dominated phase.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

We present a 1200-micron image of the Great Observatories Origin Deep Survey North (GOODS-N) field, obtained with the Max Planck Millimeter Bolometer array (MAMBO) on the IRAM 30-m telescope. The survey covers a contiguous area of 287 square arcmin to a near-uniform noise level of ~0.7mJy/beam. After Bayesian flux deboosting, a total of 30 sources are recovered (>=3.5sigma). An optimal combination of our 1200-micron data and an existing 850-micron image from the Submillimetre Common-User Bolometer Array (SCUBA) yielded 33 sources (>=4sigma). We combine our GOODS-N sample with those obtained in the Lockman Hole and ELAIS-N2 fields (Scott et al. 2002; Greve et al. 2004) in order to explore the degree of overlap between 1200-micron- and 850-micron-selected galaxies (hereafter SMGs), finding no significant difference between their 850-micron to 1200-micron flux density distributions. However, a noise-weighted stacking analysis yields a significant detection of the 1200-micron-blank SCUBA sources, whereas no significant 850-micron signal is found for the 850-micron-blank MAMBO sources. The hypothesis that the 850/1200-micron flux density distribution of SCUBA sources is also representative of the MAMBO population is rejected at the ~4sigma level, via Monte Carlo simulations. Therefore, although the populations overlap, galaxies selected at 850 and 1200micron are different, and there is compelling evidence for a significant 1200-micron-detected population which is not recovered at 850micron. These are submm drop-outs (SDOs), with S_850/S_1200 = 0.7-1.7, requiring very cold dust or unusual spectral energy distributions (T_d ~ 10K; beta ~ 1), unless SDOs reside beyond the redshift range observed for radio-identified SMGs, i.e. at z > 4.

We present a comparison between the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 and 850 microns in the Lockman Hole East with a deep Spitzer Space Telescope survey at 3.6-24 microns conducted in Guaranteed Time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer-identified submm galaxies at 1<z<1.5 appear to reside in overdensities when the density field is smoothed at 0.5-2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation - galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses <10^11 Msun, this peaks at lower redshifts than the ostensible z~2.2 maximum for submm point sources, adding to the growing consensus for ``downsizing'' in star formation. Our surveys are also consistent with ``downsizing'' in mass assembly. Both the mean star formation rates <dM/dt> and specific star formation rates <(1/M)dM/dt> are in striking disagreement with some semi-analytic predictions from the Millenium simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

We present Spitzer/IRS mid-infrared spectra for 15 gravitationally lensed, 24 micron--selected galaxies, and combine the results with 4 additional very faint galaxies with IRS spectra in the literature. The median intrinsic 24 micron flux density of the sample is 130 microJy, enabling a systematic survey of the spectral properties of the very faint 24 micron sources that dominate the number counts of Spitzer cosmological surveys. Six of the 19 galaxy spectra (32%) show the strong mid-IR continuua expected of AGN; X-ray detections confirm the presence of AGN in three of these cases, and reveal AGNs in two other galaxies. These results suggest that nuclear accretion may contribute more flux to faint 24 micron--selected samples than previously assumed. Almost all the spectra show some aromatic (PAH) emission features; the measured aromatic flux ratios do not show evolution from z=0. In particular, the high S/N mid-IR spectrum of SMM J163554.2+661225 agrees remarkably well with low--redshift, lower--luminosity templates. We compare the rest-frame 8 micron and total infrared luminosities of star--forming galaxies, and find that the behavior of this ratio with total IR luminosity has evolved modestly from z=2 to z=0. Since the high aromatic--to--continuum flux ratios in these galaxies rule out a dominant contribution by AGN, this finding implies systematic evolution in the structure and/or metallicity of infrared sources with redshift. It also has implications for the estimates of star forming rates inferred from 24 micron measurements, in the sense that at z ~2, a given observed frame 24 micron luminosity corresponds to a lower bolometric luminosity than would be inferred from low-redshift templates of similar luminosity at the corresponding rest wavelength.

We present mid/far-infrared photometry of nine FeLoBAL QSOs, taken using the Spitzer space telescope. All nine objects are extremely bright in the infrared, with rest-frame 1-1000 micron luminosities comparable to those of Ultraluminous Infrared Galaxies. Furthermore, a significant fraction of the infrared emission from many, and possibly all of the sample is likely to arise from star formation, with star formation rates of order several hundred solar masses per year. We combine these results with previous work to propose that FeLoBALs mark galaxies and QSOs in which an extremely luminous starburst is approaching its end, and in which a rapidly accreting supermassive black hole is in the last stages of casting off its dust cocoon. FeLoBAL signatures in high redshift QSOs and galaxies may thus be an efficient way of selecting sources at a critical point in their evolution.

We present measurements of the spatial clustering of ultraluminous infrared galaxies in two redshift intervals, 1.5<z<2.0 and 2<z<3. Both samples cluster strongly, with r_0=14.40+/-1.99 h^-1 Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1 Mpc for the 1.5<z<2.0 sample, making them among the most biased galaxies at these epochs. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13 solar masses. We infer that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0.

In a search for evidence of the short wavelength increment in the Sunyaev-Zel'dovich (SZ) effect, we have analyzed archival galaxy cluster data from the Sub-millimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope, resulting in the most complete pointed survey of clusters at 850 microns to date. SCUBA's 850 microns passband overlaps the peak of the SZ increment. The sample consists of 44 galaxy clusters in the range 0 < z < 1.3. Maps of each of the clusters have been made and sources have been extracted; as an ancillary product we generate the most thorough galaxy cluster point source list yet from SCUBA. Seventeen of these clusters are free of obvious AGN and have data deep enough to provide interesting measurements of the expected SZ signal. Specialized analysis techniques are employed to extract the SZ effect signal from these SCUBA data, including using SCUBA's short wavelength band as an atmospheric monitor and fitting the long wavelength channel to a model of the spatial distribution of each cluster's SZ effect. By explicitly excising the exact cluster centre from our analysis we demonstrate that emission from galaxies within the cluster does not contaminate our measurement. The SZ amplitudes from our measurements are consistently higher than the amplitudes inferred from low frequency measurements of the SZ decrement.

We have imaged a $\sim$6 arcminute$^2$ region in the Bo\"otes Deep Field using the 350 $\mu$m-optimised second generation Submillimeter High Angular Resolution Camera (SHARC II), achieving a peak 1$\sigma$ sensitivity of $\sim$5 mJy. We detect three sources above 3$\sigma$, and determine a spurious source detection rate of 1.09 in our maps. In the absence of $5\sigma$ detections, we rely on deep 24 $\mu$m and 20 cm imaging to deduce which sources are most likely to be genuine, giving two real sources. From this we derive an integral source count of 0.84$^{+1.39}_{-0.61}$ sources arcmin$^{-2}$ at $S>13$ mJy, which is consistent with 350 $\mu$m source count models that have an IR-luminous galaxy population evolving with redshift. We use these constraints to consider the future for ground-based short-submillimetre surveys.

We present 350micron observations of 36 ultraluminous infrared galaxies (ULIRGs) at intermediate redshifts (0.089 <= z <= 0.926) using the Submillimeter High Angular Resolution Camera II (SHARC-II) on the Caltech Submillimeter Observatory (CSO). In total, 28 sources are detected at S/N >= 3, providing the first flux measurements longward of 100micron for a statistically significant sample of ULIRGs in the redshift range of 0.1 < z < 1.0. Combining our 350micron flux measurements with the existing IRAS 60 and 100micron data, we fit a single-temperature model to the spectral energy distribution (SED), and thereby estimate dust temperatures and far-IR luminosities. Assuming an emissivity index of beta = 1.5, we find a median dust temperature and far-IR luminosity of Td = 42.8+-7.1K and log(Lfir/Lsolar) = 12.2+-0.5, respectively. The far-IR/radio correlation observed in local star-forming galaxies is found to hold for ULIRGs in the redshift range 0.1 < z < 0.5, suggesting that the dust in these sources is predominantly heated by starbursts. We compare the far-IR luminosities and dust temperatures derived for dusty galaxy samples at low and high redshifts with our sample of ULIRGs at intermediate redshift. A general Lfir-Td relation is observed, albeit with significant scatter, due to differing selection effects and variations in dust mass and grain properties. The relatively high dust temperatures observed for our sample compared to that of high-z submillimeter-selected starbursts with similar far-IR luminosities suggest that the dominant star formation in ULIRGs at moderate redshifts takes place on smaller spatial scales than at higher redshifts.

We report the results of pointed observations of the prototypical ultra-luminous infrared galaxy (ULIRG) Arp 220 at 850 microns using the polarimeter on the SCUBA instrument on the James Clerk Maxwell Telescope. We find a Bayesian 99 per cent confidence upper limit on the polarized emission for Arp 220 of 1.54 per cent, averaged over the 15 arcsec beam-size. Arp 220 can serve as a proxy for other, more distant such galaxies. This upper limit constrains the magnetic field geometry in Arp 220 and also provides evidence that polarized ULIRGs will not be a major contaminant for next-generation cosmic microwave background polarization measurements.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

We present two-dimensional spectroscopy covering the rest-frame wavelengths of strong optical emission lines in six luminous submillimetre galaxies at z=1.3--2.5. Together with HST ACS and NICMOS imaging we map the dynamics and morphologies of these systems on scales from 4--11kpc. Four of the systems show multiple components in their spatially-resolved spectra with average velocity offsets of ~180km/s across 8kpc in projection. From the ensemble properties from our survey and the literature, we estimate the typical dynamical masses of bright sub-mm galaxies as (5+/-3) x10^{11}Mo. This is similar to recent estimates of their stellar masses -- suggesting that the dynamics of the central regions of these galaxies are baryon dominated, with a substantial fraction of those baryons in stars by the epoch of observation. Combining our dynamical mass estimates with stellar luminosities for this population we investigate whether submillimetre galaxies can evolve onto the Faber-Jackson relation for local ellipticals. Adopting a typical lifetime of tau_{burst}~300Myr for the submillimetre-luminous phase we find that the stellar populations of sub-mm galaxies should fade to place them on the Faber-Jackson relation, at M_K ~ -25.1. Furthermore, using the same starburst lifetime we correct the observed space density of submillimetre galaxies for the duty cycle to derive a volume density of the progenitors of ~1x10^{-4} Mpc^{-3}. This is consistent with the space density of local luminous early-type galaxies with M_K ~ -25.1, indicating that submillimetre galaxies can evolve onto the scaling relations observed for local early-type galaxies, and the observed population at z~2 is then sufficient to account for the formation of the whole population of >~3 L^{*} ellipticals seen at z~0.

We present SEDs, Spitzer colours, and IR luminosities for 850 micron selected galaxies in the GOODS-N field. Using the deep Spitzer Legacy images and new data and reductions of the VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submm sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. We find that in most cases the 850 micron emission is dominated by a single 24 micron source. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ULIRGs of the same luminosity and therefore appear to be cooler. The SEDs of submm galaxies are also different from those of their high redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Using 24 micron. 850 micron and 1.4 GHz observations, we fit templates that span the mid-IR through radio to derived the integrated IR luminosity of the submm galaxies. By themselves, 24 micron and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model independent estimate of the redshift. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts show strong evidence for an active galactic nucleus dominating the IR luminosity.

We present measurements of the spatial clustering of galaxies with stellar masses >10^11Msun, infrared luminosities >10^12 Lsun, and star formation rates >200Msun per year in two redshift intervals; 1.5<z<2.0 and 2<z<3. Both samples cluster very strongly, with spatial correlation lengths of r_0=14.40+/-1.99 h^-1Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1Mpc for the 1.5<z<2.0 sample. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13Msun, making them among the most biased galaxies at these epochs. We infer, from this and previous results, that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1, both starbursts and AGN. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0. We conclude that ULIRGs at z>1 signpost stellar buildup in galaxies that will reside in clusters at z=0, with ULIRGs at increasing redshifts signposting the buildup of stars in galaxies that will reside in increasingly rich clusters.

How does the submm galaxy population detected by SCUBA fit into galaxy evolution? How do these rare starbursting systems, which contribute significantly to high redshift star-formation, relate to other galaxy populations? Deep radio observations have been most useful for studying these systems, but still leave a significant fraction of the population unidentified. Now with the deep IRAC and MIPS images from the GOODS Spitzer Legacy program and a re-analysis of the deep radio data, we are able to identify counterparts for a large fraction of SCUBA galaxies in GOODS-N. All of these counterparts are detected by Spitzer. Given the vast multi-wavelength data in this field, we can study the spectral energy distributions (SEDs) of these systems and determine what is fueling their intense infrared luminosities. A rest-frame composite optical-to-radio SED for all spectroscopically identified submm sources shows that the average SCUBA galaxy is consistent with models of ultraluminous starburst galaxies, although cooler than those observed locally. Because of this, the submm flux alone consistently overestimates IR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The wide range of 24/850 micron flux ratios as a function of redshift indicates the presence of strong mid-IR features, to be confirmed with deep IRS spectroscopy. The IRAC colours of the submm systems provide useful redshift constraints, since, at these redshifts, IRAC samples the stellar bump. The Spitzer photometry of this large sample of submm galaxies has allowed us to put constraints on many of the outstanding issues in submm astronomy.

We present high resolution 890 micron images of two 20 mJy submillimeter galaxies, SMMJ123711+622212 and MIPS~J142824.0+352619, obtained using the Submillimeter Array (SMA). Using submillimeter interferometric observations with an angular resolution of 2.5", the coordinates of these high redshift sources are determined with an accuracy of 0.2". The new SMA data on SMMJ123711+622212 reveal an unresolved submm source offset to the east by 0.8" from an optical galaxy found in deep HST images, suggesting either a large galaxy with a dusty central region, or an interacting galaxy system. The SMA image of hyper-luminous (LFIR = 3.2 x 10^{13} Lsun) source MIPS~J142824.0+352619 provides a firm upper limit to the source size of < 1.2. This constraint provides evidence that the foreground lens is only weakly affecting the observed high FIR luminosity.

We present low-resolution (64 < R < 124) mid-infrared (8--38 micron) Spitzer/IRS spectra of two z~1.3 ultraluminous infrared galaxies (LFIR~10^13) discovered in a Spitzer/MIPS survey of the Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). MIPS J142824.0+352619 is a bright 160 micron source with a large infrared-to-optical flux density ratio and a possible lensing amplification of <~10. The 6.2, 7.7, 11.3, and 12.8 micron PAH emission bands in its IRS spectrum indicate a redshift of z~1.3. The large equivalent width of the 6.2 micron PAH feature indicates that at least 50% of the mid-infrared energy is generated in a starburst, an interpretation that is supported by a large [NeII]/[NeIII] ratio and a low upper limit on the X-ray luminosity. SST24 J142827.19+354127.71 has the brightest 24 micron flux (10.55 mJy) among optically faint (R > 20) galaxies in the NDWFS. Its mid-infrared spectrum lacks emission features, but the broad 9.7 micron silicate absorption band places this source at z~1.3. Given this redshift, SST24 J142827.19+354127.71 has among the largest rest-frame 5 micron luminosities known. The similarity of its SED to those of known AGN-dominated ULIRGs and its lack of either PAH features or large amounts of cool dust indicate that the powerful mid-infrared emission is dominated by an active nucleus rather than a starburst. Our results illustrate the power of the IRS in identifying massive galaxies in the ``redshift desert'' and in discerning their power sources. Because they are bright, MIPS J142824.0+352619 (pending future observations to constrain its lensing amplification) and SST24 J142827.19+354127.71 are useful z>1 templates of a high luminosity starburst and AGN, respectively.

Using the SHARC-II camera at the Caltech Submillimeter Observatory to obtain 350 micron images of sources detected with the MIPS instrument on Spitzer, we have discovered a remarkable object at z=1.325+/-0.002 with an apparent Far-Infrared luminosity of 3.2(+/-0.7) x 10^13 Lsun. Unlike other z>1 sources of comparable luminosity selected from mid-IR surveys, MIPS J142824.0+352619 lacks any trace of AGN activity, and is likely a luminous analog of galaxies selected locally by IRAS, or at high redshift in the submillimeter. This source appears to be lensed by a foreground elliptical galaxy at z=1.034, although the amplification is likely modest (~10). We argue that the contribution to the observed optical/Near-IR emission from the foreground galaxy is small, and hence are able to present the rest-frame UV through radio Spectral Energy Distribution of this galaxy. Due to its unusually high luminosity, MIPS J142824.0+352619 presents a unique chance to study a high redshift dusty starburst galaxy in great detail.

We analyze deep X-ray, optical and mid-infrared Spitzer observations of the CDF-N/GOODS-N region to study 13 submillimeter-detected galaxies (SMGs) with spectroscopic redshifts (median z=2.2). We find a correlation between the estimated stellar and X-ray luminosity, implying that masses of the black holes may be related to the stellar masses of their host galaxies. Although the submillimeter emission implies that these galaxies are undergoing an epoch of intense star-formation, the Spitzer data reveal a massive stellar population already in place. These stellar masses are then compared to previously published black hole mass estimates derived from the X-ray luminosities under the assumption of Eddington-limit accretion. We find that the black hole masses for our high-redshift sample are approximately 1-2 orders of magnitude smaller than galaxies of comparable stellar mass in the local Universe. Although our estimates of black hole masses will increase if the accretion is sub-Eddington, and our stellar masses will decrease if we assume a much younger stellar population or a different initial mass function, we find that only through a combination of effects is it possible to shift the high redshift galaxies such that they lie on the local relation. This suggests that the black holes need to grow substantially between z=2.2 and the present-day, with much of the black hole growth occurring after the current obscured, far-infrared luminous phase of activity which is likely associated with the formation of the spheroid. This interpretation supports a scenario where SMGs pass through a subsequent accretion-dominated phase, where they would appear as optically bright quasars. (ABRIDGED)

Low-ionization-state far-infrared (FIR) emission lines may be useful diagnostics of star-formation activity in young galaxies, and at high redshift may be detectable from the ground. In practice, however, very little is known concerning how strong such line emission might be in the early Universe. We attempted to detect the 158 micron [C II] line from a lensed galaxy at z = 4.926 using the Caltech Submillimeter Observatory. This source is an ordinary galaxy, in the sense that it shows high but not extreme star formation, but lensing makes it visible. Our analysis includes a careful consideration of the calibrations and weighting of the individual scans. We find only modest improvement over the simpler reduction methods, however, and the final spectrum remains dominated by systematic baseline ripple effects. We obtain a 95 per cent confidence upper limit of 33 mJy for a 200 km/s full width at half maximum line, corresponding to an unlensed luminosity of 1x10^9 L_sun for a standard cosmology. Combining this with a marginal detection of the continuum emission using the James Clerk Maxwell Telescope, we derive an upper limit of 0.4 per cent for the ratio of L_CII/L_FIR in this object.

We present a new sub-mm Super-map in the HDF-North region (GOODS-North field), containing 40 statistically robust sources at 850 microns. This map contains additional data, and several new sources, including one of the brightest blank-sky extragalactic sub-mm sources ever detected. We have used the ACS HST images and ground-based near-IR observations from GOODS to develop a systematic approach for counterpart identification. 72 per cent of our sources with optical coverage have a unique optical counterpart using our new techniques for counterpart identification, and an additional 18 per cent have more than one possibility that meet our criteria in the ACS images. We have found a much higher ERO rate than other sub-mm surveys, due to the increased depth in the optical images. The median photometric redshift (and quartile range), from optical and near-infrared data, is 1.7 (1.3-2.5) for the radio-detected sub-mm sources, and rises to 2.3 (1.3-2.7) for the radio-undetected sub-sample. We find interesting correlations between the 850 micron flux and both the i magnitude and the photometric redshift, from which there appears to be an absence of high redshift faint counterparts to the lower flux density SCUBA sources. While the quantitative morphologies span a range of values, in general the sub-mm galaxies show larger sizes and a higher degree of asymmetry than other galaxy populations at the same redshifts.

We present optical and near-infrared integral field spectroscopy of the SCUBA galaxy SMM J163650.43+405734.5 (ELAIS N2 850.4) at z=2.385. We combine Ly-alpha and H-alpha emission line maps and velocity structure with high resolution HST ACS and NICMOS imaging to probe the complex dynamics of this vigorous star-burst galaxy. The imaging data shows a complex morphology, consisting of at least three components separated by ~1'' (8kpc) in projection. When combined with the H-alpha velocity field from UKIRT UIST IFU observations we identify two components whose redshifts are coincident with the systemic redshift, measured from previous CO observations, one of which shows signs of AGN activity. A third component is offset by 220+/-50km/s from the systemic velocity. The total star formation rate of the whole system (estimated from the narrow-line H-alpha and uncorrected for reddening) is 340+/-50Mo/yr. The Ly-alpha emission mapped by the GMOS IFU covers the complete galaxy and is offset by +270+/-40km/s from the systemic velocity. This velocity offset is comparable to that seen in rest-frame UV-selected galaxies at similar redshifts and usually interpreted as a star-burst driven wind. The extended structure of the Ly-alpha emission suggests that this wind is not a nuclear phenomenon, but is instead a galactic scale outflow. Our observations suggest that the vigorous activity in N2 850.4 is arising as a result of an interaction between at least two dynamically-distinct components, resulting in a strong starburst, a starburst-driven wind and actively-fuelled AGN activity. [abridged]

We present an 850-micron map and list of candidate sources in a 70 arcmin^2 sub-area of the Groth Strip observed using SCUBA. We initially detect 7 candidate sources with signal-to-noise ratios (SNRs) between 3.0 and 3.5 and 4 candidate sources with SNR > 3.5. Simulations suggest that on average in a map this size one expects 1.6 false positive sources for SNR > 3.5 and 4.5 for 3 < SNR < 3.5. Flux boosting in maps is a well known effect and we have developed a simple Bayesian prescription for estimating the unboosted flux distribution and used this method to determine the best flux estimates of our sources. This method is easily adapted for any other modest signal-to-noise survey in which there is prior knowledge of the source counts. We performed follow-up photometry in an attempt to confirm or reject 5 of our source candidates. We confirm the reality of 2 of the SCUBA sources, although at lower levels than suggested in the map and we find that the photometry results are consistent with and confirm the de-boosted map fluxes. Our final candidate source list contains 3 sources, including the 2 confirmed detections and 1 further candidate source with SNR > 3.5 which has a reasonable chance of being real.

We investigate the photometric redshift distribution, colours and morphologies of a large sample of submillimetre sources in the GOODS-North region. With the depth achieved by the deep ACS images, optical counterparts have been found for all the radio-detected sub-mm sources, and we have used several techniques to identify counterparts to the radio-undetected sources. The sub-mm counterparts have a median photometric redshift, from optical and near-IR photometry, of 1.9. Our sample contains a much higher fraction of EROs than other sub-mm surveys, and larger angular sizes than the average galaxy at all redshifts. The sub-mm galaxy counterparts are consistent with massive star-forming galaxies at high redshift.

We present radio, optical and X-ray detected counterparts to the sub-mm sources found using SCUBA in the Hubble Deep Field North region (GOODS-N). A new counterpart identification statistic is developed to identify properties of galaxies detected at other wavelengths that can be used to aid counterpart identification. We discriminate between criteria that can be used to pre-select sub-mm bright objects, and those that identify the counterpart to a known sub-mm object. Optically faint galaxies detected in the deepest 1.4 GHz radio continuum maps are the only effective way of pre-selecting SCUBA galaxies, and radio sources are the best way to identify counterparts to known sub-mm detections. Looking at radio spectral indices, only the steeper sources (indicative of star formation) are detected in the sub-mm. Although we find several X-ray identifications, we show that deep Chandra images do not contribute to counterpart identifications, since in all cases they are already detected in the more easily obtained VLA radio maps. We also find find no evidence for clustering between Chandra and SCUBA sources in this field. For a known SCUBA position, the reddest source tends to be the correct association, although we can find no cut on colour, magnitude, or clustering property that efficiently pre-selects for SCUBA sources. 15 micron ISO sources are statistically detected by SCUBA, but the limiting mid-IR flux is not low enough to provide useful constraints. We present postage stamp strips for each SCUBA detection in separate bands from X-ray to radio, providing direct visual evidence that approximately half of the sub-mm sources in this field remain unidentified, despite an abundance of deep multi-wavelength data.

We present a new SCUBA image of the cluster MS0451.6-0305, which exhibits strong, extended sub-mm flux at 850 microns. The most striking feature in the map is an elongated region of bright sub-mm emission, with a flux density of >10 mJy over several beam-sizes. This region is apparently coincident with a previously known optical arc (which turns out to be a strongly lensed Lyman Break Galaxy at z=2.911), as well as with a newly identified multiply-imaged ERO (Extremely Red Object) pair predicted to be at a similar, if not identical redshift. By combing a detailed lensing model with deep images from HST, Chandra, CFHT, JCMT, and spectra from the VLT, we conclude that both the strongly lensed optical arc and ERO systems have properties consistent with known sub-mm emitters. Using a simple model for the two sources, we estimate that the multiply-lensed EROs contribute the majority of the flux in the SCUBA lensed arc. Correcting for the lensing amplification, we estimate that the inherent 850 micron fluxes for both objects are ~0.4 mJy. If the LBG and ERO pair are truly at the same redshift, then they are separated by only ~10 kpc in the source plane, and hence constitute an interacting system at z~2.9. Higher angular resolution observations in sub-mm/mm will permit us to more accurately separate the contribution from each candidate, and better understand the nature of this system.

The cosmologically significant population of dusty galaxies detected at sub-mm and mm wavelengths by SCUBA and MAMBO have played a central part in many of the discussions at this workshop. Obtaining their redshifts is critical in understanding their role in galaxy formation and evolution, but this is a challenging task. Significant progress has been made by obtaining spectroscopic redshifts of the subset of sub-mm galaxies detected in deep radio images, but other groups advocate a photometric redshift approach. Such methods rely on understanding the Spectral Energy Distribution of the galaxies under study, but these are not well characterized for even the well-studied nearby Infrared luminous population. Observations at short sub-mm wavelengths, previously difficult until the commissioning of SHARC-II at the CSO, can be used to better constrain the dust properties of these galaxies.

We present new follow-up observations of the sub-mm luminous Ly_alpha-emitting object in the SSA22 z=3.09 galaxy overdensity, referred to as `Blob 1' by Steidel et al.(2000). In particular we discuss high resolution Hubble Space Telescope optical imaging, Owens Valley Radio Observatory spectral imaging, Keck spectroscopy, VLA 20cm radio continuum imaging, and Chandra X-ray observations. We also present a more complete analysis of the existing James Clerk Maxwell Telescope sub-mm data. We detect several optical continuum components which may be associated with the core of the submillimeter emitting region. A radio source at the position of one of the HST components (22:17:25.94, +00:12:38.9) identifies it as the likely counterpart to the submillimeter source. We also tentatively detect the CO(4-3) molecular line, centered on the radio position. We use the CO(4-3) intensity to estimate a limit on the gas mass for the system. The optical morphology of sources within the Ly_alpha cloud appears to be filamentary, while the optical source identified with the radio source has a dense knot which may be an AGN or compact starburst. We obtain a Keck-LRIS spectrum of this object, despite its faintness (R=26.8). The spectrum reveals weak Ly_alpha emission, but no other obvious features, suggesting that the source is not an energetic AGN (or that it is extremely obscured). We use non-detections in deep Chandra X-ray images to constrain the nature of the `Blob'. Although conclusive evidence regarding the nature of the object remains hard to obtain at this redshift, the evidence presented here is at least consistent with a dust-obscured AGN surrounded by a starburst situated at the heart of this giant Ly_alpha cloud.

We have detected the Sunyaev-Zel'dovich (SZ) increment at 850 microns in two galaxy clusters (Cl 0016+16 and MS 1054.4-0321) using SCUBA (Sub-millimetre Common User Bolometer Array) on the James Clerk Maxwell Telescope. Fits to the isothermal \beta model yield a central Compton y parameter of (2.2 \pm 0.7) x 10^{-4} and a central 850 micron flux of \Delta I_{0} = 2.2 \pm 0.7 mJy/beam in Cl 0016. This can be combined with decrement measurements to infer y = (2.38 \pm_{0.34}^{0.36}) x 10^{-4} and v_{pec} = 400 \pm_{1400}^{1900} km/s. In MS 1054 we find a peak 850 micron flux of \Delta I_{0} = 2.0 \pm 1.0 mJy/beam and y = (2.0 \pm 1.0) x 10^{-4}. To be successful such measurements require large chop throws and non-standard data analysis techniques. In particular, the 450 micron data are used to remove atmospheric variations in the 850 micron data. An explicit annular model is fit to the SCUBA difference data in order to extract the radial profile, and separately fit to the model differences to minimize the effect of correlations induced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can be used to measure the SZ increment in massive, compact galaxy clusters.

We investigate the emission of sub-millimetre-wave radiation from galaxies in a 165 square arcminute region surrounding the Hubble Deep Field North. The data were obtained from dedicated observing runs from our group and others using the SCUBA camera on the James Clerk Maxwell Telescope, and combined using techniques specifically developed for low signal-to-noise source recovery. The resulting `Super-map' is derived from about 60 shifts of JCMT time, taken in a variety of observing modes and chopping strategies, and combined here for the first time. At 850 micron we detect 19 sources at >4 sigma, including 5 not previously reported. We also list an additional 15 sources between 3.5 and 4.0 sigma (where 2 are expected by chance). The 450 micron map contains 5 sources at >4 sigma. We present a new estimate of the 850 micron and 450 micron source counts. The number of sub-mm galaxies we detect account for approximately 40% of the 850 micron sub-mm background, and we show that mild extrapolations can reproduce it entirely. A clustering analysis fails to detect any significant signal in this sample of SCUBA detected objects. A companion paper describes the multiwavelength properties of the sources.

We present results from JCMT sub-mm observations of sources selected from the {\sl ISO} FIRBACK (Far-IR BACKground) survey, along with UKIRT near-IR imaging of a sub-sample. This gives valuable insight into the brightest $\sim$10% of galaxies which contribute to the Cosmic Infrared Background (CIB). We estimate the photometric redshifts and luminosities of these sources by fitting their Spectral Energy Distributions (SEDs). The data appear to show a bimodal galaxy distribution, with normal star-forming galaxies at $z\simeq0$, and a much more luminous population at $z\sim0.4$--0.9. These are similar to the ultraluminous infrared galaxies which are found to evolve rapidly with redshift in other surveys. The detectability threshold of FIRBACK biases the sample away from much higher redshift ($z\stackrel{>}{_{\sim}}1.5$) objects. Nevertheless, the handful of $z\sim0.5$ sources which we identify are likely to be the low-$z$ counterparts of the typically higher-$z$ sources found in blank field sub-mm observations. This sub-sample, being much more nearby than the average SCUBA galaxies, has the virtue of being relatively easy to study in the optical. Hence their detailed investigation could help elucidate the nature of the sub-mm bright galaxies.

(abridged) We present multiwavelength observations for a large sample of microjansky radio sources detected in ultradeep 1.4GHz maps centered on the Hubble Deep Field-North (HDF-N) and the Hawaii Survey Fields SSA13 and SSA22. Our spectroscopic redshifts for 169 radio sources reveal a flat median redshift distribution, and these sources are hosted by similarly luminous optical L* galaxies, regardless of redshift. This is a serious concern for radio estimates of the local star formation rate density, as a substantial fraction of the ultraviolet luminosity density is generated by sub-L* galaxies at low redshifts. From our submillimeter measurements for 278 radio sources, we find error-weighted mean 850micron fluxes of 1.72$\pm$0.09 mJy for the total sample, 2.37$\pm$0.13 mJy for the optically-faint (I>23.5) subsample, and 1.04$\pm$0.13 mJy for the optically-bright (I<23.5) subsample. We significantly (>3\sigma) detect in the submillimeter 50 of the radio sources, 38 with I>23.5. Spectroscopic redshifts for three of the I<23.5 submillimeter-detected radio sources are in the range z=1.0-3.4, and all show AGN signatures. Using only the submillimeter mapped regions we find that 69\pm9% of the submillimeter-detected radio population are at I>23.5. We also find that 66\pm7% of the S850>5 mJy (>4\sigma) sources are radio-identified. We find that millimetric redshift estimates at low redshifts are best made with a FIR template intermediate between a Milky Way type galaxy and a starburst galaxy, and at high redshifts with an Arp220 template.

Using the 850 micron SCUBA camera on the JCMT and a scanning technique different from other sub-mm surveys, we have obtained a 125 square arcminute map centered on the Hubble Deep Field. The one-sigma sensitivity to point sources is roughly 3 mJy and thus our map probes the brighter end of the sub-mm source counts. We find 6 sources with a flux greater than about 12 mJy (>4 sigma) and, after a careful accounting of incompleteness and flux bias, estimate the integrated density of bright sources N(>12 mJy)= 164 (+77/-58) per square degree (68 per cent confidence bounds).

We present a sample of 43 submillimeter sources detected (at >3 sigma), drawn from our program to follow-up optically faint radio sources with SCUBA. These sources already have associated radio and in many cases optical identifications, and many are also detected at 450 microns. We compare these with 12 submillimeter sources drawn from the literature, which were discovered in blank field mapping campaigns, but also have radio detections. We then use this total sample (55 sources) to study and model the evolution of dusty galaxies. A correlation is observed in the sub-mm/radio color-magnitude diagram, which can be modeled by strong luminosity evolution. The selection effects of the radio/optical pre-selection technique are determined from the models, and a corrected redshift distribution is constrained using a range of model assumptions. The temperature/redshift effects on the 450 microns detected subset of our sample are studied in relation to the models, and prospects for improved measurements in the shorter sub-mm wavelength windows (450 microns and 350 microns) are explored.

Summary and discussion of some projects to use SCUBA to target sources selected at other wavebands, as well as to find new sub-mm galaxies in `blank fields': FIRBACK galaxies; Lyman break galaxies and `the Blob'; HDF flanking fields and the Groth Strip; survey of lensing cluster fields.

We have used the Submillimetre Common-User Bolometer Array (SCUBA) detector on the James Clerk Maxwell Telescope (JCMT) to measure bright sub-mm emission associated with a recently discovered extensive (>100/h kpc) and highly luminous, `blob' of Ly-alpha emission at z=3.09. The blob lies within a known large overdensity of optical sources in the z=3.07-3.11 range, and is centered on a locally overdense peak within this region. The best explanation for the copious sub-mm emission is a dust obscured continuum source, which may produce the ionizing flux for the Ly-alpha cloud. Cooling gas explanations are plausible but excessively complicated, and the 450/850 micron ratio rules out a significant fraction of the signal arising from the Sunyaev-Zel'dovich increment. At least two additional ~10 mJy sub-mm detections in the SCUBA map, with a surface density significantly higher than in blank field surveys, suggests that they may be associated with the z=3.09 structure. A SCUBA `photometry' observation of a second nearby Ly-alpha blob tentatively detects a weaker sub-mm counterpart.

We summarize the main results from our SCUBA survey of Lyman-break galaxies (LBGs) at z~3. Analysis of our sample of LBGs reveals a mean flux of S850=0.6$\pm$0.2 mJy, while simple models of emission based on the UV properties predict a mean flux about twice as large. Known populations of LBGs are expected to contribute flux to the weak sub-mm source portion of the far-IR background, but are not likely to comprise the bright source (S850>5 mJy) end of the SCUBA-detected source count. The detection of the LBG, Westphal-MM8, at 1.9 mJy suggests that deeper observations of individual LBGs in our sample could uncover detections at similar levels, consistent with our UV-based predictions. By the same token, many sub-mm selected sources with S850<2 mJy could be LBGs. The data are also consistent with the FarIR/$\beta$ relation holding at z=3.

We present an 11 X 11 arcminute map centred on the Hubble Deep Field taken at 850 microns with the SCUBA camera on the JCMT. The map has an average one-sigma sensitivity to point sources of about 2.3 mJy and thus probes the brighter end of the sub-mm source counts. We find 7 sources with a flux greater than 9 mJy (roughly 4 sigma), and therefore estimate N(>9 mJy)= 208 (+90/-72) per degree. This result is consistent with work from other groups, but improves the statistics at the bright end, and is suggestive of a steepening of the counts.

Recent submillimetre surveys have revealed a population of dusty, high redshift sources of great cosmological significance for understanding dust-enshrouded star formation in distant galaxies, and for determining the origin of the far-IR background. In this paper, we analyze nine rich cluster fields mapped at 850 and 450 microns with the SCUBA array on the James Clerk Maxwell telescope. Lensing models of the clusters are developed in order to derive accurate source counts for our sample. VLA maps of the same clusters are used to help constrain the redshift distribution of our SCUBA detections. Implications for high redshift galaxies and for the far-IR background are discussed. We also provide limits on distributed dust produced by cooling flows in these clusters.

We have carried out targetted sub-mm observations as part of a programme to explore the connection between the rest-frame UV and far-IR properties of star-forming galaxies at high redshift, which is currently poorly understood. On the one hand the Lyman break technique is very effective at selecting z~3 galaxies. On the other hand `blank field' imaging in the sub-mm seems to turn up sources routinely, amongst which some are star forming galaxies at similar redshifts. Already much work has been done searching for optical identifications of objects detected using the SCUBA instrument. Here we have taken the opposite approach, performing sub-mm photometry for a sample of Lyman break galaxies whose UV properties imply high star formation rates. The total signal from our Lyman break sample is undetected in the sub-mm, at an RMS level of 0.5 mJy, which implies that the population of Lyman break galaxies does not constitute a large part of the recently detected blank-field sub-mm sources. However, our one detection suggests that with reasonable SCUBA integrations we might expect to detect those few LBGs that are far-IR brightest.

The FIRBACK (Far InfraRed BACKground) survey represents the deepest extensive 170 micron images obtained by the ISO satellite. The sources detected comprise about 10% of the Cosmic IR Background (CIB) seen by COBE, and, importantly, were observed at a wavelength near the peak of the CIB. Detailed follow-up of these sources should help pin down the redshifts, Hubble types, and other properties of the galaxies which constitute this background. We have used the Submillimetre Common-User Bolometer Array (SCUBA) instrument on the James Clerk Maxwell Telescope (JCMT) to search for sub-mm emission from a sample of 10 galaxies which are bright at 170 microns and for which we had accurate radio positions. Statistically we detect this sample at both 450 microns (~4 sigma) and 850 microns (~ 7 sigma); individual objects are convincingly detected in four cases at 850 microns, and one case at 450 microns. Fits to simple spectral energy distributions suggest a range of low to moderate redshifts, perhaps z=0 to 1.5 for this FIRBACK sub-sample.

We present a detailed mm-wave and optical study of the gravitational lens system B1933+503, discovered by Sykes et al.(1998) in the radio. This object is probably the most complex lens system known, with 10 lensed components within a radius of one arcsecond. It is potentially important as a probe of the Hubble constant, although no optical counterpart has thus far been observed down to I=24.2. We have obtained new sub-millimetre detections at 450, 850 and 1350 microns. We have also constrained the possible dust emission from the proposed foreground lensing galaxy using a K-band adaptive optics image and CO(5-4) measurements. A lensing model is constructed, taking the foreground elliptical galaxy at z=0.755 as the lensing mass. From this we derive a scenario from which to model the sub-millimetre emission. Several arguments then point to the source in the B1933+503 system lying above a redshift of 2. We speculate that unlensed relatives of this source may constitute a sizable fraction of the 850 micron source counts.

The SCUBA instrument on the James Clerk Maxwell Telescope has already had an impact on cosmology by detecting relatively large numbers of dusty galaxies at high redshift. Apart from identifying well-detected sources, such data can also be mined for information about fainter sources and their correlations, as revealed through low level fluctuations in SCUBA maps. As a first step in this direction we analyse a small SCUBA data-set as if it were obtained from a Cosmic Microwave Background (CMB) differencing experiment. This enables us to place limits on CMB anisotropy at 850 microns. Expressed as Q_{flat}, the quadrupole expectation value for a flat power spectrum, the limit is 152 microKelvin at 95 per cent confidence, corresponding to C_0^{1/2} < 355 microKelvin for a Gaussian autocorrelation function, with a coherence angle of about 20--25 arcsec; These results could easily be reinterpretted in terms of any other fluctuating sky signal. This is currently the best limit for these scales at high frequency, and comparable to limits at similar angular scales in the radio. Even with such a modest data-set, it is possible to put a constraint on the slope of the SCUBA counts at the faint end, since even randomly distributed sources would lead to fluctuations. Future analysis of sky correlations in more extensive data-sets ought to yield detections, and hence additional information on source counts and clustering.