arXiv:1008.2754v2 [astro-ph.HE] 25 Jan 2011 h yeI N20b.Te hwdisrs ie( time rise its showed 2010a; studied They (2010) al. 2002bj. al. et SN (Perets et Ib Poznanski literature type 2010). the the al. et in Poznanski discussed been a)adfloigdcy( decay following and day) ..Bntie l 05 e i.1.Nvrhls,it Nevertheless, 1). M Fig. of see luminosity 2005; peak al. a et reached Benetti e.g. n ats N19b-ieeet ( events 1991bg-like SN fastest and n vnfse hntefit(yial M (typically faint the than faster even and atrta yeI N swl srglrtp aSNe Ia type regular as well as SNe Ib ( type than faster h o jce as on o ohtpso N (SN SNe of types both for found solar mass, and few events, envi- ejected 2005E-like old a low SN of The the least class the evolution. at for gener- slower found of explo- much ronments are show core-collapse and (stripped) he- the spectra, masses, massive of from their of evidence originate in sion showing to hydrogen) thought i.e. not ally SNe, (but Ib lium Type markable. rich. are and known) a o an N h oa nerdms nisejecta its in mass inferred most at total be The to SN. found was faint a not was vlto ( evolution N ie rmS 02ja hyso an eklumi- peak faint show they M as (typically 2002bj nosities these SN However, from environments. differ SNe old with associated likely N20E,as nerdt have prototypical to the inferred (with also events Ib 2005E), type SN of class a described ∆ at abr t,Psdn,C 10,USA 91101, CA Pasadena, St., Barbara Santa Israel 69978, Aviv 7610 Rehovot Science, of Institute Weizmann Israel The Physics, USA of 02138, MA Cambridge, St., Garden 60 [email protected] rpittpstuigL using 2018 typeset 22, Preprint October version Draft eety aiu euirfs-vligSehave SNe fast-evolving peculiar various Recently, uhpoete o uenveo yeI r re- are Ib type of supernovae for properties Such m 4 3 2 1 15 bevtre fteCrei nttto fWsigo,8 Washington, of Institution Carnegie the Tel- of Observatories University, Tel-Aviv Astronomy, and Physics of School Faculty Astrophysics, and Physics Particle of Astrophysics, Department for Center Harvard-Smithsonian fellow, CfA NEEGN LS FBIH,FS-VLIGSPROA WITH SUPERNOVAE FAST-EVOLVING BRIGHT, OF CLASS EMERGING AN ( B hwdahg eklmnst,wt bouemgiueM magnitude absolute with luminosity, peak high a showed hwdeiec o eylwms feet ( ejecta of mass of low time very rise for a evidence had showed It curve. light nnab aaissget httert fteepcla S peculiar discove these the of helium-s Finally, rate a the headings: 2002bj. from that Subject SN expected suggests of be origin no nearby might the in with as env SNe. be 1939B, old these to elliptical SN for The suggested an t of progenitors 1939B. (in of dwarf spectrum SN 1939B white the analysis SN of old our and spectrum support M31) from c strongly the of light inferred time bulge fast ejecta first the similarly the of showing for mass 1939B, present low SN the and discuss 1885A SN events, aa .Perets B. Hagai ) eetaayi fsproa(N 02jrvae hti w it that revealed 2002bj (SN) of analysis recent A ∆ < m 1 15 . 7 ( ..Bntie l 05 e i.1), Fig. see 2005; al. et Benetti e.g. ; B 1. ) B ∼ A INTRODUCTION T E − ∼ 1 ∼ tl mltajv 11/10/09 v. emulateapj style X . 9 0 lhuhters iei un- is time rise the although , 1 ht dwarfs white . alsBadenes Carles , 15 uenve eea uenve niiul(85,1939 (1885A, individual supernovae: – general supernovae: ∆ 15 m a) hyhv slower a have they mag), M 15 B ⊙ ( < B eese l (2010a) al. et Perets . − ∼ 1 1 ) . 7 M ∼ < 18 ⊙ 3 ∆ . < B 5 . feet,and ejecta, of 2 m − ∼ a much was ) 7 a,ie it i.e. mag, rf eso coe 2 2018 22, October version Draft 15 2,3 asfloe yadcieo .5mgprdyi -ad and B-band, in day per mag 0.25 of decline a by followed days ( 17 arArcavi Iair , B ) mag) < < 13 < 2 Abstract 0, 7 ; 0 . 15 n 99)smlrt 02ji h bevdB/V-band. 1885A observed the (SNe in curves 2002bj to light similar fast-evolving 1939B) and but for bright case with the not is this However, spectra; 2002bj. rich SN 2010a,b). calcium al. their et and Perets environments support their stellar (including properties old scenario independent detonation helium-shell several WD and the SNe) (8 2010). ple al. et Waldman 2010; models al. these from et how- expected SNe, (Shen that these with of comparable under is lu- time ever, 2005E, rise lower SN The while to conditions. slowly similar 2010), appropriate possibly al. more are et decay models (Waldman minosity luminosities, 2002bj SN similar and light of than at expected some (2010) models, peak the al. model- SN which et that helium-detonation Waldman curve their show from light (2010), (2010) curves detailed al. Kasen et & However, Woosley Shen by 2002bj. ing SN in vnfse hnosre o N20E(idtne al. ’.Ia’ et (Bildsten called 2005E SN evolve (so e.g. for curves fast observed light 2007), mass of than their ejected faster suggested production even initially 2010). low They al. the et with Kasen SNe). Shen studied & SNe Woosley 2007; (2007) al. evolving 2010; et al. Bildsten et al. Bildsten 1987; et Waldman al. et 2010; Iben (e.g. bn&Ttkv19;Lve&Ant 95.Vari- deflagration) ejecting helium (including SNe, 1995). models 1987; faint Arnett al. sub- these & predicted et Livne Iben of helium-accreting 1986; 1991; ants al. Tutukov et of & (Woosley Iben a ago predicted explosion decades theoretically few was the WDs mass Chandrasekhar from al. et Perets SNe 2010; 2010). al. al. et et Poznanski Waldman (WDs; 2010a; 2010b; dwarfs al. white mecha- et on suggested Perets detonations One helium-shell is them. differ- nism for suggest model(s) therefore explosion 2002bj), ent SN and SNe 2005E-like M 2 ohaD Simon D. Joshua , ntefloigw rsn h td ftoohrSNe other two of study the present we following the In sam- larger a 2005E, SN to similar SNe of case the In h xsec ftemncerhelium-detonation thermonuclear of existence The ⊙ .Hr edsustoadtoa historical additional two discuss we Here ). R ei tlat12%o l SNe. all of % 1-2 least at is Ne − ∼ yo l h bevdfs-vligSNe fast-evolving observed the all of ry efidn la vdnefrhlu in helium for evidence clear no find We sa paetyuiu yeI N It SN. Ib type unique apparently an as eS 85 enn nM1 and M31, in remnant 1885A SN he ∆ 18 re n o jce ass We masses. ejected low and urves eldtnto nawiedwarf, white a on detonation hell m rcso trfrainactivity), formation of traces . rnet fbt N18A(in 1885A SN both of ironments 5 15 u neteeyfast-evolving extremely an but , ( B 4 ) n vsa Gal-yam Avishay and > 2 osbya ata observed as fast as possibly , ,20b)–(stars:) – 2002bj) B, O-ASEJECTA LOW-MASS < 1 M ⊙ fmaterial of 2 2

straint on the SN ejecta. These events also provide us −20 with further data on the frequency of such SNe. Finally, the lack of clear evidence for helium in the spectrum of −19 SN 1939B could give an additional, although puzzling, SN 1939B clue to the origin of these SNe and the processes involved SN 2002bj in their explosions. −18 Type Ia SNe

2. OBSERVATIONAL DATA −17 SN 1885A 2.1. SN 2002bj This recently studied SN was shown to be an unusually −16 fast-evolving type Ib SN (Poznanski et al. 2010). Here we reproduce some of its observations to serve as context −15 for the discussion of the two additional fast-evolving SNe

Absolute B−Band Peak Luminosity [mag] SN 2005E described below. −14 0.5 1 1.5 2 2.5 3 3.5 ∆ m (B) [mag] 2.2. SN 1885A 15 SN 1885A (S Andromeda), discovered in the bulge of Fig. 1.— ∆m15(B) for various SNe. Data shown for type Ia SNe (× signs; taken from Benetti et al. 2005, table 1), in- M31, was the first SN reported in the modern era of as- cluding faint fast SNe similar to SN 1991bg (e.g. SNe with tronomy. Its close proximity (the second closest after SN ∆m15(B) > 1.7), as well as estimated ∆m15(B) for SN 2005E, in- 1987A) enabled many observations of its light curve to ◦ terpolated from its light curve ( sign; Perets et al. 2010a). These be made following its discovery in August 1885 (Fig. 2; can be compared with the ∆m15(B) measured for SN 2002bj (⋆ sign; Poznanski et al. 2010) and SN 1939B ( sign), and the esti- based on de Vaucouleurs & Corwin 1985, data extracted mated ∆m15(B) for SN 1885A, which has a range of possible peak from the published light curve figure) as well as direct ob- magnitudes (solid line, connecting the △ signs). The specific peak servations of its SNR in recent years (Fesen et al. 1989, magnitude and ∆m15(B) values for SNe SN 1885A are unknown, as no direct B-band data are available. Limiting values for the 2007). Unfortunately, quantitative spectroscopic mea- peak magnitude are obtained under two possible assumptions; (1) surements of SNe did not begin until a few decades after a constant B−V = 1.32 in the 15 days epoch post-peak (bottom tri- the detection of SN 1885A. Nevertheless, descriptions of angle), corresponding to the inferred values for SN 1885A assuming the SN colors, and more importantly several prism mea- a similar evolution as found for sub-luminous SN 1991bg-like SNe (van den Bergh 2002; Garnavich et al. 2004) and suggested from surements, provide some information on its spectral fea- visual observations (de Vaucouleurs & Corwin 1985), or (2) a con- tures. stant B −V ≃ 0 (top triangle) corresponding to a similar evolution to that observed for SN 2002bj. The photographic magnitudes of SN 1939B are taken to be similar to B-magnitudes (Bessell 2005). 2.3. SN 1939B Error bars are not shown; typical error bars on ∆m15(B) are ±0.1; typical error bars on peak luminosity are ±0.1, not including error SN 1939B, discovered in the elliptical galaxy NGC bars from distance modulus (typically ±0.15). Detailed discussion 4621 in the Virgo cluster, was one of the first SNe re- on the error estimates and the ∆m15(B)-peak luminosity relation ported by Zwicky. It was observed by several people can be found in Altavilla et al. (2004). who provided a detailed light curve (Fig. 2, based on data from Leibundgut et al. 1991, that combined several These historical SNe were previously suggested to orig- independent but consistent data sources; re-analysis of inate from helium-shell detonations (Chevalier & Plait some of the original plates by Schaefer & Girard 1999 1988), but were not considered in the study of SN 2002bj. was also found to be consistent). Although mentioned in In this study we would like to remedy this and give ad- the original discovery report (Zwicky & Minkowki 1939; ditional data to complement the picture. Adams & Seahes 1940), a spectrum taken by Minkowski The data for these SNe are of course less complete than three weeks post maximum was never published. With could be obtained with modern instruments. However, the invaluable help of the Carnegie Observatories archive we do not find evidence of any systematic errors or biases, we have been able to locate the original spectrum plate, suggesting that the data quality is sufficient to reveal the and present the reduced spectrum for the first time (Fig. intrinsic properties of the SNe. 3). In the absence of any documentation, we assume The SNe we discuss provide additional clues to the that the observational setup was similar to that described origin of the fastest evolving SNe. In particular, they in other contemporaneous works (e.g. Minkowski 1939, provide direct evidence for old, non-star forming envi- 1942). A trace was extracted from both the SN spec- ronments for at least some of these SNe (SNe 1885A and trum and the arc spectra by following a line of maximum 1939B), likely excluding young massive progenitors and light. The wavelength solution was fitted using Legendre pointing to low-mass, old progenitors such as WDs. They polynomials, yielding residuals of less than 1 Å. Because also provide further evidence for the low ejecta mass no flux calibration is available for the data, the general (much lower than the Chandrasekhar mass) of such SNe, shape of the spectrum was normalized to most resemble and possibly the importance of other elements and/or the continuum of SN2004eo at 21 days (Fig. 3), enabling processes beside 56Ni decay in driving the SN light curves a direct comparison of their spectral features. at early epochs. In addition, the low ejecta mass inferred for SN 1885A is consistent with the low X-ray luminos- 3. PROPERTIES ity of the observed SN remnant (SNR) in M31, as we 3.1. Light curves will show, thus providing a light-curve-independent con- 3

−20 10 SN1885A (V) Regular SN Ia template (V) −19 9 SN2002bj (upper limits, unfilterred) SN2002bj (unfilterred) SN 2002bj +7 −18 SN 2002bj (B,V) 8 SN 1939B (B) He I λ5876 He I λ6678 −17 SN 1991bg (V) 7 SN 2008D +19 −16 6

−15 5 SN 1994I +18 −14 4 SN 1939B +22

Absolute M (mag) −13 3 Normalized Flux + Constant −12 2 SN 2004eo +21 −11 1 −10 −20 5 30 55 80 105 130 155 180 205 230 0 Time [day] 4000 4500 5000 5500 6000 6500 7000 7500 Rest Wavelength (A)˚

−20 Fig. 3.— The spectrum of SN1939B, 22 days after peak. It is similar to both the type Ia SN 2004eo (from Pastorello et al. 2007) and the type Ic SN 1994I (Filippenko et al. 1995) at a sim- −19 ilar phase. The type Ib SN 2008D (from Mazzali et al. 2008) is shown for comparison, with helium lines (blue-shifted by 10,500 −18 km/s) marked. SN 1939B shows no clear evidence for helium, with an apparent absence of the weak He I 6678 absorption. Never- theless, the feature possibly corresponding to Na D at 5890Å + −17 5896Å could alternatively be interpreted as a 5876Å He line. The spectrum of SN 2002bj at an earlier epoch (7 days) is also shown −16 (from Poznanski et al. 2010; after continuum subtraction).

Absolute M (mag) −15 days. All three SNe decay fast, faster than SN 1991bg- like sub-luminous SNe, with ∆m15(B) of ∼ 2.2, 2.6 and −14 ∼ 3.2 for SNe 1885A, 1939B and 2002bj, respectively, where we assume a ∼constant B − V for SN 1885A dur- −13 −15 −10 −5 0 5 10 15 20 25 30 35 40 ing the 15 days post-peak (as seen for SN 2002bj, and Time [day] inferred for SN 1885A, e.g. de Vaucouleurs & Corwin 1985; van den Bergh 2002), as only V-band data are Fig. 2.— The absolute photographic (pg)/B-band and visual available for this SN. (pv)/V-band light curves of SNe 1885A (pv), 1939B (pg) and 2002bj (B; V band is almost coincident with B). Both full light The light curves of SNe 1885A and 1939B were fol- curves (a) and zoomed in (b) light curves are shown. Note the lowed for a long period and show a break at ∼ 60 similar light curve shapes of these 3 events, declining rapidly from and ∼ 20 days, respectively, with a power-law-like de- a bright peak, in contrast to both normal and faint/fast SNe Ia. cay afterward6. The slope of SN 1885A at this epoch The light curves of SNe 1885A and 1939B are visual light curves, as − compiled by de Vaucouleurs & Corwin (1985, extracted from their (0.028 0.029 mag per day) is similar to that observed published light curve figure) and Leibundgut et al. (1991) for SNe for type Ia SNe at late epochs (see the non-stretched SN 1885A and 1939B, respectively. The light curve of SN 2002bj (V Ia template light curve in Fig. 2, and the light curve of band, and unfiltered points) is taken from Poznanski et al. (2010). SN 1991bg). The light curve of SN 1939B, however, de- Error bars for other data points are unknown, besides those esti- − mated by de Vaucouleurs & Corwin for the first data points of SN cays more slowly after the break (0.022 0.024 mag per 1885A. We show for comparison the V-band light curve of a regu- day). lar non-stretched SN Ia template (Nugent et al. 2002) normalized The late time exponential decay behavior of SN 1885A to have a peak magnitude of -19. Also shown is the V-band light may suggest a similar origin as for regular type Ia SNe at curve of SN 1991bg (using data from Leibundgut et al. 1993), the 56 prototypical example for faint and fast evolving type Ia SNe. this epoch, namely the radioactive decay of Co. The re- semblance between SNe 1885A and 1991bg at this epoch The light curves of SNe 1885A, 1939B and 2002bj (as is in stark contrast to their early behavior, showing a large difference between their peak luminosities. In type well as SN 1991bg) are shown in Fig. 25. All three light curves peak at luminosities typical of type Ia SNe, but 6 Note that some photographic light curves of SNe show a lev- show a very fast rise to maximum and fast decay. The eling at late epochs, possibly due to background contamination upper limit for the rise time of SN 2002bj is < 7 days; (Leibundgut et al. 1991). However, when discussing this effect Lei- SN 1885A rises at quite a similar timescale of ∼ 6 days. bundgut et al. (1991) specify several such leveling SNe, but do not ∼ − consider SN 1939B to be one of these. Moreover, SN 1939B does The rise time of SN 1939B is somewhat longer, 7 8 not level out to some possible background level as might be ex- pected by an interpretation of background light [as indeed seen 5 Data from the NASA/IPAC extragalactic database (NED) in the leveling of SNe specifically mentioned by Leibundgut et al. were used for distance modulus and Galactic extinction (based on (1991) in that context], but continues to fainter magnitudes, e.g. Schlegel et al. 1998) of the host galaxies of the SNe; distance mod- the latest point of SN 1939B is more than two magnitudes fainter ulus error bars are of the order of ±0.15 mag, these are not shown. than the position of the light curve break. Similarly, SN 1885A The light curves are not corrected for extinction inside the host does not show signs of leveling, and continuously becomes much galaxies. fainter at late times after the break. 4

Ia SNe the peak luminosity is thought to be driven by less massive ejecta of SN 1939B, in which case the ob- the energy deposition from 56Ni decay, whereas the late served peak luminosity of SN 1939B is even lower than decay is thought to be driven by its product (56Co). it would have been if all the γ-rays were trapped. The At early epochs the γ-rays (emitted as 56Ni decays to inconsistency between the peak and the tail ratios for SN 56Co and then to 56Fe) are still trapped in the ejecta 1939B and SN 1991bg which we discuss might therefore and deposit their energy there through successive Comp- be even larger than we calculated. We conclude that one ton scatterings. At later times, however, and especially requires a different/additional mechanism to drive the for the case of low ejecta mass (as is likely the case for light curve at an early epoch. One possible explanation the fast evolving SNe) the column density of the ejecta is that SN 1885A had other radioactive elements with may may not be thick enough to trap the γ-rays and short lifetime which contributed to the early light curve they may escape without depositing their energy into generation (e.g. Shen et al. 2010; Waldman et al. 2010). the ejecta. In principle, this could explain the difference Note that an additional data point at day 20 (post- between the luminosities of the SNe at different epochs, maximum) exists for SN 2002bj in bands other than V as SN 1885A seems to have much lower mass of ejecta. (see figure 1 in Poznanski et al. 2010), possibly suggest- However, when estimating the difference in the fraction ing a steeper slope than SN 1885A and SN 1939B at of γ-rays absorbed inside the ejecta we find that it could this epoch. This is seen in the B and R bands, which only account for a small fraction of the observed differ- likely bracket the V band to a similar magnitude. The ence, as we now show. I band behavior, however, shows a less steep slope, and The electron (Thompson) optical depth is: is quite unique; we see an R − I ≃ 1.2 for SN 2002bj at this epoch, which is much larger than found for any 3 Mej Z σT type of SNe (compare with Poznanski et al. 2002). The τe = neRσT = 2 color evolution of SN 1885A and SN 1939B is unknown, 4π mp A R − as no filtered data are available. Color descriptions sum- M t 2 ∼ 9 ej (1) marized by de Vaucouleurs & Corwin (1985) suggest a 0.16M⊙  15 day value of B − V ∼ 1.3 for SN 1885A at peak, much larger than the B − V ∼ 0 observed for SN 2002bj. where Z is the atomic number, A is the mass number, mp is mass of proton, σT is the Thompson cross-section 3.2. SN 1885A remnant and R ∼ 6(t/day) AU is the radius at time t. There- fore the optical depth at the epoch of peak luminos- SN 1885A is one of the very rare cases where both the ity is sufficient to trap most of the γ-rays. Assuming SN and its SNR could be observed. After several failed a similar energy source from radioactive decay, we ex- attempts to find optical, radio, and X-ray emission at the pect the ratio between the luminosities of SNe 1885A location of SN 1885A, the SNR was finally discovered in and SN 1991bg to reflect the ratio between their 56Ni absorption against the bulge of M31 (Fesen et al. 1989) and later studied in detail using HST (Fesen et al. 2007). ejecta mass, MNi. Next, following Kasliwal et al. (2010) we use a fitting formula (as given in Kulkarni 2005; Equa- The lack of direct emission from this 125 yr old SNR tion 47) to estimate the fraction of γ-rays that are effec- is remarkable. More than a decade of monitoring X- tively absorbed inside the ejecta, η(τ ). At late times ray binaries in the central region of M31 with Chandra e (Williams et al. 2006; Li et al. 2009) has only yielded an (100 days) we find η to be ∼ 0.011 M⊙ and ∼ 0.057 ∼ × 34 ⊙ ⊙ upper limit on the X-ray luminosity of 6 10 erg M for SN ejecta of 0.2 M (corresponding to the mass −1 estimate derived below for SN 1885A) and 1.4 M⊙ (for s from the SNR, assuming a Raymond-Smith model a typical type Ia SN), respectively. At this epoch the with kbT=0.4 KeV (Zhiyuan Li, private communication, radiation is dominated by the energy deposition from 2010). the radioactive decay of 56Co rather than the 56Ni de- The thermal X-ray emission from a young SNR is de- cay (which has a half life time of only ∼ 6 days), and termined mainly by the density of the ambient medium therefore Lrad ≃ (0.965η +0.035)LCo where LCo is the and the mass and kinetic energy of the SN ejecta (e.g., radioactive power released by the decay of 56Co (fol- Hamilton & Sarazin 1984; Badenes et al. 2003). Li et al. lowing Kasliwal et al. 2010, where the kinetic energy of (2009) have analyzed the properties of the diffuse X-ray emission in the central arcminute of M31 (which includes positrons is 3.5 % of LCo; Sollerman et al. 2002). We find the position of SN 1885A), and find an average density of that the difference in γ-rays absorption between the two −3 SNe changes their relative luminosities by only a factor ∼ 0.1 cm . This gas density is low, as befits the bulge of (0.965η(1885A)+0.035)/(0.965η(1991bg) + 0.035) ≃ of a spiral galaxy, but only an order of magnitude lower 0.045/0.09 = 0.5. As can be seen in Fig. 2 the ratio than the ambient medium density around the Tycho SNR between the SN luminosities at peak vs. late epoch can (Badenes et al. 2006), and only a factor 2-3 lower than not be accounted for by this factor. Even if the effective around the remnant of SN 1006 (Raymond et al. 2007), optical depth is ten times larger than the simple Thomp- which are prototypical Type Ia SNRs. We used the code son optical depth we use (e.g. Swartz et al. 1995 find described in Badenes et al. (2003, 2005) to calculate the few to ten times larger depth) the ratio we find changes thermal X-ray luminosity of SNR models with an age only by order unity (0.4 instead of 0.5). We note that of 125 yrs, assuming different ambient medium densities the estimated Thompson optical depth we use is calcu- and SN ejecta configurations. For all canonical Type Ia ∼ ≃ lated at the bottom of the ejecta, where as not all of the SN ejecta (i.e., explosions with Mejecta 1.4 M⊙ MCh γ-rays is emitted from there. It is therefore possible that and kinetic energy of at least a few ×1050 erg; see Fig. a large fraction of the γ-rays may not be trapped even 4 for the specific models used), our calculations predict at peak luminosity. This issue is more likely to affect the X-ray luminosities in excess of 1038 ergs s−1 at the value 5

1042 b30_3d_768 seems to resemble both that of regular type Ia and type DDTe Ic SNe, e.g. SN 2004eo (Ia) and SN 1994I (Ic; see fig. 3 1040 SCH for comparison) as well as faint SNe such as SN 1986G

] (spectrum not shown). We find no evidence of hydrogen, -1 and no clear evidence for helium lines in the spectrum. 1038 Naturally, the recently studied SN 2002bj has much better spectral observations than the historical SNe. Pos- sible helium lines are seen in its spectra and may show 1036 Density Estimate (Li et al.) it to be a type Ib SN, however, Kasliwal et al. (2010) X-ray Upper Limit suggest a different interpretation, making this SN more Luminosity [erg s 34 similar to a type Ia or a type Ic SN. Note also that this 10 SN shows other similarities with type Ia SNe, mostly in Allowed Region its absorption feature near 6150 Å (rest frame), usually 1032 attributed to Si II (Poznanski et al. 2010). 10-25 10-24 10-23 AM Density [g cm-3] 3.4. Inferred masses Calculations based on the fast light curve of SN 2002bj (Poznanski et al. 2010), resulted in ∼ 0.15 M⊙ of mass in the ejecta. We can use a similar treatment, follow- Fig. 4.— Predicted thermal X-ray luminosities from the shocked ing Arnett (1982), by relating the mass, kinetic energy, ejecta in several SNR models at an age of 125 yr, plotted as a func- and the timescale of the SN light curve, to estimate the tion of ambient medium density. The models are: b30_3d_768 ejecta mass in SN 1939B. The ejecta mass of a given su- . . × 51 (from Travaglio et al. 2004; Mej = 1 4 M⊙; Ek = 0 5 10 pernova can be estimated using its light curve and the ergs; red line), delayed-detonation (DDTe; from Badenes et al. 51 observed ejecta velocities. The expansion velocity, v, of 2003; M = 1.4 M⊙; Ek = 1.0 × 10 ergs; blue line), and sub- ej a SN is proportional to (E /M )1/2, where E is the Chandrasekhar (SCH; also from Badenes et al. 2003; Mej = 1.0 kin ej kin 51 M⊙; Ek = 1.0 × 10 ergs; orange line). The horizontal line is kinetic energy and Mej is the ejected mass, while the typ- ∝ 3 1/4 the upper limit on the X-ray luminosity of the SNR (Zhiyuan Li, ical duration of a SN light curve is td (Mej/Ekin) private communication, 2010). The vertical dashed line is the es- timate of the density in the circumnuclear region of M31, which (Arnett 1982). Combining these equations and assuming includes the location of SN 1885A from Li et al. (2009). that two objects have the same opacity, we have 3 2 of the ambient medium density estimated by Li et al. v1 t1 37 −1 Ekin 1/Ekin 2 = (2) (2009), and well above 10 ergs s even for densities a , , v  t  factor two lower. Such high luminosities are inconsistent 2 2 with the observed X-ray upper limit, even after correct- and 2 ing for absorption. Only Type Ia SN models with signif- v1 t1 icantly lower ejected masses, such as sub-Chandrasekhar Mej,1/Mej,2 = . (3) v t  explosions, might be able to reproduce X-ray fluxes low 2 2 enough to match the observations. It is of course possible Given the spectral similarity of SN 1939B (line widths) that SN 1885A is in a local ’bubble’ with a density much to SNe 1986G and 2004eo, we take the photospheric ve- lower than the average value found by Li et al. (2009), locities to be similar to these. We adopt a velocity of 104 but taken at face value, the lack of direct emission from km s−1 (given the similarity to SN 1986G and SN 2004eo, the remnant of SN 1885A seems hard to reconcile with which were found to have typical SN Ia photospheric ve- the possibility that its birth event was a canonical Type locities of 9.36×103 and 10.52×103 km s−1, respectively). Ia SN, and suggests much lower ejecta mass and energy. The timescale for the rise and fall of SN 1939B is com- parable to that of SN 2002bj, although slightly longer 3.3. Spectra (where Poznanski et al. 2010 find the rise time to be 3 The available spectra of the fast-evolving SNe are times shorter than typical SNe Ia). Taking a rise time of shown in Fig. 3. A good quantitative spectrum of trise ∼ 7 days, and comparing these values with an aver- SN 1885A is not available. A historical reconstruc- age type Ia SN with tIa−rise ∼ 17.5 days, and taking the tion based on observations made with prisms is de- following typical velocities of 104 km s−1 (Hayden et al. scribed by de Vaucouleurs & Corwin (1985). The spec- 2010), total ejecta mass of 1.4 M⊙ and total energy of tral data suggest that it is not a type II SN, but given 1051 erg, we find the ejecta mass of SN 1939B to be 50 its poor quality, it cannot be used to conclusively clas- ∼ 0.22 M⊙ with a kinetic energy of 1.6 × 10 (where sify SN 1885A as a specific type I SN. We do note that similar opacities are assumed). No data are available for de Vaucouleurs & Corwin (1985) and Pastorello et al. the photospheric velocities of SN 1885A and similar di- (2008) mentioned the possibility of a type Ib/c identifica- rect estimate can not be done. However, if we take the tion (the latter suggesting SN 1885A was a massive star expansion velocity of its SNR, 1.24 ± 0.14 × 104 km s−1 exploding in a He-rich circumstellar medium, though the (Fesen et al. 2007) to be the typical ejecta velocity, we stellar environment of this SN is likely old). However, find the total ejecta mass to be ∼ 0.2 M⊙ and the kinetic detailed analysis by Wheeler & Harkness (1990) argue energy is of 2.2 × 1050 erg (taking a similar rise time as against any evidence for helium lines. We conclude that SN 2002bj of ∼ 6 days); this is consistent with our analy- SN 1885A most resembles a type Ia or a type Ic SN. sis of the SNR in section 3.2, showing that the energetics The spectrum of SN 1939B, taken ∼ 22 days after peak, and ejecta mass are likely below those of canonical type 6

Ia SNe. We conclude that all three SNe seem to have a TABLE 1 low ejecta mass, 0.1 − 0.3 M⊙ and low kinetic energies, Fraction rate of bright fast evolving SNe 1 − 3 × 1050 erg. Fast evolving Total Fraction (%) Ref. 3.5. Environments 1 (1939B) 5 20 ± 20 Zwicky 1942 1 (1939B) 30 3.3 ± 3.3 Katgert & Oort 1967 Both SNe 1885A and 1939B are found in old environ- 1 (1885A) 45∗1 2.2 ± 2.2 Tammann 1970 ments. Multi-band photometric observations of the en- 2 (1885A, 1939B) 96 2.1 ± 1.4 Tammann et al. 1994 vironment of SN 1885A (Olsen et al. 2006) in the bulge 0 110 < 0.9 ± 0.9 Cappellaro et al. 1997 of M31 found it to be composed of an old stellar popu- 1 (2002bj) 103∗∗ 1 ± 1 Poznanski et al. 2010 ∼ lation, with a characteristic age of 10 Gyr. SN 1939B ∗The sample is only of Sb and Sc galaxies occurred in the elliptical galaxy NGC 4621. HST im- ∗∗They reported SN 2002bj relative to 31 type Ia SNe, the total ages of the galaxy (WF4 detector, with the F555W and number of SNe, 103, is estimated given the fractional rate of type ∼ . F814W filters and an exposure time of 1050 s, taken on Ia SNe to be 0 3 (Li et al. 2010). th Feb 5 , 1995) and of the location of SN 1939B were used rapid decay. These results point to a low mass (

McWilliam, Doug Mink and Bob Kirshner for helping CfA and BIKURA(FIRST) fellow. A.G. is supported by with locating, reducing and analyzing of the original grants from the Israeli Science Foundation, an EU FP7 spectroscopic plate of SN 1939B. We would also like to Marie Curie IRG Fellowship, a research grant from the thank the anonymous referee, and the second referee, Peter and Patricia Gruber Awards, the Weizmann Min- David Branch, for their helpful comments. HBP is a erva program and the Benoziyo Center for Astrophysics.

REFERENCES Adams, W. S. & Seahes, F. H. 1940, Mount Wilson Observatory Livne, E. & Arnett, D. 1995, ApJ, 452, 62 Annual Report, 12, 1 Minkowski, R. 1939, ApJ, 89, 156 Altavilla, G. et al. 2004, MNRAS, 349, 1344 —. 1942, ApJ, 96, 199 Arnett, W. D. 1982, ApJ, 253, 785 Nugent, P., Kim, A., & Perlmutter, S. 2002, PASP, 114, 803 Badenes, C. et al. 2003, ApJ, 593, 358 Olsen, K. A. G. et al. 2006, AJ, 132, 271 —. 2005, ApJ, 624, 198 Pastorello, A. et al. 2008, MNRAS, 389, 113 —. 2006, ApJ, 645, 1373 Perets, H. B. et al. 2010a, Nature, 465, 322 Benetti, S. et al. 2005, ApJ, 623, 1011 —. 2010b, arXiv:1012.0570 Bessell, M. S. 2005, ARA&A, 43, 293 Poznanski, D. et al. 2002, PASP, 114, 833 Bildsten, L. et al. 2007, ApJL, 662, L95 —. 2010, Science, 327, 58 Cappellaro, E. et al. 1997, A&A, 322, 431 Raymond, J. C. et al. 2007, ApJ, 659, 1257 Chevalier, R. A. & Plait, P. C. 1988, ApJ, 331, L109 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, de Vaucouleurs, G. & Corwin, Jr., H. G. 1985, ApJ, 295, 287 525 Fesen, R. A., Saken, J. M., & Hamilton, A. J. S. 1989, ApJ, 341, Shen, K. J. et al. 2010, ApJ, 715, 767 L55 Sollerman, J., et al. 2002, A&A, 386, 944 Fesen, R. A. et al. 2007, ApJ, 658, 396 Swartz, D. A., Sutherland, P. G., & Harkness, R. P. 1995, ApJ, Garnavich, P. M. et al. 2004, ApJ, 613, 1120 446, 766 Schaefer, B. E., & Girard, T. M. 1999, Anni Mirabiles, 69 Tammann, G. A. 1970, A&A, 8, 458 Hamilton, A. J. S. & Sarazin, C. L. 1984, ApJ, 284, 601 Tammann, G. A., Loeffler, W., & Schroeder, A. 1994, ApJS, 92, Hayden, B. T. et al. 2010, ApJ, 712, 350 487 Howell, D. A. 2001, ApJ, 554, L193 Travaglio, C. et al. 2004, A&A, 425, 1029 Iben, Jr., I. & Tutukov, A. V. 1991, ApJ, 370, 615 van den Bergh, S. 2002, AJ, 123, 2045 Iben, I. J. et al. 1987, ApJ, 317, 717 Van Dyk, S. D. et al. 1999, AJ, 118, 2331 Kasliwal, M. M. et al. 2010, ApJ, 723, L98 Waldman, R. et al. 2010, arXiv:1010.1522. Katgert, P. & Oort, J. H. 1967, Bull. Astron. Inst. Netherlands, Wheeler, J. C. & Harkness, R. P. 1990, Reports on Progress in 19, 239 Physics, 53, 1467 Kulkarni, S. R. 2005, ArXiv:astro-ph/0510256 Williams, B. F. et al. 2006, ApJ, 643, 356 Leibundgut, B. et al. 1991, A&AS, 89, 537 Woosley, S. E., Taam, R. E., & Weaver, T. A. 1986, ApJ, 301, 601 —. 1993, AJ, 105, 301 Woosley, S. E., & Kasen, D. 2010, arXiv:1010.5292 Li, W. et al. 2010, arXiv:1006.4613 Zwicky, F. & Minkowki, R. 1939, IAU Circ., 774, 1 Li, Z., Wang, Q. D., & Wakker, B. P. 2009, MNRAS, 397, 148 Zwicky, F. 1942, ApJ, 96, 28