Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model is constructed from SNe Ia spectral properties and spectrophotometric data from The Nearby Supernova Factory collaboration. In a first step, a PCA-like method is used on spectral features measured at maximum light, which allows us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties are used to extract the average extinction curve. Third, an interpolation using Gaussian Processes facilitates using data taken at different epochs during the lifetime of a SN Ia and then projecting the data on a fixed time grid. Finally, the three steps are combined to build the SED model as a function of time and wavelength. This is the SUGAR model. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity, the second is correlated with their calcium lines. The addition of these parameters, as well as the high quality the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like ZTF, LSST or WFIRST.

Type Ia supernova cosmology depends on the ability to fit and standardize observations of supernova magnitudes with an empirical model. We present here a series of new models of Type Ia Supernova spectral time series that capture a greater amount of supernova diversity than possible with the models that are currently customary. These are entitled SuperNova Empirical MOdels (\textsc{SNEMO}\footnote{https://snfactory.lbl.gov/snemo}). The models are constructed using spectrophotometric time series from $172$ individual supernovae from the Nearby Supernova Factory, comprising more than $2000$ spectra. Using the available observations, Gaussian Processes are used to predict a full spectral time series for each supernova. A matrix is constructed from the spectral time series of all the supernovae, and Expectation Maximization Factor Analysis is used to calculate the principal components of the data. K-fold cross-validation then determines the selection of model parameters and accounts for color variation in the data. Based on this process, the final models are trained on supernovae that have been dereddened using the Fitzpatrick and Massa extinction relation. Three final models are presented here: \textsc{SNEMO2}, a two-component model for comparison with current Type~Ia models; \textsc{SNEMO7}, a seven component model chosen for standardizing supernova magnitudes which results in a total dispersion of $0.100$~mag for a validation set of supernovae, of which $0.087$~mag is unexplained (a total dispersion of $0.113$~mag with unexplained dispersion of $0.097$~mag is found for the total set of training and validation supernovae); and \textsc{SNEMO15}, a comprehensive $15$ component model that maximizes the amount of spectral time series behavior captured.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SNIa) for cosmology, we have statistically classified a large sample of nearby SNeIa into those located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNeIa in predominantly younger environments are DY=0.163\pm0.029 mag (5.7 sigma) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNeIa in hosts with lower or higher total stellar masses is smaller at DM=0.119\pm0.032 mag (4.5 sigma), for the same set of SNeIa. When fit simultaneously, the environment age offset remains very significant, with DY=0.129\pm0.032 mag (4.0 sigma), while the global stellar mass step is reduced to DM=0.064\pm0.029 mag (2.2 sigma). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Standardization using only the SNeIa in younger environments reduces the total dispersion from 0.142\pm0.008 mag to 0.120\pm0.010 mag. We show that as environment ages evolve with redshift a strong bias on measurement of the dark energy equation of state parameters can develop. Fortunately, data to measure and correct for this effect is likely to be available for many next-generation experiments. [abstract shorten]

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of 2.8+0.6-0.5 (1.10 +- 0.23 mag)---compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from LambdaCDM simulations---making it the most amplified SN Ia discovered behind a galaxy cluster.

Type Ia Supernovae (SNe Ia) are widely used to measure the expansion of the Universe. To perform such measurements the luminosity and cosmological redshift ($z$) of the SNe Ia have to be determined. The uncertainty on $z$ includes an unknown peculiar velocity, which can be very large for SNe Ia in the virialized cores of massive clusters. We determine which SNe Ia exploded in galaxy clusters. We then study how the correction for peculiar velocities of host galaxies inside the clusters improves the Hubble residuals. Using 145 SNe Ia from the Nearby Supernova Factory we found 11 candidates for membership in clusters. To estimate the redshift of a cluster we applied the bi-weight technique. Then, we use the galaxy cluster redshift instead of the host galaxy redshift to construct the Hubble diagram. For SNe Ia inside galaxy clusters the dispersion around the Hubble diagram when peculiar velocities are taken into account is smaller in comparison with a case without peculiar velocity correction, with a $wRMS=0.130\pm0.038$ mag instead of $wRMS=0.137\pm0.036$ mag. The significance of this improvement is 3.58 $\sigma$. If we remove the very nearby Virgo cluster member SN2006X ($z<0.01$) from the analysis, the significance decreases to 1.34 $\sigma$. The peculiar velocity correction is found to be highest for the SNe Ia hosted by blue spiral galaxies, with high local specific star formation rate and smaller stellar mass, seemingly counter to what might be expected given the heavy concentration of old, massive elliptical galaxies in clusters. As expected, the Hubble residuals of SNe Ia associated with massive galaxy clusters improve when the cluster redshift is taken as the cosmological redshift of the SN. This fact has to be taken into account in future cosmological analyses in order to achieve higher accuracy for cosmological redshift measurements. Here we provide an approach to do so.

Context. Observations of Type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions {\lambda}(uNi), {\lambda}(uTi), {\lambda}(uSi) and {\lambda}(uCa). Two of these span the Ca H&K {\lambda}{\lambda} 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. (1) The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 $\pm$ 0.011 mag RMS. A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. (2) Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of lightcurve power. (3) The uCa flux further improves standardization, yielding a 0.086 $\pm$ 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach {\chi}$^2$/dof $\sim$ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.

The scientific yield of current and future optical surveys is increasingly limited by systematic uncertainties in the flux calibration. This is the case for Type Ia supernova (SN Ia) cosmology programs, where an improved calibration directly translates into improved cosmological constraints. Current methodology rests on models of stars. Here we aim to obtain flux calibration that is traceable to state-of-the-art detector-based calibration. We present the SNIFS Calibration Apparatus (SCALA), a color (relative) flux calibration system developed for the SuperNova Integral Field Spectrograph (SNIFS), operating at the University of Hawaii 2.2 m (UH 88) telescope. By comparing the color trend of the illumination generated by SCALA during two commissioning runs, and to previous laboratory measurements, we show that we can determine the light emitted by SCALA with a long-term repeatability better than 1%. We describe the calibration procedure necessary to control for system aging. We present measurements of the SNIFS throughput as estimated by SCALA observations. The SCALA calibration unit is now fully deployed at the UH\,88 telescope, and with it color-calibration between 4000 {\AA} and 9000 {\AA} is stable at the percent level over a one-year baseline.

Correction of Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning 3,300 - 9,200 {\AA}, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E(B-V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (RV , RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands, and compare the 5780 {\AA} band with the dust-to-band ratio for the Milky Way. Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the Milky Way. Furthermore we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al. (1989), O'Donnell (1994), and Fitzpatrick (1999), and find the predictions of Fitzpatrick (1999) fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6$\pm$1.1 Mpc. We compare this result with distance measurements in the literature.

While recent supernova cosmology research has benefited from improved measurements, current analysis approaches are not statistically optimal and will prove insufficient for future surveys. This paper discusses the limitations of current supernova cosmological analyses in treating outliers, selection effects, shape- and color-standardization relations, unexplained dispersion, and heterogeneous observations. We present a new Bayesian framework, called UNITY (Unified Nonlinear Inference for Type-Ia cosmologY), that incorporates significant improvements in our ability to confront these effects. We apply the framework to real supernova observations and demonstrate smaller statistical and systematic uncertainties. We verify earlier results that SNe Ia require nonlinear shape and color standardizations, but we now include these nonlinear relations in a statistically well-justified way. This analysis was primarily performed blinded, in that the basic framework was first validated on simulated data before transitioning to real data. We also discuss possible extensions of the method.

We introduce a method for identifying "twin" Type Ia supernovae, and using them to improve distance measurements. This novel approach to Type Ia supernova standardization is made possible by spectrophotometric time series observations from the Nearby Supernova Factory (SNfactory). We begin with a well-measured set of supernovae, find pairs whose spectra match well across the entire optical window, and then test whether this leads to a smaller dispersion in their absolute brightnesses. This analysis is completed in a blinded fashion, ensuring that decisions made in implementing the method do not inadvertently bias the result. We find that pairs of supernovae with more closely matched spectra indeed have reduced brightness dispersion. We are able to standardize this initial set of SNfactory supernovae to 0.083 +/- 0.012 magnitudes, implying a dispersion of 0.072 +/- 0.010 magnitudes in the absence of peculiar velocities. We estimate that with larger numbers of comparison SNe, e.g, using the final SNfactory spectrophotometric dataset as a reference, this method will be capable of standardizing high-redshift supernovae to within 0.06-0.07 magnitudes. These results imply that at least 3/4 of the variance in Hubble residuals in current supernova cosmology analyses is due to previously unaccounted-for astrophysical differences among the supernovae

NICMOS 2 observations are crucial for constraining distances to most of the existing sample of z > 1 SNe Ia. Unlike the conventional calibration programs, these observations involve long exposure times and low count rates. Reciprocity failure is known to exist in HgCdTe devices and a correction for this effect has already been implemented for high and medium count-rates. However observations at faint count-rates rely on extrapolations. Here instead, we provide a new zeropoint calibration directly applicable to faint sources. This is obtained via inter-calibration of NIC2 F110W/F160W with WFC3 in the low count-rate regime using z ~ 1 elliptical galaxies as tertiary calibrators. These objects have relatively simple near-IR SEDs, uniform colors, and their extended nature gives superior signal-to-noise at the same count rate than would stars. The use of extended objects also allows greater tolerances on PSF profiles. We find ST magnitude zeropoints (after the installation of the NICMOS cooling system, NCS) of 25.296 +- 0.022 for F110W and 25.803 +- 0.023 for F160W, both in agreement with the calibration extrapolated from count-rates 1,000 times larger (25.262 and 25.799). Before the installation of the NCS, we find 24.843 +- 0.025 for F110W and 25.498 +- 0.021 for F160W, also in agreement with the high-count-rate calibration (24.815 and 25.470). We also check the standard bandpasses of WFC3 and NICMOS 2 using a range of stars and galaxies at different colors and find mild tension for WFC3, limiting the accuracy of the zeropoints. To avoid human bias, our cross-calibration was "blinded" in that the fitted zeropoint differences were hidden until the analysis was finalized.

We present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova (SLSN) from the Dark Energy Survey (DES). We briefly discuss the data and search algorithm used to find this event in the first year of DES operations, and outline the spectroscopic data obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm its redshift (z = 0.663 +/- 0.001 based on the host-galaxy emission lines) and likely spectral type (type I). Using this redshift, we find M_U_peak = -21.05 +0.10 -0.09 for the peak, rest-frame U-band absolute magnitude, and find DES13S2cmm to be located in a faint, low metallicity (sub-solar), low stellar-mass host galaxy (log(M/M_sun) = 9.3 +/- 0.3); consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest declining tails (beyond +30 days rest frame past peak), and is the faintest at peak. Moreover, we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of only 0.2-0.3 magnitudes between +25 and +30 days after peak (rest frame) depending on redshift range studied; this could be important for 'standardising' such supernovae, as is done with the more common type Ia. We fit the bolometric light curve of DES13S2cmm with two competing models for SLSNe-I - the radioactive decay of 56Ni, and a magnetar - and find that while the magnetar is formally a better fit, neither model provides a compelling match to the data. Although we are unable to conclusively differentiate between these two physical models for this particular SLSN-I, further DES observations of more SLSNe-I should break this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 days in the rest frame of the supernova.

Using three magnified Type Ia supernovae (SNe Ia) detected behind CLASH clusters, we perform a first pilot study to see whether standardizable candles can be used to calibrate cluster mass maps created from strong lensing observations. Such calibrations will be crucial when next generation HST cluster surveys (e.g. FRONTIER) provide magnification maps that will, in turn, form the basis for the exploration of the high redshift Universe. We classify SNe using combined photometric and spectroscopic observations, finding two of the three to be clearly of type SN Ia and the third probable. The SNe exhibit significant amplification, up to a factor of 1.7 at $\sim5\sigma$ significance (SN-L2). We conducted this as a blind study to avoid fine tuning of parameters, finding a mean amplification difference between SNe and the cluster lensing models of $0.09 \pm 0.09^{stat} \pm 0.05^{sys}$ mag. This impressive agreement suggests no tension between cluster mass models and high redshift standardized SNe Ia. However, the measured statistical dispersion of $\sigma_{\mu}=0.21$ mag appeared large compared to the dispersion expected based on statistical uncertainties ($0.14$). Further work with the supernova and cluster lensing models, post unblinding, reduced the measured dispersion to $\sigma_{\mu}=0.12$. An explicit choice should thus be made as to whether SNe are used unblinded to improve the model, or blinded to test the model. As the lensed SN samples grow larger, this technique will allow improved constraints on assumptions regarding e.g. the structure of the dark matter halo.

We present measurements of the rates of high-redshift Type Ia supernovae derived from the Subaru/XMM-Newton Deep Survey (SXDS). We carried out repeat deep imaging observations with Suprime-Cam on the Subaru Telescope, and detected 1040 variable objects over 0.918 deg$^2$ in the Subaru/XMM-Newton Deep Field. From the imaging observations, light curves in the observed $i'$-band are constructed for all objects, and we fit the observed light curves with template light curves. Out of the 1040 variable objects detected by the SXDS, 39 objects over the redshift range $0.2 < z < 1.4$ are classified as Type Ia supernovae using the light curves. These are among the most distant SN Ia rate measurements to date. We find that the Type Ia supernova rate increase up to $z \sim 0.8$ and may then flatten at higher redshift. The rates can be fitted by a simple power law, $r_V(z)=r_0(1+z)^\alpha$ with $r_0=0.20^{+0.52}_{-0.16}$(stat.)$^{+0.26}_{-0.07}$(syst.)$\times 10^{-4} {\rm yr}^{-1}{\rm Mpc}^{-3}$, and $\alpha=2.04^{+1.84}_{-1.96}$(stat.)$^{+2.11}_{-0.86}$(syst.).

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

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional mass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{\alpha}), we obtain \alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction \alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

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.

We present spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

We present data from the Supernova Cosmology Project for five high redshift Type Ia supernovae (SNe Ia) that were obtained using the NICMOS infrared camera on the Hubble Space Telescope. We add two SNe from this sample to a rest-frame I-band Hubble diagram, doubling the number of high redshift supernovae on this diagram. This I-band Hubble diagram is consistent with a flat universe (Omega_Matter, Omega_Lambda= 0.29, 0.71). A homogeneous distribution of large grain dust in the intergalactic medium (replenishing dust) is incompatible with the data and is excluded at the 5 sigma confidence level, if the SN host galaxy reddening is corrected assuming R_V=1.75. We use both optical and infrared observations to compare photometric properties of distant SNe Ia with those of nearby objects. We find generally good agreement with the expected color evolution for all SNe except the highest redshift SN in our sample (SN 1997ek at z=0.863) which shows a peculiar color behavior. We also present spectra obtained from ground based telescopes for type identification and determination of redshift.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

Using deep two-band imaging from the Hubble Space Telescope, we measure the color-magnitude relations (CMR) of E/S0 galaxies in a set of 9 optically-selected clusters principally from the Red-Sequence Cluster Survey (RCS) at 0.9 < z < 1.23. We find that the mean scatter in the CMR in the observed frame of this set of clusters is 0.049 +/- 0.008, as compared to 0.031 +/- 0.007 in a similarly imaged and identically analyzed X-ray sample at similar redshifts. Single-burst stellar population models of the CMR scatter suggest that the E/S0 population in these RCS clusters truncated their star-formation at z~1.6, some 0.9 Gyrs later than their X-ray E/S0 counterparts which were truncated at z~2.1. The notion that this is a manifestation of the differing evolutionary states of the two populations of cluster galaxies is supported by comparison of the fraction of bulge-dominated galaxies found in the two samples which shows that optically-selected clusters contain a smaller fraction of E/S0 galaxies at the their cores

We present a study of the morphological fractions and color-magnitude relation in the most distant X-ray selected galaxy cluster currently known, XMMXCS J2215.9-1738 at z=1.46, using a combination of optical imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys, and infrared data from the Multi-Object Infrared Camera and Spectrograph, mounted on the 8.2m Subaru telescope. We find that the morphological mix of the cluster galaxy population is similar to clusters at z~1: approximately ~62% of the galaxies identified as likely cluster members are ellipticals or S0s; and ~38% are spirals or irregulars. We measure the color-magnitude relations for the early type galaxies, finding that the slope in the z_850-J relation is consistent with that measured in the Coma cluster, some ~9 Gyr earlier, although the uncertainty is large. In contrast, the measured intrinsic scatter about the color-magnitude relation is more than three times the value measured in Coma, after conversion to rest frame U-V. From comparison with stellar population synthesis models, the intrinsic scatter measurements imply mean luminosity weighted ages for the early type galaxies in J2215.9-1738 of ~3 Gyr, corresponding to the major epoch of star formation coming to an end at z_f = 3-5. We find that the cluster exhibits evidence of the `downsizing' phenomenon: the fraction of faint cluster members on the red sequence expressed using the Dwarf-to-Giant Ratio (DGR) is 0.32+/-0.18. This is consistent with extrapolation of the redshift evolution of the DGR seen in cluster samples at z < 1. In contrast to observations of some other z > 1 clusters, we find a lack of very bright galaxies within the cluster.

We present observations of SCP 06F6, an unusual optical transient discovered during the Hubble Space Telescope Cluster Supernova Survey. The transient brightened over a period of ~100 days, reached a peak magnitude of ~21.0 in both i_775 and z_850, and then declined over a similar timescale. There is no host galaxy or progenitor star detected at the location of the transient to a 3 sigma upper limit of i_775 = 26.4 and z_850 = 26.1, giving a corresponding lower limit on the flux increase of a factor of ~120. Multiple spectra show five broad absorption bands between 4100 AA and 6500 AA and a mostly featureless continuum longward of 6500 AA. The shape of the lightcurve is inconsistent with microlensing. The transient's spectrum, in addition to being inconsistent with all known supernova types, is not matched to any spectrum in the Sloan Digital Sky Survey (SDSS) database. We suggest that the transient may be one of a new class.

The recent robust and homogeneous analysis of the world's supernova distance-redshift data, together with cosmic microwave background and baryon acoustic oscillation data, provides a powerful tool for constraining cosmological models. Here we examine particular classes of scalar field, modified gravity, and phenomenological models to assess whether they are consistent with observations even when their behavior deviates from the cosmological constant Lambda. Some models have tension with the data, while others survive only by approaching the cosmological constant, and a couple are statistically favored over LCDM. Dark energy described by two equation of state parameters has considerable phase space to avoid Lambda and next generation data will be required to constrain such physics.

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.

We present a new measurement of the volumetric rate of Type Ia supernova up to a redshift of 1.7, using the Hubble Space Telescope (HST) GOODS data combined with an additional HST dataset covering the North GOODS field collected in 2004. We employ a novel technique that does not require spectroscopic data for identifying Type Ia supernovae (although spectroscopic measurements of redshifts are used for over half the sample); instead we employ a Bayesian approach using only photometric data to calculate the probability that an object is a Type Ia supernova. This Bayesian technique can easily be modified to incorporate improved priors on supernova properties, and it is well-suited for future high-statistics supernovae searches in which spectroscopic follow up of all candidates will be impractical. Here, the method is validated on both ground- and space-based supernova data having some spectroscopic follow up. We combine our volumetric rate measurements with low redshift supernova data, and fit to a number of possible models for the evolution of the Type Ia supernova rate as a function of redshift. The data do not distinguish between a flat rate at redshift > 0.5 and a previously proposed model, in which the Type Ia rate peaks at redshift >1 due to a significant delay from star-formation to the supernova explosion. Except for the highest redshifts, where the signal to noise ratio is generally too low to apply this technique, this approach yields smaller or comparable uncertainties than previous work.

We present Keck diffraction limited H-band photometry of a z~1.3 Type Ia supernova (SN) candidate, first identified in a Hubble Space Telescope (HST) search for SNe in massive high redshift galaxy clusters. The adaptive optics (AO) data were obtained with the Laser Guide Star facility during four observing runs from September to November 2005. In the analysis of data from the observing run nearest to maximum SN brightness, the SN was found to have a magnitude H=23.9 +/- 0.14 (Vega). We present the H-band (approximately rest-frame R) light curve and provide a detailed analysis of the AO photometric uncertainties. By constraining the aperture correction with a nearby (4" separation) star we achieve 0.14 magnitude photometric precision, despite the spatially varying AO PSF.