arxiv: v2 [astro-ph] 1 Jul 2009

Transcription

1 DRAFT VERSION NOVEMBER 1, 2018 Preprint typeset using LATEX style emulateapj v. 11/26/04 A NEW CLASS OF LUMINOUS TRANSIENTS AND A FIRST CENSUS OF THEIR MASSIVE STELLAR PROGENITORS TODD A. THOMPSON 1,2,3, JOSÉ L. PRIETO 1,2, K. Z. STANEK 1,2, MATTHEW D. KISTLER 2,4, JOHN F. BEACOM 1,2,4, & CHRISTOPHER S. KOCHANEK 1,2 Draft version November 1, 2018 arxiv: v2 [astro-ph] 1 Jul 2009 ABSTRACT The progenitors of SN 2008S and the 2008 luminous transient in NGC 300 were deeply dust-enshrouded massive stars, with extremely red mid-infrared colors and relatively low bolometric luminosities ( L ). The transients were optically faint compared to normal core-collapse supernovae, with peak absolute visual magnitudes of 13 M V 15, and their spectra exhibit narrow Balmer and [Ca II] emission lines. These events are unique among transient-progenitor pairs and hence constitute a new class. Additional members of this class may include the M85 transient, SN 1999bw, 2002bu, and others. Whether they are true supernovae or bright massive-star eruptions, we argue that their rate is of order 20% of the core-collapse supernova rate in star-forming galaxies. This fact is remarkable in light of the observation that a very small fraction of all massive stars in any one galaxy, at any moment, have the infrared colors of the progenitors of SN 2008S and the NGC 300 transient. We show this by extracting mid-infrared and optical luminosity, color, and variability properties of massive stars in M33 using archival imaging. We find that the fraction of massive stars with colors consistent with the progenitors of SN 2008S and the NGC 300 transient is In fact, only 10 similar objects exist in M33 (and perhaps 1) all of which lie at the luminous red extremum of the AGB sequence. That these transients are simultaneously relatively common with respect to supernovae, while their progenitors are remarkably rare compared to the massive star population, implies that the dust-enshrouded phase is a short-lived phenomenon in the lives of many massive stars. This shrouded epoch can occur only in the last 10 4 yr before explosion, be it death or merely eruption. We discuss the implications of this finding for the evolution and census of low-mass massive stars (i.e., 8 12 M ), and we connect it with theoretical discussions of electron-capture supernovae near this mass range. Other potential mechanisms, including the explosive birth of massive white dwarfs, and massive star outbursts are also discussed. A systematic census with (warm) Spitzer of galaxies in the local universe (D 10 Mpc) for analogous progenitors would significantly improve our knowledge of this channel to massive stellar explosions, and potentially to others with obscured progenitors. Subject headings: supernovae: general, individual (SN 2008S, 1999bw, 2002bu) surveys: stars, galaxies stars: evolution 1. INTRODUCTION Identifying the progenitors of core-collapse supernovae, the outbursts of Luminous Blue Variables (LBVs), and other massive-star transients is essential for understanding the physics, demographics, variability, evolution, and end-states of massive stars. The problem of identifying the progenitors of bright transients from massive stars is difficult, and traditionally limited to serendipitous archival imaging of nearby galaxies in the optical and near-infrared (e.g., Van Dyk et al. 2003; Smartt et al. 2004; Li et al. 2007; see the extensive summary in Smartt et al. 2009). Progenitor searches are complemented by statistical studies of supernova environments within their host galaxies, which provide indirect evidence for associations between certain types of supernovae and broad classes of progenitors (e.g., James & Anderson 2006; Kelly et al. 2008; Prieto et al. 2008; Anderson & James 2008). A much more direct method for understanding the relation between types of massive stars and their transients is to cata- 1 Department of Astronomy, The Ohio State University, 140 W. 18th Ave., Columbus, OH 43210; thompson, prieto, kstanek, ckochanek@astronomy.ohio-state.edu. 2 Center for Cosmology & Astro-Particle Physics, The Ohio State University, 191 W. Woodruff Ave., Columbus, OH Alfred P. Sloan Fellow 4 Department of Physics, The Ohio State University, 191 W. Woodruff Ave., Columbus, OH 43210; kistler, beacom@mps.ohio-state.edu. log all of the massive stars in the local universe (D 10 Mpc) before explosion. While surveys for bright optical transients in the local universe are well-developed (e.g., Li et al. 2001), a fairly complete census of the massive stars in nearby galaxies has only recently been proposed and undertaken (Massey et al. 2006; Kochanek et al. 2008). Despite the technical challenges required by the depth, area, and cadence of the observations, these surveys are critical for our understanding of the one-to-one correspondence between massive stars and their end-states, whether they are successful or failed explosions (e.g., Kochanek et al. 2008). The long-term promise of these surveys is to produce a catalog within which the characteristics of progenitors of future supernovae are listed preexplosion, as in the case of SN 1987A (West et al. 1987; Menzies et al. 1987). They will provide an essential mechanism for understanding the direct causal mapping between individual progenitor types and their transients (see, e.g., Gal-Yam et al. 2007). Here, we describe a new link in this causal mapping. The discovery by Prieto et al. (2008a) and Prieto (2008) of the dust-obscured progenitors of the luminous outbursts in NGC 6946 (SN 2008S; Arbour & Boles 2008) and in NGC 300 (Monard 2008) with Spitzer, opens up qualitatively new possibilities in the study of the connection between massive stars and their explosions. We show that these discoveries allow us to make a strong and perhaps unprecedented

