This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We measure the imprint of primordial baryon acoustic oscillations (BAO) in the correlation function of Ly$\alpha$ absorption in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS) in Data Release 14 (DR14) of the Sloan Digital Sky Survey (SDSS)-IV. In addition to 179,965 spectra with absorption in the Lyman-$\alpha$ (Ly$\alpha$) region, we use, for the first time, Ly$\alpha$ absorption in the Lyman-$\beta$ region of 56,154 spectra. We measure the Hubble distance, $D_H$, and the comoving angular diameter distance, $D_M$, relative to the sound horizon at the drag epoch $r_d$ at an effective redshift $z=2.34$. Using a physical model of the correlation function outside the BAO peak, we find $D_H(2.34)/r_d=8.86\pm 0.29$ and $D_M(2.34)/r_d=37.41\pm 1.86$, within 1$\sigma$ from the flat-$\Lambda$CDM model consistent with CMB anisotropy measurements. With the addition of polynomial "broadband" terms, the results remain within one standard deviation of the CMB-inspired model. Combined with the quasar-Ly$\alpha$ cross-correlation measurement presented in a companion paper Blomqvist19, the BAO measurements at $z=2.35$ are within 1.7$\sigma$ of the predictions of this model.

We present a measurement of the baryon acoustic oscillation (BAO) scale at redshift $z=2.35$ from the three-dimensional correlation of Lyman-$\alpha$ (Ly$\alpha$) forest absorption and quasars. The study uses 266,590 quasars in the redshift range $1.77<z<3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 14 (DR14). The sample includes the first two years of observations by the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS), providing new quasars and re-observations of BOSS quasars for improved statistical precision. Statistics are further improved by including Ly$\alpha$ absorption occurring in the Ly$\beta$ wavelength band of the spectra. From the measured BAO peak position along and across the line of sight, we determined the Hubble distance $D_{H}$ and the comoving angular diameter distance $D_{M}$ relative to the sound horizon at the drag epoch $r_{d}$: $D_{H}(z=2.35)/r_{d}=9.20\pm 0.36$ and $D_{M}(z=2.35)/r_{d}=36.3\pm 1.8$. These results are consistent at $1.5\sigma$ with the prediction of the best-fit spatially-flat cosmological model with the cosmological constant reported for the Planck (2016) analysis of cosmic microwave background anisotropies. Combined with the Ly$\alpha$ auto-correlation measurement presented in a companion paper, the BAO measurements at $z=2.34$ are within $1.7\sigma$ of the predictions of this model.

We present the first attempt at measuring the baryonic acoustic oscillations (BAO) in the large scale cross-correlation between the magnesium-II doublet (MgII) flux transmission field and the position of quasar and galaxy tracers. The MgII flux transmission continuous field at $0.3 < z < 1.6$ is measured from 500,589 quasar spectra obtained in the Baryonic Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS). The position of 246,697 quasar tracers and 1,346,776 galaxy tracers are extracted from the Sloan Digital Sky Survey (SDSS) I, II, BOSS and eBOSS catalogs. In addition to measuring the cosmological BAO scale and the biased matter density correlation, this study allows tests and improvements to cosmological Lyman-$\alpha$ analyses. A feature consistent with that of the BAO is detected at a significance of $\Delta \chi^{2} = 7.25$. The measured MgII linear transmission bias parameters are $b_{\mathrm{MgII(2796)}} (z = 0.59) = (-6.82 \pm 0.54) \, \times 10^{-4}$ and $b_{\mathrm{MgII(2804)}} (z = 0.59) = (-5.55 \pm 0.46) \, \times 10^{-4}$, and the MgI bias is $b_{\mathrm{MgI(2853)}} (z = 0.59) = (-1.48 \pm 0.24) \, \times 10^{-4}$. Their redshift evolution is characterized by the power-law index: $\gamma_{\mathrm{Mg}} = 3.36 \pm 0.46$. These measurements open a new window towards using BAO from flux transmission at $z < 1$ in the final eBOSS sample and in the upcoming sample from the Dark Energy Spectroscopic Instrument.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

