A non-local-thermodynamic-equilibrium (NLTE) level population model of the first and second ionisation stages of iron, nickel and cobalt is used to fit a sample of XShooter optical + near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia). From the ratio of the NIR lines to the optical lines limits can be placed on the temperature and density of the emission region. We find a similar evolution of these parameters across our sample. Using the evolution of the Fe II 12$\,$570$\,\mathring{A}\,$to 7$\,$155$\,\mathring{A}\,$line as a prior in fits of spectra covering only the optical wavelengths we show that the 7200$\,\mathring{A}\,$feature is fully explained by [Fe II] and [Ni II] alone. This approach allows us to determine the abundance of Ni II$\,$/$\,$Fe II for a large sample of 130 optical spectra of 58 SNe Ia with uncertainties small enough to distinguish between Chandrasekhar mass (M$_{\text{Ch}}$) and sub-Chandrasekhar mass (sub-M$_{\text{Ch}}$) explosion models. We conclude that the majority (85$\%$) of normal SNe Ia have a Ni/Fe abundance that is in agreement with predictions of sub-M$_{\text{Ch}}$ explosion simulations of $\sim Z_\odot$ progenitors. Only a small fraction (11$\%$) of objects in the sample have a Ni/Fe abundance in agreement with M$_{\text{Ch}}$ explosion models.

We present a comprehensive dataset of optical and near-infrared photometry and spectroscopy of type~Ia supernova (SN) 2016hnk, combined with integral field spectroscopy (IFS) of its host galaxy, MCG -01-06-070, and nearby environment. Properties of the SN local environment are characterized by means of single stellar population synthesis applied to IFS observations taken two years after the SN exploded. SN 2016hnk spectra are compared to other 1991bg-like SNe Ia, 2002es-like SNe Ia, and Ca-rich transients. In addition, abundance stratification modelling is used to identify the various spectral features in the early phase spectral sequence and the dataset is also compared to a modified non-LTE model previously produced for the sublumnious SN 1999by. SN 2016hnk is consistent with being a sub-luminous (M$_{\rm B}=-16.7$ mag, s$_{\rm BV}$=0.43$\pm$0.03), highly reddened object. IFS of its host galaxy reveals both a significant amount of dust at the SN location, as well as residual star formation and a high proportion of old stellar populations in the local environment compared to other locations in the galaxy, which favours an old progenitor for SN 2016hnk. Inspection of a nebular spectrum obtained one year after maximum contains two narrow emission lines attributed to the forbidden [Ca II] $\lambda\lambda$7291,7324 doublet with a Doppler shift of 700 km s$^{-1}$. Based on various observational diagnostics, we argue that the progenitor of SN 2016hnk was likely a near Chandrasekhar-mass ($M_{\rm Ch}$) carbon-oxygen white dwarf that produced 0.108 $M_\odot$ of $^{56}$Ni. Our modeling suggests that the narrow [Ca II] features observed in the nebular spectrum are associated with $^{48}$Ca from electron capture during the explosion, which is expected to occur only in white dwarfs that explode near or at the $M_{\rm Ch}$ limit.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

While Chandrasekhar-mass ($M_{\rm Ch}$) models with a low $^{56}\rm{Ni}$ yield can match the peak luminosities of fast-declining, 91bg-like Type Ia supernovae (SNe Ia), they systematically fail to reproduce their faster light-curve evolution. Here we illustrate the impact of a low ejecta mass on the radiative display of low-luminosity SNe Ia, by comparing a sub-$M_{\rm Ch}$ model resulting from the pure central detonation of a C-O White Dwarf (WD) to a $M_{\rm Ch}$ delayed-detonation model with the same $^{56}\rm{Ni}$ yield of 0.12 M$_{\odot}$. Our sub-$M_{\rm Ch}$ model from a 0.90 M$_{\odot}$ WD progenitor has a $\sim$5 day shorter rise time in the integrated UV-optical-IR (uvoir) luminosity, as well as in the $B$-band, and a $\sim$20 per cent higher peak uvoir luminosity ($\sim$1 mag brighter peak $M_B$). This sub-$M_{\rm Ch}$ model also displays bluer maximum-light colours due to the larger specific heating rate, and larger post-maximum uvoir and $B$-band decline rates. The luminosity decline at nebular times is also more pronounced, reflecting the enhanced escape of gamma rays resulting from the lower density of the progenitor WD. The deficit of stable nickel in the innermost ejecta leads to a notable absence of forbidden lines of [Ni II] in the nebular spectra. In contrast, the $M_{\rm Ch}$ model displays a strong line due to [Ni II] 1.939 $\mu\rm{m}$, which could in principle serve to distinguish between different progenitor scenarios. Our sub-$M_{\rm Ch}$ model offers an unprecedented agreement with optical and near-infrared observations of the 91bg-like SN 1999by, making a strong case for a WD progenitor significantly below the Chandrasekhar-mass limit for this event and other low-luminosity SNe Ia.

We present multiband photometry of 60 spectroscopically-confirmed supernovae (SN): 39 SN II/IIP, 19 IIn, one IIb and one that was originally classified as a IIn but later as a Ibn. Forty-six have only optical photometry, six have only near infrared (NIR) photometry and eight have both optical and NIR. The median redshift of the sample is 0.016. We also present 192 optical spectra for 47 of the 60 SN. All data are publicly available. There are 26 optical and two NIR light curves of SN II/IIP with redshifts z > 0.01, some of which may give rise to useful distances for cosmological applications. All photometry was obtained between 2000 and 2011 at the Fred Lawrence Whipple Observatory (FLWO), via the 1.2m and 1.3m PAIRITEL telescopes for the optical and NIR, respectively. Each SN was observed in a subset of the $u'UBVRIr'i'JHK_s$ bands. There are a total of 2932 optical and 816 NIR light curve points. Optical spectra were obtained using the FLWO 1.5m Tillinghast telescope with the FAST spectrograph and the MMT Telescope with the Blue Channel Spectrograph. Our photometry is in reasonable agreement with other samples from the literature. Comparison with Pan-STARRS shows that two-thirds of our individual star sequences have weighted-mean V offsets within $\pm$0.02 mag. In comparing our standard-system SN light curves with common Carnegie Supernova Project objects using their color terms, we found that roughly three-quarters have average differences within $\pm$0.04 mag. The data from this work and the literature will provide insight into SN II explosions, help with developing methods for photometric SN classification, and contribute to their use as cosmological distance indicators.

(Original) Recent high-resolution spectra of the Type Ia SN 2006X have revealed the presence of time-variable and blueshifted Na I D features, interpreted by Patat et al. as originating in circumstellar material within the progenitor system. The variation seen in SN 2006X induces relatively large changes in the total Na I D equivalent width ($\Delta\rm{EW}\approx 0.5 \unicode{x212B}$ in just over two weeks), that would be detectable at lower resolutions. We have used a large data set comprising 2400 low-resolution spectra of 450 Type Ia supernovae (SNe Ia) obtained by the CfA Supernova Program to search for variable Na I D features. Out of the 31 SNe Ia (including SN 2006X) in which we could have detected similar EW variations, only one other (SN 1999cl) shows variable Na I D features, with an even larger change over a similar ~10-day timescale ($\Delta\rm{EW} = 1.66 \pm 0.21 \unicode{x212B}$). Interestingly, both SN 1999cl and SN 2006X are the two most highly-reddened objects in our sample, raising the possibility that the variability is connected to dusty environments. (Erratum) The large variation in the Na I D equivalent width observed in SN 1999cl results in fact from a measurement error. Our new measurements show that the EW variation is significantly lower, at $0.43 \pm 0.14 \unicode{x212B}$. While the EW variation remains statistically significant (3.1$\sigma$ different from zero), it is now below the detection threshold of 0.5 $\unicode{x212B}$ derived from the Monte Carlo simulations published in the original paper. As a result, SN 1999cl should no longer be considered as an object displaying variable Na I D lines in our study. The fraction of SNe Ia in our sample displaying Na I D lines thus goes from $\sim$6% (2/31) in the original study to $\sim$3% (1/31) in the revised analysis, SN 2006X being the only SN Ia in our sample with variable Na I D lines.

