We present 90, 140, and 268GHz sub-arcminute resolution imaging of the Sunyaev-Zel'dovich effect (SZE) in MACSJ0717.5+3745. Our 90GHz SZE data result in a sensitive, 34uJy/bm map at 13" resolution using MUSTANG. Our 140 and 268GHz SZE imaging, with resolutions of 58" and 31" and sensitivities of 1.8 and 3.3mJy/beam respectively, was obtained using Bolocam. We compare these maps to a 2-dimensional pressure map derived from Chandra X-ray observations. Our MUSTANG data confirm previous indications from Chandra of a pressure enhancement due to shock-heated, >20keV gas immediately adjacent to extended radio emission seen in low-frequency radio maps. The MUSTANG data also detect pressure substructure that is not well-constrained by the X-ray data in the remnant core of a merging subcluster. We find that the small-scale pressure enhancements in the MUSTANG data amount to ~2% of the total pressure measured in the 140GHz Bolocam observations. The X-ray template also fails on larger scales to accurately describe the Bolocam data, particularly at the location of a subcluster known to have a high line of sight optical velocity (~3200km/s). Our Bolocam data are adequately described when we add an additional component - not described by a thermal SZE spectrum - coincident with this subcluster. Using flux densities extracted from our model fits, and marginalizing over the temperature constraints for the region, we fit a thermal+kinetic SZE spectrum to our data and find the subcluster has a best-fit line of sight proper velocity of 3600+3440/-2160km/s. This agrees with the optical velocity estimates for the subcluster. The probability of velocity<0 given our measurements is 2.1%. Repeating this analysis using flux densities measured non-parametrically results in a 3.4% probability of a velocity<=0. We note that this tantalizing result for the kinetic SZE is on resolved, subcluster scales.

The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in the $K_{\rm s}$--3\,GHz range. Our method simultaneously fits the galaxies, the system foreground, and the leakage of flux from galaxies located in masked areas and corrects for an "over-fitting" effect not previously accounted for in stacking methods. We explore the contribution to the CIB as a function of galaxy survey wavelength and depth. We find high contributions to the CIB with the deep $r$ ($m_{\rm AB} \le 26.5$), $K_{\rm s}$ ($m_{\rm AB} \le 24.0$) and 3.6\,$\mu$m ($m_{\rm AB} \le 25.5$) catalogues. We combine these three deep catalogues and find a total CIB contributions of 10.5 $\pm$ 1.6, 6.7 $\pm$ 1.5 and 3.1 $\pm$ 0.7\,nWm$^{-2}$sr$^{-1}$ at 250, 350 and 500\,$\mu$m, respectively. Our CIB estimates are consistent with recent phenomenological models, prior based SPIRE number counts and with (though more precise than) the diffuse total measured by FIRAS. Our results raise the interesting prospect that the CIB contribution at $\lambda \le 500\,\mu$m from known galaxies has converged. Future large-area surveys like those with the Large Synoptic Survey Telescope are therefore likely to resolve a substantial fraction of the population responsible for the CIB at 250\,$\mu$m $\leq \lambda \leq$ 500\,$\mu$m.

We present measurements of the power spectra of cosmic infrared background (CIB) and cosmic microwave background (CMB) fluctuations in six frequency bands. Maps at the lower three frequency bands, 95, 150, and 220 GHz (3330, 2000, 1360 $\mu$m) are from the South Pole Telescope, while the upper three frequency bands, 600, 857, and 1200 GHz (500, 350, 250 $\mu$m) are observed with Herschel/SPIRE. From these data, we produce 21 angular power spectra (six auto- and fifteen cross-frequency) spanning the multipole range $600 \le \ell \le 11,000$. Our measurements are the first to cross-correlate measurements near the peak of the CIB spectrum with maps at 95 GHz, complementing and extending the measurements from Planck Collaboration et al. (2014) at 218, 550, and 857 GHz. The observed fluctuations originate largely from clustered, infrared-emitting, dusty star-forming galaxies, the CMB, and to a lesser extent radio galaxies, active galactic nuclei, and the Sunyaev-Zel'dovich effect.

We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (ICM), and HST imaging to derive models of each galaxy cluster's mass density. The galaxy clusters gravitationally lens the background CIB, which produced an on-average reduction in brightness towards the galaxy clusters' centers after the brightest galaxies were subtracted. We corrected for this deficit, which was between 5-25% of the 270 GHz SZ effect signal within $R_{2500}$. Using the SZ effect measurements, along with the X-ray constraint on $T_e$, we measured each galaxy cluster's average line of sight (LOS) velocity $v_z$ within $R_{2500}$, with a median per-cluster uncertainty of +-700 km/s. We found an ensemble-mean <$v_z$> of 430+-210 km/s, and an intrinsic cluster-to-cluster scatter $\sigma_{int}$ of 470+-340 km/s. We also obtained maps of $v_z$ over each galaxy cluster's face with an angular resolution of 70". All four galaxy clusters previously identified as having a merger oriented along the LOS showed an excess variance in these maps at a significance of 2-4$\sigma$, indicating an internal $v_z$ rms of $\gtrsim$1000 km/s. None of the six galaxy clusters previously identified as relaxed or plane of sky mergers showed any such excess variance.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Astrophysical measurements away from the 1 AU orbit of Earth can enable several astrophysical science cases that are challenging or impossible to perform from Earthbound platforms, including: building a detailed understanding of the extragalactic background light throughout the electromagnetic spectrum; measurements of the properties of dust and ice in the inner and outer solar system; determinations of the mass of planets and stellar remnants far from luminous stars using gravitational microlensing; and stable time-domain astronomy. Though potentially transformative for astrophysics, opportunities to fly instrumentation capable of these measurements are rare, and a mission to the distant solar system that includes instrumentation expressly designed to perform astrophysical science, or even one primarily for a different purpose but capable of precise astronomical investigation, has not yet been flown. In this White Paper, we describe the science motivations for this kind of measurement, and advocate for future flight opportunities that permit intersectional collaboration and cooperation to make these science investigations a reality.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

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

We quantify the star formation (SF) in the inner cores ($\mathcal{R}$/$R_{200}$$\leq$0.3) of 24 massive galaxy clusters at 0.2$\lesssim$$z$$\lesssim$0.9 observed by the $Herschel$ Lensing Survey and the Cluster Lensing and Supernova survey with $Hubble$. These programmes, covering the rest-frame ultraviolet to far-infrared regimes, allow us to accurately characterize stellar mass-limited ($\mathcal{M}_{*}$$>$$10^{10}$ $M_{\odot}$) samples of star-forming cluster members (not)-detected in the mid- and/or far-infrared. We release the catalogues with the photometry, photometric redshifts, and physical properties of these samples. We also quantify the SF displayed by comparable field samples from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We find that in intermediate-$z$ cluster cores, the SF activity is suppressed with respect the field in terms of both the fraction ($\mathcal{F}$) of star-forming galaxies (SFG) and the rate at which they form stars ($\mathcal{SFR}$ and $s\mathcal{SFR} = \mathcal{SFR}/\mathcal{M}_{*}$). On average, the $\mathcal{F}$ of SFGs is a factor $\sim$$2$ smaller in cluster cores than in the field. Furthermore, SFGs present average $\mathcal{SFR}$ and $s\mathcal{SFR}$ typically $\sim$0.3 dex smaller in the clusters than in the field along the whole redshift range probed. Our results favour long time-scale quenching physical processes as the main driver of SF suppression in the inner cores of clusters since $z$$\sim$0.9, with shorter time-scale processes being very likely responsible for a fraction of the missing SFG population.

The outer solar system provides a unique, quiet vantage point from which to observe the universe around us, where measurements could enable several niche astrophysical science cases that are too difficult to perform near Earth. NASA's New Horizons mission comprises an instrument package that provides imaging capability from ultraviolet (UV) to near-infrared (near-IR) wavelengths with moderate spectral resolution located beyond the orbit of Pluto. A carefully designed survey with New Horizons can optimize the use of expendable propellant and the limited data telemetry bandwidth to allow several measurements, including a detailed understanding of the cosmic extragalactic background light; studies of the local and extragalactic UV background; measurements of the properties of dust and ice in the outer solar system; confirmation and characterization of transiting exoplanets; determinations of the mass of dark objects using gravitational microlensing; and rapid follow-up of transient events. New Horizons is currently in an extended mission designed to focus on Kuiper Belt science that will conclude in 2021. The astrophysics community has a unique, generational opportunity to use this mission for astronomical observation at heliocentric distances beyond 50 au in the next decade. In this paper, we discuss the potential science cases for such an extended mission, and provide an initial assessment of the most important operational requirements and observation strategies it would require. We conclude that New Horizons is capable of transformative science, and that it would make a valuable and unique asset for astrophysical science that is unlikely to be replicated in the near future.

