We present SALT2X, an extension of the SALT2 model for SN Ia supernova light curves. SALT2X separates the light-curve-shape parameter x1 into an x1r and x1f for the rise and fall portions of the light curve. Using the Joint Lightcurve Analysis (JLA) SN sample, we assess the importance of the rising and falling portions of the light curve for cosmological standardization using a modified version of the Unified Nonlinear Inference for Type Ia cosmologY (UNITY) framework. We find strong evidence that x1r has a stronger correlation with peak magnitude than x1f. We see evidence that standardizing on the rise affects the color standardization relation, and reduces the size of the host-galaxy standardization and the unexplained ("intrinsic") luminosity dispersion. Since SNe Ia generally rise more quickly than they decline, a faster observing cadence in future surveys will be necessary to maximize the gain from this work, and to continue to explore the impacts of decoupling the rising and falling portions of SN Ia light curves.

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

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

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.

B-band light-curve rise times for eight unusually well-observed nearby Type Ia supernovae (SNe) are fitted by a newly developed template-building algorithm, using light-curve functions that are smooth, flexible, and free of potential bias from externally derived templates and other prior assumptions. From the available literature, photometric BVRI data collected over many months, including the earliest points, are reconciled, combined, and fitted to a unique time of explosion for each SN. On average, after they are corrected for light-curve decline rate, three SNe rise in 18.81 +- 0.36 days, while five SNe rise in 16.64 +- 0.21 days. If all eight SNe are sampled from a single parent population (a hypothesis not favored by statistical tests), the rms intrinsic scatter of the decline-rate-corrected SN rise time is 0.96 +0.52 -0.25 days -- a first measurement of this dispersion. The corresponding global mean rise time is 17.44 +- 0.39 days, where the uncertainty is dominated by intrinsic variance. This value is ~2 days shorter than two published averages that nominally are twice as precise, though also based on small samples. When comparing high-z to low-z SN luminosities for determining cosmological parameters, bias can be introduced by use of a light-curve template with an unrealistic rise time. If the period over which light curves are sampled depends on z in a manner typical of current search and measurement strategies, a two-day discrepancy in template rise time can bias the luminosity comparison by ~0.03 magnitudes.

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

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