We present SMA observations at resolutions from 0.35 to 3 arcseconds of a sample of 34 candidate high redshift dusty star forming galaxies (DSFGs). These sources were selected from the HerMES Herschel survey catalogues to have SEDs rising from 250 to 350 to 500$\mu$m, a population termed 500-risers. We detect counterparts to 24 of these sources, with four having two counterparts. We conclude that the remaining ten sources that lack detected counterparts are likely to have three or more associated sources which blend together to produce the observed Herschel source. We examine the role of lensing, which is predicted to dominate the brightest (F500 $>$ 60mJy) half of our sample. We find that while lensing plays a role, at least 35% of the bright sources are likely to be multiple sources rather than the result of lensing. At fainter fluxes we find a blending rate comparable to, or greater than, the predicted 40%. We determine far-IR luminosities and star formation rates for the non-multiple sources in our sample and conclude that, in the absence of strong lensing, our 500-risers are very luminous systems with L$_{FIR} > 10^{13}$L$_{\odot}$ and star formation rates $> 1000$M$_{\odot}$/yr.

We aim to present 70 spectra of 68 new high-redshift type Ia supernovae (SNeIa) measured at ESO's VLT during the final two years of operation (2006-2008) of the Supernova Legacy Survey (SNLS). We use the full five year SNLS VLT spectral set to investigate a possible spectral evolution of SNeIa populations with redshift and study spectral properties as a function of lightcurve fit parameters and the mass of the host-galaxy. Reduction and extraction are based on both IRAF standard tasks and our own reduction pipeline. Redshifts are estimated from host-galaxy lines whenever possible or alternatively from supernova features. We used the spectrophotometric SNIa model SALT2 combined with a set of galaxy templates that model the host-galaxy contamination to assess the type Ia nature of the candidates. We identify 68 new SNeIa with redshift ranging from z=0.207 to z=0.98 (<z>=0.62). Each spectrum is presented individually along with its best-fit SALT2 model. The five year dataset contains 209 spectra corresponding to 192 SNeIa identified at the VLT. We also publish the redshifts of other candidates (host galaxies or other transients) whose spectra were obtained at the same time as the spectra of live SNe Ia. Using the full VLT SNeIa sample, we build composite spectra around maximum light with cuts in color, lightcurve shape parameter ('stretch'), host-galaxy mass and redshift. We find that high-z SNeIa are bluer, brighter and have weaker intermediate mass element absorption lines than their low-z counterparts at a level consistent with what is expected from selection effects. We also find a flux excess in the range [3000-3400] A for SNeIa in low mass host-galaxies or with locally blue U-V colors, and suggest that the UV flux (or local color) may be used in future cosmological studies as a third standardization parameter in addition to stretch and color.

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.

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

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

We present observations of four rapidly rising (t_{rise}~10d) transients with peak luminosities between those of supernovae (SNe) and superluminous SNe (M_{peak}~-20) - one discovered and followed by the Palomar Transient Factory (PTF) and three by the Supernova Legacy Survey (SNLS). The light curves resemble those of SN 2011kl, recently shown to be associated with an ultra-long-duration gamma ray burst (GRB), though no GRB was seen to accompany our SNe. The rapid rise to a luminous peak places these events in a unique part of SN phase space, challenging standard SN emission mechanisms. Spectra of the PTF event formally classify it as a Type II SN due to broad Halpha emission, but an unusual absorption feature, which can be interpreted as either high velocity Halpha (though deeper than in previously known cases) or Si II (as seen in Type Ia SNe), is also observed. We find that existing models of white dwarf detonations, CSM interaction, shock breakout in a wind (or steeper CSM) and magnetar spindown can not readily explain the observations. We consider the possibility that a "Type 1.5 SN" scenario could be the origin of our events. More detailed models for these kinds of transients and more constraining observations of future such events should help better determine their nature.

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

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.

The Galileon model is a modified gravity model that can explain the late-time accelerated expansion of the Universe. In a previous work, we derived experimental constraints on the Galileon model with no explicit coupling to matter and showed that this model agrees with the most recent cosmological data. In the context of braneworld constructions or massive gravity, the Galileon model exhibits a disformal coupling to matter, which we study in this paper. After comparing our constraints on the uncoupled model with recent studies, we extend the analysis framework to the disformally coupled Galileon model and derive the first experimental constraints on that coupling, using precise measurements of cosmological distances and the growth rate of cosmic structures. In the uncoupled case, with updated data, we still observe a low tension between the constraints set by growth data and those from distances. In the disformally coupled Galileon model, we obtain better agreement with data and favour a non-zero disformal coupling to matter at the $2.5\sigma$ level. This gives an interesting hint of the possible braneworld origin of Galileon theory.

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

Galaxies' rest-frame ultraviolet (UV) properties are often used to directly infer the degree to which dust obscuration affects the measurement of star formation rates. While much recent work has focused on calibrating dust attenuation in galaxies selected at rest-frame ultraviolet wavelengths, locally and at high-$z$, here we investigate attenuation in dusty, star-forming galaxies (DSFGs) selected at far-infrared wavelengths. By combining multiwavelength coverage across 0.15--500\,$\mu$m in the COSMOS field, in particular making use of {\it Herschel} imaging, and a rich dataset on local galaxies, we find a empirical variation in the relationship between rest-frame UV slope ($\beta$) and ratio of infrared-to-ultraviolet emission ($L_{\rm IR}/L_{\rm UV}\equiv\,IRX$) as a function of infrared luminosity, or total star formation rate, SFR. Both locally and at high-$z$, galaxies above SFR$\gt$50\,M$_\odot$\,yr$^{-1}$ deviate from the nominal $IRX-\beta$ relation towards bluer colors by a factor proportional to their increasing IR luminosity. We also estimate contamination rates of DSFGs on high-$z$ dropout searches of $\ll1$\%\ at $z\lt4-10$, providing independent verification that contamination from very dusty foreground galaxies is low in LBG searches. Overall, our results are consistent with the physical interpretation that DSFGs, e.g. galaxies with $>50$\,M$_\odot$\,yr$^{-1}$, are dominated at all epochs by short-lived, extreme burst events, producing many young O and B stars that are primarily, yet not entirely, enshrouded in thick dust cocoons. The blue rest-frame UV slopes of DSFGs are inconsistent with the suggestion that most DSFGs at $z\sim2$ exhibit steady-state star formation in secular disks.

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

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

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

We present 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 present cosmological constraints from a joint analysis of type Ia supernova (SN Ia) observations obtained by the SDSS-II and SNLS collaborations. The data set includes several low-redshift samples (z<0.1), all 3 seasons from the SDSS-II (0.05 < z < 0.4), and 3 years from SNLS (0.2 <z < 1) and totals \ntotc spectroscopically confirmed type Ia supernovae with high quality light curves. We have followed the methods and assumptions of the SNLS 3-year data analysis except for the following important improvements: 1) the addition of the full SDSS-II spectroscopically-confirmed SN Ia sample in both the training of the SALT2 light curve model and in the Hubble diagram analysis (\nsdssc SNe), 2) inter-calibration of the SNLS and SDSS surveys and reduced systematic uncertainties in the photometric calibration, performed blindly with respect to the cosmology analysis, and 3) a thorough investigation of systematic errors associated with the SALT2 modeling of SN Ia light-curves. We produce recalibrated SN Ia light-curves and associated distances for the SDSS-II and SNLS samples. The large SDSS-II sample provides an effective, independent, low-z anchor for the Hubble diagram and reduces the systematic error from calibration systematics in the low-z SN sample. For a flat LCDM cosmology we find Omega_m=0.295+-0.034 (stat+sys), a value consistent with the most recent CMB measurement from the Planck and WMAP experiments. Our result is 1.8sigma (stat+sys) different than the previously published result of SNLS 3-year data. The change is due primarily to improvements in the SNLS photometric calibration. When combined with CMB constraints, we measure a constant dark-energy equation of state parameter w=-1.018+-0.057 (stat+sys) for a flat universe. Adding BAO distance measurements gives similar constraints: w=-1.027+-0.055.

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3.

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 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/

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

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.

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

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.

We present spectra and lightcurves of SNLS 06D4eu and SNLS 07D2bv, two hydrogen-free superluminous supernovae discovered by the Supernova Legacy Survey. At z = 1.588, SNLS 06D4eu is the highest redshift superluminous SN with a spectrum, at M_U = -22.7 is one of the most luminous SNe ever observed, and gives a rare glimpse into the restframe ultraviolet where these supernovae put out their peak energy. SNLS 07D2bv does not have a host galaxy redshift, but based on the supernova spectrum, we estimate it to be at z ~ 1.5. Both supernovae have similar observer-frame griz lightcurves, which map to restframe lightcurves in the U-band and UV, rising in ~ 20 restframe days or longer, and declining over a similar timescale. The lightcurves peak in the shortest wavelengths first, consistent with an expanding blackbody starting near 15,000 K and steadily declining in temperature. We compare the spectra to theoretical models, and identify lines of C II, C III, Fe III, and Mg II in the spectrum of SNLS 06D4eu and SCP 06F6, and find that they are consistent with an expanding explosion of only a few solar masses of carbon, oxygen, and other trace metals. Thus the progenitors appear to be related to those suspected for SNe Ic. A high kinetic energy, 10^52 ergs, is also favored. Normal mechanisms of powering core- collapse or thermonuclear supernovae do not seem to work for these supernovae. We consider models powered by 56Ni decay and interaction with circumstellar material, but find that the creation and spin-down of a magnetar with a period of 2ms, magnetic field of 2 x 10^14 Gauss, and a 3 solar mass progenitor provides the best fit to the data.

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.

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 present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

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

The Galileon model is a modified gravity theory that may provide an explanation for the accelerated expansion of the Universe. This model does not suffer from instabilities or ghost problems (normally associated with higher-order derivative theories), restores local General Relativity -- thanks to the Vainshtein screening effect -- and predicts late time acceleration of the expansion. In this paper, we derive a new definition of the Galileon parameters that allows us to avoid having to choose initial conditions for the Galileon field, and then test this model against precise measurements of the cosmological distances and the rate of growth of cosmic structures. We observe a small tension between the constraints set by growth data and those from distances. However, we find that the Galileon model remains consistent with current observations and is still competitive with the \Lambda CDM model, contrary to what was concluded in recent publications.

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.

