The InfraRed Imaging Spectrograph (IRIS) is the first-light client instrument for the Narrow Field Infrared Adaptive Optics System (NFIRAOS) on the Thirty Meter Telescope (TMT). IRIS includes three natural guide star (NGS) On-Instrument Wavefront Sensors (OIWFS) to measure tip/tilt and focus errors in the instrument focal plane. NFIRAOS also has an internal natural guide star wavefront sensor, and IRIS and NFIRAOS must precisely coordinate the motions of their wavefront sensor positioners to track the locations of NGSs while the telescope is dithering (offsetting the telescope to cover more area), to avoid a costly re-acquisition time penalty. First, we present an overview of the sequencing strategy for all of the involved subsystems. We then predict the motion of the telescope during dithers based on finite-element models provided by TMT, and finally analyze latency and jitter issues affecting the propagation of position demands from the telescope control system to individual motor controllers.

Infrared Imaging Spectrograph (IRIS) is the first light instrument for the Thirty Meter Telescope (TMT) that consists of a near-infrared (0.84 to 2.4 micron) imager and integral field spectrograph (IFS) which operates at the diffraction-limit utilizing the Narrow-Field Infrared Adaptive Optics System (NFIRAOS). The imager will have a 34 arcsec x 34 arcsec field of view with 4 milliarcsecond (mas) pixels. The IFS consists of a lenslet array and slicer, enabling four plate scales from 4 mas to 50 mas, multiple gratings and filters, which in turn will operate hundreds of individual modes. IRIS, operating in concert with NFIRAOS will pose many challenges for the data reduction system (DRS). Here we present the updated design of the real-time and post-processing DRS. The DRS will support two modes of operation of IRIS: (1) writing the raw readouts sent from the detectors and performing the sampling on all of the readouts for a given exposure to create a raw science frame; and (2) reduction of data from the imager, lenslet array and slicer IFS. IRIS is planning to save the raw readouts for a given exposure to enable sophisticated processing capabilities to the end users, such as the ability to remove individual poor seeing readouts to improve signal-to-noise, or from advanced knowledge of the point spread function (PSF). The readout processor (ROP) is a key part of the IRIS DRS design for writing and sampling of the raw readouts into a raw science frame, which will be passed to the TMT data archive. We discuss the use of sub-arrays on the imager detectors for saturation/persistence mitigation, on-detector guide windows, and fast readout science cases (< 1 second).

The InfraRed Imaging Spectrograph (IRIS) is a first-light instrument for the Thirty Meter Telescope (TMT) that will be used to sample the corrected adaptive optics field by NFIRAOS with a near-infrared (0.8 - 2.4 $\mu$m) imaging camera and Integral Field Spectrograph (IFS). In order to understand the science case specifications of the IRIS instrument, we use the IRIS data simulator to characterize photometric precision and accuracy of the IRIS imager. We present the results of investigation into the effects of potential ghosting in the IRIS optical design. Each source in the IRIS imager field of view results in ghost images on the detector from IRIS's wedge filters, entrance window, and Atmospheric Dispersion Corrector (ADC) prism. We incorporated each of these ghosts into the IRIS simulator by simulating an appropriate magnitude point source at a specified pixel distance, and for the case of the extended ghosts redistributing flux evenly over the area specified by IRIS's optical design. We simulate the ghosting impact on the photometric capabilities, and found that ghosts generally contribute negligible effects on the flux counts for point sources except for extreme cases where ghosts coalign with a star of $\Delta$m$>$2 fainter than the ghost source. Lastly, we explore the photometric precision and accuracy for single sources and crowded field photometry on the IRIS imager.

We have used the Submillimeter Array at 860$\,\mu$m to observe the brightest SCUBA-2 sources in 4$\,$deg$^{2}$ of the Cosmology Legacy Survey. We have targeted 75 of the brightest single-dish SCUBA-2 850$\,\mu$m sources down to $S_{850}\,{\approx}\,8\,$mJy, achieving an average synthesized beam of 2.4$^{\prime\prime}$ and an average rms of $\sigma_{860}\,{=}\,1.5\,$mJy in our primary beam-corrected maps. We searched our maps for $4\sigma$ peaks, corresponding to $S_{860}\,{\gtrsim}\,6\,$mJy sources, and detected 59 single galaxies and three pairs of galaxies. We include in our study 28 archival observations, bringing our sample size to 103 bright single-dish submillimetre sources with interferometric follow-up. We compute the cumulative and differential number counts of our sample, finding them to overlap with previous single-dish survey number counts within the uncertainties, although our cumulative number count is systematically lower than the parent SCUBA-2 cumulative number count by $24\,{\pm}\,6$ per cent between 11 and 15$\,$mJy. We estimate the probability that a ${\gtrsim}\,10\,$mJy single-dish submillimetre source resolves into two or more galaxies with similar flux densities, causing a significant change in the number counts, to be about 15 per cent. Assuming the remaining 85 per cent of the targets are ultra-luminous starburst galaxies between $z\,{=}\,2$-3, we find a likely volume density of ${\gtrsim}\,400\,$M$_{\odot}\,$yr$^{-1}$ sources to be ${\sim}\,3^{+0.7}_{-0.6}\,{\times}\,10^{-7}\,$Mpc$^{-3}$. We show that the descendants of these galaxies could be ${\gtrsim}\,4\,{\times}\,10^{11}\,$M$_{\odot}$ local quiescent galaxies, and that about 10 per cent of their total stellar mass would have formed during these short bursts of star-formation.