The faster light-curve evolution of low-luminosity Type Ia supernovae (SNe Ia) suggests that they could result from the explosion of white dwarf (WD) progenitors below the Chandrasekhar mass ($M_{\rm Ch}$). Here we present 1D non-LTE time-dependent radiative transfer simulations of pure central detonations of carbon-oxygen WDs with a mass ($M_\rm{tot}$) between 0.88 M$_{\odot}$ and 1.15 M$_{\odot}$, and a $^{56}\rm{Ni}$ yield between 0.08 M$_{\odot}$ and 0.84 M$_{\odot}$. Their lower ejecta density compared to $M_{\rm Ch}$ models results in a more rapid increase of the luminosity at early times and an enhanced $\gamma$-ray escape fraction past maximum light. Consequently, their bolometric light curves display shorter rise times and larger post-maximum decline rates. Moreover, the higher $M(^{56}\rm{Ni})/M_\rm{tot}$ ratio at a given $^{56}\rm{Ni}$ mass enhances the temperature and ionization level in the spectrum-formation region for the less luminous models, giving rise to bluer colours at maximum light and a faster post-maximum evolution of the $B-V$ colour. For sub-$M_{\rm Ch}$ models fainter than $M_B\approx -18.5$ mag at peak, the greater bolometric decline and faster colour evolution lead to a larger $B$-band post-maximum decline rate, $\Delta M_{15}(B)$. In particular, all of our previously-published $M_{\rm Ch}$ models (standard and pulsational delayed detonations) are confined to $\Delta M_{15}(B) < 1.4$ mag, while the sub-$M_{\rm Ch}$ models with $M_\rm{tot}\lesssim 1$ M$_{\odot}$ extend beyond this limit to $\Delta M_{15}(B)\approx 1.65$ mag for a peak $M_B\approx -17$ mag, in better agreement with the observed width-luminosity relation (WLR). Regardless of the precise ignition mechanism, these simulations suggest that fast-declining SNe Ia at the faint end of the WLR could result from the explosion of WDs whose mass is significantly below the Chandrasekhar limit.

Fast-declining Type Ia supernovae (SN Ia) separate into two categories based on their bolometric and near-infrared (NIR) properties. The peak bolometric luminosity ($\mathrm{L_{max}}$), the phase of the first maximum relative to the optical, the NIR peak luminosity and the occurrence of a second maximum in the NIR distinguish a group of very faint SN Ia. Fast-declining supernovae show a large range of peak bolometric luminosities ($\mathrm{L_{max}}$ differing by up to a factor of $\sim$ 8). All fast-declining SN Ia with $\mathrm{L_{max}} < 0.3 \cdot$ 10$^{43}\mathrm{erg s}^{-1}$ are spectroscopically classified as 91bg-like and show only a single NIR peak. SNe with $\mathrm{L_{max}} > 0.5 \cdot$ 10$^{43}\mathrm{erg s}^{-1}$ appear to smoothly connect to normal SN Ia. The total ejecta mass (M$_{ej}$) values for SNe with enough late time data are $\lesssim$1 $M_{\odot}$, indicating a sub-Chandrasekhar mass progenitor for these SNe.

The ESSENCE survey discovered 213 Type Ia supernovae at redshifts 0.1 < z < 0.81 between 2002 and 2008. We present their R and I-band photometry, measured from images obtained using the MOSAIC II camera at the CTIO 4 m Blanco telescope, along with rapid-response spectroscopy for each object. We use our spectroscopic follow-up observations to determine an accurate, quantitative classification and a precise redshift. Through an extensive calibration program we have improved the precision of the CTIO Blanco natural photometric system. We use several empirical metrics to measure our internal photometric consistency and our absolute calibration of the survey. We assess the effect of various potential sources of systematic bias on our measured fluxes, and we estimate that the dominant term in the systematic error budget from the photometric calibration on our absolute fluxes is ~1%.

The increase in the number of Type Ia supernovae (SNe\,Ia) has demonstrated that the population shows larger diversity than has been assumed in the past. The reasons (e.g. parent population, explosion mechanism) for this diversity remain largely unknown. We have investigated a sample of SNe\,Ia near-infrared light curves and have correlated the phase of the second maximum with the bolometric peak luminosity. The peak bolometric luminosity is related to the time of the second maximum (relative to the {\it B} light curve maximum) as follows : $L_{max}(10^{43} erg s^{-1}) = (0.039 \pm 0.004) \times t_2(J)(days) + (0.013 \pm 0.106)$. $^{56}$Ni masses can be derived from the peak luminosity based on Arnett's rule, which states that the luminosity at maximum is equal to instantaneous energy generated by the nickel decay. We check this assumption against recent radiative-transfer calculations of Chandrasekhar-mass delayed detonation models and find this assumption is valid to within 10\% in recent radiative-transfer calculations of Chandrasekhar-mass delayed detonation models. The $L_{max}$ vs. $t_2$ relation is applied to a sample of 40 additional SNe\,Ia with significant reddening ($E(B-V) >$ 0.1 mag) and a reddening-free bolometric luminosity function of SNe~Ia is established. The method is tested with the $^{56}$Ni mass measurement from the direct observation of $\gamma-$rays in the heavily absorbed SN 2014J and found to be fully consistent. Super-Chandrasekhar-mass explosions, in particular SN\,2007if, do not follow the relations between peak luminosity and second IR maximum. This may point to an additional energy source contributing at maximum light. The luminosity function of SNe\,Ia is constructed and is shown to be asymmetric with a tail of low-luminosity objects and a rather sharp high-luminosity cutoff, although it might be influenced by selection effects.

We present 1D non-local thermodynamic equilibrium (non-LTE) time-dependent radiative-transfer simulations of a Chandrasekhar-mass delayed-detonation model which synthesizes 0.51 Msun of 56Ni, and confront our results to the Type Ia supernova (SN Ia) 2002bo over the first 100 days of its evolution. Assuming only homologous expansion, this same model reproduces the bolometric and multi-band light curves, the secondary near-infrared (NIR) maxima, and the optical and NIR spectra. The chemical stratification of our model qualitatively agrees with previous inferences by Stehle et al., but reveals significant quantitative differences for both iron-group and intermediate-mass elements. We show that +/-0.1 Msun (i.e., +/-20 per cent) variations in 56Ni mass have a modest impact on the bolometric and colour evolution of our model. One notable exception is the U-band, where a larger abundance of iron-group elements results in less opaque ejecta through ionization effects, our model with more 56Ni displaying a higher near-UV flux level. In the NIR range, such variations in 56Ni mass affect the timing of the secondary maxima but not their magnitude, in agreement with observational results. Moreover, the variation in the I, J, and K_s magnitudes is less than 0.1 mag within ~10 days from bolometric maximum, confirming the potential of NIR photometry of SNe Ia for cosmology. Overall, the delayed-detonation mechanism in single Chandrasekhar-mass white dwarf progenitors seems well suited for SN 2002bo and similar SNe Ia displaying a broad Si II 6355 A line. Whatever multidimensional processes are at play during the explosion leading to these events, they must conspire to produce an ejecta comparable to our spherically-symmetric model.

