A review of type Ia supernova spectra

SN 2011fe was the nearest and best-observed type Ia supernova in a generation, and brought previous incomplete datasets into sharp contrast with the detailed new data. In retrospect, documenting spectroscopic behaviors of type Ia supernovae has been more often limited by sparse and incomplete temporal sampling than by consequences of signal-to-noise ratios, telluric features, or small sample sizes. As a result, type Ia supernovae have been primarily studied insofar as parameters discretized by relative epochs and incomplete temporal snapshots near maximum light. Here we discuss a necessary next step toward consistently modeling and directly measuring spectroscopic observables of type Ia supernova spectra. In addition, we analyze current spectroscopic data in the parameter space defined by empirical metrics, which will be relevant even after progenitors are observed and detailed models are refined.


Introduction
The transient nature of extragalactic type Ia supernovae (SN Ia) prevent studies from conclusively singling out unobserved progenitor configurations (Roelofs et al. 2008;Li et al. 2011b;Kilic et al. 2013). It remains fairly certain that the progenitor system of SN Ia comprises at least one compact C+O white dwarf (Chandrasekhar 1957;Nugent et al. 2011;Bloom et al. 2012). However, how the state of this primary star reaches a critical point of disruption continues to elude astronomers. This is particularly so given that less than ∼ 15% of locally observed white dwarfs have a mass a few 0.1M greater than a solar mass; very few systems near the formal Chandrasekhar-mass limit 5 , M Ch ≈ 1.38 M (Vennes 1999;Liebert et al. 2005;Napiwotzki et al. 2005;Parthasarathy et al. 2007;Napiwotzki et al. 2007).
Thus far observational constraints of SN Ia have been inconclusive in distinguishing between the following three separate theoretical considerations about possible progenitor scenarios. Along side perturbations in the critical mass limit or masses of the progenitors, e.g., from rotational support (Mueller & Eriguchi 1985;Yoon & Langer 2005;Chen & Li 2009;Hachisu et al. 2012;Tornambé & Piersanti 2013) or variances of white dwarf (WD) populations (van Kerkwijk et al. 2010;Dan et al. 2013), the primary WD may reach the critical point by accretion of material from a low-mass, radially-confined secondary star (Whelan & Iben 1973;Nomoto & Sugimoto 1977;Hayden et al. 2010a;Bianco et al. 2011;Bloom et al. 2012;Hachisu et al. 2012;Wheeler 2012;Mazzali et al. 2013;Chen et al. 2013), and/or through one of several white dwarf merger scenarios with a close binary companion (Webbink 1984;Iben & Tutukov 1984;Paczynski 1985;Thompson 2011;Wang et al. 2013a;Pakmor et al. 2013). In addition, the presence (or absence) of circumstellar material may not solely rule out particular progenitor systems as now both single-and double-degenerate systems are consis-tent with having polluted environments prior to the explosion (Shen et al. 2013;Phillips et al. 2013).
Meanwhile, and within the context of a well-observed spectroscopically normal SN 2011fe, recent detailed models and spectrum synthesis along with SN Ia rates studies, a strong case for merging binaries as the progenitors of normal SN Ia has surfaced (c.f., van Kerkwijk et al. 2010;Li et al. 2011c;Blondin et al. 2012;Chomiuk 2013;Dan et al. 2013;Moll et al. 2013;Maoz et al. 2013;Johansson et al. 2014). However, because no progenitor system has ever been connected to any SN Ia, most observational constraints and trends are difficult to robustly impose on a standard model picture for even a single progenitor channel; the SN Ia problem is yet to be confined for each SN Ia subtype.
As for restricting SN Ia subtypes to candidate progenitor systems: (i) observed "jumps" between mean properties of SN Ia subtypes signify potential differences of progenitors and/or explosion mechanisms, (ii) the dispersions of individual subtypes are thought to arise from various abundance, density, metallicity, and/or temperature enhancements of the original progenitor system's post-explosion ejecta tomography, and (iii) "transitional-type" SN Ia complicate the already similar overlap of observed SN Ia properties (Nugent et al. 1995;Lentz et al. 2000;Benetti et al. 2005;Branch et al. 2009;Höflich et al. 2010;Wang et al. 2012Wang et al. , 2013cDessart et al. 2013a). Moreover, our physical understanding of all observed SN Ia subclasses remains based entirely on interpretations of idealized explosion models that are so far constrained and evaluated by "goodness of fit" comparisons to incomplete observations, particularly for SN Ia spectra at all epochs.
By default, spectra have been a limiting factor of supernova studies due to associated observational consequences, e.g., impromptu transient targets, variable intrinsic peak luminosities, a sparsity of complete datasets in wavelength and time, insufficient signalto-noise ratios, and the ever-present obstacle of spectroscopic line blending (Payne- Gaposchkin & Whipple 1940). Subsequently, two frequently relied upon empirical quantifiers of SN Ia spectroscopic diversity have been the rate at which rest-frame 6100Å absorption minima shift redward vis-à-vis projected Doppler velocities of the absorbing Si-rich material Wang et al. 2009a) and absorption strength measurements (a.k.a. pseudo equivalent widths; pEWs) of several lines of interest (see Hachinger et al. 2006;Silverman et al. 2012b;Blondin et al. 2012). Together these classification schemes moreor-less describe the same events by two interconnected parameter spaces (i.e. flux and expansion velocities, Branch et al. 2009;Foley & Kasen 2011;Blondin et al. 2012) that are dependent on a multi-dimensional array of physical properties. Naturally, the necessary next step for supernova studies alike is the development of prescriptions for the physical diagnosis of spectroscopic behaviors (see §2.2 and Kerzendorf & Sim 2014).
The fact that certain subsets of normal SN Ia constitute a near homogenous group of intrinsically bright events has led to their use as standardizable distance indicators (Kowal 1968;Elias et al. 1985a;Branch & Tammann 1992;Riess et al. 1999;Perlmutter et al. 1999;Schmidt 2004;Mandel et al. 2011;Maeda et al. 2011;Sullivan et al. 2011a;Hicken et al. 2012). However, this same attribute of homogeneity remains the greatest challenge in the individual study of SN Ia given that the time-evolving spectrum of a supernova is unique unto itself from the earliest to the latest epochs.
Because SN Ia are invaluable tools for both cosmology and understanding progenitor populations, a multitude of large scale surveys, searches, and observing campaigns 6 are continually being carried out with regularly improved precision. Subsequently, this build-up of competing resources has also resulted in an ever growing number of new and important discoveries, with less than complete information for each. In fact, with so many papers published each year on various aspects of SN Ia, it can be difficult to keep track of new results and important developments, including the validity of past and present theoretical explosion simulations and their related observational interpretations (see Maoz et al. 2013 for the latest).
Here we compile some of the discussions on spectroscopic properties of SN Ia from the past decade of published works. In §2 we overview the most common means for studying SN Ia: light curves ( §2.1), spectra ( §2.2), and detailed explosion models ( §2.3). In particular, we overview how far the well-observed SN 2011fe has progressed the degree of confidence associated with reading highly blended SN Ia spectra. Issues of SN Ia diversity are discussed in §3. Next, in §4 we recall several SN Ia that have made up the bulk of recent advances in uncovering the extent of their properties and peculiarities (see also the Appendix for a guide of some recent events). Finally, in §5 we summarize and conclude with some observational lessons of SN 2011fe.

Light curves
The interaction between the radiation field and the ejecta can be interpreted to zeroth order with the bolometric light curve. For SN Ia, the rise and fall of the light curve is said to be "powered" by 56 Ni produced in the explosion (Colgate & McKee 1969;Arnett 1982;Khokhlov et al. 1993;Mazzali et al. 1998;Pinto & Eastman 2000a;Stritzinger & Leibundgut 2005). Additional sources are expected to contribute to the overall luminosity behavior at various epochs 7 .  Timmes et al. 2003;Meng et al. 2011), interaction with circumstellar material (CSM, see Quimby et al. 2006b;Patat et al. 2007;Simon et al. 2007;Kasen 2010;Hayden et al. 2010a;Sternberg et al. 2011;Foley et al. 2012a;Förster et al. 2012;Shen et al. 2013;Silverman et al. 2013d; or an enshrouding C+O envelope Taubenberger et al. For example, Nomoto et al. (2003) has suggested that the variation of the carbon mass fraction in the C+O WD (C/O), or the variation of the initial WD mass, causes the diversity of SN Ia brightnesses (see Höflich et al. 2010). Similarly, Meng et al. (2011) argue that C/O and progenitor metallicity, Z, are intimately related for a fixed WD mass, and particularly for high metallicities given that it results in lower 3α burning rates plus an increased reduction of carbon via 12 C(α,γ) 16 O. For Z > Z (∼0.02), Meng et al. (2011) find that both C/O and Z have an approximately equal influence on 56 Ni production since, for a given WD mass, high progenitor metallicities (a greater abundance of species heavier than oxygen) and low C/O abundances (low carbon-rich fuel assuming a singledegenerate scenario) result in a low 56 Ni yield and subsequently dimmer SN Ia. For near solar metallicities or less, the carbon mass fraction plays a dominant role in 56 Ni production (Timmes et al. 2003). This then suggests that the average C/O ratio in the final state of the progenitor is an important physical cause, in addition to metallicity, for the observed width-luminosity relationship (WLR 8 ) of normal SN Ia light curves (Umeda et al. 1999a;Timmes et al. 2003;Nomoto et al. 2003;Bravo et al. 2010;Meng et al. 2011).
At the same time, the observed characteristics of SN Ia light curves and spectra can be fairly matched by adopting radial and/or axial shifts in the distribution of 56 Ni, possibly due to a delayed-and/or pulsational-detonation-like explosion mechanism (see Khokhlov 1991b;Hoflich et al. 1995;Baron et al. 2008;Bravo et al. 2009;Maeda et al. 2010b; Baron et al. 2012;Dessart et al. 2013a) or a merger scenario (e.g., Dan et al. 2013;Moll et al. 2013). Central ignition densities are also expected to play a secondary role in the form of the WLR since they are dependent upon the accretion rate of H and/or He-rich material and cooling time Höflich et al. 2010;Meng et al. 2010;Krueger et al. 2010;Sim et al. 2013), in addition to the spin-down timescales for differentially 2013), differences in total progenitor system masses (Hachisu et al. 2012;Pakmor et al. 2013;Chen et al. 2013), and directional dependent aspects of binary configurations (e.g., Blondin et al. 2011;Moll et al. 2013). 8 A WLR is followed when a SN Ia has a proportionately broader light curve for its intrinsic brightness at maximum light (Phillips 1993). Phillips et al. (1999) later extended this correlation by incorporating measurement of the extinction via late time B − V color measurements and B − V and V − I measurements at maximum light (see also Germany et al. 2004;Prieto et al. 2006b). Because lights curves of faint SN Ia evolve promptly before 15 days post-maximum light, light curve shape measurements are better suited for evaluating the light curve "stretch" (Conley et al. 2008).
rotating WDs (Hachisu et al. 2012;Tornambé & Piersanti 2013). Generally, discerning which of these factors dominate the spectrophotometric variation from one SN Ia to another remains a challenging task . As a result, astronomers are still mapping a broad range of SN Ia characteristics and trends ( §3).
Meanwhile, cosmological parameters determined by SN Ia light curves depend on an accurate comparison of nearby and distant events 9 . For distant and therefore redshifted SN Ia, a "K-correction" converts an observed magnitude to that which would be observed in the rest frame in another bandpass filter, allowing for the comparison of SN Ia brightnesses at various redshifts (Hogg et al. 2002). Consequently, K-corrections require the spectral energy distribution (SED) of the SN Ia and depend on SN Ia broad-band colors and the diversity of spectroscopic features . While some light curve fitters take a K-correction-less approach (e.g., Guy et al. 2005Guy et al. , 2007Conley et al. 2008), an SED is still required. A spectral template time series dataset is usually used since there exists remarkable homogeneity in the observed optical spectra of "normal" SN Ia (e.g., Hsiao et al. 2007).
Unfortunately there do remain poorly understood differences regarding spectroscopic feature strengths and inferred expansion velocities for these and other types of thermonuclear supernovae (see §2.2 and §3). At best, the spectroscopic diversity of SN Ia has been determined to be multidimensional Benetti et al. 2005;Branch et al. 2009;Wang et al. 2009a). Verily, SN Ia diversity studies require numerous large spectroscopic datasets in order to subvert many complex challenges faced when interpreting the data and extracting both projected Doppler velocities and "feature strength" measurements. However, studies that seek to primarily utilize SN Ia broad band luminosities need only collect a handful of sporadically sampled spectra in order to type the supernova event as a bona fide SN Ia. We note that interests in precision cosmology conflict at this point with the study of SN Ia. This is primarily because obtaining UBVRI photometry for hundreds of events is cheaper than collecting complete spectroscopy for a lesser number of SN Ia at various redshifts.
Nevertheless, the brightness decline rate in the Bband during the first 15 rest-frame days post-maximum light, ∆m 15 (B ), has proven useful for all SN Ia surveys. Phillips (1993) noted that ∆m 15 (B ) is well correlated with the intrinsic luminosity, a.k.a. the widthluminosity relationship. Previously, Khokhlov et al. (1993) did predict the existence of a WLR given that the light curve shape is sensitive to the time-dependent state of the ejected material.  recently utilized multidimensional time-dependent Monte Carlo radiative transfer calculations of Chandrasekhar-mass SN Ia models to access the physical relationship between the luminosity and light curve decline rate. They found that the WLR is largely a consequence of the radiative transfer inherent to SN Ia atmospheres, whereby the ionization evolution of iron redirects flux red ward and is hastened for dimmer and/or cooler SN Ia. Woosley et al. (2007) later explored the diversity of SN Ia light curves using a grid of 130 one-dimensional models. They concluded that a WLR is satisfied when SN Ia burn ∼ 1.1 M of material, with iron-group elements extending out to ∼ 8000 km s −1 .
Broadly speaking, the shape of the WLR is fundamentally influenced by the ionization evolution of iron group elements . However, since broad band luminosities are the sum of a supernova SED per wavelength interval, details of SN Ia diversity risk being "blurred out" for large samples of SN Ia. Therefore, decoding the spectra of all SN Ia subtypes, in addition to indirectly constraining detailed explosion models by the WLR, is of vital importance since variable signatures of iron-peak elements (IPEs) blend themselves within an SED typically populated by relatively strong features of overlapping signatures of intermediate-mass elements (IMEs).