We introduce QuasarNET, a deep convolutional neural network that performs classification and redshift estimation of astrophysical spectra with human-expert accuracy. We pose these two tasks as a \emph{feature detection} problem: presence or absence of spectral features determines the class, and their wavelength determines the redshift, very much like human-experts proceed. When ran on BOSS data to identify quasars through their emission lines, QuasarNET defines a sample $99.51\pm0.03$\% pure and $99.52\pm0.03$\% complete, well above the requirements of many analyses using these data. QuasarNET significantly reduces the problem of line-confusion that induces catastrophic redshift failures to below 0.2\%. We also extend QuasarNET to classify spectra with broad absorption line (BAL) features, achieving an accuracy of $98.0\pm0.4$\% for recognizing BAL and $97.0\pm0.2$\% for rejecting non-BAL quasars. QuasarNET is trained on data of low signal-to-noise and medium resolution, typical of current and future astrophysical surveys, and could be easily applied to classify spectra from current and upcoming surveys such as eBOSS, DESI and 4MOST.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present a fully consistent catalog of local and global properties of host galaxies of 882 Type Ia supernov\ae\ (SNIa) that were selected based on their light-curve properties, spanning the redshift range $0.01 < z < 1.\text{}$ This catalog corresponds to a preliminary version of the compilation sample and includes Supernova Legacy Survey (SNLS) 5-year data, Sloan Digital Sky Survey (SDSS), and low-redshift surveys. We measured low- and moderate-redshift host galaxy photometry in SNLS and SDSS images and used spectral energy distribution (SED) fitting techniques to derive host properties such as stellar mass and $U-V$ rest-frame colors; the latter are an indicator of the luminosity-weighted age of the stellar population in a galaxy. We also estimated the local observed fluxes at the supernova location within a proper distance radius of 3 kpc, and transposed them into local $U-V$ rest-frame colors. Selecting SNIa based on host photometry quality, we then performed cosmological fits using local color as a third standardization variable, for which we split the sample at the median value. We find a local color step as significant as the maximum mass step effect. Correcting for the maximum mass step correction, we still find a significant local color effect, which shows that additional information is provided by the close environment of SNIa. Departures from the initial choices were investigated, and we discuss the possible implications for cosmology. This will be of tremendous importance for the forthcoming SNIa surveys, and in particular for the Large Synoptic Survey Telescope (LSST), for which uncertainties on the dark energy equation of state will be comparable to the effects reported here.

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.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We present a measurement of the volumetric rate of superluminous supernovae (SLSNe) at z~1, measured using archival data from the first four years of the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). We develop a method for the photometric classification of SLSNe to construct our sample. Our sample includes two previously spectroscopically-identified objects, and a further new candidate selected using our classification technique. We use the point-source recovery efficiencies from Perrett et.al. (2010) and a Monte Carlo approach to calculate the rate based on our SLSN sample. We find that the three identified SLSNe from SNLS give a rate of 91 (+76/-36) SNe/Yr/Gpc^3 at a volume-weighted redshift of z=1.13. This is equivalent to 2.2 (+1.8/-0.9) x10^-4 of the volumetric core collapse supernova rate at the same redshift. When combined with other rate measurements from the literature, we show that the rate of SLSNe increases with redshift in a manner consistent with that of the cosmic star formation history. We also estimate the rate of ultra-long gamma ray bursts (ULGRBs) based on the events discovered by the Swift satellite, and show that it is comparable to the rate of SLSNe, providing further evidence of a possible connection between these two classes of events. We also examine the host galaxies of the SLSNe discovered in SNLS, and find them to be consistent with the stellar-mass distribution of other published samples of SLSNe.

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.

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