We present a multi-wavelength analysis of galaxies selected at 450 and 850um from the deepest SCUBA-2 observations in the Extended Groth Strip (EGS) field, which have an average depth of sigma_450=1.9 and sigma_850=0.46 mJy/beam over ~70 sq. arcmin. The final sample comprises 95 sources: 56 (59%) are detected at both wavelengths, 31 (33%) are detected only at 850um, and 8 (8%) are detected only at 450um. We identify counterparts for 75% of the whole sample. The redshift distributions of the 450 and 850um samples peak at different redshifts with median values of z=1.66 +/- 0.18 and z=2.30 +/- 0.20, respectively. However, the two populations have similar IR luminosities, SFRs, and stellar masses, with mean values of (1.5 +/- 0.2) x 10^12 L_sun, (150 +/- 20) M_sun/yr, and (9.0 +/- 0.6) x 10^10 M_sun, respectively. This places most of our sources (>85%) on the high-mass end of the `main-sequence' of star-forming galaxies. Exploring the IR excess vs UV-slope (IRX-beta) relation we find that the most luminous galaxies are consistent with the Meurer law, while the less luminous galaxies lie below this relation. Using the results of a two-dimensional modelling of the HST H_160-band imaging, we derive a median Sersic index of n=1.4 +0.3 -0.1 and a median half-light radius of R_1/2 = 4.8 +/ 0.4 kpc. Based on a visual-like classification in the same band, we find that the dominant component for most of the galaxies at all redshifts is a disk-like structure, although there is a transition from irregular disks to disks with a spheroidal component at z~1.4, which morphologically supports the scenario of SMGs as progenitors of massive elliptical galaxies.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $\sigma_{450}=1.2$ mJy/beam and $\sigma_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

Submillimetre galaxies (SMGs) are among the most luminous dusty galaxies in the Universe, but their true nature remains unclear; are SMGs the progenitors of the massive elliptical galaxies we see in the local Universe, or are they just a short-lived phase among more typical star-forming galaxies? To explore this problem further, we investigate the clustering of SMGs identified in the SCUBA-2 Cosmology Legacy Survey. We use a catalogue of submillimetre ($850\mu$m) source identifications derived using a combination of radio counterparts and colour/IR selection to analyse a sample of 914 SMGs in the UKIDSS Ultra Deep Survey (UDS), making this the largest high redshift sample of these galaxies to date. Using angular cross-correlation techniques, we estimate the halo masses for this large sample of SMGs and compare them with passive and star-forming galaxies selected in the same field. We find that SMGs, on average, occupy high-mass dark matter halos (M$_{\text{halo}} >10^{13}$M$_{\odot}$) at redshifts $z > 2.5$, consistent with being the progenitors of massive quiescent galaxies in present-day galaxy clusters. We also find evidence of downsizing, in which SMG activity shifts to lower mass halos at lower redshifts. In terms of their clustering and halo masses, SMGs appear to be consistent with other star-forming galaxies at a given redshift.

IRIS (InfraRed Imaging Spectrograph) is the diffraction-limited first light instrument for the Thirty Meter Telescope (TMT) that consists of a near-infrared (0.84 to 2.4 $\mu$m) imager and integral field spectrograph (IFS). The IFS makes use of a lenslet array and slicer for spatial sampling, which will be able to operate in 100's of different modes, including a combination of four plate scales from 4 milliarcseconds (mas) to 50 mas with a large range of filters and gratings. The imager will have a field of view of 34$\times$34 arcsec$^{2}$ with a plate scale of 4 mas with many selectable filters. We present the preliminary design of the data reduction system (DRS) for IRIS that need to address all of these observing modes. Reduction of IRIS data will have unique challenges since it will provide real-time reduction and analysis of the imaging and spectroscopic data during observational sequences, as well as advanced post-processing algorithms. The DRS will support three basic modes of operation of IRIS; reducing data from the imager, the lenslet IFS, and slicer IFS. The DRS will be written in Python, making use of open-source astronomical packages available. In addition to real-time data reduction, the DRS will utilize real-time visualization tools, providing astronomers with up-to-date evaluation of the target acquisition and data quality. The quicklook suite will include visualization tools for 1D, 2D, and 3D raw and reduced images. We discuss the overall requirements of the DRS and visualization tools, as well as necessary calibration data to achieve optimal data quality in order to exploit science cases across all cosmic distance scales.

The InfraRed Imaging Spectrograph (IRIS) will be a first-light client instrument for the Narrow Field Infrared Adaptive Optics System (NFIRAOS) on the Thirty Meter Telescope. IRIS includes three configurable tip/tilt (TT) or tip/tilt/focus (TTF) On-Instrument Wavefront Sensors (OIWFS). These sensors are positioned over natural guide star (NGS) asterisms using movable polar-coordinate pick-off arms (POA) that patrol an approximately 2-arcminute circular field-of-view (FOV). The POAs are capable of colliding with one another, so an algorithm for coordinated motion that avoids contact is required. We have adopted an approach in which arm motion is evaluated using the gradient descent of a scalar potential field that includes an attractive component towards the goal configuration (locations of target stars), and repulsive components to avoid obstacles (proximity to adjacent arms). The resulting vector field is further modified by adding a component transverse to the repulsive gradient to avoid problematic local minima in the potential. We present path planning simulations using this computationally inexpensive technique, which exhibit smooth and efficient trajectories.

The InfraRed Imaging Spectrograph (IRIS) will be the first light adaptive optics instrument on the Thirty Meter Telescope (TMT). IRIS is being built by a collaboration between Caltech, the University of California, NAOJ and NRC Herzberg. We present novel aspects of the Support Structure, Rotator and On-Instrument Wavefront Sensor systems being developed at NRC Herzberg. IRIS is suspended from the bottom port of the Narrow Field Infrared Adaptive Optics System (NFIRAOS), and provides its own image de-rotation to compensate for sidereal rotation of the focal plane. This arrangement is a challenge because NFIRAOS is designed to host two other science instruments, which imposes strict mass requirements on IRIS. We have been tasked with keeping the instrument mass under seven tonnes which has resulted in a mass reduction of 30 percent for the support structure and rotator compared to the most recent IRIS designs. To accomplish this goal, while still being able to withstand earthquakes, we developed a new design with composite materials. As IRIS is a client instrument of NFIRAOS, it benefits from NFIRAOS's superior AO correction. IRIS assists this correction by sensing low-order aberrations with an On-Instrument Wavefront Sensor (OIWFS). The OIWFS consists of three independently positioned natural guide star wavefront sensors that patrol a 2-arcminute field of view. We expect tip-tilt measurements from faint stars within the IRIS imager focal plane will further stabilize the delivered image quality. We describe how the use of On-Detector Guide Windows (ODGWs) in the IRIS imager can be incorporated into the AO correction. Finally, we present our strategies for acquiring and tracking sources with this complex AO system, and for mitigating and measuring the various potential sources of image blur and misalignment due to properties of the mechanical structure and interfaces. (Abridged)

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

Thermal emission from cold dust grains in giant molecular clouds can be used to probe the physical properties, such as density, temperature and emissivity in star-forming regions. We present the SCUBA-2 shared-risk observations at 450 $\mu$m and 850 $\mu$m of the Orion A molecular cloud complex taken at the James Clerk Maxwell Telescope (JCMT). Previous studies showed that molecular emission lines can contribute significantly to the measured fluxes in those continuum bands. We use the HARP $^{12}$CO J=3-2 integrated intensity map for Orion A in order to evaluate the molecular line contamination and its effects on the SCUBA-2 maps. With the corrected fluxes, we have obtained a new spectral index $\alpha$ map for the thermal emission of dust in the well-known integral-shaped filament. Furthermore, we compare a sample of 33 sources, selected over the Orion A molecular cloud complex for their high $^{12}$CO J=3-2 line contamination, to 27 previously identified clumps in OMC-4. This allows us to quantify the effect of line contamination on the ratio of 850 $\mu$m to 450 $\mu$m flux densities and how it modifies the deduced spectral index of emissivity $\beta$ for the dust grains. We also show that at least one Spitzer-identified protostellar core in OMC-5 has a $^{12}$CO J=3-2 contamination level of 16 %. Furthermore, we find the strongest contamination level (44 %) towards a young star with disk near OMC-2. This work is part of the JCMT Gould Belt Legacy Survey.

We describe a search for submillimeter emission in the vicinity of one of the most distant, luminous galaxies known, HerMES FLS3 at z=6.34, exploiting it as a signpost to a potentially biased region of the early Universe, as might be expected in hierarchical structure formation models. Imaging to the confusion limit with the innovative, wide-field submillimeter bolometer camera, SCUBA-2, we are sensitive to colder and/or less luminous galaxies in the surroundings of HFLS3. We use the Millennium Simulation to illustrate that HFLS3 may be expected to have companions if it is as massive as claimed, but find no significant evidence from the surface density of SCUBA-2 galaxies in its vicinity, or their colors, that HFLS3 marks an over-density of dusty, star-forming galaxies. We cannot rule out the presence of dusty neighbours with confidence, but deeper 450-um imaging has the potential to more tightly constrain the redshifts of nearby galaxies, at least one of which likely lies at z>~5. If associations with HFLS3 can be ruled out, this could be taken as evidence that HFLS3 is less biased than a simple extrapolation of the Millennium Simulation may imply. This could suggest either that it represents a rare short-lived, but highly luminous, phase in the evolution of an otherwise typical galaxy, or that this system has suffered amplification due to a foreground gravitational lens and so is not as intrinsically luminous as claimed.

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on subparsec to 100 parsec scales, leading to the formation of starless cores, is not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 $\mu$m maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to a large-scale filament shape model for Lupus I, or relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 $\mu$m with the Balloon-borne Large Aperture Telescope for Polarimetry (BLASTPol) during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics - including secondary filaments that often run orthogonally to the primary filament - and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core.

We mapped five massive star forming regions with the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT). Temperature and column density maps are obtained from the SCUBA-2 450 and 850 $\mu$m images. Most of the dense clumps we find have central temperatures below 20 K with some as cold as 8 K, suggesting that they have no internal heating due to the presence of embedded protostars. This is surprising, because at the high densities inferred from these images and at these low temperatures such clumps should be unstable, collapsing to form stars and generating internal heating. The column densities at the clump centres exceed 10$^{23}$ cm$^{-2}$, and the derived peak visual extinction values are from 25-500 mag for $\beta$ = 1.5-2.5, indicating highly opaque centres. The observed cloud gas masses range from $\sim$ 10 to 10$^{3}$ M$_{\odot}$. The outer regions of the clumps follow an $r^{-2.36\pm0.35}$ density distribution and this power-law structure is observed outside of typically 10$^{4}$ AU. All these findings suggest that these clumps are high-mass starless clumps and most likely contain high-mass starless cores.

We analyse new SCUBA-2 submillimeter and archival SPIRE far-infrared imaging of a z=1.62 cluster, Cl0218.3-0510, which lies in the UKIDSS/UDS field of the SCUBA-2 Cosmology Legacy Survey. Combining these tracers of obscured star formation activity with the extensive photometric and spectroscopic information available for this field, we identify 31 far-infrared/submillimeter-detected probable cluster members with bolometric luminosities >1e12 Lo and show that by virtue of their dust content and activity, these represent some of the reddest and brightest galaxies in this structure. We exploit Cycle-1 ALMA submillimeter continuum imaging which covers one of these sources to confirm the identification of a SCUBA-2-detected ultraluminous star-forming galaxy in this structure. Integrating the total star-formation activity in the central region of the structure, we estimate that it is an order of magnitude higher (in a mass-normalised sense) than clusters at z~0.5-1. However, we also find that the most active cluster members do not reside in the densest regions of the structure, which instead host a population of passive and massive, red galaxies. We suggest that while the passive and active populations have comparable near-infrared luminosities at z=1.6, M(H)~-23, the subsequent stronger fading of the more active galaxies means that they will evolve into passive systems at the present-day which are less luminous than the descendants of those galaxies which were already passive at z~1.6 (M(H)~-20.5 and M(H)~-21.5 respectively at z~0). We conclude that the massive galaxy population in the dense cores of present-day clusters were already in place at z=1.6 and that in Cl0218.3-0510 we are seeing continuing infall of less extreme, but still ultraluminous, star-forming galaxies onto a pre-existing structure.

MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.

We report the serendipitous detection of a significant overdensity of Herschel-SPIRE 250 micron sources in the vicinity of MRC1138-26. We use an adaptive kernel density estimate to quantify the significance, including a comparison with other fields. The overdensity has a size of ~3.5-4' and stands out at ~5sigma with respect to the background estimate. No features with similar significance were found in four extragalactic control fields: GOODS-North, Lockman, COSMOS and UDS. The chance of having a similar overdensity in a field with the same number but randomly distributed sources is less than 2%. The clump is also visible as a low surface brightness feature in the Planck 857 GHz map. We detect 76 sources at 250 micron (with a signal-to-noise ratio greater than 3), in a region of 4' radius; 43 of those are above a flux density limit of 20 mJy. This is a factor of 3.6 in excess over the average in the four control fields, considering only the sources above 20 mJy. We also find an excess in the number counts of sources with 250 micron flux densities between 30 and 40 mJy, compared to deep extragalactic blank-field number counts. Assuming a fixed dust temperature (30 K) and emissivity (beta=1.5) a crude, blackbody-derived redshift distribution, zBB, of the detected sources is significantly different from the distributions in the control fields and exhibits a significant peak at zBB ~ 1.5, although the actual peak redshift is highly degenerate with the temperature. We tentatively suggest, based on zBB and the similar S250/S350 colours of the sources within the peak, that a significant fraction of the sources in the clump may be at a similar redshift. Since the overdensity lies ~7' south of the z=2.16 Spiderweb protocluster MRC1138-26, an intriguing possibility (that is presently unverifiable given the data in hand) is that it lies within the same large-scale structure.(abridged)

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 Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 {\mu}m. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is an instrument operating on the 15-m James Clerk Maxwell Telescope, nominally consisting of 5120 bolometers in each of two simultaneous imaging bands centred over 450 and 850 um. The camera is operated by scanning across the sky and recording data at a rate of 200 Hz. As the largest of a new generation of multiplexed kilopixel bolometer cameras operating in the (sub)millimetre, SCUBA-2 data analysis represents a significant challenge. We describe the production of maps using the Sub-Millimetre User Reduction Facility (SMURF) in which we have adopted a fast, iterative approach to map-making that enables data reduction on single, modern, high-end desktop computers, with execution times that are typically shorter than the observing times. SMURF is used in an automated setting, both at the telescope for real-time feedback to observers, as well as for the production of science products for the JCMT Science Archive at the Canadian Astronomy Data Centre. Three detailed case studies are used to: (i) explore convergence properties of the map-maker using simple prior constraints (Uranus -- a point source); (ii) achieve the white-noise limit for faint point-source studies (extragalactic blank-field survey of the Lockman Hole); and (iii) demonstrate that our strategy is capable of recovering angular scales comparable to the size of the array footprint (approximately 5 arcmin) for bright extended sources (star-forming region M17).

SCUBA-2 is a 10000-bolometer submillimetre camera on the James Clerk Maxwell Telescope (JCMT). The instrument commissioning was completed in September 2011, and full science operations began in October 2011. To harness the full potential of this powerful new astronomical tool, the instrument calibration must be accurate and well understood. To this end, the algorithms for calculating the line-of-sight opacity have been improved, and the derived atmospheric extinction relationships at both wavebands of the SCUBA-2 instrument are presented. The results from over 500 primary and secondary calibrator observations have allowed accurate determination of the flux conversion factors (FCF) for the 850 and 450 micron arrays. Descriptions of the instrument beam-shape and photometry methods are presented. The calibration factors are well determined, with relative calibration accuracy better than 5 per cent at 850 microns and 10 per cent at 450 microns, reflecting the success of the derived opacity relations as well as the stability of the performance of the instrument over several months. The sample-size of the calibration observations and accurate FCFs have allowed the determination of the 850 and 450 micron fluxes of several well-known submillimetre sources, and these results are compared with previous measurements from SCUBA.

SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 microns, the vast increase in pixel count means that SCUBA-2 maps the sky 100-150 times faster than the previous SCUBA instrument. In this paper we present an overview of the instrument, discuss the physical characteristics of the superconducting detector arrays, outline the observing modes and data acquisition, and present the early performance figures on the telescope. We also showcase the capabilities of the instrument via some early examples of the science SCUBA-2 has already undertaken. In February 2012, SCUBA-2 began a series of unique legacy surveys for the JCMT community. These surveys will take 2.5 years and the results are already providing complementary data to the shorter wavelength, shallower, larger-area surveys from Herschel. The SCUBA-2 surveys will also provide a wealth of information for further study with new facilities such as ALMA, and future telescopes such as CCAT and SPICA.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

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 Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

SCUBA-2 is the largest submillimetre wide-field bolometric camera ever built. This 43 square arc-minute field-of-view instrument operates at two wavelengths (850 and 450 microns) and has been installed on the James Clerk Maxwell Telescope on Mauna Kea, Hawaii. SCUBA-2 has been successfully commissioned and operational for general science since October 2011. This paper presents an overview of the on-sky performance of the instrument during and since commissioning in mid-2011. The on-sky noise characteristics and NEPs of the 450 and 850 micron arrays, with average yields of approximately 3400 bolometers at each wavelength, will be shown. The observing modes of the instrument and the on-sky calibration techniques are described. The culmination of these efforts has resulted in a scientifically powerful mapping camera with sensitivities that allow a square degree of sky to be mapped to 10 mJy/beam rms at 850 micron in 2 hours and 60 mJy/beam rms at 450 micron in 5 hours in the best weather.

The AzTEC 1.1 mm survey of the SCUBA HAlf Degree Extragalactic Survey (SHADES) fields is the largest (0.7 deg2) blank-field millimetre-wavelength survey undertaken to date at a resolution of ~18" and a depth of ~1 mJy. We have used the deep optical-to-radio multi-wavelength data in the SHADES Lockman Hole East and SXDF/UDS fields to obtain galaxy identifications for ~64% (~80% with tentative identifications) of the 148 AzTEC-SHADES 1.1 mm sources, exploiting deep radio and 24 um data complemented by methods based on 8 um flux-density and red optical-infrared (i-K) colour. This unusually high identification rate can be attributed to the relatively bright millimetre-wavelength flux-density threshold, combined with the relatively deep supporting multi-frequency data. We have further exploited the optical-mid-infrared-radio data to derive a ~60% (~75% with tentative identifications) complete redshift distribution for the AzTEC-SHADES sources, yielding a median redshift of z~2.2, with a high-redshift tail extending to at least z~4. Despite the larger area probed by the AzTEC survey relative to the original SCUBA SHADES imaging, the redshift distribution of the AzTEC sources is consistent with that displayed by the SCUBA sources, and reinforces tentative evidence that the redshift distribution of mm/sub-mm sources in the Lockman Hole field is significantly different from that found in the SXDF/UDS field. Comparison with simulated surveys of similar scale extracted from semi-analytic models based on the Millennium simulation indicates that this is as expected if the mm/sub-mm sources are massive (M>10^11 Mo) star-forming galaxies tracing large-scale structures over scales of 10-20 Mpc. This confirms the importance of surveys covering several square degrees (as now underway with SCUBA2) to obtain representative samples of bright (sub)mm-selected galaxies.

We investigate the properties (e.g. star formation rate, dust attentuation, stellar mass and metallicity) of a sample of infrared luminous galaxies at z \sim 1 via near-IR spectroscopy with Subaru-FMOS. Our sample consists of Herschel SPIRE and Spitzer MIPS selected sources in the COSMOS field with photometric redshifts in the range 0.7 < z-phot < 1.8, which have been targeted in 2 pointings (0.5 sq. deg.) with FMOS. We find a modest success rate for emission line detections, with candidate H{\alpha} emission lines detected for 57 of 168 SPIRE sources (34 per cent). By stacking the near-IR spectra we directly measure the mean Balmer decrement for the H{\alpha} and H{\beta} lines, finding a value of <E(B-V)> = 0.51\pm0.27 for <LIR> = 10^12 Lsol sources at <z> = 1.36. By comparing star formation rates estimated from the IR and from the dust uncorrected H{\alpha} line we find a strong relationship between dust attenuation and star formation rate. This relation is broadly consistent with that previously seen in star-forming galaxies at z ~ 0.1. Finally, we investigate the metallicity via the N2 ratio, finding that z ~ 1 IR-selected sources are indistinguishable from the local mass-metallicity relation. We also find a strong correlation between dust attentuation and metallicity, with the most metal-rich IR-sources experiencing the largest levels of dust attenuation.

The submillimeter opacity of dust in the diffuse Galactic interstellar medium (ISM) has been quantified using a pixel-by-pixel correlation of images of continuum emission with a proxy for column density. We used three BLAST bands at 250, 350, and 500 \mu m and one IRAS at 100 \mu m. The proxy is the near-infrared color excess, E(J-Ks), obtained from 2MASS. Based on observations of stars, we show how well this color excess is correlated with the total hydrogen column density for regions of moderate extinction. The ratio of emission to column density, the emissivity, is then known from the correlations, as a function of frequency. The spectral distribution of this emissivity can be fit by a modified blackbody, whence the characteristic dust temperature T and the desired opacity \sigma_e(1200) at 1200 GHz can be obtained. We have analyzed 14 regions near the Galactic plane toward the Vela molecular cloud, mostly selected to avoid regions of high column density (N_H > 10^{22} cm^-2) and small enough to ensure a uniform T. We find \sigma_e(1200) is typically 2 to 4 x 10^{-25} cm^2/H and thus about 2 to 4 times larger than the average value in the local high Galactic latitude diffuse atomic ISM. This is strong evidence for grain evolution. There is a range in total power per H nucleon absorbed (re-radiated) by the dust, reflecting changes in the interstellar radiation field and/or the dust absorption opacity. These changes affect the equilibrium T, which is typically 15 K, colder than at high latitudes. Our analysis extends, to higher opacity and lower T, the trend of increasing opacity with decreasing T that was found at high latitudes. The recognition of changes in the emission opacity raises a cautionary flag because all column densities deduced from dust emission maps, and the masses of compact structures within them, depend inversely on the value adopted.

The source counts of galaxies discovered at sub-millimetre and millimetre wavelengths provide important information on the evolution of infrared-bright galaxies. We combine the data from six blank-field surveys carried out at 1.1 mm with AzTEC, totalling 1.6 square degrees in area with root-mean-square depths ranging from 0.4 to 1.7 mJy, and derive the strongest constraints to date on the 1.1 mm source counts at flux densities S(1100) = 1-12 mJy. Using additional data from the AzTEC Cluster Environment Survey to extend the counts to S(1100) ~ 20 mJy, we see tentative evidence for an enhancement relative to the exponential drop in the counts at S(1100) ~ 13 mJy and a smooth connection to the bright source counts at >20 mJy measured by the South Pole Telescope; this excess may be due to strong lensing effects. We compare these counts to predictions from several semi-analytical and phenomenological models and find that for most the agreement is quite good at flux densities > 4 mJy; however, we find significant discrepancies (>3sigma) between the models and the observed 1.1 mm counts at lower flux densities, and none of them are consistent with the observed turnover in the Euclidean-normalised counts at S(1100) < 2 mJy. Our new results therefore may require modifications to existing evolutionary models for low luminosity galaxies. Alternatively, the discrepancy between the measured counts at the faint end and predictions from phenomenological models could arise from limited knowledge of the spectral energy distributions of faint galaxies in the local Universe.

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 {\mu}m (Spitzer/MIPS); 70, 100, and 160{\mu}m (Herschel/PACS); 250, 350, 500{\mu}m (BLAST); and 870{\mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

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 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 constrain the evolution of the rest-frame far-infrared (FIR) luminosity function out to high redshift, by combining several pieces of complementary information provided by the deep Balloon-borne Large-Aperture Submillimeter Telescope surveys at 250, 350 and 500 micron, as well as other FIR and millimetre data. Unlike most other phenomenological models, we characterise the uncertainties in our fitted parameters using Monte Carlo Markov Chains. We use a bivariate local luminosity function that depends only on FIR luminosity and 60-to-100 micron colour, along with a single library of galaxy spectral energy distributions indexed by colour, and apply simple luminosity and density evolution. We use the surface density of sources, Cosmic Infrared Background (CIB) measurements and redshift distributions of bright sources, for which identifications have been made, to constrain this model. The precise evolution of the FIR luminosity function across this crucial range has eluded studies at longer wavelengths (e.g., using SCUBA and MAMBO) and at shorter wavelengths (e.g., Spitzer), and should provide a key piece of information required for the study of galaxy evolution. Our adoption of Monte Carlo methods enables us not only to find the best-fit evolution model, but also to explore correlations between the fitted parameters. Our model-fitting approach allows us to focus on sources of tension coming from the combination of data-sets. We specifically find that our choice of parameterisation has difficulty fitting the combination of CIB measurements and redshift distribution of sources near 1 mm. Existing and future data sets will be able to dramatically improve the fits, as well as break strong degeneracies among the models. [abridged]

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 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 identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

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.

We have carried out the first general submillimeter analysis of the field towards GRSMC 45.46+0.05, a massive star forming region in Aquila. The deconvolved 6 deg^2 (3\degree X 2\degree) maps provided by BLAST in 2005 at 250, 350, and 500 micron were used to perform a preliminary characterization of the clump population previously investigated in the infrared, radio, and molecular maps. Interferometric CORNISH data at 4.8 GHz have also been used to characterize the Ultracompact HII regions (UCHIIRs) within the main clumps. By means of the BLAST maps we have produced an initial census of the submillimeter structures that will be observed by Herschel, several of which are known Infrared Dark Clouds (IRDCs). Our spectral energy distributions of the main clumps in the field, located at ~7 kpc, reveal an active population with temperatures of T~35-40 K and masses of ~10^3 Msun for a dust emissivity index beta=1.5. The clump evolutionary stages range from evolved sources, with extended HII regions and prominent IR stellar population, to massive young stellar objects, prior to the formation of an UCHIIR.The CORNISH data have revealed the details of the stellar content and structure of the UCHIIRs. In most cases, the ionizing stars corresponding to the brightest radio detections are capable of accounting for the clump bolometric luminosity, in most cases powered by embedded OB stellar clusters.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

SCUBA-2 is the largest submillimetre array camera in the world and was commissioned on the James Clerk Maxwell Telescope (JCMT) with two arrays towards the end of 2009. A period of shared-risks observing was then completed and the full planned complement of 8 arrays, 4 at 850 microns and 4 at 450 microns, are now installed and ready to be commissioned. SCUBA-2 has 10,240 bolometers, corresponding to a data rate of 8 MB/s when sampled at the nominal rate of 200 Hz. The pipeline produces useful maps in near real time at the telescope and often publication quality maps in the JCMT Science Archive (JSA) hosted at the Canadian Astronomy Data Centre (CADC).

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey at 250, 350, and 500 micron, combined with rich datasets for this well-studied region. Our primary goal is to investigate the early stages of high mass star formation. We have detected 184 compact sources in various stages of evolution across all three BLAST bands. From their well-constrained spectral energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced from their position in the luminosity-mass (L-M) diagram. The BLAST spectral coverage, near the peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be seen by earlier shorter-wavelength surveys like IRAS. We detect thermal emission from infrared dark clouds and investigate the phenomenon of cold ``starless cores" more generally. Spitzer images of these cold sources often show stellar nurseries, but these potential sites for massive star formation are ``starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: some cold massive compact sources might never form a cluster containing massive stars; and clusters with massive stars might not have an identifiable compact cold massive precursor.

We carry out a multi-wavelength study of individual galaxies detected by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and identified at other wavelengths, using data spanning the radio to the ultraviolet (UV). We develop a Monte Carlo method to account for flux boosting, source blending, and correlations among bands, which we use to derive deboosted far-infrared (FIR) luminosities for our sample. We estimate total star-formation rates for BLAST counterparts with z < 0.9 by combining their FIR and UV luminosities. Star formation is heavily obscured at L_FIR > 10^11 L_sun, z > 0.5, but the contribution from unobscured starlight cannot be neglected at L_FIR < 10^11 L_sun, z < 0.25. We assess that about 20% of the galaxies in our sample show indication of a type-1 active galactic nucleus (AGN), but their submillimeter emission is mainly due to star formation in the host galaxy. We compute stellar masses for a subset of 92 BLAST counterparts; these are relatively massive objects, with a median mass of ~10^11 M_sun, which seem to link the 24um and SCUBA populations, in terms of both stellar mass and star-formation activity. The bulk of the BLAST counterparts at z<1 appear to be run-of-the-mill star-forming galaxies, typically spiral in shape, with intermediate stellar masses and practically constant specific star-formation rates. On the other hand, the high-z tail of the BLAST counterparts significantly overlaps with the SCUBA population, in terms of both star-formation rates and stellar masses, with observed trends of specific star-formation rate that support strong evolution and downsizing.

We present a joint analysis of the overlapping BLAST 250, 350, 500um, and LABOCA 870um observations (from the LESS survey) of the ECDF-S. Out to z~3, the BLAST filters sample near the peak wavelength of far-infrared (FIR) emission from galaxies (rest-frame wavelengths ~60--200um), primarily produced by dust heated through absorption in star-forming clouds. However, identifying counterparts to individual BLAST peaks is very challenging, given the large beams (FWHM 36--60"). In contrast, the 870um observations have a significantly smaller 19" FWHM beam, and are sensitive to higher redshifts (z~1--5, and potentially beyond) due to the more favourable negative K-correction. We use the LESS data, as well as deep Spitzer and VLA imaging, to identify 118 individual sources that produce significant emission in the BLAST bands. We characterize the temperatures and FIR luminosities for a subset of 69 sources which have well-measured submm SEDs and redshift measurements out to z~3. For flux-limited sub-samples in each BLAST band, and a dust emissivity index \beta=2.0, we find a median temperature T=30K (all bands) as well as median redshifts: z=1.1 (interquartile range 0.2--1.9) for S250 > 40mJy; z=1.3 (interquartile range 0.6--2.1) for S350 > 30mJy; and z=1.6 (interquartile range 1.3--2.3) for S500 > 20mJy. Taking into account selection effects for our survey (a bias toward detecting lower-temperature galaxies), we find no evidence for evolution in the local FIR-temperature correlation out to z~2.5. Comparing with star-forming galaxy SED templates, about 8% of our sample appears to exhibit significant excesses in the radio and/or mid-IR, consistent with those sources harbouring an AGN. We describe the following techniques in two appendices: our `matched filter' for identifying sources in the presence of point-source confusion; and our approach for identifying counterparts using likelihood ratios.