We present 645 optical spectra of 73 supernovae (SNe) of Types IIb, Ib, Ic, and broad-lined Ic. All of these types are attributed to the core collapse of massive stars, with varying degrees of intact H and He envelopes before explosion. The SNe in our sample have a mean redshift <cz> = 4200 km/s. Most of these spectra were gathered at the Harvard-Smithsonian Center for Astrophysics (CfA) between 2004 and 2009. For 53 SNe, these are the first published spectra. The data coverage range from mere identification (1-3 spectra) for a few SNe to extensive series of observations (10-30 spectra) that trace the spectral evolution for others, with an average of 9 spectra per SN. For 44 SNe of the 73 SNe presented here, we have well-determined dates of maximum light to determine the phase of each spectrum. Our sample constitutes the most extensive spectral library of stripped-envelope SNe to date. We provide very early coverage (as early as 30 days before V-band max) for photospheric spectra, as well as late-time nebular coverage when the innermost regions of the SNe are visible (as late as 2 years after explosion, while for SN1993J, we have data as late as 11.6 years). This data set has homogeneous observations and reductions that allow us to study the spectroscopic diversity of these classes of stripped SNe and to compare these to SNe associated with gamma-ray bursts. We undertake these matters in follow-up papers.

We explore the physics of SN Ia light curves and spectra using the 1-D non-LTE time-dependent radiative-transfer code CMFGEN. Rather than adjusting ejecta properties to match observations, we select as input one "standard" 1-D Chandrasekhar-mass delayed-detonation hydrodynamical model, and then explore the sensitivity of radiation and gas properties on radiative-transfer modeling assumptions. The correct computation of SN Ia radiation is not exclusively a solution to an "opacity problem", characterized by the treatment of a large number of lines. It is also key to treat important atomic processes consistently. Besides handling line blanketing in non-LTE, we show that including forbidden line transitions of metals is increasingly important for the temperature and ionization of the gas beyond maximum light. Non-thermal ionization and excitation are also critical since they affect the color evolution and the Delta-M15 of our model. While impacting little the bolometric luminosity, a more complete treatment of decay routes leads to enhanced line blanketing, e.g., associated with 48Ti in the U and B bands. Overall, we find that SN Ia radiation properties are influenced in a complicated way by the atomic data we employ, so that obtaining converged results is a challenge. We nonetheless obtain a good match to the golden standard type Ia SN 2005cf in the optical and near-IR, from 5 to 60d after explosion, suggesting that assuming spherical symmetry is not detrimental to SN Ia radiative-transfer modeling at these times. Multi-D effects no doubt matter, but they are perhaps less important than accurately treating non-LTE processes [abridged].

Observations of SN 2011fe at early times reveal an evolution analogous to a fireball model of constant color. In contrast, our unmixed delayed detonations of Chandrasekhar-mass white dwarfs (DDC series) exhibit a faster brightening concomitant with a shift in color to the blue. In this paper, we study the origin of these discrepancies. We find that strong chemical mixing largely resolves the photometric mismatch at early times, but it leads to an enhanced line broadening that contrasts, for example, with the markedly narrow SiII6355A line of SN 2011fe. We also explore an alternative configuration with pulsational-delayed detonations (PDDEL model series). Because of the pulsation, PDDEL models retain more unburnt carbon, have little mass at high velocity, and have a much hotter outer ejecta after the explosion. The pulsation does not influence the inner ejecta, so PDDEL and DDC models exhibit similar radiative properties beyond maximum. However, at early times, PDDEL models show bluer optical colors and a higher luminosity, even for weak mixing. Their early-time radiation is derived primarily from the initial shock-deposited energy in the outer ejecta rather than radioactive decay heating. Furthermore, PDDEL models show short-lived CII lines, reminiscent of SN 2013dy. They typically exhibit lines that are weaker, narrower, and of near-constant width, reminiscent of SN 2011fe. In addition to multi-dimensional effects, varying configurations for such ``pulsations" offer a source of spectral diversity amongst SNe Ia. PDDEL and DDC models also provide one explanation for low- and high-velocity gradient SNe Ia.

The high metal content and fast expansion of supernova (SN) Ia ejecta lead to considerable line overlap in their optical spectra. Uncertainties in composition and ionization further complicate the process of line identification. In this paper, we focus on the 5900A emission feature seen in SN Ia spectra after bolometric maximum, a line which in the last two decades has been associated with [CoIII]5888A or NaID. Using non-LTE time-dependent radiative-transfer calculations based on Chandrasekhar-mass delayed-detonation models, we find that NaID line emission is extremely weak at all post-maximum epochs. Instead, we predict the presence of [CoIII]5888A after maximum in all our SN Ia models, which cover a range from 0.12 to 0.87Msun of 56Ni. We also find that the [CoIII]5888A forbidden line is present within days of bolometric maximum, and strengthens steadily for weeks thereafter. Both predictions are confirmed by observations. Rather than trivial taxonomy, these findings confirm that it is necessary to include forbidden-line transitions in radiative-transfer simulations of SNe Ia, both to obtain the correct ejecta cooling rate and to match observed optical spectra.

The delayed-detonation explosion mechanism applied to a Chandrasekhar-mass white dwarf offers a very attractive model to explain the inferred characteristics of Type Ia supernovae (SNe Ia). The resulting ejecta are chemically stratified, have the same mass and roughly the same asymptotic kinetic energy, but exhibit a range in 56Ni mass. We investigate the contemporaneous photometric and spectroscopic properties of a sequence of delayed-detonation models, characterized by 56Ni masses between 0.18 and 0.81 Msun. Starting at 1d after explosion, we perform the full non-LTE, time-dependent radiative transfer with the code CMFGEN, with an accurate treatment of line blanketing, and compare our results to SNe Ia at bolometric maximum. Despite the 1D treatment, our approach delivers an excellent agreement to observations. We recover the range of SN Ia luminosities, colours, and spectral characteristics from the near-UV to 1 micron, for standard as well as low-luminosity 91bg-like SNe Ia. Our models predict an increase in rise time to peak with increasing 56Ni mass, from ~15 to ~21d, yield peak bolometric luminosities that match Arnett's rule to within 10%, and reproduce the much smaller scatter in near-IR magnitudes compared to the optical. We reproduce the morphology of individual spectral features, the stiff dependence of the R(Si) spectroscopic ratio on 56Ni mass, and the onset of blanketing from TiII/ScII in low-luminosity SNe Ia with a 56Ni mass <0.3 Msun. We find that ionization effects, which often dominate over abundance variations, can produce high-velocity features in CaII lines, even in 1D. Distinguishing between different SN Ia explosion mechanisms is a considerable challenge but the results presented here provide additional support to the viability of the delayed-detonation model.

We present non-LTE time-dependent radiative-transfer simulations of pair-instability supernovae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet (WR) star rotation-free progenitors born in the mass range 160-230Msun, at 10^-4 Zsun. Although subject to uncertainties in convection and stellar mass-loss rates, our initial conditions come from physically-consistent models that treat evolution from the main-sequence, the onset of the pair-production instability, and the explosion phase. With our set of input models characterized by large 56Ni and ejecta masses, and large kinetic energies, we recover qualitatively the Type II-Plateau, II-peculiar, and Ib/c light-curve morphologies, although they have larger peak bolometric luminosities (~10^9 to 10^10 Lsun) and a longer duration (~200d). We discuss the spectral properties for each model during the photospheric and nebular phases, including Balmer lines in II-P and II-pec at early times, the dominance of lines from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthening of metal line blanketing thereafter. Having similar He-core properties, all models exhibit similar post-peak spectra that are strongly blanketed by FeII and FeI lines, characterized by red colors, and that arise from photospheres/ejecta with a temperature of <4000K. Combined with the modest line widths after bolometric peak, these properties contrast with those of known super-luminous SNe suggesting that PISNe are yet to be discovered. Being reddish, PISNe will be difficult to observe at high redshift except when they stem from RSG explosions, in which case they could be used as metallicity probes and distance indicators.