2 2 Thompson et al connection between a dust-enshrouded sub-population of massive stars and a new relatively common class of bright transients. The argument presented in this paper can be summarized in four points: 1. The transients SN 2008S and NGC constitute a class. Both have peak absolute V -band magnitude M V 14±1 ( 2 3 magnitudes fainter than normal core-collapse supernovae; e.g., Richardson et al. 2002), strong evidence for internal extinction in their spectra, narrow emission lines (similar to low-luminosity Type IIn supernovae; e.g., Filippenko 1997), and progenitors that are optically obscured and deeply dust-embedded (dust-reprocessed emission giving blackbody temperatures 500 K), and with bolometric luminosities of L. In addition, they show little infrared variability on a few-year baseline. The details of this unique class of progenitor-transient pairs and its members, as well as a comparison with other classes of optical transients are presented in 2. An in-depth search for analogs to the progenitors of SN 2008S and NGC 300 in M33 is presented in Transients of this type are relatively common with respect to core-collapse supernovae. A total of 22 core-collapse supernovae or supernova-like transients have been discovered within 10 Mpc since Sixteen were confirmed supernovae, two were LBV eruptions, one was a Type IIn supernova (SN 2002bu in NGC 4242, D 8 Mpc; Puckett & Gauthier 2002), whose relatively low peak magnitude (M V 15; Hornoch 2002) suggests a close similarity with the remaining two bright transients (see 2.5), which are of primary interest in this paper: SN 2008S (D 5.6 Mpc) and NGC 300 (D 1.9 Mpc), whose physical nature is uncertain. They may be either true (but optically sub-luminous) supernovae, or a new class of massive star eruptions. Taken at face value, these numbers imply that the rate of SN 2008S-like transients is of order 10% of the supernova rate. 6 Because of incompleteness, the true rate is likely higher. We discuss the frequency of these events in detail in The progenitors of this class are extremely rare among massive stars at any moment, in any star-forming galaxy. Although the bolometric luminosities of the progenitors of SN 2008S and NGC 300 are unremarkable for massive stars ( L ), their colors put them in a class consisting of less than 10 4 of all massive stars ([3.6] [4.5]µm color 2.0 and 2.7, respectively). In a mid-infrared (MIR) color-magnitude diagram (see Figs. 1 & 2), these progenitors lie at the extremum of the AGB sequence in both luminosity and MIR color. We refer to them as extreme AGB (EAGB) stars throughout this work. Because of their relatively low bolometric luminosities, they are not η Carinae, cool hypergiant, or classical LBV analogs (see, e.g., Smith 2008). They are thus distinct from the supernova impostors, produced by bright outbursts of optically-luminous LBVs (e.g., SN 1997bs; Van Dyk et al. 2000, 2003). In 4 we present results from a comprehensive survey of M33 for massive stars with 5 Throughout this paper we denote the transient in NGC 300 as NGC 300 (e.g., the progenitor of NGC 300 ) unless we specifically refer to the host galaxy. 6 Throughout this paper we use supernova rate to mean the core-collapse supernova rate unless we specifically mention the contribution from Type Ia supernovae. properties similar (in bolometric luminosity, obscuration, and variability) to the progenitors of SN 2008S and NGC 300. We find remarkably few. We compare with MIR surveys of the LMC, SMC, and NGC Conclusion: A large fraction of all massive stars undergo a dust-enshrouded phase within 10 4 yr of explosion. This is the most natural explanation for the facts of points (2) and (3) above. If these transients have a rate comparable to the supernova rate ( 20%; point 2), then the timescale for the obscured phase is determined by the ratio of the number of dust-obscured massive stars relative to the entire population ( 10 4 ; see 4, 5) times the average lifetime of a massive star ( 10 7 yr). Importantly, from the rarity of SN 2008Slike progenitors alone ( 10 4 of massive stars; point 3), one would naively expect a comparable fraction of supernovae to have progenitors of this type, if the dust-obscured phase occurs at a random time in the life of a massive star. However, the relative frequency of these explosions (point 2) shows that this phase must come in the last 10 4 yr, just before explosion. Thus, there must be a causal relation between the occurrence of the highly dust-enshrouded phase and eruption. These points, together with a discussion of the theory of the evolution of massive stars, the potential connection with electron-capture supernovae, white dwarf birth, and other hypotheses for the physical mechanism of 2008S-like explosions, as well as a call for a more comprehensive Spitzer survey for analogous sources within D 10 Mpc, are presented in THE CLASS We start by listing the objects we view as likely to represent this new class of SN 2008S-like transients and progenitors. The two objects that define the class SN 2008S ( 2.1) and NGC 300 ( 2.2) are unique among transient-progenitor pairs. The progenitors are relatively low-luminosity, have low variability, and are deeply dust-embedded on 100 AU scales. The transients are low-luminosity, with spectra exhibiting both narrow Balmer lines (similar to low-luminosity IIn s and impostors), and [Ca II] in emission, and have rapidly decaying lightcurves compared to IIP supernovae. The transient in M85 ( 2.3) and SN 1999bw ( 2.4) may also be members of the same class, but we cannot confirm the existence of a dust-obscured progenitor similar to SN 2008S/NGC 300. In 2.5, we contrast this class with other peculiar outbursts, such as the supernova impostors and low-luminosity Type IIP supernovae, and we note a number of other transients that are not excluded as members of the SN 2008S-like class SN 2008S SN 2008S in NGC 6946 (D 5.6 Mpc; Sahu et al. 2006) was discovered February 1.79 UT (Arbour & Boles 2008). Because of the presence of narrow Balmer lines (FWHM 1000 km s 1 ) it was initially classified as a young Type IIn supernova. Stanishev et al. (2008) reported strong Na D absorption with equivalent width of 5 Å, indicating a high degree of internal extinction. On the basis of its relatively low luminosity (M V 14) and peculiar spectrum (including the presence of strong and narrow [Ca II] 730 nm doublet in emission), Steele et al. (2008) proposed that SN 2008S was a supernova impostor such as SN 1997bs (Van Dyk et al. 2000). Botticella et al. (2009) present the late-time lightcurve of SN 2008S, which shows evidence for a power-law time dependence with