The near-infrared background between 0.5 $\mu$m to 2 $\mu$m contains a wealth of information related to radiative processes in the universe. Infrared background anisotropies encode the redshift-weighted total emission over cosmic history, including any spatially diffuse and extended contributions. The anisotropy power spectrum is dominated by undetected galaxies at small angular scales and diffuse background of Galactic emission at large angular scales. In addition to these known sources, the infrared background also arises from intra-halo light (IHL) at $z < 3$ associated with tidally-stripped stars during galaxy mergers. Moreover, it contains information on the very first galaxies from the epoch of reionization (EoR). The EoR signal has a spectral energy distribution (SED) that goes to zero near optical wavelengths due to Lyman absorption, while other signals have spectra that vary smoothly with frequency. Due to differences in SEDs and spatial clustering, these components may be separated in a multi-wavelength-fluctuation experiment. To study the extent to which EoR fluctuations can be separated in the presence of IHL, extra-galactic and Galactic foregrounds, we develop a maximum likelihood technique that incorporates a full covariance matrix among all the frequencies at different angular scales. We apply this technique to simulated deep imaging data over a 2$\times$100 deg$^2$ sky area from 0.75 $\mu$m to 5 $\mu$m in 9 bands and find that such a "frequency tomography" can successfully reconstruct both the amplitude and spectral shape for representative EoR, IHL and the foreground signals.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

Intensity mapping provides a unique means to probe the epoch of reionization (EoR), when the neutral intergalactic medium was ionized by the energetic photons emitted from the first galaxies. The [CII] 158$\mu$m fine-structure line is typically one of the brightest emission lines of star-forming galaxies and thus a promising tracer of the global EoR star-formation activity. However, [CII] intensity maps at $6 \lesssim z \lesssim 8$ are contaminated by interloping CO rotational line emission ($3 \leq J_{\rm upp} \leq 6$) from lower-redshift galaxies. Here we present a strategy to remove the foreground contamination in upcoming [CII] intensity mapping experiments, guided by a model of CO emission from foreground galaxies. The model is based on empirical measurements of the mean and scatter of the total infrared luminosities of galaxies at $z < 3$ and with stellar masses $M_{*} > 10^{8}\,\rm M_{\rm \odot}$ selected in $K$-band from the COSMOS/UltraVISTA survey, which can be converted to CO line strengths. For a mock field of the Tomographic Ionized-carbon Mapping Experiment (TIME), we find that masking out the "voxels" (spectral-spatial elements) containing foreground galaxies identified using an optimized CO flux threshold results in a $z$-dependent criterion $m^{\rm AB}_{\rm K} \lesssim 22$ (or $M_{*} \gtrsim 10^{9} \,\rm M_{\rm \odot}$) at $z < 1$ and makes a [CII]/CO$_{\rm tot}$ power ratio of $\gtrsim 10$ at $k=0.1$ $h$/Mpc achievable, at the cost of a moderate $\lesssim 8\%$ loss of total survey volume.

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).

We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light.

We present a demonstration of delensing the observed cosmic microwave background (CMB) B-mode polarization anisotropy. This process of reducing the gravitational-lensing generated B-mode component will become increasingly important for improving searches for the B modes produced by primordial gravitational waves. In this work, we delens B-mode maps constructed from multi-frequency SPTpol observations of a 90 deg$^2$ patch of sky by subtracting a B-mode template constructed from two inputs: SPTpol E-mode maps and a lensing potential map estimated from the $\textit{Herschel}$ $500\,\mu m$ map of the CIB. We find that our delensing procedure reduces the measured B-mode power spectrum by 28% in the multipole range $300 < \ell < 2300$; this is shown to be consistent with expectations from theory and simulations and to be robust against systematics. The null hypothesis of no delensing is rejected at $6.9 \sigma$. Furthermore, we build and use a suite of realistic simulations to study the general properties of the delensing process and find that the delensing efficiency achieved in this work is limited primarily by the noise in the lensing potential map. We demonstrate the importance of including realistic experimental non-idealities in the delensing forecasts used to inform instrument and survey-strategy planning of upcoming lower-noise experiments, such as CMB-S4.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

We present physical properties [redshifts (z), star-formation rates (SFRs) and stellar masses (Mstar)] of bright (S850>4mJy) submm galaxies in the ~2deg2 COSMOS and UDS fields selected with SCUBA-2/JCMT. We complete the galaxy identification process for all (~2000) S/N>3.5 850-um sources, but focus our scientific analysis on a high-quality sub-sample of 651 S/N>4 sources with complete multi-wavelength coverage including 1.1-mm imaging. We check the reliability of our identifications, and the robustness of the SCUBA-2 fluxes by revisiting the recent ALMA follow-up of 29 sources in our sample. Considering >4mJy ALMA sources, our identification method has a completeness of ~86 per cent with a reliability of ~92 per cent, and only ~15-20 per cent of sources are significantly affected by multiplicity (when a secondary component contributes >1/3 of the primary source flux). The impact of source blending on the 850-um source counts as determined with SCUBA-2 is modest; scaling the single-dish fluxes by ~0.9 reproduces the ALMA source counts. For our final SCUBA-2 sample we find median z=2.40+0.10-0.04, SFR=287+-6Moyr-1, and log(Mstar/Mo)=11.12+-0.02 (the latter for 349/651 sources with optical identifications). These properties clearly locate bright submm galaxies on the high-mass end of the 'main sequence' of star-forming galaxies out to z~6, suggesting that major mergers are not a dominant driver of the high-redshift submm-selected population. Their number densities are also consistent with the evolving galaxy stellar mass function. Hence, the submm galaxy population is as expected, albeit reproducing the evolution of the main sequence of star-forming galaxies remains a challenge for theoretical models/simulations.

The Extragalactic Background Light (EBL) captures the total integrated emission from stars and galaxies throughout the cosmic history. The amplitude of the near-infrared EBL from space absolute photometry observations has been controversial and depends strongly on the modeling and subtraction of the Zodiacal light foreground. We report the first measurement of the diffuse background spectrum at 0.8-1.7 um from the CIBER experiment. The observations were obtained with an absolute spectrometer over two flights in multiple sky fields to enable the subtraction of Zodiacal light, stars, terrestrial emission, and diffuse Galactic light. After subtracting foregrounds and accounting for systematic errors, we find the nominal EBL brightness, assuming the Kelsall Zodiacal light model, is 42.7+11.9/-10.6 nW/m2/sr at 1.4 um. We also analyzed the data using the Wright Zodiacal light model, which results in a worse statistical fit to the data and an unphysical EBL, falling below the known background light from galaxies at <1.3 um. Using a model-independent analysis based on the minimum EBL brightness, we find an EBL brightness of 28.7+5.1/-3.3 nW/m2/sr at 1.4 um. While the derived EBL amplitude strongly depends on the Zodiacal light model, we find that we cannot fit the spectral data to Zodiacal light, Galactic emission, and EBL from solely integrated galactic light from galaxy counts. The results require a new diffuse component, such as an additional foreground or an excess EBL with a redder spectrum than that of Zodiacal light.

The cosmic optical background is an important observable that constrains energy production in stars and more exotic physical processes in the universe, and provides a crucial cosmological benchmark against which to judge theories of structure formation. Measurement of the absolute brightness of this background is complicated by local foregrounds like the Earth's atmosphere and sunlight reflected from local interplanetary dust, and large discrepancies in the inferred brightness of the optical background have resulted. Observations from probes far from the Earth are not affected by these bright foregrounds. Here we analyze data from the Long Range Reconnaissance Imager (LORRI) instrument on NASA's New Horizons mission acquired during cruise phase outside the orbit of Jupiter, and find a statistical upper limit on the optical background's brightness similar to the integrated light from galaxies. We conclude that a carefully performed survey with LORRI could yield uncertainties comparable to those from galaxy counting measurements.