Much uncertainty surrounds the origin of super-luminous supernovae (SNe). Motivated by the discovery of the Type Ic SN2007bi, we study its proposed association with a pair-instability SN (PISN). We compute stellar-evolution models for primordial ~200Msun stars, simulating the implosion/explosion due to the pair-production instability, and use them as inputs for detailed non-LTE time-dependent radiative-transfer simulations that include non-local energy deposition and non-thermal processes. We retrieve the basic morphology of PISN light curves from red-supergiant, blue-supergiant, and Wolf-Rayet (WR) star progenitors. Although we confirm that a progenitor 100Msun helium core (PISN model He100) fits well the SN2007bi light curve, the low ratios of its kinetic energy and 56Ni mass to the ejecta mass, similar to standard core-collapse SNe, conspire to produce cool photospheres, red spectra subject to strong line blanketing, and narrow line profiles, all conflicting with SN2007bi observations. He-core models of increasing 56Ni-to-ejecta mass ratio have bluer spectra, but still too red to match SN2007bi, even for model He125 -- the effect of 56Ni heating is offset by the associated increase in blanketing. In contrast, the delayed injection of energy by a magnetar represents a more attractive alternative to reproduce the blue, weakly-blanketed, and broad-lined spectra of super-luminous SNe. The extra heat source is free of blanketing and is not explicitly tied to the ejecta. Experimenting with a ~9Msun WR-star progenitor, initially exploded to yield a ~1.6B SN Ib/c ejecta but later influenced by tunable magnetar-like radiation, we produce a diversity of blue spectral morphologies reminiscent of SN2007bi, the peculiar Type Ib SN2005bf, and super-luminous SN2005ap-like events.

We present 2603 spectra of 462 nearby Type Ia supernovae (SN Ia) obtained during 1993-2008 through the Center for Astrophysics Supernova Program. Most of the spectra were obtained with the FAST spectrograph at the FLWO 1.5m telescope and reduced in a consistent manner, making data set well suited for studies of SN Ia spectroscopic diversity. We study the spectroscopic and photometric properties of SN Ia as a function of spectroscopic class using the classification schemes of Branch et al. and Wang et al. The width-luminosity relation appears to be steeper for SN Ia with broader lines. Based on the evolution of the characteristic Si II 6355 line, we propose improved methods for measuring velocity gradients, revealing a larger range than previously suspected, from ~0 to ~400 km/s/day considering the instantaneous velocity decline rate at maximum light. We find a weaker and less significant correlation between Si II velocity and intrinsic B-V color at maximum light than reported by Foley et al., owing to a more comprehensive treatment of uncertainties and host galaxy dust. We study the extent of nuclear burning and report new detections of C II 6580 in 23 early-time spectra. The frequency of C II detections is not higher in SN Ia with bluer colors or narrower light curves, in conflict with the recent results of Thomas et al. Based on nebular spectra of 27 SN Ia, we find no relation between the FWHM of the iron emission feature at ~4700 A and Dm15(B) after removing the two low-luminosity SN 1986G and SN 1991bg, suggesting that the peak luminosity is not strongly dependent on the kinetic energy of the explosion for most SN Ia. Finally, we confirm the correlation of velocity shifts in some nebular lines with the intrinsic B-V color of SN Ia at maximum light, although several outliers suggest a possible non-monotonic behavior for the largest blueshifts.

We compare models for Type Ia supernova (SN Ia) light curves and spectra with an extensive set of observations. The models come from a recent survey of 44 two-dimensional delayed-detonation models computed by Kasen, Roepke & Woosley (2009), each viewed from multiple directions. The data include optical light curves of 251 SNe Ia and 2231 low-dispersion spectra from the Center for Astrophysics, plus data from the literature. The analysis uses standard techniques employed by observers, including MLCS2k2, SALT2, and SNooPy for light-curve analysis, and the Supernova Identification (SNID) code of Blondin & Tonry for spectroscopic comparisons to assess how well the models match the data. We show that the models that match observed spectra best lie systematically on the observed width-luminosity relation. Conversely, we reject six models with highly asymmetric ignition conditions and a large amount (>1 M_sun) of synthesized 56Ni that yield poor matches to observed SN Ia spectra. More subtle features of the comparison include the general difficulty of the models to match the U-band flux at early times, caused by a hot ionized ejecta that affect the subsequent redistribution of flux at longer wavelengths. We examine ways in which the asymptotic kinetic energy of the explosion affects both the predicted velocity and velocity gradient in the Si II and Ca II lines. Models with an asymmetric distribution of 56Ni are found to result in a larger variation of photometric and spectroscopic properties with viewing angle, regardless of the initial ignition setup. We discuss more generally whether highly anisotropic ignition conditions are ruled out by observations, and how detailed comparisons between models and observations involving both light curves and spectra can lead to a better understanding of SN Ia explosion mechanisms.

We report the first detection of asymmetry in a supernova (SN) photosphere based on SN light echo (LE) spectra of Cas A from the different perspectives of dust concentrations on its LE ellipsoid. New LEs are reported based on difference images, and optical spectra of these LEs are analyzed and compared. After properly accounting for the effects of finite dust-filament extent and inclination, we find one field where the He I and H alpha features are blueshifted by an additional ~4000 km/s relative to other spectra and to the spectra of the Type IIb SN 1993J. That same direction does not show any shift relative to other Cas A LE spectra in the Ca II near-infrared triplet feature. We compare the perspectives of the Cas A LE dust concentrations with recent three-dimensional modeling of the SN remnant (SNR) and note that the location having the blueshifted He I and H alpha features is roughly in the direction of an Fe-rich outflow and in the opposite direction of the motion of the compact object at the center of the SNR. We conclude that Cas A was an intrinsically asymmetric SN. Future LE spectroscopy of this object, and of other historical SNe, will provide additional insight into the connection of explosion mechanism to SN to SNR, as well as give crucial observational evidence regarding how stars explode.

The light echo systems of historical supernovae in the Milky Way and local group galaxies provide an unprecedented opportunity to reveal the effects of asymmetry on observables, particularly optical spectra. Scattering dust at different locations on the light echo ellipsoid witnesses the supernova from different perspectives and the light consequently scattered towards Earth preserves the shape of line profile variations introduced by asymmetries in the supernova photosphere. However, the interpretation of supernova light echo spectra to date has not involved a detailed consideration of the effects of outburst duration and geometrical scattering modifications due to finite scattering dust filament dimension, inclination, and image point-spread function and spectrograph slit width. In this paper, we explore the implications of these factors and present a framework for future resolved supernova light echo spectra interpretation, and test it against Cas A and SN 1987A light echo spectra. We conclude that the full modeling of the dimensions and orientation of the scattering dust using the observed light echoes at two or more epochs is critical for the correct interpretation of light echo spectra. Indeed, without doing so one might falsely conclude that differences exist when none are actually present.

[Abridged] We investigate the use of a wide variety of spectroscopic measurements to determine distances to low-redshift Type Ia supernovae (SN Ia). We consider linear models for predicting distances to SN Ia using light-curve width and color parameters (determined using the SALT2 light-curve fitter) and a spectroscopic indicator, and evaluate the resulting Hubble diagram scatter using a cross-validation procedure. We confirm the ability of spectral flux ratios alone at maximum light to reduce the scatter of Hubble residuals by ~10% with respect to the standard combination of light-curve width and color. When used in combination with the SALT2 color parameter, the color-corrected flux ratio R^c(6420/5290) at maximum light leads to an even lower scatter, although the improvement has low statistical significance (<2 sigma) given the size of our sample (26 SN Ia). We highlight the importance of an accurate relative flux calibration and the failure of this method for highly-reddened objects. Comparison with synthetic spectra from 2D delayed-detonation explosion models shows that the correlation of R(6630/4400) with SN Ia absolute magnitudes can be largely attributed to intrinsic color variations and not to reddening by dust in the host galaxy. We consider flux ratios at other ages, as well as the use of pairs of flux ratios, revealing the presence of small-scale intrinsic spectroscopic variations in the iron-group dominated absorption features around ~4300 A and ~4800 A. The best flux ratio overall is the color-corrected R^c(4610/4260) at t=-2.5d from maximum light, which leads to ~30% lower scatter with respect to the standard combination of light-curve width and color. We examine other spectroscopic indicators related to line-profile morphology, but none appear to lead to a significant improvement over the standard light-curve width and color parameters.