3 A New Class of Transients 3 a slope indicative of being powered by the radioactive decay 56 Co. Following the suggestion by Prieto et al. (2008, hereafter P08a) and the discussion presented in 5, Botticella et al. argue that SN 2008S may have been an electron-capture supernova. P08a used a deep pre-explosion archival Large Binocular Telescope (LBT) image of NGC 6946 to put stringent upper limits on the optical emission from the site of the supernova (3σ limits of M U > 4.8, M B > 4.3, and M V > 3.8). The failure to detect a progenitor in the optical led P08a to examine archival Spitzer IRAC imaging of NGC 6946 from the SINGS Legacy Survey (Kennicutt et al. 2003) for the progenitor. A point source at the location of SN 2008S was detected at 4.5, 5.8, and 8µm, but undetected at 3.6µm and 24µm, leading to a lower limit on the [3.6] [4.5] color of 2. The bestfitting blackbody temperature to the SED was 440 K. 7 The integrated luminosity was L, consistent with a relatively low-mass massive-star progenitor with a zero-age main sequence (ZAMS) mass of M 10 M. The progenitor SED is shown in the left panel of Figure 5. Simple arguments suggest that the obscuration was circumstellar, with an optical depth at visual wavelengths larger than a few on a physical scale of order 150 AU (P08a). The explosion itself was also serendipitously observed by Spitzer 5 days after discovery (Wesson et al. 2008), with an infrared luminosity 50 times larger than the progenitor, suggesting a substantial amount of dust-reprocessing (see Fig. 2 of P08c). The lack of variability of the progenitor on a 3 yr timescale argues that the obscuring medium was either a continuous wind with a steady photosphere or a (implausibly?) slowly ( 10 km s 1 ) expanding shell NGC 300 A luminous optical transient (M V 13) in the nearby galaxy NGC 300 (D 1.9 Mpc; Gieren et al. 2005) was discovered by Monard (2008), and reported by Berger & Soderberg (2008). The latter used archival imaging with ACS/WFC onboard the Hubble Space Telescope (HST) to put very tight limits on the optical emission (see Fig. 5). These limits led them to suggest that the NGC 300 transient was analogous to the outburst V838 Monocerotis (Bond et al. 2003). Similar to the case of SN 2008S, Prieto (2008, hereafter P08b) discovered a deeply-embedded source in archival Spitzer imaging (PI: R. Kennicutt) at the location of the transient. The source was detected in all IRAC bands, as well as in MIPS 24µm, and had a [3.6] [4.5] color of 2.7. The SED implies a 330 K blackbody (as with SN 2008S, a blackbody is a fairly poor fit to the SED) with a bolometric luminosity of L. This finding confirms the massive stellar origin of the NGC 300 transient, and is consistent with a relatively low-mass massive star (see Fig. 5). Importantly, depending on the details of stellar models for ZAMS masses in the range of 10 M, at fixed final bolometric luminosity, the inferred initial progenitor mass may be multiply-valued and a luminosity of L can imply a 5, 8 or 11 M progenitor (see Smartt et al. 2009, their Fig. 2; 5.2). 8 The luminosity and blackbody temperature of the progenitor of NGC 300 suggest an obscuring medium with physical 7 A blackbody provides a rather poor fit to the SED, perhaps indicating an interesting grain size distribution in the obscuring material (P08a). 8 P08b assumed standard singly-valued masses, and given the measured luminosity of the progenitor of NGC 300 quoted ( 10 5 L ) they inferred an initial progenitor mass of M. scale of order 300 AU. The deep limits on the optical emission from the progenitor with HST suggest an optical depth at V considerably larger than unity ( 8 10). 9 The fact that the transient in NGC 300 and SN 2008S were similar, both in their luminosities (both relatively faint with respect to typical supernovae with M V 14) and spectra (with narrow Balmer lines and strong [Ca II] 730 nm doublet in emission; Bond et al. 2008), and that their progenitors were similar (highly dust-obscured, relatively modest bolometric luminosities) led P08b to propose that NGC 300 and SN 2008S share a common origin: the explosion whether supernova or outburst of a massive star enshrouded in its own dust. As we show below, the fact that their progenitors are so rare among massive stars implies that just two events (although, see 2.3 & 2.4) are sufficient to define a class M85 The optical transient in the Virgo galaxy M85 (NGC 4382) was discovered in early 2006 by the Lick Observatory Supernova Search Team (KAIT; Li et al. 2001) and discussed in Kulkarni et al. (2007) and Pastorello et al. (2007) and may also be a member of the class defined by SN 2008S and NGC 300. The transient had peak R-band absolute magnitude of 12 with a day plateau, similar to lowluminosity Type IIP supernovae. Optical limits constrain the progenitor to be less than 7 M or highly obscured (Ofek et al. 2008; see also Pastorello et al. 2007). 10 The low optical luminosity, plateau, and redward spectral evolution of the transient led Kulkarni et al. (2007) to propose that it was analogous to the outburst V838 Mon. In contrast, Pastorello et al. (2007) argued that it was a low-luminosity Type IIP supernova (see 2.5). Prieto et al. (2008c, hereafter P08c) again searched archival Spitzer imaging and discovered an infrared source at the site of the optical transient, but, by chance, taken 8.8 days before the discovery in the optical by KAIT. The source was detected in all IRAC bands and was associated with the transient itself, and not the progenitor, since archival Spitzer imaging from 2004 (PID 3649; PI P. Côté) does not reveal a point source at the location of the M85 transient. Using these images, we derive 3σ upper limits of L and L at 3.6µm and 4.5 µm, respectively, for the progenitor. The transient was also detected in the mid-ir seven months after the initial discovery by Rau et al. (2007), the fluxes having decreased by a factor of 5 over that time. The bright infrared transient discovered by P08c is adequately fit by a blackbody with temperature of 800 K with luminosity L. The optical and NIR photometry of Kulkarni et al. (2007) indicate a second component to the SED with a blackbody temperature of 3900 K and with a luminosity of L. The cooler re-radiated dust emission arises from a region of order AU in physical scale. Assuming that the optical emission did not vary in the 8.8 days between the IR discovery and the optical discovery, the ratio of the power in these two blackbody components implies that the optical depth at V (τ V ) is less than unity. The physical scale of the obscuring medium indicates that it is likely circumstellar, and the result of a mass-loaded wind. In addition, the luminosity of the transient suggests that any pre-existing dust within 9 Assuming a spherical homogeneous medium and a Galactic dust-to-gas ratio, this optical depth at V requires a total mass of obscuring material of roughly M on a scale of 300 AU. 10 Ofek et al. (2008) quote a limiting absolute magnitude at F850LP(z) with HST of > 6.2.