We report the discovery of a redshift 1.71 supernova in the GOODS North field. The Hubble Space Telescope (HST) ACS spectrum has almost negligible contamination from the host or neighboring galaxies. Although the rest frame sampled range is too blue to include any Si ii line, a principal component analysis allows us to confirm it as a Type Ia supernova with 92% confidence. A recent serendipitous archival HST WFC3 grism spectrum contributed a key element of the confirmation by giving a host-galaxy redshift of 1.713 +/- 0.007. In addition to being the most distant SN Ia with spectroscopic confirmation, this is the most distant Ia with a precision color measurement. We present the ACS WFC and NICMOS 2 photometry and ACS and WFC3 spectroscopy. Our derived supernova distance is in agreement with the prediction of LambdaCDM.

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

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

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

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

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 study the effects of current systematic errors in Type Ia supernova (SN Ia) measurements on dark energy (DE) constraints using current data from the Supernova Legacy Survey (SNLS). We consider how SN systematic errors affect constraints from combined SN Ia, baryon acoustic oscillations (BAO), and cosmic microwave background (CMB) data, given that SNe Ia still provide the strongest constraints on DE but are arguably subject to more significant systematics than the latter two probes. We focus our attention on the temporal evolution of DE described in terms of principal components (PCs) of the equation of state, though we examine a few of the more common, simpler parametrizations as well. We find that the SN Ia systematics degrade the total generalized figure of merit (FoM), which characterizes constraints in multi-dimensional DE parameter space, by a factor of two to three. Nevertheless, overall constraints obtained on more than five PCs are very good even with current data and systematics. We further show that current constraints are robust to allowing for the finite detection significance of the BAO feature in galaxy surveys.

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.

Deeper understanding of the properties of dark energy via SNIa surveys, and to a large extent other methods as well, will require unprecedented photometric precision. Laboratory and solar photometry and radiometry regularly achieve precisions on the order of parts in ten thousand, but photometric calibration for non-solar astronomy presently remains stuck at the percent or greater level. We discuss our project to erase this discrepancy, and our steps toward achieving laboratory-level photometric precision for surveys late this decade. In particular, we show near-field observations of the balloon-borne light source we are presently testing, in addition to previous work with a calibrated laser source presently in low-Earth orbit. Our technique is additionally applicable to microwave astronomy. Observation of gravitational waves in the polarized CMB will similarly require unprecedented polarimetric and radiometric precision, and we briefly discuss our plans for a calibrated microwave source above the atmosphere as well.

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.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate (SNR_Ia) as a function of redshift for the first four years of data from the Canada-France-Hawaii Telescope (CFHT) Supernova Legacy Survey (SNLS). This analysis includes 286 spectroscopically confirmed and more than 400 additional photometrically identified SNe Ia within the redshift range 0.1<z<1.1. The volumetric SNR_Ia evolution is consistent with a rise to z~1.0 that follows a power-law of the form (1+z)^alpha, with alpha=2.11+/-0.28. This evolutionary trend in the SNLS rates is slightly shallower than that of the cosmic star-formation history over the same redshift range. We combine the SNLS rate measurements with those from other surveys that complement the SNLS redshift range, and fit various simple SN Ia delay-time distribution (DTD) models to the combined data. A simple power-law model for the DTD (i.e., proportional to t^-beta) yields values from beta=0.98+/-0.05 to beta=1.15+/-0.08 depending on the parameterization of the cosmic star formation history. A two-component model, where SNR_Ia is dependent on stellar mass (Mstellar) and star formation rate (SFR) as SNR_Ia(z)=AxMstellar(z) + BxSFR(z), yields the coefficients A=1.9+/-0.1 SNe/yr/M_solar and B=3.3+/-0.2 SNe/yr/(M_solar/yr). More general two-component models also fit the data well, but single Gaussian or exponential DTDs provide significantly poorer matches. Finally, we split the SNLS sample into two populations by the light curve width (stretch), and show that the general behavior in the rates of faster-declining SNe Ia (0.8<s<1.0) is similar, within our measurement errors, to that of the slower objects (1.0<s<1.3) out to z~0.8.

We use the wide-field capabilities of the 2dF fibre positioner and the AAOmega spectrograph on the Anglo-Australian Telescope (AAT) to obtain redshifts of galaxies that hosted supernovae during the first three years of the Supernova Legacy Survey (SNLS). With exposure times ranging from 10 to 60 ksec per galaxy, we were able to obtain redshifts for 400 host galaxies in two SNLS fields, thereby substantially increasing the total number of SNLS supernovae with host galaxy redshifts. The median redshift of the galaxies in our sample that hosted photometrically classified Type Ia supernovae (SNe Ia) is 0.77, which is 25% higher than the median redshift of spectroscopically confirmed SNe Ia in the three-year sample of the SNLS. Our results demonstrate that one can use wide-field fibre-fed multi-object spectrographs on 4m telescopes to efficiently obtain redshifts for large numbers of supernova host galaxies over the large areas of sky that will be covered by future high-redshift supernova surveys, such as the Dark Energy Survey.

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

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

We calculate the average stretch-corrected rise-time of type Ia supernovae (SNe Ia) in the Supernova Legacy Survey. We use the aggregate lightcurves of spectroscopic and photometrically identified SNe Ia to fit the rising part of the lightcurve with a simple quadratic model. We obtain a lightcurve shape corrected, i .e. stretch-corrected, fiducial rise-time of 17.02^{+0.18}_{-0.28} (stat) days. The measured rise-time differs from an earlier finding by the SNLS (Conley et al. 2006) due to the use of different SN Ia templates. We compare it to nearby samples using the same methods and find no evolution in the early part of the lightcurve of SNe Ia up to z=1. We search for variations among different populations, particularly subluminous objects, by dividing the sample in stretch. Bright and slow decliners (s>1.0) have consistent stretch-corrected rise-times compared to fainter and faster decliners (0.8<s<1.0); they are shorter by 0.57^{+0.47}_{-0.50} (stat) days. Subluminous SNe Ia (here defined as objects with s<0.8), although less constrained, are also consistent, with a rise-time of 18.03^{+0.81}_{-1.37} (stat) days. We study several systematic biases and find that the use of different fiducial templates may affect the average rise-time but not the intrinsic differences between populations. Based on our results, we estimate that subluminous SNe Ia are powered by 0.05-0.35 solar masses of radioactive nickel synthesized in the explosion. Our conclusions are the same for the single-stretch and two-stretch parameterizations of the lightcurve.

We present a sample of 485 photometrically identified Type Ia supernova candidates mined from the first three years of data of the CFHT SuperNova Legacy Survey (SNLS). The images were submitted to a deferred processing independent of the SNLS real-time detection pipeline. Light curves of all transient events were reconstructed in the g_M, r_M, i_M and z_M filters and submitted to automated sequential cuts in order to identify possible supernovae. Pure noise and long-term variable events were rejected by light curve shape criteria. Type Ia supernova identification relied on event characteristics fitted to their light curves assuming the events to be normal SNe Ia. The light curve fitter SALT2 was used for this purpose, assigning host galaxy photometric redshifts to the tested events. The selected sample of 485 candidates is one magnitude deeper than that allowed by the SNLS spectroscopic identification. The contamination by supernovae of other types is estimated to be 4%. Testing Hubble diagram residuals with this enlarged sample allows us to measure the Malmquist bias due to spectroscopic selections directly. The result is fully consistent with the precise Monte Carlo based estimate used to correct SN Ia distance moduli in the SNLS 3-year cosmological analyses. This paper demonstrates the feasibility of a photometric selection of high redshift supernovae with known host galaxy redshifts, opening interesting prospects for cosmological analyses from future large photometric SN Ia surveys.

While it is generally accepted that Type Ia supernovae are the result of the explosion of a carbon-oxygen White Dwarf accreting mass in a binary system, the details of their genesis still elude us, and the nature of the binary companion is uncertain. Kasen (2010) points out that the presence of a non-degenerate companion in the progenitor system could leave an observable trace: a flux excess in the early rise portion of the lightcurve caused by the ejecta impact with the companion itself. This excess would be observable only under favorable viewing angles, and its intensity depends on the nature of the companion. We searched for the signature of a non-degenerate companion in three years of Supernova Legacy Survey data by generating synthetic lightcurves accounting for the effects of shocking and comparing true and synthetic time series with Kolmogorov-Smirnov tests. Our most constraining result comes from noting that the shocking effect is more prominent in rest-frame B than V band: we rule out a contribution from white dwarf-red giant binary systems to Type Ia supernova explosions greater than 10% at 2 sigma, and than 20% at 3 sigma level.

We present observational constraints on the nature of dark energy using the Supernova Legacy Survey three year sample (SNLS3) of Guy et al. (2010) and Conley et al. (2011). We use the 472 SNe Ia in this sample, accounting for recently discovered correlations between SN Ia luminosity and host galaxy properties, and include the effects of all identified systematic uncertainties directly in the cosmological fits. Combining the SNLS3 data with the full WMAP7 power spectrum, the Sloan Digital Sky Survey luminous red galaxy power spectrum, and a prior on the Hubble constant H0 from SHOES, in a flat universe we find omega_m=0.269+/-0.015 and w=-1.061+0.069-0.068 -- a 6.5% measure of the dark energy equation-of-state parameter w. The statistical and systematic uncertainties are approximately equal, with the systematic uncertainties dominated by the photometric calibration of the SN Ia fluxes -- without these calibration effects, systematics contribute only a ~2% error in w. When relaxing the assumption of flatness, we find omega_m=0.271+/-0.015, omega_k=-0.002+/-0.006, and w=-1.069+0.091-0.092. Parameterizing the time evolution of w as w(a)=w_0+w_a(1-a), gives w_0=-0.905+/-0.196, w_a=-0.984+1.094-1.097 in a flat universe. All of our results are consistent with a flat, w=-1 universe. The size of the SNLS3 sample allows various tests to be performed with the SNe segregated according to their light curve and host galaxy properties. We find that the cosmological constraints derived from these different sub-samples are consistent. There is evidence that the coefficient, beta, relating SN Ia luminosity and color, varies with host parameters at >4sigma significance (in addition to the known SN luminosity--host relation); however this has only a small effect on the cosmological results and is currently a sub-dominant systematic.