Spectra
Supernova spectra detail information about the explosion and its local environment. To isolate and extract physical details (and determine their order of influence), several groups have invested greatly in advancing the computation of synthetic spectra for SN Ia, particularly during the early phases of homologous expansion (e.g., Hauschildt & Baron 1999;Kasen et al. 2002;Thomas et al. 2002;Höflich et al. 2002;Sauer et al. 2006;Kasen et al. 2006;Jeffery & Mazzali 2007;Sim et al. 2010a; Thomas et al. 2011a;Hoffmann et al. 2013;Pauldrach et al. 2013;Kerzendorf & Sim 2014). Although, even the basic facets of the supernova radiation environment serve as obstacles for timely computations of physically accurate, statistically representative, and robustly certain synthetic spectra (e.g., consequences of expansion).
It is the time-dependent interaction of the radiation field with the expanding material that complicates drawing conclusions about the explosion physics from the observations 10 . In a sense, there are two stages during which direct (and accessible) information about the progenitor system is driven away from being easily discernible within the post-explosion spectra: explosive nucleosynthesis and radiation transport 11 . That is to say, the ability to reproduce both the observed light curve and spectra, as well as the range of observed characteristics among SN Ia, is essential towards validating and/or restricting any explosion model for a given subtype.
Moreover, this assumes the sources of observed spectroscopic signatures in all varieties of SN Ia are known a priori, which is not necessarily the case given the immense volume of actively contributing atomic line transitions and continuum processes (Baron et al. , 1996Kasen et al. 2008;Bongard et al. 2008;Sauer et al. 2008). In fact, several features throughout the spectra have been either tentatively associated with a particular blend of atomic lines or identified with a multiple of conflicting suggestions (e.g., forbidden versus permitted lines at late or "nebular" transitional phases, see Bowers et al. 1997;Branch et al. 2005;Friesen et al. 2012;Dessart et al. 2013b). Meanwhile others are simply misidentified or unresolved due to the inherent high degeneracy of solutions and warrant improvements to 10 There is general consensus that the observed spectroscopic diversity of most SN Ia are influenced by: different configurations of 56 Ni produced in the events (Colgate & McKee 1969;Arnett 1982;Khokhlov et al. 1993;Baron et al. 2012), their effective temperatures (Nugent et al. 1995), density profiles and the amount of IPEs present within the outermost layers of ejecta (Hatano et al. 1999a;Baron et al. 2006;Hachinger et al. 2012), global symmetries of Si-rich material (Thomas et al. 2002), departures from spherical symmetry for Ca and Si-rich material at high velocities (Wang et al. 2007;Kasen et al. 2009;Maeda et al. 2010a;Maund et al. 2013;Moll et al. 2013;Dessart et al. 2013a), efficiencies of flux redistribution Jack et al. 2012), the radial extent of stratified material resulting from a detonation phase , host galaxy dust (Tripp & Branch 1999;Childress et al. 2013a), and the metallicity of the progenitors Lentz et al. 2000;Timmes et al. 2003;Howell et al. 2009;Bravo et al. 2010;Jackson et al. 2010;Wang et al. 2013c). 11 Some relevant obstacles include: a high radiation energy density in a low matter density environment, radiative versus local collisional processes (non-LTE conditions) and effects (Baron et al. 1996), time-dependent effects and the dominance of line over continuum opacity (Pinto & Eastman 2000a,b), and relativistic flows as well as GR effects on line profiles (Chen et al. 2007;Knop et al. 2009). In addition, the entire light emission is powered by decay-chain γ-rays, interactions with CSM, and is influenced by positrons, fast electrons, and Auger electrons in later phases (Kozma & Fransson 1992;Seitenzahl et al. 2009). the models for further study (e.g., Na I versus [Co III]; Dessart et al. 2013b).
For example, the debate over whether or not hydrogen and/or helium are detected in some early Ibc spectra has been difficult to navigate on account of the wavelength separation of observed weak features and the number of plausible interpretations (Deng et al. 2000;Branch et al. 2002b;Anupama et al. 2005a;Elmhamdi et al. 2006;Parrent et al. 2007;Ketchum et al. 2008;Soderberg et al. 2008;James & Baron 2010;Benetti et al. 2011;Chornock et al. 2011;Dessart et al. 2012;Milisavljevic et al. 2013a,b;Takaki et al. 2013). Historically, the term "conspicuous" has defined whether or not a supernova belongs to a particular spectroscopic class. By way of illustration, photographic spectrograms of type II events reveal conspicuous emission bands of hydrogen while type I events do not (Minkowski 1941). With the advent of CCD cameras in modern astronomy, it has been determined that 6300Å absorption features (however weak) in the early spectra of some type Ibc supernovae are often no less conspicuous than 6100Å Si II λ6355 absorption features in SN Ia spectra, where some 6300Å features produced by SN Ibc may be due to Si II and/or higher velocity Hα (Filippenko 1988;Filippenko et al. 1990;Filippenko 1992). That is, while SN Ibc are of the type I class, they do not necessarily lack hydrogen and/or helium within their outer-most layers of ejecta, hence the conservative definition of type I supernovae as "hydrogen/helium-poor" events.
This conundrum of which ion signatures construct each observed spectral feature rests proportionately on the signal-to-noise ratio (S/N) of the data. However, resolving this spectroscopic dilemma is primarily dependent on the wavelength and temporal coverage of the observations and traces back to the pioneering work of McLaughlin (1963) who studied spectra of the type Ib supernova, SN 1954A, in NGC 4214 (Wellmann 1955Branch 1972;Blaylock et al. 2000;Casebeer et al. 2000). Contrary to previous interpretations that supernova spectra were the result of broad, overlapping emission features (Gaposchkin 1936;Humason 1936;Baade 1936;Walter & Strohmeier 1937;Minkowski 1939;Payne-Gaposchkin & Whipple 1940;Zwicky 1942;Baade et al. 1956), it was D. B. McLaughlin who first began to repeatedly entertain the idea that "absorption-like" fea-  Pereira et al. (2013). We have normalized and over-plotted each spectrum at the 6100Å P Cygni profile in order to show the relative locations of all ill-defined features as they evolve with the expansion of the ejecta. The quoted rise-time to maximum light (dashed black) is from Mazzali et al. (2013).
tures were present 12 in regions that "lacked emission" (McLaughlin 1959(McLaughlin , 1960(McLaughlin , 1963. The inherent difficulties in reading supernova spectra and the history of uncertain line identifications for both conspicuous and concealed absorption signatures are almost as old as the supernova field itself (Payne- Gaposchkin & Whipple 1940;Dessart et al. 2013b). Still, spectroscopic intuitions can only evolve as far as the data allow. Therefore it is both appropriate and informative to recall the progression of early discussions on the spectra of supernovae, during which spectroscopic designations of type I and type II were first introduced: There appears to be a general opinion that the evidence concerning the spectrum of the most luminous nova of modern times was so contradictory that conclusions as to its spectra nature are impossible. This view is expressed, for example, by Miss Cannon: "With the testimony apparently so conflicting, it is difficult to form any conception of the class of this spectrum" (Gaposchkin 1936).
It also seems ill advised to conclude anything regarding the distribution of temperature in super-novae from the character 12 Admittedly Minkowski (1941) had previously mentioned "absorptions and broad emission bands are developed [in the spectra of supernovae]." Although, this was primarily within the context of early epoch observations that revealed a featureless, blue continuum: "Neither absorptions nor emission bands can be definitely seen but some emission is suspected in the region of Hα" (Minkowski 1940).
of their visible spectra as long as a satisfactory explanation of some of the most important features of these spectra is completely lacking (Zwicky 1936).
The spectrum is not easy to interpret, as true boundaries of the wide emission lines are difficult to determine (Humason 1936).
Those [emission] bands with distinct maxima and a fairly sharp redward or violetward edge, excepting edges due to a drop in plate spectral sensitivity, may give an indication of expansion velocity (Popper 1937).
Instead of the typical pattern of broad, diffuse emissions dominated by a band about 4600Å, it appeared like a continuum with a few deep and several shallow absorption-like minima. Two of the strongest "absorption lines," when provisionally interpreted as λλ4026, 4472 He I, give velocities near −5000 km s −1 [...] The author is grateful to N. U. Mayall and R. Minkowski for the use of spectrograms, and for helpful discussions. However, this does not imply agreement with the author's interpretations (McLaughlin 1959).
It is hardly necessary to emphasize in detail the difficulties of establishing the correct interpretation of a spectrum which may reflect unusual chemical composition, whose features may represent emission, absorption, or both mixed, and whose details are too ill-defined to admit precise measures of wavelengths (Minkowski 1963).
Given that our general understanding of blended spectral lines remains in a continual state of improvement, the frequently recurrent part of "the supernova problem" is pairing observed features with select elements of the periodic table (Hummer 1976;Axelrod 1980;Jeffery & Branch 1990;Hatano et al. 1999b;Branch et al. 2000). In fact, it was not until nearly a half-century after Minkowski (1963), with the discovery and prompt spectroscopic follow-up of SN 2011fe (Figure 1 and §4.1) that the loose self-similarity of SN Ia time series spectra from the perceived beginning of the event to near maximum light was roundly confirmed (Nugent et al. 2011, see also Garavini et al. 2005;Foley et al. 2012c;Silverman et al. 2012d;Childress et al. 2013c;Zheng et al. 2013).
While SN 2011fe may not have revealed a direct confirmation on its progenitor system (Li et al. 2011b), daily spectroscopic records at optical wavelengths were finally achieved, establishing the most efficient approach for observing ill-defined features over time (Pereira et al. 2013). This is important given that UV to NIR line identifications of all observed complexes are highly time-dependent, are sensitive to most physically relevant effects, continuously vary between subtypes, and rely on minimal constraint for all observed events 13 .
Even so, this rarely attainable observing strategy does not necessarily illuminate nor eliminate all degeneracies in spectral feature interpretations. However the advantage of complimentary high frequency follow-up observations is that the spectrum solution associated with any proposed explosion scenario can at least be consistently tested and constrained by the observed rapid changes over time ("abundance tomography" goals, e.g., Hauschildt & Baron 1999;Stehle et al. 2005;Sauer et al. 2006;Kasen et al. 2006;Pauldrach et al. 2013). It then follows that hundreds of well-observed spectrophotometric datasets serve to carve out the characteristic information, f (λ; t), for each SN Ia between subtypes, in addition to establishing the perceived boundaries of the SN Ia diversity problem (see Fig. 11 of Blondin et al. 2012 for this concept at maximum light).
For supernovae in general, Figure 1 also serves as a reminder that all relative strengths evolve continuously over time, where entire features are always redshifting across wavelength (line velocity space) during the rise and fall in brightness. A corollary of this situation is that prescriptions for taking measurements of spectroscopic behaviors (whereby interpretations rely on a subjective "goodness of fit") and robustly associating with any number of physical causes do not exist. Instead there are two primary means for interpreting SN Ia spectra and taking measurements of features for the purposes of extracting physical properties.
13 See Foley et al. (2012b) for "The First Maximum-light Ultraviolet through Near-infrared Spectrum of a Type Ia Supernova." Fig. 2 Top: A schematic representation of how an assumed spherically sharp and embedded photosphere amounts to a pure line-resonance P Cygni profile under the conditions of Sobolev line transfer within a geometry of Absorbing, Emitting, and Occulted regions of material (Jeffery & Branch 1990;Branch et al. 2005). The approximate photospheric velocity, v phot , is proportional to the blue ward shift of an unblended absorption minimum. Bottom: Application of the above P Cygni diagram to SN Ia spectra in terms of which species dominate and what other species are known to influence the temporal behavior (Bongard et al. 2008), each of which are constrainable from complete spectroscopic coverage. For each series of spectra, the black line in bold represents maximum light.
Indirect analysis assumes a detailed explosion model and is primarily tasked with assessing the accuracy and flaws of the model. Direct analysis seeks to manually measure via spectrum synthesis where one can either assume an initial post-explosion ejecta composition or give up abundance information altogether to assess the associated uncertainties and consequences of supernova line blending via purposeful high parameterizations. For the latter of these direct inference methods, the conclusions about spectroscopic interpretations−which are supported by remnants of inconsistencies throughout the literature−are summarized as follows.
For the most part, particularly at early epochs and as far as anyone can tell with current limiting datasets, the features in SN Ia spectra are due to IMEs and IPEs formed by resonance scattering of continuum and decay-chain photons, and have P Cygni-type profiles overall (Pskovskii 1969;Mustel 1971;Branch & Patchett 1973;Kirshner et al. 1973a; see Figure 2). Emission components peak at or near the rest wavelength and absorption components are blue-shifted according to the opacity profile of matter at and above the photospheric line forming region. The combination of these effects can often lead to "trumped" emission features (Jeffery & Branch 1990), giving SN Ia spectra their familiar shapes.
Essentially all relevant atomic species (isotope plus ionization state) are present somewhere within the ejecta, each with its own 3-dimensional abundance profile. At optical wavelengths, conditions and abundance tomographies of the ejecta maintain the dominance of select singly−triply ionized subsets of C+O, IMEs, and IPEs (Hatano et al. 1999b). From shortly after the onset of the explosion to around the time of maximum light, the optical−NIR spectrum of a normal SN Ia consists of a continuum level with superimposed features that are primarily consistent with strong permitted lines of ions such as O I, Mg II, Si II, Si III, S II, Ca II, Fe II, Fe III, and trace signatures of C I and C II Thomas et al. 2007;Bongard et al. 2008;Nugent et al. 2011;Parrent et al. 2012;Hsiao et al. 2013;Mazzali et al. 2013;Dessart et al. 2013a). After the pre-maximum light phase, blends of Fe II (in addition to other IPEs) begin to dominate or influence the temporal behavior of many optical−NIR features over timescales from weeks to months (see  and references therein).
With the above mentioned approximated view of line formation in mind (Figure 2), the real truth is that the time-dependent state of the ejecta and radiation field at all locations dictates how the material presence within the line forming regions will be imparted onto the spectral continuum, i.e. the radiation field and the matter are said to be "coupled." With the additional condition of near-relativistic expansion velocities (∼0.1c), line identifications themselves can also be thought of as coupled to the abundance tomography of ejected material, which includes the projected Doppler velocities spanned by the recipe of absorbing material. Subsequently, while spectra can be used for constraining limits of some model parameters, it comes with a cost of certainty on account of natural uncertainties imparted by the large expansion velocities and associated expansion opacities.
As an exercise in this point, in Figure 3 we have constructed an early epoch set of toy model line profiles that are representative of normal SN Ia line identification procedures (e.g., Branch et al. 2005;Parrent et al. 2011) and over-plot them with an early optical−NIR spectrum (the observed outermost layers, sans UV) of SN 2011fe. We summarize the take away points of Figure 3 as follows.
• Even without considering weak contributions, at no place along the (UV−) optical−NIR spectrum is any observed feature removed from being due to less than 2 sources (more precisely, see also Bongard et al. 2008). That is, under the basic assumptions of pure resonance line scattering and homologous expansion (Figure 2), all features are complex blends of at least 2+ ions and are universally influenced by multiple regions of emitting and/or absorbing material (e.g., "high[-er] velocity" and "photospheric velocity" intervals of material, see also Marion et al. 2013). • For supernovae, the components of the spectrum are most easily constrained via spectrum synthesis, and subsequently measurable (not the converse), when the bounds of wavelength coverage, λ a and λ b , are between ∼2000−3500 and 12000Å, respectively. If λ b < 7500−9500Å, then the velocities and relative strengths of several physically relevant ions (e.g., C I, O I, Mg II, and Ca II) are said to be devoid of useful constraint and provide a null (or uncertain) measurement for every other overlapping spectral line signature (i.e. all features). That is, in order to viably "identify" and measure a single feature, the entire spectrum must be reproduced. While empirical measurements of certain absorption features are extremely useful for identifying trends in the observed behavior of SN Ia, these methods do not suffice to measure the truest underlying atomic recipe and its time-dependent behavior, much less the "strength" of contributing lines (e.g., multiple velocity components of Si (Hsiao et al. 2013;Pereira et al. 2013). Calculations are based on an optical set of photospheric phase spectra (see Parrent et al. 2012) and are true-to-scale. Bands of color are intended to show overlap between lines under the simplified however informative assumption of permitted line scattering under homologous expansion. Some of the weaker lines have not been highlighted for clarity.
2008) or the use of standardized relative strength parameters (e.g., Childress et al. 2013b). • Therefore, as in Figure 2, employing stacked Doppler velocity scaled time series spectra provides useful and timely first-order comparative estimates for when (epoch) and where (projected Doppler velocity) contributing ions appear, disappear, and span as the photospheric region recedes inward over time.
We speak on this only to point out that even simple questions−particularly for homogeneous SN Ia−are awash in detection/non-detection ambiguities. However, it should be noted that a powerful exercise in testing uncertain line identifications and resolving complex blends can be done, in part, without the use of additional synthetic spectrum calculations. That is, by comparing a single observed spectrum to that of other well-observed SN Ia, where the analysis of the latter offers a greater context for interpretation than the single spectrum itself, one can deduce whether or not a "mystery" absorption feature is common to most SN Ia in general. On the other hand, if a matching absorption feature is not found, then one can infer the presence of either a newly identified, compositionally consistent ion or the unblended line of an already accounted for species (resulting from forbidden line emission, non-LTE effects, and/or when line strengths or expansion velocities differ between subtypes). Given also the intrinsic dispersion of expansion opacities between SN Ia, it is likely that an "unidentified" feature is that of a previously known ion at higher and/or lower velocities. It is this interplay between expansion opacities and blended absorption features that keep normal and some peculiar SN Ia within the description of a homogenous set of objects, however different they may appear.
In fact, when one compares the time series spectra of a broad sample of SN Ia subtypes, however blended, there is little room for degeneracy among plausible ion assignments (sans IPEs, e.g., Fe II versus Cr II during post-maximum phases). In other words, there exists a unique set of ions, common to most SN Ia atmospheres, that make up the resulting spectrum, where differences in subtype are associated with differences in temperature and/or the abundance tomography of the outermost layers (Tanaka et al. 2008). The atomic species listed in Figure 3 do not so much represent a complete account of the composition, or the "correct" answer, as they are consistent with the subsequent time evolution of the spectrum toward maximum light, and therefore serve to construct characteristic standards for direct comparative diversity assessments.
Said another way, it is the full time series dataset that enables the best initial spectrum solution hypothesis, which can be further tested and refined for the approximate measurement of SN Ia features (Branch et al. 2007a). Therefore, this idea of a unique set of ions remains open since−with current limiting datasets−species with minimal constraint and competing line transfer processes can be ambiguously present 14 , even for data with an infinite S/N (i.e. sources with few strong lines, or lines predominately found blue ward of ∼6100Å, e.g., C III, O III, Si IV, Fe I, Co II, Ni II). One can still circumvent these uncertainties of direct analysis by either using dense time series observations (e.g., Parrent et al. 2012) or by ruling out spurious inferred detections by including adjacent wavelength regions into the spectroscopic analysis (UV−optical−NIR; see Foley et al. 2012b;Hsiao et al. 2013;Mazzali et al. 2013).

Models
A detailed account of SN Ia models is beyond the scope of our general review of SN Ia spectra (for the latest discussions, see Wang & Han 2012;Nomoto et al. 2013;Hillebrandt et al. 2013;Calder et al. 2013;Maoz et al. 2013). However, in order to understand the context by which observations are taken and synthetic comparisons made, here we only mention the surface layer of matters relating to observed spectra. For some additional recent modeling work, see Fryer & Diehl (2008) Moll et al. (2013), and .
Realistic models are not yet fully ready because of the complicated physical conditions in the binary stellar evolution that leads up to an expanding SN Ia atmosphere. For instance, the explosive conditions of the SN Ia problem take place over a large dynamic range of relevant length-scales (R W D ∼ 1R ⊕ and flamethicknesses of ∼ 0.1 cm; Timmes & Woosley 1992;, involve turbulent flames that are fundamentally multi-dimensional (Khokhlov 1995(Khokhlov , 2000Reinecke et al. 2002a,b;Gamezo et al. 2003Gamezo et al. , 2005Seitenzahl et al. 2013), and consist of uncertainties in both the detonation velocity (Domínguez & Khokhlov 2011) and certain nuclear reaction rates, especially 12 C+ 12 C (Bravo et al. 2011, however see also Bravo & Martínez-Pinedo 2012;Chen et al. 2013).
Most synthetic spectra are angle-averaged representations of higher-dimensional detailed models. Overall, 14 See Fig. 9 of Stritzinger et al. (2013) to see clear detections of permitted Co II lines in the NIR spectra of the peculiar and faint SN 2010ae. the observed spectra of normal SN Ia have differed less amongst themselves than that of some detailed models compared to the data of normal SN Ia. This is not from a lack of efforts, but is simply telling of the inherent difficulty of the problem and limiting assumptions and interests of various calculations. Kasen et al. (2008) reviewed previous work done of N-dimensional SN Ia models and presented the first high-resolution 3D calculation of a SN Ia spectrum at maximum light. Their results are still in a state of infancy, however they represent the first step toward the ultimate goal of SN Ia modeling, i.e. to trace observed SN Ia properties and infer the details of the progenitor and its subsequent disruption by comparing 3D model spectra and light curves of 3D explosion simulations with the best observed temporal datasets.
Still, progress has been made in understanding general observed properties of SN Ia and their relation to predictions of simulated explosion models. For example, one-dimensional (1D) numerical models of SN Ia have been used in the past to test the possible explosion mechanisms such as subsonic flame or supersonic detonation models, as well as conjoined delayed-detonations (e.g., Arnett 1968;Nomoto et al. 1984;Lentz et al. 2001a). The one-dimensional models disfavor the route of a pure thermonuclear detonation as the mechanism to explain most SN Ia events (Hansen & Wheeler 1969;Arnett 1969;Axelrod 1980). Such a mechanism produces mostly 56 Ni and almost none of the IMEs observed in the spectra of all SN Ia (e.g., Branch et al. 1982;Filippenko 1997;Pastorello et al. 2007a).
However, one-dimensional models have shown that a detonation can produce intermediate mass elements if it propagates through a Chandrasekhar-mass WD that has pre-expanded during an initial deflagration stage (Khokhlov 1991a;Yamaoka et al. 1992;Khokhlov et al. 1993;Arnett & Livne 1994a,b;Wheeler et al. 1995;Hoflich et al. 1995;Khokhlov et al. 1997). To their advantage, these deflagration-to-detonation transition (DDT) and pulsating delayed-detonation (PDD) models are able to reproduce the observed characteristics of SN Ia, however not without the use of an artificially-set transition density between stages of burning (Khokhlov 1991b;Hoflich et al. 1995;Lentz et al. 2001a,b;Baron et al. 2008;Bravo et al. 2009;Dessart et al. 2013a). Subsequently, a bulk of the efforts within the modeling community has been the pursuit of conditions or mechanisms which cause the burning front to naturally transition from a sub-sonic deflagration to a super-sonic detonation, e.g., gravitationally confined detonations (Jordan et al. 2009), prompt detonations of merging WDs, a.k.a. "peri-mergers" .
With the additional possibility that the effectively burned portion of the progenitor is enclosed or obscured by some body of circumstellar or envelope/disk of material (see Sternberg et al. 2011;Foley et al. 2012a;Förster et al. 2012;Scalzo et al. 2012;Silverman et al. 2013d;Dan et al. 2013;Dessart et al. 2013a;Moll et al. 2013), the intrinsically multidimensional nature of the explosion itself is also expected to manifest signatures of asymmetric plumes of burned material and pockets of unburned material within a spheroidal debris field of flexible asymmetries (see Khokhlov 1995;Niemeyer & Hillebrandt 1995;Gamezo et al. 2004;Wang & Wheeler 2008;Patat et al. 2009;Kasen et al. 2009). Add to this the degeneracy of SN Ia flux behaviors, i.e. colors are sensitive to dust/CSM extinction and intrinsic dispersions in the same direction (Tripp & Branch 1999), whether large or small redshift-color dependencies (Saha et al. 1999;Jha et al. 1999;Parodi et al. 2000;Wang et al. 2008a;Goobar 2008;Wang et al. 2009a;Foley & Kasen 2011;Mohlabeng & Ralston 2013), and we find the true difficulty in constraining SN Ia models. Blondin et al. (2013) recently presented and discussed the photometric and spectroscopic properties at maximum light of a sequence of 1D DDT explosion models, with ranges of synthesized 56 Ni masses between 0.18 and 0.81 M . In addition to showing broad consistencies with the diverse array of observed SN Ia properties, the synthetic spectra of Blondin et al. (2013) predict weaker absorption features of unburned oxygen (O I λ7774) at maximum light, in proportion to the amount of 56 Ni produced. This is to be expected (Hoflich et al. 1995), however constraints on the remaining amount of unburned material, in addition to its temporal behavior, are more readily seen during the earliest epochs (within the outermost layers of ejecta) via C II λ6580 and O I λ7774 (Thomas et al. 2007;Parrent et al. 2011;Nugent et al. 2011). Consequently, temporal spectrum calculations of detailed explosion models are needed for the purposes of understanding why the properties of SN Ia are most divergent well before maximum light Dessart et al. 2013a).
Nucleosynthesis in two-dimensional (2D) delayed detonation models of SN Ia were explored by Maeda et al. (2010a). In particular, they focused on the distribution of species in an off-center DDT model and found the abundance tomography to be stratified, with an inner region of 56 Ni surrounded by an off-center shell of electron-capture elements (e.g., Fe 54 , Ni 58 ). Later, Maeda et al. (2010b) investigated the late time emission profiles associated with this off-center inner-shell of material within several observed SN Ia and found a correlation between possible nebular-line Doppler shifts along the line-of-sight and the rate-of-decline of Si II velocities at earlier epochs. Their interpretation is to suggest that some SN Ia subtypes may represent two different hemispheres of the "same" SN Ia (LVG vs. HVG subtypes; see §3.2). Moreover, the findings of Maeda et al. (2010b) and Maund et al. (2010b) remain largely consistent with the additional early and late time observations of the well-observed SN 2011fe (Smith et al. 2011;McClelland et al. 2013) and those of larger SN Ia samples (Blondin et al. 2012;Silverman et al. 2013a). However, even the results of Maeda et al. (2010b) and others that rely on spectroscopic measurements at all epochs are not without reservation given that late time emission profiles are subject to more than line-shifts due to Doppler velocities and ionization balance (Bongard et al. 2008;Friesen et al. 2012). Seitenzahl et al. (2013) presented 14 3-dimensional (3D) high resolution Chandrasekhar-mass delayeddetonations that produce a range of 56 Ni (depending on the location of ignition points) between ∼ 0.3 and 1.1 M . For this set of models, unburned carbon extends down to 4000 km s −1 while oxygen is not present below 10,000 km s −1 . Seitenzahl et al. (2013) conclude that if delayed-detonations are to viably produce normal SN Ia brightnesses, the region of ignition cannot be far off-center so as to avoid the over-production of 56 Ni. As noted by Seitenzahl et al. (2013), these models warrant tests via spectrum synthesis given their 3D nature and possible predictive relations to the WLR, spectropolarimetry, and C+O "footprints" Baron et al. 2003;Thomas et al. 2007;Wang & Wheeler 2008). Dessart et al. (2013a) recently compared synthetic light curves and spectra of a suite of DDT and PDD models. Based on comparisons to SN 2002bo and SN 2011fe, two SN Ia of different spectroscopic subtypes, and based on poor to moderate agreement between recent DDT models and observed SN Ia diversity (Blondin et al. 2011), Dessart et al. (2013a convincingly argue that these two SN Ia varieties (LVG vs. HVG, as above) are dissimilar enough to be explained by different explosion scenarios and/or progenitor systems (Wang et al. 2013c). For SN Ia in general, delineating spectroscopic diversity has been a difficult issue Branch et al. 2009), and has only recently been made clear with the belated release of decades-worth of unpublished data (Blondin et al. 2012;Silverman et al. 2012c).

Spectroscopic Diversity of SN Ia
Observationally and particularly at optical wavelengths, SN Ia increase in brightness over ∼13 to 23 days before reaching maximum light (t rise = 17.38 ± 0.17; Hayden et al. 2010b). However, it is not until ∼1 year later that the period of observation is said to be "complete." From the time of the explosion our perspective as outside observers begins at the outermost layers if the SN Ia is caught early enough. In the approximate sense, this is because the line-forming region (the "photosphere") recedes as the ejecta expand outward, which in turn means that the characteristic information for each explosion mechanism and progenitor channel is specified by the temporal spectrophotometric attributes of the "inner" and "outer" layers of freshly synthesized and remaining primordial material. In addition, because the expanding material cools as it expands, the net flux of photons samples different layers (of different states and distributions) over time. And since the density profile of the material roughly declines from the center outward, significant changes within the spectra for an individual SN Ia take place daily before or near maximum light, and weekly to monthly thereafter.
Documenting the breadth of temporal spectroscopic properties for each SN Ia is not only useful for theoretical purposes, but is also necessary for efficiently typing and estimating the epoch of newly found possible supernova candidates before they reach maximum light. Several supernova identification tools have been made that allow for fair estimates of both subtype and epoch (e.g., SNID; Blondin & Tonry 2007, Gelato;Harutyunyan et al. 2008, Superfit;Howell et al. 2005). In addition, the spectroscopic goodness-of-fit methods of  allow one to find the "nearest neighbors" of any particular SN Ia within a sample of objects, enabling the study of so called "transitional subtype" SN Ia (those attributed with contrasting characteristics of two or more subtypes).

Data
One of the major limitations of spectroscopic studies has been data quality. For example, the signal-tonoise ratio, S/N, of a spectrum signifies the quality across wavelength and is usually moderate to high for high-z events. Similarly, and at least for low-z SN Ia, there should exist a quantity that specifies the density of spectra within a time series dataset. We suggest S/N •(P) ≡ the number of continual follow-up spectra / the mean number of nights passed between exposures • (total number of spectra prior to maximum light). In Figure 4 we apply this quantity to literature data.
An ideal dataset consisting of 25 spectra during the first 25 days post-explosion would yield S/N •(P) = 25 (16) (e.g., SN 2011fe), whereas a dataset of spectra Continual spectroscopic follow-up efficiencies for the most "well-observed" SN Ia at early phases (not counting multiple spectra per day). Some of the values reported may be slightly lower for instances of unpublished data. Dashed lines represent the upper-limit spectroscopic efficiencies and peak number of pre-maximum light spectra for one and two day follow-up cadences during the first 25 days post-explosion. See §3.1.
at days −12, −10, −7, −4, +0, +3, +8, +21, +48, +119 (a common occurrence) would be said to have S/N •(P) = 3.3 (4) plus follow-up at days +21, +48, and +119. By including the total number of spectra prior to maximum light in parentheses, we are anticipating those cases where S/N = 1, but with P = 3, e.g., a dataset with days −12, −9, and −6 observed. It may serve a purpose to also add second and third terms to this quantity that take into account the number of post-maximum light and late time spectra.
Regardless of moniker and definition, a quantity that specifies the density of spectra observed during the earliest epochs would aid in determining, quantitatively, which datasets are most valuable for various SN Ia diversity studies. Clearly such a high follow-up rate for slow-evolving events (e.g., SN 2009dc) or events caught at maximum light are not as imperative. However, when SN Ia are found and typed early, a high S/N ensures no loss of highly time sensitive information, e.g., when high velocity features and C+O signatures dissipate. Since most datasets are less than ideal for detailed temporal inspections of many events (by default), astronomers have instead relied upon comparative studies ( §3.2); those that maximize sample sizes by prioritizing the most commonly available spectroscopic observables, e.g., line velocities of 6100Å absorption minima near maximum light.
Another limitation of spectroscopic studies has been the localized release of all published data. The Online Supernova Spectrum Archive (SuSpect 15 ; Richardson et al. 2001) carried the weight of addressing data foraging during the past decade, collecting a total of 867 SN Ia spectra (1741 SN spectra in all). Many of these were either at the request of or donation to SuSpect, while some other spectra were digitized from original publications in addition to original photographic plates (Casebeer et al. 1998. Prior to and concurrent with SuSpect, D. Jeffery managed a collection of SUpernova spectra PENDing further analysis (SUSPEND 16 ).
With the growing need for a manageable influx of data, the Weizmann Interactive Supernova Data Repository (WISeREP 17 ; Yaron & Gal-Yam 2012) has since served as a replacement and ideal central data hub, and has increased the number of SN Ia spectra to 7661 (with 7933 publicly available SN spectra out of 13,334 in all). We encourage all groups to upload published data to WISeREP, whether or not made available elsewhere.

Samples
By far the largest data releases occurred during the past five years, and are available on WISeREP and their affiliated archives. Matheson et al. (2008)  These samples provide a substantial improvement and crux by which to explore particular issues of SN Ia diversity. However, the remaining limitation is that our routine data collection efforts continue to yield several thousand SN Ia with few to several spectra by which to dissect and compare SN Ia atmospheres.

Comparisons of "Well-Observed" SN Ia
Given that both quantitative and qualitative spectrum comparisons are at the heart of SN Ia diversity studies, in Figure 5 - Figure 13 we plot spectroscopic temporal snapshots for as many "well-observed" SN Ia as are currently available on WISeREP (Tables 1 and 2). Because the decline parameter, ∆m 15 (B ), remains a useful parameter for probing differences of synthesized 56 Ni mass, properties of the ejecta, limits of CSM interaction, etc., we have loosely ordered the spectra with increasing ∆m 15 (B ) (top-down) based on average values found throughout the literature (Tables 3−6) and M B (peak) considerations for cases that are reported as having the same ∆m 15 (B ). The spectra have been normalized with respect to 6100Å line profiles in order to amplify relative strengths of the remaining features (see caption of Figure 5). We also denote the spectroscopic subtype for each object in color in order to show the overlap of these properties between particular SN Ia subclasses (see §3.2 and Blondin et al. 2012).
By inspection, the collected spectra show how altogether different and similar SN Ia (both odd and normal varieties) have come to be since nearly 32 years ago. With regard to the recent modeling of Blondin et al. (2013) and their accompanying synthetic spectra, we plot the spectra in Figure 5 - Figure 11 in the flux-representation of λ 2 F λ for ease of future comparisons. These juxtapositions should reveal the severity 4000 5000 6000 7000 8000 9000 Rest Wavelength [Å]  Figure 8 for comparison, as it is not a prototypical SN Ia, however appearing similar to SN 1991bg-like events (Sullivan et al. 2010;Kasliwal et al. 2012). ! "#    of the SN Ia diversity problem as well as the future of promising studies and work that lie ahead.

Deciphering 21 st Century SN Ia Subtypes
Observationally, the whole of SN Ia are hetero-, homogeneous events (Oke & Searle 1974;Filippenko 1997); some of the observed differences in their spectra are clear, while other suspected differences are small enough to fall below associable certainty. Because of this, observational studies have concentrated on quantitatively organizing a mapping between the most peculiar and normal events. In this section our aim is to review SN Ia subtypes. In all, three observational classification schemes will be discussed Wang et al. 2009a), as well as the recent additions of so-called over-and sub-luminous events (see Scalzo et al. 2012, Foley et al. 2013, Silverman et al. 2013d and references therein). For other relatively new and truly peculiar subclasses of supernova transients, we refer the reader to Shen et al. (2010), Kasliwal et al. (2012) and references therein.

Benetti et al. (2005) Classification
Understanding the origin of the WLR is a key issue for understanding the diversity of SN Ia light curves and spectra, as well as their use as cosmological distance indicators. Brighter SN Ia with broader light curves tend to occur in late-type spiral galaxies, while dimmer, faster declining SN Ia are preferentially located in an older stellar population and thus the age and/or metallicity of the progenitor system may be relevant factors affecting SN Ia properties (Hamuy et al. 1995;Howell 2001;Pan et al. 2013, see also Hicken et al. (2009a)). With this in mind, Benetti et al. (2005) studied the observational properties of 26 well-observed SN Ia (e.g., SN 1984A, 1991T, 1991bg, 1994D) with the intent of exploring SN Ia diversity. Based on the observed projected Doppler velocity evolution from the spectra 18 , in conjunction with characteristics of the light curve (M B , ∆m 15 ), Benetti et al. (2005) considered three different groups of SN Ia: (1) "FAINT" SN 1991bg-likes, (2) "low velocity gradient" (LVG) SN 1991T/1994D-likes, and (3) "high velocity gradient" (HVG) SN 1984A-like events. The velocity gradient here is based on the time-evolution of 6100 ("6150") A absorption minima as inferred from Si II λ6355 line 18 The velocity gradient−the mean velocity decline rate ∆v/∆t−of a particular absorption minimum (e.g.,v Si II ) has been redefined to be measured over a fixed phase range [t 0 , t 1 ] (Blondin et al. 2012).
velocities. Overall, HVG SN Ia have higher mean expansion velocities than FAINT and LVG SN Ia, while LVG SN Ia are brighter than FAINT and HVG SN Ia on average (Silverman et al. 2012b;Blondin et al. 2012). Given an apparent separation of SN Ia subgroups from this sample of 26 objects, Benetti et al. (2005) considered it as evidence that LVG, HVG, and FAINT classifications signify three distinct kinds of SN Ia. Branch et al. (2005Branch et al. ( , 2007bBranch et al. ( , 2009) published a series of papers based on systematic, comprehensive, and comparative direct analysis of normal and peculiar SN Ia spectra at various epochs with the parameterized supernova synthetic spectrum code, SYNOW 19 (Fisher 2000;Branch et al. 2007a). From the systematic analysis of 26 spectra of SN 1994D, Branch et al. (2005 infer a compositional structure that is radially stratified, overall. In addition, several features are consistent with being due to permitted lines well into the late postmaximum phases (∼120 days, see Friesen et al. 2012). Another highlight of this work is that, barring the usual short comings of the model, SYNOW is shown to provide a necessary consistency in the direct quantification of spectroscopic diversity (Branch et al. 2007a). Consequently, the SYNOW model has been useful for assessing the basic limits of a spectroscopic "goodness of fit" ( In their second paper of the series on comparative direct analysis of SN Ia spectra,  studied the spectra of 24 SN Ia close to maximum light. Based on empirical pEW measurements of 5750, 6100Å absorption features, in addition to spectroscopic modeling with SYNOW,  organized SN Ia diversity by four spectroscopic patterns: (1) "Core-Normal" (CN) SN 1994D-likes, (2) "Broad-line" (BL), where one of the most extreme cases is SN 1984A, (3) "Cool" (CL) SN 1991bg-likes, and (4) "Shallow-Silicon" (SS) SN 1991T-likes. In this manner, a particular SN Ia is defined by its spectroscopic similarity to one or more SN Ia prototype via 5750, 6100Å features. These spectroscopic subclasses also materialized from analysis of pre-maximum light spectra (Branch et al. 2007b).

Branch et al. (2006) Classification
The overlap between both Benetti et al. (2005) and  classifications schemes comes by comparing Table 1 in Benetti et al. (2005) to Table  1 of , and it reveals the following SN Ia descriptors: HVG−BL, LVG−CN, LVG−SS, and FAINT−CL. This holds true throughout the subsequent literature (Branch et al. 2009;Folatelli et al. 2012;Blondin et al. 2012;Silverman et al. 2012b).
In contrast with Benetti et al. (2005) who interpreted FAINT, LVG, and HVG to correspond to the "discrete grouping" of distinctly separate SN Ia origins among these subtypes,  found a continuous distribution of properties between the four subclasses defined above. We should point out that this classification scheme of  is primarily tied to the notion that SN Ia spectroscopic diversity is related to the temperature sequence found by Nugent et al. 1995. That is, despite the contrast with Benetti et al. (2005) (continuous versus discrete subgrouping of SN Ia), so far these classifications say more about the state of the ejecta than the various number of possible progenitor systems and/or explosion mechanisms (see also Dessart et al. 2013a). Furthermore, the existence of "transitional" subtype events support this notion (e.g., SN 2004eo, 2006bt, 2009ig, 2001ay, and PTF10ops; see appendix). Branch et al. (2009) later analyzed a larger sample of SN Ia spectra. They found that SN 1991bg-likes are not a physically distinct subgroup (Doull & Baron 2011), and that there are probably many SN 1999aalike events (A.5) that similarly may not constitute a physically distinct variety of SN Ia.
With regard to the fainter variety of SN Ia, Doull & Baron (2011) made detailed comparative analysis of spectra of peculiar SN 1991bg-likes. They also studied the intermediates, such as SN 2004eo (A.23), and discussed the spectroscopic subgroup distribution of SN Ia. The CL SN Ia are dim, undergo a rapid decline in luminosity, and produce significantly less 56 Ni than normal SN Ia. They also have an unusually deep and wide trough in their spectra around 4200Å suspected as due to Ti II (Filippenko et al. 1992b), in addition to a relatively strong 5750Å absorption (due to more than Si II λ5972; see Bongard et al. 2008). Doull & Baron (2011) analyzed the spectra of SN 1991bg, 1997cn, 1999by, and 2005bl using SYNOW, and found this group of SN Ia to be fairly homogeneous, with many of the blue spectral features well fit by Fe II.

Wang et al. (2009a) Classification
Based on the maximum light expansion velocities inferred from Si II λ6355 absorption minimum line velocities, Wang et al. (2009a) studied 158 SN Ia, separating them into two groups called "high velocity" (HV) and "normal velocity" (NV). This classification scheme is similar to those previous of Benetti et al. (2005) and , where NV and HV SN Ia are akin to LVG−CN and HVG−BL SN Ia, respectively. That is, while the subtype notations differ among authors, memberships between these classification schemes are roughly equivalent (apart from outliers such as the HV-CN SN 2009ig, see Blondin et al. 2012).
Explicitly, Benetti et al. (2005) and Wang et al. (2009a) subclassifications are based on empirically estimated mean expansion velocities near maximum light (±4 days; ±500 − 2000 km s −1 ) of 6100Å features produced by an assumed single broad component of Si II. The notion of a single photospheric layer, much less a single-epoch snapshot, does not realistically account for the multilayered nature of spectrum formation (Bongard et al. 2008), its subsequent evolution postmaximum light Scalzo et al. 2012), and potential relations to line-of-sight considerations (Maeda et al. 2010b;Blondin et al. 2011;Moll et al. 2013). In the strictest sense of SN Ia sub-classification, "normal" refers to both of these subtypes since they differ foremost by a continuum of inferred mean expansion velocities and the extent of expansion opacities, simultaneously.
Furthermore, note from a sample of 13 LVG and 8 HVG SN Ia that Benetti et al. (2005) found 10 v Si II (km s −1 day −1 ) 67 (±7) and 75 v Si II 125 (±20) for each, respectively. Similarly, and from a sample of 14 LVG and 29 HVG SN Ia, Silverman et al. (2012b) report that 10 v Si II 445 (±50) and 15 v Si II 290 (±140) for LVG and HVG events, respectively. Additionally, the pEW measurements of 5750, 6100Å absorption features (among others) are seen to share a common convergence in observed values Hachinger et al. 2006;Blondin et al. 2012;Silverman et al. 2012b). The continually consistent overlap between the measured properties for these two SN Ia "subtypes" implies that the notion of a characteristic separation value forv Si II ∼ 70 km s −1 day −1 (including the inferred maximum light separation velocity, v 0 12,000 km s −1 ) is still devoid of any physical significance beyond overlapping bimodal distributions of LVG−CN and HVG−BL SN Ia properties (see §5.3 of Silverman et al. 2012b, §5.2 of Blondin et al. 2012, andSilverman et al. 2012a). Rather, a continuum of empirically measured properties exists between the extremities of these two particular historically-based SN Ia classes (e.g., SN 1984A and1994D). Given also the natural likelihood for a physical continuum between NV and HV subgroups, considerable care needs to be taken when concluding on underlying connections to progenitor systems from under-observed, early epoch snapshots of blended 6100Å absorption minima.
Hence, the primary obstacle within SN Ia diversity studies has been that it is not yet clear if the expanse of all observed characteristics of each subtype has been fully charted. For the observed properties of normal SN Ia, it is at least true thatv Si II resides between 10−445 km s −1 day −1 , with a median value of ∼ 60−120 km s −1 day −1 Blondin et al. 2012;Silverman et al. 2012a), while the rise to peak B -band brightness ranges from 16.3 to 19 days Mazzali et al. 2013).
Recently, Wang et al. (2013c) applied this NV and HV subgrouping to 123 "Branch normal" SN Ia with known positions within their host galaxies and report that HV SN Ia more often inhabit the central and brighter regions of their hosts than NV SN Ia. This appears to suggest that a supernova with "higher velocities at maximum light" is primarily a consequence of a progenitor with larger than solar metallicities, or that PDD/HVG SN Ia are primarily found within the galactic distribution of DDT/LVG SN Ia (c.f. Blondin et al. 2011Blondin et al. , 2012Dessart et al. 2013a). This is seemingly in contrast to interpretations of Maeda et al. (2010b) who propose, based on both early epoch and late time considerations, that LVG and HVG SN Ia are possibly one in the same event where the LVG-to-HVG transition is ascribed to an off-center ignition.
While it is true that increasing the C+O layer metallicity can affect the blueshift of the 6100Å absorption feature−in addition to lower temperatures and increased UV line-blocking−this is not primarily responsible for the shift in 6100Å absorption minima (Lentz et al. 2001a,b), where the dependence of this effect is not easily decoupled from changes in the temperature structure . However, it is also worthwhile to point out that, while the early epoch spectra of SN 2011fe (a NV event) are consistent with a DDTlike composition with a sub-solar C+O layer metallicity ("W7 + ," Mazzali et al. 2013) and a PDD-like composition (Dessart et al. 2013a), the outermost layers of SN 2010jn (a HV event; A.41) are practically void of unburned material and subsequently already overabundant in synthesized metals for progenitor metallicity to be well determined (Hachinger et al. 2013). Therefore, discrepancies between NV and HV SN Ia must still be largely dependent on more than a single parameter, e.g. differences in explosion mechanisms (Dessart et al. 2013a;Moll et al. 2013), where progenitor metallicity is likely to be only one of several factors influencing the dispersions of each subgroup Höflich et al. 2010;Wang et al. 2012).
It should be acknowledged again that metallicitydependent aspects of stellar evolution are expected to contribute, in part, to the underlying variance of holistic SN Ia characteristics. However thus far, the seen discrepancies from metallicities share similarly uncertain degrees of influence as for asymmetry and lineof-sight considerations of ejecta-CSM interactions for a wide variety of SN Ia Kasen et al. 2003;Leloudas et al. 2013). Similar to this route of interpretation for SN Ia subtypes are active galactic nuclei and the significance of the broad absorption line quasi-stellar objects (BALQSOs, see de Kool & Begelman 1995;Becker et al. 1997;Elvis 2000;Branch et al. 2002a;Hamann & Sabra 2004;Casebeer et al. 2008;Leighly et al. 2009;Elvis 2012).