4 4 Thompson et al AU would have been destroyed during the explosion. Given the fact that the optical depth to the source scales as r 1 in a freely-expanding wind, it is not implausible that τ V to the progenitor was a factor of larger before explosion. These estimates suggest a pre-explosion obscuring medium similar in its gross properties to the SN 2008S and NGC 300 progenitors. The M85 transient also showed narrow Balmer lines in emission, as well as some Fe II lines, similar to SN 2008S and NGC 300. Because of the strong evidence for obscuration of the progenitor, as evidenced by the bright IR transient, and the similarity of the spectra, P08c proposed that these outbursts share a common origin and that their obscured progenitors may give rise a new class of 2008S-like transients. We emphasize that because the character of the progenitor is not known (except for the optical and IR limits), the connection to SN 2008S- and NGC 300-like events is plausible rather than certain. Nevertheless, if the M85 transient was associated with an embedded massive star, the IR limits we derive are consistent with the luminosities derived for the 2008S and NGC 300 progenitors SN 1999bw The Lick Observatory Supernova Search reported in April 1999 the discovery of a possible supernova in the galaxy NGC 3198 (Li 1999). The optical spectrum of the transient, dominated by narrow Balmer lines in emission (Garnavich et al. 1999; Filippenko et al. 1999), and its low V-band absolute magnitude at maximum of 13 (D 13.7 Mpc; Freedman et al. 2001) led Li et al. (2002) to propose that this transient was an LBV-like outburst. Like SN 2008S and NGC 300, its spectrum showed [Ca II] in emission. Additionally, an infrared source coincident with the optical position of the transient was detected in archival Spitzer imaging obtained with IRAC by the SINGS Legacy Survey five years after the discovery of SN 1999bw (Sugerman et al. 2004). The source was detected in all IRAC bands and the SED was well-fit by a 450 K blackbody with an integrated luminosity of L, which translates into a blackbody scale of 300 AU. 11 We have checked archival IRAC images obtained in December 2005 (PID 20320; PI B. Sugerman), 1.5 yr after the detection in the SINGS images, and we confirm that the MIR source is indeed the transient, since the fluxes have declined by a factor of more than 3 during this time. The combination of a low optical luminosity at maximum, an optical spectrum dominated by narrow Balmer lines in emission, the presence of [Ca II] emission, and a luminous infrared emission detected with Spitzer, make SN 1999bw similar to SN 2008S, NGC 300, and the transient in M85. However, as in the case of M85, we emphasize that because there is no information on the progenitor we cannot be sure that SN 1999bw was of the same class as SN 2008S and NGC The Connection to Other Transients As we discuss in detail in 4, perhaps the primary distinguishing characteristic of this class of transients is their deeply embedded progenitors. Since we are unable to confirm the presence of such progenitors for the M85 transient or SN 1999bw, we are unable to make a direct analogy with SN 2008S and NGC 300, and instead rely on the fact that the tran- 11 The luminosity and blackbody scale have been adjusted to the distance employed here. sients themselves provides strong evidence for obscuration on few-hundred AU scales. In our effort to understand which cataloged transients might belong to the SN 2008S/NGC 300-like class we have examined archival imaging of many recent supernovae, as well as archetypal peculiar supernovae, including supernova impostors, LBV outbursts, and low-luminosity Type IIP supernovae. Here, we provide a brief discussion in an effort to orient the reader. As we mentioned in 1 (point 2), the low-luminosity Type IIn SN 2002bu is an interesting transient that may be a member of the class defined by SN 2008S and NGC 300. We checked archival Spitzer data of the host galaxy NGC 4242 (PID 69; PI G. Fazio) taken 2 years after discovery and we find a bright infrared point source detected in all IRAC bands within 0. 6 of the supernova position. Also, two epochs of MIPS data (PID 40204; PI R. Kennicutt), obtained in 2008, 6 years after explosion, and separated by just 6 days, reveal a 24µm source at the position of the supernova. This is qualitatively similar to the case of M85 and SN 1999bw, but because only a single post-explosion IRAC epoch exists, and because the two MIPS epochs are separated by such a short time, we are unable to definitively confirm that the MIR source is associated with SN 2002bu. SN 1997bs is an intriguing example of an object that does not fit into this class, although its peak absolute magnitude, lightcurve, color, and some spectral features are comparable to SN 2008S and NGC 300. In this case, a luminous unobscured progenitor has been identified in the optical (M V 7) and the transient itself has been argued to be the outburst of an LBV (Van Dyk et al. 2003). Of interest is the fact that no object has been subsequently identified in the optical at the site of the transient (Van Dyk 2005). We have checked archival Spitzer data of the host galaxy (NGC 3627) obtained in 2004 ( 7 years after discovery) by SINGS and we do not detect a bright MIR source at the site of SN 1997bs, in contrast with SN 2008S, M85, and SN 1999bw. The event SN 2003gm is also interesting in this context since it had photometric and spectroscopic evolution similar to SN 1997bs, and also showed an optically luminous progenitor (M V 7.5; Maund et al. 2006). The fact that both SN 1997bs and 2003gm had bright unobscured progenitors is our primary reason for excluding them from the class defined by SN 2008S and NGC 300. Historical LBV eruptions in nearby galaxies that have been initially classified as supernovae are also worth mentioning here. These include SN 1954J (e.g., Smith et al. 2000) and SN 2002kg (e.g., Maund et al. 2006; Van Dyk et al. 2006) in NGC 2403, SN 1961V (e.g., Humphreys 2005) in NGC 1058, and SN 2000ch in NGC 3432 (Wagner et al. 2004). As in the cases of SN 1997bs and SN 2003gm, a very important common difference between these objects and SN 2008S or NGC 300 is that they all had optically luminous progenitors with absolute magnitudes 6, consistent with originating from very massive stars. These transients also have other properties that are not consistent with SN 2008S-like explosions: (1) their peak absolute magnitudes range between 18 M V 9, (2) the transient timescales vary widely from a few days to years, and (3) they are not luminous MIR sources in archival Spitzer data. Other optical transients classified as sub-luminous Type IIP supernovae that might potentially fall into the class defined by SN 2008S and NGC 300 include SN 1994N, SN 1997D, SN 1999eu, 1999br, 2001dc, 2003Z, and 2005cs (e.g., Pa-

5 A New Class of Transients 5 storello et al. 2004, 2006). We note, however, that there are fundamental differences in the spectra of low-luminosity Type IIP SNe compared with SN 2008S-like transients. In particular, sub-luminous Type IIP SNe show Balmer lines with strong P-Cygni absorption profiles and velocities of a few thousand km s 1, as is observed in more luminous Type IIP SNe. This stands in sharp contrast with the Balmer lines in SN 2008S-like transients, which are fairly narrow (FWHM 1000 km s 1 ) and which do not show strong P-Cygni absorption features. In this way, SN 2008S-like transients most closely resemble the spectra of low-luminosity Type IIn SNe and LBV outbursts. In addition to the very interesting case of 2002bu, there were five other low-luminosity transients classified as impostors or unknown that could have been LBV outbursts, but for which no progenitor has been identified, and which might be 2008S-like: NGC 4656, SN 2001ac, 2006bv, 2006fp, and 2007sv (see 3). However, as we have emphasized and as the sources discussed above imply, the properties of the progenitor cannot be deduced from the character of the optical outburst alone (e.g., contrast SN 1997bs and SN 2003gm with SN 2008S). Thus, in order to understand the causal mapping between progenitor and explosion, a census of the progenitors must first be completed. What is clearly needed is a comprehensive survey for bright MIR sources in all nearby galaxies ( 10 Mpc) with (warm) Spitzer, analogous to the survey proposed by Kochanek et al. (2008) in the optical. In the next section we discuss our search for deeply-embedded progenitors in M33. We discuss a more complete census in RATES The absolute rate of transients analogous to NGC 300 and SN 2008S is uncertain. Current samples of supernovae over the last decade in the local volume within 10, 20, and 30 Mpc allow us to make only a rough estimate of the true rate. A systematic transient search in the local volume is crucial to solidify these numbers. Nevertheless, we estimate that 2008S-like transients occur with a frequency equivalent to 20% of the Type II supernova rate. We discuss the observed rates in 3.1 below, and then we enumerate several arguments suggesting that the sample of SN 2008S-like transients may be highly incomplete in the local universe ( 3.2) Observed Counts D 10 Mpc In addition to SN 2008S and NGC 300, 20 other corecollapse supernovae or supernova-like transients have been discovered within 10 Mpc since Sixteen were confirmed supernovae; they are SN 1999em, 1999ev, 1999gi, 1999gq, 2002hh, 2002ap, 2003gd, 2004am, 2004dj, 2004et, 2005af, 2005at, 2005cs, 2007gr, 2008bk, and 2008ax. Two were bona fide LBV eruptions (SN 2000ch, Wagner et al. 2004; 2002kg, Schwartz et al. 2003; Weis & Bomans 2005; Maund et al. 2006; Van Dyk et al. 2006). One was a Type IIn supernova potentially of the SN 2008S class (SN 2002bu in NGC 4242, D 8 Mpc; Puckett & Gauthier 2002; Hornoch 2002). Finally, the transient in NGC 4656 also had some of the spectral characteristics of low-luminosity IIn supernovae (e.g., narrow Hα in emission), but only reached an absolute magnitude of 11.5, (Rich et al. 2005; Elias-Rosa et al. 2005). 12 Throughout the discussion here we exclude Type Ia supernovae. Taken at face value, with no correction for incompleteness, these numbers imply that 2/22 9% or 3/22 14% (including SN 2002bu) of all optically bright transients are SN 2008S-like. Removing the two bona fide LBV outbursts (2000ch and 2002kg) for comparison with the supernova sample proper, if SN 2008S and NGC 300 are supernovae, they represent 10% and 15% (again, with SN 2002bu) of the sample D 20 Mpc A similar exercise can be carried out within the larger volume of 20 Mpc. With a peak absolute magnitude of 14, SN 2008S-like transients would have an apparent magnitude of 17.5 at D = 20 Mpc, without including a correction for extinction intrinsic to the transient. In fact, SN 2008S and NGC 300 had A V 1.2 (P08a) and A V (Bond et al. 2008), respectively. Using the Smartt et al. (2009) compilation, we find 29 IIP, 4 IIb, 15 Ib/c, 2 IIn (1998S & 2002bu), and 2 IIL supernovae in the last decade. There are 6 classified as LBV eruptions/impostors (1999bw, 2000ch, 2002kg, 2003gm, 2007sv, and NGC ), but only 2000ch and 2002kg have strong evidence for an LBV progenitor. Whether the yellow supergiant progenitor of 2003gm survived the explosion, and hence whether 2003gm was in fact a supernova, has not yet been definitively established (Maund et al. 2006). Similarly, although SN 2007sv, which reached an absolute magnitude of 14, bears some similarity to 1997bs, the nature of its progenitor has not been established (Duszanowicz et al. 2007; Harutyunyan et al. 2007). Finally, there are 3 events whose nature is classified in Smartt et al. (2009) as unknown: M85, NGC 300, and SN 2008S. With no correction for incompleteness, and taking only NGC 300 and SN 2008S, this compilation implies an overall rate of 2/61 3.3% within 20 Mpc. Including M85 and 2002bu in the sample of 2008S analogs doubles the rate. Including 1999bw, 2007sv, and NGC 4656 brings the overall observed rate of SN 2008S-like transients to 10% within D 20 Mpc D 30 Mpc With a limiting magnitude of 18 19, amateur and professional supernova surveys could in principle find SN 2008Slike transients with M V 14 to a luminosity distance of Mpc (again assuming no extinction). Objects of interest discovered over the last 10 years with faint absolute magnitudes in this distance range, and with the spectral characteristics of IIn supernovae analogous to SN 2008S include 2001ac, 2006bv, and 2006fp. A total of 92 core-collapse supernovae and 7 LBV eruptions appear in the recent compilation of Smartt et al. (2009). Three events M85, NGC 300, and SN 2008S are classified as unknown. Of the 7 LBV eruptions, only 2000ch, 2002kg, and 2003gm have LBV-like progenitors. The remainder 1999bw, 2001ac, 2006fp, and 2007sv have little or no progenitor information, spectra that resemble IIn s, and may be SN 2008S analogs. Given the number of confirmed NGC 300 and SN 2008S-like analogs (just 2) and those suspected of belonging to this class (M85, 1999bw, 2002bu), as well as those tradiationally labelled LBV outbursts, but with no strong confirmation (4656, 2001ac, 2006bv, 2006fp, 2007sv) 13 Note that NGC 4656 is not included in Smartt et al. (2009).