We present optical and near-infrared photometry, as well as ground-based optical spectra and Hubble Space Telescope ultraviolet spectra, of the Type Ia supernova (SN) 2001ay. At maximum light the Si II and Mg II lines indicated expansion velocities of 14,000 km/sec, while Si III and S II showed velocities of 9,000 km/sec There is also evidence for some unburned carbon at 12,000 km/sec. SN 2001ay exhibited a decline-rate parameter Delta m_15(B) = 0.68 \pm 0.05 mag; this and the B-band photometry at t > +25 d past maximum make it the most slowly declining Type Ia SN yet discovered. Three of four super-Chandrasekhar-mass candidates have decline rates almost as slow as this. After correction for Galactic and host-galaxy extinction, SN 2001ay had M_B = -19.19 and M_V = -19.17 mag at maximum light; thus, it was not overluminous in optical bands. In near-infrared bands it was overluminous only at the 2-sigma level at most. For a rise time of 18 d (explosion to bolometric maximum) the implied Ni-56 yield was (0.58 \pm 0.15)/alpha M_Sun, with alpha = L_max/E_Ni probably in the range 1.0 to 1.2. The Ni-56 yield is comparable to that of many Type Ia supernovae. The "normal" Ni-56 yield and the typical peak optical brightness suggest that the very broad optical light curve is explained by the trapping of the gamma rays in the inner regions.

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

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

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

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

We combine high redshift Type Ia supernovae from the first 3 years of the Supernova Legacy Survey (SNLS) with other supernova (SN) samples, primarily at lower redshifts, to form a high-quality joint sample of 472 SNe (123 low-$z$, 93 SDSS, 242 SNLS, and 14 {\it Hubble Space Telescope}). SN data alone require cosmic acceleration at >99.9% confidence, including systematic effects. For the dark energy equation of state parameter (assumed constant out to at least $z=1.4$) in a flat universe, we find $w = -0.91^{+0.16}_{-0.20}(\mathrm{stat}) ^{+0.07}_{-0.14} (\mathrm{sys})$ from SNe only, consistent with a cosmological constant. Our fits include a correction for the recently discovered relationship between host-galaxy mass and SN absolute brightness. We pay particular attention to systematic uncertainties, characterizing them using a systematics covariance matrix that incorporates the redshift dependence of these effects, as well as the shape-luminosity and color-luminosity relationships. Unlike previous work, we include the effects of systematic terms on the empirical light-curve models. The total systematic uncertainty is dominated by calibration terms. We describe how the systematic uncertainties can be reduced with soon to be available improved nearby and intermediate-redshift samples, particularly those calibrated onto USNO/SDSS-like systems.

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.

The rate evolution of subluminous Type Ia Supernovae is presented using data from the Supernova Legacy Survey. This sub-sample represents the faint and rapidly-declining light-curves of the observed supernova Ia (SN Ia) population here defined by low stretch values (s<0.8). Up to redshift z=0.6, we find 18 photometrically-identified subluminous SNe Ia, of which six have spectroscopic redshift (and three are spectroscopically-confirmed SNe Ia). The evolution of the subluminous volumetric rate is constant or slightly decreasing with redshift, in contrast to the increasing SN Ia rate found for the normal stretch population, although a rising behaviour is not conclusively ruled out. The subluminous sample is mainly found in early-type galaxies with little or no star formation, so that the rate evolution is consistent with a galactic mass dependent behavior: $r(z)=A\times M_g$, with $A=(1.1\pm0.3)\times10^{-14}$ SNe per year and solar mass.

We report results from the Supernova Photometric Classification Challenge (SNPCC), a publicly released mix of simulated supernovae (SNe), with types (Ia, Ibc, and II) selected in proportion to their expected rate. The simulation was realized in the griz filters of the Dark Energy Survey (DES) with realistic observing conditions (sky noise, point-spread function and atmospheric transparency) based on years of recorded conditions at the DES site. Simulations of non-Ia type SNe are based on spectroscopically confirmed light curves that include unpublished non-Ia samples donated from the Carnegie Supernova Project (CSP), the Supernova Legacy Survey (SNLS), and the Sloan Digital Sky Survey-II (SDSS-II). A spectroscopically confirmed subset was provided for training. We challenged scientists to run their classification algorithms and report a type and photo-z for each SN. Participants from 10 groups contributed 13 entries for the sample that included a host-galaxy photo-z for each SN, and 9 entries for the sample that had no redshift information. Several different classification strategies resulted in similar performance, and for all entries the performance was significantly better for the training subset than for the unconfirmed sample. For the spectroscopically unconfirmed subset, the entry with the highest average figure of merit for classifying SNe~Ia has an efficiency of 0.96 and an SN~Ia purity of 0.79. As a public resource for the future development of photometric SN classification and photo-z estimators, we have released updated simulations with improvements based on our experience from the SNPCC, added samples corresponding to the Large Synoptic Survey Telescope (LSST) and the SDSS, and provided the answer keys so that developers can evaluate their own analysis.

We present photometric properties and distance measurements of 252 high redshift Type Ia supernovae (0.15 < z < 1.1) discovered during the first three years of the Supernova Legacy Survey (SNLS). These events were detected and their multi-colour light curves measured using the MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), by repeatedly imaging four one-square degree fields in four bands. Follow-up spectroscopy was performed at the VLT, Gemini and Keck telescopes to confirm the nature of the supernovae and to measure their redshifts. Systematic uncertainties arising from light curve modeling are studied, making use of two techniques to derive the peak magnitude, shape and colour of the supernovae, and taking advantage of a precise calibration of the SNLS fields. A flat LambdaCDM cosmological fit to 231 SNLS high redshift Type Ia supernovae alone gives Omega_M = 0.211 +/- 0.034(stat) +/- 0.069(sys). The dominant systematic uncertainty comes from uncertainties in the photometric calibration. Systematic uncertainties from light curve fitters come next with a total contribution of +/- 0.026 on Omega_M. No clear evidence is found for a possible evolution of the slope (beta) of the colour-luminosity relation with redshift.

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.

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

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.

GMOS optical long-slit spectroscopy at the Gemini-North telescope was used to classify targets from the Supernova Legacy Survey (SNLS) from July 2005 and May 2006 - May 2008. During this time, 95 objects were observed. Where possible the objects' redshifts (z) were measured from narrow emission or absorption features in the host galaxy spectrum, otherwise they were measured from the broader supernova features. We present spectra of 68 confirmed or probable SNe Ia from SNLS with redshifts in the range 0.17 \leq z \leq 1.02. In combination with earlier SNLS Gemini and VLT spectra, we used these new observations to measure pseudo-equivalent widths (EWs) of three spectral features - CaII H&K, SiII and MgII - in 144 objects and compared them to the EWs of low-redshift SNe Ia from a sample drawn from the literature. No signs of changes with z are seen for the CaII H&K and MgII features. Systematically lower EW SiII is seen at high redshift, but this can be explained by a change in demographics of the SNe Ia population within a two-component model combined with an observed correlation between EW SiII and photometric lightcurve stretch.

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

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

The Supernova Legacy Survey (SNLS) has produced a high-quality, homogeneous sample of Type Ia supernovae (SNe Ia) out to redshifts greater than z=1. In its first four years of full operation (to June 2007), the SNLS discovered more than 3000 transient candidates, 373 of which have been confirmed spectroscopically as SNe Ia. Use of these SNe Ia in precision cosmology critically depends on an analysis of the observational biases incurred in the SNLS survey due to the incomplete sampling of the underlying SN Ia population. This paper describes our real-time supernova detection and analysis procedures, and uses detailed Monte Carlo simulations to examine the effects of Malmquist bias and spectroscopic sampling. Such sampling effects are found to become apparent at z~0.6, with a significant shift in the average magnitude of the spectroscopically confirmed SN Ia sample towards brighter values for z>0.75. We describe our approach to correct for these selection biases in our three-year SNLS cosmological analysis (SNLS3), and present a breakdown of the systematic uncertainties involved.

We present a method using the SALT2 light curve fitter to determine the redshift of Type Ia supernovae in the Supernova Legacy Survey (SNLS) based on their photometry in g', r', i' and z'. On 289 supernovae of the first three years of SNLS data, we obtain a precision $\sigma_{\Delta z/(1+z)} = 0.022$ on average up to a redshift of 1.0, with a higher precision of 0.016 for z<0.45 and a lower one of 0.025 for z>0.45. The rate of events with $|\Delta z|/(1+z)>0.15$ (catastrophic errors) is 1.4%. Both the precision and the rate of catastrophic errors are better than what can be currently obtained using host galaxy photometric redshifts. Photometric redshifts of this precision may be useful for future experiments which aim to discover up to millions of supernovae Ia but without spectroscopy for most of them.

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

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.

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 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 have publicly released a blinded mix of simulated SNe, with types (Ia, Ib, Ic, II) selected in proportion to their expected rate. The simulation is realized in the griz filters of the Dark Energy Survey (DES) with realistic observing conditions (sky noise, point spread function and atmospheric transparency) based on years of recorded conditions at the DES site. Simulations of non-Ia type SNe are based on spectroscopically confirmed light curves that include unpublished non-Ia samples donated from the Carnegie Supernova Project (CSP), the Supernova Legacy Survey (SNLS), and the Sloan Digital Sky Survey-II (SDSS-II). We challenge scientists to run their classification algorithms and report a type for each SN. A spectroscopically confirmed subset is provided for training. The goals of this challenge are to (1) learn the relative strengths and weaknesses of the different classification algorithms, (2) use the results to improve classification algorithms, and (3) understand what spectroscopically confirmed sub-sets are needed to properly train these algorithms. The challenge is available at www.hep.anl.gov/SNchallenge, and the due date for classifications is May 1, 2010.