We present near-infrared (0.8-1.8 $\mu$m) spectra of 105 bright (${m_{J}}$ $<$ 10) stars observed with the low resolution spectrometer on the rocket-borne Cosmic Infrared Background Experiment (CIBER). As our observations are performed above the earth's atmosphere, our spectra are free from telluric contamination, which makes them a unique resource for near-infrared spectral calibration. Two-Micron All Sky Survey (2MASS) photometry information is used to identify cross-matched stars after reduction and extraction of the spectra. We identify the spectral types of the observed stars by comparing them with spectral templates from the Infrared Telescope Facility (IRTF) library. All the observed spectra are consistent with late F to M stellar spectral types, and we identify various infrared absorption lines.

Intensity mapping is now becoming a useful tool to study the large-scale structure of the universe through spatial variations in the integrated emission from galaxies and the intergalactic medium. We study intensity mapping of the H-alpha 6563AA, [OIII]5007AA, [OII]3727AA and H-beta 4861AA lines at 0.8<z<5.2. The mean intensities of these four emission lines are estimated using the observed luminosity functions (LFs), cosmological simulations, and the star formation rate density (SFRD) derived from observations at z<5. We calculate the intensity power spectra and consider the foreground contamination of other lines at lower redshifts. We use the proposed NASA small explorer SPHEREx (the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) as a case study for the detectability of the intensity power spectra of the four emission lines. We also investigate the cross correlation with the 21-cm line probed by CHIME (the Canadian Hydrogen Intensity Mapping Experiment), Tianlai experiment and SKA (the Square Kilometer Array) at 0.8<z<2.4. We find both the auto and cross power spectra can be well measured for the H-alpha, [OIII] and [OII] lines at z<3, while it is more challenging for the H-beta line. Finally, we estimate the constraint on the SFRD from intensity mapping, and find we can reach accuracy higher than 7% at z<4, which is better than usual measurements using the LFs of galaxies.

Measurement of the gas velocity distribution in galaxy clusters provides insight into the physics of mergers, through which large scale structures form in the Universe. Velocity estimates within the intracluster medium (ICM) can be obtained via the Sunyaev-Zel'dovich (SZ) effect, but its observation is challenging both in term of sensitivity requirement and control of systematic effects, including the removal of contaminants. In this paper we report resolved observations, at 150 and 260 GHz, of the SZ effect toward the triple merger MACS J0717.5+3745 (z=0.55), using data obtained with the NIKA camera at the IRAM 30m telescope. Assuming that the SZ signal is the sum of a thermal (tSZ) and a kinetic (kSZ) component and by combining the two NIKA bands, we extract for the first time a resolved map of the kSZ signal in a cluster. The kSZ signal is dominated by a dipolar structure that peaks at -5.1 and +3.4 sigma, corresponding to two subclusters moving respectively away and toward us and coincident with the cold dense X-ray core and a hot region undergoing a major merging event. We model the gas electron density and line-of-sight velocity of MACS J0717.5+3745 as four subclusters. Combining NIKA data with X-ray observations from XMM-Newton and Chandra, we fit this model to constrain the gas line-of-sight velocity of each component, and we also derive, for the first time, a velocity map from kSZ data (i.e. that is model-dependent). Our results are consistent with previous constraints on the merger velocities, and thanks to the high angular resolution of our data, we are able to resolve the structure of the gas velocity. Finally, we investigate possible contamination and systematic effects with a special care given to radio and submillimeter galaxies. Among the sources that we detect with NIKA, we find one which is likely to be a high redshift lensed submillimeter galaxy.

We present multi-wavelength detections of nine candidate gravitationally-lensed dusty star-forming galaxies (DSFGs) selected at 218GHz (1.4mm) from the ACT equatorial survey. Among the brightest ACT sources, these represent the subset of the total ACT sample lying in Herschel SPIRE fields, and all nine of the 218GHz detections were found to have bright Herschel counterparts. By fitting their spectral energy distributions (SEDs) with a modified blackbody model with power-law temperature distribution, we find the sample has a median redshift of $z=4.1^{+1.1}_{-1.0}$ (68 per cent confidence interval), as expected for 218GHz selection, and an apparent total infrared luminosity of $\log_{10}(\mu L_{\rm IR}/{\rm L}_\odot) = 13.86^{+0.33}_{-0.30}$, which suggests that they are either strongly lensed sources or unresolved collections of unlensed DSFGs. The effective apparent diameter of the sample is $\sqrt{\mu}d= 4.2^{+1.7}_{-1.0}$kpc, further evidence of strong lensing or multiplicity, since the typical diameter of dusty star-forming galaxies is $1.0$--$2.5$ kpc. We emphasize that the effective apparent diameter derives from SED modelling without the assumption of optically thin dust (as opposed to image morphology). We find that the sources have substantial optical depth ($\tau = 4.2^{+3.7}_{-1.9}$) to dust around the peak in the modified blackbody spectrum ($\lambda_{\rm obs} \le 500$ $\mu$m), a result that is robust to model choice.

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.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

Mega-pixel charge-integrating detectors are common in near-IR imaging applications. Optimal signal-to-noise ratio estimates of the photocurrents, which are particularly important in the low-signal regime, are produced by fitting linear models to sequential reads of the charge on the detector. Algorithms that solve this problem have a long history, but can be computationally intensive. Furthermore, the cosmic ray background is appreciable for these detectors in Earth orbit, particularly above the Earth's magnetic poles and the South Atlantic Anomaly, and on-board reduction routines must be capable of flagging affected pixels. In this paper we present an algorithm that generates optimal photocurrent estimates and flags random transient charge generation from cosmic rays, and is specifically designed to fit on a computationally restricted platform. We take as a case study the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer (SPHEREx), a NASA Small Explorer astrophysics experiment concept, and show that the algorithm can easily fit in the resource-constrained environment of such a restricted platform. Detailed simulations of the input astrophysical signals and detector array performance are used to characterize the fitting routines in the presence of complex noise properties and charge transients. We use both Hubble Space Telescope Wide Field Camera-3 and Wide-field Infrared Survey Explorer to develop an empirical understanding of the susceptibility of near-IR detectors in low earth orbit and build a model for realistic cosmic ray energy spectra and rates. We show that our algorithm generates an unbiased estimate of the true photocurrent that is identical to that from a standard line fitting package, and characterize the rate, energy, and timing of both detected and undetected transient events.

The extragalactic background light (EBL) is comprised of the cumulative radiation from all galaxies and active galactic nuclei over the cosmic history. In addition to point sources, EBL also contains information from diffuse sources of radiation. The angular power spectra of the near-infrared intensities could contain additional signals and a complete understanding of the nature of the IR background is still lacking in the literature. Here we explore the constraints that can be placed on particle decays, especially candidate dark matter models involving axions that trace dark matter halos of galaxies. Axions with a mass around a few eV will decay via two photons with wavelengths in the near-IR band, and will leave a signature in the IR background intensity power spectrum. Using recent power spectra measurements from the Hubble Space Telescope (HST) and Cosmic Infrared Background Experiment (CIBER), we find that the 0.6 to 1.6 micron power spectra can be explained by axions with masses around 4 eV. The total axion abundance Omega_a~0.05, and it is comparable to the baryon density of the Universe. The suggested mean axion mass and abundance are not ruled out by existing cosmological observations. Interestingly, the axion model with a mass distribution is preferred by the data, which cannot be explained by the standard quantum chromodynamics (QCD) theory and needs further discussion.

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

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

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

We report on the temporal and spatial fluctuations in the atmospheric brightness in the narrow band between Meinel emission lines at 1191.3 nm using an R=320 near-infrared instrument. We present the instrument design and implementation, followed by a detailed analysis of data taken over the course of a night from Table Mountain Observatory. The absolute sky brightness at this wavelength is found to be 5330 +/- 30 nW m^-2 sr^-1, consistent with previous measurements of the inter-band airglow at these wavelengths. This amplitude is larger than simple models of the continuum component of the airglow emission at these wavelengths, confirming that an extra emissive or scattering component is required to explain the observations. We perform a detailed investigation of the noise properties of the data and find no evidence for a noise component associated with temporal instability in the inter-line continuum. This result demonstrates that in several hours of ~100s integrations the noise performance of the instrument does not appear to significantly degrade from expectations, giving a proof of concept that near-IR line intensity mapping may be feasible from ground-based sites.

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