Recent observations have revealed that some Type Ia supernovae exhibit narrow, time-variable Na I D absorption features. The origin of the absorbing material is controversial, but it may suggest the presence of circumstellar gas in the progenitor system prior to the explosion, with significant implications for the nature of the supernova progenitors. We present the third detection of such variable absorption, based on six epochs of high-resolution spectroscopy of the Type Ia supernova SN 2007le from Keck and the HET. The data span ~3 months, from 5 days before maximum light to 90 days after maximum. We find that one component of the Na D absorption lines strengthened significantly with time, indicating a total column density increase of ~2.5 x 10^12 cm^-2. The changes are most prominent after maximum light rather than at earlier times when the UV flux from the SN peaks. As with SN 2006X, we detect no change in the Ca II H&K lines over the same time period, rendering line-of-sight effects improbable and suggesting a circumstellar origin for the absorbing material. Unlike the previous two SNe exhibiting variable absorption, SN 2007le is not highly reddened (E_B-V = 0.27 mag), also pointing toward circumstellar rather than interstellar absorption. Photoionization models show that the data are consistent with a dense (10^7 cm^-3) cloud or clouds of gas located ~0.1 pc from the explosion. These results broadly support the single-degenerate scenario previously proposed to explain the variable absorption, with mass loss from a nondegenerate companion star responsible for providing the circumstellar gas. We also present tentative evidence for narrow Halpha emission associated with the SN, which will require followup observations at late times to confirm. [abridged]

We present extensive early photometric (ultraviolet through near-infrared) and spectroscopic (optical and near-infrared) data on supernova (SN) 2008D as well as X-ray data analysis on the associated Swift/X-ray transient (XRT) 080109. Our data span a time range of 5 hours before the detection of the X-ray transient to 150 days after its detection, and detailed analysis allowed us to derive constraints on the nature of the SN and its progenitor; throughout we draw comparisons with results presented in the literature and find several key aspects that differ. We show that the X-ray spectrum of XRT 080109 can be fit equally well by an absorbed power law or a superposition of about equal parts of both power law and blackbody. Our data first established that SN 2008D is a spectroscopically normal SN Ib (i.e., showing conspicuous He lines), and show that SN 2008D had a relatively long rise time of 18 days and a modest optical peak luminosity. The early-time light curves of the SN are dominated by a cooling stellar envelope (for \Delta t~0.1- 4 day, most pronounced in the blue bands) followed by 56^Ni decay. We construct a reliable measurement of the bolometric output for this stripped-envelope SN, and, combined with estimates of E_K and M_ej from the literature, estimate the stellar radius R_star of its probable Wolf-Rayet progenitor. According to the model of Waxman et al. and of Chevalier & Fransson, we derive R_star^{W07}= 1.2+/-0.7 R_sun and R_star^{CF08}= 12+/-7 R_sun, respectively; the latter being more in line with typical WN stars. Spectra obtained at 3 and 4 months after maximum light show double-peaked oxygen lines that we associate with departures from spherical symmetry, as has been suggested for the inner ejecta of a number of SN Ib cores.

We combine the CfA3 supernova Type Ia (SN Ia) sample with samples from the literature to calculate improved constraints on the dark energy equation of state parameter, w. The CfA3 sample is added to the Union set of Kowalski et al. (2008) to form the Constitution set and, combined with a BAO prior, produces 1+w=0.013 +0.066/-0.068 (0.11 syst), consistent with the cosmological constant. The CfA3 addition makes the cosmologically-useful sample of nearby SN Ia between 2.6 and 2.9 times larger than before, reducing the statistical uncertainty to the point where systematics play the largest role. We use four light curve fitters to test for systematic differences: SALT, SALT2, MLCS2k2 (R_V=3.1), and MLCS2k2 (R_V=1.7). SALT produces high-redshift Hubble residuals with systematic trends versus color and larger scatter than MLCS2k2. MLCS2k2 overestimates the intrinsic luminosity of SN Ia with 0.7 < Delta < 1.2. MLCS2k2 with R_V=3.1 overestimates host-galaxy extinction while R_V=1.7 does not. Our investigation is consistent with no Hubble bubble. We also find that, after light-curve correction, SN Ia in Scd/Sd/Irr hosts are intrinsically fainter than those in E/S0 hosts by 2 sigma, suggesting that they may come from different populations. We also find that SN Ia in Scd/Sd/Irr hosts have low scatter (0.1 mag) and reddening. Current systematic errors can be reduced by improving SN Ia photometric accuracy, by including the CfA3 sample to retrain light-curve fitters, by combining optical SN Ia photometry with near-infrared photometry to understand host-galaxy extinction, and by determining if different environments give rise to different intrinsic SN Ia luminosity after correction for light-curve shape and color.

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

- Constraining the cosmological parameters and understanding Dark Energy have tremendous implications for the nature of the Universe and its physical laws. - The pervasive limit of systematic uncertainties reached by cosmography based on Cepheids and Type Ia supernovae (SNe Ia) warrants a search for complementary approaches. - Type II SNe have been shown to offer such a path. Their distances can be well constrained by luminosity-based or geometric methods. Competing, complementary, and concerted efforts are underway, to explore and exploit those objects that are extremely well matched to next generation facilities. Spectroscopic follow-up will be enabled by space- based and 20-40 meter class telescopes. - Some systematic uncertainties of Type II SNe, such as reddening by dust and metallicity effects, are bound to be different from those of SNe Ia. Their stellar progenitors are known, promising better leverage on cosmic evolution. In addition, their rate - which closely tracks the ongoing star formation rate - is expected to rise significantly with look- back time, ensuring an adequate supply of distant examples. - These data will competitively constrain the dark energy equation of state, allow the determination of the Hubble constant to 5%, and promote our understanding of the processes involved in the last dramatic phases of massive stellar evolution.

We have obtained 1087 NIR (JHKs) measurements of 21 SNe Ia using PAIRITEL, nearly doubling the number of well-sampled NIR SN Ia light curves. These data strengthen the evidence that SNe Ia are excellent standard candles in the NIR, even without correction for optical light-curve shape. We construct fiducial NIR templates for normal SNe Ia from our sample, excluding only the three known peculiar SNe Ia: SN 2005bl, SN 2005hk, and SN 2005ke. The H-band absolute magnitudes in this sample of 18 SNe Ia have an intrinsic rms of only 0.15 mag with no correction for light-curve shape. We found a relationship between the H-band extinction and optical color excess of AH=0.2E(B-V). This variation is as small as the scatter in distance modulus measurements currently used for cosmology based on optical light curves after corrections for light-curve shape. Combining the homogeneous PAIRITEL measurements with 23 SNe Ia from the literature, these 41 SNe Ia have standard H-band magnitudes with an rms scatter of 0.16 mag. The good match of our sample with the literature sample suggests there are few systematic problems with the photometry. We present a nearby NIR Hubble diagram that shows no correlation of the residuals from the Hubble line with light-curve properties. Future samples that account for optical and NIR light-curve shapes, absorption, spectroscopic variation, or host-galaxy properties may reveal effective ways to improve the use of SNe Ia as distance indicators. Since systematic errors due to dust absorption in optical bands remain the leading difficulty in the cosmological use of supernovae, the good behavior of SN Ia NIR light curves and their relative insensitivity to reddening make these objects attractive candidates for future cosmological work.