We report on work to increase the number of well-measured Type Ia supernovae (SNe Ia) at high redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe Ia that were discovered during 2001 are presented. Additionally, for the two SNe with z>1, we present ground-based J-band photometry from Gemini and the VLT. These are among the most distant SNe Ia for which ground based near-IR observations have been obtained. We add these six SNe Ia together with other data sets that have recently become available in the literature to the Union compilation (Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the most important ones being the refitting of all light curves with the SALT2 fitter and an improved handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 compilation. The flat concordance LambdaCDM model remains an excellent fit to the Union2 data with the best fit constant equation of state parameter w=-0.997^{+0.050}_{-0.054} (stat) ^{+0.077}_{-0.082} (stat+sys\ together) for a flat universe, or w=-1.035^{+0.055}_{-0.059} (stat)^{+0.093}_{-0.097} (stat+sys together) with curvature. We also present improved constraints on w(z). While no significant change in w with redshift is detected, there is still considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data.

(Abridged) Precision cosmology with Type Ia supernovae (SNe Ia) makes use of the fact that SN Ia luminosities depend on their light-curve shapes and colours. Using Supernova Legacy Survey (SNLS) and other data, we show that there is an additional dependence on the global characteristics of their host galaxies: events of the same light-curve shape and colour are, on average, 0.08mag (~4.0sigma) brighter in massive host galaxies (presumably metal-rich) and galaxies with low specific star-formation rates (sSFR). SNe Ia in galaxies with a low sSFR also have a smaller slope ("beta") between their luminosities and colours with ~2.7sigma significance, and a smaller scatter on SN Ia Hubble diagrams (at 95% confidence), though the significance of these effects is dependent on the reddest SNe. SN Ia colours are similar between low-mass and high-mass hosts, leading us to interpret their luminosity differences as an intrinsic property of the SNe and not of some external factor such as dust. If the host stellar mass is interpreted as a metallicity indicator, the luminosity trends are in qualitative agreement with theoretical predictions. We show that the average stellar mass, and therefore the average metallicity, of our SN Ia host galaxies decreases with redshift. The SN Ia luminosity differences consequently introduce a systematic error in cosmological analyses, comparable to the current statistical uncertainties on parameters such as w. We show that the use of two SN Ia absolute magnitudes, one for events in high-mass (metal-rich) galaxies, and one for events in low-mass (metal-poor) galaxies, adequately corrects for the differences. Cosmological fits incorporating these terms give a significant reduction in chi^2 (3.8-4.5sigma). We conclude that future SN Ia cosmological analyses should use a correction of this (or similar) form to control demographic shifts in the galaxy population.

This paper exploits the gravitational magnification of SNe Ia to measure properties of dark matter haloes. The magnification of individual SNe Ia can be computed using observed properties of foreground galaxies and dark matter halo models. We model the dark matter haloes of the galaxies as truncated singular isothermal spheres with velocity dispersion and truncation radius obeying luminosity dependent scaling laws. A homogeneously selected sample of 175 SNe Ia from the first 3-years of the Supernova Legacy Survey (SNLS) in the redshift range 0.2 < z < 1 is used to constrain models of the dark matter haloes associated with foreground galaxies. The best-fitting velocity dispersion scaling law agrees well with galaxy-galaxy lensing measurements. We further find that the normalisation of the velocity dispersion of passive and star forming galaxies are consistent with empirical Faber-Jackson and Tully-Fisher relations, respectively. If we make no assumption on the normalisation of these relations, we find that the data prefer gravitational lensing at the 92 per cent confidence level. Using recent models of dust extinction we deduce that the impact of this effect on our results is very small. We also investigate the brightness scatter of SNe Ia due to gravitational lensing. The gravitational lensing scatter is approximately proportional to the SN Ia redshift. We find the constant of proportionality to be B = 0.055 +0.039 -0.041 mag (B < 0.12 mag at the 95 per cent confidence level). If this model is correct, the contribution from lensing to the intrinsic brightness scatter of SNe Ia is small for the SNLS sample.

The observed brightness of Type Ia supernovae is affected by gravitational lensing caused by the mass distribution along the line of sight, which introduces an additional dispersion into the Hubble diagram. We look for evidence of lensing in the SuperNova Legacy Survey 3-year data set. We investigate the correlation between the residuals from the Hubble diagram and the gravitational magnification based on a modeling of the mass distribution of foreground galaxies. A deep photometric catalog, photometric redshifts, and well established mass luminosity relations are used. We find evidence of a lensing signal with a 2.3 sigma significance. The current result is limited by the number of SNe, their redshift distribution, and the other sources of scatter in the Hubble diagram. Separating the galaxy population into a red and a blue sample has a positive impact on the significance of the signal detection. On the other hand, increasing the depth of the galaxy catalog, the precision of photometric redshifts or reducing the scatter in the mass luminosity relations have little effect. We show that for the full SuperNova Legacy Survey sample (~400 spectroscopically confirmed Type Ia SNe and ~200 photometrically identified Type Ia SNe), there is an 80% probability of detecting the lensing signal with a 3 sigma significance.

We have combined the large SN Ia database of the Canada-France-Hawaii Telescope Supernova Legacy Survey and catalogs of galaxies with photometric redshifts, VLA 1.4 GHz radio sources, and Spitzer infrared sources. We present eight SNe Ia in early-type host galaxies which have counterparts in the radio and infrared source catalogs. We find the SN Ia rate in subsets of radio and infrared early-type galaxies is ~1-5 times the rate in all early-type galaxies, and that any enhancement is always <~ 2 sigma. Rates in these subsets are consistent with predictions of the two component "A+B" SN Ia rate model. Since infrared properties of radio SN Ia hosts indicate dust obscured star formation, we incorporate infrared star formation rates into the "A+B" model. We also show the properties of SNe Ia in radio and infrared galaxies suggest the hosts contain dust and support a continuum of delay time distributions for SNe Ia, although other delay time distributions cannot be ruled out based on our data.

We use multi-wavelength, matched aperture, integrated photometry from GALEX, SDSS and the RC3 to estimate the physical properties of 166 nearby galaxies hosting 168 well-observed Type Ia supernovae (SNe Ia). Our data corroborate well-known features that have been seen in other SN Ia samples. Specifically, hosts with active star formation produce brighter and slower SNe Ia on average, and hosts with luminosity-weighted ages older than 1 Gyr produce on average more faint, fast and fewer bright, slow SNe Ia than younger hosts. New results include that in our sample, the faintest and fastest SNe Ia occur only in galaxies exceeding a stellar mass threshhold of ~10^10 M_sun, indicating that their progenitors must arise in populations that are older and/or more metal rich than the general SN Ia population. A low host extinction sub-sample hints at a residual trend in peak luminosity with host age, after correcting for light-curve shape, giving the appearance that older hosts produce less-extincted SNe Ia on average. This has implications for cosmological fitting of SNe Ia and suggests that host age could be useful as a parameter in the fitting. Converting host mass to metallicity and computing 56Ni mass from the supernova light curves, we find that our local sample is consistent with a model that predicts a shallow trend between stellar metallicity and the 56Ni mass that powers the explosion, but we cannot rule out the absence of a trend. We measure a correlation between 56Ni mass and host age in the local universe that is shallower and not as significant as that seen at higher redshifts. The details of the age -- 56Ni mass correlations at low and higher redshift imply a luminosity-weighted age threshhold of ~3 Gyr for SN Ia hosts, above which they are less likely to produce SNe Ia with 56Ni masses above ~0.5 M_sun. (Abridged)

We present 139 spectra of 124 Type Ia supernovae (SNeIa) that were observed at the ESO/VLT during the first three years of the Canada-France-Hawai Telescope (CFHT) Supernova Legacy Survey (SNLS). This homogeneous data set is used to test for redshift evolution of SNeIa spectra, and will be used in the SNLS 3rd year cosmological analyses. Spectra have been reduced and extracted with a dedicated pipeline that uses photometric information from deep CFHT Legacy Survey (CFHT-LS) reference images to trace, at sub-pixel accuracy, the position of the supernova on the spectrogram as a function of wavelength. It also separates the supernova and its host light in 60% of cases. The identification of the supernova candidates is performed using a spectrophotometric SNIa model. A total of 124 SNeIa, roughly 50% of the overall SNLS spectroscopic sample, have been identified using the ESO/VLT during the first three years of the survey. Their redshifts range from z=0.149 to z=1.031. The average redshift of the sample is z=0.63+/-0.02. This constitutes the largest SNIa spectral set to date in this redshift range. The spectra are presented along with their best-fit spectral SNIa model and a host model where relevant. In the latter case, a host subtracted spectrum is also presented. We produce average spectra for pre-maximum, maximum and post-maximum epochs for both z<0.5 and z>=0.5 SNeIa. We find that z<0.5 spectra have deeper intermediate mass element absorptions than z>= 0.5 spectra. The differences with redshift are consistent with the selection of brighter and bluer supernovae at higher redshift.

We present the photometric calibration of the Supernova Legacy Survey (SNLS) fields. The SNLS aims at measuring the distances to SNe Ia at (0.3<z<1) using MegaCam, the 1 deg^2 imager on the Canada-France-Hawaii Telescope (CFHT). The uncertainty affecting the photometric calibration of the survey dominates the systematic uncertainty of the key measurement of the survey, namely the dark energy equation of state. The photometric calibration of the SNLS requires obtaining a uniform response across the imager, calibrating the science field stars in each survey band (SDSS-like ugriz bands) with respect to standards with known flux in the same bands, and binding the calibration to the UBVRI Landolt standards used to calibrate the nearby SNe from the literature necessary to produce cosmological constraints. The spatial non-uniformities of the imager photometric response are mapped using dithered observations of dense stellar fields. Photometric zero-points against Landolt standards are obtained. The linearity of the instrument is studied. We show that the imager filters and photometric response are not uniform and publish correction maps. We present models of the effective passbands of the instrument as a function of the position on the focal plane. We define a natural magnitude system for MegaCam. We show that the systematics affecting the magnitude-to-flux relations can be reduced if we use the spectrophotometric standard star BD +17 4708 instead of Vega as a fundamental flux standard. We publish ugriz catalogs of tertiary standards for all the SNLS fields.

The Carnegie Supernova Project (CSP) is designed to measure the luminosity distance for Type Ia supernovae (SNe Ia) as a function of redshift, and to set observational constraints on the dark energy contribution to the total energy content of the Universe. The CSP differs from other projects to date in its goal of providing an I-band {rest-frame} Hubble diagram. Here we present the first results from near-infrared (NIR) observations obtained using the Magellan Baade telescope for SNe Ia with 0.1 < z < 0.7. We combine these results with those from the low-redshift CSP at z <0.1 (Folatelli et al. 2009). We present light curves and an I-band Hubble diagram for this first sample of 35 SNe Ia and we compare these data to 21 new SNe Ia at low redshift. These data support the conclusion that the expansion of the Universe is accelerating. When combined with independent results from baryon acoustic oscillations (Eisenstein et al. 2005), these data yield Omega_m = 0.27 +/- 0.0 (statistical), and Omega_DE = 0.76 +/- 0.13 (statistical) +/- 0.09 (systematic), for the matter and dark energy densities, respectively. If we parameterize the data in terms of an equation of state, w, assume a flat geometry, and combine with baryon acoustic oscillations, we find that w = -1.05 +/- 0.13 (statistical) +/- 0.09 (systematic). The largest source of systematic uncertainty on w arises from uncertainties in the photometric calibration, signaling the importance of securing more accurate photometric calibrations for future supernova cosmology programs. Finally, we conclude that either the dust affecting the luminosities of SNe Ia has a different extinction law (R_V = 1.8) than that in the Milky Way (where R_V = 3.1), or that there is an additional intrinsic color term with luminosity for SNe Ia independent of the decline rate.

We use three years of data from the Supernova Legacy Survey (SNLS) to study the general properties of core-collapse and type Ia supernovae. This is the first such study using the "rolling search" technique which guarantees well-sampled SNLS light curves and good efficiency for supernovae brighter than $i^\prime\sim24$. Using host photometric redshifts, we measure the supernova absolute magnitude distribution down to luminosities $4.5 {\rm mag}$ fainter than normal SNIa. Using spectroscopy and light-curve fitting to discriminate against SNIa, we find a sample of 117 core-collapse supernova candidates with redshifts $z<0.4$ (median redshift of 0.29) and measure their rate to be larger than the type Ia supernova rate by a factor $4.5\pm0.8(stat.) \pm0.6 (sys.)$. This corresponds to a core-collapse rate at $z=0.3$ of $[1.42\pm 0.3(stat.) \pm0.3(sys.)]\times10^{-4}\yr^{-1}(h_{70}^{-1}\Mpc)^{-3}$.

The answers to fundamental science questions in astrophysics, ranging from the history of the expansion of the universe to the sizes of nearby stars, hinge on our ability to make precise measurements of diverse astronomical objects. As our knowledge of the underlying physics of objects improves along with advances in detectors and instrumentation, the limits on our capability to extract science from measurements is set, not by our lack of understanding of the nature of these objects, but rather by the most mundane of all issues: the precision with which we can calibrate observations in physical units. We stress the need for a program to improve upon and expand the current networks of spectrophotometrically calibrated stars to provide precise calibration with an accuracy of equal to and better than 1% in the ultraviolet, visible and near-infrared portions of the spectrum, with excellent sky coverage and large dynamic range.

In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0<z<4. Aside from cosmology, SNe Ia regulate galactic and cluster chemical evolution, inform stellar evolution, and are laboratories for extreme physics. Essential probes of SNe Ia in these contexts include spectroscopy from the UV to the IR, X-ray cluster and SN remnant observations, spectropolarimetry, and advanced theoretical studies. While there are an abundance of discovery facilities planned, there is a deficit of follow-up resources. Living in the systematics era demands deep understanding rather than larger statistics. NOAO ReSTAR initiative to build 2-4m telescopes would provide necessary follow-up capability. Finally, to fully exploit LSST, well-matched wide-field spectroscopic capabilities are desirable.

Timmes, Brown & Truran found that metallicity variations could theoretically account for a 25% variation in the mass of 56Ni synthesized in Type Ia supernovae (SNe Ia), and thus account for a large fraction of the scatter in observed SN Ia luminosities. Higher-metallicity progenitors are more neutron-rich, producing more stable burning products relative to radioactive 56Ni. We develop a new method for estimating bolometric luminosity and 56Ni yield in SNe Ia and use it to test the theory with data from the Supernova Legacy Survey. We find that the average 56Ni yield does drop in SNe Ia from high metallicity environments, but the theory can only account for 7%--10% of the dispersion in SN Ia 56Ni mass, and thus luminosity. This is because the effect is dominant at metallicities significantly above solar, whereas we find that SN hosts have predominantly subsolar or only moderately above-solar metallicities. We also show that allowing for changes in O/Fe with the metallicity [Fe/H] does not have a major effect on the theoretical prediction of Timmes, Brown & Truran, so long as one is using the O/H as the independent variable. Age may have a greater effect than metallicity -- we find that the luminosity weighted age of the host galaxy is correlated with 56Ni yield, and thus more massive progenitors give rise to more luminous explosions. This is hard to understand if most SNe Ia explode when the primaries reach the Chandrasekhar mass. Finally, we test the findings of Gallagher et al., that the residuals of SNe Ia from the Hubble diagram are correlated with host galaxy metallicity, and we find no such correlation.

Massive stars undergo a violent death when the supply of nuclear fuel in their cores is exhausted, resulting in a catastrophic "core-collapse" supernova. Such events are usually only detected at least a few days after the star has exploded. Observations of the supernova SNLS-04D2dc with the Galaxy Evolution Explorer space telescope reveal a radiative precursor from the supernova shock before the shock reached the surface of the star and show the initial expansion of the star at the beginning of the explosion. Theoretical models of the ultraviolet light curve confirm that the progenitor was a red supergiant, as expected for this type of supernova. These observations provide a way to probe the physics of core-collapse supernovae and the internal structures of their progenitor stars

We present SiFTO, a new empirical method for modeling type Ia supernovae (SNe Ia) light curves by manipulating a spectral template. We make use of high-redshift SN observations when training the model, allowing us to extend it bluer than rest frame U. This increases the utility of our high-redshift SN observations by allowing us to use more of the available data. We find that when the shape of the light curve is described using a stretch prescription, applying the same stretch at all wavelengths is not an adequate description. SiFTO therefore uses a generalization of stretch which applies different stretch factors as a function of both the wavelength of the observed filter and the stretch in the rest-frame B band. We compare SiFTO to other published light-curve models by applying them to the same set of SN photometry, and demonstrate that SiFTO and SALT2 perform better than the alternatives when judged by the scatter around the best fit luminosity distance relationship. We further demonstrate that when SiFTO and SALT2 are trained on the same data set the cosmological results agree.

For the first time the cross-correlation between type Ia supernova host galaxies and surrounding field galaxies is measured using the Supernova Legacy Survey sample. Over the z=0.2 to 0.9 redshift range we find that supernova hosts are correlated an average of 60% more strongly than similarly selected field galaxies over the 3-100 arcsec range and about a factor of 3 more strongly below 10 arcsec. The correlation errors are empirically established with a jackknife analysis of the four SNLS fields. The hosts are more correlated than the field at a significance of 99% in the fitted amplitude and slope, with the point-by-point difference of the two correlation functions having a reduced $\chi^2$ for 8 degrees of freedom of 4.3, which has a probability of random occurrence of less than 3x10^{-5}. The correlation angle is 1.5+/-0.5 arcsec, which deprojects to a fixed co-moving correlation length of approximately 6.5+/- 2/h mpc. Weighting the field galaxies with the mass and star formation rate supernova frequencies of the simple A+B model produces good agreement with the observed clustering. We conclude that these supernova clustering differences are primarily the expected outcome of the dependence of supernova rates on galaxy masses and stellar populations with their clustering environment.

We present a new compilation of Type Ia supernovae (SNe Ia), a new dataset of low-redshift nearby-Hubble-flow SNe and new analysis procedures to work with these heterogeneous compilations. This ``Union'' compilation of 414 SN Ia, which reduces to 307 SNe after selection cuts, includes the recent large samples of SNe Ia from the Supernova Legacy Survey and ESSENCE Survey, the older datasets, as well as the recently extended dataset of distant supernovae observed with HST. A single, consistent and blind analysis procedure is used for all the various SN Ia subsamples, and a new procedure is implemented that consistently weights the heterogeneous data sets and rejects outliers. We present the latest results from this Union compilation and discuss the cosmological constraints from this new compilation and its combination with other cosmological measurements (CMB and BAO). The constraint we obtain from supernovae on the dark energy density is $\Omega_\Lambda= 0.713^{+0.027}_{-0.029} (stat)}^{+0.036}_{-0.039} (sys)}$, for a flat, LCDM Universe. Assuming a constant equation of state parameter, $w$, the combined constraints from SNe, BAO and CMB give $w=-0.969^{+0.059}_{-0.063}(stat)^{+0.063}_{-0.066} (sys)$. While our results are consistent with a cosmological constant, we obtain only relatively weak constraints on a $w$ that varies with redshift. In particular, the current SN data do not yet significantly constrain $w$ at $z>1$. With the addition of our new nearby Hubble-flow SNe Ia, these resulting cosmological constraints are currently the tightest available.

The Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS) has created a large homogeneous database of intermediate redshift (0.2 < z < 1.0) type Ia supernovae (SNe Ia). The SNLS team has shown that correlations exist between SN Ia rates, properties, and host galaxy star formation rates. The SNLS SN Ia database has now been combined with a photometric redshift galaxy catalog and an optical galaxy cluster catalog to investigate the possible influence of galaxy clustering on the SN Ia rate, over and above the expected effect due to the dependence of SFR on clustering through the morphology-density relation. We identify three cluster SNe Ia, plus three additional possible cluster SNe Ia, and find the SN Ia rate per unit mass in clusters at intermediate redshifts is consistent with the rate per unit mass in field early-type galaxies and the SN Ia cluster rate from low redshift cluster targeted surveys. We also find the number of SNe Ia in cluster environments to be within a factor of two of expectations from the two component SNIa rate model.

We analyze the mean rest-frame ultraviolet (UV) spectrum of Type Ia Supernovae (SNe Ia) and its dispersion using high signal-to-noise Keck-I/LRIS-B spectroscopy for a sample of 36 events at intermediate redshift (z=0.5) discovered by the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). We introduce a new method for removing host galaxy contamination in our spectra, exploiting the comprehensive photometric coverage of the SNLS SNe and their host galaxies, thereby providing the first quantitative view of the UV spectral properties of a large sample of distant SNe Ia. Although the mean SN Ia spectrum has not evolved significantly over the past 40% of cosmic history, precise evolutionary constraints are limited by the absence of a comparable sample of high quality local spectra. Within the high-redshift sample, we discover significant UV spectral variations and exclude dust extinction as the primary cause by examining trends with the optical SN color. Although progenitor metallicity may drive some of these trends, the variations we see are much larger than predicted in recent models and do not follow expected patterns. An interesting new result is a variation seen in the wavelength of selected UV features with phase. We also demonstrate systematic differences in the SN Ia spectral features with SN light curve width in both the UV and the optical. We show that these intrinsic variations could represent a statistical limitation in the future use of high-redshift SNe Ia for precision cosmology. We conclude that further detailed studies are needed, both locally and at moderate redshift where the rest-frame UV can be studied precisely, in order that future missions can confidently be planned to fully exploit SNe Ia as cosmological probes [ABRIDGED].