We present an implementation of the iterative flux-conserving Lucy-Richardson (L-R) deconvolution method of image restoration for maps produced by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). We have analyzed its performance and convergence extensively through simulations and cross-correlations of the deconvolved images with available highresolution maps. We present new science results from two BLAST surveys, in the Galactic regions K3-50 and IC 5146, further demonstrating the benefits of performing this deconvolution. We have resolved three clumps within a radius of 4.'5 inside the star-forming molecular cloud containing K3-50. Combining the well-resolved dust emission map with available multi-wavelength data, we have constrained the Spectral Energy Distributions (SEDs) of five clumps to obtain masses (M), bolometric luminosities (L), and dust temperatures (T). The L-M diagram has been used as a diagnostic tool to estimate the evolutionary stages of the clumps. There are close relationships between dust continuum emission and both 21-cm radio continuum and 12CO molecular line emission. The restored extended large scale structures in the Northern Streamer of IC 5146 have a strong spatial correlation with both SCUBA and high resolution extinction images. A dust temperature of 12 K has been obtained for the central filament. We report physical properties of ten compact sources, including six associated protostars, by fitting SEDs to multi-wavelength data. All of these compact sources are still quite cold (typical temperature below ~ 16 K) and are above the critical Bonner-Ebert mass. They have associated low-power Young Stellar Objects (YSOs). Further evidence for starless clumps has also been found in the IC 5146 region.