Additional Peculiar SN Ia Subtypes
Spectroscopically akin to some luminous SS SN Ia are a growing group of events thought to be "twice as massive," aka super-Chandrasekhar candidates (SCC, Howell et al. 2006;Jeffery et al. 2006;Hillebrandt et al. 2007;Hicken et al. 2007;Maeda et al. 2009;Chen & Li 2009;Yamanaka et al. 2009a;Scalzo et al. 2010;Tanaka et al. 2010;Yuan et al. 2010;Silverman et al. 2011;Taubenberger et al. 2011;Kamiya et al. 2012;Scalzo et al. 2012;Hachinger et al. 2012;Yamanaka et al. 2013). Little is known about this particular class of over-luminous events, which is partly due to there having been only a handful of events studied. Thus far, SCC SN Ia are associated with metal-poor environments (Childress et al. 2011;Khan et al. 2011a). Spectroscopically, the differences that set these events apart from normal SN Ia are fairly weak Si II/Ca II signatures and strong C II absorption features relative to the strength of Si II lines. Most other features are comparable in relative strengths to those of normal SN Ia, if not muted by either toplighting or effects of CSM interaction Leloudas et al. 2013), and are less blended overall due to lower mean expansion velocities. In addition, there is little evidence to suggest that SCC SN Ia spectra consist of contributions from physically separate high velocity regions of material ( 4000 km s −1 above photospheric). This range of low expansion velocities (∼5000−18,000 km s −1 ), in conjunction with larger than normal C II absorption signatures, are difficult to explain with some M Ch explosion models Kamiya et al. 2012, however see also Hachisu et al. 2012;Dessart et al. 2013a;Moll et al. 2013 for related discussions). Silverman et al. (2013d) recently searched the BSNIP and PTF datasets, in addition to the literature sample, and compiled a list of 16 strongly CSM interacting SN Ia (referred to as "Ia-CSM" events). These supernovae obtain their name from a conspicuous signature of narrow hydrogen emission atop a weaker hydrogen P Cygni profile that together are superimposed on a loosely identifiable SS-like SN Ia spectrum (Aldering et al. 2006;Prieto et al. 2007;Leloudas et al. 2013). Apart from exhibiting similar properties to the recent PTF11kx ( §4.5.1) and SN 2005gj ( §4.5.3), Silverman et al. (2013d) find that SN Ia-CSM have a range of peak absolute magnitudes (−21.3 ≤ M R ≤ −19), are a spectroscopically homogenous class, and all reside in late-type spiral and irregular host-galaxies.
As for peculiar sub-luminous events, Narayan et al. (2011) and Foley et al. (2013) discussed the heterogeneity of the SN 2002cx-like subclass of SN Ia. Consisting of around 25 members spectroscopically similar to SN 2002cx , these new events generally have lower maximum light velocities spanning from 2000 to 8000 km s −1 and a range of peak luminosities that are typically lower than those of FAINT SN Ia (−14.2 to −18.9). In addition, this class of objects have "hot" temperature structures and−in contrast to SN Ia that follow the WLR−have low luminosities for their light curve shape. This suggests a distinct origin, such as a failed deflagration of a C+O white dwarf (Foley et al. 2009;Jordan et al. 2012;Kromer et al. 2013a) or double detonations of a sub-Chandrasekhar mass white dwarf with non-degenerate helium star companion Sim et al. 2012;Wang et al. 2013a). It is estimated that for every 100 SN Ia, there are 31 +17 −13 peculiar SN 2002cx-like objects in a given volume (Foley et al. 2013).