6 6 Thompson et al the overall rate ranges from 2/102 2% to 10/102 10% when measured within 30 Mpc Arguments for Incompleteness & Some Implications To summarize the previous subsections, conservatively taking NGC 300 and SN 2008S as the only examples of their type ever observed (that is, excluding all other low-luminosity transients), the observed rate is 9%, 3%, and 2% within 10, 20, and 30 Mpc volumes, respectively, with respect to all bright optical transients when averaged over the last 10 years. It is difficult to estimate the degree of uncertainty in these numbers since the surveys that find local supernovae are a combination of professional and amateur, with complicated and unquantitified selection functions for transient identification. Most surveys responsible for transient discoveries in the local universe do not have detailed descriptions of completeness in the literature. On the contrary, the large majority of the transients in the local volume are discovered by amateurs. 14 Nevertheless, it is possible to make an estimate of completeness that gives a sense of how large the correction to the rate of SN 2008S-like transients may be. As an example, B. Monard typically quotes ±0.2 magnitude (unfiltered) photometric errors on discovery observations of SNe detected at magnitude (e.g., Monard 2006). This photometric error translates to a signal-to-noise ratio of S/N 5. In order to estimate a lower limit on the incompleteness, we can compare this value for S/N with the mean detection efficiency for Type Ia supernovae in the SDSS-II Supernova Survey, which employs a well-tested photometric pipeline that uses difference imaging to subtract off the host galaxy (see Dilday et al. 2008). Their Figure 7 shows that for S/N 5, the detection efficiency is 0.5 in Sloan gri. The detection efficiency drops to for S/N 2, corresponding to a photometric error of ±0.4 magnitudes. Momentarily ignoring the difference in the shape of Type Ia lightcurves with respect to those of 2008S and NGC 300, the results of Dilday et al. (2008) imply that surveys achieving limiting magnitudes of are of order 10% complete for SN 2008S-like transients with M V 14 at Mpc (S/N 1 2). Thus, a factor of 10 correction should be applied to the 30 Mpc sample in from the Smartt et al. (2009) catalog for the true rate of low-luminosity IIn supernovae like SN 2008S. Within 20 Mpc the correction for incompletness is likely a factor of 5, and within 10 Mpc the incompleteness is probably closer to a factor of 2. Similar corrections should be applied to the observed rate of true LBV eruptions. There is another argument for incompleteness at the factor of 2 level within 10 Mpc for SN 2008S-like transients. As summarized by Horiuchi et al. (2008) (their Section II.B), the observed rate of supernovae of all types within 30 Mpc yields a ratio of Type Ia supernovae to core-collapse supernovae that is significantly in excess of the cosmic ratio measured at high redshift (0.5 z 1; Dahlen et al. 2008). Indeed, the ratio within 30 Mpc is large enough that we would expect to have seen several Type Ia supernovae within 10 Mpc in the last 10 years. Yet, none have been found. This fact implies that the ratio of Ia to core-collapse supernovae has been over-estimated within 30 Mpc because normal core-collapse supernovae are intrinsically fainter than Type Ia s. Thus, 14 Of the 13, 17 and 13 local (D 30 Mpc) supernovae (including Type Ia s) discovered in 2006, 2007, and 2008, only 1, 4 and 4 were discovered by professional searches. the sample of normal core-collapse supernovae is incomplete within 30 Mpc at the factor of level, even though these objects typically have peak absolute visual magnitudes of M V 16 to 18. Naively, analogs to SN 2008S with M V 14 would be 6 times more incomplete than normal supernovae at D 30 Mpc. Of course, the actual incompleteness correction depends on the overall extinction correction for the transient population and on the cadence of the observations since the lightcurve declines much faster for 2008S-like transients than for supernova of type IIP. This implied nearly order-of-magnitude incompleteness correction at D 30 Mpc for 2008S-like events strongly indicates that the sample is incomplete at order unity within D 10 Mpc. Additionally, a plot of the discovery rate of all supernovae (Ia s included) within 30 Mpc over the last 10 years shows an increasing trend, super-poisson variance, and strong dependence on the results and observing strategy of a single survey (LOSS; Li et al. 2000). There is also an asymmetry between the rate of discovery in the northern and the southern sky in excess of a simple extrapolation of star formation from the catalog of Karachentsev et al. (2004). Finally, there is an unquantified bias against small star-forming galaxies in the local universe. These points further solidify the case that the normal core-collapse supernova rate is incomplete, which implies that we are missing 2008S analogs in abundance in the local universe. Taking yet another angle on the question of overall rate, we may consider the a posteriori statistics of the events SN 2008S and NGC 300 themselves. Taking the overall supernova rate as 1 2 yr 1 within 10 Mpc implies a probability of of seeing two events in a single year if the overall rate of 2008S-like transients is 2% of the supernova rate. We consider this uncomfortably small. Similarly, if we were to take the true SN 2008S-like transient rate to be 100% of local supernova rate, we would be forced to explain the fact that only 2such events have been seen within 10 Mpc in the last 10 years. Given the discussion of incompleteness above, an overall rate of 20% with respect to supernovae gives a reasonable chance of seeing two in one year and of seeing only a handful on a 10 year baseline. Taken together, these arguments imply that the sample of transients in the local universe when averaged over the last 10 years is highly incomplete. We suggest that the incompleteness correction is a factor of 2, 5, and 10 for M V 14 transients at 10 ( 16.0 mag), 20 ( 17.5 mag), and 30 Mpc ( 18.4 mag), even before accounting for the potentially higher average extinction of these transients relative to normal supernovae. These estimates are consistent with the Richardson et al. (2002), who conclude that low-luminosity supernovae with M B 15 may constitute more than 20% of the overall supernova population (for related arguments, see Schaefer 1996; Hatano, Fisher, & Branch 1997; Pastorello et al. 2004). There are several immediate implications: 1. The true rate of SN 2008S-like transients is 20% of the core-collapse supernova rate. However, we emphasize that lower and higher values at the factor of 2 level are not excluded until a more thorough census has been made. 2. The true rate of massive star eruptions (LBV-like and otherwise) is similarly incomplete. The observed rate of LBV eruptions within 10 and 20 Mpc in and implies that they are 1 3 times more common than SN 2008S-like transients. Thus, the true rate of