We present results from a survey carried out by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) on a 9 deg^2 field near the South Ecliptic Pole at 250, 350 and 500 {\mu}m. The median 1{\sigma} depths of the maps are 36.0, 26.4 and 18.4 mJy, respectively. We apply a statistical method to estimate submillimeter galaxy number counts and find that they are in agreement with other measurements made with the same instrument and with the more recent results from Herschel/SPIRE. Thanks to the large field observed, the new measurements give additional constraints on the bright end of the counts. We identify 132, 89 and 61 sources with S/N>4 at 250, 350, 500 {\mu}m, respectively and provide a multi-wavelength combined catalog of 232 sources with a significance >4{\sigma} in at least one BLAST band. The new BLAST maps and catalogs are available publicly at http://blastexperiment.info.

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.

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.

We present a study of the cosmic infrared background, which is a measure of the dust obscured activity in all galaxies in the Universe. We venture to isolate the galaxies responsible for the background at 1mm; with spectroscopic and photometric redshifts we constrain the redshift distribution of these galaxies. We create a deep 1.16mm map (sigma ~ 0.5mJy) by combining the AzTEC 1.1mm and MAMBO 1.2mm datasets in GOODS-N. This combined map contains 41 secure detections, 13 of which are new. By averaging the 1.16mm flux densities of individually undetected galaxies with 24um flux densities > 25uJy, we resolve 31--45 per cent of the 1.16mm background. Repeating our analysis on the SCUBA 850um map, we resolve a higher percentage (40--64 per cent) of the 850um background. A majority of the background resolved (attributed to individual galaxies) at both wavelengths comes from galaxies at z > 1.3. If the ratio of the resolved submillimeter to millimeter background is applied to a reasonable scenario for the origins of the unresolved submillimeter background, 60--88 per cent of the total 1.16mm background comes from galaxies at z > 1.3.

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.

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.