We report measurements of the Diffuse Galactic Light (DGL) spectrum in the near-infrared, spanning the wavelength range 0.95-1.65 {\mu}m by the Cosmic Infrared Background ExpeRiment (CIBER). Using the low-resolution spectrometer (LRS) calibrated for absolute spectro-photometry, we acquired long-slit spectral images of the total diffuse sky brightness towards four high-latitude fields spread over four sounding rocket flights. To separate the DGL spectrum from the total sky brightness, we correlated the spectral images with a 100 {\mu}m intensity map, which traces the dust column density in optically thin regions. The measured DGL spectrum shows no resolved features and is consistent with other DGL measurements in the optical and at near-infrared wavelengths longer than 1.8 {\mu}m. Our result implies that the continuum is consistently reproduced by models of scattered starlight in the Rayleigh scattering regime with a few large grains.

Extragalactic background light (EBL) anisotropy traces variations in the total production of photons over cosmic history, and may contain faint, extended components missed in galaxy point source surveys. Infrared EBL fluctuations have been attributed to primordial galaxies and black holes at the epoch of reionization (EOR), or alternately, intra-halo light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report new EBL anisotropy measurements from a specialized sounding rocket experiment at 1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely explained by IHL emission. The measured fluctuations are associated with an EBL intensity that is comparable to the background from known galaxies measured through number counts, and therefore a substantial contribution to the energy contained in photons in the cosmos.

We present a statistical analysis of the millimeter-wavelength properties of 1.4 GHz-selected sources and a detection of the Sunyaev-Zel'dovich (SZ) effect associated with the halos that host them. The Atacama Cosmology Telescope (ACT) has conducted a survey at 148 GHz, 218 GHz and 277 GHz along the celestial equator. Using samples of radio sources selected at 1.4 GHz from FIRST and NVSS, we measure the stacked 148, 218 and 277 GHz flux densities for sources with 1.4 GHz flux densities ranging from 5 to 200 mJy. At these flux densities, the radio source population is dominated by active galactic nuclei (AGN), with both steep and flat spectrum populations, which have combined radio-to-millimeter spectral indices ranging from 0.5 to 0.95, reflecting the prevalence of steep spectrum sources at high flux densities and the presence of flat spectrum sources at lower flux densities. The thermal SZ effect associated with the halos that host the AGN is detected at the 5$\sigma$ level through its spectral signature. When we compare the SZ effect with weak lensing measurements of radio galaxies, we find that the relation between the two is consistent with that measured by Planck for local bright galaxies. We present a detection of the SZ effect in some of the lowest mass halos (average $M_{200}\approx10^{13}$M$_{\odot}h_{70}^{-1}$) studied to date. This detection is particularly important in the context of galaxy evolution models, as it confirms that galaxies with radio AGN also typically support hot gaseous halos. With Herschel observations, we show that the SZ detection is not significantly contaminated by dust. We show that 5 mJy$<S_{1.4}<$200 mJy radio sources contribute $\ell(\ell+1)C_{\ell}/(2\pi)=0.37\pm0.03\mu$K$^2$ to the angular power spectrum at $\ell=3000$ at 148 GHz, after accounting for the SZ effect associated with their host halos.

We present the first cross-correlation measurement between Sloan Digital Sky Survey (SDSS) Type 1 quasars and the cosmic infrared background (CIB) measured by Herschel. The distribution of the quasars at 0.15<z<3.5 covers the redshift range where we expect most of the CIB to originate. We detect the sub-mm emission of the quasars, which dominates on small scales, as well as correlated emission from dusty star-forming galaxies (DSFGs) dominant on larger scales. The mean sub-mm flux densities of the DR7 quasars (median redshift <z>=1.4) is $11.1^{+1.6}_{-1.4}$, $7.1^{+1.6}_{-1.3}$ and $3.6^{+1.4}_{-1.0}$ mJy at 250, 350 and 500 microns, respectively, while the mean sub-mm flux densities of the DR9 quasars (<z>=2.5) is $5.7^{+0.7}_{-0.6}$, $5.0^{+0.8}_{-0.7}$ and $1.8^{+0.5}_{-0.4}$ mJy. We find that the correlated sub-mm emission includes both the emission from satellite DSFGs in the same halo as the central quasar and the emission from DSFGs in separate halos (correlated with the quasar-hosting halo). The amplitude of the one-halo term is ~10 times smaller than the sub-mm emission of the quasars, implying the the satellites have a lower star-formation rate than the quasars. The satellite fraction for the DR7 quasars is $0.008^{+0.008}_{-0.005}$ and the host halo mass scale for the central and satellite quasars is $10^{12.36\pm0.87}$ M$_{\odot}$ and $10^{13.60\pm0.38}$ M$_{\odot}$, respectively. The satellite fraction of the DR9 quasars is $0.065^{+0.021}_{-0.031}$ and the host halo mass scale for the central and satellite quasars is $10^{12.29\pm0.62}$ M$_{\odot}$ and $10^{12.82\pm0.39}$ M$_{\odot}$, respectively. Thus, the typical halo environment of the SDSS Type 1 quasars is found to be similar to that of DSFGs, which supports the generally accepted view that dusty starburst and quasar activity are evolutionarily linked.

We study the far-infrared properties of 498 Lyman Alpha Emitters (LAEs) at z=2.8, 3.1 and 4.5 in the Extended Chandra Deep Field-South, using 250, 350 and 500 micron data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and 870 micron data from the LABOCA ECDFS Submillimeter Survey (LESS). None of the 126, 280 or 92 LAEs at z=2.8, 3.1 and 4.5, respectively, are individually detected in the far-infrared data. We use stacking to probe the average emission to deeper flux limits, reaching $1\sigma$ depths of ~0.1 to 0.4 mJy. The LAEs are also undetected at $\ge3\sigma$ in the stacks, although a $2.5\sigma$ signal is observed at 870 micron for the z=2.8 sources. We consider a wide range of far-infrared spectral energy distributions (SEDs), including a M82 and an Sd galaxy template, to determine upper limits on the far-infrared luminosities and far-infrared-derived star-formation rates of the LAEs. These star-formation rates are then combined with those inferred from the Ly$\alpha$ and UV emission to determine lower limits on the LAEs Ly$\alpha$ escape fraction ($f_{\rm esc}($Ly$\alpha$)). For the Sd SED template, the inferred LAEs $f_{\rm esc}($Ly$\alpha$) are $\gtrsim30%$ ($1\sigma$) at z=2.8, 3.1 and 4.5, which are all significantly higher than the global $f_{\rm esc}($Ly$\alpha$) at these redshifts. Thus, if the LAEs $f_{\rm esc}($Ly$\alpha$) follows the global evolution then they have warmer far-infrared SEDs than the Sd galaxy template. The average and M82 SEDs produce lower limits on the LAE $f_{\rm esc}($Ly$\alpha$) of ~10 to 20% ($1\sigma$), all of which are slightly higher than the global evolution of $f_{\rm esc}($Ly$\alpha$) but consistent with it at the 2 to 3$\sigma$ level.

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

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