We present the results of spectroscopic observations from the ESSENCE high-redshift supernova (SN) survey during its first four years of operation. This sample includes spectra of all SNe Ia whose light curves were presented by Miknaitis et al. (2007) and used in the cosmological analyses of Davis et al. (2007) and Wood-Vasey et al. (2007). The sample represents 273 hours of spectroscopic observations with 6.5 - 10-m-class telescopes of objects detected and selected for spectroscopy by the ESSENCE team. We present 174 spectra of 156 objects. Combining this sample with that of Matheson et al. (2005), we have a total sample of 329 spectra of 274 objects. From this, we are able to spectroscopically classify 118 Type Ia SNe. As the survey has matured, the efficiency of classifying SNe Ia has remained constant while we have observed both higher-redshift SNe Ia and SNe Ia farther from maximum brightness. Examining the subsample of SNe Ia with host-galaxy redshifts shows that redshifts derived from only the SN Ia spectra are consistent with redshifts found from host-galaxy spectra. Moreover, the phases derived from only the SN Ia spectra are consistent with those derived from light-curve fits. By comparing our spectra to local templates, we find that the rate of objects similar to the overluminous SN 1991T and the underluminous SN 1991bg in our sample are consistent with that of the local sample. We do note, however, that we detect no object spectroscopically or photometrically similar to SN 1991bg. Although systematic effects could reduce the high-redshift rate we expect based on the low-redshift surveys, it is possible that SN 1991bg-like SNe Ia are less prevalent at high redshift.

The spectral properties of type Ia supernovae in the ultraviolet (UV) are investigated using the early-time spectra of SN 2001ep and SN 2001eh obtained using the Hubble Space Telescope (HST). A series of spectral models is computed with a Monte Carlo spectral synthesis code, and the dependence of the UV flux on the elemental abundances and the density gradient in the outer layers of the ejecta is tested. A large fraction of the UV flux is formed by reverse fluorescence scattering of photons from red to blue wavelengths. This process, combined with ionization shifts due to enhanced line blocking, can lead to a stronger UV flux as the iron-group abundance in the outer layers is increased, contrary to previous claims.

Double-peaked oxygen lines in the nebular spectra of two peculiar Type Ib/c Supernovae (SN Ib/c) have been interpreted as off-axis views of a GRB-jet or unipolar blob ejections. Here we present late-time spectra of 8 SN IIb, Ib and Ic and show that this phenomenon is common and should not be so firmly linked to extraordinary explosion physics. The line profiles are most likely caused by ejecta expanding with a torus- or disk-like geometry. Double-peaked oxygen profiles are not necessarily the indicator of a mis-directed GRB jet.

We report the discovery and orbit determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging dataset. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the Solar System ecliptic (-21 deg < beta < -5 deg), reaching limiting magnitudes per observation of I~23.1 and R~23.7. We examine several of the newly detected objects in detail, including 2003 UC_414 which orbits entirely between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the Solar System. 2003 SS_422 and 2007 TA_418 have high eccentricities and large perihelia, making them candidate members of an outer class of trans-Neptunian objects. We also report a new member of the ''extended'' or ''detached'' scattered disk, 2004 VN_112, and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ~2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early Solar System.

We present multiepoch spectra of 13 high-redshift Type Ia supernovae (SNe Ia) drawn from the literature, the ESSENCE and SNLS projects, and our own separate dedicated program on the ESO Very Large Telescope. We use the Supernova Identification (SNID) code of Blondin & Tonry to determine the spectral ages in the supernova rest frame. Comparison with the observed elapsed time yields an apparent aging rate consistent with the 1/(1+z) factor (where z is the redshift) expected in a homogeneous, isotropic, expanding universe. These measurements thus confirm the expansion hypothesis, while unambiguously excluding models that predict no time dilation, such as Zwicky's "tired light" hypothesis. We also test for power-law dependencies of the aging rate on redshift. The best-fit exponent for these models is consistent with the expected 1/(1+z) factor.

We present 432 low-dispersion optical spectra of 32 Type Ia supernovae (SNe Ia) that also have well-calibrated light curves. The coverage ranges from 6 epochs to 36 epochs of spectroscopy. Most of the data were obtained with the 1.5m Tillinghast telescope at the F. L. Whipple Observatory with typical wavelength coverage of 3700-7400A and a resolution of ~7A. The earliest spectra are thirteen days before B-band maximum; two-thirds of the SNe were observed before maximum brightness. Coverage for some SNe continues almost to the nebular phase. The consistency of the method of observation and the technique of reduction makes this an ideal data set for studying the spectroscopic diversity of SNe Ia.

We report the successful identification of the type of the supernova responsible for the supernova remnant SNR 0509-675 in the Large Magellanic Cloud (LMC) using Gemini spectra of surrounding light echoes. The ability to classify outbursts associated with centuries-old remnants provides a new window into several aspects of supernova research and is likely to be successful in providing new constraints on additional LMC supernovae as well as their historical counterparts in the Milky Way Galaxy (MWG). The combined spectrum of echo light from SNR 0509-675 shows broad emission and absorption lines consistent with a supernova (SN) spectrum. We create a spectral library consisting of 26 SNe Ia and 6 SN Ib/c that are time-integrated, dust-scattered by LMC dust, and reddened by the LMC and MWG. We fit these SN templates to the observed light echo spectrum using $\chi^2$ minimization as well as correlation techniques, and we find that overluminous 91T-like SNe Ia with $\dm15<0.9$ match the observed spectrum best.

We analyze the Type II Plateau supernovae (SN II-P) 2005cs and 2006bp with the non-LTE model atmosphere code CMFGEN. We fit 13 spectra in the first month for SN 2005cs and 18 for SN 2006bp. {\sl Swift} ultraviolet photometry and ground-based optical photometry calibrate each spectrum. Our analysis shows both objects were discovered less than 3 days after they exploded, making these the earliest SN II-P spectra ever studied. They reveal broad and very weak lines from highly-ionized fast ejecta with an extremely steep density profile. We identify He{\sc ii} 4686\AA emission in the SN 2006bp ejecta. Days later, the spectra resemble the prototypical Type II-P SN 1999em, which had a supergiant-like photospheric composition. Despite the association of SN 2005cs with possible X-ray emission, the emergent UV and optical light comes from the photosphere, not from circumstellar emission. We surmise that the very steep density fall-off we infer at early times may be a fossil of the combined actions of the shock wave passage and radiation driving at shock breakout. Based on tailored CMFGEN models, the direct-fitting technique and the Expanding Photosphere Method both yield distances and explosion times that agree within a few percent. We derive a distance to NGC 5194, the host of SN 2005cs, of 8.9$\pm$0.5 Mpc and 17.5$\pm$0.8 Mpc for SN 2006bp in NGC 3953. The luminosity of SN 2006bp is 1.5 times that of SN 1999em, and 6 times that of SN 2005cs. Reliable distances to Type II-P supernovae that do not depend on a small range in luminosity provide an independent route to the Hubble Constant and improved constraints on other cosmological parameters.

We present the first large-scale effort of creating composite spectra of high-redshift type Ia supernovae (SNe Ia) and comparing them to low-redshift counterparts. Through the ESSENCE project, we have obtained 107 spectra of 88 high-redshift SNe Ia with excellent light-curve information. In addition, we have obtained 397 spectra of low-redshift SNe through a multiple-decade effort at Lick and Keck Observatories, and we have used 45 UV spectra obtained by HST/IUE. The low-redshift spectra act as a control sample when comparing to the ESSENCE spectra. In all instances, the ESSENCE and Lick composite spectra appear very similar. The addition of galaxy light to the Lick composite spectra allows a nearly perfect match of the overall spectral-energy distribution with the ESSENCE composite spectra, indicating that the high-redshift SNe are more contaminated with host-galaxy light than their low-redshift counterparts. This is caused by observing objects at all redshifts with the same slit width, which corresponds to different projected distances. After correcting for the galaxy-light contamination, subtle differences in the spectra remain. We have estimated the systematic errors when using current spectral templates for K-corrections to be ~0.02 mag. The variance in the composite spectra give an estimate of the intrinsic variance in low-redshift maximum-light SN spectra of ~3% in the optical and growing toward the UV. The difference between the maximum light low and high-redshift spectra constrain SN evolution between our samples to be < 10% in the rest-frame optical.