Submillimetre continuum radiation allows us to probe cold objects, particularly the earliest, dusty phases of star formation, high-redshift galaxies and circumstellar disks. The submillimetre window gives a unique view of the physical and dynamical conditions in the neutral and molecular interstellar medium. In the next decade a combination of wide-field surveys with single-dish telescopes and targeted follow-up with ALMA and other facilities should enable rapid progress in answering questions about the origins of planetary systems, stars and galaxies.

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 Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLAST-Pol) is a suborbital mapping experiment designed to study the role played by magnetic fields in the star formation process. BLAST-Pol is the reconstructed BLAST telescope, with the addition of linear polarization capability. Using a 1.8 m Cassegrain telescope, BLAST-Pol images the sky onto a focal plane that consists of 280 bolometric detectors in three arrays, observing simultaneously at 250, 350, and 500 um. The diffraction-limited optical system provides a resolution of 30'' at 250 um. The polarimeter consists of photolithographic polarizing grids mounted in front of each bolometer/detector array. A rotating 4 K achromatic half-wave plate provides additional polarization modulation. With its unprecedented mapping speed and resolution, BLAST-Pol will produce three-color polarization maps for a large number of molecular clouds. The instrument provides a much needed bridge in spatial coverage between larger-scale, coarse resolution surveys and narrow field of view, and high resolution observations of substructure within molecular cloud cores. The first science flight will be from McMurdo Station, Antarctica in December 2010.

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.

We use new large area far infrared maps ranging from 65 - 500 microns obtained with the AKARI and the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) missions to characterize the dust emission toward the Cassiopeia A supernova remnant (SNR). Using the AKARI high resolution data we find a new "tepid" dust grain population at a temperature of ~35K and with an estimated mass of 0.06 solar masses. This component is confined to the central area of the SNR and may represent newly-formed dust in the unshocked supernova ejecta. While the mass of tepid dust that we measure is insufficient by itself to account for the dust observed at high redshift, it does constitute an additional dust population to contribute to those previously reported. We fit our maps at 65, 90, 140, 250, 350, and 500 microns to obtain maps of the column density and temperature of "cold" dust (near 16 K) distributed throughout the region. The large column density of cold dust associated with clouds seen in molecular emission extends continuously from the surrounding interstellar medium to project on the SNR, where the foreground component of the clouds is also detectable through optical, X-ray, and molecular extinction. At the resolution available here, there is no morphological signature to isolate any cold dust associated only with the SNR from this confusing interstellar emission. Our fit also recovers the previously detected "hot" dust in the remnant, with characteristic temperature 100 K.

We discuss high-angular resolution observations of ammonia toward four candidate high-mass starless cores (HMSCs). The cores were identified by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) during its 2005 survey of the Vulpecula region where 60 compact sources were detected simultaneously at 250, 350, and 500 micron. Four of these cores, with no IRAS-PSC or MSX counterparts, were observed with the NRAO Very Large Array (VLA) in the NH3(1,1) and (2,2) spectral lines. Our observations indicate that the four cores are cold (Tk <~ 14K) and show a filamentary and/or clumpy structure. They also show a significant velocity substructure within ~1km/s. The four BLAST cores appear to be colder and more quiescent than other previously observed HMSC candidates, suggesting an earlier stage of evolution.

We present observations at 250, 350, and 500 um of the nearby galaxy cluster Abell 3112 (z=0.075) carried out with BLAST, the Balloon-borne Large Aperture Submillimeter Telescope. Five cluster members are individually detected as bright submillimetre sources. Their far-infrared SEDs and optical colours identify them as normal star-forming galaxies of high mass, with globally evolved stellar populations. They all have B-R colours of 1.38+/-0.08, transitional between the blue, active population and the red, evolved galaxies that dominate the cluster core. We stack to determine the mean submillimetre emission from all cluster members, which is determined to be 16.6+/-2.5, 6.1+/-1.9, and 1.5+/-1.3 mJy at 250, 350, and 500 um, respectively. Stacking analyses of the submillimetre emission of cluster members reveal trends in the mean far-infrared luminosity with respect to cluster-centric radius and Ks-band magnitude. We find that a large fraction of submillimetre emission comes from the boundary of the inner, virialized region of the cluster, at cluster-centric distances around R_500. Stacking also shows that the bulk of the submillimetre emission arises in intermediate-mass galaxies (L<L*), with Ks magnitude ~1 mag fainter than the giant ellipticals. The results and constraints obtained in this work will provide a useful reference for the forthcoming surveys to be conducted on galaxy clusters by Herschel.

The Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) carried out a 250, 350 and 500 micron survey of the galactic plane encompassing the Vela Molecular Ridge, with the primary goal of identifying the coldest dense cores possibly associated with the earliest stages of star formation. Here we present the results from observations of the Vela-D region, covering about 4 square degrees, in which we find 141 BLAST cores. We exploit existing data taken with the Spitzer MIPS, IRAC and SEST-SIMBA instruments to constrain their (single-temperature) spectral energy distributions, assuming a dust emissivity index beta = 2.0. This combination of data allows us to determine the temperature, luminosity and mass of each BLAST core, and also enables us to separate starless from proto-stellar sources. We also analyze the effects that the uncertainties on the derived physical parameters of the individual sources have on the overall physical properties of starless and proto-stellar cores, and we find that there appear to be a smooth transition from the pre- to the proto-stellar phase. In particular, for proto-stellar cores we find a correlation between the MIPS24 flux, associated with the central protostar, and the temperature of the dust envelope. We also find that the core mass function of the Vela-D cores has a slope consistent with other similar (sub)millimeter surveys.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has recently surveyed ~8.7 deg^2 centered on GOODS-South at 250, 350, and 500 microns. In Dye et al. (2009) we presented the catalogue of sources detected at 5-sigma in at least one band in this field and the probable counterparts to these sources in other wavebands. In this paper, we present the results of a redshift survey in which we succeeded in measuring redshifts for 82 of these counterparts. The spectra show that the BLAST counterparts are mostly star-forming galaxies but not extreme ones when compared to those found in the Sloan Digital Sky Survey. Roughly one quarter of the BLAST counterparts contain an active nucleus. We have used the spectroscopic redshifts to carry out a test of the ability of photometric redshift methods to estimate the redshifts of dusty galaxies, showing that the standard methods work well even when a galaxy contains a large amount of dust. We have also investigated the cases where there are two possible counterparts to the BLAST source, finding that in at least half of these there is evidence that the two galaxies are physically associated, either because they are interacting or because they are in the same large-scale structure. Finally, we have made the first direct measurements of the luminosity function in the three BLAST bands. We find strong evolution out to z=1, in the sense that there is a large increase in the space-density of the most luminous galaxies. We have also investigated the evolution of the dust-mass function, finding similar strong evolution in the space-density of the galaxies with the largest dust masses, showing that the luminosity evolution seen in many wavebands is associated with an increase in the reservoir of interstellar matter in galaxies.

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

We describe the application of a statistical method to estimate submillimeter galaxy number counts from confusion limited observations by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). Our method is based on a maximum likelihood fit to the pixel histogram, sometimes called 'P(D)', an approach which has been used before to probe faint counts, the difference being that here we advocate its use even for sources with relatively high signal-to-noise ratios. This method has an advantage over standard techniques of source extraction in providing an unbiased estimate of the counts from the bright end down to flux densities well below the confusion limit. We specifically analyse BLAST observations of a roughly 10 sq. deg. map centered on the Great Observatories Origins Deep Survey South (GOODS-S) field. We provide estimates of number counts at the three BLAST wavelengths, 250, 350, and 500 microns; instead of counting sources in flux bins we estimate the counts at several flux density nodes connected with power-laws. We observe a generally very steep slope for the counts of about -3.7 at 250 microns and -4.5 at 350 and 500 microns, over the range ~0.02-0.5 Jy, breaking to a shallower slope below about 0.015 Jy at all three wavelengths. We also describe how to estimate the uncertainties and correlations in this method so that the results can be used for model-fitting. This method should be well-suited for analysis of data from the Herschel satellite.