We report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km/s (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit line-of-sight velocity of v_z = +3450 +- 900 km/s, with (1 - Prob[v_z > 0]) = 1.3 x 10^-5 (4.2 sigma away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v_z = +2550 +- 1050 km/s, with (1 - Prob[v_z > 0]) = 2.2 x 10^-3 (2.9 sigma).

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

The potential for Planck to detect clusters of dusty, star-forming galaxies at z greater than 1 is tested by examining the Herschel-SPIRE images of Planck Early Release Compact Source Catalog (ERCSC) sources lying in fields observed by the HerMES survey. Of the 16 Planck sources that lie in the roughly 90 sq. deg. examined, we find that twelve are associated with single bright Herschel sources. The remaining four are associated with overdensities of Herschel sources, making them candidate clusters of dusty, starforming galaxies. We use complementary optical and NIR data for these clumps to test this idea, and find evidence for the presence of galaxy clusters in all four cases. We use photometric redshifts and red sequence galaxies to estimate the redshifts of these clusters, finding that they range from 0.8 to 2.3. These redshifts imply that the Herschel sources in these clusters, which contribute to the detected Planck flux, are forming stars very rapidly, with typical total cluster star formation rates greater than 1000Msun per yr. The high redshift clusters discovered in these observations are used to constrain the epoch of cluster galaxy formation, finding that the galaxies in our clusters are 1 to 1.5 Gy old at z about 1 to 2. Prospects for the discovery of further clusters of dusty galaxies are discussed, using not only all sky Planck surveys, but also deeper, smaller area, Herschel surveys.

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

We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin^2) K-selected sources (K_AB < 24.0), split according to their rest-frame U - V vs. V - J colors into 72,216 star-forming and 9,034 quiescent galaxies, on maps from Spitzer/MIPS (24um), Herschel/SPIRE (250, 350, 500um), Herschel/PACS (100, 160um), and AzTEC (1100um). The fraction of the CIB resolved by our catalog is (69 $\pm$ 15)% at 24um, (78 $\pm$ 17)% at 70um, (58 $\pm$ 13)% at 100um, (78 $\pm$ 18)% at 160um, (80 $\pm$ 17)% at 250um, (69 $\pm$ 14)% at 350um, (65 $\pm$ 12)% at 500um, and (45 $\pm$ 8)% at 1100um. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at $\lambda$ < 200um is sourced predominantly from galaxies at z < 1, while at $\lambda$ > 200um the bulk originates from 1 < z < 2. Galaxies with stellar masses log(M/ M_sun)=9.5-11 are responsible for the majority of the CIB, with those in the log(M/ M_sun)=9.5-10 contributing mostly at $\lambda$ < 250um, and those in the log(M/ M_sun)=10.5-11 bin dominating at $\lambda$ > 350um. The contribution from galaxies in the log(M/ M_sun)=9.0-9.5 and log(M/ M_sun)=11.0-12.0 stellar mass bins contribute the least, both of order 5%, although the highest stellar-mass bin is a significant contributor to the luminosity density at z > 2. The luminosities of the galaxies responsible for the CIB shifts from a combination of "normal" and luminous infrared galaxies (LIRGs) at $\lambda$ < 160um, to LIRGs at 160um < $\lambda$ < 500um, to finally LIRGs and ultra-luminous infrared galaxies (ULIRGs) at $\lambda$ > 500um. Stacking analyses were performed with SIMSTACK (available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html) which accounts for possible biases due to clustering.

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

We study the link between observed ultraviolet luminosity, stellar mass, and dust attenuation within rest-frame UV-selected samples at z~ 4, 3, and 1.5. We measure by stacking at 250, 350, and 500 um in the Herschel/SPIRE images from the HerMES program the average infrared luminosity as a function of stellar mass and UV luminosity. We find that dust attenuation is mostly correlated with stellar mass. There is also a secondary dependence with UV luminosity: at a given UV luminosity, dust attenuation increases with stellar mass, while at a given stellar mass it decreases with UV luminosity. We provide new empirical recipes to correct for dust attenuation given the observed UV luminosity and the stellar mass. Our results also enable us to put new constraints on the average relation between star formation rate and stellar mass at z~ 4, 3, and 1.5. The star formation rate-stellar mass relations are well described by power laws (SFR~ M^0.7), with the amplitudes being similar at z~4 and z~3, and decreasing by a factor of 4 at z~1.5 at a given stellar mass. We further investigate the evolution with redshift of the specific star formation rate. Our results are in the upper range of previous measurements, in particular at z~3, and are consistent with a plateau at 3<z<4. Current model predictions (either analytic, semi-analytic or hydrodynamic) are inconsistent with these values, as they yield lower predictions than the observations in the redshift range we explore. We use these results to discuss the star formation histories of galaxies in the framework of the Main Sequence of star-forming galaxies. Our results suggest that galaxies at high redshift (2.5<z<4) stay around 1 Gyr on the Main Sequence. With decreasing redshift, this time increases such that z=1 Main Sequence galaxies with 10^8<M_*/Msun<10^10 stay on the Main Sequence until z=0.

Gravitational lensing of the cosmic microwave background generates a curl pattern in the observed polarization. This "B-mode" signal provides a measure of the projected mass distribution over the entire observable Universe and also acts as a contaminant for the measurement of primordial gravity-wave signals. In this Letter we present the first detection of gravitational lensing B modes, using first-season data from the polarization-sensitive receiver on the South Pole Telescope (SPTpol). We construct a template for the lensing B-mode signal by combining E-mode polarization measured by SPTpol with estimates of the lensing potential from a Herschel-SPIRE map of the cosmic infrared background. We compare this template to the B modes measured directly by SPTpol, finding a non-zero correlation at 7.7 sigma significance. The correlation has an amplitude and scale-dependence consistent with theoretical expectations, is robust with respect to analysis choices, and constitutes the first measurement of a powerful cosmological observable.

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

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

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3{\sigma}). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

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

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