We present an algorithm to identify the type of an SN spectrum and to determine its redshift and age. This algorithm, based on the correlation techniques of Tonry & Davis, is implemented in the Supernova Identification (SNID) code. It is used by members of ongoing high-redshift SN searches to distinguish between type Ia and type Ib/c SNe, and to identify "peculiar" SNe Ia. We develop a diagnostic to quantify the quality of a correlation between the input and template spectra, which enables a formal evaluation of the associated redshift error. Furthermore, by comparing the correlation redshifts obtained using SNID with those determined from narrow lines in the SN host galaxy spectrum, we show that accurate redshifts (with a typical error less than 0.01) can be determined for SNe Ia without a spectrum of the host galaxy. Last, the age of an input spectrum is determined with a typical 3-day accuracy, shown here by using high-redshift SNe Ia with well-sampled light curves. The success of the correlation technique confirms the similarity of some SNe Ia at low and high redshifts. The SNID code, which is available to the community, can also be used for comparative studies of SN spectra, as well as comparisons between data and models.

Spectra and light curves of SN 2006gz show the strongest signature of unburned carbon and one of the slowest fading light curves ever seen in a type Ia event (Delta m_15 = 0.69 +/- 0.04). The early-time Si II velocity is low, implying it was slowed by an envelope of unburned material. Our best estimate of the luminosity implies M_V = -19.74 and the production of ~ 1.2 M_sun of 56Ni. This suggests a super-Chandrasekhar mass progenitor. A double degenerate merger is consistent with these observations.

Motivated by the recognition that variation in the optical transmission of the atmosphere is probably the main limitation to the precision of ground-based CCD measurements of celestial fluxes, we review the physical processes that attenuate the passage of light through the Earth's atmosphere. The next generation of astronomical surveys, such as PanSTARRS and LSST, will greatly benefit from dedicated apparatus to obtain atmospheric transmission data that can be associated with each survey image. We review and compare various approaches to this measurement problem, including photometry, spectroscopy, and LIDAR. In conjunction with careful measurements of instrumental throughput, atmospheric transmission measurements should allow next-generation imaging surveys to produce photometry of unprecedented precision. Our primary concerns are the real-time determination of aerosol scattering and absorption by water along the line of sight, both of which can vary over the course of a night's observations.

The first cosmological results from the ESSENCE supernova survey (Wood-Vasey et al. 2007) are extended to a wider range of cosmological models including dynamical dark energy and non-standard cosmological models. We fold in a greater number of external data sets such as the recent Higher-z release of high-redshift supernovae (Riess et al. 2007) as well as several complementary cosmological probes. Model comparison statistics such as the Bayesian and Akaike information criteria are applied to gauge the worth of models. These statistics favor models that give a good fit with fewer parameters. Based on this analysis, the preferred cosmological model is the flat cosmological constant model, where the expansion history of the universe can be adequately described with only one free parameter describing the energy content of the universe. Among the more exotic models that provide good fits to the data, we note a preference for models whose best-fit parameters reduce them to the cosmological constant model.

We describe the implementation and optimization of the ESSENCE supernova survey, which we have undertaken to measure the equation of state parameter of the dark energy. We present a method for optimizing the survey exposure times and cadence to maximize our sensitivity to the dark energy equation of state parameter w=P/rho c^2 for a given fixed amount of telescope time. For our survey on the CTIO 4m telescope, measuring the luminosity distances and redshifts for supernovae at modest redshifts (z~0.5 +- 0.2) is optimal for determining w. We describe the data analysis pipeline based on using reliable and robust image subtraction to find supernovae automatically and in near real-time. Since making cosmological inferences with supernovae relies crucially on accurate measurement of their brightnesses, we describe our efforts to establish a thorough calibration of the CTIO 4m natural photometric system. In its first four years, ESSENCE has discovered and spectroscopically confirmed 102 type Ia SNe, at redshifts from 0.10 to 0.78, identified through an impartial, effective methodology for spectroscopic classification and redshift determination. We present the resulting light curves for the all type Ia supernovae found by ESSENCE and used in our measurement of w, presented in Wood-Vasey et al, 2007.

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

We report early photospheric-phase observations of the Type IIP Supernova (SN) 2005cs obtained by Swift's Ultraviolet-Optical and X-Ray Telescopes. Observations started within two days of discovery and continued regularly for three weeks. During this time the V-band magnitude remained essentially constant, while the UV was initially bright but steadily faded until below the brightness of an underlying UV-bright HII region. This UV decay is similar to SNe II observed by the International Ultraviolet Explorer. UV grism spectra show the P-Cygni absorption of MgII 2798A, indicating a photospheric origin of the UV flux. Based on non-LTE model atmosphere calculations with the CMFGEN code, we associate the rapid evolution of the UV flux with the cooling of the ejecta, the peak of the spectral energy distribution (SED) shifting from ~700A on June 30th to ~1200A on July 5th. Furthermore, the corresponding recombination of the ejecta, e.g., the transition from FeIII to FeII, induces considerable strengthening of metal line-blanketing at and above the photosphere, blocking more effectively this fading UV flux. SN2005cs was not detected in X-rays, and the upper limit to the X-ray luminosity yields a limit to the mass loss rate of the progenitor of about 10^-5 solar masses per year. Overall, Swift presents a unique opportunity to capture the early and fast evolution of Type II SNe in the UV, providing additional constraints on the reddening, the SED shortward of 4000A, and the ionization state and temperature of the photon-decoupling regions.

We present an algorithm to identify the types of supernova spectra, and determine their redshift and phase. This algorithm, based on the correlation techniques of Tonry & Davis, is implemented in the SuperNova IDentification code (SNID). It is used by members of the ESSENCE project to determine whether a noisy spectrum of a high-redshift supernova is indeed of type Ia, as opposed to, e.g., type Ib/c. Furthermore, by comparing the correlation redshifts obtained using SNID with those determined from narrow lines in the supernova host galaxy spectrum, we show that accurate redshifts (with a typical error < 0.01) can be determined for SNe Ia for which a spectrum of the host galaxy is unavailable. Last, the phase of an input spectrum is determined with a typical accuracy of ~3 days.

We present a set of 11 type Ia supernova (SN Ia) lightcurves with dense, pre-maximum sampling. These supernovae (SNe), in galaxies behind the Large Magellanic Cloud (LMC), were discovered by the SuperMACHO survey. The SNe span a redshift range of z = 0.11 - 0.35. Our lightcurves contain some of the earliest pre-maximum observations of SNe Ia to date. We also give a functional model that describes the SN Ia lightcurve shape (in our VR-band). Our function uses the "expanding fireball" model of Goldhaber et al. (1998) to describe the rising lightcurve immediately after explosion but constrains it to smoothly join the remainder of the lightcurve. We fit this model to a composite observed VR-band lightcurve of three SNe between redshifts of 0.135 to 0.165. These SNe have not been K-corrected or adjusted to account for reddening. In this redshift range, the observed VR-band most closely matches the rest frame V-band. Using the best fit to our functional description of the lightcurve, we find the time between explosion and observed VR-band maximum to be 17.6+-1.3(stat)+-0.07(sys) rest-frame days for a SN Ia with a VR-band Delta m_{-10} of 0.52mag. For the redshifts sampled, the observed VR-band time-of-maximum brightness should be the same as the rest-frame V-band maximum to within 1.1 rest-frame days.

We present early photometric and spectroscopic data on the afterglow of GRB 060218 and report the evolution of the underlying supernova 2006aj. Our data span a time-range of 4 days to 22 days after the GRB and clearly establish that SN 2006aj is a fast-evolving broad-lined Type Ic SN with an extremely short rise-time (~ 10 days) and a large optical luminosity (M_V = -18.7 mag). The SN properties are deduced well since the GRB afterglow does not contribute a significant amount to the total light output. The spectra show broad lines indicative of large expansion velocities, but are better matched by those of SN 2002ap and SN 1997ef (that are not associated with a GRB) than those of the proto-typical GRB-related SN 1998bw. We refine the redshift estimate to z = 0.0335 +/- 0.00007. The host-galaxy is a low-metallicity dwarf galaxy (with M_V \~ -16.0 mag), similar to host-galaxies of other GRB-associated SNe.