SN Ia Subtype Summary
In Figure 14 we plot average literature values of M B (peak), ∆m 15 (B ), and V peak (Si II λ6355) versus one another for all known SN Ia subtypes. For M B (peak) versus ∆m 15 (B ), the WLR is apparent. We have included the brightest SN 2002cx-likes (Foley et al. 2013) for reference, as these events are suspected as having separate origins from the bulk of normal SN Ia (Hillebrandt et al. 2013). We have not included Ia-CSM events given that estimates of expansion velocities and luminosities, without detailed modeling, are obscured by CSM interaction. However, it suffices to say for Figure 14 that Ia-CSM are nearest to SS and SCC SN Ia in both projected Doppler velocities and peak M R brightness (Silverman et al. 2013d). At a separate end of these SN Ia diversity plane(s), V peak (Si II λ6355) versus ∆m 15 (B ) further separates FAINT−CL SN Ia and peculiar events away from the pattern between SCC/SN 1991T-like over-luminous SN Ia and normal subtypes, where the former tend to be slowdecliners (i.e. typically brighter) with slower average velocities.  Blondin et al. (2012), and additional points discussed in Pakmor et al. (2013). Bottom: Expansion velocities at maximum light (± 3 days; via Si II λ6355 line velocities) versus ∆m15(B ). All subtypes have been tagged in accordance with the same color-scheme as in Figure 5 - Figure 13. Included for reference are the brightest, peculiar SN 2002cx-likes (light blue circles). Outliers for each subtype have been labeled for clarity and reference. We also plot mean values for the SCC and CL subtypes (larger circles), and include the mean values for SS, CN, and BL SN Ia (large diamonds) as reported by Blondin et al. (2012).
To summarize the full extent of SN Ia subtypes in terms of the qualitative luminosity and expansion velocity patterns, in Figure 15 we have outlined how SN Ia relate to one another thus far (for quantitative assessments, see Blondin et al. 2012;Silverman et al. 2012b;Folatelli et al. 2013). Broadly speaking, the red ward evolution of SN Ia features span low to high rates of decline for a large range of luminosities. Shallow Silicon and Super-Chandrasekhar Candidate SN Ia are by far the brightest, while Ia-CSM SN exhibit bright Hα emission features. These "brightest" SN Ia also show low to moderate expansion velocities andv Si II . From BL to CN to SS/SSC SN Ia, mean peak absolute brightnesses scale up with an overall decrease in maximum light line velocities. Meanwhile, CL SN Ia fall between low velocity and high velocity gradients, but lean toward HVG SN Ia in terms of their photospheric velocity evolution. Comparatively, peculiar SN 2002cx-like and other sub-luminous events are by far the largest group of thermonuclear outliers.
Obtaining observations of SN Ia that lie outside the statistical norm is important for gauging the largest degree by which SN Ia properties diverge in nature. However, just as imperative for the cause remains filling the gaps of observed SN Ia properties (e.g. color evolution) with well-observed SN Ia. This is especially true for those SN Ia most similar to one another, aka "nearest neighbors" , and transitional-type SN Ia.

Signatures of C+O Progenitor Material
If the primary star of most SN Ia is a C+O WD, and if the observed range of SN Ia properties is primarily due to variances in the ejected mass or abundances of material synthesized in the explosion (e.g., 56 Ni), then this should also be reflected in the remaining amount of carbon and oxygen if M Ch is a constant parameter (see Maeda et al. 2010a;Blondin et al. 2013;Dessart et al. 2013a). On the other hand, if one assumes that the progenitor system is the merger of two stars (Webbink 1984;Iben & Tutukov 1984;Pakmor et al. 2013;Moll et al. 2013) or a rapidly rotating WD (Hachisu et al. 2012)−both of which are effectively obscured by an amorphous region and/or disk of C+O material−then the properties of C and O absorption features will be sensitive to the interplay between ejecta and the remaining unburned envelope (see Livio & Pringle 2011).
Oxygen absorption features (unburned plus burned ejecta) are often present as O I λ7774 in the premaximum spectra of SN Ia ( Figure 5). They may exhibit similar behavior to those seen in SN 2011fe ( §4.1),

Sub-luminous
Over-luminous Spectroscopic detections of carbon-rich material have been documented since the discovery of SN 1990N (see Leibundgut et al. 1991;Jeffery et al. 1992;Branch et al. 2007b;Tanaka et al. 2008) and have been primarily detected as singly ionized in the optical spectra of LVG−CN SN Ia (Parrent et al. 2011). However, NIR spectra of some SN Ia subtypes have been suspected of harboring C I absorption features Hsiao et al. 2013, see also Marion et al. 2006Marion et al. , 2009a, while C III has been tentatively identified in "hotter" SS/SN 1991T-like SN Ia (Hatano et al. 2002;Garavini et al. 2004;Chornock et al. 2006).
Observations of the over-luminous SCC SN 2003fg suggested the presence of a larger than normal C II λ6580 absorption signature (Howell et al. 2006). Later in 2006, with the detection of a conspicuous C II λ6580 absorption "notch" in the early epoch observations of the normal SN Ia, SN 2006D, Thomas et al. (2007) reconsidered the question of whether or not spectroscopic signatures of carbon were a ubiquitous property of all or at least some SN Ia subtypes.
As follow-up investigations, Parrent et al. (2011) andFolatelli et al. (2012) presented studies of carbon features in SN Ia spectra, particularly those of C II λλ6580, 7234 (which are easier to confirm than λλ4267, 4745). However weak, conspicuous 6300Å absorption features were reported in several SN Ia spectra obtained during the pre-maximum phase. It was shown that most of the objects that exhibit clear signatures are of the LVG−CN SN Ia subtype, while HVG−BL SN Ia may either be void of conspicuous signatures due to severe line blending, or lack carbon altogether, the latter of which is consistent with DDT models (e.g., Hachinger et al. 2013) and could also be partially due to increased progenitor metallicities Meng et al. 2011;Milne et al. 2013). This requires further study and spectrum synthesis from detailed models. Thomas et al. (2011b) presented additional evidence of unburned carbon at photospheric velocities from observations of 5 SN Ia obtained by the Nearby Supernova Factory. Detections were based on the presence of relatively strong C II 6300Å absorption signatures in multiple spectra of each SN, supported by automated fitting with the SYNAPPS code (Thomas et al. 2011a). They estimated that at least 22 +10 -6 % of SN Ia exhibit spectroscopic C II signatures as late as day −5, i.e. carbon features, whether or not present in all SN Ia, are not often seen even as early as day −5.
Folatelli et al. (2012) later searched through the Carnegie Supernova Project (CSP) sample and found at least 30% of the objects show an absorption feature that can be attributed to C II λ6580.  searched for carbon in the BSNIP sample and found that ∼ 11% of the SN Ia studied show carbon absorption features, while ∼ 25% show some indication of weak 6300Å absorption. From their sample, they find that if the spectra of SN Ia are obtained before day −5, then the detection percentage is higher than ∼ 30%. Recently it has also been confirmed that "carbon-positive" SN Ia tend to have bluer near-UV colors than those without conspicuous C II λ6580 signatures (Thomas et al. 2011b;Milne et al. 2013).  estimate the range of carbon masses in normal SN Ia ejecta to be (2 − 30) x 10 −3 M . For SN 2006D, Thomas et al. 2007 estimated 0.007 M of carbon between 10,000 and 14,000 km s −1 as a lower limit. Thomas et al. 2007 also note that the most vigorous model of Röpke et al. (2006) left behind 0.085 M of carbon in the same velocity interval. However, we are not aware of any subsequent spectrum synthesis for this particular model that details the state of an associated 6300Å signature.
In the recent detailed study on SN Ia spectroscopic diversity, Blondin et al. (2012) searched for signatures of C II λ6580 in a sample of 2603 spectra of 462 nearby SN Ia and found 23 additional "carbon-positive" SN Ia. Given that seven of the nine CN SN Ia reported by Blondin et al. (2012) with spectra prior to day −10 clearly exhibit signatures of C II, and that ∼30−40% of the SN Ia within their sample are of the CN subtype, it is likely that at least 30−40% of all SN Ia leave behind some amount of carbon-rich material, spanning velocities between 8000 − 18,000 km s −1 (Parrent et al. 2011;Pereira et al. 2013;Cartier et al. 2013).
Considering the volume-limited percentage of Branch normal SN Ia estimated by Li et al. (2011c), roughly 50% or more are expected to contain detectable carbonrich material in the outermost layers. If this is true, then it implies that explosion scenarios that do not naturally leave behind at least a detectable amount (pEWs ∼ 5−25Å) of unprocessed carbon can only explain half of all SN Ia or less (sans considerations of Ia-CSM progenitors, subtype-ejecta hemisphere dualities, and effects of varying metallicities; see below).
Historically, time-evolving signatures of C II λλ6580, 7234 from the computed spectra of some detailed models have not revealed themselves to be consistent with the current interpretations of the observations. This could be due to an inadequate lower-extent of carbon within the models (Thomas et al. 2007;Parrent et al. 2011;Blondin et al. 2012) or the limits of the resolution for the computed spectra (Blondin et al. 2011).
It should be noted that 6300Å features are present in the non-LTE pre-maximum light spectra of Lentz et al. (2000) who assessed metallicity effects on the spectrum for a pure deflagration model (see their Fig. 7). Overall, Lentz et al. (2000) find that an increase in C+O layer metallicities results in a decreased flux (primarily UV) in addition to a blue ward shift of absorption minima (primarily the Si II 6100Å feature). While Lentz et al. (2000) did not discuss whether or not the weak 6300Å absorption signatures are due to C II λ6580, it is likely the case given that an increase in C+O layer metallicities is responsible for the seen decrease in the strength of the 6300Å feature. However, it should be emphasized that the "strength" of this supposed C II λ6580 feature appears to be a consequence of how Lentz et al. (2000) renormalized abundances for metallicity enhancements in the C+O layer. In other words, even though the preponderance of normal SN Ia with detectable C II λ6580 notches are of the NV subgroup, the fact that HV SN Ia are thus far void of 6300Å notches does not imply robust consistency with the idea that nearest neighbor HV SN Ia properties are solely the result of a progenitor with relatively higher metallicities (Lentz et al. 2001b). Such a claim would need to be verified by exploring a grid of models with accompanying synthetic spectra.
Additionally, carbon absorption features could signify an origin that is separate from explosion nucleosynthesis if most SN Ia are the result of a merger. For example, Moll et al. (2013) recently presented angleaveraged synthetic spectra for a few "peri-merger" detonation scenarios. In particular, they find a causal connection between "normal" C II λ6580 signatures and the secondary star for both sub-and super-Chandrasekhar mass cases (c.f. Hicken et al. 2007;Zheng et al. 2013;Dessart et al. 2013a). With constraints from UV spectra ) and high velocity features ( §3.4), peri-megers can be used to explore the expanse of their spectroscopic influence within the broader picture of SN Ia diversity.
Coincident with understanding the relevance of remaining carbon-rich material is the additional goal of grasping the spectroscopic role of species that arise from carbon-burning below the outermost layers, e.g., magnesium . While signatures of Mg II λλ4481, 7890 are frequently observed at optical wavelengths during the earliest phases prior to maximum light, these wavelength regions undergo severe line-blending compared to the NIR signatures of Mg II. Consequently, Mg II λ10927 has served as a better investment for measuring the lower regional extent and conditions during which a DDT is thought to have taken place (e.g., Rudy et al. 2002;Marion et al. 2003Marion et al. , 2006Marion et al. , 2009aHsiao et al. 2013, however see our §4.1).

High Velocity (>16,000 km s −1 ) Features
The spectra of many SN Ia have shown evidence for high-velocity absorption lines of the Ca II NIR triplet (IR3) in addition to an often concurrent signature of high-velocity Si II λ6355 (Mazzali et al. 2005a,b). Most recently, high velocity features (HVFs) have also been seen in SN 2009ig (Foley et al. 2012c), SN 2012fr (Childress et al. 2013c, and the SN 2000cx-like SN 2013bh (Silverman et al. 2013c). Overall, HVFs are more common before maximum light, display a rich diversity of behaviors (Childress et al. 2013b), tend to be concurrent with polarization signatures (Leonard et al. 2005;Tanaka et al. 2010), and may be due to an intrinsically clumpy distribution of material Kasen et al. 2003;Thomas et al. 2004;Tanaka et al. 2006;Hole et al. 2010). Maund et al. (2010b) showed that the Si II λ6355 line velocity decline rate,v Si II , is correlated with the polarization of the same line at day −5, p Si II , and is consistent with an asymmetric distribution of IMEs. This interpretation is also complimentary with a previous finding thatv Si II is correlated with v neb , the apparent Doppler line shift of [Fe II] 7155 emitted from the "core" during late times (Maeda et al. 2010b;Silverman et al. 2013a). For the recent SN 2012fr, high velocity features of Ca II IR3 and Si II λ6355 at day −11 show concurrent polarization signatures that decline in strength during post-maximum light phases (Maund et al. 2013).
As for the origin of HVFs, they may be the result of abundance and/or density enhancements due to the presence of a circumstellar medium (Gerardy et al. 2004;Quimby et al. 2006b). If abundance enhancements are responsible, it could be explained by an overabundant, outer region of Si and Ca synthesized during a pre-explosion simmering phase (see Piro 2011 andZingale et al. 2011). On the other hand, the HVFs in LVG SN Ia spectra could indicate the presence of an opaque disk. For example, it is plausible that HVFs are due to magnetically induced merger outflows (Ji et al. 2013, pending abundance calculations of a successful detonation), or interaction with a tidal tail and/or secondary star (e.g., Moll et al. 2013).
Most recently, Childress et al. (2013b) studied 58 low-z SN Ia (z < 0.03) with well-sampled light curves and spectra near maximum light in order to access potential relationships between light curve decline rates and empirical relative strength measurements of Si II and Ca II HVFs. They find a consistent agreement with Maguire et al. (2012) in that the Ca II velocity profiles assume a variety of characteristics for a given ∆m 15 (B ) solely because of the overlapping presence of HVFs. In addition, Childress et al. (2013b) show for their sample that the presence of HVFs is not strongly related to the overall intrinsic (B − V ) max colors. It is also seen that SN Ia with ∆m 15 (B ) > 1.4 continue to be void of conspicuous HVFs, while the strength of HVFs in normal SN Ia is generally larger for objects with broader light curves. Finally, and most importantly, the strength of HVFs at maximum light does not uniquely characterize HVF pre-maximum light behavior.
Notably, Silverman et al. (2013a) find no correlation between nebular velocity and ∆m 15 (B ), and for a given light-curve shape there is a large range of observed nebular velocities. Similarly Blondin et al. (2012) found no relation between the FWHM of late time 4700Å iron emission features and ∆m 15 (B ). This implies the peak brightness of these events do not translate toward uniquely specifying their late time characteristics, however the data do indicate a correlation between observed (B − V) max and this particular measure of line-of-sight nebular velocities.
We should also note that while HVG SN Ia do not clearly come with HVFs in the same sense as for LVG SN Ia, the entire 6100Å absorption feature for HVG SN Ia spans across velocity intervals for HVFs detected in LVG SN Ia. This makes it difficult to regard LVG and HVG subtypes as two separately distinct explosion scenarios. Instead we can only conclude that HVFs are a natural component of all normal SN Ia, whether conspicuously separate from a photospheric region or concealed as an extended region of absorbing material in the radial direction.

Empirical Diversity Diagnostics
The depth ratio between 5750 and 6100Å absorption features, R("Si II") (Nugent et al. 1995), has been found to correlate with components of the WLR. In addition, Benetti et al. (2005) find a rich diversity of R(Si II) pre-maximum evolution among LVG and HVG SN Ia.
As for some observables not directly related to the decline rate parameter, Patat et al. (1996) studied a small sample of well observed SN Ia and found no apparent correlation between the blue-shift of the 6100 A absorption feature at the time of maximum and ∆m 15 (B ). Similarly, Hatano et al. (2000) showed that R(Si II) does not correlate well with v 10 (Si II), the photospheric velocity derived from the Si II λ6355 Doppler line velocities 10 days after maximum light. This could arise from two or more explosion mechanisms, however Hatano et al. (2000) note that their interpretation is "rudimentary" on account of model uncertainties and the limited number of temporal datasets available at that time. In the future, it would be worthwhile to reaccess these trends with the latest detailed modeling. Hachinger et al. (2006) made empirical measurements of spectroscopic feature pEWs, flux ratios, and projected Doppler velocities for 28 well-observed SN Ia, which include LVG, HVG, and FAINT subtypes. For normal LVG SN Ia they find similar observed maximum light velocities (via Si II λ6355; ∼ 9000−10,600 km s −1 ). Meanwhile the HVG SN Ia in their sample revealed a large spread of maximum light velocities (∼ 10,300−12,500 km s −1 ), regardless of the value of ∆m 15 (B ). This overlap in maximum light velocities implies a natural continuum between LVG and HVG SN Ia (enabling unification through asymmetrical contexts; Maeda et al. 2010b;Maund et al. 2010b). They also note that FAINT SN Ia tend to show slightly smaller velocities at B -band maximum for larger values of ∆m 15 (B ), however no overreaching trend of maximum light expansion velocities from LVG to HVG to FAINT SN Ia was apparent from this particular sample of SN Ia. Hachinger et al. (2006) did find several flux ratios to correlate with ∆m 15 (B). In particular, they confirm that the flux ratio, R("S II λ5454, 5640"/"Si II λ5972"), is a fairly reliable spectroscopic luminosity indicator in addition to R(Si II). Hachinger et al. (2006) conclude that these and other flux ratio comparisons are the result of changes in relative abundances across the three main SN Ia subtypes. In a follow-up investigation, Hachinger et al. (2009) argue that the correlation with luminosity is a result of ionization balance, where dimmer objects tend to have a larger value of R(Si II). Silverman et al. 2012b later studied correlations between these and other flux ratios of SN Ia from the BSNIP sample and find evidence to suggest that CSM-associated events tend to have larger 6100Å blue-shifts in addition to broader absorption features at the time of maximum light (see also Arsenijevic 2011;Foley et al. 2012a). Altavilla et al. (2009) studied the R(Si II) ratio and expansion velocities of intermediate-redshift supernovae. They find that the comparison of intermediateredshift SN Ia spectra with high S/N spectra of nearby SN Ia do not reveal significant differences in the optical features and the expansion velocities derived from the Si II and Ca II lines that are within the range observed for nearby SN Ia. This agreement is also found in the color and decline of the light curve (see also Mohlabeng & Ralston 2013).
While the use of empirically determined singleparameter descriptions of SN Ia have proved to be useful in practice, they do not fully account for the observed diversity of SN Ia Benetti et al. 2004;Pignata et al. 2004). With regard to SN Ia diversity, it should be reemphasized that special care needs to be taken with the implementation of flux ratios and pEWs. Detailed modeling is needed when attempting to draw connections between solitary characteristics of the observed spectrum and the underlying radiative environment, where a photon-ray's route crosses many radiative contributions that form the spectrum's various shapes, from UV to IR wavelengths. For example, the relied upon 5750, 6100Å features used for R(Si II) have been shown to be influenced by more than simply Si II, as well as from more locations (and therefore various temperatures) than a single region of line formation (Bongard et al. 2008). In fact, it is likely that a number of effects are at play, e.g., line blending and phase evolution effects. Furthermore, v 10 (Si II) is a measure of the 6100Å absorption minimum during a phase of intense line blending with no less than Fe II, which imparts a bewildering array of lines throughout the optical bands (Baron et al. , 1996. Still, parameters such as R(Si II) have served as useful tools for SN Ia diversity studies in that they often correlate with luminosity  and are relatively accessible empirical measurements for large samples of under-observed SN Ia. A detailed study on the selection of global spectral indicators can be found in Bailey et al. 2009. 3.6 The Adjacent Counterparts of Optical Wavelengths 3.6.1 Ultraviolet Spectra SN Ia are known as relatively "weak emitters" at UV wavelengths (< 3500Å; Panagia 2007). It has been shown that UV flux deficits are influenced by lineblanketing effects from IPEs within the outermost layers of ejecta Hachinger et al. 2013;Mazzali et al. 2013), overall higher expansion velocities (Foley & Kasen 2011;Wang et al. 2012), progenitor metallicity Lentz et al. 2000;Sim et al. 2010b), viewing angle effects (e.g., Blondin et al. 2011), or a combination of these ). Although, it is not certain which of these play the dominant role(s) in controlling UV flux behaviors among all SN Ia. For SN 1990N and SN 1992A, two extensively studied SN Ia, pre-maximum light UV observations were made and presented by Leibundgut et al. (1991) and Kirshner et al. (1993), respectively. These observations revealed their expected sensitivity to the source temperature and opacity at UV wavelengths.
It was not until recently when a larger UV campaign of high S/N, multi-epoch spectroscopy of distant SN Ia was presented and compared to that of local SN Ia (Ellis et al. 2008, see also Milne et al. 2013). Most notably, Ellis et al. (2008) found a larger intrinsic dispersion of UV properties than could be accounted for by the span of effects seen in the latest models (e.g., metallicity of the progenitor, see Höflich et al. 1998;Lentz et al. 2000).
As a follow-up investigation, Cooke et al. (2011) utilized and presented data from the STIS spectrograph onboard Hubble Space Telescope (HST) with the intent of studying near-UV, near-maximum light spectra (day −0.32 to day +4) of nearby SN Ia. Between a high-z and low-z sample, they find a noticeable difference between the mean UV spectrum of each, suggesting that the cause may be related to different metallicities between the statistical norm of each sample. Said another way, their UV observations suggest a plausible measure of two different populations of progenitors (or constituent scenarios) that could also be dependent on the metallically thereof, including potentially larger dependencies such as variable 56 Ni mass and line blanketing due to enhanced burning within the outermost layers (Marion et al. 2013). It should be noted, despite the phase selection criterion invoked by Cooke et al. (2011), it may not be enough to simply designate a phase range in order to avoid phase evolution effects (see Fig. 7 of Childress et al. 2013b).
In order to confirm spectroscopic trends at UV wavelengths, a better method of selection will be necessary as the largest UV difference found by Cooke et al. (2011) and Maguire et al. (2012) between the samples overlaps with the Si II, Ca II H&K absorption features (3600−3900Å), i.e. a highly blended feature that is far too often a poorly understood SN Ia variable, both observationally (across subtype and phase) and theoretically, within the context of line formation at UV−NIR wavelengths (Mazzali 2000;Kasen et al. 2003;Thomas et al. 2004;Foley 2012;Marion et al. 2013;Childress et al. 2013b). While it is true that different radiative processes dominate within different wavelength regions, there are a multitude of explanations for such a difference between the Si II−Ca II blend near 3700Å. Furthermore, the STIS UV spectra do not offer a look at either the state of the 6100Å absorption feature (is it completely photospheric?−the answer requires spectrum synthesis even for maximum light phases), nor is it clear if the same is true for Ca II in the NIR where high-velocity components thereof are most easily discernible .
It is important to further reemphasize that the timedependent behavior of the sum total of radiative processes that generate a spectrum from a potentially axially-asymmetric (and as of yet unknown) progenitor system and explosion mechanism are not well understood, much less easily decipherable with an only recently obtained continuous dataset for how the spectrum itself evolves over time at optical wavelengths 20 . Which is only to say, given the current lack of certain predictability between particular observational characteristics of SN Ia (e.g., spectroscopic phase transition times), time series observations at UV wavelengths would offer a beneficial route for the essential purposes of hand-selecting the 'best' spectrum comparisons in order to ensure a complete lack of phase evolution effects.
Recently Rather, the inclusion of asymmetry into a standard model picture of SN Ia should be a relevant part of their observed diversity (e.g., Kasen et al. 2009;Blondin et al. 2011).
More recently, Mazzali et al. (2013) obtained 10 HST UV−NIR spectra of SN 2011fe, spanning −13.5 to +41 days relative to B -band maximum. They analyzed the data along side spectrum synthesis results from three explosion models, namely a 'fast deflagration' (W7), a low-energy delayed-detonation (WS15DD1; Iwamoto et al. 1999), and a third model treated as an intermediary between the outer-layer density profiles of the other two models ("W7 + "). From the seen discrepancies between W7 and WS15DD1 during the early pre-maximum phase, in addition to optical flux excess for W7 and a mismatch in observed velocities for WS15DD1, Mazzali et al. (2013) conclude that their modified W7 + model is able to provide a better fit to the data because of the inclusion of a high velocity tail of low density material. In addition, and based on a spectroscopic rise time of ∼ 19 days, Mazzali et al. (2013) infer a ∼ 1.4 day period of optical quiescence after the explosion (see Piro & Nakar 2013;Chomiuk 2013).