We describe the far-infrared (FIR; rest-frame 8--1000\mu m) properties of a sample of 443 H\alpha-selected star-forming galaxies in the COSMOS and UDS fields detected by the HiZELS imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select H\alpha\ (and [OII] if available) emitters in a narrow redshift slice at z = 1.47+/-0.02. We use a stacking approach in Spitzer, Herschel (from PEP and HerMES surveys) and AzTEC images to describe their typical FIR properties. We find that HiZELS galaxies with observed H\alpha\ luminosities of ~ 10^{8.1-9.1} Lo have bolometric FIR luminosities of typical LIRGs, L_FIR ~ 10^{11.48+/-0.05} Lo. Combining the H\alpha\ and FIR luminosities, we derive median SFR = 32+/-5 Mo/yr and H\alpha\ extinctions of A(H\alpha) = 1.0+/-0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to A(H\alpha) relations and the little or no evolution up to z = 1.47. For HiZELS galaxies, we provide an empirical parametrisation of the SFR as a function of (u-z)_rest colours and 3.6\mu m photometry. We find that the observed H\alpha\ luminosity is a dominant SFR tracer when (u-z)_rest ~< 0.9 mag or when 3.6\mu m photometry > 22 mag (Vega) or when M* < 10^9.7 Mo. We do not find any correlation between the [OII]/H\alpha\ and FIR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions. The luminosity-limited HiZELS sample tends to lie above of the so-called `main sequence' for star-forming galaxies, especially at low M*. This work suggests that obscured star formation is linked to the assembly of M*, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

We have developed a near-infrared spectrometer designed to measure the absolute intensity of the Solar 854.2 nm CaII Fraunhofer line, scattered by interplanetary dust, in the Zodiacal light spectrum. Based on the known equivalent line width in the Solar spectrum, this measurement can derive the Zodiacal brightness, testing models of the Zodiacal light based on morphology that are used to determine the extragalactic background light in absolute photometry measurements. The spectrometer is based on a simple high-resolution tipped filter placed in front of a compact camera with wide-field refractive optics to provide the large optical throughput and high sensitivity required for rocket-borne observations. We discuss the instrument requirements for an accurate measurement of the absolute Zodiacal light brightness, the measured laboratory characterization, and the instrument performance in flight.

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

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

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

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

We use a temperature map of the cosmic microwave background (CMB) obtained using the South Pole Telescope at 150 GHz to construct a map of the gravitational convergence to z ~ 1100, revealing the fluctuations in the projected mass density. This map shows individual features that are significant at the ~ 4 sigma level, providing the first image of CMB lensing convergence. We cross-correlate this map with Herschel/SPIRE maps covering 90 square degrees at wavelengths of 500, 350, and 250 microns. We show that these submillimeter-wavelength (submm) maps are strongly correlated with the lensing convergence map, with detection significances in each of the three submm bands ranging from 6.7 to 8.8 sigma. We fit the measurement of the cross power spectrum assuming a simple constant bias model and infer bias factors of b=1.3-1.8, with a statistical uncertainty of 15%, depending on the assumed model for the redshift distribution of the dusty galaxies that are contributing to the Herschel/SPIRE maps.

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

We have constructed an extended halo model (EHM) which relates the total stellar mass and star-formation rate (SFR) to halo mass (M_h). An empirical relation between the distribution functions of total stellar mass of galaxies and host halo mass, tuned to match the spatial density of galaxies over 0<z<2 and the clustering properties at z~0, is extended to include two different scenarios describing the variation of SFR on M_h. We also present new measurements of the redshift evolution of the average SFR for star-forming galaxies of different stellar mass up to z=2, using data from the Herschel Multi-tiered Extragalactic Survey (HerMES) for infrared-bright galaxies. Combining the EHM with the halo accretion histories from numerical simulations, we trace the stellar mass growth and star-formation history in halos spanning a range of masses. We find that: (1) The intensity of the star-forming activity in halos in the probed mass range has steadily decreased from z~2 to 0; (2) At a given epoch, halos in the mass range between a few times 10^{11} M_Sun and a few times 10^{12} M_Sun are the most efficient at hosting star formation; (3) The peak of SFR density shifts to lower mass halos over time; (4) Galaxies that are forming stars most actively at z~2 evolve into quiescent galaxies in today's group environments, strongly supporting previous claims that the most powerful starbursts at z~2 are progenitors of today's elliptical galaxies.

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 study the far-infrared (IR) and sub-millimeter properties of a sample of ultraviolet (UV) selected galaxies at z\sim1.5. Using stacking at 250, 350 and 500 um from Herschel Space Observatory SPIRE imaging of the COSMOS field obtained within the HerMES key program, we derive the mean IR luminosity as a function of both UV luminosity and slope of the UV continuum beta. The IR to UV luminosity ratio is roughly constant over most of the UV luminosity range we explore. We also find that the IR to UV luminosity ratio is correlated with beta. We observe a correlation that underestimates the correlation derived from low-redshift starburst galaxies, but is in good agreement with the correlation derived from local normal star-forming galaxies. Using these results we reconstruct the IR luminosity function of our UV-selected sample. This luminosity function recovers the IR luminosity functions measured from IR selected samples at the faintest luminosities (Lir ~ 10^{11} L_sun), but might underestimate them at the bright-end (Lir > 5.10^{11} L_sun). For galaxies with 10^{11}<Lir/L_sun<10^{13}, the IR luminosity function of a UV selection recovers (given the differences in IR-based estimates) 52-65 to 89-112 per cent of the star-formation rate density derived from an IR selection. The cosmic star-formation rate density derived from this IR luminosity function is 61-76 to 100-133 per cent of the density derived from IR selections at the same epoch. Assuming the latest Herschel results and conservative stacking measurements, we use a toy model to fully reproduce the far IR luminosity function from our UV selection at z\sim 1.5. This suggests that a sample around 4 magnitudes deeper (i.e. reaching u \sim 30 mag) and a large dispersion of the IR to UV luminosity ratio are required.

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.

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs.

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.

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.

Using far-infrared imaging from the "Herschel Lensing Survey", we derive dust properties of spectroscopically-confirmed cluster member galaxies within two massive systems at z~0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30K). Several sub-LIRG (L_IR < 10^11 L_sun) Bullet Cluster members are much warmer (T_dust > 37K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR-luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (> 3sigma), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e. preferentially located in the cluster periphery, although the galaxy hosts tend towards lower stellar masses (M_* < 10^10 M_sun). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover a SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR-luminosity or dust mass than the currently observed normal star-forming population.

We have developed and characterized an imaging instrument to measure the spatial properties of the diffuse near-infrared extragalactic background light in a search for fluctuations from z > 6 galaxies during the epoch of reionization. The instrument is part of the Cosmic Infrared Background Experiment (CIBER), designed to observe the extragalactic background light above the Earth's atmosphere during a suborbital sounding rocket flight. The imaging instrument incorporates a 2x2 degree field of view, to measure fluctuations over the predicted peak of the spatial power spectrum at 10 arcminutes, and 7"x7" pixels, to remove lower redshift galaxies to a depth sufficient to reduce the low-redshift galaxy clustering foreground below instrumental sensitivity. The imaging instrument employs two cameras with \Delta \lambda / \lambda ~0.5 bandpasses centered at 1.1 and 1.6 microns to spectrally discriminate reionization extragalactic background fluctuations from local foreground fluctuations. CIBER operates at wavelengths where the electromagnetic spectrum of the reionization extragalactic background is thought to peak, and complements fluctuations measurements by AKARI and Spitzer at longer wavelengths. We have characterized the instrument in the laboratory, including measurements of the sensitivity, flat-field response, stray light performance, and noise properties. Several modifications were made to the instrument following a first flight in 2009 February. The instrument performed to specifications in subsequent flights in 2010 July and 2012 March, and the scientific data are now being analyzed.

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

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

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

We present far-infrared (FIR) analysis of 68 Brightest Cluster Galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500um), we calculate the obscured star formation rate (SFR). 22(+6.2,-5.3)% of the BCGs are detected in the far-infrared, with SFR= 1-150 M_sun/yr. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Ha emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2x10^11 L_sun), only a small (<0.4 mag) reddening correction is required for SFR(Ha) to agree with SFR_FIR. The relatively low Ha extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate form normal stellar mass loss.

Absolute spectrophotometric measurements of diffuse radiation at 1 \mu m to 2 \mu m are crucial to our understanding of the radiative content of the Universe from nucleosynthesis since the epoch of reionization, the composition and structure of the Zodiacal dust cloud in our solar system, and the diffuse galactic light arising from starlight scattered by interstellar dust. The Low Resolution Spectrometer (LRS) on the rocket-borne Cosmic Infrared Background Experiment (CIBER) is a \lambda / \Delta \lambda \sim 15-30 absolute spectrophotometer designed to make precision measurements of the absolute near-infrared sky brightness between 0.75 \mu m < \lambda < 2.1 \mu m. This paper presents the optical, mechanical and electronic design of the LRS, as well as the ground testing, characterization and calibration measurements undertaken before flight to verify its performance. The LRS is shown to work to specifications, achieving the necessary optical and sensitivity performance. We describe our understanding and control of sources of systematic error for absolute photometry of the near-infrared extragalactic background light.

The Cosmic Infrared Background Experiment (CIBER) is a suite of four instruments designed to study the near infrared (IR) background light from above the Earth's atmosphere. The instrument package comprises two imaging telescopes designed to characterize spatial anisotropy in the extragalactic IR background caused by cosmological structure during the epoch of reionization, a low resolution spectrometer to measure the absolute spectrum of the extragalactic IR background, and a narrow band spectrometer optimized to measure the absolute brightness of the Zodiacal light foreground. In this paper we describe the design and characterization of the CIBER payload. The detailed mechanical, cryogenic, and electrical design of the system are presented, including all system components common to the four instruments. We present the methods and equipment used to characterize the instruments before and after flight, and give a detailed description of CIBER's flight profile and configurations. CIBER is designed to be recoverable and has flown twice, with modifications to the payload having been informed by analysis of the first flight data. All four instruments performed to specifications during the second flight, and the scientific data from this flight are currently being analyzed.

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.

The atomic CII fine-structure line is one of the brightest lines in a typical star-forming galaxy spectrum with a luminosity ~ 0.1% to 1% of the bolometric luminosity. It is potentially a reliable tracer of the dense gas distribution at high redshifts and could provide an additional probe to the era of reionization. By taking into account of the spontaneous, stimulated and collisional emission of the CII line, we calculate the spin temperature and the mean intensity as a function of the redshift. When averaged over a cosmologically large volume, we find that the CII emission from ionized carbon in individual galaxies is larger than the signal generated by carbon in the intergalactic medium (IGM). Assuming that the CII luminosity is proportional to the carbon mass in dark matter halos, we also compute the power spectrum of the CII line intensity at various redshifts. In order to avoid the contamination from CO rotational lines at low redshift when targeting a CII survey at high redshifts, we propose the cross-correlation of CII and 21-cm line emission from high redshifts. To explore the detectability of the CII signal from reionization, we also evaluate the expected errors on the CII power spectrum and CII-21 cm cross power spectrum based on the design of the future milimeter surveys. We note that the CII-21 cm cross power spectrum contains interesting features that captures physics during reionization, including the ionized bubble sizes and the mean ionization fraction, which are challenging to measure from 21-cm data alone. We propose an instrumental concept for the reionization CII experiment targeting the frequency range of $\sim$ 200 to 300 GHz with 1, 3 and 10 meter apertures and a bolometric spectrometer array with 64 independent spectral pixels with about 20,000 bolometers.

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

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

Ultraviolet emission from the first generation of stars in the Universe ionized the intergalactic medium in a process which was completed by z~6; the wavelength of these photons has been redshifted by (1+z) into the near infrared today and can be measured using instruments situated above the Earth's atmosphere. First flying in February 2009, the Cosmic Infrared Background Experiment (CIBER) comprises four instruments housed in a single reusable sounding rocket borne payload. CIBER will measure spatial anisotropies in the extragalactic IR background caused by cosmological structure from the epoch of reionization using two broadband imaging instruments, make a detailed characterization of the spectral shape of the IR background using a low resolution spectrometer, and measure the absolute brightness of the Zodical light foreground with a high resolution spectrometer in each of our six science fields. This paper presents the scientific motivation for CIBER and details of its first two flights, including a review of the published scientific results from the first flight and an outlook for future reionization science with CIBER data.

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 present a catalog of compact sources derived from the QUaD Galactic Plane Survey. The survey covers ~800 square degrees of the inner galaxy (|b|<4 degrees) in Stokes I, Q and U parameters at 100 and 150 GHz, with angular resolution 5 and 3.5 arcminutes. 505 unique sources are identified in I, of which 239 are spatially matched between frequency bands, with 50 (216) detected at 100 (150) GHz alone; 182 sources are identified as ultracompact HII (UCHII) regions. Approximating the distribution of total intensity source fluxes as a power-law, we find a slope of $\gamma_{S,100}=-1.8\pm0.4$ at 100 GHz, and $\gamma_{S,150}=-2.2\pm0.4$ at 150 GHz. Similarly, the power-law index of the source two-point angular correlation function is $\gamma_{\theta,100}=-1.21\pm0.04$ and $\gamma_{\theta,150}=-1.25\pm0.04$. The total intensity spectral index distribution peaks at $\alpha_{I}\sim0.25$, indicating that dust emission is not the only source of radiation produced by these objects between 100 and 150 GHz; free-free radiation is likely significant in the 100 GHz band. Four sources are detected in polarized intensity P, of which three have matching counterparts in I. Three of the polarized sources lie close to the galactic center, Sagittarius A*, Sagittarius B2 and the Galactic Radio Arc, while the fourth is RCW 49, a bright HII region. An extended polarized source, undetected by the source extraction algorithm on account of its $\sim0.5^{\circ}$ size, is identified visually, and is an isolated example of large-scale polarized emission oriented distinctly from the bulk galactic dust polarization.

We present a survey of ~800 square degrees of the galactic plane observed with the QUaD telescope. The primary product of the survey are maps of Stokes I, Q and U parameters at 100 and 150 GHz, with spatial resolution 5 and 3.5 arcminutes respectively. Two regions are covered, spanning approximately 245-295 and 315-5 degrees in galactic longitude l, and -4<b<+4 degrees in galactic latitude b. At 0.02 degree square pixel size, the median sensitivity is 74 and 107 kJy/sr at 100 GHz and 150 GHz respectively in I, and 98 and 120 kJy/sr for Q and U. In total intensity, we find an average spectral index of 2.35+/-0.01 (stat) +/- 0.02 (sys) for |b|<1 degree, indicative of emission components other than thermal dust. A comparison to published dust, synchrotron and free-free models implies an excess of emission in the 100 GHz QUaD band, while better agreement is found at 150 GHz. A smaller excess is observed when comparing QUaD 100 GHz data to WMAP 5-year W band; in this case the excess is likely due to the wider bandwidth of QUaD. Combining the QUaD and WMAP data, a two-component spectral fit to the inner galactic plane (|b|<1 degree) yields mean spectral indices of -0.32 +/- 0.03 and 2.84 +/- 0.03; the former is interpreted as a combination of the spectral indices of synchrotron, free-free and dust, while the second is attributed largely to the thermal dust continuum. In the same galactic latitude range, the polarization data show a high degree of alignment perpendicular to the expected galactic magnetic field direction, and exhibit mean polarization fraction 1.38+/-0.08 (stat) +/-0.1 (sys) % at 100 GHz and 1.70 +/- 0.06 (stat) +/- 0.1 (sys) % at 150GHz. We find agreement in polarization fraction between QUaD 100 GHz and WMAP W band, the latter giving 1.1 +/- 0.4 %.

The surface density of submillimeter (sub-mm) galaxies as a function of flux, usually termed the source number counts, constrains models of the evolution of the density and luminosity of starburst galaxies. At the faint end of the distribution, direct detection and counting of galaxies are not possible. However, gravitational lensing by clusters of galaxies allows detection of sources which would otherwise be too dim to study. We have used the largest catalog of sub-mm-selected sources along the line of sight to galaxy clusters to estimate the faint end of the 850 micron number counts; integrating to S = 0.10 mJy the equivalent flux density at 850 microns is v I_{v} = 0.24 +/- 0.03 nW/m^2/sr. This provides a lower limit to the extragalactic far-infrared background and is consistent with direct estimates of the full intensity from the FIRAS. The results presented here can help to guide strategies for upcoming surveys carried out with single dish sub-mm instruments.

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

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

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

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

Centaurus (Cen) A represents one of the best candidates for an isolated, compact, highly polarized source that is bright at typical cosmic microwave background (CMB) experiment frequencies. We present measurements of the 4 degree by 2 degree region centered on Cen A with QUaD, a CMB polarimeter whose absolute polarization angle is known to 0.5 degrees. Simulations are performed to assess the effect of misestimation of the instrumental parameters on the final measurement and systematic errors due to the field's background structure and temporal variability from Cen A's nuclear region are determined. The total (Q, U) of the inner lobe region is (1.00 +/- 0.07 (stat.) +/- 0.04 (sys.), -1.72 +/- 0.06 +/- 0.05) Jy at 100 GHz and (0.80 +/- 0.06 +/- 0.06, -1.40 +/- 0.07 +/- 0.08) Jy at 150 GHz, leading to polarization angles and total errors of -30.0 +/- 1.1 degrees and -29.1 +/- 1.7 degrees. These measurements will allow the use of Cen A as a polarized calibration source for future millimeter experiments.

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.

Interplanetary dust (IPD) scatters solar radiation which results in the zodiacal light that dominates the celestial diffuse brightness at optical and near-infrared wavelengths. Both asteroid collisions and cometary ejections produce the IPD, but the relative contribution from these two sources is still unknown. The Low Resolution Spectrometer (LRS) onboard the Cosmic Infrared Background Experiment (CIBER) observed the astrophysical sky spectrum between 750 and 2100 nm over a wide range of ecliptic latitude. The resulting zodiacal light spectrum is redder than the solar spectrum, and shows a broad absorption feature, previously unreported, at approximately 900 nm, suggesting the existence of silicates in the IPD material. The spectral shape of the zodiacal light is isotropic at all ecliptic latitudes within the measurement error. The zodiacal light spectrum, including the extended wavelength range to 2500 nm using IRTS data, is qualitatively similar to the reflectance of S-type asteroids. This result can be explained by the proximity of S-type asteroidal dust to Earth's orbit, and the relativily high albedo of asteridal dust compared with cometary dust.

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

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

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

We evaluate the contribution of cosmic microwave background (CMB) polarization spectra to cosmological parameter constraints. We produce cosmological parameters using high-quality CMB polarization data from the ground-based QUaD experiment and demonstrate for the majority of parameters that there is significant improvement on the constraints obtained from satellite CMB polarization data. We split a multi-experiment CMB dataset into temperature and polarization subsets and show that the best-fit confidence regions for the LCDM 6-parameter cosmological model are consistent with each other, and that polarization data reduces the confidence regions on all parameters. We provide the best limits on parameters from QUaD EE/BB polarization data and we find best-fit parameters from the multi-experiment CMB dataset using the optimal pivot scale of k_p=0.013 Mpc-1 to be {omch2, ombh2, H_0, A_s, n_s, tau}= {0.113, 0.0224, 70.6, 2.29 times 10^-9, 0.960, 0.086}.

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

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

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

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

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

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

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

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

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

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

We present an improved analysis of the final dataset from the QUaD experiment. Using an improved technique to remove ground contamination, we double the effective sky area and hence increase the precision of our CMB power spectrum measurements by ~30% versus that previously reported. In addition, we have improved our modeling of the instrument beams and have reduced our absolute calibration uncertainty from 5% to 3.5% in temperature. The robustness of our results is confirmed through extensive jackknife tests and by way of the agreement we find between our two fully independent analysis pipelines. For the standard 6-parameter LCDM model, the addition of QUaD data marginally improves the constraints on a number of cosmological parameters over those obtained from the WMAP experiment alone. The impact of QUaD data is significantly greater for a model extended to include either a running in the scalar spectral index, or a possible tensor component, or both. Adding both the QUaD data and the results from the ACBAR experiment, the uncertainty in the spectral index running is reduced by ~25% compared to WMAP alone, while the upper limit on the tensor-to-scalar ratio is reduced from r < 0.48 to r < 0.33 (95% c.l). This is the strongest limit on tensors to date from the CMB alone. We also use our polarization measurements to place constraints on parity violating interactions to the surface of last scattering, constraining the energy scale of Lorentz violating interactions to < 1.5 x 10^{-43} GeV (68% c.l.). Finally, we place a robust upper limit on the strength of the lensing B-mode signal. Assuming a single flat band power between l = 200 and l = 2000, we constrain the amplitude of B-modes to be < 0.57 micro-K^2 (95% c.l.).

We present measurements of the cosmic microwave background (CMB) radiation temperature anisotropy in the multipole range 2000<ell<3000 from the QUaD telescope's second and third observing seasons. After masking the brightest point sources our results are consistent with the primary LCDM expectation alone. We estimate the contribution of residual (un-masked) radio point sources using a model calibrated to our own bright source observations, and a full simulation of the source finding and masking procedure. Including this contribution slightly improves the chi^2. We also fit a standard SZ template to the bandpowers and see no strong evidence of an SZ contribution, which is as expected for sigma_8 approx 0.8.

We describe the QUaD experiment, a millimeter-wavelength polarimeter designed to observe the Cosmic Microwave Background (CMB) from a site at the South Pole. The experiment comprises a 2.64 m Cassegrain telescope equipped with a cryogenically cooled receiver containing an array of 62 polarization-sensitive bolometers. The focal plane contains pixels at two different frequency bands, 100 GHz and 150 GHz, with angular resolutions of 5 arcmin and 3.5 arcmin, respectively. The high angular resolution allows observation of CMB temperature and polarization anisotropies over a wide range of scales. The instrument commenced operation in early 2005 and collected science data during three successive Austral winter seasons of observation.

The Cosmic Infrared Background ExpeRiment (CIBER) is a rocket-borne absolute photometry imaging and spectroscopy experiment optimized to detect signatures of first-light galaxies present during reionization in the unresolved IR background. CIBER-I consists of a wide-field two-color camera for fluctuation measurements, a low-resolution absolute spectrometer for absolute EBL measurements, and a narrow-band imaging spectrometer to measure and correct scattered emission from the foreground zodiacal cloud. CIBER-I was successfully flown on February 25th, 2009 and has one more planned flight in early 2010. We propose, after several additional flights of CIBER-I, an improved CIBER-II camera consisting of a wide-field 30 cm imager operating in 4 bands between 0.5 and 2.1 microns. It is designed for a high significance detection of unresolved IR background fluctuations at the minimum level necessary for reionization. With a FOV 50 to 2000 times largerthan existing IR instruments on satellites, CIBER-II will carry out the definitive study to establish the surface density of sources responsible for reionization.

We constrain parity-violating interactions to the surface of last scattering using spectra from the QUaD experiment's second and third seasons of observations by searching for a possible systematic rotation of the polarization directions of CMB photons. We measure the rotation angle due to such a possible "cosmological birefringence" to be 0.55 deg. +/- 0.82 deg. (random) +/- 0.5 deg. (systematic) using QUaD's 100 and 150 GHz TB and EB spectra over the multipole range 200 < l < 2000, consistent with null, and constrain Lorentz violating interactions to < 2^-43 GeV (68% confidence limit). This is the best constraint to date on electrodynamic parity violation on cosmological scales.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

We report results from the second and third seasons of observation with the QUaD experiment. Angular power spectra of the Cosmic Microwave Background are derived for both temperature and polarization at both 100 GHz and 150 GHz, and as cross frequency spectra. All spectra are subjected to an extensive set of jackknife tests to probe for possible systematic contamination. For the implemented data cuts and processing technique such contamination is undetectable. We analyze the difference map formed between the 100 and 150 GHz bands and find no evidence of foreground contamination in polarization. The spectra are then combined to form a single set of results which are shown to be consistent with the prevailing LCDM model. The sensitivity of the polarization results is considerably better than that of any previous experiment -- for the first time multiple acoustic peaks are detected in the E-mode power spectrum at high significance.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

In this paper we present a parameter estimation analysis of the polarization and temperature power spectra from the second and third season of observations with the QUaD experiment. QUaD has for the first time detected multiple acoustic peaks in the E-mode polarization spectrum with high significance. Although QUaD-only parameter constraints are not competitive with previous results for the standard 6-parameter LCDM cosmology, they do allow meaningful polarization-only parameter analyses for the first time. In a standard 6-parameter LCDM analysis we find the QUaD TT power spectrum to be in good agreement with previous results. However, the QUaD polarization data shows some tension with LCDM. The origin of this 1 to 2 sigma tension remains unclear, and may point to new physics, residual systematics or simple random chance. We also combine QUaD with the five-year WMAP data set and the SDSS Luminous Red Galaxies 4th data release power spectrum, and extend our analysis to constrain individual isocurvature mode fractions, constraining cold dark matter density, alpha(cdmi)<0.11 (95 % CL), neutrino density, alpha(ndi)<0.26 (95 % CL), and neutrino velocity, alpha(nvi)<0.23 (95 % CL), modes. Our analysis sets a benchmark for future polarization experiments.

QUaD is a bolometric CMB polarimeter sited at the South Pole, operating at frequencies of 100 and 150 GHz. In this paper we report preliminary results from the first season of operation (austral winter 2005). All six CMB power spectra are presented derived as cross spectra between the 100 and 150 GHz maps using 67 days of observation in a low foreground region of approximately 60 square degrees. This data is a small fraction of the data acquired to date. The measured spectra are consistent with the LCDM cosmological model. We perform jackknife tests which indicate that the observed signal has negligible contamination from instrumental systematics. In addition by using a frequency jackknife we find no evidence for foreground contamination.

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

The Sub-millimetre Common User Bolometer Array (SCUBA) at the JCMT employs a chopping and nodding observation technique to remove variations in the atmospheric signal and improve the long term stability of the instrument. In order to understand systematic effects in SCUBA data, we have analysed single-nod time streams from across SCUBA's lifetime, and present an analysis of an ubiquitous optical pick-up signal connected to the pointing of the secondary mirror. This pick-up is usually removed to a high level by subtracting data from two nod positions and is therefore not obviously present in most reduced SCUBA data. However, if the nod cancellation is not perfect this pick-up can swamp astronomical signals. We discuss various methods which have been suggested to account for this type of imperfect cancellation, and also examine the impact of this pick-up on observations with future bolometric cameras at the JCMT.

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

We present estimates of the confusion limit of measurements of the Sunyaev-Zel'dovich effect (SZ) increment near its peak intensity based on Monte Carlo simulations of sub-millimetre observations of galaxy clusters. Specifically, we evaluate the contribution of gravitational lensing of high redshift background sources by a target cluster to its inferred SZ signal. We find that the background confusion limit is 0.6 mJy per beam without any lensing. With gravitational lensing, the confusion limit increases to 0.9 mJy per beam. Removing bright sources in the gravitationally lensed fields decreases the confusion limit to 0.7 mJy per beam, essentially background levels. These limits are found by removing exceptionally poor fits to the expected SZ profile, which are due to the different structure of the SZ and lensing angular profiles. We conclude that experiments designed to measure the SZ increment must have high enough angular resolution to identify background point sources while surveying enough of the sky to determine the SZ profile. Observations of the SZ increment with SCUBA, which has sufficient sensitivity and resolution for the task, should be straight forward.