Using archival data of low-redshift (z < 0.01) Type Ia supernovae (SN Ia) and recent observations of high-redshift (0.16 < z <0.64; Matheson et al. 2005) SN Ia, we study the "uniformity'' of the spectroscopic properties of nearby and distant SN Ia. We find no difference in the measures we describe here. In this paper, we base our analysis solely on line-profile morphology, focusing on measurements of the velocity location of maximum absorption (vabs) and peak emission (vpeak). We find that the evolution of vabs and vpeak for our sample lines (Ca II 3945, Si II 6355, and S II 5454, 5640) is similar for both the low- and high-redshift samples. We find that vabs for the weak S II 5454, 5640 lines, and vpeak for S II 5454, can be used to identify fast-declining [dm15 > 1.7] SN Ia, which are also subluminous. In addition, we give the first direct evidence in two high-z SN Ia spectra of a double-absorption feature in Ca II 3945, an event also observed, though infrequently, in low-redshift SN Ia spectra (6/22 SN Ia in our local sample). We report for the first time the unambiguous and systematic intrinsic blueshift of peak emission of optical P-Cygni line profiles in Type Ia spectra, by as much as 8000 km/s. All the high-z SN Ia analyzed in this paper were discovered and followed up by the ESSENCE collaboration, and are now publicly available.

The proper usage of Type Ia supernovae (SNe Ia) as distance indicators has revolutionized cosmology, and added a new dominant component to the energy density of the Universe, dark energy. Following the discovery and confirmation era, the currently ongoing SNe Ia surveys aim to determine the properties of the dark energy. ESSENCE is a five year ground-based supernova survey aimed at finding and characterizing 200 SNe Ia in the redshift domain z=[0.2-0.8]. The goal of the project is to put constraints on the equation of state parameter, w, of the dark energy with an accuracy of <10%. This paper presents these ongoing efforts in the context of the current developments in observational cosmology.

We present broad-band light curves of nine supernovae ranging in redshift from 0.5 to 0.8. The supernovae were discovered as part of the ESSENCE project, and the light curves are a combination of Cerro Tololo 4-m and Hubble Space Telescope (HST) photometry. On the basis of spectra and/or light-curve fitting, eight of these objects are definitely Type Ia supernovae, while the classification of one is problematic. The ESSENCE project is a five-year endeavor to discover about 200 high-redshift Type Ia supernovae, with the goal of tightly constraining the time average of the equation-of-state parameter [w = p/(rho c^2)] of the "dark energy." To help minimize our systematic errors, all of our ground-based photometry is obtained with the same telescope and instrument. In 2003 the highest-redshift subset of ESSENCE supernovae was selected for detailed study with HST. Here we present the first photometric results of the survey. We find that all but one of the ESSENCE SNe have slowly declining light curves, and the sample is not representative of the low-redshift set of ESSENCE Type Ia supernovae. This is unlikely to be a sign of evolution in the population. We attribute the decline-rate distribution of HST events to a selection bias at the high-redshift edge of our sample and find that such a bias will infect other magnitude-limited SN Ia searches unless appropriate precautions are taken.

We present the results of spectroscopic observations of targets discovered during the first two years of the ESSENCE project. The goal of ESSENCE is to use a sample of ~200 Type Ia supernovae (SNe Ia) at moderate redshifts (0.2 < z < 0.8) to place constraints on the equation of state of the Universe. Spectroscopy not only provides the redshifts of the objects, but also confirms that some of the discoveries are indeed SNe Ia. This confirmation is critical to the project, as techniques developed to determine luminosity distances to SNe Ia depend upon the knowledge that the objects at high redshift are the same as the ones at low redshift. We describe the methods of target selection and prioritization, the telescopes and detectors, and the software used to identify objects. The redshifts deduced from spectral matching of high-redshift SNe Ia with low-redshift SNe Ia are consistent with those determined from host-galaxy spectra. We show that the high-redshift SNe Ia match well with low-redshift templates. We include all spectra obtained by the ESSENCE project, including 52 SNe Ia, 5 core-collapse SNe, 12 active galactic nuclei, 19 galaxies, 4 possibly variable stars, and 16 objects with uncertain identifications.

We use a new technique to extract the spectrum of a supernova from that of the contaminating background of its host galaxy, and apply it to the specific case of high-redshift Type Ia supernova (SN Ia) spectroscopy. The algorithm is based on a two-channel iterative technique employing the Richardson-Lucy restoration method and is implemented in the IRAF code 'specinholucy'. We run the code both on simulated (SN Ia at z=0.5 embedded in a bright host galaxy) and observed (SNe Ia at various phases up to z=0.236) data taken with VLT+FORS1 and show the advantages of using such a deconvolution technique in comparison with less elaborate methods. This paper is motivated by the need for optimal supernova spectroscopic data reduction in order to make meaningful comparisons between the low and high-redshift SN Ia samples. This may reveal subtle evolutionary and systematic effects that could depend on redshift and bias the cosmological results derived from comparisons of local and high-z SNe Ia in recent years. We describe the various aspects of the extraction in some detail as guidelines for the first-time user and present an optimal observing strategy for successful implementation of this method in future high-z SN Ia spectroscopic follow-up programmes.

We present results from the Hubble Higher-z Supernova Search, the first space-based open field survey for supernovae (SNe). In cooperation with the Great Observatories Origins Deep Survey, we have used the Hubble Space Telescope with the Advanced Camera for Surveys to cover 300 square arcmin in the area of the Chandra Deep Field South and the Hubble Deep Field North on five separate search epochs (separated by 45 day intervals) to a limiting magnitude of z'=26. These deep observations have allowed us to discover 42 SNe in the redshift range 0.2 < z < 1.6. As these data span a large range in redshift, they are ideal for testing the validity of Type Ia supernova progenitor models with the distribution of expected ``delay times,'' from progenitor star formation to SN Ia explosion, and the SN rates these models predict. Through a Bayesian maximum likelihood test, we determine which delay-time models best reproduce the redshift distribution of SNe Ia discovered in this survey. We find that models that require a large fraction of ``prompt'' (less than 2 Gyr) SNe Ia poorly reproduce the observed redshift distribution and are rejected at 95% confidence. We find that Gaussian models best fit the observed data for mean delay times in the range of 3 to 4 Gyr.

We present photometric and spectroscopic observations of 23 high redshift supernovae spanning a range of z=0.34-1.03, 9 of which are unambiguously classified as Type Ia. These supernovae were discovered during the IfA Deep Survey, which began in September 2001 and observed a total of 2.5 square degrees to a depth of approximately m=25-26 in RIZ over 9-17 visits, typically every 1-3 weeks for nearly 5 months, with additional observations continuing until April 2002. We give a brief description of the survey motivations, observational strategy, and reduction process. This sample of 23 high-redshift supernovae includes 15 at z>0.7, doubling the published number of objects at these redshifts, and indicates that the evidence for acceleration of the universe is not due to a systematic effect proportional to redshift. In combination with the recent compilation of Tonry et al. (2003), we calculate cosmological parameter density contours which are consistent with the flat universe indicated by the CMB (Spergel et al. 2003). Adopting the constraint that Omega_total = 1.0, we obtain best-fit values of (Omega_m, Omega_Lambda)=(0.33, 0.67) using 22 SNe from this survey augmented by the literature compilation. We show that using the empty-beam model for gravitational lensing does not eliminate the need for Omega_Lambda > 0. Experience from this survey indicates great potential for similar large-scale surveys while also revealing the limitations of performing surveys for z>1 SNe from the ground.