We forecast dark energy constraints that could be obtained from a new large sample of Type Ia supernovae where those at high redshift are acquired with the Euclid space mission. We simulate a three-prong SN survey: a z<0.35 nearby sample (8000 SNe), a 0.2<z<0.95 intermediate sample (8800 SNe), and a 0.75<z<1.55 high-z sample (1700 SNe). The nearby and intermediate surveys are assumed to be conducted from the ground, while the high-z is a joint ground- and space-based survey. This latter survey, the "Dark Energy Supernova Infra-Red Experiment" (DESIRE), is designed to fit within 6 months of Euclid observing time, with a dedicated observing program. We simulate the SN events as they would be observed in rolling-search mode by the various instruments, and derive the quality of expected cosmological constraints. We account for known systematic uncertainties, in particular calibration uncertainties including their contribution through the training of the supernova model used to fit the supernovae light curves. Using conservative assumptions and a 1-D geometric Planck prior, we find that the ensemble of surveys would yield competitive constraints: a constant equation of state parameter can be constrained to sigma(w)=0.022, and a Dark Energy Task Force figure of merit of 203 is found for a two-parameter equation of state. Our simulations thus indicate that Euclid can bring a significant contribution to a purely geometrical cosmology constraint by extending a high-quality SN Hubble diagram to z~1.5. We also present other science topics enabled by the DESIRE Euclid observations

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 present a combined photometric calibration of the SNLS and the SDSS supernova survey, which results from a joint effort of the SDSS and the SNLS collaborations. We deliver fluxes calibrated to the HST spectrophotometric star network for large sets of tertiary stars that cover the science fields of both surveys in all photometric bands. We also cross-calibrate directly the two surveys and demonstrate their consistency. For each survey the flat-fielding is revised based on the analysis of dithered star observations. The calibration transfer from the HST spectrophotometric standard stars to the multi-epoch tertiary standard star catalogs in the science fields follows three different paths: observations of primary standard stars with the SDSS PT telescope; observations of Landolt secondary standard stars with SNLS MegaCam instrument at CFHT; and direct observation of faint HST standard stars with MegaCam. In addition, the tertiary stars for the two surveys are cross-calibrated using dedicated MegaCam observations of stripe 82. This overlap enables the comparison of these three calibration paths and justifies using their combination to improve the calibration accuracy. Flat-field corrections have improved the uniformity of each survey as demonstrated by the comparison of photometry in overlapping fields: the rms of the difference between the two surveys is 3 mmag in gri, 4 mmag in z and 8 mmag in u. We also find a remarkable agreement (better than 1%) between the SDSS and the SNLS calibration in griz. The cross-calibration and the introduction of direct calibration observations bring redundancy and strengthen the confidence in the resulting calibration. We conclude that the surveys are calibrated to the HST with a precision of about 0.4% in griz. This precision is comparable to the external uncertainty affecting the color of the HST primary standard stars.

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.

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 study the pseudo equivalent width of the Si II 4000 feature of Type Ia supernovae in the redshift range 0.0024 < z < 0.634. We find that this spectral indicator correlates with the lightcurve color excess (SALT2 c) as well as previously defined spectroscopic subclasses (Branch types) and the evolution of the Si II 6150 velocity, i.e., the so called velocity gradient. Based on our study of 55 objects from different surveys, we find indications that the Si II 4000 spectral indicator could provide important information to improve cosmological distance measurements with Type Ia supernovae.

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 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 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 present a forecast of dark energy constraints that could be obtained from a large sample of distances to Type Ia supernovae detected and measured from space. We simulate the supernova events as they would be observed by a EUCLID-like telescope with its two imagers, assuming those would be equipped with 4 visible and 3 near infrared swappable filters. We account for known systematic uncertainties affecting the cosmological constraints, including those arising through the training of the supernova model used to fit the supernovae light curves. Using conservative assumptions and Planck priors, we find that a 18 month survey would yield constraints on the dark energy equation of state comparable to the cosmic shear approach in EUCLID: a variable two-parameter equation of state can be constrained to ~0.03 at z~0.3. These constraints are derived from distances to about 13,000 supernovae out to z=1.5, observed in two cones of 10 and 50 deg^2. These constraints do not require measuring a nearby supernova sample from the ground. Provided swappable filters can be accommodated on EUCLID, distances to supernovae can be measured from space and contribute to obtain the most precise constraints on dark energy properties.

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.

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.

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

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

We combine measurements of weak gravitational lensing from the CFHTLS-Wide survey, supernovae Ia from CFHT SNLS and CMB anisotropies from WMAP5 to obtain joint constraints on cosmological parameters, in particular, the dark energy equation of state parameter w. We assess the influence of systematics in the data on the results and look for possible correlations with cosmological parameters. We implement an MCMC algorithm to sample the parameter space of a flat CDM model with a dark-energy component of constant w. Systematics in the data are parametrised and included in the analysis. We determine the influence of photometric calibration of SNIa data on cosmological results by calculating the response of the distance modulus to photometric zero-point variations. The weak lensing data set is tested for anomalous field-to-field variations and a systematic shape measurement bias for high-z galaxies. Ignoring photometric uncertainties for SNLS biases cosmological parameters by at most 20% of the statistical errors, using supernovae only; the parameter uncertainties are underestimated by 10%. The weak lensing field-to-field variance pointings is 5%-15% higher than that predicted from N-body simulations. We find no bias of the lensing signal at high redshift, within the framework of a simple model. Assuming a systematic underestimation of the lensing signal at high redshift, the normalisation sigma_8 increases by up to 8%. Combining all three probes we obtain -0.10<1+w<0.06 at 68% confidence (-0.18<1+w<0.12 at 95%), including systematic errors. Systematics in the data increase the error bars by up to 35%; the best-fit values change by less than 0.15sigma. [Abridged]

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.