7 A New Class of Transients 7 FIG. 1. Mid-infrared color-magnitude diagram for M33. Absolute magnitude at 4.5µm is plotted versus [3.6] [4.5] color for all detected sources (3σ limits denoted by dashed lines; see 4.1; Table 1). For comparison, the positions of the progenitors of NGC 300 (filled square) and SN 2008S (filled triangle; lower limit) are also plotted (P08ab). Stars analogous to these progenitors are exceedingly rare. Compare with Figure 2. The main sequence, AGB, and EAGB stars are clearly visible. Note the spur in the data extending to fainter M 4.5 and redder color at [3.6] [4.5] 1 and M , originating at the red extremum of the AGB population. To our knowledge, this is the first time that such a feature has been identified in a MIR CMD. massive LBV eruptions is likely 20 60% of the corecollapse supernova rate The observed rate of core-collapse supernovae is incomplete at the factor of 2 level for D 30 Mpc. 4. A FIRST CENSUS Because of the implied frequency of events similar to SN 2008S and NGC 300 ( 1, point 2) and the interesting character of their progenitors, we searched for analogous sources in archival Spitzer imaging of nearby galaxies. Our goal was to identify the underlying sub-population of massive stars from which these progenitors emerge, to characterize their properties and frequency, and to catalog them for future study. The key characteristics of the progenitors of SN 2008S and NGC 300 are that they are optically-obscured and deeply em- 15 We thank the anonymous referee for pointing this out. bedded, with very red MIR colors, that their bolometric luminosities are indicative of relatively low-mass massive stars (L and L, respectively), and that the several epochs on NGC 6946 revealed that the progenitor of SN 2008S was not highly variable in the 3 years before explosion (P08a). We discuss the variability of the progenitor of NGC 300 in 4.4 below based on just two pre-explosion epochs. Although we are only able to derive a lower limit to its RMS variation at 4.5 µm, like the progenitor of SN 2008S, we find that it too is consistent with a low level of variability. For a first census, we searched for the deepest archival Spitzer observations of a nearby relatively massive bright starforming galaxy, with already extant optical catalogs. The Triangulum galaxy M33 is a perfect test case. It has an absolute B-band magnitude of M B 19.2, a distance of 0.96 Mpc (distance modulus µ = 24.92; Bonanos et al. 2006), and it has extensive optical (e.g., Hartman et al. 2006; Massey et al. 2006, hereafter M06), Hα (Massey et al. 2007, hereafter M07), and MIR and FIR imaging (McQuinn et al. 2007, here-

8 8 Thompson et al FIG. 2. Same as Figure 1, but focused on the red and bright region of interest. Here, the [3.6] [4.5] > 1.5 and M 4.5 < 10 selection for extremely red and bright objects is shown explicitly by the dotted lines, as are the 18 EAGB sources in M33 that meet these criteria (large open circles; see 4.2). We denote them S1-S18, ordered by color. The reddest source, S1, is shown in Figure 3. Most of these sources are highly variable (see Figs. 7, A1, & A2). SEDs are shown in Figure 5 (see Table 2). The EAGB population is clearly distinct from the optically-luminous LBV candidates from the catalog of M07, shown by the open triangles (see 4; Table 3). Most of the LBV candidates are considerably more bolometrically luminous and much less variable at 4.5µm than the EAGB sample (see Figs. 5, 7, B1, & B2). A number of cool hyper-giants such as VY CMa, NML Cyg, and M33 Var A, as well as η-carina, are shown for comparison. The brightest LBV candidate is M33 Var C. The LBV candidates appear to be bimodal in MIR color (right panel of Fig. 5). after Mc07). A similarly rich dataset exists for several other local galaxies, including the Magellanic clouds (e.g., Blum et al. 2006; Bolatto et al. 2007) and M31 (e.g., Mould et al. 2008). An analysis similar to that described below will be the subject of future work (but, see 4.5). In this Section, we present the MIR color-magnitude and color-color diagrams for all cataloged point sources in M33 obtained from multi-epoch archival Spitzer/IRAC (Fazio et al. 2004) observations (PI R. Gehrz; PID 5). This dataset allows us to search for and identify stars analogous to the progenitors of SN 2008S and NGC 300. We also present a variability study for both the reddest sources we find in M33 (all extreme-agb EAGB stars) and for the optically-selected LBV candidates from M07, which are detected in the MIR imaging. We find that the population of the reddest sources those most likely to be true analogs of the progenitors of SN 2008S and NGC 300 is completely distinct from the population of LBV candidates in the primary metrics of color, magnitude, and variability. Indeed, we find very few sources with the properties of the SN 2008S and NGC 300 progenitors. Section 4.1 describes our procedure for extracting point sources and the resulting catalog. In 4.2 we present the color-magnitude diagram for M33 and we discuss the reddest sources vis á vis the optically-selected LBV candidates from M07. In 4.3 and 4.4 we discuss their SEDs and variability, respectively. Finally, in 4.5 we discuss a preliminary search for similar sources in NGC 300, the LMC, and the SMC Catalog We coadded six epochs of mid-ir imaging of M33 obtained between Jan. 9, 2004 and Feb. 4, 2006 with IRAC ( µm; see Mc07 for details of the observing program). We produced the coadds in all four IRAC channels from the flux cal-