Aims: We present a quantitative study of a new data set of high redshift Type Ia supernovae spectra, observed at the Gemini telescopes during the first 34 months of the Supernova Legacy Survey. During this time 123 supernovae candidates were observed, of which 87 have been identified as SNe Ia at a median redshift of z=0.720. Spectra from the entire second year of the survey and part of the third year (59 total SNe candidates with 46 confirmed SNe Ia) are published here for the first time. The spectroscopic measurements made on this data set are used determine if these distant SNe comprise a population similar to those observed locally. Methods: Rest-frame equivalent width and ejection velocity measurements are made on four spectroscopic features. Corresponding measurements are presented for a set of 167 spectra from 24 low-z SNe Ia from the literature. Results: We show that there exists a sample at high redshift with properties similar to nearby SNe. No significant difference was found between the distributions of measurements at low and high redsift for three of the features. The fourth feature displays a possible difference that should be investigated further. Correlations between Type Ia SNe properties and host galaxy morphology were also found to be similar at low and high z, and within each host galaxy class we see no evidence for redshift-evolution in SN properties. A new correlation between SNe Ia peak magnitude and the equivalent width of SiII absorption is presented. We demonstrate that this correlation reduces the scatter in SNe Ia luminosity distances in a manner consistent with the lightcurve shape-luminosity corrections that are used for Type Ia SNe cosmology. Conclusions: We show that this new sample of SNLS SNe Ia has spectroscopic properties similar to nearby objects. (Abridged)

Recent studies indicate that Type Ia supernovae (SNe Ia) consist of two groups - a "prompt" component whose rates are proportional to the host galaxy star formation rate, whose members have broader lightcurves and are intrinsically more luminous, and a "delayed" component whose members take several Gyr to explode, have narrower lightcurves, and are intrinsically fainter. As cosmic star formation density increases with redshift, the prompt component should begin to dominate. We use a two-component model to predict that the average lightcurve width should increase by 6% from z=0-1.5. Using data from various searches we find an 8.1% +/- 2.7% increase in average lightcurve width for non-subluminous SNe Ia from z=0.03 - 1.12, corresponding to an increase in the average intrinsic luminosity of 12%. To test whether there is any bias after supernovae are corrected for lightcurve shape we use published data to mimic the effect of population evolution and find no significant difference in the measured dark energy equation of state parameter, w. However, future measurements of time-variable w will require standardization of SN Ia magnitudes to 2% up to z=1.7, and it is not yet possible to assess whether lightcurve shape correction works at this level of precision. Another concern at z=1.5 is the expected order of magnitude increase in the number of SNe Ia that cannot be calibrated by current methods.

We examine recent evidence from the luminosity-redshift relation of Type Ia Supernovae (SNe Ia) for the $\sim 3 \sigma$ detection of a ``Hubble bubble'' -- a departure of the local value of the Hubble constant from its globally averaged value \citep{Jha:07}. By comparing the MLCS2k2 fits used in that study to the results from other light-curve fitters applied to the same data, we demonstrate that this is related to the interpretation of SN color excesses (after correction for a light-curve shape-color relation) and the presence of a color gradient across the local sample. If the slope of the linear relation ($\beta$) between SN color excess and luminosity is fit empirically, then the bubble disappears. If, on the other hand, the color excess arises purely from Milky Way-like dust, then SN data clearly favors a Hubble bubble. We demonstrate that SN data give $\beta \simeq 2$, instead of the $\beta \simeq 4$ one would expect from purely Milky-Way-like dust. This suggests that either SN intrinsic colors are more complicated than can be described with a single light-curve shape parameter, or that dust around SN is unusual. Disentangling these possibilities is both a challenge and an opportunity for large-survey SN Ia cosmology.

We present a measurement of the distant Type Ia supernova rate derived from the first two years of the Canada -- France -- Hawaii Telescope Supernova Legacy Survey. We observed four one-square degree fields with a typical temporal frequency of <Delta t> ~ 4 observer-frame days over time spans of from 158 to 211 days per season for each field, with breaks during full moon. We used 8-10 meter-class telescopes for spectroscopic followup to confirm our candidates and determine their redshifts. Our starting sample consists of 73 spectroscopically verified Type Ia supernovae in the redshift range 0.2 < z < 0.6. We derive a volumetric SN Ia rate of r_V(<z>=0.47) = 0.42^{+0.13}_{-0.09} (systematic) +- 0.06 (statistical) X 10^-4 yr^-1 Mpc^3, assuming h = 0.7, Omega_m = 0.3 and a flat cosmology. Using recently published galaxy luminosity functions derived in our redshift range, we derive a SN Ia rate per unit luminosity of r_L(<z>=0.47) = 0.154^{+0.048}_{-0.033} (systematic) ^{+0.039}_{-0.031} (statistical) SNu. Using our rate alone, we place an upper limit on the component of SN Ia production that tracks the cosmic star formation history of 1 SN Ia per 10^3 M_sun of stars formed. Our rate and other rates from surveys using spectroscopic sample confirmation display only a modest evolution out to z=0.55.

We quantify the effect of supernova Type Ia peculiar velocities on the derivation of cosmological parameters. The published distant and local Ia SNe used for the Supernova Legacy Survey first-year cosmology report form the sample for this study. While previous work has assumed that the local SNe are at rest in the CMB frame (the No Flow assumption), we test this assumption by applying peculiar velocity corrections to the local SNe using three different flow models. The models are based on the IRAS PSCz galaxy redshift survey, have varying beta = Omega_m^0.6/b, and reproduce the Local Group motion in the CMB frame. These datasets are then fit for w, Omega_m, and Omega_Lambda using flatness or LambdaCDM and a BAO prior. The chi^2 statistic is used to examine the effect of the velocity corrections on the quality of the fits. The most favored model is the beta=0.5 model, which produces a fit significantly better than the No Flow assumption, consistent with previous peculiar velocity studies. By comparing the No Flow assumption with the favored models we derive the largest potential systematic error in w caused by ignoring peculiar velocities to be Delta w = +0.04. For Omega_Lambda, the potential error is Delta Omega_Lambda = -0.04 and for Omega_m, the potential error is Delta Omega_m < +0.01. The favored flow model (beta=0.5) produces the following cosmological parameters: w = -1.08 (+0.09,-0.08), Omega_m = 0.27 (+0.02,-0.02) assuming a flat cosmology, and Omega_Lambda = 0.80 (+0.08,-0.07) and Omega_m = 0.27 (+0.02,-0.02) for a w = -1 (LambdaCDM) cosmology.

We develop a method to measure the strength of the absorption features in Type Ia supernova (SN Ia) spectra and use it to make a quantitative comparison between the spectra of Type Ia supernovae at low and high redshifts. In this case study, we apply the method to 12 high-redshift (0.212 < z < 0.912) SNe Ia observed by the Supernova Cosmology Project . Through measurements of the strengths of these features and of the blueshift of the absorption minimum in Ca II H&K, we show that the spectra of the high-redshift SNe Ia are quantitatively similar to spectra of nearby SNe Ia (z < 0.15). One supernova in our high redshift sample, SN 2002fd at z=0.279, is found to have spectral characteristics that are associated with peculiar SN 1991T/SN 1999aa-like supernovae.

With the advent of large dedicated Type Ia supernova (SN Ia) surveys, K-corrections of SNe Ia and their uncertainties have become especially important in the determination of cosmological parameters. While K-corrections are largely driven by SN Ia broad-band colors, it is shown here that the diversity in spectral features of SNe Ia can also be important. For an individual observation, the statistical errors from the inhomogeneity in spectral features range from 0.01 (where the observed and rest-frame filters are aligned) to 0.04 (where the observed and rest-frame filters are misaligned). To minimize the systematic errors caused by an assumed SN Ia spectral energy distribution (SED), we outline a prescription for deriving a mean spectral template time series which incorporates a large and heterogeneous sample of observed spectra. We then remove the effects of broad-band colors and measure the remaining uncertainties in the K-corrections associated with the diversity in spectral features. Finally, we present a template spectroscopic sequence near maximum light for further improvement on the K-correction estimate. A library of ~600 observed spectra of ~100 SNe Ia from heterogeneous sources is used for the analysis.

We present an empirical model of Type Ia supernovae spectro-photometric evolution with time. The model is built using a large data set including light-curves and spectra of both nearby and distant supernovae, the latter being observed by the SNLS collaboration. We derive the average spectral sequence of Type Ia supernovae and their main variability components including a color variation law. The model allows us to measure distance moduli in the spectral range 2500-8000 A with calculable uncertainties, including those arising from variability of spectral features. Thanks to the use of high-redshift SNe to model the rest-frame UV spectral energy distribution, we are able to derive improved distance estimates for SNe Ia in the redshift range 0.8<z<1.1. The model can also be used to improve spectroscopic identification algorithms, and derive photometric redshifts of distant Type Ia supernovae.

We present a progress report on a project to derive the evolution of the volumetric supernova Type Ia rate from the Supernova Legacy Survey. Our preliminary estimate of the rate evolution divides the sample from Neill et al. (2006) into two redshift bins: 0.2 < z < 0.4, and 0.4 < z < 0.6. We extend this by adding a bin from the sample analyzed in Sullivan et al. (2006) in the range 0.6 < z < 0.75 from the same time period. We compare the derived trend with previously published rates and a supernova Type Ia production model having two components: one component associated closely with star formation and an additional component associated with host galaxy mass. Our observed trend is consistent with this model, which predicts a rising SN Ia rate out to at least z=2.