Infrared Light Curves and Spectra
By comparing absolute magnitudes at maximum of two dozen SN Ia, Krisciunas (2005) argue that SN Ia can be best used as standard candles at NIR wavelengths (which was also suggested by Elias et al. 1985a,b), even without correction for optical light curve shape. Wood-Vasey et al. (2008) later confirmed this to be the case from the analysis of 1087 near-IR (JHK) measurements of 21 SN Ia. Based on their data and data from the literature, they derive absolute magnitudes of 41 SN Ia in the H -band with rms scatter of 0.16 magnitudes. Folatelli et al. (2010) find a weak dependence of J -band luminosities on the decline rate from 9 NIR datasets, in addition to V −J corrected J -band magnitudes with a dispersion of 0.12 magnitudes. Mandel et al. (2011) constructed a statistical model for SN Ia light curves across optical and NIR passbands and find that near-IR luminosities enable the most ideal use of SN Ia as standard candles, and are less sensitive to dust extinction as well. Kattner et al. (2012) analyzed the standardizability of SN Ia in the near-IR by investigating the correlation between observed peak near-IR absolute magnitude and post-maximum ∆m 15 (B ). They confirm that there is a bimodal distribution in the near-IR absolute magnitudes of fast-declining SN Ia (Krisciunas et al. 2009) and suggest that applying a correction to SN Ia peak luminosities for decline rate is likely to be beneficial in the J and H bands, making SN Ia more precise distance indicators in the IR than at optical wavelengths (Barone- Nugent et al. 2012).
While optical spectra of SN Ia have received a great deal of attention in the recent past, infrared datasets (e.g., Kirshner et al. 1973b;Meikle et al. 1996;Bowers et al. 1997;Rudy et al. 2002;Höflich et al. 2002) are either not obtained, or are not observed at the same epochs or rate as their optical counterparts. This has only recently begun to change. Thus far, the largest NIR datasets can be found in Marion et al. (2003) and Marion et al. (2009a). Marion et al. (2003) obtained NIR spectra (0.8−2.5 µm) of 12 normal SN Ia, with fairly early coverage. Later, Marion et al. (2009a) presented and studied a catalogue of NIR spectra (0.7 − 2.5 µm) of 41 additional SN Ia. In all, they report an absence of conspicuous signatures of hydrogen and helium in the spectra, and no indications of carbon via C I λ10693 (however, see our §4.1). For an extensive review on IR observations, we refer the reader to Phillips (2012).