9 A New Class of Transients V FIG. 3. Image showing the point source S1 at 3.6 µm, 4.5 µm, as well as V and I band (M06). This is the reddest of the 18 sources in [3.6] [4.5] color selected in Figure 1. It is optically obscured (see Fig. 5). S1 is also highly variable in color and 4.5 µm magnitude (see Figs. 6 & 7). ibrated mosaics provided by the Spitzer Science Center (post- BCD data). Our final mosaics cover an area of ( pixels, with 1. 2/pix) centered on M33, this is approximately within R 25 ( 35 ; de Vaucouleurs et al. 1991). We performed source detection and PSF fitting photometry on the coadds using the DAOPHOT/ALLSTAR package (Stetson 1992). The PSF magnitudes obtained with ALLSTAR were transformed to Vega-calibrated magnitudes using simple zero point shifts derived from aperture photometry (using a 12 radius), performed in the original images, of bright and isolated stars in each band. We estimate errors in the photometric transformations to the Vega system of 0.04 mag in 3.6 µm, 0.05 mag in 4.5 µm, 0.07 mag in 5.8 µm, and 0.07 mag in 8.0µm channel. For the 3.6 and 4.5µm bands, the detection limits (3σ) in the coadds are [3.6] 18.9 mag and [4.5] = 18.2 mag, respectively. The sample becomes incomplete at [3.6] [4.5] > 17.1, approximately 0.5 magnitudes deeper than Mc07. A total of 80,000 sources are detected in either 3.6 or 4.5 µm. We cross-matched the two catalogs using a 0.5 pixel (0. 6) matching radius to obtain a final catalog with 53,200 individual sources detected at both 3.6 and 4.5µm (see Table 1). There are several reasons for producing new point source and variability (see 4.4) catalogs, given the already existing catalog from Mc07: (1) because we are specifically looking for objects with extreme colors, we wanted to be able to relax the criterion for point source detection in all IRAC bands (in the Figures that follow, all sources of interest are detected with 3σ confidence); (2) we wanted to be able to derive our own upper limits in each band for the same reason; (3) we wanted full IRAC SEDs, whereas the catalog of Mc07 does not provide data at 5.8µm; (4) we wanted to combine the images used by Mc07 with a sixth archival epoch (PI R. Gehrz; PID 5); (5) for the reddest stars, we wanted to derive full sixepoch lightcurves for a more complete measure of variability. The resulting MIR color-magnitude diagram looks different from that in Mc07. Of the 18 sources we discuss extensively below, only 2 were detected at both 3.6 and 4.5 µm in the catalog of Mc07. I 4.2. The Color-Magnitude Diagram The primary result of this effort on M33 was the production of the MIR color-magnitude diagram (CMD), shown in Figure 1, which shows the [3.6] [4.5] color for all the sources detected at both 3.6 and 4.5µm as a function of absolute magnitude at 4.5µm, M 4.5. The dashed line marks the 3σ completeness limit in this plane. The Spitzer colors and magnitudes of the progenitors of SN 2008S and NGC 300 are shown for comparison (filled triangle and filled square, respectively; P08a, P08b). Note that there are remarkably few objects inhabiting the bright and (very) red region of the CMD. Among the massive stars in M33 (see 5), 18, 186, and 567 have both M 4.5 < 10 and [3.6] [4.5] color larger than 1.5, 1.0, and 0.7, respectively, which correspond to blackbody temperatures of 500, 700, and 1000 K, respectively. A total of 2264 point sources are detected with [3.6] [4.5] 0.7 and M 4.5 < 10. Figure 2 shows an expanded view of the brightest MIR sources. In addition to the progenitors of SN 2008S and NGC 300, we include several well-studied LBVs (e.g., η-carina), and a number of cool hyper-giants (VY CMa, NML Cyg, and Var A in M33) for comparison. The MIR magnitudes and colors for η-carina, VY CMa, and NML Cyg were synthesized from ISO spectra (Sloan et al. 2003). The magnitudes of M33 Var A were obtained from Humphreys et al. (2006). Also included are the 16 sources matched between the catalog of MIR sources presented here and the LBV sample of M07, obtained from narrow-band Hα imaging, using a 0. 6 matching radius. The larger circles within the dotted lines show the EAGB stars, which we discuss in detail below. The primary point of Figures 1 and 2 is to show that there are very few massive stars with the colors and MIR luminosities of the progenitors of SN 2008S and NGC 300. Although it is difficult to identify a quantitative criterion for inclusion in the class defined by the progenitors of SN 2008S and NGC 300, it is clear from Figure 1 that the number of analogs in color and magnitude is very small with respect to the total number of massive stars in M33. For example, if to be included as an analog to the SN 2008S and NGC 300 progenitors we require that M 4.5 be brighter than or equal to the progenitor of NGC 300 and we require that the color be redder than the lower limit on the progenitor of SN 2008S, we find a single source. If we require M 4.5 < 10 and color redder than SN 2008S, we find just two sources. Casting the net more widely, for the purpose of having a sample larger than one or two objects and in an effort to be conservative, we use M 4.5 < 10 and [3.6] [4.5]> 1.5 to identify a sample of 18 EAGB stars (large open circles). We discuss the spectra and variability properties of these sources in 4.3 and 4.4. Our choice of the cuts M 4.5 < 10 and [3.6] [4.5] > 1.5 to identify objects of interest is somewhat arbitrary. Because our argument in this paper relies on the fact that analogs to the progenitors of SN 2008S and NGC 300 are intrinsically rare ( 1 & 5), this issue deserves discussion. The magnitude limit is straightforward: it is meant to select objects that have bolometric luminosities indicative of massive stars ( 8 M ). In 4.3 below, we show that M 4.5 < 10 is conservative; only half of the 18 sources selected have bolometric luminosities large enough to be massive stars. Note that had we required M 4.5 to be brighter than or equal to the NGC 300 progenitor, we would exclude 8 of our 18 sources (see Fig. 5). Our goal was to not miss any deeply-embedded massive stars and the criterion M 4.5 < 10 accomplishes that goal.