Abstract
It is well known that massive stars (M > 8 M☉) evolve up to the collapse of the stellar core, resulting in most cases in a supernova (SN) explosion. Their heterogeneity is related mainly to different configurations of the progenitor star at the moment of the explosion and to their immediate environments. We present photometry and spectroscopy of SN 2010bt, which was classified as a Type IIn SN from a spectrum obtained soon after discovery and was observed extensively for about 2 months. After the seasonal interruption owing to its proximity to the Sun, the SN was below the detection threshold, indicative of a rapid luminosity decline. We can identify the likely progenitor with a very luminous star (log L/L☉ ≈ 7) through comparison of Hubble Space Telescope images of the host galaxy prior to explosion with those of the SN obtained after maximum light. Such a luminosity is not expected for a quiescent star, but rather for a massive star in an active phase. This progenitor candidate was later confirmed via images taken in 2015 (∼5 yr post-discovery), in which no bright point source was detected at the SN position. Given these results and the SN behavior, we conclude that SN 2010bt was likely a Type IIn SN and that its progenitor was a massive star that experienced an outburst shortly before the final explosion, leading to a dense H-rich circumstellar environment around the SN progenitor.
Export citation and abstract BibTeX RIS
1. Introduction
Type II "narrow" (IIn) supernovae (SNe) are a heterogeneous subset of core-collapse SNe (CC-SNe; Schlegel 1990). According to Li et al. (2011) and Smith et al. (2011a), they represent 9% of all CC-SNe and seem to preferentially occur in small and late-type spiral galaxies.
This subclass of objects can be distinguished from other types of CC-SNe by their spectral appearance. On the other hand, they can also be confused with nonterminal outbursts of very massive stars. The broad absorption lines typical of many SNe are weak or absent in SNe IIn throughout their evolution. Instead, they show strong, narrow Balmer emission components (with FWHM intensity ranging from a few tens to a few hundreds of km s−1) atop broader emission (which can have intermediate-velocity components with FWHM ≈ 1000 km s−1, as well as broad components with FWHM of a few thousands of km s−1; see, e.g., Filippenko 1997). The narrow lines are thought to arise from the surrounding circumstellar material (CSM) ionized by the shock-interaction emission (e.g., Chugai & Danziger 1994). Their light curves, however, exhibit a wide range of properties (see, e.g., Smith 2017). This diversity is related to the mass-loss history during the evolution of massive stars.
It was suggested that the progenitors of a fraction of these interacting SNe are massive stars in a luminous blue variable (LBV) phase. These stars are among the most luminous (Mbol < −9.6 mag) and massive (>50 M☉) stars in late-type galaxies (e.g., Humphreys & Davidson 1994). The evidence arises from the identification of the progenitor stars of a few SNe IIn, such as SN 2005gl (Gal-Yam et al. 2007; Gal-Yam & Leonard 2009), SN 2009ip (Smith et al. 2010, 2014; Foley et al. 2011), and SN 2015bh (Elias-Rosa et al. 2016; Thöne et al. 2017). A luminous, blue point-like source, originating from either a young cluster or a single star, was first proposed as the progenitor of SN 2010jl by Smith et al. (2011b). Subsequently, Fox et al. (2017) argued that this source was most likely a massive young cluster, although they did not rule out the possibility of a very luminous progenitor star obscured by dust (see also Dwek et al. 2017). However, the nature of these progenitors is not fully clear, and other precursor channels have been proposed (see, e.g., Mauerhan & Smith (2012) in the case of the SN 1998S, or Smith (2014) for a larger review of this topic).
As mentioned above, the powerful nonterminal outbursts of very massive stars can mimic genuine SN explosions, in terms of energetics and spectral appearance (both show incipient narrow lines of hydrogen in emission); thus, they are usually referred to as "SN impostors" (e.g., SN 1997bs, Van Dyk et al. 2000; SN 2000ch, Pastorello et al. 2010; see also the general discussions in Smith et al. 2011c; Van Dyk & Matheson 2012). Discriminating between SN impostors and SNe IIn is challenging, and for several objects both possibilities remain viable (see, e.g., the case of SN 2007sv; Tartaglia et al. 2015).
SN 2010bt was discovered in NGC 7130 on 2010 April 17.11 (UT dates are used throughout this paper) at an unfiltered magnitude of 15.9 and confirmed on 2010 April 18.14 at a magnitude of 15.8 (Monard 2010; Figure 1). The presence of Balmer emission with multicomponent profiles in the spectrum from 2010 April 18.39 led to the classification of SN 2010bt as an SN IIn, a few days after discovery (Turatto et al. 2010).
NGC 7130 (or IC 5135; z = 0.01617 ; morphological type Sa pec) was classified as a Seyfert 2 galaxy by Phillips et al. (1983)—i.e., an active galactic nucleus (AGN) with obscuring material that precludes a direct view of its nuclear region. Throughout the paper, we adopt a distance to NGC 7130 of 65.4 ± 4.6 Mpc (μ = 34.08 ± 0.15 mag), resulting from the recession velocity of the galaxy (de Vaucouleurs et al. 1976) corrected for Local Group infall into the Virgo Cluster (Mould et al. 2000), vVir = 4771 ± 17 km s−1 (z = 0.016), assuming H0 = 73 km s−1 Mpc−1 (values taken from NED).
This paper presents and discusses our photometric and spectroscopic observations of SN 2010bt in Section 2. In Section 3, we compare pixel by pixel four sets of Hubble Space Telescope (HST) data, and we discuss the nature of the progenitor star in Section 4. We conclude in Section 6.
2. The Nature of SN 2010bt
2.1. Photometry
Optical BVRI images of SN 2010bt were obtained with the 1.3 m Small and Moderate Aperture Research Telescope System (SMARTS)+ANDICAM at Cerro Tololo Inter-American Observatory (CTIO), the 6.5 m Magellan Clay Telescope+LDSS-3 at Las Campanas Observatory, and the 3.58 m New Technology Telescope (NTT)+EFOSC2 at the European Southern Observatory (ESO) of La Silla Observatory, all located in Chile. The data were obtained thanks to a collaboration between American and European researchers. We also include in our data set unfiltered data from the Bronberg Observatory (South Africa) and archival images from the 1.0 m Jacobus Kapteyn Telescope at Roque de Los Muchachos Observatory (Spain). More information about the telescopes and instruments used in this follow-up campaign can be found in Table 1.
Table 1. Basic Information about Telescopes and Instruments Used
Table Keya | Telescope | Instrument | Pixel Scale | Location |
---|---|---|---|---|
(arcsec pixel−1) | ||||
ANDICAM | 1.3 m SMARTSb | ANDICAM | 0.37 | CTIO,c Chile |
BO | 0.3 m Meade RCX400 f/8 telescope | SBIG ST8-XME CCD | 1.37 | Bronberg Observatory, South Africa |
EFOSC2 | 3.6 m New Technology Telescope | EFOSC2 | 0.12 | ESO,d La Silla Obs., Chile |
HAWK-I | 8.2 m Very Large Telescope UT4 | HAWK-I | 0.11 | ESO, La Silla Obs., Chile |
HST_WFPC2 | 2.4 m HST | WFPC2 | 0.05e | ⋯ |
HST_ACS/HRC | 2.4 m HST | ACS/HRC | 0.03 | ⋯ |
HST_ACS/WFC | 2.4 m HST | ACS/WFC | 0.05 | ⋯ |
HST_WFC3 | 2.4 m HST | WFC3/UVIS | 0.04 | ⋯ |
JAG | 1.0 m Jacobus Kapteyn Telescope | JAGf | 0.33 | ORM,g La Palma, Spain |
LDSS-3 | 6.5 m Magellan Clay Telescope | LDSS-3 | 0.19 | Las Campanas Obs., Chile |
SOFI | 3.6 m New Technology Telescope | SOFI | 0.29 | ESO, La Silla Obs., Chile |
Spitzer | 0.8 m Spitzer Space Telescope | IRAC | 0.60 | ⋯ |
Swift | 0.3 m Ritchey–Chretien UV/optical Telesc. | Swift | 0.50 | ⋯ |
Notes.
aSee Tables 3–5. bSmall and Moderate Aperture Research Telescope System. cCerro Tololo Inter-American Observatory. dEuropean Southern Observatory. eWFPC2 contains four chips. The SN 2010bt field was observed with the Planetary Camera (005 pixel−1). fJKT Acquisition and Guidance Unit. gObservatorio del Roque de los Muchachos.Download table as: ASCIITypeset image
The photometric observations were processed following the standard recipe in iraf for CCD images (trimming, overscan, bias, and flat-field corrections). Due to the location of SN 2010bt in NGC 7130, contamination of the SN photometry by the host-galaxy light was a serious problem. We therefore used the template-subtraction technique to remove this background and hence to measure more accurately the SN magnitudes. The template images of NGC 7130 were obtained with the 1.3 m SMARTS telescope+ANDICAM at CTIO on 2011 August 03. Each SN image was registered geometrically and photometrically with its corresponding template using a dedicated pipeline (SNOoPY). This consists of a collection of python scripts calling standard iraf tasks (through pyraf) and other specific analysis tools, in particular sextractor, for source extraction and star/galaxy separation. The instrumental magnitudes of the SN and the reference stars in the SN field were measured in the subtracted images (produced with hotpants) using the point-spread function (PSF) fitting technique with the daophot package. In order to calibrate the instrumental magnitudes to the standard photometric vegamag system, we used the magnitudes and colors of 15 local sequence stars in the SN field (Figure 1 and Table 2). The unfiltered data were transformed to the Johnson–Cousins R band, for which the effective wavelength is similar to the natural instrumental band defined by the CCD quantum efficiency of the detector that was used.
Table 2. Magnitudes of the Local Comparison Starsa
Star | B | V | R | I | J | H | K |
---|---|---|---|---|---|---|---|
(mag) | (mag) | (mag) | (mag) | (mag) | (mag) | (mag) | |
1 | 19.86 (0.07) | 19.04 (0.05) | 18.51 (0.03) | 17.96 (0.04) | ⋯ | ⋯ | ⋯ |
2 | 20.24 (0.07) | 18.70 (0.05) | 17.60 (0.07) | 16.19 (0.04) | 14.81 (0.03) | 14.14 (0.02) | 13.85 (0.04) |
3 | 21.28 (0.07) | 19.66 (0.04) | 18.30 (0.04) | 16.51 (0.05) | 14.84 (0.03) | 14.20 (0.04) | 13.92 (0.05) |
4 | 18.68 (0.04) | 17.80 (0.04) | 17.20 (0.03) | 16.62 (0.03) | 15.91 (0.07) | 15.40 (0.09) | 15.38 (0.13) |
5 | 18.13 (0.04) | 17.40 (0.03) | 16.92 (0.02) | 16.51 (0.03) | 15.90 (0.07) | 15.54 (0.08) | 15.43 (0.14) |
6 | 19.23 (0.04) | 18.28 (0.04) | 17.65 (0.04) | 17.10 (0.03) | 16.56 (0.11) | 15.96 (0.11) | 15.52 (0.17) |
7 | 17.31 (0.03) | 16.49 (0.03) | 16.00 (0.03) | 15.54 (0.03) | 14.86 (0.03) | 14.27 (0.04) | 14.29 (0.05) |
8 | 19.61 (0.05) | 19.15 (0.06) | 18.89 (0.04) | 18.51 (0.04) | ⋯ | ⋯ | ⋯ |
9 | 20.63 (0.04) | 19.19 (0.04) | 18.35 (0.04) | 17.18 (0.04) | 15.94 (0.06) | 15.32 (0.07) | 15.33 (0.13) |
10 | 19.97 (0.06) | 18.68 (0.03) | 17.87 (0.03) | 17.05 (0.02) | 15.96 (0.08) | 15.40 (0.07) | 15.29 (0.13) |
12 | 17.08 (0.03) | 16.46 (0.02) | 16.12 (0.04) | 15.73 (0.02) | 15.26 (0.04) | 14.84 (0.05) | 14.76 (0.09) |
12 | 20.12 (0.04) | 19.18 (0.05) | 18.62 (0.05) | 18.12 (0.07) | ⋯ | ⋯ | ⋯ |
13 | 20.97 (0.03) | 19.58 (0.01) | 18.76 (0.02) | 17.79 (0.04) | 16.65 (0.11) | 16.20 (0.15) | 15.38 (0.15) |
14 | 18.99 (0.04) | 18.47 (0.04) | 18.11 (0.04) | 17.80 (0.03) | ⋯ | ⋯ | ⋯ |
15 | 15.99 (0.04) | 14.98 (0.03) | 14.38 (0.06) | 13.83 (0.05) | 12.93 (0.03) | 12.41 (0.02) | 12.33 (0.03) |
Note.
aQuoted uncertainties are 1σ.Download table as: ASCIITypeset image
Near-infrared (NIR) observations were also obtained with the 3.6 m NTT+SOFI at the ESO Observatory of La Silla. We include in this work the data obtained at the 8.2 m unit telescope UT4 of the Very Large Telescope+HAWK-I at the ESO Observatory of Cerro Paranal (PI: E. Cappellaro, 083.D-0259(A)), presented by Miluzio et al. (2013). The data were reduced using standard procedures for the VLT+HAWK-I data (see Miluzio et al. 2013, for details on the data reduction). For the NIR images we used two sets of images of the field as templates to subtract from the SN images, depending on the instrumentation, e.g., images from ESO NTT+SOFI taken on 2004 October 04 (PI: P. Lira, 074.B-0375(A)) and images from ESO VLT UT4+HAWK-I on 2011 April 24 (PI: E. Cappellaro, 083.D-0259(A)). The instrumental JHKs photometry was calibrated using Two Micron All Sky Survey (2MASS18 ) stars in the field.
When the transient was not detected, upper limits were adopted corresponding to 2.5 times the standard deviation in the background. Uncertainty estimates were obtained through an artificial star experiment, combined (in quadrature) with the PSF fit error returned by daophot and the propagated uncertainties from the photometric calibration.
The final vegamag calibrated magnitudes of SN 2010bt are listed in Tables 3 and 4. The single epoch with magnitudes in the Sloan system was transformed to the vegamag system and included in Table 3 by employing the relations given by Blanton & Roweis (2007).
Table 3. Optical Johnson–Cousins Photometry of SN 2010bt (vegamag)a
Date | MJD | Phaseb | U | B | V | R | I | Instrument Key |
---|---|---|---|---|---|---|---|---|
(days) | (mag) | (mag) | (mag) | (mag) | (mag) | |||
20010813 | 52,134.0 | −3169.1 | ⋯ | ⋯ | ⋯ | >18.0 | ⋯ | JAG |
20100417 | 55,303.1 | 0.0 | ⋯ | ⋯ | ⋯ | 16.00 (0.41) | ⋯ | BO |
20100418 | 55,304.1 | 1.0 | ⋯ | ⋯ | ⋯ | 16.08 (0.18) | ⋯ | BO |
20100418 | 55,304.4 | 1.3 | ⋯ | 17.00 (0.08) | 16.33 (0.08) | 16.05 (0.10) | 15.79 (0.12) | EFOSC2 |
20100420 | 55,306.4 | 3.3 | ⋯ | 17.11 (0.03) | 16.54 (0.03) | 16.26 (0.04) | 15.88 (0.04) | ANDICAM |
20100420 | 55,306.9 | 3.8 | 16.88 (0.14) | 17.33 (0.15) | 16.58 (0.17) | ⋯ | ⋯ | Swift |
20100422 | 55,308.1 | 5.0 | 17.09 (0.17) | 17.35 (0.16) | 16.42 (0.15) | ⋯ | ⋯ | Swift |
20100423 | 55,309.4 | 6.2 | ⋯ | 17.42 (0.08) | 16.81 (0.07) | 16.50 (0.05) | 16.17 (0.04) | ANDICAM |
20100424 | 55,310.2 | 7.0 | 17.26 (0.19) | 17.58 (0.18) | 17.09 (0.24) | ⋯ | ⋯ | Swift |
20100426 | 55,312.8 | 9.7 | 17.71 (0.26) | 17.75 (0.22) | ⋯ | ⋯ | ⋯ | Swift |
20100427 | 55,313.3 | 10.2 | ⋯ | 17.89 (0.18) | 17.16 (0.16) | 16.73 (0.09) | 16.39 (0.14) | ANDICAM |
20100428 | 55,314.7 | 11.6 | 17.71 (0.26) | 17.77 (0.20) | 17.25 (0.26) | ⋯ | ⋯ | Swift |
20100430 | 55,316.0 | 12.9 | >18.3 | ⋯ | ⋯ | ⋯ | ⋯ | Swift |
20100430 | 55,316.0 | 12.9 | ⋯ | 18.23 (0.28) | 17.46 (0.30) | ⋯ | ⋯ | Swift |
20100430 | 55,316.4 | 13.3 | ⋯ | 18.11 (0.15) | 17.36 (0.13) | 16.89 (0.10) | 16.66 (0.09) | ANDICAM |
20100501 | 55,317.9 | 14.7 | ⋯ | ⋯ | >17.8 | ⋯ | ⋯ | Swift |
20100502 | 55,318.9 | 15.7 | 18.15 (0.35) | 18.16 (0.27) | ⋯ | ⋯ | ⋯ | Swift |
20100503 | 55,319.4 | 16.3 | ⋯ | 18.24 (0.14) | 17.52 (0.12) | 16.95 (0.09) | 16.70 (0.10) | ANDICAM |
20100504 | 55,320.9 | 17.8 | 18.27 (0.39) | 18.46 (0.34) | 17.69 (0.38) | ⋯ | ⋯ | Swift |
20100506 | 55,322.1 | 19.0 | ⋯ | ⋯ | ⋯ | 17.18 (0.14) | ⋯ | BO |
20100506 | 55,322.2 | 19.1 | 18.10 (0.35) | 18.27 (0.30) | ⋯ | ⋯ | ⋯ | Swift |
20100506 | 55,322.2 | 19.1 | ⋯ | ⋯ | >17.7 | ⋯ | ⋯ | Swift |
20100507 | 55,323.4 | 20.2 | ⋯ | 18.56 (0.06) | 17.79 (0.05) | 17.30 (0.04) | 16.96 (0.06) | ANDICAM |
20100508 | 55,324.5 | 21.4 | >18.3 | >18.6 | ⋯ | ⋯ | ⋯ | Swift |
20100510 | 55,326.3 | 23.2 | ⋯ | 18.98 (0.13) | 17.98 (0.11) | 17.52 (0.06) | 17.22 (0.06) | ANDICAM |
20100511 | 55,327.1 | 24.0 | ⋯ | ⋯ | ⋯ | >16.7 | ⋯ | BO |
20100518 | 55,334.4 | 31.3 | ⋯ | 19.58 (0.11) | 18.59 (0.09) | 17.98 (0.06) | 17.69 (0.07) | ANDICAM |
20100519 | 55,335.4 | 32.3 | ⋯ | 19.60 (0.08) | 18.62 (0.07) | 18.01 (0.06) | 17.68 (0.06) | ANDICAM |
20100523 | 55,339.3 | 36.2 | ⋯ | 19.84 (0.14) | 18.84 (0.12) | 18.19 (0.08) | 17.91 (0.12) | ANDICAM |
20100524 | 55,340.1 | 37.0 | ⋯ | ⋯ | ⋯ | >16.5 | ⋯ | BO |
20100524 | 55,340.3 | 37.2 | ⋯ | 19.89 (0.14) | 18.74 (0.12) | 18.23 (0.11) | 17.95 (0.12) | ANDICAM |
20100526 | 55,342.3 | 39.2 | ⋯ | 20.10 (0.12) | 19.04 (0.11) | 18.34 (0.07) | 18.03 (0.08) | ANDICAM |
20100531 | 55,347.4 | 44.2 | ⋯ | 20.47 (0.13) | 19.22 (0.12) | 18.48 (0.07) | 18.16 (0.09) | ANDICAM |
20100604 | 55,351.3 | 48.2 | ⋯ | 20.79 (0.13) | 19.48 (0.11) | 18.63 (0.10) | 18.24 (0.07) | ANDICAM |
20100613 | 55,360.1 | 57.0 | ⋯ | ⋯ | ⋯ | >16.6 | ⋯ | BO |
20100915 | 55,454.2 | 151.1 | ⋯ | >19.8 | >19.8 | >18.9 | >19.7 | EFOSC2 |
20101004 | 55,473.1 | 170.0 | ⋯ | >20.4 | >20.0 | >19.3 | >18.7 | EFOSC2 |
20101009 | 55,478.4 | 175.3 | ⋯ | ⋯ | 22.50 (0.05) | ⋯ | ⋯ | HST_ACS/WFC |
20101028 | 55,497.2 | 194.1 | ⋯ | >21.3 | >21.0 | >20.2 | >19.7 | EFOSC2 |
20150917 | 57,283.0 | 1979.9 | ⋯ | >18.6 | ⋯ | >18.5 | ⋯ | LDSS-3 |
20151014 | 57,309.0 | 2005.9 | ⋯ | ⋯ | >24.4 | ⋯ | ⋯ | HST_WFC3 |
Notes.
aQuoted uncertainties are 1σ. bPhases are relative to the discovery, MJD = 55,303.1.Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
Table 4. Near-infrared Photometry of SN 2010bt (vegamag)a
Date | MJD | Phaseb | J | H | K | Instrument Key |
---|---|---|---|---|---|---|
(days) | (mag) | (mag) | (mag) | |||
20100525 | 55,341.0 | 37.9 | ⋯ | ⋯ | 16.69 (0.23) | HAWK-I |
20100606 | 55,353.0 | 49.9 | ⋯ | ⋯ | 17.51 (0.31) | HAWK-I |
20100707 | 55,384.0 | 80.9 | ⋯ | ⋯ | 18.49 (0.22) | HAWK-I |
20100917 | 55,456.1 | 153.0 | >18.7 | >17.4 | >18.9 | SOFI |
20101029 | 55,498.1 | 195.0 | >17.2 | >15.4 | ⋯ | SOFI |
Notes.
aQuoted uncertainties are 1σ. bPhases are relative to the discovery, MJD = 55,303.1.Download table as: ASCIITypeset image
Space-based optical and ultraviolet (UV) data taken with the Ultraviolet/Optical Telescope (UVOT) on board Swift complement our ground-based photometry. The calibrated SN images in the UVOT-Vega system were obtained from SOUSA (Swift's Optical/Ultraviolet Supernova Archive; Brown et al. 2014). Upper limits corresponding to three times the standard deviation in the background were estimated when the transient was not detected. The Swift photometry for SN 2010bt is reported in Table 5.
Table 5. Swift UV Photometry of SN 2010bt (vegamag)a
Date | MJD | Phaseb | UVW1 | UVM2 | UVW2 |
---|---|---|---|---|---|
(days) | (mag) | (mag) | (mag) | ||
20100420 | 55,306.9 | 3.8 | >18.5 | >18.8 | >18.9 |
20100422 | 55,308.1 | 5.0 | 18.20 (0.34) | ⋯ | ⋯ |
20100422 | 55,308.2 | 5.0 | ⋯ | >18.8 | >18.9 |
20100424 | 55,310.2 | 7.0 | 18.48 (0.42) | ⋯ | ⋯ |
20100424 | 55,310.2 | 7.0 | ⋯ | >18.7 | >18.9 |
20100426 | 55,312.8 | 9.7 | >18.5 | ⋯ | >18.4 |
20100428 | 55,314.7 | 11.6 | >18.5 | >18.3 | >18.9 |
20100430 | 55,316.0 | 12.9 | >18.6 | >18.9 | >19.0 |
20100502 | 55,318.9 | 15.7 | >18.6 | >18.9 | >19.0 |
20100504 | 55,320.9 | 17.8 | >18.5 | >18.8 | >18.9 |
20100506 | 55,322.2 | 19.1 | >18.5 | >18.9 | >18.9 |
20100508 | 55,324.5 | 21.4 | >18.5 | >18.9 | >18.9 |
Notes.
aQuoted uncertainties are 1σ. bPhases are relative to the discovery, MJD = 55,303.1.Download table as: ASCIITypeset image
Finally, HST observed the SN 2010bt field at four epochs: 1994, 2003, 2010, and 2015 (see Sections 3 and 4 for more details). The vegamag magnitudes of the progenitor candidate and transient were obtained using Dolphot (see Table 10).
SN 2010bt was observed during the first ∼50 (optical) to 80 (NIR) days after discovery. Subsequently, it came too close to the Sun's direction and was lost. The telescopes could point to the field again around 2 months after the last observation in the NIR; however, at that time the SN was no longer visible (at least with ground-based telescopes). Instead, it was detected in a F606W (∼V-) band HST exposure obtained ∼175 days after discovery. The UBVRIJHK light curves, including the UBV-UVOT data, are shown in Figure 2. In the figure, phase is relative to the discovery date on 2010 April 17.10, or MJD 55,303.1. However, in the following, and in order to better compare SN 2010bt with other SNe, we will refer to the V maximum computed as the discovery date minus 10 days (see Section 2.3 for more details).
Download figure:
Standard image High-resolution imageAs one can see in Figure 2, SN 2010bt was discovered after maximum light. The light curves show a decline similar to those of the rapidly declining SNe II studied by Anderson et al. (2014; with an average V-band light-curve decline >1.3 mag in the first 50 days after peak). In fact, the UBV light curves of SN 2010bt show constant decline rates of 4.3, 4.0, and 3.3 mag/50 days,19 respectively, while the redder light curves (RI) slightly flatten out in brightness around day 30 after SN discovery. The K-band light curve was also slowly declining (2.0 mag/50 days, considering the interval from ∼38 to 80 days after the discovery date). SN 2010bt was only detected in two epochs in the Swift UVW1 band, with a similar decline to that of the bluer light curves (∼3.4 mag/50 days).
Figure 3 shows the evolution of the absolute V magnitude of SN 2010bt, compared with the SNe IIn SN 1996al (Benetti et al. 2016) and SN 1998S (Fassia et al. 2000; Liu et al. 2000; Pozzo et al. 2004) and the SN IIn/SN impostors SN 2009ip (Fraser et al. 2013; Mauerhan et al. 2013; Pastorello et al. 2013; Margutti et al. 2014) and SN 2015bh (Elias-Rosa et al. 2016). The comparison SNe have been corrected for extinction using published estimates and assuming the Cardelli et al. (1989) extinction law (see also Table 6). SN 2010bt exhibits a rapid decline at early times, similar to the decline from the "b" event of SN 2009ip and slightly faster than SN 1996al. This decline rate is also reminiscent of rapidly declining SNe II or even of SNe IIb (see, e.g., Li et al. 2011; Anderson et al. 2014). The absolute V magnitude at maximum of SN 2010bt was −19 mag or brighter, consistent with both the typical V-band peak magnitudes of SNe IIn (MV = −18.4 ± 1.0 mag; Kiewe et al. 2012) and those of SNe II-P/L (MV = − 16.89 ± 0.98 mag; Galbany et al. 2016).
Download figure:
Standard image High-resolution imageTable 6. Properties of the Comparison Supernovae
SN | Host Galaxy | Redshift | Distancea | E(B − V)tot | Vmax Date | Source |
---|---|---|---|---|---|---|
(Mpc) | (mag) | (MJD) | ||||
1996al | NGC 7689 | 0.007 | 22.9 | 0.110 | 50,265.0 | a |
1998S | NGC 3877 | 0.002 | 15.7 | 0.219 | 50,890.5 | b |
2009ip | NGC 7259 | 0.006 | 25.0 | 0.019 | 56,203.5 | c |
2015bh | NGC 2770 | 0.007 | 29.3 | 0.208 | 57,167.0 | d |
2010bt | NGC 7130 | 0.016 | 65.4 (4.6) | 0.40 (0.14) | 55,293.1b | This work |
Notes.
aDistances have been scaled to H0 = 73 km s−1 Mpc−1. bWe have assumed V maximum date = discovery date (MJD 55,303.1) − 10 days. Sources: (a) Benetti et al. (2016); (b) Fassia et al. (2000), Leonard et al. (2000), Pozzo et al. (2004), NED; (c) Pastorello et al. (2013), Fraser et al. (2013), Mauerhan et al. (2013), Margutti et al. (2014); (d) Elias-Rosa et al. (2016).Download table as: ASCIITypeset image
The early-time (B − V)0, , and color curves of SN 2010bt (see Figure 4) exhibit rapid evolution from blue to red, similar to the color evolution of the other SNe IIn, yet different from SN 1996al.
Download figure:
Standard image High-resolution imageWe have computed the pseudobolometric light curve of SN 2010bt by integrating the flux at different wavelengths derived from the extinction-corrected optical apparent magnitudes over the sparse observations in the UVW1 through K bands. Fluxes were measured considering only the epochs when V-band observations were available. When photometric measurements in one band at given epochs were not available, the flux was estimated by interpolating magnitudes from epochs close in time or, when necessary, by extrapolating the missing photometry assuming a constant color. We estimated the pseudobolometric flux at each epoch by integrating the spectral energy distribution (SED) using the trapezoidal rule and assuming zero flux outside the integration boundaries. Finally, the effective fluxes were converted to luminosities using the adopted distance to the SN (see Section 1). The errors in the bolometric luminosity include the uncertainties in the distance estimate, the extinction, and the apparent magnitudes. Note that the UVW1 band provides about 13% of the luminosity.
In Figure 5 we present the pseudobolometric light curve of SN 2010bt, as well as those of SN 1996al, SN 1998S, SN 2009ip, and SN 2015bh computed with the same prescriptions (for SN 1996al and SN 1998S we have no observations in UV bands). As shown earlier, the luminosity decline of SN 2010bt is similar to that of SN 2009ip, while it was more luminous than all other SNe in the sample, with the possible exception of SN 1998S. Considering the first epoch in the V band for SN 2010bt, its peak may have reached a luminosity >1.3 × 1043 erg s−1 (SN 1998S had a luminosity at peak of 1.6 × 1043 erg s−1).
Download figure:
Standard image High-resolution imageAssuming that the tail of SN 2010bt followed the radioactive 56Co decay with full trapping, and considering the explosion date to be between 15 and 50 days before the discovery date (see Section 2.3) and the HST detection in F606W (∼V) at 175.3 days (from discovery), we can roughly estimate ≤0.03 ± 0.01 M☉ for the 56Ni mass, using the formula given by Hamuy (2003). This value is at the low end in comparison with typical values for CC-SNe, 0.001–0.3 M☉ (Hamuy 2003). The 56Ni mass of SN 2010bt is also in agreement with what was estimated for the interacting SN 1996al, SN 2009ip, and SN 2015bh (MNi ≤ 0.02, ≤0.08, and ≤0.04 M⊙, respectively; Fraser et al. 2013; Margutti et al. 2014; Smith et al. 2014; Benetti et al. 2016; Elias-Rosa et al. 2016).
2.2. Spectroscopy
Four low-resolution optical spectra of SN 2010bt were obtained on 2010 April 18 (classification spectrum), September 15 and 16, and October 28 with the 3.58 m NTT+EFOSC2 at ESO of La Silla (Chile). Only the first spectrum was taken with the slit along the parallactic angle to avoid differential flux losses (Filippenko 1982). For the other spectra it was necessary to point to a nearby bright star and then rotate the slit to position the SN 2010bt site inside the aperture. Basic information on our spectra is reported in Table 7.
Table 7. Log of Spectroscopic Observations of SN 2010bt
UT Date | MJD | Phasea | Instrument | Grism/Grating + slit | Spectral Range | Resolution |
---|---|---|---|---|---|---|
(days) | Key | (Å) | (Å) | |||
20100418 | 55304.4 | 1.3 | EFOSC2 | gm11+gm16+15 | 3400–10,300 | 14 |
20100915 | 55454.1 | 151.0 | EFOSC2 | gm11+gm16+10 | 3400–10,300 | 14 |
20100916 | 55455.1 | 152.0 | EFOSC2 | gm11+10 | 3400–7500 | 14 |
20101028 | 55497.1 | 194.0 | EFOSC2 | gm11+10 | 3400–7500 | 14 |
Note.
aPhases are relative to the date of discovery, MJD = 55,303.1.Download table as: ASCIITypeset image
All spectra were reduced following standard procedures with iraf routines. The 2D spectroscopic frames were debiased and flat-fielded, before the optimized extraction (Horne 1986) of the 1D spectra. Wavelength calibration was accomplished with the help of arc-lamp exposures obtained the same night. Small adjustments estimated from night-sky lines in the object frames were applied. The spectra were flux-calibrated using the well-exposed continua of spectrophotometric standard stars (Oke 1990; Hamuy et al. 1992, 1994). An atmospheric extinction correction was applied using tabulated extinction coefficients for the ESO-La Silla Observatory. The strongest telluric absorption bands were removed using the standard-star spectra. Finally, the flux of each spectrum was cross-checked against the photometry.
Figure 6 shows the sequence of optical spectra of SN 2010bt,20 and in Figures 7 and 8 we compare some of these spectra with those of SNe IIn SN 1996al (Benetti et al. 2016) and SN 1998S (Leonard et al. 2000) and the transients SN 2009ip (Fraser et al. 2013; Pastorello et al. 2013) and SN 2015bh (Elias-Rosa et al. 2016) at similar epochs. All of the spectra have been corrected for extinction and deredshifted using values from the literature (see also Table 6).
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageThe early-time spectrum (Figure 6, top panel) exhibits a blue continuum, with relatively weak spectral features, except for the strong Balmer emission lines with complex yet relatively narrow profiles. SN 2010bt does not show signs of the blue pseudocontinuum that characterizes some of the most energetic SNe IIn such as SN 1988Z (Turatto et al. 1993; Kiewe et al. 2012). The pseudocontinuum is a "blue excess" arising from the overlap of emission lines when the expanding ejecta of the SN interact with the surrounding CSM. The lack of a pseudocontinuum indicates that for SN 2010bt the continuum is thermal with a blackbody-like behavior.
Balmer emission lines show P Cygni absorption components with expansion velocities of 4000–3500 km s−1, estimated from their absorption minima. A blend of He i λ5876 and Na i D is also visible in the spectrum, with the broad emission centered at ∼5950 Å and with a very broad FWHM ≈ 290 Å. The photospheric temperature at this epoch, estimated by fitting the SED of SN 2010bt with a blackbody function, is around 12,900 K (see also Table 8). Given this temperature and the pseudobolometric luminosity of SN 2010bt at that epoch, a radius of the photosphere of ∼8 × 1014 cm is inferred.
Table 8. Main Parameters Inferred from Spectra of SN 2010bt
UT Date | MJD | Phasea | Temp.b | Radiusc | FWHMHα,nard | FWHMHα,br | ||
---|---|---|---|---|---|---|---|---|
(days) | (K) | (1014 cm) | (km s−1) | (km s−1) | (km s−1) | (1039 erg s−1) | ||
20100418 | 55,304.4 | 11.3 | 12900 | 8 (2) | 750 (100) | 12,800 (1600) | 3000 (800) | 50 (7) |
FWHMBlue | FWHMRest | |||||||
(km s−1) | (km s−1) | |||||||
20100915 | 55,454.1 | 161.0 | ⋯ | ⋯ | 1500 (400) | 1050 (300) | ⋯ | 3 (1) |
20100916 | 55,455.1 | 162.0 | ⋯ | ⋯ | 1550 (400) | 750 (500) | ⋯ | 2 (1) |
20101028 | 55,497.1 | 204.0 | ⋯ | ⋯ | 1450 (400) | 1100 (350) | ⋯ | 1 (0.4) |
Notes.
aPhases are relative to the assumed V maximum date = discovery date (MJD 55,303.1) − 10 days. bWe consider a conservative uncertainty in the temperature of about ±500 K. cWe have propagated the uncertainties from the Stefan–Boltzmann equation. dThe velocities are computed from the decomposition of the Hα profile.Download table as: ASCIITypeset image
Consistent with the photometry, the Hα profile of SN 2010bt resembles those of the interacting transients SN 2009ip and SN 2015bh and shows differences with SN 1996al and SN 1998S (see the close-up view of the Hα profiles in Figures 7 and 8). The blackbody temperature shaping the continuum, though, is similar to that of SN 1998S, as are the strengths of the He i lines.
Download figure:
Standard image High-resolution imageThe extractions of the late-time spectra (>150 days) were arduous, given their low signal-to-noise ratio (S/N). Note that these spectra are contemporaneous with the upper photometric limits measured for SN 2010bt. The only distinct features are those of Hβ and Hα, together with typical residual emission lines, such as [N ii] λλ6548, 6584 and [S ii] λλ6717, 6731, from neighboring H ii regions. Despite the low S/N of the late-time SN 2010bt spectra, we notice some differences comparing its Hα profile with that of the other SNe. None of the transients share the same profile shape, with that of SN 2010bt being double peaked yet narrower than that of the other objects (see the Hα profiles magnified in Figure 8).
2.2.1. Hydrogen Feature Decomposition
The Hα line profile of SN 2010bt appears to consist of multiple components with evident changes in morphology between early and late phases. We have decomposed the line into a Lorentzian profile for the narrow component and two Gaussian functions for the emission and absorption components of the broader P Cygni profile at the early epochs, and into two Lorentzian profiles for the oldest phases. To do this, we have used a Python script for least-squares minimization.
Figure 9 presents the Hα emission line decomposition for three of our four epochs. The procedure independently fits the parameters, while the uncertainties were estimated using a bootstrap resampling technique, varying randomly the flux of each pixel according to a normal distribution having variance equal to the noise of the continuum. The velocities of the different gas components are listed in Table 8.
Download figure:
Standard image High-resolution imageThe relatively narrow Hα emission at the first epoch (phase ∼ 11.3 days) is resolved with an FWHM of ∼750 km s−1 (after correction for instrumental resolution), while the broader component has FWHM ≈ 12,800 km s−1. In addition, the blue side of the Hα profile is "absorbed" by a P Cygni component with minimum at ∼3000 km s−1. The broader component is redshifted, probably caused by electron scattering as the Hα line photons diffuse through dense CSM ahead of the shock. The observed luminosity of Hα was estimated from the integrated flux of the entire line to be 5 × 1040 erg s−1.
At later phases (>150 days), the Hα profile is well reproduced with two Lorentzian components called "Blue" and "Rest," given that they are centered at wavelengths of ∼6542 and ∼6565 Å, respectively. Clearly visible is a broad feature of Hα in the 2D images of the late-time SN 2010bt spectra (see Figure 10). It is also evident that these spectra are likely contaminated by neighboring H ii regions. In fact, a residual emission line centered at ∼6588 Å, most likely corresponding to [N ii] λ6584, is present. Consequently, we have considered this line in the decomposition of Hα by adding a third Lorentzian component. The FWHMs of the Blue and Rest components at all late epochs are relatively constant, at ∼1500 and ∼1000 km s−1, respectively. The total luminosity of the Hα line at these phases decreases to 3 × 1039 erg s−1.
Download figure:
Standard image High-resolution image2.3. Explosion Date of SN 2010bt
Monard (2010) reported that no sources were visible at the SN position in images taken on 2009 December 22.80 (limit >17.8 mag, or MR ≳ −17.3 mag, for this particular case). Nothing was also visible in a HAWK-I image taken on 2009 July 26.60 (limit K > 19.0 mag; in this case MK ≳ −15.2 mag; Miluzio et al. 2013). Therefore, the SN explosion occurred less than 115 days before discovery.
A more accurate date of the explosion is difficult to estimate in the case of SN 2010bt. The main complications are the peculiar behaviors of both its photometric and spectroscopic evolution.
The light curves (see Section 2.1) show that the SN was discovered after maximum light. Assuming that the V-band light curve of SN 2010bt has a similar decay slope to that of SN 2009ip (see Section 2.1 and Figure 3), the former would have been discovered no more than 10 days after maximum light. This estimate is also confirmed by the behavior of the color curves. Considering that SNe IIn generally have a rise time of >5 days (Ofek et al. 2014), SN 2010bt would have exploded at least 15 days before the discovery. On the other hand, and according to the template-fitting code GELATO, acceptable matches with SNe IIn at phases between 10 and 50 days after explosion are found for the early-time SN 2010bt spectrum (although the best fit is at 10 days).
Thus, in view of all these uncertainties, we are able to determine only that the explosion occurred 15–50 days before discovery.
3. Identification of the Progenitor Candidate
In order to search for a possible SN progenitor, we isolated archival21 HST images of NGC 7130 taken with the Wide Field and Planetary Camera 2 (WFPC2) in F606W (∼V; 500 s exposure) on 1994 August 23 by program GO-5479 (PI: M. Malkan) and with the High Resolution Channel (HRC) of the Advanced Camera for Surveys (ACS) in F330W (∼u; 1200 s exposure) on 2003 May 08 by program GO-9379 (PI: H. Schmitt). We worked with drizzled images downloaded from the Hubble Legacy Archive.22 These images have been resampled to a uniform grid to correct for geometric distortions.
We performed relative astrometry by geometrically transforming the pre-explosion images to match those taken after the explosion. We first used the ground-based, post-explosion CTIO/SMARTS image taken on 2010 April 23 (with seeing 15) to approximately locate the position of the SN in the WFPC2 images. Then, we confirmed the identification of this candidate through high-resolution HST ACS images using the Wide Field Channel (WFC, ∼005 pixel−1; a pair of observations of 40 and 80 s exposures), obtained through the F606W filter on 2010 October 09, as part of our Target-of-Opportunity program GO-11575 (PI: S. Van Dyk). The individual exposures of this trigger were drizzled to produce a final mosaic. However, as we showed in Section 2.1, the SN was already quite faint at the time of these observations ( mag; see Table 3), leaving some ambiguity in its identification. Therefore, we used the HST post-explosion image to identify the progenitor in the pre-SN WFPC2 and ACS/HRC images, as follows.
(i) Initially, we obtained a precise SN position, at α = 21h48m2028, δ = −34°57'177 (J2000.0), with rms uncertainties of ∼01, based on nine point-like sources as fiducials in the V-band NTT+EFOSC2 image taken on 2010 April 18. We measured the centroids of the fiducial sources with the iraf tasks imexamine. We adopted the 2MASS Point Source Catalog as the astrometric grid and used the iraf task ccmap to obtain the solution. At that point, we were able to confirm the SN position in the HST post-explosion mosaic, with rms uncertainty <006, and located the faint SN.
(ii) Then, we achieved high-precision relative astrometry between the pre-explosion WFPC2 F606W (we worked only with the drizzled PC image, since the SN site is located on this chip, with 0045 pixel−1) and ACS/HRC F330W images (with ∼0025 pixel−1) and the post-explosion ACS/WFC image, by geometrically transforming the former pair to match the latter. We used 17 point-like sources in common between the three sets of images and the iraf task geomap for the transformation.
The positions (and their uncertainties) for the SN and the progenitor candidate are obtained by averaging the measurements from two centroiding methods, the task daofind within iraf/daophot and imexamine within iraf. As a result, in the pre-SN F606W image we identify an object very close to the SN position, which we consider to be the progenitor candidate. This same object is faintly detected at this position in the F303W image; see Figures 11(a) and (b). The pixel position for the SN transformed into the pre-SN F606W image is [242.83, 363.79], while the candidate position is [243.08, 363.96]; for the pre-SN F330W image, these are [348.36, 317.88] and [348.92, 318.02], respectively. The differences between the SN and the progenitor candidate positions, compared with the total estimated astrometric uncertainty, are given in Table 9. The latter uncertainty was calculated as a quadrature sum of the uncertainties in the SN and progenitor candidate positions and the rms uncertainty in the geometric transformation.
Download figure:
Standard image High-resolution imageTable 9. SN 2010bt and the Progenitor Candidate Position Comparison
F330W | F606W | F555W | |
---|---|---|---|
(α/δ) | (α/δ) | (α/δ) | |
Total uncertainty (mas) | 10/10 | 9/8 | 10/10 |
Diff. position SN/candidate (mas) | 28/7 | 13/8 | ⋯ |
Note. Uncertainties are 1σ.
Download table as: ASCIITypeset image
Table 10. Brightnessa of SN 2010bt and the Progenitor Candidate
F330W | F606W | F555W | |
---|---|---|---|
(mag) | (mag) | (mag) | |
Progenitor candidate | 24.29 (0.23)b | 22.23 (0.02)c | ⋯ |
SN 2010bt | ⋯ | 22.50 (0.05)d | ≳24.4e |
Notes.
aMagnitude uncertainties are 1σ. bImage taken with HST+ACS/HRC on 2003 May 08. cImage taken with HST+WFPC2 on 1994 August 23. dImage taken with HST+ACS/WFC on 2010 October 09. eImage taken with HST+WFC3/UVIS on 2015 October 14.Download table as: ASCIITypeset image
From the results in Table 9, it can be seen that the difference between the SN position and the position of the progenitor candidate is slightly larger in right ascension for both bands than the measurement uncertainties; the agreement is much better in declination. However, the differences in position are within 3σ of the total uncertainties. We therefore consider the candidate as the SN progenitor star and attribute the larger offset in right ascension (about 28 and 13 mas in F330W and F606W, respectively) to the complex background in the bright spiral arm of the host galaxy on which the progenitor candidate is located. We emphasize that no other point-like source exists within this 3σ radius from the progenitor candidate position (see Figure 11(c)).
The field was observed again on 2015 October 14 as part of our Target-of-Opportunity program GO-13683 (PI: S. Van Dyk). On this occasion we obtained six dithered images of 100 s exposure with the HST Wide Field Camera 3 (WFC3) UVIS channel (∼004 pixel−1) in F555W band (hereafter labeled as "late-time HST images"). Using relative astrometry, we matched the observations of SN 2010bt before and after the explosion with this late-time drizzled HST image. We did not detect any luminous point source at the SN position in the F555W image (Figure 11(d)), confirming the identification of the progenitor (see next section for further details).
4. The Nature of the Progenitor Candidate
As discussed in Section 2.2, the early-time spectrum of SN 2010bt has a rather blue continuum, yet it also exhibits relatively strong, narrow Na i D absorption, suggesting the simultaneous presence of dust and extinction suffered by the SN. Measuring the equivalent width (EW) of the blended Na i D line at the host-galaxy redshift (z0 = 0.016) from the early-time optical spectrum of SN 2010bt (∼1.8 Å), we can attempt to estimate . Any relationship between EW and tends to break down for SNe with moderate to high reddening. We obtain a large dispersion in the results following the relations of different authors. Specifically, we estimate and 0.9 mag (i.e., AV ≈ 0.8 and 2.7 mag, respectively—assuming the Cardelli et al. [1989] reddening law with updated wavelengths and RV = 3.1) following Turatto et al. (2003) and mag (i.e., AV ≈ 5.0 mag) using Poznanski et al. (2012). Because of the large differences obtained with the EW(Na i D) versus relations, we decided not to use this method to estimate the extinction toward SN 2010bt. Phillips et al. (2013) have also cautioned that the EW(Na i D) versus relationship is associated with large scatter.
Thus, in order to estimate a consistent value for the extinction (AV) toward SN 2010bt, we considered different methods based on comparisons of the object's SED and luminosity with those of other "standard" SNe IIn. It should be noted that in fact SNe IIn exhibit a wide diversity in both photometry and spectra (see, e.g., Kiewe et al. 2012; Taddia et al. 2013). We first matched simultaneously the intrinsic , , and color curves of SN 2010bt with those of other SNe IIn (see Figure 4), finding an extinction AV = 1.47 ± 0.31 mag. We also compared the early-time optical SED of SN 2010bt (we did not used the late-time spectra, given their comparatively poor S/N) with those of the interacting SN 1998S, SN 2009ip, and SN 2015bh at similar epochs. Spectra of the comparison SNe were first corrected for redshift and extinction and then scaled to the distance of SN 2010bt. The average of the good matches in all cases is AV = 0.98 ± 0.33 mag. We adopt the uncertainty-weighted average value, mag (i.e., mag), as the extinction toward SN 2010bt.
We measured the brightness of the progenitor candidate using Dolphot, finding mF606W = 22.23 ± 0.02 mag (also reported in Table 9). Dolphot reported that the progenitor candidate has an "object type" flag of "1," which means that the source is likely stellar. Additionally, the object had a "sharpness" and a "crowding" parameter23 ∼0, further indicating that the detected source is point-like and "clean" (Dolphin 2000).
Dolphot also detected the candidate in F330W, although the object type is "2," indicating that the source is a star too faint for a PSF determination; the other two parameters indicated that this object is an otherwise good and clean source. The magnitude measured (via aperture photometry) in this filter is mF330W = 24.29 ± 0.23 mag. Note that the HST images taken with F606W were obtained in 1994 and those using F330W in 2003.
We obtained mF606W = 22.50 ± 0.05 mag for the SN, running Dolphot on the HST images taken on 2010 (assuming that the F606W bandpass is approximately Johnson V, we include this measurement in Table 3). The SN was therefore found to be slightly fainter than the progenitor candidate.
Finally, we ran Dolphot on the 2015 late-time HST images of the SN field. By adding artificial stars at the SN position, we obtained conservatively an upper brightness limit for a single point source with a sigma threshold of 3.0 of mF555W ≳ 24.4 mag. That is, the previously identified source had disappeared, being at least 2 mag fainter than in the observations in 1994. The precise match between the HST images taken before and after the SN explosion, the fit quality given by Dolphot, and the disappearance of the progenitor star 5 yr after the explosion tell us that this star is likely the progenitor of SN 2010bt. Note that the 2015 HST images were obtained with WFC3, which has better throughput than WFPC2. SN 2009ip was also observed at late time phases (>1100 days after the maximum of its brightest event 2012b), at which time the luminosity of the transient was probably still dominated by CSM interaction, and was found to be marginally fainter than its quiescent progenitor (Thöne et al. 2015; Smith et al. 2016; Graham et al. 2017). Unfortunately, like SN 2009ip, then, this apparent fading of SN 2010bt 15 yr after the nominally SN-like event cannot alone prove that it was a terminal explosion.
Correcting by the total extinction and distance assumed for the SN (see the beginning of this section and Section 1), we find that the absolute magnitude of the progenitor was mag and mag, while the SN was at −12.72 ± 0.42 mag (F606W) in 2010 and ≳−11 mag (F555W) in 2015.
Adopting zero bolometric correction and that mF606W ≈ V, we find that the progenitor candidate has a bolometric luminosity of log(L/L☉) ≈ 7.1 ± 0.2. Such a luminosity is too large for a star in quiescence, and thus it is more likely that the progenitor was observed in eruption at that epoch.24 Curiously, this value of luminosity is evocative of η Carinae during its Great Eruption (Humphreys et al. 1999; Rest et al. 2012; see also Figure 12).
Download figure:
Standard image High-resolution imageIt is well known that the surrounding medium of a star can affect the peak luminosity and spectra during the explosion. Indeed, SN 2010bt was classified as Type IIn owing to the presence of narrow circumstellar H emission. A possible explanation for the fact that the SN light faded so drastically in just a few months after the explosion could be the result of possible formation of dust in the ejected material, which could have obscured the SN light (or the "remnant" if it survived the explosion). In this case, we should expect at late times a strong IR excess, as in the case of SN 1998S (Pozzo et al. 2004) or SN 2005ip (e.g., Smith et al. 2009; Fox et al. 2010). Therefore, we analyzed several sets of IR images obtained with the InfraRed Array Camera (IRAC) on board the Spitzer Space Telescope, before (2004 October 30; 3.6, 4.5, 5.8, and 8 μm channels; PI: J. M. Mazzarella; Program ID 3672) and after (2012–2015; 3.6, and 4.5 μm channels; PIs D. B. Sanders, O. D. Fox, D. Stern; Program IDs 80089, 90031, 10098, respectively) the SN explosion. We worked with the post-basic calibrated data (PBCD) provided by the Spitzer Heritage Archive,25 which are already fully co-added and calibrated. Neither the progenitor candidate nor SN 2010bt was detected in any Spitzer channel (see Figure 13). We therefore attempted to constrain the dust formation at the SN site by measuring the variance of the integrated flux in an area of 3 × 3 pixels around the progenitor/SN position using MOPEX. Following the "recipe" advised by the Spitzer team, we obtained the upper limits reported in Table 11.
Download figure:
Standard image High-resolution imageTable 11. Spitzer Flux Upper Limits at the SN 2010bt Site
Date | MJD | Phasea | 3.6 μm | 4.5 μm | 5.8 μm | 8.0 μm | Program ID |
---|---|---|---|---|---|---|---|
(days) | (μJy) | (μJy) | (μJy) | (μJy) | |||
20041030 | 53,308.2 | −1994.9 | <195.5 | <196.2 | <749.2 | <2560.0 | 3672 |
20121221 | 56,282.7 | 979.5 | <202.6 | <209.7 | ⋯ | ⋯ | 80089 |
20121221 | 56,282.9 | 979.8 | <169.3 | ⋯ | ⋯ | ⋯ | 90031 |
20130123 | 56,315.5 | 1012.4 | <205.3 | ⋯ | ⋯ | ⋯ | 90031 |
20130718 | 56,491.4 | 1188.3 | <177.3 | ⋯ | ⋯ | ⋯ | 90031 |
20130820 | 56,524.2 | 1221.1 | <238.7 | ⋯ | ⋯ | ⋯ | 90031 |
20131224 | 56,650.4 | 1347.3 | <223.1 | ⋯ | ⋯ | ⋯ | 90031 |
20140131 | 56,688.6 | 1385.5 | <210.1 | ⋯ | ⋯ | ⋯ | 90031 |
20140727 | 56,865.9 | 1562.8 | <193.8 | ⋯ | ⋯ | ⋯ | 90031 |
20140829 | 56,898.4 | 1595.3 | <197.5 | ⋯ | ⋯ | ⋯ | 90031 |
20150103 | 57,025.9 | 1722.8 | <237.6 | <177.5 | ⋯ | ⋯ | 10098 |
Note.
aPhases are relative to the discovery date, MJD = 55,303.1.Download table as: ASCIITypeset image
As one can see, there is no evidence of IR excess emission from years 2004 through 2014. The upper limits on the IR emission from SN 2010bt tell us that there are no detections larger than ∼9 × 1040 erg s−1 and ∼7 × 1040 erg s−1 at 3.6 and 4.5 μm, respectively. Thus, it is quite likely that the drop in luminosity of SN 2010bt is not caused by dust formation, but rather by the end of the ejecta interaction with the CSM.
5. Was SN 2010bt a Nonterminal Explosion?
As discussed by many authors, the dividing line between SN impostors and real SNe is not clear, with both possibilities viable for several objects. Both families of transients share observable similarities; however, the big difference is in the progenitor star fate—the star is destroyed in an SN explosion, but it still remains after the outburst in the case of an SN impostor. It is also hard to find a unique origin of SN 2010bt. We discuss below the observables of this transient.
- 1.Both early and late spectra show spectra dominated by strong multicomponent Balmer emission lines, indicators that define the interacting transients. The Hα profile of SN 2010bt at early time resembles more those of the interacting transients, such as SN 2009ip and SN 2015bh, than those of confirmed SNe IIn SN 1996al and SN 1998S (see Section 2.2). Instead, at late phases, none of the comparison transients share the same Hα profile as SN 2010bt. This is understandable since the CSM may evolve in a different way for each object, leading to diversity among the transients.
- 2.Photometrically (see Section 2.1), SN 2010bt exhibits a rapid decrease in luminosity after maximum, which is reminiscent of the SN 2009ip-like objects, and is somehow faster than other SNe. However, the luminosity at maximum of SN 2010bt is higher than that of SN 2009ip and most consistent with the peak luminosity of bright SNe IIn.
- 3.The flattening of the late-time light curve seems slower than that expected for interacting SNe; however, it is still at the low end of the typical 56Ni mass range for CC-SNe. It is, however, in agreement with what was estimated for SN 2009ip-like SNe.
- 4.
- 5.Unfortunately, except for the estimated magnitude of the progenitor star, we have not found any other data indicative of a possible pre-SN eruptive event.
- 6.We estimated an absolute magnitude for the progenitor star of mag (see Section 4), which is similar to that derived for the SN in 2010, −12.7 mag (F606W). That luminosity is more like that of a star in outburst, if we compare with SN impostors, such as UGC 2773-OT (Smith et al. 2010), SNhunt248 (Kankare et al. 2015), or PSN J09132750+7627410 (Tartaglia et al. 2016), which were all between −13 and −14 mag, or objects in outburst, such as SN 2009ip (e.g., Pastorello et al. 2013), SN 2015bh (Elias-Rosa et al. 2016; Thöne et al. 2017), or η Carinae during its Great Eruption (Humphreys et al. 1999; Rest et al. 2012). Note also that the progenitor luminosity found for the SN IIn 2005gl (Gal-Yam et al. 2007; Gal-Yam & Leonard 2009), or some SN 2009ip-like transients (Smith et al. 2010, 2014; Foley et al. 2011; Elias-Rosa et al. 2016; Pastorello et al. 2018), is around ∼−10 mag.
- 7.Several years after the discovery, the progenitor of SN 2010bt seems to have vanished, indicating that the star may have finally exploded as an SN (see Section 4). This has been inferred for interacting SNe, such as SN 2005gl, but not yet for any member of the SN 2009ip-like family. As discussed by Pastorello et al. (2018), SN 2009ip-like transients share similarities with SN 2005gl, supporting the possible terminal explosion scenario for all.
In short, the SN 2010bt observations have not helped us clarify whether the transient is an SN or an impostor, with the exception of the very late-time HST data that may show that the progenitor had vanished. It makes the terminal explosion scenario plausible for this event.
6. Summary
SN 2010bt was classified as an SN IIn from an optical spectrum taken not long after the explosion. The observational campaign was interrupted after ∼2 months, owing to the SN's proximity to the Sun in the sky, and was again continued just a few months later. By that time, the SN had become much fainter or was nearly undetectable, which is unusual behavior for an SN in general. SN detection at late phases was obtained only through images taken with HST, at a luminosity of ∼−12.7 mag in the F606W band.
Comparing HST images of the host galaxy prior to the explosion and those of the SN at late times, high-precision relative astrometry allowed us to identify the likely SN progenitor star with magnitude ∼−13 in the F606W band, which is almost certainly the luminosity of a star in eruption. At first we found that the brightness of the SN nearly 5 months after discovery was somewhat smaller than that estimated for the progenitor candidate. Still, no source more luminous than −11 mag in the F555W band was found at the position of the SN 5 yr after the SN discovery.
Overall, we have not found a unique explanation for the chain of events observed for SN 2010bt. In the following we will list our best understanding of the transient's observables.
- 1.The SN 2010bt progenitor, identified in the HST pre-SN images, was in outburst ( mag; log(L/L⊙) ≈ 7) at the moment of the observations in 1994. Unfortunately, no information about other possible pre-SN eruptive events has been found.
- 2.Some time thereafter, a powerful (terminal or nonterminal) outburst occurs, resulting in SN 2010bt.
- 3.The ejecta interacted with a compact shell created during eruptions or heavy mass loss from the massive progenitor star prior to the explosion. This shock/CSM interaction led to SN 2010bt reaching a quite luminous MV ≥ −19 mag.
- 4.A probably patchy photosphere is located in the external symmetric CSM, where the CSM–ejecta interaction is taking place, since the early-time observed Hα profile is quite symmetric and a weak broad component is visible.
- 5.The CSM+ejecta recombined quickly, leading to the observed rapid decline of the SN 2010bt light curve.
- 6.The ejecta continued to propagate into more distant CSM. This probably had an asymmetric geometry, since the late-time Hα emission showed a double-peaked profile.
- 7.Over time, the CSM–ejecta interaction became far less strong; consequently, we saw no further trace of SN 2010bt (at least not brighter than −11 mag) in 2015.
In conclusion, we confirm the identification of the SN 2010bt progenitor in outburst and favor the scenario in which SN 2010bt was a genuine SN surrounded by dense CSM with a complex geometry.
N.E.-R. thanks Avet Harutyunyan for his help and acknowledges the hospitality of the Institut de Ciències de l'Espai at the Autonomous University of Barcelona's Campus, where much of this work was written. Support for this work was provided by NASA/HST through grants GO-11575, GO-13684, GO-14668, and AR-14295 from the Space Telescope Science Institute (STScI), which is operated by AURA, Inc., under National Aeronautics and Space Administration (NASA) contract NAS5-26555. S.B. and M.T. acknowledge partial financial support by the PRIN-INAF 2017 (project "Towards the SKA and CTA Era: Discovery, Localization, and Physics of Transient Sources"). A.V.F. is grateful for generous financial assistance from the Christopher R. Redlich Fund, the TABASGO Foundation, the Miller Institute for Basic Research in Science (UC Berkeley), and US NSF grant AST-1211916. G.P. acknowledges support provided by the Millennium Institute of Astrophysics (MAS) through grant IC120009 of the Programa Iniciativa Cientíifica Milenio del Ministerio de Economía, Fomento y Turismo de Chile. L.G. was supported in part by NSF grant AST-1311862. N.S. was supported in part by NSF grants AST-1312221 and AST-1515559. This research is based in part on observations made with the NASA/ESA Hubble Space Telescope and obtained from the Hubble Legacy Archive, which is a collaboration between STScI/NASA, the Space Telescope European Coordinating Facility (ST-ECF/ESA), and the Canadian Astronomy Data Centre (CADC/NRC/CSA); the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA (support was provided by NASA through an award issued by JPL/Caltech); the SMARTS Consortium 1.3 m telescope located at Cerro Tololo Inter-American Observatory (CTIO), Chile; and the New Technology Telescope at the European Southern Observatory–La Silla Observatory. Data in this work have been taken in the framework of the European supernova collaboration involved in ESO-NTT large program 184.D-1140 led by Stefano Benetti. Many of the spectra used here for comparison were obtained from the Padova-Asiago Supernova Archive (ASA). This work has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with NASA, and data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by NASA and the NSF.
Facilities: Swift - Swift Gamma-Ray Burst Mission, Spitzer - Spitzer Space Telescope satellite, ING:Kapteyn (JAG) - Isaac Newton Group's Jacobus Kapteyn Telescope, CTIO:1.3m (ANDICAM) - Cerro Tololo Inter-American Observatory's 1.3 meter Telescope, HST(WFPC2) - Hubble Space Telescope satellite, ACS - , WFC3 - , ESO:3.6m (EFOSC2 - , SOFI) - , Magellan:Clay (LDSS-3) - Magellan II Landon Clay Telescope, VLT:Yepun (HAWK-I) - Very Large Telescope (Yepun).
Software: IRAF (Tody 1986, 1993), SNOoPy (http://sngroup.oapd.inaf.it/snoopy.html), SExtractor (Bertin & Arnouts 1996), HOTPANTS (Becker 2015), Dolphot (Dolphin 2000), GELATO (https://gelato.tng.iac.es), MOPEX (http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/mopex/).
Footnotes
- 17
NED, NASA/IPAC Extragalactic Database; http://www.ipac.caltech.edu/.
- 18
- 19
Considering the interval from ∼1 to 50 days after the discovery date.
- 20
Our spectra are available on WiseREP (Yaron & Gal-Yam 2012).
- 21
- 22
- 23
The "sharpness" parameter indicates the reliability that a detected source is indeed point-like. It is considered a "good star" when this parameter is between −0.3 and +0.3 and a "perfectly fit star" when the value is 0. The "crowding" parameter is defined in magnitude and describes the measured quality brightness of a star. Isolated stars have a crowding value of zero. This value increases as more stars surround the star under study, contaminating the measurement. Hence, a low value increases the certainty of the measured stars.
- 24
Assuming mag, the bolometric luminosity of SN 2010bt's progenitor is log(L/L☉) ≈ 6.8 ± 0.1. This luminosity is still more appropriate for a star in eruption.
- 25