Drawing Conclusions about SN Ia Diversity from SN Ia Rates Studies
It has long been perceived that a supernova's local environment, rate of occurrence, and host galaxy properties (e.g., WD population) serve as powerful tools for uncovering solutions to SN Ia origins (Zwicky 1961;Hamuy et al. 1995 Despite this broad extent of the progenitor problem, measurements of the total cosmic SN Ia rate, R SN Ia (z), can be made to gauge the general underlying behavior of actively contributing systems ). Further insight into how various progenitor populations impart their signature onto R SN Ia (z) comes about by considering which scenarios lead to a "prompt" (or a "tardy") stellar demise (Scannapieco & Bildsten 2005;Mannucci 2005). Whether or not mergers involve both a ("prompt") helium-burning or ("tardy") degenerate secondary star remains to be seen (Woods et al. 2011;Hillebrandt et al. 2013;Dan et al. 2013 and references therein). Because brighter SN Ia prefer younger, metalpoor galaxies, and a linear relation exists between the SN Ia light curve shape and gas-phase metallicity, the principle finding has been that the rate of the universally prompt component is proportional to the star formation rate of the host galaxy, whereas the second delayed component's rate is proportional to the stellar mass of the galaxy (Sullivan et al. 2006;Howell et al. 2007;Sullivan et al. 2010;Zhang 2011;Pan et al. 2013). The SN Ia galaxy morphology study of Hicken et al. (2009a) has since progressed this discussion of linking certain observed SN Ia properties with their individual environments. Overall, the trend of brighter/dimmer SN Ia found in younger/older hosts remains, however now with indications that a continuous distribution of select SN Ia subtypes exist in multiple host galaxy morphologies and projected distances within each host.
To understand the full form of R SN Ia (z), taking into account the delay time distribution (DTD) for every candidate SN Ia system is necessary (see Bonaparte et al. 2013;Claeys et al. 2014).  find that the DTD peaks prior to 2.2 Gyr and has a long tail out to ∼ 10 Gyr. They conclude that a DTD with a power-law t −1.2 starting at time t = 400 Myr to a Hubble time can satisfy both constraints of observed cluster SN rates and iron-to-stellar mass ratios, implying that that half to a majority of all SN Ia events occur within one Gyr of star formation (see also Strolger et al. 2010;Meng et al. 2011).
In general, the DTD may be the result of binary mergers (Ruiter et al. 2009;Toonen et al. 2012;Nelemans et al. 2013) and/or a single-degenerate scenario (Hachisu et al. 2008(Hachisu et al. , 2012Chen et al. 2013), but with the consideration that evidence for delay times as short as 100 Myr have been inferred from SN remnants in the Magellanic Clouds (Badenes et al. 2009;Maoz & Badenes 2010). From a recent comparison of low/highz SN Ia rate measurements and DTDs of various binary population synthesis models, Graur et al. (2013) argue that single-degenerate systems are ruled out between 1.8 < z < 2.4. Overall, their results support the existence of a double-degenerate progenitor channel for SN Ia if the the number of double-degenerate systems predicted by binary population synthesis models can be "aptly" increased .
However, initial studies have primarily focussed upon deriving the DTD without taking into account the possible effects of stellar metallicity on the SN Ia rate in a given galaxy. Given that lower metallicity stars leave behind higher mass WD stars (Umeda et al. 1999b;Timmes et al. 2003), Kistler et al. (2013) and Meng et al. (2011) argue that the effects of metallicity may serve to significantly alter the SN Ia rate (see also Pan et al. 2013). In fact, models that include the effects of metallicity (e.g., Kistler et al. 2013) find similar consistencies with the observed R SN Ia (z). Notably, recent spectroscopic studies do indicate a stronger preference of low-metallicity hosts for super-Chandrasekhar candidate SN Ia Childress et al. 2011), which may just as well be explained by low metallicity single-degenerate systems (Hachisu et al. 2012). While there are not enough close binary WD systems in our own galaxy that would result in SCC DD scenarios (Parthasarathy et al. 2007), sub-Chandrasekhar merging binaries may be able to account for discrepancies in the observed rate of SN Ia (Badenes & Maoz 2012;Kromer et al. 2013b).
Although, we wish to remind the reader that since spectrophotometry of SN Ia so far offer the best visual insight into these distant extragalactic events, and because there is no clear consensus on the origin of their observed spectrophotometric diversity, there is no clear certainty as to what distribution of progenitor scenarios connect with any kind of SN Ia since none have been observed prior to the explosion. Furthermore, whether or not brighter or dimmer SN Ia "tend to" correlate with any property of their hosts does not alleviate the discussion down to one or two progenitor systems (e.g., singleversus double-degenerate systems) since the most often used tool for probing SN Ia diversity over all distance scales, i.e. the "stretch" of a light curve, does not necessarily uniquely determine the spectroscopic subtype. Rather, such correlations reveal the degree of an underlying effect from samples of uncertain and unknown SN Ia subtype biases, i.e. dust extinction in star formation galaxies and progenitor ages also evolve along galaxy mass sequences (Childress et al. 2013a) and the redshift-color evolution of SN Ia remains an open issue (Mohlabeng & Ralston 2013;Pan et al. 2013;Wang et al. 2013b).
While it is important to consider the full redshift range over which various hierarchies of progenitor and subtype sequences may dominate over others, such studies are rarely able to incorporate spectroscopic diversity as input (a "serendipitous" counter-example being Krughoff et al. 2011). This is relevant given that the landscape of SN Ia spectroscopic diversity has not yet been seen to be void of line-of-sight discrepancies for all progenitor scenarios (particularly so for double degenerate detonations/mergers, e.g., Shen et al. 2013;Pakmor et al. 2013;Moll et al. 2013;. Ultimately, robust theories should be able to connect spectroscopic subtypes with individual or dual instances of particular progenitor systems, which requires detailed spectroscopic modeling.
Thus, the consensus as to how many progenitor channels contribute to SN Ia populations is still unclear.
Broadly speaking, there are likely to be no less than two to three progenitor scenarios for normal SN Ia so long as single-degenerate systems remain viable (Hachisu et al. 2012), if not restricted to explaining Ia-CSM SN alone (see Han & Podsiadlowski 2006;Silverman et al. 2013d;Leloudas et al. 2013). Given also a low observed frequency of massive white dwarfs and massive doubledegenerate binaries near the critical mass limit with orbital periods short enough to merge within a Hubble time, some normal SN Ia are still perceived as originating from single-degenerate systems (Parthasarathy et al. 2007). Meanwhile, some portion of events may also be the result of a core-degenerate merger (Soker et al. 2013), while some merger phenomena are possibly accelerated within triple systems (Thompson 2011;Kushnir et al. 2013;Dong et al. 2014). It likewise remains unclear whether or not some double-degenerate mergers predominately result in the production of a neutron star instead of a SN Ia (Saio & Nomoto 1985;Nomoto & Kondo 1991;Piersanti et al. 2003;Saio & Nomoto 2004;Dan et al. 2013;Tauris et al. 2013). At present, separately distinct origins for spectroscopically similar SN Ia cannot be ruled out by even one discovery of a progenitor system; the spectroscopic diversity is currently too great and too poorly understood to confirm without greater unanimity among explosion models and uniformity in data collection efforts.

Some Recent SN Ia
During the past decade, several normal, interesting, and peculiar SN Ia have been discovered. For example, the recent SN 2009ig, 2011fe, and 2012fr are nearby SN Ia that were discovered extremely young with respect to the onset of the explosion (Nugent et al. 2011;Foley et al. 2012c;Childress et al. 2013c) and have been extensively studied at all wavelengths, yielding a clearer understanding of the time-dependent behavior of SN spectroscopic observations, in addition to a better context by which to compare. Below we briefly summarize some of the highlighted discoveries during the most recent decade, during which it has revealed a greater diversity of SN Ia than was previously known. In the appendix we provide a guide to the recent literature of other noteworthy SN Ia discoveries. We emphasize that these sections are not meant to replace reading the original publications, and are only summarized here as a navigation tool for the reader to investigate further.

SN 2011fe in M101
Thus far, the closest spectroscopically normal SN Ia in the past 25 years, SN 2011fe (PTF11kly), has provided a great amount of advances, including testing SN Ia distance measurement methods (Matheson et al. 2012;Vinkó et al. 2012;Lee & Jang 2012). For example, the early spectroscopy of SN 2011fe showed a clear and certain time-evolving signature of high-velocity oxygen that varied on time scales of hours, indicating sizable overlap between C+O, Si, and Ca-rich material and newly synthesized IMEs within the outermost layers (Nugent et al. 2011). Parrent et al. (2012) carried out analysis of 18 spectra of SN 2011fe during its first month. Consequently, they were able to follow the evolution of C II λ6580 absorption features from near the onset of the explosion until they diminished after maximum light, providing strong evidence for overlapping regions of burned and unburned material between ejection velocities of at least 10,000 and 16,000 km s −1 . At the same time, the evolution of a 7400Å absorption feature experienced a declining Doppler-shift until 5 days post-maximum light, with O I λ7774 line velocities ranging 11,500 to 21,000 km s −1 (Nugent et al. 2011). Parrent et al. (2012) concluded that incomplete burning (in addition to progenitor scenarios) is a relevant source of spectroscopic diversity among SN Ia (Tanaka et al. 2008;Maeda et al. 2010a). Pereira et al. (2013) presented high quality spectrophotometric observations of SN 2011fe, which span from day −15 to day +97, and discussed comparisons to other observations made by Brown et al. (2012), Richmond & Smith (2012), Vinkó et al. (2012), and Munari et al. (2013). From an observed peak bolometric luminosity of 1.17 ± 0.04 x 10 43 erg s −1 , they estimate SN 2011fe to have produced between ∼ 0.44 ± 0.08 − 0.53 ± 0.11 M of 56 Ni.
By contrast, Pastorello et al. (2007a) and Wang et al. (2009b) estimate 56 Ni production for the normal SN 2005cf (A.26) to be ∼ 0.7 M . It is also interesting to note that SN 2011fe and the fast-declining SN 2004eo produced similar amount of radioactive nickel, however lower for SN 2004eo (∼ 0.4 M ; Mazzali et al. 2008). Pereira et al. (2013) also made comparisons between SN 2011fe, a SNFactory normal SN Ia (SNF20080514-002) and the broad-lined HV-CN SN 2009ig (Foley et al. 2012c). Pereira et al. (2013) note similarities (sans the UV) and notable contrast with respect to high-velocity features, respectively. Pereira et al. (2013) calculatedv Si II for SN 2011fe to be ∼ 60 (± 3) km s −1 day −1 , near the high end of low-velocity gradient SN Ia events (see Benetti et al. 2005;Blondin et al. 2012). Given their high S/N, time series dataset, Pereira et al. (2013) were also able to place tighter constraints on the velocities over which C II λ6580 is observed to be present in SN 2011fe.
They conclude that C II is present down to at least as low as 8000 km s −1 , which is 2000 km s −1 lower than that estimated by Parrent et al. (2012), and is also ∼ 4000−6000 km s −1 (or more) lower than what is predicted by some past and presently favored SN Ia abundance models (e.g., W7; Nomoto et al. 1984, and the delayed detonations of Röpke et al. 2012). Hsiao et al. (2013) presented and discussed NIR time series spectra of SN 2011fe that span between day −15 and day +17. In particular, they report a detection of C I λ10693 on the blue-most side of a blended Mg II λ10927 absorption feature at roughly the same velocities and epochs as C II λ6580 found by Parrent et al. (2011) and Pereira et al. (2013), which itself is blended on its blue-most side with Si II λ6355. While searches and studies of C I λ10693 are extremely useful for understanding the significance of C-rich material from normal to cooler sub-luminous SN Ia within the greater context of all C I, C II, C III, O I absorption features (C III for "hotter" SN 1991T-likes), blended C I λ10693 absorption shoulders are certainly no more (nor no less) useful for probing lower velocity boundaries than C II λ6580 absorption notches. This is especially true given that C I λ10693 absorption features are blended from the red-most side (lower velocities) by the neighboring Mg II line, which will only serve to obscure the lower velocity information of the C I profile for the non-extreme cases (e.g., SN 1999by, Höflich et al. 2002. Hsiao et al. (2013) used the observed temporal behavior, and later velocity-plateau, of Mg II λ10927 to estimate a lower extent of ∼ 11,200 km s −1 for carbonburning products within SN 2011fe. Given that this in contrast to the refined lower extent of C II at ∼ 8000 km s −1 by Pereira et al. (2013), this could imply (i.e. assuming negligible temperature differences and/or non-LTE effects) that either some unburned material has been churned below the boundary of carbon-burning products via turbulent instabilities Gamezo et al. , 2004 and/or the distribution of emitting and absorbing carbon-rich material is truly globally lopsided (Kasen et al. 2009;Maeda et al. 2010b;Blondin et al. 2011), and may indicate the remains of a degenerate secondary star .
Detailed studies of this nearby, normal, and unreddenned SN 2011fe have given strong support for doubledegenerate scenarios (assuming low environmental abundances of hydrogen) and have placed strong constraint on single-degenerate scenarios, i.e. MS and RG companion stars have been strongly constrained for SN 2011fe (see Shappee et al. 2013 Bloom et al. (2012) confirm that the primary star was a compact star (R * 0.1 R , c.f., Bloom et al. 2012;Piro & Nakar 2013;Chomiuk 2013). From the lack of evidence for an early shock outbreak (Kasen 2010;Nakar & Sari 2012), non-detections of radio and X-ray emissions (Horesh et al. 2012;Chomiuk et al. 2012;Margutti et al. 2012), non-detections of narrow Ca II H&K or Na D lines or pre-existing dust that could be associated with the event (Patat et al. 2013;Johansson et al. 2013), and low upper-limits on hydrogen-rich gas (Lundqvist et al. 2013), the paucity of evidence for an environment dusted in CSM from a non-degenerate secondary strongly supports the double degenerate scenario for SN 2011fe. Plus, this inferred ambient environment is consistent with that of recent merger simulations (Dan et al. 2012), and could signal an avenue of interpretation for signatures of carbon-rich material as well (Branch et al. 2005;Dan et al. 2013;Moll et al. 2013). Specifically, the remaining amount of carbonrich material predicted by some explosion models may already be accounted for, and more so than would be required by the existence of low velocity detections of C I and C II. If this turns out to be the case, spectroscopic signatures of both C and O could tap into understanding (i) the sizes of merger C+O common envelopes, (ii) potential downward mixing effects between the envelope and the underlying ejecta, and/or (iii) test theories on possible asymmetries of C+O material within the post-explosion ejecta of the primary and secondary stars (Livio & Pringle 2011), which is expected to depend on the degree of coalescence .
Of course, this all rests on the assumptions that (i) the surrounding environment of a single degenerate scenario just prior to the explosion ought to be contaminated with some amount of CSM, above which it would be detected (Justham 2011;Brown et al. 2012), and (ii) the surrounding environment of a merger remains relatively "clean" (Shen et al. 2013;. In this instance, and assuming similarly above that current DDT-like models roughly fit the outcome of the explosion, absorption signatures of C (+ HV O I) may point to super-massive single-degenerate progenitors with variable enclosed envelopes and/or disks of material (e.g., Yoon & Langer 2004Hachisu et al. 2012;Scalzo et al. 2012;Tornambé & Piersanti 2013;Dan et al. 2013) or sub-Chandrasekhar mass "perimergers" for resolve (see Moll et al. 2013 and references therein).

SN 2009ig in NGC 1015
Foley et al. (2012c) obtained well-sampled, early UV and optical spectra of SN 2009ig as it was discovered 17 hr after the event (Kleiser et al. 2009;Navasardyan et al. 2009). SN 2009ig is found to be a normal SN Ia, rising to B -band maximum in ∼ 17.3 days. From the earliest spectra, Foley et al. (2012c) find Si II λ6355 line velocities around 23,000 km s −1 , which is exceptionally high for such a spectroscopically normal SN Ia (see also Blondin et al. 2012). SN 2009ig possess either an overall shallower density profile than other CN SN Ia, or a buildup of IMEs is present at high velocities. Marion et al. (2013) recently analyzed the photospheric to post-maximum light phase spectra of SN 2009ig, arguing for the presence of additional highvelocity absorption signatures from not only Si II, Ca II, but also Si III, S II and Fe II. Whether or not two separate but compositionally equal regions of line formation is a ubiquitous property of similar SN Ia remains to be seen. However, it should not be unlikely for primordial amounts of said atomic species to be present (in addition to singly-ionized silicon and calcium) on account of possible density and/or abundance enhancements within the outermost layers (Thomas et al. 2004;Mazzali et al. 2005b,a). For example, simmering effects during convective phases prior to the explosion may be responsible for dredging up IMEs later seen as HVFs, which would give favorability to single-degenerate progenitor scenarios (see Piro 2011;Zingale et al. 2011). Similarly, it is worthwhile to access the versatility of mergers in producing high-velocity features.

SN 2012fr in NGC 1365
Childress et al. (2013c) report on their time series spectroscopic observations of SN 2012fr (Klotz et al. 2012;Childress et al. 2012;Buil 2012), complete with 65 spectra that cover between ∼ 15 days before and 40 days after it reached a peak B -band brightness of −19.3. In addition to the simultaneous spectropolarimetric observations of Maund et al. (2013), the early to maximum light phase spectra of SN 2012fr reveal one of the clearest indications that SN Ia of similar type (e.g., SN 1994D, 2001el, 2009ig, 2011fe, and many others; Mazzali05a) tend to have two distinctly separate regions of Si-, Ca-based material that differ by a range of separation velocities (Childress et al. 2013b). Childress et al. (2013c) and Maund et al. (2013) discussed the various interpretations that have been presented in the past, however no firm conclusions on the origin of HVFs could be realized given the uncertainties of current explosion models. Despite this, the most recent advance toward understanding HVFs is the continual detection of polarization signatures due to the high-velocity Si II and Ca II absorption features, indicating a departure from a radially stratified, spherically symmetric geometry at some layer near or above the "photospheric region" of IMEs.

Over-luminous SN 2003fg (SNLS-03D3bb)
SN 2003fg was discovered as part of the Supernova Legacy Survey (SNLS); z = 0.2440 (Howell et al. 2006). Its peak absolute magnitude was estimated to be −19.94 in V -band, placing SN 2003fg completely outside the M V -distribution of normal low-z SN Ia (2.2 times brighter). Assuming Arnett's rule, such a high luminosity corresponds to ∼ 1.3 M of 56 Ni, which would be in conflict with SN 2003fg's spectra since only ∼ 60% of a Chandrasekhar pure detonation ends up as radioactive nickel (Steinmetz et al. 1992, however see also Pfannes et al. 2010). Given also the lower mean expansion velocities, this builds upon the picture of a super-Chandrasekhar mass progenitor for SN 2003fg and others like it (Howell et al. 2006;Jeffery et al. 2006). Yoon & Langer (2005) proposed the formation of super-Chandrasekhar mass WD stars as a result of rapid rotation. Pfannes et al. (2010) later reworked these models and found that the "prompt" detonation of a super-Chandrasekhar mass WD produces enough nickel, as well as a remainder of IMEs in the outer layers (in contrast to Steinmetz et al. 1992), to explain over-luminous SN Ia. Hachisu et al. (2012) added to this model by taking into account processes of binary evolution. Namely, with the inclusion of mass-striping, optically thick winds of a differentially rotating primary star, Hachisu et al. (2012) find three critical mass ranges that are each separated according to the spin-down time of the accreting WD. All three of these single-degenerate scenarios may explain a majority of events from sub-luminous to over-luminous SN Ia. So far no super-Chandrasekhar mass WD stars that would result in a SN Ia have been found in the sample of known WD stars in our Galaxy (Saffer et al. 1998, see also Kilic et al. 2012). However, this does not so much rule out super-Chandrasekhar mass models as it suggests that these systems are rare in the immediate vicinity within our own galaxy. Hillebrandt et al. (2007) proposed an alternative scenario involving only a Chandrasekhar-mass WD progenitor to explain the SN 2003fg event. They demon-strated that an off-center explosion of a Chandrasekharmass WD could explain the super-bright SN Ia. However, in this off-center explosion model it is not easy to account for the high Ni mass in the outer layers, in addition to the special viewing direction. Yamanaka et al. (2009a) presented early phase optical and NIR observations for SN 2009dc (Puckett et al. 2009;Harutyunyan et al. 2009;Marion et al. 2009b;Nicolas & Prosperi 2009). From the peak V -band absolute magnitude they conclude that SN 2009dc belongs to the most luminous class of SN Ia (∆m 15 (B ) = 0.65), and estimate the 56 Ni mass to be 1.2 to 1.6 M . Based on the JHK photometry Yamanaka et al. (2009a) also find SN 2009dc had an unusually high NIR luminosity with enhanced fading after ∼ day +200 (Maeda et al. 2009;Silverman et al. 2011;Taubenberger et al. 2011). The spectra of SN 2009dc also show strong, long lasting 6300Å absorption features (until ∼ two weeks postmaximum light) Based on the observed spectropolarimetric indicators, in combination with photometric and spectroscopic properties, Tanaka et al. (2010) similarly conclude that the progenitor mass of SN 2009dc was of super-Chandrasekhar origin and that the explosion geometry was globally spherically symmetric, with a clumpy distribution of IMEs. Silverman et al. (2011) presented an analysis of 14 months of observations of SN 2009dc and estimate a rise-time of ∼ 23 days and ∆m 15 (B ) = 0.72. They find a lower limit of the peak bolometric luminosity ∼ 2.4 x 10 43 erg s −1 and caution that the actual value is likely almost 40% larger. Based on the high luminosity and low mean expansion velocities of SN 2009dc, Silverman et al. (2011 derive a mass of more than 2M for the white dwarf progenitor and a 56 Ni mass of ∼ 1.4 to 1.7 M . Taubenberger et al. (2011) find the minimum 56 Ni mass to be 1.8 M , assuming the smallest possible rise-time of 22 days, and the ejecta mass to be 2.8 M . Taubenberger et al. (2013) compared photometric and spectroscopic observations of normal and SCC SN Ia at late epochs, including SN 2009dc, and find a large diversity of properties spanning through normal, SS, and SCC SN Ia. In particular the decline in the light curve "radioactive tail" for SCC SN Ia is larger than normal, along with weaker than normal [Fe III] emission in the nebular phase spectra. Taubenberger et al. (2013) argue that the weak [Fe III] emission is indicative of an ejecta environment with higher than normal densities. Previously, Hachinger et al. (2012) carried out spectroscopic modeling for SN 2009dc and discussed the model alternatives, such as a 2 M rotating WD, a core-collapse SN, and a CSM interaction scenario. Overall, Hachinger et al. (2012) found the interaction scenario to be the most promising in that it does not require the progenitor to be supermassive. Taubenberger et al. (2013) furthered this discussion in conjunction with their late time comparisons and conclude that the models of Hachinger et al. (2012) do not simultaneously match the peak brightness and decline of SN 2009dc (see also Yamanaka et al. 2013). Following the interaction scenario of Hachinger et al. (2012), Taubenberger et al. (2013) propose a non-violent merger model that produces ∼1M of 56 Ni and is enshrouded by ∼0.6−0.7M of C+O-rich material. In order to reconcile the low 56 Ni production, Taubenberger et al. (2013) note that additional luminosity from interaction with CSM is required during the first two months post-explosion. Further support for CSM interaction comes from the observed suppression of the double peak in the I -band, which is thought to arise from a breaking of ejecta stratification in the outermost layers (Kasen 2006;Kamiya et al. 2012).

Over-luminous SN 2009dc in UGC 10064
It is not yet clear if SN 2003fg, 2006gz (A.32), 2007if (A.34), and SN 2009dc are the result of a single super-Chandrasekhar mass WD star, given that even in our galaxy there is no observational evidence for the existence of such a system. Likewise, there is no direct observational evidence for the presence of very rapidly rotating massive WD stars, either single WDs or in binary systems as well. In fact, no double-degenerate close binary systems with a total mass amounting to super-Chandrasekhar mass configurations that can merge in Hubble-time have been found (Parthasarathy et al. 2007). Therefore, our current understanding of the origin of over-luminous SN Ia is limited, and more observations are needed. For example, progress has been made with the recent discovery of 24 merging WD systems via the extremely low mass Survey (see Kilic et al. 2012 and references therein), however it is unclear if any are systems that would produce a normal SN Ia. They obtained photometric and spectroscopic observations which revealed it to be unique among all observed SN Ia. Li et al. (2003) described SN 2002cx as having SN 1991T-like pre-maximum spectrum, a SN 1991bglike luminosity, and expansion velocities roughly half those of normal SN Ia.

The Peculiar SN 2002cx-like Class of SN
Photometrically, SN 2002cx has a broad peak in the R-band, a plateau phase in the I -band, and a slow late time decline. The B − V color evolution are described as nearly normal, while the V − R and V − I colors are redder than normal. Spectra of SN 2002cx during early phases evolve rapidly and are dominated by lines from IMEs and IPEs, but the features are weak overall. In addition, emission lines are present around 7000Å during post-maximum light phases, while the late time nebular spectrum shows narrow lines of iron and cobalt. Jha et al. (2006a) presented late time spectroscopy of SN 2002cx, which includes spectra at 227 and 277 days post-maximum light. They considered it as a prototype of a new subclass of SN Ia. The spectra do not appear to be dominated by the forbidden emission lines of iron, which is not expected during the "nebular phase," where instead they find a number of P Cygni profiles of Fe II at exceptionally low expansion velocities of ∼ 700 km s −1 (Branch et al. 2004a). A tentative identification of O I λ7774 is also reported for SN 2002cx, suggesting the presence of oxygen-rich material. Currently, it is difficult to explain all the observed photometric and spectroscopic properties of SN 2002cx using the standard SN Ia models (see Foley et al. 2013). However, the spectral characteristics of SN 2002cx support pure deflagration or failed-detonation models that leave behind a bound remnant instead of delayed detonations (Jordan et al. 2012;Kromer et al. 2013a;Hillebrandt et al. 2013). Phillips et al. (2007) presented extensive multi-color photometry and optical spectroscopy of SN 2005hk (Quimby et al. 2005). Sahu et al. (2008) also studied the spectrophotometric evolution SN 2005hk, covering pre-maximum phase to around 400 days after the event. These datasets reveal that SN 2005hk is nearly identical in its observed properties to SN 2002cx. Both supernovae exhibited high ionization SN 1991T-like premaximum light spectra but with low peak luminosities like that of SN 1991bg. The spectra reveal that SN 2005hk, like SN 2002cx, has expansion velocities that are roughly half those of typical SN Ia.

SN 2005hk in UGC 00272
The R and I -band light curves of both supernovae are also peculiar for not displaying the secondary maximum observed for normal SN Ia. Phillips et al. (2007) constructed a bolometric light curve from 15 days before to 60 days after B -band maximum. They conclude that the shape and exceptionally low peak luminosity of the bolometric light curve, low expansion velocities, and absence of a secondary maximum in the NIR light curves are in reasonable agreement with model calculations of a three-dimensional deflagration that produces 0.25 M of 56 Ni. Note however that the low amount of continuum polarization observed for SN 2005hk (∼ 0.2%−0.4%) is far too similar to that of more normal SN Ia to serve as an explanation for the spectroscopic peculiarity of SN 2005hk, and possibly other SN 2002cx-like events (Chornock et al. 2006;Maund et al. 2010a).

SN 2009ku
SN 2009ku was discovered by Pan-STARS-1 as a SN Ia belonging to the peculiar SN 2002cx class. Narayan et al. (2011) studied SN 2009ku and find that while its multi-band light curves are similar to that of SN 2002cx, they are slightly broader and have a later rise to gband maximum. Its peak brightness was found to be M V = −18.4 and the ejecta velocity at 18 days after maximum brightness was found to be ∼ 2000 km s −1 . Spectroscopically, SN 2009ku is similar to SN 2008ha (A.37). Narayan et al. (2011) note that the high luminosity and low ejecta velocity for SN 2009ku is not in agreement with the trend seen for SN 2002cx class of SN Ia. The spectroscopic and photometric characteristics of SN 2009ku indicate that the SN 2002cx class of SN Ia are not homogeneous, and that the SN 2002cx class of events may have a significant dispersion in their progenitor population and/or explosion physics (see also Kasliwal et al. 2012 for differences between this class and sub-luminous "calcium-rich" transients). PTF11kx. Using time series, high-resolution optical spectra, they find direct evidence supporting a singledegenerate progenitor system based on several narrow, temporal (∼ 65 km s −1 ) spectroscopic features of the hydrogen Balmer series, He I, Na I, Ti II, and Fe II. In addition, and for the first time, PTF11kx observations reveal strong, narrow, highly time-dependent Ca II absorption features that change from saturated absorption signatures to emission lines within ∼ 40 days. Dilday et al. (2012) considered the details of these observations and concluded that the complex CSM environment that enshrouds PTF11kx is strongly indicative of mass loss or "outflows," prior to the onset of the explosion of the progenitor system. Other SN Ia have exhibited narrow, temporal Na D lines before (e.g., SN 2006X, 2007le;see Simon et al. 2009;Patat et al. 2009Patat et al. , 2010Patat et al. , 2011Sternberg et al. 2011), but none have been reported as having signatures of these particular ions, which are clearly present in the high-resolution spectra of PTF11kx. On the whole, and during the earliest epochs, Dilday et al. (2012) find that the underlying SN Ia spectroscopic component of PTF11kx most closely resembles that of SN 1991T (Filippenko et al. 1992a;Gómez & López 1998) and 1999aa (Garavini et al. 2004. As for the late time phases, Silverman et al. (2013b) studied spectroscopic observations of PTF11kx from 124 to 680 days post-maximum light and find that its nebular phase spectra are markedly different from those of normal SN Ia. Specifically, the late time spectra of PTF11kx are void of the strong cobalt and iron emission features typically seen in other SN 1991T/1999aa-like and normal SN Ia events (e.g., Ruiz-Lapuente & Lucy 1992;Salvo et al. 2001;Branch et al. 2003;Stehle et al. 2005;Kotak et al. 2005;McClelland et al. 2013;Silverman et al. 2013a). For the most part, the late time spectra of PTF11kx are seen to be dominated by broad (FWHM ∼ 2000 km s −1 ) Hα emission and strong Ca II emission features that are superimposed onto a relatively blue, overly luminous continuum level that may be serving to wash out the underlying SN Ia spectroscopic information. Silverman et al. (2013b) note that the Hα emission increases in strength for ∼1 yr before decreasing. In addition, from the absence of strong Hβ, He I, and O I emission, as well as a larger than normal late time luminosity, Silverman et al. (2013b) conclude that PTF11kx indeed interacted with some form of CSM material; possibly of multiply thin shells, shocked into radiative modes of collisional excitation as the SN ejecta overtakes the slower-moving CSM. However, it should be noted that it is not yet clear if the CSM originates from a single-degenerate scenario or a H-rich layer of material that is ejected prior to a double-degenerate merger event (Shen et al. 2013, see also Soker et al. 2013).

SN 2002ic
Hamuy et al. (2003) Hamuy et al. (2003) argue that the broad component arose from ejecta−CSM interaction. By day +48, they find that the spectrum is similar to that of SN 1990N. Hamuy et al. (2003) argue that the progenitor system contained a massive AGB star, associated with a few solar masses of hydrogen-rich CSM. Kotak et al. (2004) obtained the first high resolution, high S/N spectrum of SN 2002ic. The resolved Hα line has a P Cygni-type profile, indicating the presence of a dense, slow-moving outflow (about 100 km s −1 ). They also find a relatively large and unusual NIR excess and argue that this is the result of an infrared light-echo originating from the presence of CSM. They estimate the mass of CSM to be more than 0.3 M , produced by a progenitor mass loss rate greater than 10 −4 M yr −1 . For the progenitor, Kotak et al. (2004) favor a single-degenerate system with a post-AGB companion star.
Wood-Vasey et al. (2004) obtained pre-maximum and late time photometry of SN 2002ic and find that a non-SN Ia component of the light curve becomes pronounced about 20 days post-explosion. They suggest the non-SN Ia component to be due to heating from a shock interaction between SN ejecta and CSM. Wood-Vasey et al. (2004) also suggest that the progenitor system consisted of a WD and an AGB star in the protoplanetary nebula phase. Wood-Vasey et al. (2004) and Sokoloski et al. (2006) proposed that a nova shell ejected from a recurrent nova progenitor system, creating the evacuated region around the explosion center of SN 2002ic. They suggest that the periodic shell ejections due to nova explosions on a WD sweep up the slow wind from the binary companion, creating density variations and instabilities that lead to structure in the circumstellar medium. This type of phenomenon may occur in SN Ia with recurrent nova progenitors, however Schaefer (2011) recently reported on an ongoing observational campaign of recurrent novae (RN) orbital period changes between eruptions. For at least two objects (CI Aquilae and U Scorpii), he finds that the RN lose mass, thus making RN unlikely progenitors for SN Ia.
Nearly one year after the explosion, Wang et al. (2004) found that the supernova had become fainter overall, but Hα emission had brightened and broadened compared to earlier observations. From their spectropolarimetry observations, Wang et al. (2004) find that hydrogen-rich matter is asymmetrically distributed. Likewise, Deng et al. (2004) also found evidence of a hydrogen-rich asymmetric circumstellar medium. From their observations of SN 2002ic, Wang et al. (2004) conclude that the event took place within a "dense, clumpy, disk-like" circumstellar medium. They suggest that the star responsible for SN 2002ic could either be a post-AGB star or WD companion (see also Hachisu et al. 1999;Han & Podsiadlowski 2006).

SN 2005gj
Similar to SN 2002ic, Aldering et al. (2006) argue that SN 2005gj is a SN Ia in a massive circumstellar envelope (see also Prieto et al. 2007), which is located in a low metallicity host galaxy with a significant amount of star formation. Their first spectrum of SN 2005gj shows a blue continuum level with broad and narrow Hα emission. Subsequent spectra reveal muted SN Ia features combined with broad and narrow Hγ, Hβ, Hα and He I λλ5876, 7065 in emission, where high resolution spectra reveal narrow P Cygni profiles. An inverted P Cygni profile for [O III] λ5007 was also detected, indicating top-lighting effects from CSM interaction . From their early photometry of SN 2005gj, Aldering et al. (2006) find that the interaction between the supernova ejecta and CSM was much weaker for SN 2002ic. Notably, both Aldering et al. (2006) and Prieto et al. (2007) agree that a SN 1991T-like spectrum can account for many of the observed profiles with an assumed increase in continuum radiation from interaction with the hydrogen-rich material. Aldering et al. (2006) also find that the light curve and measured velocity of the unshocked CSM imply mass loss as recent as 1998. This is in contrast to SN 2002ic, for which an inner cavity in the circumstellar matter was inferred (Wood-Vasey et al. 2004 all exhibit large CSM interactions and are from lowluminosity hosts.
Consistent with this interpretation for CSM interactions is the recent report by Fox & Filippenko (2013) that a NIR re-brightening, possibly due to emission from "warm" dust, took place at late times for both SN 2002ic and 2005gj. Notably, and in contrast to SN 2002ic, Fox & Filippenko (2013 find that the mid-IR luminosity of SN 2005gj increased to ∼ twice its early epoch brightness.

Summary and Concluding Remarks
Observations of a significant number of SN Ia during the last two decades have enabled us to document a larger expanse of their physical properties which is manifested through spectrophotometric diversity. While in general SN Ia have long been considered a homogeneous class, they do exhibit up to 3.5 mag variations in the peak luminosity, whereas "normal" SN Ia dispersions are ∼1 mag, and constitute several marginally distinct subtypes (Blondin et al. 2012;Scalzo et al. 2012;Silverman et al. 2013d;Foley et al. 2013;Dessart et al. 2013a). Consequently, the use of normal SN Ia for cosmological purposes depends on empirical calibration methods (e.g., Bailey et al. 2009), where one of the most physically relevant methods is the use of the widthluminosity relation (Phillips 1993;Phillips et al. 1999).
Understanding the physics and origin of the widthluminosity-relationship of SN Ia light curves is an important aspect in the modeling of SN Ia (Khokhlov et al. 1993;Lentz et al. 2000;Timmes et al. 2003;Nomoto et al. 2003;Kasen et al. 2009;Meng et al. 2011;Blondin et al. 2011). Brighter SN Ia often have broad light curves that decline slowly after peak brightness. Slightly less bright or dimmer SN Ia have narrower and relatively rapidly declining light curves. In addition, several SN Ia do not follow the width-luminosity-relationship (e.g., SN 2001ay, 2004dt, 2010jn, SCC, CL and SN 2002cx-like SN Ia), which reinforces the notion that a significant number of physically relevant factors influence the diversity of SN Ia overall (see Wang et al. 2012;Baron et al. 2012).
Despite the ever increasing number of caught-early supernovae, our perspective on their general properties and individual peculiarities undergoes a continual convergence toward a set of predictive standards with which models must be seen to comply. The most recent observational example is that of SN 2012fr (Maund et al. 2013;Childress et al. 2013c), a normal/lowvelocity-gradient SN Ia that has been added to the growing list of similar SN Ia that exhibit stark evidence for a distinctly separate region of "high-velocity" material (>16,000 km s −1 ). While the origin of high velocity features in the spectra of SN Ia is not well understood, it is concurrent with polarization signatures in most cases which implies some amount of ejecta density asymmetries (e.g., Kasen et al. 2003;Wang & Wheeler 2008;Smith et al. 2011;Maund et al. 2013). Furthermore, since understanding the temporal behavior of high velocity Si II/Ca II depends on knowing the same for the photospheric component, studies that focus on velocity gradients and potential velocity-plateaus of the photospheric component could make clearer the significance of the physical separation between these two regions of material (see Patat et al. 1996;Kasen et al. 2003;Tanaka et al. 2008;Foley et al. 2012c;Parrent et al. 2012;Scalzo et al. 2012;Childress et al. 2013c;Marion et al. 2013). However, it is at least certain that all viable models that encompass "normal" SN Ia conditions must account for the range of properties related to velocity evolution (see Blondin et al. 2012 and references therein), the occasionally observed however potentially under-detected signatures of C+O material at both low and high velocities (Thomas et al. 2007;Parrent et al. 2011;Thomas et al. 2011b;Folatelli et al. 2012;Piro & Nakar 2013;Mazzali et al. 2013), a high-velocity region of either clumps or an amorphous plumage of opaque Si-, Ca-based material (Gamezo et al. 2004;Leonard et al. 2005;Wang et al. 2007;Maund et al. 2010b;Piro 2011), and the supposed blue/red-shift of nebular lines emitted from the inner IPE-rich material (Maeda et al. 2010b;Mc-Clelland et al. 2013;Silverman et al. 2013a).
For at least normal SN Ia, there remain two viable explosion channels (with a few sub-and super-M Ch subchannels) regardless of the hierarchical dominance of each at various redshifts and/or ages of galactic constituents (c.f., Röpke et al. 2012;Hachisu et al. 2012;Seitenzahl et al. 2013;Pakmor et al. 2013;Moll et al. 2013;Claeys et al. 2014). Also, it may or may not be the case that some SN Ia are 2+ subtypes viewed upon from various lines of sight amidst variable CSM interaction (Maeda et al. 2010b;Foley et al. 2012a;Scalzo et al. 2012;Leloudas et al. 2013;Dan et al. 2013;Moll et al. 2013;Dessart et al. 2013a). However, with the current lack of complete observational coverage in wavelength, time, and mode (i.e. spectrophotometric and spectropolarimetric observations) for all SN Ia subtypes and "well-observed" events, there is a limit for how much constraint can be placed on many of the proposed explosion models and progenitor scenarios. That is, despite observational indications for and theoretical consistencies with the supposition of multiple progenitor channels, the observed diversity of SN Ia does not yet necessitate that each spectrophotometric subtype be from a distinctly separate explosive binary scenario than that of others within the SN Ia family of observed events; particularly so for normal SN Ia.
For the purposes of testing the multifaceted predictions of theoretical explosion models, time series spectroscopic observations of SN Ia serve to visualize the post-explosion material of an unknown progenitor system. For example during the summer of 2011 astronomers bore witness to SN 2011fe, the best observed normal type Ia supernova of the modern era. The prompt discovery and follow-up of this nearby event uniquely allowed for a more complete record of observed properties than all previous well-observed events. More specifically, the full range (in wavelength and time) of rapid spectroscopic changes was documented with continual day-to-day follow-up into the object's postmaximum light phases and well beyond. However, the observational side of visualizing other SN Ia remains inefficient without the logistical coordination of many telescope networks (e.g., LCOGT;Brown et al. 2013), telescopes large enough to make nearly all SN "nearby" in terms of improved signal-to-noise ratios (e.g., The Thirty Meter Telescope, The Giant Magellan Telescope), or a space-based facility dedicated to the study of such time sensitive UV−optical−NIR transients.
Existing SN Ia surveys are currently acting toward optimizing a steady flow of discoveries, while other programs have produced a significant number of publicly available spectra (Richardson et al. 2001;Matheson et al. 2008;Blondin et al. 2012;Yaron & Gal-Yam 2012;Silverman et al. 2012c;Folatelli et al. 2013). However, for the longterm future we believe it is imperative to begin a discussion of a larger (digital) network of international collaboration by way of (data-) cooperative competition like that done for both The Large Hadron Collider Experiment and Fermi Lab's Tevatron, with multiple competing experiments centered about mutual goals and mutual resources. Otherwise we feel the simultaneous collection of even very high quality temporal datasets by multiple groups will continue to create an inefficient pursuit of over-observing the most high profile event(s) of the year with a less than complete dataset.
Such observational pursuits require an increasingly focused effort toward observing bright and nearby events. For example, 206 supernovae were reported in 1999 and 67 were brighter than 18th magnitude while only three reached ∼ 13 magnitude 21 . By 2012 the number of found supernovae increased to 1045 while 78 were brighter than 16th magnitude and five brighter than 13th magnitude. This clearly indicates that supernovae caught early are more prevalent than ∼ 15 years ago and it is worthwhile for multiple groups to continually increase collaborative efforts for the brightest events. Essentially this could be accomplished without interfering with spectrum-limited high[er]-z surveys by considering a distance threshold ( 10−30 Mpc) as part of the public domain. Additionally, surveys that corroborate the immediate release of discoveries would further increase the number of well-observed events and 21 Quoted from the archives page of http://www. rochesterastronomy.org/snimages/. could be supplemented and sustained with staggered observations given that there are two celestial hemispheres, unpredictable weather patterns, and caughtearly opportunities nearly every week during active surveying.
In conclusion, to extract details of the spectroscopic behavior for all SN Ia subtypes, during all phases, larger samples of well-observed events are essential, beginning from as close to the onset of the explosion as possible (e.g., SN 1999ac, 2009ig, 2011fe, 2012cg, 2012fr), where SN Ia homogeneity diverges the most (see Zheng et al. 2013 for the most recent instance in SN 2013dy). Nearcontinuous temporal observations are most important for at least the first 1−2 months post-explosion and biweekly to monthly follow-up thereafter for ∼ 1 year. SN Ia spectra are far too complicated to do so otherwise. Even normal SN Ia deserve UV−optical−IR spectroscopic follow-up at a 1:1 to 2:1 ratio between days passed and spectrum taken, whenever possible, given that fine differences between normal SN Ia detail the variance in explosion mechanism parameters and initial conditions of their unobserved progenitor systems. It is through such observing campaigns that the true diversity to the underlying nature of SN Ia events will be better understood.

Acknowledgements
This work was supported in part by NSF grant AST-0707704, and US DOE Grant DE-FG02-07ER41517, and by SFB 676, GRK 1354 from the DFG. Support for Program number HST-GO-12298.05-A was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. We wish to acknowledge the use of the Kurucz & Bell (1995)  This review was made possible by collaborative discussions at the 2011 UC-UC-HIPACC International As-troComputing Summer School on Computational Explosive Astrophysics. We would like to thank Dan Kasen and Peter Nugent for organizing the program and providing the environment for a productive summit, and we hope that such programs for summer learning opportunities will continue in the future. We MP is thankful to Prof. David Branch, Prof. Eddie Baron, Prof. Ramanath Cowsik, Prof. Shoken Miyama, Prof. Masahiko Hayashi, Prof. Yoichi Takeda, Prof. Wako Aoki, Prof. Ajit Kembhavi, Prof. Kandaswamy Subramanian, and Prof. T. Padamanabhan for their kind support, encouragement, and hospitality. JTP would like to thank the University of Oklahoma Supernova Group, Rollin Thomas, and Alicia Soderberg for several years of support and many enlightening discussions on reading supernova spectra. JTP wishes to acknowledge helpful discussions with B. Dilday, R. A. Fesen, R. Foley, M. L. Graham, D. A. Howell, G. H. Marion, D. Milisavljevic, P. Milne, D. Sand, and S. Valenti, as well as S. Perlmutter for an intriguing conversation on "characteristic information of SN Ia" at the 221st American Astronomical Society Meeting in Long Beach, CA. JTP is also indebted to Natalie Buckley-Medrano for influential comments on the text and figures presented here.

SN Name
References SN 2006X Wang et al. 2008bSN 2006bt Hicken et al. 2009bFoley et al. 2010bSN 2006gz Hicken et al. 2007Scalzo et al. 2010SN 2007ax Kasliwal et al. 2008SN 2007if Scalzo et al. 2010SN 2007qd McClelland et al. 2010SN 2008J Taddia et al. 2012SN 2008ha Foley et al. 2009SN 2009dc Yamanaka et al. 2009bScalzo et al. 2010;Silverman et al. 2011;Taubenberger et al. 2011SN 2009ig Foley et al. 2012cSN 2009ku Narayan et al. 2011SN 2009nr Khan et al. 2011bTsvetkov et al. 2011PTF09dav Sullivan et al. 2011bSN 2010jn Hachinger et al. 2013SN 2011fe Richmond & Smith 2012Munari et al. 2013;Pereira et al. 2013SN 20011iv Foley et al. 2012bSN 2012cg Silverman & Filippenko 2012Munari et al. 2013SN 2012fr Childress et al. 2013c A Some recent SN Ia, Continued A.1 Peculiar SN 1997br in ESO 576-G40 Li et al. (1999) presented observations of the peculiar SN 1991T-like, SN 1997br (Bao Supernova Survey et al. 1997. Hatano et al. (2002) analyzed the spectra of SN 1997br and raised the question of whether or not Fe III and Ni III features in the early spectra are produced by 54 Fe and 58 Ni rather than by 56 Fe and 56 Ni. In addition, Hatano et al. (2002) discussed the issue of SN 1991T-like events as more powerful versions of normal SN Ia, rather than a physically distinct subgroup of events.
A.2 SN 1997cn in NGC 5490 Turatto et al. (1998) studied the faint SN 1997cn, which is located in an elliptical host galaxy Turatto et al. 1997). Like SN 1991bg, spectra of SN 1997cn show a deep Ti II trough between 4000 and 5000Å, strong Ca II IR3 absorption features, a large R("Si II"), and slow mean expansion velocities.
A.3 SN 1997ff and other "farthest known" SN Ia With a redshift of z = 1.7, SN 1997ff was the most distant SN Ia discovered at that time (Riess et al. 2001;Benítez et al. 2002). There have been ∼110 high-z (1 < z < 2) SN Ia discoveries since SN 1997ff (Riess et al. 2004(Riess et al. , 2007Suzuki et al. 2012), with the most recent and and one of the most distant SN Ia known being "SN UDS10Wil" at z = 1.914 (Jones et al. 2013). With these and future observations of "highest-z" SN Ia, constraints on DTD timescales (Strolger et al. 2010;Graur et al. 2013) and dark energy (Rubin et al. 2013) will certainly improve.
A.4 SN 1998aq in NGC 3982 Branch et al. (2003) used SYNOW to study 29 optical spectra of the normal SN 1998aq (Hurst et al. 1998), covering 9 days before to 241 days after maximum light (days −9 and +241, respectively). Notably, they find evidence for C II down to 11,000 km s −1 , ∼ 3000 km s −1 below the cutoff of carbon in the pure deflagration model, W7 (Nomoto et al. 1984).
A.5 SN 1999aa in NGC 2595 From day −11 to day +58, Garavini et al. (2004) obtained 25 optical spectra of SN 1999aa (Armstrong & Schwartz 1999). While SN 1999aa appears SN 1991Tlike, Garavini et al. (2004) note that the Ca II absorption feature strengths are between those of SN 1991T (SS) and the SN 1990N (CN), along with a phase transition to normal SN Ia characteristics that sets in earlier than SN 1991T. Subsequently, they suggest SN 1999aa to be a link between SN 1991T-likes and spectroscopically normal SN Ia. Evidence of carbon-rich material is also found in SN 1999aa; decisively as C II λ6580, tentatively as C III λ4649 (see Parrent et al. 2011). A schematic representation of their SYNOW fitting results is also presented, showing the inference of Co II, Ni II, and Ni III during the pre-maximum phases. These results deserve further study from more detailed models.
A.6 SN 1999ac in NGC 2841 Between day −15 and day +42, Garavini et al. (2005) obtained spectroscopic observations of the unusual SN 1999ac (Modjaz et al. 1999). The pre-maximum light spectra are similar to that of SN 1999aa-like, while appearing spectroscopically normal during later epochs. Garavini et al. (2005) find evidence of a fairly conspicuous, heavily blended C II λ6580 feature in the day −15 spectrum with approximate ejection velocities >16,000 km s −1 . By day −9, the C II absorption feature is weak or absent amidst blending with the neighboring 6100Å feature. This alone indicates that studies cannot fully constrain SN Ia models without spectra prior to day −10.
A.7 SN 1999aw in a low luminosity host galaxy Strolger et al. (2002) find SN 1999aw to be a luminous, slow-declining SN Ia, similar to 1999aa-like events. Strolger et al. (2002) derive a peak luminosity of 1.51 x 10 43 and a 56 Ni mass of 0.76 M .
A.8 SN 1999by in NGC 2841 Vinkó et al. (2001) presented and discussed the first three pre-maximum light spectra of SN 1999by (Arbour et al. 1999), where they find it to be a sub-luminous SN Ia similar to SN 1991bg (Filippenko et al. 1992b;Leibundgut et al. 1993;Turatto et al. 1996), SN 1992K (Hamuy et al. 1994), SN 1997cn (Turatto et al. 1998), and SN 1998de (Modjaz et al. 2001); in addition, the list of sub-luminous SN Ia include SN 1957A, 1960H, 1971I, 1980I, 1986G (see Branch et al. 1993Doull & Baron 2011 and references therein) and several other recently discovered under-luminous SN Ia (Howell 2001;McClelland et al. 2010;Hachinger et al. 2009). Premaximum spectra of SN 1999by show relatively strong features due to O, Mg, and Si, which are due to explosive carbon burning. In addition, blue wavelength regions reveal spectra dominated by Ti II and some other IPEs. Meanwhile, Höflich et al. (2002) studied the infrared spectra of SN 1999by, covering from day −4 to day +14. Post-maximum spectra show features which can be attributed to incomplete Si burning, while further support for incomplete burning comes from the detection of a pre-maximum C II absorption feature (Garnavich et al. 2004). Höflich et al. (2002) analyzed the spectra through the construction of an extended set of delayed detonation models covering the entire range of normal to sub-luminous SN Ia. They estimate the 56 Ni mass for SN 1999by to be on the order of 0.1 M . Garnavich et al. (2004) obtained UBVRIJHK light curves of SN 1999by. From the photometry of SN 1999by, the recent Cepheid distance to NGC 2841 (Macri et al. 2001), and minimal dust extinction along the line-ofsight, Garnavich et al. (2004) derive a peak absolute magnitude of M B = −17.15.
In order to assess the role of material asymmetries as being responsible for the observed peculiarity of subluminous SN Ia, Howell et al. (2001) obtained polarization spectra of SN 1999by near maximum light. They find relatively low levels of polarization (0.3%−0.8%), however significant enough to be consistent with a 20% departure from spherical symmetry (Maund et al. 2010b).
A.9 SN 1999ee in IC 5179 From day −10 to day +53, Stritzinger et al. (2002) obtained well-sampled UBVRIz light curves of SN 1999ee . They find the B -band light curve is broader than normal SN Ia, however sitting toward the over-luminous end of SN Ia peak brightnesses, with M B = −19.85 ± 0.28 and ∆m 15 = 0.94. Hamuy et al. (2002) obtained optical and infrared spectroscopy of SN 1999ee between day −9 and day +42. Before maximum light, the spectra of SN 1999ee are normal, with relatively strong Si II 6100Å absorption, however within the SS subtype (Branch et al. 2009). Hamuy et al. (2002) compared the infrared spectra of SN 1999ee to that of other SN Ia out to 60 days post-explosion, and find similar characteristics for SN 1999ee and 1994D (Meikle et al. 1996).
A.10 SN 2000E in NGC 6951 Valentini et al. (2003) obtained UBVRIJHK photometry and optical spectra of SN 2000E, which is located in a spiral galaxy. Optical spectra were obtained from 6 days before B -band maximum to 122 days after Bband maximum. The photometric observations span 230+ days, starting at day −16. The photometric light curves of SN 2000E are similar to other SS SN Ia, however SN 2000E is classified as a slowly declining, spectroscopically "normal" SN Ia similar to SN 1990N. Valentini et al. (2003) estimate the 56 Ni mass to be 0.9 M from the bolometric light curve. A.11 SN 2000cx in NGC 524 One of the brightest supernovae observed in the year 2000 was the peculiar SN 2000cx, located in an S0 galaxy (Yu et al. 2000;Li et al. 2001;Candia et al. 2003). It was classified as a SN Ia with a spectrum resembling that of the peculiar SN 1991T (Chornock et al. 2000). Sollerman et al. (2004) obtained late time BVRIJH light curves of SN 2000cx covering 360 to 480 days after maximum. During these epochs, they find relatively constant NIR magnitudes, indicating the increasing importance (with time) of the NIR contribution to the bolometric light curve. Branch et al. (2004b) decomposed the photosphericphase spectra of SN 2000cx with SYNOW. Apart from confirming HVFs of Ca II IR3 (which are consistent with primordial abundances; Thomas et al. 2004), Branch et al. (2004b) also find HVFs of Ti II. They attribute the odd behavior of SN 2000cx's B -band light curve to the time-dependent behavior of these highly line blanketing Ti II absorption signatures. Branch et al. (2004b) find an absorption feature near 4530Å in the spectra of SN 2000cx that can be tentatively associated with Hβ at high velocities, however this feature is more likely due to C III λ4649 or S II/Fe II instead (see Parrent et al. 2011 and references therein). Rudy et al. (2002) obtained 0.8−2.5 µm spectra of SN 2000cx at day −7 and day −8 before maximum light. From the λ10926 line of Mg II, they find that carbon-burning has taken place up to ∼25,000 km s −1 . Given the SS subtype nature of SN 2000cx, the early epoch IR spectra of Rudy et al. (2002) are valuable for comparison with other SN Ia IR datasets.
A.12 SN 2001V in NGC 3987 Vinkó et al. (2003) presented photometry of SN 2001V . They find that SN 2001V is overluminous, relative to the majority of SN Ia. Spectroscopic observations, spanning from day −14 to day +106, can be found in Matheson et al. (2008) and reveal it to be a SS SN Ia, consistent with its observed brightness.
A.13 Slowly declining SN 2001ay in IC 4423 Krisciunas et al. (2011) obtained optical and near infrared photometry, and optical and UV spectra of SN 2001ay (Swift et al. 2001). They find maximum light Si II and Mg II line velocities of ∼ 14,000 km s −1 , with Si III and S II near 9,000 km s −1 . SN 2001ay is one of the most slowly declining SN Ia. However, a ∆m 15 (B ) = 0.68 is odd given it is not over-luminous like SCC SN Ia slow decliners. In fact, the 56 Ni yield of 0.58 is comparable to that of many normal SN Ia. Baron et al. (2012) note this apparent WLR violation is related to a decrease in γ-ray trapping deeper within the ejecta due to an overall outward shift of 56 Ni, thus creating a fast rise in brightness followed by a slow decline caused by enhanced heating of the outer regions of material, which is a consequence of the larger expansion opacities.
A.14 SN 2001el in NGC 1448 Krisciunas et al. (2003) obtained well-sampled UBVRI-JHK light curves of the nearby (about 18 Mpc) and normal SN 2001el (Monard et al. 2001), from day −11 to day +142. Because Krisciunas et al. (2003) obtained UBVRI and JHK light curves, they were able to measure a true optical−NIR reddening value (A V = 0.57 mag along the line-of-sight) for the first time. Mattila et al. (2005) obtained early time high resolution and low resolution optical spectra of SN 2001el. They estimate the mass loss rate (assuming 10−50 km s −1 wind velocities) from the progenitor system of SN 2001el to be no greater than 9 x 10 −6 M yr −1 and 5 x 10 −5 M yr −1 , respectively. The low resolution spectrum was obtained 400 days after maximum light with no apparent signatures of hydrogen Balmer lines. High velocity Ca II was detected out to 34,000 km s −1 , while the 6100Å absorption feature is suspected of harboring high velocity Si II (see also Kasen et al. 2003).

A.15 SN 2002bo in NGC 3190
Between day −13 and day +102, Benetti et al. (2004) collected optical and NIR spectra and photometry of the BL SN 2002bo (Cacella et al. 2002;Krisciunas et al. 2004). Estimates on host galaxy extinction from Na D equivalent width measurements are consistent with the inferred color excess determined by comparison to the Lira relation (Lira 1995;Riess et al. 1996;Phillips et al. 1999). From the time-evolution of the 6100Å absorption feature, Benetti et al. (2004) find that SN 2002bo is an intermediary between the BL SN 1984A and the CN SN 1994D. Benetti et al. (2004) also discuss SN 2002bo argue that some of the IME high velocity material may be primordial, while most is produced during the explosion and possibly by prolonged burning in a delayed detonation. This interpretation is also consistent with a lack of any clear signatures of unburned carbon. Stehle et al. (2005) studied the abundance stratification by fitting a series of spectra with a Monte Carlo code and found that the elements synthesized in different stages of burning are not completely mixed within the ejecta. In the case of SN 2002bo, they derived the total mass of 56 Ni to be 0.52 M . Similar to SN 2001ay's fast rise to maximum light , Stehle et al. (2005) attribute SN 2002bo's fast rise to outward mixing of 56 Ni.

A.16 SN 2002cv in NGC 3190
The NIR photometry of SN 2002cv reveal an obscured SN Ia (Di Paola et al. 2002); more than 8 magnitudes of visual extinction. Both optical and NIR spectroscopy indicate SN 2002cv is most similar to SN 1991T (Meikle et al. 2002;Filippenko et al. 2002). It should also be noted that the SS SN 2002cv and the BL SN 2002bo share the same host galaxy. Elias-Rosa et al. (2008) obtained and analyzed VRIJHK photometry, in addition to a sampling of optical and NIR spectroscopy near and after maximum light, and find a best fit value for the ratio between inferred extinction and reddening, R V = 1.59 ± 0.07 whereas 3.1 is often assumed for normal SN Ia (however see Tripp 1998;Astier et al. 2006;Krisciunas et al. 2006;Folatelli et al. 2010 (Hutchings & Li 2002). The dataset presented make it one of the most well-observed SN 1984A-like SN Ia and is a valuable tool for the discussion of SN Ia diversity.

A.18 SN 2002er in UGC 10743
From day −11 to day +215, Kotak et al. (2005) carried out spectroscopic follow-up for the reddened, CN SN 2002er (Wood-Vasey et al. 2002c. By contrast with the photometric behavior seen for SN 1992A, 1994D, and 1996X, SN 2002er stands out for its slightly delayed second peak in the I -band and similarly for V and R-bands. Pignata et al. (2004); Kotak et al. (2005) estimated the mass of 56 Ni to be on the order of 0.6 to 0.7 M , where the uncertainty in the exact distance to SN 2002er was the primary limitation.

A.19 SN 2003du in UGC 09391
For 480 days, and starting from day −13, Stanishev et al. (2007) monitored the CN SN 2003du. From modeling of the bolometric light curve, Stanishev et al. (2007) estimate the mass of 56 Ni to between 0.6 and 0.8 M . Like other normal SN Ia, the early spectra of SN 2003du contain HVFs of Ca II and a 6100Å feature that departs from being only due to photospheric Si II, suggesting either a distinctly separate region of HV Si II or the radial extension of opacities from below. Tanaka et al. (2011) studied the chemical composition distribution in the ejecta of SN 2003du by modeling a one year extended time series of optical spectra. Tanaka et al. (2011) do not find SN 2003du to be as fully mixed as a some 3D deflagration models. Specifically, from their modeling Tanaka et al. (2011) that the a core of stable IPEs supersedes 56 Ni out to ∼ 3000 km s −1 ( 0.2 in mass coordinate). Atop this 0.65 M of 56 Ni are layers of IMEs, while the outermost layers consist of oxygen, some silicon, and no more than 0.016 M of carbon above 10,500 km s −1 .
A.20 SN 2003gs in NGC 936 Krisciunas et al. (2009) obtained near-maximum to late time optical and NIR observations of SN 2003gs, offering a chance to study the post-maximum light bolometric behavior of a fast declining SN Ia that was subluminous at optical wavelengths, but of standard luminosity in NIR bands at maximum light. Krisciunas et al. (2009) find ∆m 15 (B ) = 1.83 ± 0.02 and discuss comparisons to other fast decliners; namely, SN 2003hv (Leloudas et al. 2009), SN 2004gs (Folatelli et al. 2010Contreras et al. 2010), SN 2005bl (Taubenberger et al. 2008Folatelli et al. 2010;Wood-Vasey et al. 2008), and SN 2005ke, 2006gt, and 2006mr (Folatelli et al. 2010Contreras et al. 2010). In particular, Krisciunas et al. (2009) note that, in contrast to normal and over-luminous SN Ia, the delay in the time of J -band maximum from that of the B -band, for fast decliners, is inversely proportional to the peak NIR magnitude. Furthermore, they discussed the possibility for two subsets of FAINT−CL fast decliners; those that do and do not show a J -band peak before the B -band (see also Kattner et al. 2012);respectively, SN 1986G, 2003gs, 2003hv, and 2006gt, and SN 1991bg, 1999by, 2005bl, 2005ke, and 2006mr. Krisciunas et al. (2009 conclude that differences in NIR opacity within the outer layers are responsible for dissimilar γ-ray trapping, and therefore longer J -band than B -band diffusion times for FAINT−CL SN Ia that are fainter in the NIR. However, the origin (differences of explosion mechanism and/or progenitor systems) for this apparent 'bimodal' difference in NIR opacity is not clear.
For SN 2003gs, Krisciunas et al. (2009 used UBVRI-JHK photometry and Arnett's rule (Arnett 1982, but also see Stritzinger & Leibundgut 2005) to estimate 0.25 M of 56 Ni was produced during the explosion. As for the optical spectra, SN 2003gs is similar to SN 2004eo (A.23) in that it is found to have absorption signatures that are consistent with its photometric characteristics; a FAINT−CL SN Ia with a larger than normal R(Si II) and the presence of Ti II features near 4000−4500Å.
A.21 SN 2003hv in NGC 1201 Leloudas et al. (2009) studied SN 2003hv out to very late phases (day +786). Notably, this seemingly spectroscopically normal SN Ia has ∆m 15 (B ) = 1.61, while the late time light curves show a deficit in flux that follow the decay of radioactive 56 Co, assuming full and instantaneous positron trapping. Leloudas et al. (2009) consider this as possibly due to a redistribution of flux for SN 2003hv (a.k.a. an infrared catastrophe, see Axelrod 1980) from a dense clumping of inner material, and would also explain the flat-topped nebular emission lines (Motohara et al. 2006;Gerardy et al. 2007). Mazzali et al. (2011) also studied the nebular spectrum of 2003hv and consider it as a non-standard event. They note that its late time flux deficit, compared to normal SN Ia, could be due to SN 2003hv having a lower mean density structure, possibly consistent with a sub-Chandrasekhar mass origin. Motohara et al. (2006) presented NIR Subaru Telescope spectra of SN 2003du, 2003hv, and 2005W during their late phase evolution (+200 days post-maximum light). For both SN 2003du and 2003hv, they find a flat-topped [Fe II] λ16440 emission feature that is blue-shifted by ∼ 2000 km s −1 from the SN rest frame (FWHM ∼ 4000 km s −1 ). Motohara et al. (2006) further argue that the [Fe II] emission would be rounded on top if the neutron-rich Fe-peak isotopes produced in the explosion were thoroughly mixed with the surrounding distribution of 56 Ni; for SN 2003du and 2003hv they suggest that this is not the case. In fact, they find that SN 1991T and 2005W (see their Fig. 1), at least, may represent instances where the inner most regions have been thoroughly mixed.
Similarly, Gerardy et al. (2007) looked to address the nature of the thermonuclear burning front by utilizing late time (+135 days) mid-IR (5.2−15.2 µm) Spitzer Space Telescope spectra of SN 2003hv and 2005df. In particular, Gerardy et al. (2007) find direct evidence in SN 2005df for a small inner zone of nickel that is surrounded by 56 Co and an asymmetric shell-like structure of Ar. While it is not clear why a supposed initial deflagration phase of a DDT explosion mechanism produces little to no mixing for these two SN Ia, the observations of Gerardy et al. (2007) give strong support for a stratified abundance tomography like those seen in DDT-like models; the various species of material are restricted to radially confined zones, which is inconsistent with the large-scale mixing that is expected to occur in 3D deflagration models. This is also in agreement with X-ray observations of the Tycho supernova remnant (Badenes et al. 2006) in addition to optical and UV line resonance absorption imaging of SNR 1885 in M31 (Fesen et al. 2007).
A.22 SN 2004dt in NGC 0799 Wang et al. (2006) and Altavilla et al. (2007) studied the early spectral evolution of SN 2004dt from more than a week before optical maximum, when line profiles show matter moving at velocities as high as 25,000 km s −1 . The variation of the polarization across some Si II lines approaches 2%, making SN 2004dt one of the most highly polarized SN Ia observed and an outlier in the polarization-nebular velocity plane (Maund et al. 2010b). In contrast with the polarization associated with Si II, Wang et al. (2006) find that the strong 7400 A O I−Mg II absorption complex shows little or no polarization signature. Wang et al. (2006) conclude this is due to a spherical geometry of oxygen-rich material encompassing a lumpy distribution of IMEs.

A.25 Under-luminous SN 2005bl in NGC 4070
Both Taubenberger et al. (2008) and Hachinger et al. (2009) studied the sub-luminous SN 2005bl with observations made between day −6 and day +66, and carried out spectral analysis ("abundance tomography") of SN 2005bl (Morrell et al. 2005). They find it to be one of incomplete burning similar to SN 1991bg and 1999by. Compared to SN 1991bg, a noteworthy difference of SN 1999by is the likely presence of carbon C II in pre-maximum spectra (Taubenberger et al. 2008), whereas C I λ10691 is also clearly detected in NIR spectra . However, this is likely a biased comparison to SN 1991bg given that the earliest spectrum obtained was on day −1 (potentially too late to detect unburned material via C II λ6580). To our knowledge no conspicuous C I λ10691 absorption features have been documented for other SN Ia. For example, C I λ10691 is present in the pre-maximum spectra of SN 2011fe but it is not a conspicuous signature. Similarly, pre-maximum spectra of SN 2005bl show less conspicuous detections of C I and C II but still indicate low burning efficiency with a significant amount of leftover unburned material (Taubenberger et al. 2008). Hachinger et al. (2009) suggest that a detonation at low pre-expanded densities is responsible for the abundance stratification of IMEs seen in the spectra of SN 2005bl. This would also explain the remaining carbon-rich material seen for some CL SN Ia when caught early enough.
A.26 SN 2005cf in MCG-01-39-003 Wang et al. (2009b) studied UV−optical−NIR observations of the normal SN 2005cf (see also Pastorello et al. 2007b). During the early evolution of the spectrum, HVFs of Ca II and Si II are found to be present above 18,000 km s −1 (confirming observations of Garavini et al. 2007). Gall et al. (2012) studied the NIR spectra of SN 2005cf at epochs from day −10 to day +42, which show clear signatures of Co II during postmaximum phases. In addition, they attribute fluorescence emission in making the underlying shape of the SED.
A.27 SN 2005cg in a low-luminosity, star forming host Quimby et al. (2006b) presented and discussed the spectroscopic evolution and light curve of the SS SN 2005cg, which was discovered by ROTSE-IIIc. Pre-maximum spectra reveal HVFs of Ca II and Si II out to ∼ 24,000 km s −1 and Quimby et al. (2006b) find good consistency between observed and modeled Si II profiles. They interpret the steep rise in the blue wing of the Si II to be an indication of circumstellar interaction given that abundance estimates for HVFs suggest modest amounts of swept up material (∼ 10 −4 − 10 −3 M ; see Quimby et al. 2006b;. A.28 SN 2005hj Quimby et al. (2007 obtained optical spectra of the SS SN 2005hj during pre-maximum and post-maximum light phases. From a ROTSE-IIIb unfiltered light curve, SN 2005hj reached an over-luminous peak absolute magnitude of −19.6 (assuming z = 0.0574). Interestingly, the sharp and shallow 6100Å feature remains fairly stagnant at ∼ 10,600 km s −1 near and after maximum light, with a sudden decrease at later epochs. Similar to Quimby et al. (2006b), Quimby et al. (2007) find this also consistent with the interpretation that CSM is influencing spectral profiles of SN 1999aa-like SN Ia (see also Scalzo et al. 2010Scalzo et al. , 2012. A.29 SN 2006D in MCG-01-33-034 Thomas et al. (2007) obtained the spectra of the spectroscopically normal SN 2006D from day −7 to day +13. The spectra show one of the clearest signatures of carbon-rich material at photospheric velocities observed in a normal SN Ia (below 10,000 km s −1 ). The 6300Å carbon feature becomes weaker with time and disappears as the photosphere recedes and the SN reaches maximum brightness. These observations−like all SN Ia diversity studies−underscore the importance of obtaining spectra of SN Ia during all phases. If [O I] and [C I] lines are present in the spectra during post-maximum light phases at velocities below 10,000 km s −1 , this would indicate the presence of unburned matter. These particular lines have not been detected, however the absence of said signatures does not imply a complete lack of C+O material at low velocities (Baron et al. 2003;Kozma et al. 2005).
A.30 SN 2006X in M100 Wang et al. (2008b) presented UBVRI and JK light curves and optical spectroscopy of the reddened BL SN 2006X (Stockdale et al. 2006;Quimby et al. 2006a) and find high mean expansion velocities during pre-maximum light phases ( 20,000 km s −1 ). Wang et al. (2008b) suggest the observed properties of SN 2006X may be due to interaction with CSM. Yamanaka et al. (2009b) also presented and discussed the early spectral evolution. They note that the R(Si II) ratio is unusually low for such a high-velocity gradient SN Ia. However, rather than this being an indication of low effective temperature, they suggest that the low R(Si II) value is due to line-blending, likely from a higher velocity component of Si II. Both Wang et al. (2008b) and Yamanaka et al. (2009b) find the observed properties of SN 2006X to be consistent with characteristics of delayed detonation models.
Equipped with high resolution spectra of narrow Na D signatures spanning ∼100 days post-maximum light, Patat et al. (2007) infer the presence of intervening CSM and argue a mass loss history associated with SN 2006X in the decades prior to explosion. In fact, at least half of all SN Ia with narrow rest frame, blueshifted Na D absorption profiles are associated with high ejecta velocities (Sternberg et al. 2011;Foley et al. 2012a). This indicates that CSM outflows are either present in some explosion scenarios or associated with all progenitors scenarios at some point during the lead up to the explosion. Patat et al. (2009) later discussed the VLT spectropolarimetry of SN 2006X. In particular, they find that the presence of the high-velocity Ca II is coincident with a relatively high polarization signature (∼1.4%) at day −10, that diminishes by only ∼15% near maximum light, and is still present 41 days later. Patat et al. (2009) note that this day +40 detection is not seen for SN 2001el (Wang et al. 2003) or SN 2004du (Leonard et al. 2005). As for the high-velocity Si II, its polarization signature is seen to peak (∼1.1%) at day −6, drop by ∼30% near maximum, and is undetected well into the post-maximum phase. While the findings of spectropolarimetry studies of SN Ia are thought to be associated with, for example, "deflagration phase plumes" with time-dependent photospheric covering fractions, Patat et al. (2009) are unable to conclude why SN 2006X exhibits a sizable post-maximum light re-polarization signature by day +39.
A.31 SN 2006bt in CGCG 108-013 Foley et al. (2010b) obtained optical light curves and spectra of transitional CN/CL SN 2006bt (Lee & Li 2006). The B -band decline rate, ∆m 15 (B) = 1.09, is within the range that is observed for normal SN Ia, however SN 2006bt shows a larger than normal R(Si II), slightly lower mean expansion velocities, and a lack of a double peak in the I -band; CL SN 1991bg-like properties. A tentative C II 6300 feature is identified, however with a minimum at ∼ 6450Å. Foley et al. (2010b) suggest this inferred lower projected Doppler velocity could be accounted for by a clump of carbon offset from the line of sight at photospheric velocities. Because of an association within a halo population of its passive host galaxy, Foley et al. (2010b) conclude that the progenitor was also likely to be from an old population of stars.
A.32 Over-luminous SN 2006gz in IC 1277 Hicken et al. (2007) studied SN 2006gz (Prieto et al. 2006a) and estimated a peak intrinsic V -band brightness of −19.74 and ∆m 15 (B ) = 0.69, implying M( 56 Ni) ∼ 1.0−1.2 M (assuming R V = 2.1−3.1; see also Maeda et al. 2009). The spectroscopic signatures during early phases are relatively narrow on account of slightly lower mean expansion velocities. At two weeks before maximum light, Hicken et al. (2007) attributed a relatively strong 6300Å feature to C II λ6580 that diminishes in strength by day −10 (Prieto et al. 2006b). Compared to a 5Å equivalent width 6300Å absorption feature observed in the CN SN 1990N (Leibundgut et al. 1991;Jeffery et al. 1992), the absorption feature has an observed equivalent width of 25Å in the early spectra of SN 2006gz (Hicken et al. 2007). So far, spectroscopic modeling that incorporate signatures of C II λ6580 predict carbon mass fractions, X(C), that span an order of magnitude and are broadly consistent with both singleand double-degenerate scenarios. Maeda et al. (2009) Kasliwal et al. (2008) find that it shares similarities with a sub-luminous SN 2005ke Hughes et al. 2007;Patat et al. 2012) and also shows clear excess UV emission ∼ 20 days post-maximum light. Based on the small amount of synthesized 56 Ni that is inferred (0.05 − 0.09 M ), along with SN Ia-like expansion velocities near maximum light (∼9000 km s −1 ), Kasliwal et al. (2008) conclude that SN 2007ax is not compatible with a number of theoretical models that have been proposed to explain FAINT−CL SN Ia.
A.34 Over-luminous SN 2007if Scalzo et al. (2010) find that SN 2007if qualifies as a SCC SN Ia, i.e. it is over-luminous (M V = −20.4), has a slow-rise to peak brightness (t rise = 24 days), the early spectra contain signatures of stronger than normal C II, and SN 2007if resides in a low-luminosity host (M g = −14.10). Despite having a red B − V color (+0.16) at B -band maximum, signs of host reddening via Na D lines appear negligible. Utilizing Keck observations of the young metal-poor host galaxy, Childress et al. (2011) concluded that SN 2007if is likely to have originated from a young, metal-poor progenitor. From the Hα line of the host galaxy, Yuan et al. (2010) derived a redshift of 0.0736.
Based on the bolometric light curve and the sluggish Si II velocity evolution, Scalzo et al. (2010) conclude that SN 2007if was the death of a super-Chandrasekhar mass progenitor. They estimate the total mass of the system to be 2.4 M , with 1.6 M of 56 Ni, and 0.3 to 0.5 M in the form of a C+O envelope. Given the possibility that other over-luminous events could potentially stem from similar super-Chandrasekhar mass origins, Scalzo et al. (2012) searched the SNFactory sample (based on a criterion of SN 1991T/2007if-like selections) and found four additional super-Chandrasekhar mass candidates.
A.35 SN 2007on in NGC 1404 SN 2007on was found associated with the elliptical galaxy, NGC 1404 (Pollas & Klotz 2007).  reported the discovery of the progenitor of SN 2007on based on a detected X-ray source in presupernova archival X-ray images, located 0.9" ± 1.3" (later 1.15" ± 0.27"; Roelofs et al. 2008) from the position of SN 2007on within its host galaxy. However, Roelofs et al. (2008) later reevaluated the detection of the progenitor of SN 2007on and concluded that given the offset discrepancy between the X-ray source and the SN location, the probability for a connection is of order 1 percent. However, should SN Ia progenitors reveal themselves to be producers of pre-explosive X-ray sources,  suggest this would be consistent with a merger model with an accretion disc, formed from the disrupted companion star rather than an explosion immediately upon or soon after the merger of the two stars.  A.40 PTF10ops, another peculiar cross-type SN Ia Maguire et al. (2011) presented optical photometric and spectroscopic observations of a somewhat peculiar and sub-luminous SN Ia, PTF10ops (−17.77 mag). Spectroscopically, this object has been noted as belonging to the CL class of SN Ia on account of the presence of conspicuous Ti II absorption features blue ward of 5000Å, in addition to a larger than normal R(Si II) ratio (partially indicative of cooler effective temperatures). Photometrically, PTF10ops overlaps "normal" SN Ia properties in ∆m 15 (B ) (1.12 ± 0.06 mag) and its rise-time to maximum light (19.6 days). Maguire et al. (2011) estimate ∼0.17 M of 56 Ni was produced during the explosion, which is well below what is expected for LVG−CN SN Ia. Maguire et al. (2011) also note that either PTF10ops remains without a visible host galaxy, or it resides within the outskirts of a massive spiral galaxy located at least 148 kpc away, which would be consistent with a possible influence of low metallicities or an old progenitor population. Maguire et al. (2011) suggest the progenitor could have been the merger of two compact objects (Pakmor et al. 2010), however time series spectrum synthesis is needed to confirm.

A.41 SN 2010jn in NGC 2929
The BL SN 2010jn was discovered by the Palomar Transient Factory (PTF10ygu) 15 days before it reached maximum light. Hachinger et al. (2013) performed spectroscopic analysis of the photospheric phase observations and find that the outer layers of SN 2010jn are rich in iron-group elements. At such high velocities (>16,000 km s −1 ), iron-group elements have been tentatively identified in the spectra of SN Ia before (Hatano et al. 1999a) and may also be a ubiquitous property of SN Ia. However, more early epoch, time series observations are needed in order to test and confirm such claims. For SN 2010jn at least, Hachinger et al. (2013) favor a Chandrasekhar-mass delayed detonation, where the presence of iron-group elements within the outermost layers may be a consequence of outward mixing via hydrodynamical instabilities prior to or during the explosion (see Piro 2011Piro , 2012. A.42 SN 2011iv in NGC 1404 Foley et al. (2012b) presented the first maximum-light UV through NIR spectrum of a SN Ia (SN 2011iv; Drescher et al. 2011). Despite having a normal looking spectrum, SN 2011iv declined in brightness fairly quickly (∆m 15 (B ) = 1.69). Since the UV region of a SN Ia spectrum is extremely sensitive to the composition of the outer layers, they offer the potential for strong constraints as soon as observational UV spectroscopic diversity is better understood.
A.43 SN 2012cg in NGC 4424 Silverman et al. (2012d) presented early epoch observations of the nearby spectroscopically normal SN 2012cg (Kandrashoff et al. 2012;Cenko et al. 2012;Marion et al. 2012), discovered immediately after the event (∼1.5 days after). Compared to the width of other normal SN Ia B -band light curves, Silverman et al. (2012d) find that SN 2012cg's light curve relatively narrow for its peak absolute brightness, with t rise = 17.3 days (coincident photometry was also presented by Munari et al. 2013). Mean expansion velocities within 2.5 days of the event were found to be more than 14,000 km s −1 , while the earliest observations show high-velocity components of both Si II and Ca II. The C II λλ6580, 7234 absorption features were also detected very early. Johansson et al. (2013) obtained upper limits on dust emission via far infrared Herschel Space Observatory flux measurements in the vicinity of the recent and nearby SN 2011by, 2011fe, and 2012cg. From non-detections during post-maximum epochs at 70 µm and 160 µm band-passes and archival image measurements, Johansson et al. (2013) exclude dust masses 7 x 10 −3 M for SN 2011fe, and 10 −1 M for SN 2011by and 2012cg for ∼ 500 K dust temperatures, ∼ 10 17 cm dust shell radii, and peak SN bolometric luminosities of ∼ 10 9 L .