We detect correlations in the cosmic far-infrared background due to the clustering of star-forming galaxies in observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST, at 250, 350, and 500 microns. We perform jackknife and other tests to confirm the reality of the signal. The measured correlations are well fit by a power law over scales of 5-25 arcminutes, with Delta I/I = 15.1 +/- 1.7%. We adopt a specific model for submillimeter sources in which the contribution to clustering comes from sources in the redshift ranges 1.3 <= z <= 2.2, 1.5 <= z <= 2.7, and 1.7 <= z <= 3.2, at 250, 350, and 500 microns, respectively. With these distributions, our measurement of the power spectrum, P(k_theta), corresponds to linear bias parameters, b = 3.8 +/- 0.6, 3.9 +/- 0.6 and 4.4 +/- 0.7, respectively. We further interpret the results in terms of the halo model, and find that at the smaller scales, the simplest halo model fails to fit our results. One way to improve the fit is to increase the radius at which dark matter halos are artificially truncated in the model, which is equivalent to having some star-forming galaxies at z >= 1 located in the outskirts of groups and clusters. In the context of this model we find a minimum halo mass required to host a galaxy is log (M_min / M_sun) = 11.5 (+0.4/-0.1), and we derive effective biases $b_eff = 2.2 +/- 0.2, 2.4 +/- 0.2, and 2.6 +/- 0.2, and effective masses log (M_eff / M_sun) = 12.9 +/- 0.3, 12.8 +/- 0.2, and 12.7 +/- 0.2, at 250, 350, and 500 microns, corresponding to spatial correlation lengths of r_0 = 4.9, 5.0, and 5.2 +/- 0.7 h^-1 Mpc, respectively. Finally, we discuss implications for clustering measurement strategies with Herschel and Planck.

We report multi-wavelength power spectra of diffuse Galactic dust emission from BLAST observations at 250, 350, and 500 microns in Galactic Plane fields in Cygnus X and Aquila. These submillimeter power spectra statistically quantify the self-similar structure observable over a broad range of scales and can be used to assess the cirrus noise which limits the detection of faint point sources. The advent of submillimeter surveys with the Herschel Space Observatory makes the wavelength dependence a matter of interest. We show that the observed relative amplitudes of the power spectra can be related through a spectral energy distribution (SED). Fitting a simple modified black body to this SED, we find the dust temperature in Cygnus X to be 19.9 +/- 1.3 K and in the Aquila region 16.9 +/- 0.7 K. Our empirical estimates provide important new insight into the substantial cirrus noise that will be encountered in forthcoming observations.

Over the course of two flights, the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) made resolved maps of seven nearby (<25 Mpc) galaxies at 250, 350, and 500 microns. During its June 2005 flight from Sweden, BLAST observed a single nearby galaxy, NGC 4565. During the December 2006 flight from Antarctica, BLAST observed the nearby galaxies NGC 1097, NGC 1291, NGC 1365, NGC 1512, NGC 1566, and NGC 1808. We fit physical dust models to a combination of BLAST observations and other available data for the galaxies observed by Spitzer. We fit a modified blackbody to the remaining galaxies to obtain total dust mass and mean dust temperature. For the four galaxies with Spitzer data, we also produce maps and radial profiles of dust column density and temperature. We measure the fraction of BLAST detected flux originating from the central cores of these galaxies and use this to calculate a "core fraction," an upper limit on the "AGN fraction" of these galaxies. We also find our resolved observations of these galaxies give a dust mass estimate 5-19 times larger than an unresolved observations would predict. Finally, we are able to use these data to derive a value for the dust mass absorption co-efficient of kappa = 0.29 +/-0.03 m^2 kg^-1 at 250 microns. This study is an introduction to future higher-resolution and higher-sensitivity studies to be conducted by Herschel and SCUBA-II.

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

We present the 250, 350, and 500 micron detection of bright submillimeter emission in the direction of the Bullet Cluster measured by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The 500 micron centroid is coincident with an AzTEC 1.1 mm point-source detection at a position close to the peak lensing magnification produced by the cluster. However, the 250 micron and 350 micron centroids are elongated and shifted toward the south with a differential shift between bands that cannot be explained by pointing uncertainties. We therefore conclude that the BLAST detection is likely contaminated by emission from foreground galaxies associated with the Bullet Cluster. The submillimeter redshift estimate based on 250-1100 micron photometry at the position of the AzTEC source is z_phot = 2.9 (+0.6 -0.3), consistent with the infrared color redshift estimation of the most likely IRAC counterpart. These flux densities indicate an apparent far-infrared luminosity of L_FIR = 2E13 Lsun. When the amplification due to the gravitational lensing of the cluster is removed, the intrinsic far-infrared luminosity of the source is found to be L_FIR <= 10^12 Lsun, consistent with typical luminous infrared galaxies.

We present first results from an unbiased 50 deg^2 submillimeter Galactic survey at 250, 350, and 500 micron from the 2006 flight of the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The map has resolution ranging from 36 arcsec to 60 arcsec in the three submillimeter bands spanning the thermal emission peak of cold starless cores. We determine the temperature, luminosity, and mass of more than one thousand compact sources in a range of evolutionary stages and an unbiased statistical characterization of the population. From comparison with C^(18)O data, we find the dust opacity per gas mass, kappa r = 0.16 cm^2 g^(-1) at 250 micron, for cold clumps. We find that 2% of the mass of the molecular gas over this diverse region is in cores colder than 14 K, and that the mass function for these cold cores is consistent with a power law with index alpha = -3.22 +/- 0.14 over the mass range 14 M_sun < M < 80 M_sun. Additionally, we infer a mass-dependent cold core lifetime of t_c(M) = 4E6 (M/20 M_sun)^(-0.9) years - longer than what has been found in previous surveys of either low or high mass cores, and significantly longer than free fall or likely turbulent decay times. This implies some form of non-thermal support for cold cores during this early stage of star formation.

We directly measure redshift evolution in the mean physical properties (far-infrared luminosity, temperature, and mass) of the galaxies that produce the cosmic infrared background (CIB), using measurements from the Balloon-borne Large Aperture Sub-millimeter Telescope (BLAST), and Spitzer which constrain the CIB emission peak. This sample is known to produce a surface brightness in the BLAST bands consistent with the full CIB, and photometric redshifts are identified for all of the objects. We find that most of the 70 micron background is generated at z <~ 1 and the 500 micron background generated at z >~ 1. A significant growth is observed in the mean luminosity from ~ 10^9 - 10^12 L_sun, and in the mean temperature by 10 K, from redshifts 0< z < 3. However, there is only weak positive evolution in the comoving dust mass in these galaxies across the same redshift range. We also measure the evolution of the far-infrared luminosity density, and the star-formation rate history for these objects, finding good agreement with other infrared studies up to z ~1, exceeding the contribution attributed to optically-selected galaxies.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made one square degree, deep, confusion limited maps at three different bands, centered on the Great Observatories Origins Deep Survey South field. By calculating the covariance of these maps with catalogs of 24 micron sources from the Far-Infrared Deep Extragalactic Legacy Survey (FIDEL), we have determined that the total submillimeter intensities are 8.60 +/- 0.59, 4.93 +/- 0.34, and 2.27 +/- 0.20 nW m^2 sr^(-1) at 250, 350, and 500 micron, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 micron-selected galaxies generating the full intensity of the CIB. We find that the fraction of the CIB that originates from sources at z >= 1.2 increases with wavelength, with 60% from high redshift sources at 500 micron. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 micron-faint sources than it is for 24 micron-bright sources. Galaxies identified as active galactic nuclei (AGN) by their Infrared Array Camera (IRAC) colors are 1.6-2.6 times brighter than the average population at 250-500 micron, consistent with what is found for X-ray-selected AGN. BzK-selected galaxies are found to be moderately brighter than typical 24 micron-selected galaxies in the BLAST bands. These data provide high precision constraints for models of the evolution of the number density and intensity of star forming galaxies at high redshift.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 250-hour flight over Antarctica in December 2006 (BLAST06). As part of the calibration and pointing procedures, the red hypergiant star VY CMa was observed and used as the primary calibrator. Details of the overall BLAST06 calibration procedure are discussed. The 1-sigma absolute calibration is accurate to 10, 12, and 13% at the 250, 350, and 500 micron bands, respectively. The errors are highly correlated between bands resulting in much lower error for the derived shape of the 250-500 micron continuum. The overall pointing error is <5" rms for the 36, 42, and 60" beams. The performance of the optics and pointing systems is discussed.

We have identified radio and/or mid-infrared counterparts to 198 out of 350 sources detected at >=5 sigma over ~ 9 square degrees centered on the Chandra Deep Field South (CDFS) by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at 250, 350 and 500 um. We have matched 114 of these counterparts to optical sources with previously derived photometric redshifts and fitted SEDs to the BLAST fluxes and fluxes at 70 and 160 um acquired with the Spitzer Space Telescope. In this way, we have constrained dust temperatures, total far-infrared/sub-millimeter luminosities and star formation rates for each source. Our findings show that on average, the BLAST sources lie at significantly lower redshifts and have significantly lower rest-frame dust temperatures compared to submm sources detected in surveys conducted at 850 um. We demonstrate that an apparent increase in dust temperature with redshift in our sample arises as a result of selection effects. Finally, we provide the full multi-wavelength catalog of >= 5 sigma BLAST sources contained within the complete ~ 9 square degree survey area.