Aim: We present new extraction and identification techniques for supernova (SN) spectra developed within the Supernova Legacy Survey (SNLS) collaboration. Method: The new spectral extraction method takes full advantage of photometric information from the Canada-France-Hawai telescope (CFHT) discovery and reference images by tracing the exact position of the supernova and the host signals on the spectrogram. When present, the host spatial profile is measured on deep multi-band reference images and is used to model the host contribution to the full (supernova + host) signal. The supernova is modelled as a Gaussian function of width equal to the seeing. A chi-square minimisation provides the flux of each component in each pixel of the 2D spectrogram. For a host-supernova separation greater than <~ 1 pixel, the two components are recovered separately and we do not use a spectral template in contrast to more standard analyses. This new procedure permits a clean extraction of the supernova separately from the host in about 70% of the 3rd year ESO/VLT spectra of the SNLS. A new supernova identification method is also proposed. It uses the SALT2 spectrophotometric template to combine the photometric and spectral data. A galaxy template is allowed for spectra for which a separate extraction of the supernova and the host was not possible. Result: These new techniques have been tested against more standard extraction and identification procedures. They permit a secure type and redshift determination in about 80% of cases. The present paper illustrates their performances on a few sample spectra.

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.

We present multi-band light curves and distances for five type Ia supernovae at intermediate redshifts, 0.18<z<0.27. Three telescopes on the Canary Island of La Palma, INT, NOT and JKT, were used for discovery and follow-up of type Ia supernovae in the g' and r' filters. Supernova fluxes were measured by simultaneously fitting a supernova and host galaxy model to the data, and then calibrated using star catalogues from the Sloan Digital Sky Survey. The light curve peak luminosities, corrected for light curve shape and colour, are consistent with the expectations for a flat LambdaCDM universe at the 1.5-sigma level. One supernova in the sample, SN1999dr, shows surprisingly large reddening, considering that it is both located at a significant distance from the core of its host (~4 times the fitted exponential radius) and that the galaxy can be spectroscopically classified as early-type with no signs of ongoing star formation.

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)

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 present spectra of six Type Ia and two Type II supernovae obtained in June 2002 at the William Herschel Telescope during a search for Type Ia supernovae (SNIa) at intermediate redshift. Supernova type identification and phase determination are performed using a fitting technique based on a Xi2 minimization against a series of model templates. The spectra range from z=0.033 to z=0.328, including one spectroscopically underluminous SNIa at z=0.033. This set of spectra significantly increases the sample of well-observed type SNIa supernovae available in the range 0.15< z <0.35. Together with the twelve supernovae observed by our team in 1999 in the same redshift range, they form an homogeneous sample of seventeen type Ia supernovae with comparable signal-to-noise ratio and regular phase sampling in a still largely unexplored region of the redshift space.

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 study intermediate--z SNe Ia using the empirical physical diagrams which enable to learn about those SNe explosions. This information can be very useful to reduce systematic uncertainties of the Hubble diagram of SNe Ia up to high z. The study of the expansion velocities and the measurement of the ratio $\mathcal{R}$(\SiII) allow to subtype those SNe Ia as done for nearby samples. The evolution of this ratio as seen in the diagram $\mathcal{R}$(\SiII)--(t) together with $\mathcal{R}$(\SiII)$_{max}$ versus (B-V)$_{0}$ indicate consistency of the properties at intermediate z compared with local SNe. At intermediate--z, the expansion velocities of Ca II and Si II are similar to the nearby counterparts. This is found in a sample of 6 SNe Ia in the range 0.033$\leq z \leq$0.329 discovered within the {\it International Time Programme} (ITP) of {\it Cosmology and Physics with SNe Ia} during the spring of 2002. Those supernovae were identified using the 4.2m William Herschel Telescope. Two SNe Ia at intermediate z were of the cool FAINT type, one being a SN1986G--like object highly reddened. The $\mathcal{R}$(\SiII) ratio as well as subclassification of the SNe Ia beyond templates help to place SNe Ia in their sequence of brightness and to distinguish between reddened and intrinsically red supernovae. This test can be done with very high z SNe Ia and it will help to reduce systematic uncertainties due to extinction by dust. It should allow to map the high-z sample into the nearby one.