We present constraints on the dark energy equation-of-state parameter, w=P/(rho c^2), using 60 Type Ia supernovae (SNe Ia) from the ESSENCE supernova survey. We derive a set of constraints on the nature of the dark energy assuming a flat Universe. By including constraints on (Omega_M, w) from baryon acoustic oscillations, we obtain a value for a static equation-of-state parameter w=-1.05^{+0.13}_{-0.12} (stat; 1 sigma) +- 0.11 (sys) and Omega_M=0.274^{+0.033}_{-0.020} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.96. These results are consistent with those reported by the SuperNova Legacy Survey in a similar program measuring supernova distances and redshifts. We evaluate sources of systematic error that afflict supernova observations and present Monte Carlo simulations that explore these effects. Currently, the largest systematic currently with the potential to affect our measurements is the treatment of extinction due to dust in the supernova host galaxies. Combining our set of ESSENCE SNe Ia with the SuperNova Legacy Survey SNe Ia, we obtain a joint constraint of w=-1.07^{+0.09}_{-0.09} (stat; 1 sigma) +- 0.12 (sys), Omega_M=0.267^{+0.028}_{-0.018} (stat; 1 sigma) with a best-fit chi^2/DoF of 0.91. The current SN Ia data are fully consistent with a cosmological constant.

The acceleration of the expansion of the universe, and the need for Dark Energy, were inferred from the observations of Type Ia supernovae (SNe Ia). There is consensus that SNe Ia are thermonuclear explosions that destroy carbon-oxygen white dwarf stars that accrete matter from a companion star, although the nature of this companion remains uncertain. SNe Ia are thought to be reliable distance indicators because they have a standard amount of fuel and a uniform trigger -- they are predicted to explode when the mass of the white dwarf nears the Chandrasekhar mass -- 1.4 solar masses. Here we show that the high redshift supernova SNLS-03D3bb has an exceptionally high luminosity and low kinetic energy that both imply a super-Chandrasekhar mass progenitor. Super-Chandrasekhar mass SNe Ia should preferentially occur in a young stellar population, so this may provide an explanation for the observed trend that overluminous SNe Ia only occur in young environments. Since this supernova does not obey the relations that allow them to be calibrated as standard candles, and since no counterparts have been found at low redshift, future cosmology studies will have to consider contamination from such events.

We compare the rise times of nearby and distant Type Ia supernovae (SNe Ia) as a test for evolution using 73 high-redshift spectroscopically-confirmed SNe Ia from the first two years of the five year Supernova Legacy Survey (SNLS) and published observations of nearby SN. Because of the ``rolling'' search nature of the SNLS, our measurement is approximately 6 times more precise than previous studies, allowing for a more sensitive test of evolution between nearby and distant supernovae. Adopting a simple $t^2$ early-time model (as in previous studies), we find that the rest-frame $B$ rise times for a fiducial SN Ia at high and low redshift are consistent, with values $19.10^{+0.18}_{-0.17}({stat}) \pm 0.2 ({syst})$ and $19.58^{+0.22}_{-0.19}$ days, respectively; the statistical significance of this difference is only 1.4 \sg . The errors represent the uncertainty in the mean rather than any variation between individual SN. We also compare subsets of our high-redshift data set based on decline rate, host galaxy star formation rate, and redshift, finding no substantive evidence for any subsample dependence.

We present measurements of \Omega_m and \Omega_{\Lambda} from a blind analysis of 21 high-redshift supernovae using a new technique (CMAGIC) for fitting the multi-color lightcurves of Type Ia supernovae, first introduced in Wang et al. (2003). CMAGIC takes advantage of the remarkably simple behavior of Type Ia supernovae on color-magnitude diagrams, and has several advantages over current techniques based on maximum magnitudes. Among these are a reduced sensitivity to host galaxy dust extinction, a shallower luminosity-width relation, and the relative simplicity of the fitting procedure. This allows us to provide a cross check of previous supernova cosmology results, despite the fact that current data sets were not observed in a manner optimized for CMAGIC. We describe the details of our novel blindness procedure, which is designed to prevent experimenter bias. The data are broadly consistent with the picture of an accelerating Universe, and agree with a flat Universe within 1.7\sigma, including systematics. We also compare the CMAGIC results directly with those of a maximum magnitude fit to the same SNe, finding that CMAGIC favors more acceleration at the 1.6\sigma level, including systematics and the correlation between the two measurements. A fit for w assuming a flat Universe yields a value which is consistent with a cosmological constant within 1.2\sigma.

(ABRIDGED) We show that Type Ia supernovae (SNe Ia) are formed within both very young and old stellar populations, with observed rates that depend on the stellar mass and mean star-formation rates (SFRs) of their host galaxies. Models where the SN Ia rate depends solely on host galaxy stellar mass are ruled out with 99% confidence. Our analysis is based on 100 spectroscopically-confirmed SNe Ia, plus 24 photometrically-classified events, all from the Supernova Legacy Survey (SNLS) and distributed over 0.2<z<0.75. Using multi-band photometry, we estimate stellar masses and SFRs for the SN Ia host galaxies by fitting their broad-band spectral energy distributions with the galaxy spectral synthesis code, PEGASE.2. We show that the SN Ia rate per unit mass is proportional to the specific SFR of the parent galaxies -- more vigorously star-forming galaxies host more SNe Ia per unit stellar mass, broadly equivalent to the trend of increasing SN Ia rate in later-type galaxies seen in the local universe. Following earlier suggestions for a simple "two-component" model approximating the SN Ia rate, we find bivariate linear dependencies of the SN Ia rate on both the stellar masses and the mean SFRs of the host systems. We also demonstrate a dependence of distant SN Ia light-curve shapes on star-formation in the host galaxy, similar to trends observed locally. Passive galaxies, with no star-formation, preferentially host faster-declining/dimmer SNe Ia, while slower-declining/brighter events are only found in systems with ongoing star-formation. We model the light-curve width distribution in star-forming galaxies as the sum of a young component, and an old component taken from the distribution in non-star-forming galaxies.

We present the first high-redshift Hubble diagram for Type II-P supernovae (SNe II-P) based upon five events at redshift up to z~0.3. This diagram was constructed using photometry from the Canada-France-Hawaii Telescope Supernova Legacy Survey and absorption line spectroscopy from the Keck observatory. The method used to measure distances to these supernovae is based on recent work by Hamuy & Pinto (2002) and exploits a correlation between the absolute brightness of SNe II-P and the expansion velocities derived from the minimum of the Fe II 516.9 nm P-Cygni feature observed during the plateau phases. We present three refinements to this method which significantly improve the practicality of measuring the distances of SNe II-P at cosmologically interesting redshifts. These are an extinction correction measurement based on the V-I colors at day 50, a cross-correlation measurement for the expansion velocity and the ability to extrapolate such velocities accurately over almost the entire plateau phase. We apply this revised method to our dataset of high-redshift SNe II-P and find that the resulting Hubble diagram has a scatter of only 0.26 magnitudes, thus demonstrating the feasibility of measuring the expansion history, with present facilities, using a method independent of that based upon supernovae of Type Ia.

Published B and V fluxes from nearby Type Ia supernovae are fitted to light-curve templates with 4-6 adjustable parameters. Separately, B magnitudes from the same sample are fitted to a linear dependence on B-V color within a post-maximum time window prescribed by the CMAGIC method. These fits yield two independent SN magnitude estimates B_max and B_BV. Their difference varies systematically with decline rate Delta m_15 in a form that is compatible with a bilinear but not a linear dependence; a nonlinear form likely describes the decline-rate dependence of B_max itself. A Hubble fit to the average of B_max and B_BV requires a systematic correction for observed B-V color that can be described by a linear coefficient R = 2.59 +- 0.24, well below the coefficient R_B ~ 4.1 commonly used to characterize the effects of Milky Way dust. At 99.9% confidence the data reject a simple model in which no color correction is required for SNe that are clustered at the blue end of their observed color distribution. After systematic corrections are performed, B_max and B_BV exhibit mutual rms intrinsic variation equal to 0.074 +- 0.019 mag, of which at least an equal share likely belongs to B_BV. SN magnitudes measured using maximum-luminosity or CMAGIC methods show comparable rms deviations of order ~ 0.14 mag from the Hubble line. The same fit also establishes a 95% confidence upper limit of 486 km/s on the rms peculiar velocity of nearby SNe relative to the Hubble flow.

We present a method for selecting high-redshift type Ia supernovae (SNe Ia) located via rolling SN searches. The technique, using both color and magnitude information of events from only 2-3 epochs of multi-band real-time photometry, is able to discriminate between SNe Ia and core collapse SNe. Furthermore, for the SNe Ia, the method accurately predicts the redshift, phase and light-curve parameterization of these events based only on pre-maximum-light data. We demonstrate the effectiveness of the technique on a simulated survey of SNe Ia and core-collapse SNe, where the selection method effectively rejects most core-collapse SNe while retaining SNe Ia. We also apply the selection code to real-time data acquired as part of the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). During the period May 2004 to January 2005 in the SNLS, 440 SN candidates were discovered of which 70 were confirmed spectroscopically as SNe Ia and 15 as core-collapse events. For this test dataset, the selection technique correctly identifies 100% of the identified SNe II as non-SNe Ia with only a 1-2% false rejection rate. The predicted parameterization of the SNe Ia has a precision of |delta_z|/(1+z_spec)<0.09 in redshift, and +/- 2-3 rest-frame days in phase, providing invaluable information for planning spectroscopic follow-up observations. We also investigate any bias introduced by this selection method on the ability of surveys such as SNLS to measure cosmological parameters (e.g., w and omega matter), and find any effect to be negligible.

We present spectra for 14 high-redshift (0.17 < z < 0.83) supernovae, which were discovered by the Supernova Cosmology Project as part of a campaign to measure cosmological parameters. The spectra are used to determine the redshift and classify the supernova type, essential information if the supernovae are to be used for cosmological studies. Redshifts were derived either from the spectrum of the host galaxy or from the spectrum of the supernova itself. We present evidence that these supernovae are of Type Ia by matching to spectra of nearby supernovae. We find that the dates of the spectra relative to maximum light determined from this fitting process are consistent with the dates determined from the photometric light curves, and moreover the spectral time-sequence for SNe Type Ia at low and high redshift is indistinguishable. We also show that the expansion velocities measured from blueshifted CaHK are consistent with those measured for low-redshift Type Ia supernovae. From these first-level quantitative comparisons we find no evidence for evolution in SNIa properties between these low- and high-redshift samples. Thus even though our samples may not be complete, we conclude that there is a population of SNe Ia at high redshift whose spectral properties match those at low redshift.