We present results from a multi-wavelength study of 29 sources (false detection probabilities <5%) from a survey of the Great Observatories Origins Deep Survey-North field at 1.1mm using the AzTEC camera. Comparing with existing 850um SCUBA studies in the field, we examine differences in the source populations selected at the two wavelengths. The AzTEC observations uniformly cover the entire survey field to a 1-sigma depth of ~1mJy. Searching deep 1.4GHz VLA, and Spitzer 3--24um catalogues, we identify robust counterparts for 21 1.1mm sources, and tentative associations for the remaining objects. The redshift distribution of AzTEC sources is inferred from available spectroscopic and photometric redshifts. We find a median redshift of z=2.7, somewhat higher than z=2.0 for 850um-selected sources in the same field, and our lowest redshift identification lies at a spectroscopic redshift z=1.1460. We measure the 850um to 1.1mm colour of our sources and do not find evidence for `850um dropouts', which can be explained by the low-SNR of the observations. We also combine these observed colours with spectroscopic redshifts to derive the range of dust temperatures T, and dust emissivity indices $\beta$ for the sample, concluding that existing estimates T~30K and $\beta$~1.75 are consistent with these new data.

Empirical "backward" galaxy evolution models for infrared bright galaxies are constrained using multi-band infrared surveys. We developed a new Monte-Carlo algorithm for this task, implementing luminosity dependent distribution functions for the galaxies' infrared spectral energy distributions (SEDs) and for the AGN contribution, allowing for evolution of these quantities. The adopted SEDs take into account the contributions of both starbursts and AGN to the infrared emission, for the first time in a coherent treatment rather than invoking separate AGN and star-forming populations. In the first part of the paper we consider the quantification of the AGN contribution for local universe galaxies, as a function of total infrared luminosity. It is made using a large sample of LIRGs and ULIRGs for which mid-infrared spectra are available in the Spitzer archive. In the second part we present the model. Our best-fit model adopts very strong luminosity evolution, $L=L_0(1+z)^{3.4}$, up to $z=2.3$, and density evolution, $\rho=\rho_0(1+z)^2$, up to $z=1$, for the population of infrared galaxies. At higher $z$, the evolution rates drop as $(1+z)^{-1}$ and $(1+z)^{-1.5}$ respectively. To reproduce mid-infrared to submillimeter number counts and redshift distributions, it is necessary to introduce both an evolution in the AGN contribution and an evolution in the luminosity-temperature relation. Our models are in plausible agreement with current photometry-based estimates of the typical AGN contribution as a function of mid-infrared flux, and well placed to be compared to upcoming Spitzer spectroscopic results. As an example of future applications, we use our best-fitting model to make predictions for surveys with Herschel.

Submillimetre surveys during the past decade have discovered a population of luminous, high-redshift, dusty starburst galaxies. In the redshift range 1 <= z <= 4, these massive submillimetre galaxies go through a phase characterized by optically obscured star formation at rates several hundred times that in the local Universe. Half of the starlight from this highly energetic process is absorbed and thermally re-radiated by clouds of dust at temperatures near 30 K with spectral energy distributions peaking at 100 microns in the rest frame. At 1 <= z <= 4, the peak is redshifted to wavelengths between 200 and 500 microns. The cumulative effect of these galaxies is to yield extragalactic optical and far-infrared backgrounds with approximately equal energy densities. Since the initial detection of the far-infrared background (FIRB), higher-resolution experiments have sought to decompose this integrated radiation into the contributions from individual galaxies. Here we report the results of an extragalactic survey at 250, 350 and 500 microns. Combining our results at 500 microns with those at 24 microns, we determine that all of the FIRB comes from individual galaxies, with galaxies at z >= 1.2 accounting for 70 per cent of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.

The steep source counts and negative K-corrections of bright submillimetre galaxies (SMGs) suggest that a significant fraction of those observed at high flux densities may be gravitationally lensed, and that the lensing objects may often lie at redshifts above 1, where clusters of galaxies are difficult to detect through other means. In this case follow-up of bright SMGs may be used to identify dense structures along the line of sight. Here we investigate the probability for SMGs to experience strong lensing, using the latest N-body simulations and observed source flux and redshift distributions. We find that almost all high redshift sources with a flux density above 100mJy will be lensed, if they are not relatively local galaxies. We also give estimates of the fraction of sources experiencing strong lensing as a function of observed flux density. This has implications for planning follow-up observations for bright SMGs discovered in future surveys with SCUBA-2 and other instruments. The largest uncertainty in these calculations is the maximum allowed lensing amplification, which is dominated by the presently unknown spatial extent of SMGs.

We measure the local galaxy far-infrared (FIR) 60-to-100 um colour-luminosity distribution using an all-sky IRAS survey. This distribution is an important reference for the next generation of FIR--submillimetre surveys that have and will conduct deep extra-galactic surveys at 250--500 um. With the peak in dust-obscured star-forming activity leading to present-day giant ellipticals now believed to occur in sub-mm galaxies near z~2.5, these new FIR--submillimetre surveys will directly sample the SEDs of these distant objects at rest-frame FIR wavelengths similar to those at which local galaxies were observed by IRAS. We have taken care to correct for temperature bias and evolution effects in our IRAS 60 um-selected sample. We verify that our colour-luminosity distribution is consistent with measurements of the local FIR luminosity function, before applying it to the higher-redshift Universe. We compare our colour-luminosity correlation with recent dust-temperature measurements of sub-mm galaxies and find evidence for pure luminosity evolution of the form (1+z)^3. This distribution will be useful for the development of evolutionary models for BLAST and SPIRE surveys as it provides a statistical distribution of rest-frame dust temperatures for galaxies as a function of luminosity.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital experiment designed to study the process of star formation in local galaxies (including the Milky Way) and in galaxies at cosmological distances. Using a 2-m Cassegrain telescope, BLAST images the sky onto a focal plane, which consists of 270 bolometric detectors split between three arrays, observing simultaneously in 30% wide bands, centered at 250, 350, and 500 microns. The diffraction-limited optical system provides a resolution of 30" at 250 microns. The pointing system enables raster-like scans with a positional accuracy of ~30", reconstructed to better than 5" rms in post-flight analysis. BLAST had two successful flights, from the Arctic in 2005, and from Antarctica in 2006, which provided the first high-resolution and large-area (~0.8-200 deg^2) submillimeter surveys at these wavelengths. As a pathfinder for the SPIRE instrument on Herschel, BLAST shares with the ESA satellite similar focal plane technology and scientific motivation. A third flight in 2009 will see the instrument modified to be polarization-sensitive (BLAST-Pol). With its unprecedented mapping speed and resolution, BLAST-Pol will provide insights into Galactic star-forming nurseries, and give the necessary link between the larger, coarse resolution surveys and the narrow, resolved observations of star-forming structures from space and ground based instruments being commissioned in the next 5 years.

We have conducted a deep and uniform 1.1 mm survey of the GOODS-N field with AzTEC on the James Clerk Maxwell Telescope (JCMT). Here we present the first results from this survey including maps, the source catalogue, and 1.1 mm number-counts. The results presented here were obtained from a 245 sq-arcmin region with near uniform coverage to a depth of 0.96-1.16 mJy/beam. Our robust catalogue contains 28 source candidates detected with S/N >= 3.75, only 1-2 of which are expected to be spurious detections. Of these source candidates, 8 are also detected by SCUBA at 850 um in regions where there is good overlap between the two surveys. The major advantage of our survey over that with SCUBA is the uniformity of coverage. We calculate number counts using two different techniques: the first using a frequentist parameter estimation, and the second using a Bayesian method. The two sets of results are in good agreement. We find that the 1.1 mm differential number counts are well described in the 2-6 mJy range by the functional form dN/dS = N' (S'/S) exp(-S/S') with fitted parameters S' = 1.25 +/-0.38 mJy and dN/dS = 300 +/- 90 per mJy per sq-deg at 3 mJy.

We present high-resolution interferometric imaging of LH850.02, the brightest 850- and 1200-micron submillimetre (submm) galaxy in the Lockman Hole. Our observations were made at 890 micron with the Submillimetre Array (SMA). Our high-resolution submm imaging detects LH850.02 at >6-sigma as a single compact (size < 1 arcsec or < 8 kpc) point source and yields its absolute position to ~0.2-arcsec accuracy. LH850.02 has two alternative radio counterparts within the SCUBA beam (LH850.02N & S), both of which are statistically very unlikely to be so close to the SCUBA source position by chance. However, the precise astrometry from the SMA shows that the submm emission arises entirely from LH850.02N, and is not associated with LH850.02S (by far the brighter of the two alternative identifications at 24-micron). Fits to the optical-infrared multi-colour photometry of LH850.02N & S indicate that both lie at z~3.3, and are therefore likely to be physically associated. At these redshifts, the 24 micron--to--submm flux density ratios suggest that LH850.02N has an Arp220-type starburst-dominated far-IR SED, while LH850.02S is more similar to Mrk231, with less dust-enshrouded star-formation activity, but a significant contribution at 24-micron (rest-frame ~5-6 micron) from an active nucleus. This complex mix of star-formation and AGN activity in multi-component sources may be common in the high redshift ultraluminous galaxy population, and highlights the need for precise astrometry from high resolution interferometric imaging for a more complete understanding.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 100-hour flight from northern Sweden in June 2005 (BLAST05). As part of the calibration and pointing procedures, several compact sources were mapped, including solar system, Galactic, and extragalactic targets, specifically Pallas, CRL 2688, LDN 1014, IRAS 20126+4104, IRAS 21078+5211, IRAS 21307+5049, IRAS 22134+5834, IRAS 23011+6126, K3-50, W 75N, and Mrk 231. One additional source, Arp 220, was observed and used as our primary calibrator. Details of the overall BLAST05 calibration procedure are discussed here. The BLAST observations of each compact source are described, flux densities and spectral energy distributions are reported, and these are compared with previous measurements at other wavelengths. The 250, 350, and 500 um BLAST data can provide useful constraints to the amplitude and slope of the submillimeter continuum, which in turn may be useful for the improved calibration of other submillimeter instruments.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital surveying experiment designed to study the evolutionary history and processes of star formation in local galaxies (including the Milky Way) and galaxies at cosmological distances. The BLAST continuum camera, which consists of 270 detectors distributed between 3 arrays, observes simultaneously in broad-band (30%) spectral-windows at 250, 350, and 500 microns. The optical design is based on a 2m diameter telescope, providing a diffraction-limited resolution of 30" at 250 microns. The gondola pointing system enables raster mapping of arbitrary geometry, with a repeatable positional accuracy of ~30"; post-flight pointing reconstruction to ~5" rms is achieved. The on-board telescope control software permits autonomous execution of a pre-selected set of maps, with the option of manual override. In this paper we describe the primary characteristics and measured in-flight performance of BLAST. BLAST performed a test-flight in 2003 and has since made two scientifically productive long-duration balloon flights: a 100-hour flight from ESRANGE (Kiruna), Sweden to Victoria Island, northern Canada in June 2005; and a 250-hour, circumpolar-flight from McMurdo Station, Antarctica, in December 2006.

We describe a map-making method which we have developed for the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) experiment, but which should have general application to data from other submillimeter arrays. Our method uses a Maximum Likelihood based approach, with several approximations, which allows images to be constructed using large amounts of data with fairly modest computer memory and processing requirements. This new approach, Signal And Noise Estimation Procedure Including Correlations (SANEPIC), builds upon several previous methods, but focuses specifically on the regime where there is a large number of detectors sampling the same map of the sky, and explicitly allowing for the the possibility of strong correlations between the detector timestreams. We provide real and simulated examples of how well this method performs compared with more simplistic map-makers based on filtering. We discuss two separate implementations of SANEPIC: a brute-force approach, in which the inverse pixel-pixel covariance matrix is computed; and an iterative approach, which is much more efficient for large maps. SANEPIC has been successfully used to produce maps using data from the 2005 BLAST flight.

We present the first results from a new 250, 350, and 500 micron Galactic Plane survey taken with the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) in 2005. This survey's primary goal is to identify and characterize high-mass proto-stellar objects (HMPOs). The region studied here covers 4 sq. deg near the open cluster NGC 6823 in the constellation Vulpecula (l=59). We find 60 compact sources (<60'' diameter) detected simultaneously in all three bands. Their spectral energy distributions (SEDs) are constrained through BLAST, IRAS, Spitzer MIPS, and MSX photometry, with inferred dust temperatures spanning ~12-40K assuming a dust emissivity index beta=1.5. The luminosity-to-mass ratio, a distance-independent quantity, spans ~0.2-130 L_\odot M_\odot^{-1}. Distances are estimated from coincident 13CO (1->0) velocities combined with a variety of other velocity and morphological data in the literature. In total, 49 sources are associated with a molecular cloud complex encompassing NGC 6823 (distance ~2.3kpc), 10 objects with the Perseus Arm (~8.5kpc) and one object is probably in the outer Galaxy (~14kpc). Near NGC 6823, the inferred luminosities and masses of BLAST sources span ~40-10^4 L_\odot, and ~15-700 M_\odot, respectively. The mass spectrum is compatible with molecular gas masses in other high-mass star forming regions. Several luminous sources appear to be Ultra Compact HII regions powered by early B stars. However, many of the objects are cool, massive gravitationally-bound clumps with no obvious internal radiation from a protostar, and hence excellent HMPO candidates.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

Using the 100-m Green Bank Telescope, we have conducted a cm-wavelength search for CO J=1-0 line emission towards the high-redshift, far-infrared luminous object, HDF850.1 over the redshift interval 3.3<z<5.4. Despite the wealth of existing multi-wavelength observations, and the recent identification of a galaxy counterpart in deep K' band (2.2 um) imaging, an unambiguous spectroscopic redshift has not yet been obtained for this object. A far-infrared-to-radio wavelength photometric redshift technique however, predicts a ~90% probability that the redshift is in the range, 3.3<z<5.4 (equivalent to an observed redshifted CO J=1-0 emission line frequency, 26.5>nu(obs)>18.0 GHz), making HDF850.1 a potential occupent of the `high-redshift tail' of submm selected galaxies. We have also conducted a search for CO J=2-1 line emission over the narrower redshift range, 3.9<z<4.3. although we do not detect any CO line emission in this object, our limits to the CO line luminosity are in broad agreement with the median value measured in the current sample of high-redshift, submm selected objects detected in high-J CO line emission, but not sufficient to fully test the validity of the photometric redshift technique.

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]