We present distance measurements to 71 high redshift type Ia supernovae discovered during the first year of the 5-year Supernova Legacy Survey (SNLS). These events were detected and their multi-color 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 redshift. With this data set, we have built a Hubble diagram extending to z=1, with all distance measurements involving at least two bands. Systematic uncertainties are evaluated making use of the multi-band photometry obtained at CFHT. Cosmological fits to this first year SNLS Hubble diagram give the following results : Omega_M = 0.263 +/- 0.042(stat) +/- 0.032(sys) for a flat LambdaCDM model; and w = -1.023 +/- 0.090(stat) +/- 0.054(sys) for a flat cosmology with constant equation of state w when combined with the constraint from the recent Sloan Digital Sky Survey measurement of baryon acoustic oscillations.

We present spectra of twelve Type Ia supernovae obtained in 1999 at the William Herschel Telescope and the Nordic Optical Telescope during a search for Type Ia supernovae (SN Ia) at intermediate redshift. The spectra range from z=0.178 to z=0.493, including five high signal-to-noise ratio SN Ia spectra in the still largely unexplored range 0.15 < z < 0.3. Most of the spectra were obtained before or around restframe B-band maximum light. None of them shows the peculiar spectral features found in low-redshift over- or under-luminous SN Ia. Expansion velocities of characteristic spectral absorption features such as SiII at 6355 angs., SII at 5640 angs. and CaII at 3945 angs. are found consistent with their low-z SN Ia counterparts.

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 discuss the expansion velocities of different elements in the ejecta of the intermediate-z SNe Ia discovered as a part of the International Time Programme (ITP) project ``Omega and Lambda from Supernovae and the Physics of Supernova Explosions'' at the European Northern Observatory (ENO). The expansion velocities measured for each normal SNIa are found to be within the typical velocity dispersion for their epoch. Meanwhile, the subluminous SN 2002lk SiII expansion velocity is significantly higher than that of SN 1991bg shortly after maximum. The observed phase was younger in SN2002lk than in the local subluminous SNIa SN1991bg.

We study 8 supernovae discovered as part of the International Time Programme (ITP) project ``Omega and Lambda from Supernovae and the Physics of Supernova Explosions'' at the European Northern Observatory (ENO). The goal of the project is to increase the sample of intermediate redshift (0.1<z<0.4) SNe Ia for testing properties of SNe Ia along z and for enlarging the sample in the Hubble diagram up to large z.

We discuss the host galaxy morphology of the 8 SNe discovered as a part of the International Time Programme (ITP) project ``Omega and Lambda from Supernovae, and the Physics of Supernovae Explosions'' at the European Northern Observatory (ENO). Identification of the SN host galaxy types was done exploiting both imaging and spectroscopy. A peculiar SNIa at z= 0.033 is found in a spiral galaxy, as most other SNeIa with z between 0.1 and 0.4. A complete account of these studies will be given elsewhere.

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.

A semi-analytic model is proposed that couples the Press-Schechter formalism for the number of galaxies with a prescription for galaxy-galaxy interactions that enables to follow the evolution of galaxy morphologies along the Hubble sequence. Within this framework, we calculate the chemo-spectrophotometric evolution of galaxies to obtain spectral energy distributions. We find that such an approach is very successful in reproducing the statistical properties of galaxies as well as their time evolution. We are able to make predictions as a function of galaxy type: for clarity, we restrict ourselves to two categories of galaxies: early and late types that are identified with ellipticals and disks. In our model, irregulars are simply an early stage of galaxy formation. In particular, we obtain good matches for the galaxy counts and redshift distributions of sources from UV to submm wavelengths. We also reproduce the observed cosmic star formation history and the diffuse background radiation, and make predictions as to the epoch and wavelength at which the dust-shrouded star formation of spheroids begins to dominate over the star formation that occurs more quiescently in disks. A new prediction of our model is a rise in the FIR luminosity density with increasing redshift, peaking at about $z\sim 3$, and with a ratio to the local luminosity density $\rho_{L,\nu} (z = z_{peak})/ \rho_{L,\nu} (z = 0)$ about 10 times higher than that in the blue (B-band) which peaks near $z\sim 2$.