We present optical spectra of the peculiar Type Ia supernova (SN Ia) 1999ac. The data extend from -15 to +42 days with respect to B-band maximum and reveal an event that is unusual in several respects. Prior to B-band maximum, the spectra resemble those of SN 1999aa, a slowly declining event, but possess stronger SiII and CaII signatures (more characteristic of a spectroscopically normal SN). Spectra after B-band maximum appear more normal. The expansion velocities inferred from the Iron lines appear to be lower than average; whereas, the expansion velocity inferred from Calcium H and K are higher than average. The expansion velocities inferred from SiII are among the slowest ever observed, though SN 1999ac is not particularly dim. The analysis of the parameters v_10, R(SiII), dv(SiII)/dt, and d_m15 further underlines the unique characteristics of SN 1999ac. We find convincing evidence of CII 6580 in the day -15 spectrum with ejection velocity v > 16,000 km/s, but this signature disappears by day -9. This rapid evolution at early times highlights the importance of extremely early-time spectroscopy.

We present a novel technique for fitting restframe I-band light curves on a data set of 42 Type Ia supernovae (SNe Ia). Using the result of the fit, we construct a Hubble diagram with 26 SNe from the subset at 0.01< z<0.1. Adding two SNe at z~0.5 yields results consistent with a flat Lambda-dominated``concordance universe'' ($\Omega_M,\Omega_\Lambda$)=(0.25,0.75). For one of these, SN 2000fr, new near infrared data are presented. The high redshift supernova NIR data are also used to test for systematic effects in the use of SNe Ia as distance estimators. A flat, Lambda=0, universe where the faintness of supernovae at z~0.5 is due to grey dust homogeneously distributed in the intergalactic medium is disfavoured based on the high-z Hubble diagram using this small data-set. However, the uncertainties are large and no firm conclusion may be drawn. We explore the possibility of setting limits on intergalactic dust based on B-I and B-V colour measurements, and conclude that about 20 well measured SNe are needed to give statistically significant results. We also show that the high redshift restframe I-band data points are better fit by light curve templates that show a prominent second peak, suggesting that they are not intrinsically underluminous.

We present VLT FORS1 and FORS2 spectra of 39 candidate high-redshift supernovae that were discovered as part of a cosmological study using Type Ia supernovae (SNe Ia) over a wide range of redshifts. From the spectra alone, 20 candidates are spectrally classified as SNe Ia with redshifts ranging from z=0.212 to z=1.181. Of the remaining 19 candidates, 1 might be a Type II supernova and 11 exhibit broad supernova-like spectral features and/or have supernova-like light curves. The candidates were discovered in 8 separate ground-based searches. In those searches in which SNe Ia at z ~ 0.5 were targeted, over 80% of the observed candidates were spectrally classified as SNe Ia. In those searches in which SNe Ia with z > 1 were targeted, 4 candidates with z > 1 were spectrally classified as SNe Ia and later followed with ground and space based observatories. We present the spectra of all candidates, including those that could not be spectrally classified as supernova.

We present the type Ia rate measurement based on two EROS supernova search campaigns (in 1999 and 2000). Sixteen supernovae identified as type Ia were discovered. The measurement of the detection efficiency, using a Monte Carlo simulation, provides the type Ia supernova explosion rate at a redshift ~ 0.13. The result is $0.125^{+0.044+0.028}_{-0.034-0.028} h_{70}^2$ SNu where 1 SNu = 1 SN / $10^{10} L_{sun}^B$ / century. This value is compatible with the previous EROS measurement (Hardin et al. 2000), done with a much smaller sample, at a similar redshift. Comparison with other values at different redshifts suggests an evolution of the type Ia supernova rate.

We present an extensive new time-series of spectroscopic data of the peculiar SN 1999aa in NGC 2595. Our data set includes 25 optical spectra between -11 and +58 days with respect to B-band maximum light, providing an unusually complete time history. The early spectra resemble those of a SN 1991T-like object but with a relatively strong Ca H&K absorption feature. The first clear sign of Si II 6355, characteristic of Type Ia supernovae, is found at day -7 and its velocity remains constant up to at least the first month after B-band maximum light. The transition to normal-looking spectra is found to occur earlier than in SN 1991T suggesting SN 1999aa as a possible link between SN 1991T-like and Branch-normal supernovae. Comparing the observations with synthetic spectra, doubly ionized Fe, Si and Ni are identified at early epochs. These are characteristic of SN 1991T-like objects. Furthermore, in the day -11 spectrum, evidence is found for an absorption feature which could be identified as high velocity C II 6580 or H-alpha. At the same epoch C III 4648.8 at photospheric velocity is probably responsible for the absorption feature at 4500 A. High velocity Ca is found around maximum light together with Si II and Fe II confined in a narrow velocity window. Implied constraints on supernovae progenitor systems and explosion hydrodynamical models are briefly discussed.

We report measurements of $\Omega_M$, $\Omega_\Lambda$, and w from eleven supernovae at z=0.36-0.86 with high-quality lightcurves measured using WFPC-2 on the HST. This is an independent set of high-redshift supernovae that confirms previous supernova evidence for an accelerating Universe. Combined with earlier Supernova Cosmology Project data, the new supernovae yield a flat-universe measurement of the mass density $\Omega_M=0.25^{+0.07}_{-0.06}$ (statistical) $\pm0.04$ (identified systematics), or equivalently, a cosmological constant of $\Omega_\Lambda=0.75^{+0.06}_{-0.07}$ (statistical) $\pm0.04$ (identified systematics). When the supernova results are combined with independent flat-universe measurements of $\Omega_M$ from CMB and galaxy redshift distortion data, they provide a measurement of $w=-1.05^{+0.15}_{-0.20}$ (statistical) $\pm0.09$ (identified systematic), if w is assumed to be constant in time. The new data offer greatly improved color measurements of the high-redshift supernovae, and hence improved host-galaxy extinction estimates. These extinction measurements show no anomalous negative E(B-V) at high redshift. The precision of the measurements is such that it is possible to perform a host-galaxy extinction correction directly for individual supernovae without any assumptions or priors on the parent E(B-V) distribution. Our cosmological fits using full extinction corrections confirm that dark energy is required with $P(\Omega_\Lambda>0)>0.99$, a result consistent with previous and current supernova analyses which rely upon the identification of a low-extinction subset or prior assumptions concerning the intrinsic extinction distribution.

(Abridged) We present new results on the Hubble diagram of distant type Ia supernovae (SNe Ia) segregated according to the type of host galaxy. This makes it possible to check earlier evidence for a cosmological constant by explicitly comparing SNe residing in galaxies likely to contain negligible dust with the larger sample. The cosmological parameters derived from these SNe Ia hosted by presumed dust-free early-type galaxies supports earlier claims for a cosmological constant, which we demonstrate at 5 sigma significance, and the internal extinction implied is small even for late-type systems (A_B<0.2). Thus, our data demonstrate that host galaxy extinction is unlikely to systematically dim distant SNe Ia in a manner that would produce a spurious cosmological constant. We classify the host galaxies of 39 distant SNe discovered by the Supernova Cosmology Project (SCP) using the combination of HST STIS imaging, Keck ESI spectroscopy and ground-based broad-band photometry. We compare with a low-redshift sample of 25 SNe Ia. The scatter observed in the SNe Ia Hubble diagrams correlates closely with host galaxy morphology. We find the scatter in the SNe Ia Hubble diagram is smallest for SNe occurring in early-type hosts and largest for those occurring in late-type galaxies. Moreover, SNe residing in early-type hosts appear only ~0.14+/-0.09 mag brighter in their light-curve-width-corrected luminosity than those in late-type hosts, implying only a modest amount of dust extinction even in the late-type systems.

This paper provides a progress report on a collaborative program at the Las Campanas and Cerro Tololo Observatories to observe the near-IR light curves of Type Ia supernovae. We discuss how the morphologies of the JHK light curves change as a function of the decline rate. Evidence is presented which indicates that the absolute magnitudes in the H band have little or no dependence on the decline rate, suggesting that SNe Ia may be nearly perfect cosmological standard candles in the near-IR. A preliminary Hubble diagram in the H band is presented and compared with a similar diagram in V for the same objects. Finally, observations of two peculiar supernovae, 1999ac and 2001ay, are briefly discussed.

SN 1999aw was discovered during the first campaign of the Nearby Galaxies Supernova Search (NGSS) project. This luminous, slow-declining (Delta m_{15} (B) = 0.81 \pm 0.03) Type Ia supernova was noteworthy in at least two respects. First, it occurred in an extremely low luminosity host galaxy that was not visible in the template images, nor in initial subsequent deep imaging. Secondly, the photometric and spectral properties of this supernova indicate that it very likely was similar to the subclass of Type Ia supernovae whose prototype is SN 1999aa. This paper presents the BVRI and J_{s} HK_{s} lightcurves of SN 1999aw (through 100 days past maximum light), as well as several epochs of optical spectra. From these data we calculate the bolometric light curve, and give estimates of the luminosity at maximum light and the initial ^{56}Ni mass. In addition, we present deep BVI images obtained recently with the Baade 6.5-meter telescope at Las Campanas Observatory which reveal the remarkably low-luminosity host galaxy.

R-band intensity measurements along the light curve of Type Ia supernovae discovered by the Supernova Cosmology Project (SCP) are fitted in brightness to templates allowing a free parameter the time-axis width factor w = s(1+z). The data points are then individually aligned in the time-axis, normalized and K-corrected back to the rest frame, after which the nearly 1300 normalized intensity measurements are found to lie on a well-determined common rest-frame B-band curve which we call the ``composite curve''. The same procedure is applied to 18 low-redshift Calan/Tololo SNe with z < 0.11; these nearly 300 B-band photometry points are found to lie on the composite curve equally well. The SCP search technique produces several measurements before maximum light for each supernova. We demonstrate that the linear stretch factor, s, which parameterizes the light-curve timescale appears independent of z,and applies equally well to the declining and rising parts of the light curve. In fact, the B-band template that best fits this composite curve fits the individual supernova photometry data when stretched by a factor s with chi^2/DoF approx = 1, thus as well as any parameterization can, given the current data sets. The measurement of the date of explosion, however, is model dependent and not tightly constrained by the current data. We also demonstrate the 1+z light-curve time-axis broadening expected from cosmological expansion. This argues strongly against alternative explanations, such as tired light, for the redshift of distant objects.