SCUBA-2 is an innovative 10,000 pixel submillimeter camera due to be delivered to the James Clerk Maxwell Telescope in late 2006. The camera is expected to revolutionize submillimeter astronomy in terms of the ability to carry out wide-field surveys to unprecedented depths addressing key questions relating to the origins of galaxies, stars and planets. This paper presents an update on the project with particular emphasis on the laboratory commissioning of the instrument. The assembly and integration will be described as well as the measured thermal performance of the instrument. A summary of the performance results will be presented from the TES bolometer arrays, which come complete with in-focal plane SQUID amplifiers and multiplexed readouts, and are cooled to 100mK by a liquid cryogen-free dilution refrigerator. Considerable emphasis has also been placed on the operating modes of the instrument and the common-user aspect of the user interface and data reduction pipeline. These areas will also be described in the paper.

We have developed two redundant daytime star cameras to provide the fine pointing solution for the balloon-borne submillimeter telescope, BLAST. The cameras are capable of providing a reconstructed pointing solution with an absolute accuracy < 5 arcseconds. They are sensitive to stars down to magnitudes ~ 9 in daytime float conditions. Each camera combines a 1 megapixel CCD with a 200 mm f/2 lens to image a 2 degree x 2.5 degree field of the sky. The instruments are autonomous. An internal computer controls the temperature, adjusts the focus, and determines a real-time pointing solution at 1 Hz. The mechanical details and flight performance of these instruments are presented.

In this paper we apply our Monte-Carlo photometric-redshift technique, introduced in paper I (Hughes et al. 2002), to the multi-wavelength data available for 77 galaxies selected at 850um and 1.25mm. We calculate a probability distribution for the redshift of each galaxy, which includes a detailed treatment of the observational errors and uncertainties in the evolutionary model. The cumulative redshift distribution of the submillimetre galaxy population that we present in this paper, based on 50 galaxies found in wide-area SCUBA surveys, is asymmetric, and broader than those published elsewhere, with a significant high-z tail for some of the evolutionary models considered. Approximately 40 to 90 per cent of the sub-mm population is expected to have redshifts in the interval 2 < z < 4. Whilst this result is completely consistent with earlier estimates for the sub-mm galaxy population, we also show that the colours of many (< 50 per cent) individual sub-mm sources, detected only at 850um with non-detections at other wavelengths, are consistent with those of starburst galaxies that lie at extreme redshifts, z > 4. Spectroscopic confirmation of the redshifts, through the detection of rest-frame FIR--mm wavelength molecular transition-lines, will ultimately calibrate the accuracy of this technique. We use the redshift probability distribution of HDF850.1 to illustrate the ability of the method to guide the choice of possible frequency tunings on the broad-band spectroscopic receivers that equip the large aperture single-dish mm and cm-wavelength telescopes.

Future extragalactic sub-millimetre and millimetre surveys have the potential to provide a sensitive census of the level of obscured star formation in galaxies at all redshifts. While in general there is good agreement between the source counts from existing SCUBA (850um) and MAMBO (1.25mm) surveys of different depths and areas, it remains difficult to determine the redshift distribution and bolometric luminosities of the sub-millimetre and millimetre galaxy population. This is principally due to the ambiguity in identifying an individual sub-millimetre source with its optical, IR or radio counterpart which, in turn, prevents a confident measurement of the spectroscopic redshift. Additionally, the lack of data measuring the rest-frame FIR spectral peak of the sub-millimetre galaxies gives rise to poor constraints on their rest-frame FIR luminosities and star formation rates. In this paper we describe Monte-Carlo simulations of ground-based, balloon-borne and satellite sub-millimetre surveys that demonstrate how the rest-frame FIR-sub-millimetre spectral energy distributions (250-850um) can be used to derive photometric redshifts with an r.m.s accuracy of +/- 0.4 over the range 0 < z < 6. This opportunity to break the redshift deadlock will provide an estimate of the global star formation history for luminous optically-obscured galaxies [L(FIR) > 3 x 10^12 Lsun] with an accuracy of 20 per cent.

More than 150 galaxies have been detected in blank-field millimetre and sub-millimetre surveys. However the redshift distribution of sub-mm galaxies remains uncertain due to the difficulty in identifying their optical-IR counterparts, and subsequently obtaining their spectroscopic emission-line redshifts. In this paper we discuss results from a Monte-Carlo analysis of the accuracy with which one can determine redshifts from photometric measurements at sub-millimetre-FIR wavelengths. The analysis takes into account the dispersion in colours introduced by including galaxies with a distribution of SEDs, and by including photometric and absolute calibration errors associated with real observations. We present examples of the probability distribution of redshifts for individual galaxies detected in the future BLAST and Herschel/SPIRE surveys. We show that the combination of BLAST and 850um observations constrain the photometric redshifts with sufficient accuracy to pursue a program of spectroscopic follow-up with the 100m GBT.

A robust measurement of the clustering amplitude of the sub-mm population of starburst galaxies requires large-area surveys (>> 1 deg^2). The largest-format arrays subtend only 10 arcmin^2 on the sky and hence scan-mapping is a necessary observing mode. Providing realistic representations of the extragalactic sky and atmosphere, as the input to a detailed simulator of the telescope and instrument performance, allows important decisions to be made about the design of large-area fully-sampled surveys and observing strategies. In this paper we present preliminary simulations that include detector noise, time-constants and array geometry, telescope pointing errors, scan speeds and scanning angles, sky noise and sky rotation.