Recent studies of galactic nuclei suggest that most galaxies are seeded by super-massive black holes which power the central nucleus. In this picture, the proto-galactic object is likely to have undergone a very active phase during which the surrounding medium was shocked and heated up to very high temperatures. We predict the cosmological effects of this scenario in terms of the thermal and kinetic Sunyaev-Zel'dovich distortions induced on galactic scales by a population of proto-galaxies. These predictions are compared to the observational limit on the mean Compton distortion set by the COBE-FIRAS instrument. This enables us to derive tight constraints on the fraction of proto-galaxies seeded by super-massive black holes as well as on the black hole-to-spheroid mass ratio. Finally, we estimate the contribution of such a population to the angular power spectrum of the Cosmic Microwave Background temperature anisotropies on very small angular scales ($l\simeq 10^4-10^5$).

We present a simple model of the infra-red universe, based as much as possible on local observations. We model the luminosity and number evolution of disk and starburst galaxies, including the effects of dust, gas and spectral evolution. Although simple, our approach is able to reproduce observations of galaxy number counts and the infra-red and sub-millimeter extra-galactic backgrounds. It provides a useful probe of galaxy formation and evolution out to high redshift. The model demonstrates the significant role of the starburst population and predicts high star formation rates at redshifts 3 to 4, consistent with recent extinction-corrected observations of Lyman break galaxies. Starbursting galaxies are predicted to dominate the current SCUBA surveys. Their star formation is driven predominantly by strong tidal interactions and mergers of galaxies. This leads to the creation of spheroidal stellar systems, which may act as the seeds for disk formation as gas infalls. We predict the present-day baryonic mass in bulges and halos is comparable to that in disks. From observations of the extra-galactic background, the model predicts that the vast majority of star formation in the Universe occurs at z<5.

We present a semi-analytical model in which galaxy collisions and strong tidal interactions, both in the field and during the collapse phase of groups and clusters help determine galaxy morphology. From a semi-analytical analysis based on simulation results of tidal collisions (Aguilar & White 1985), we propose simple rules for energy exchanges during collisions that allow to discriminate between different Hubble types: efficient collisions result in the disruption of disks and substantial star formation, leading to the formation of elliptical galaxies; inefficient collisions allow a large gas reservoir to survive and form disks. Assuming that galaxy formation proceeds in a Omega_0=1 Cold Dark Matter universe, the model both reproduces a number of observations and makes predictions, among which are the redshifts of formation of the different Hubble types in the field. When the model is normalized to the present day abundance of X-ray clusters, the amount of energy exchange needed to produce elliptical galaxies in the field implies that they formed at z>2.5 while spiral galaxies formed at z<1.5. The model also offers a natural explanation for biasing between the various morphological types. We find that the present day morphology-density relation in the field is well reproduced under the collision hypothesis. Finally, predictions of the evolution of the various galaxy populations with redshift are made, in the field as well as in clusters.

We study the uncertainty in galaxy cluster mass estimates derived from X-ray data assuming hydrostatic equilibrium (HE) for the intra cluster gas. Using a Monte-Carlo procedure we generate a general class of mass models allowing very massive clusters. We then compute the corresponding temperature profiles via the HE equation and compare them to observational data on some clusters. We find several massive clusters that pass the observational constraints, with integrated masses varying in a quite wide range. The resulting accuracy of the mass estimates is rather poor, larger than what is generally claimed. We argue that the tight constraints on cluster masses previously obtained come from the fact that a too restricted class of mass density profiles has been investigated so far. Applying our procedure to Perseus then Coma, we find that the improvement of the observational constraints results in a quite modest improvement in the accuracy of the mass estimate. For Coma, using the best current available data, we end up with a factor of two of uncertainty in the mass within the Abell radius. This uncertainty rapidly increases at further radius.

The temperature distribution function (TDF) of X-ray clusters is derived for different normalizations of the HDM power spectrum in an Omega_0=1 universe using the statistics of peaks of a random Gaussian field. Concerning cluster formation only, the neutrino picture appears to be marginally consistent with observations for H_0=70 km/s/Mpc provided the normalization is alpha=sigma_0/delta_s=0.42 with sigma_0 the rms variance of the fluctuations linearly evolved to the present epoch and delta_s the collapse threshold required in structure formation models. However, using the constraint that the formation epoch of QSO's puts on the HDM model, it appears that this normalization is too small to be consistent with quasars formation. Allowing higher values for the Hubble constant alleviates this problem but leads to an unacceptable short age of the universe. Moreover, the constraint put by the recent measurement of the level of fluctuations in the cosmic background radiation by COBE tends to reject models with high values of H_0. We conclude that the neutrino model fails in any case to account for the whole set of available observations in a consistent way.