Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

A publishing partnership

SN 2012fr: Ultraviolet, Optical, and Near-infrared Light Curves of a Type Ia Supernova Observed within a Day of Explosion*

, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2018 May 18 © 2018. The American Astronomical Society. All rights reserved.
, , Citation Carlos Contreras et al 2018 ApJ 859 24 DOI 10.3847/1538-4357/aabaf8

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/859/1/24

Abstract

We present detailed ultraviolet, optical, and near-infrared light curves of the Type Ia supernova (SN) 2012fr, which exploded in the Fornax cluster member NGC 1365. These precise high-cadence light curves provide a dense coverage of the flux evolution from −12 to +140 days with respect to the epoch of B-band maximum (${t}_{{B}_{\max }}$). Supplementary imaging at the earliest epochs reveals an initial slow and nearly linear rise in luminosity with a duration of ∼2.5 days, followed by a faster rising phase that is well reproduced by an explosion model with a moderate amount of 56Ni mixing in the ejecta. From our analysis of the light curves, we conclude that: (i) the explosion occurred <22 hr before the first detection of the supernova, (ii) the rise time to peak bolometric (λ > 1800 Å) luminosity was 16.5 ± 0.6 days, (iii) the supernova suffered little or no host-galaxy dust reddening, (iv) the peak luminosity in both the optical and near-infrared was consistent with the bright end of normal Type Ia diversity, and (v) 0.60 ± 0.15 M of 56Ni was synthesized in the explosion. Despite its normal luminosity, SN 2012fr displayed unusually prevalent high-velocity Ca ii and Si ii absorption features, and a nearly constant photospheric velocity of the Si ii λ6355 line at ∼12,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$ that began ∼5 days before ${t}_{{B}_{\max }}$. We also highlight some of the other peculiarities in the early phase photometry and the spectral evolution. SN 2012fr also adds to a growing number of Type Ia supernovae that are hosted by galaxies with direct Cepheid distance measurements.

Export citation and abstract BibTeX RIS

1. Introduction

Type Ia supernovae (SNe Ia) are the major producers of iron in the Galaxy (Tinsley 1979) and, therefore, they are intimately tied to its chemical evolution (e.g., see McWilliam 1997). The fact that they are among the brightest and most homogenous of the supernovae (SNe) has allowed them to be used with great success as distance indicators to measure the expansion history of the universe (e.g., see Betoule et al. 2014). Nevertheless, progress has remained disappointingly slow in identifying the progenitor systems of these objects and understanding the details of their explosion mechanisms. More than 50 years ago, Hoyle & Fowler (1960) recognized that SNe Ia were thermonuclear disruptions of white dwarfs in binary systems. However, there is still considerable debate as to whether the companion is a main-sequence or giant star (single-degenerate model), or is instead another white dwarf (double-degenerate model). The collision of two white dwarfs (as opposed to their merger) has recently garnered interest as a third possible way of producing SNe Ia (Raskin et al. 2009; Rosswog et al. 2009; Thompson 2011; Kushnir et al. 2013). As to the explosion mechanism, the deflagration (subsonic burning) of a Chandrasekhar mass white dwarf that at some point transitions to a supersonic detonation (Khokhlov 1991) provides the best match to the observational properties of SNe Ia (Höflich & Khokhlov 1996; Höflich et al. 1996; Hoeflich et al. 2017), but the details of how or why this occurs are still a mystery. Meanwhile, detonations of sub-Chandrasekhar-mass white dwarfs also have many attractive properties for explaining SNe Ia (Sim et al. 2010; Shen et al. 2018). These detonations can be triggered in a number of ways, including the double detonation mechanism (e.g., Woosley & Weaver 1994) from the accretion of a helium shell (e.g., Fink et al. 2010), a detonation in an accretion stream (Guillochon et al. 2010; Dan et al. 2012) in a violent merger involving a massive WD (Pakmor et al. 2012), or even in a more long-term evolution of a merger remnant (Shen et al. 2012).

Howell (2011) emphasized that large samples of SNe Ia help to improve our understanding of the progenitors and explosion mechanisms. For example, measurements of the relative rates in different types of galaxies can be used to infer the delay time distribution, which in turn can serve as a discriminant between different progenitor scenarios (Totani et al. 2008; Maoz et al. 2010). The high precision photometry and spectroscopy of individual SNe Ia also provides an insight into the progenitors and the physics of the explosion mechanism (e.g., see Höflich et al. 2010; Hsiao et al. 2013). Photometric observations of SNe Ia at the very earliest epochs following an explosion can be used to provide important constraints on the initial radius of the primary star (Piro et al. 2010; Bloom et al. 2012), the size of its companion star (Kasen 2010; Hosseinzadeh et al. 2017), and the distribution of 56Ni in the outermost ejecta and/or the possible presence of circumstellar material (Piro & Morozova 2016).

SN 2012fr was discovered by the TAROT (Télescopes à Action Rapide pour les Objets Transitoires) collaboration in images taken with their robotic telescope at the La Silla Observatory, Chile (Klotz 2012). With J2000.0 coordinates of α = 03h33m35fs99 and δ = −36°07'37farcs7, the transient was located 3'' west and 52'' north from the center of the SBb host-galaxy NGC 1365 (see Figure 1). Within a day and a half of the discovery, SN 2012fr was spectroscopically classified as a young, normal SN Ia, which was caught well before maximum light (Childress et al. 2012). Childress et al. (2013) analyzed an extensive set of optical spectral data of SN 2012fr, and documented the presence of high-velocity Si ii and Ca ii absorption components in the first spectrum obtained more than 14 days before ${t}_{{B}_{\max }}$, which persisted for the two weeks leading up to maximum. Although the presence of such high-velocity features (HVFs) at early epochs is common, it is relatively unusual for them to remain visible at maximum light (Mazzali et al. 2005). The evolution of the photospheric velocity as inferred from the Si ii λ6355 absorption was also unusually flat. Nevertheless, the overall spectral evolution of SN 2012fr was that of a normal SN Ia. Maund et al. (2013) published spectropolarimetry at four epochs between 11 days before and 24 days after ${t}_{{B}_{\max }}$. They found that although the continuum polarization of the SN was low (<0.01%) throughout this period, the Si ii and Ca ii HVFs were highly polarized at the earliest epoch. Zhang et al. (2014) also published extensive ultraviolet and optical photometry of SN 2012fr, along with optical spectroscopy, and they conjectured that SN 2012fr might be related to the 1991T-like events but was viewed from a direction where the ejecta presented a clumpy structure. More recently, Graham et al. (2017) presented Las Cumbres Observatory early-time gri and BVRI photometry of SN 2012fr, and Childress et al. (2015) and Graham et al. (2017) published nebular-phase spectra.

Figure 1.

Figure 1. Swope telescope color image of NGC 1365 with local sequence stars and SN 2012fr indicated.

Standard image High-resolution image

Between 2004 and 2009, the Carnegie Supernova Project (CSP-I, Hamuy et al. 2006) carried out precise optical and near-infrared (NIR) photometry of more than 100 nearby SNe Ia using the facilities at the Las Campanas Observatory (LCO). A second phase of the CSP (CSP-II) was initiated in 2011 to obtain the optical and NIR light curves of ∼100 SNe Ia located further in the Hubble flow at z ≃ 0.03–0.08 (M. M. Phillips et al. 2018, in preparation). In addition, CSP-II obtained light curves of a sample of ∼100 more nearby SNe Ia at z ≤ 0.04, for which detailed NIR spectroscopic time-series observations could be made. In this paper we present comprehensive optical and NIR light curves of one of the latter objects, SN 2012fr. We also analyze the earliest images of the SN, including pre-discovery observations obtained by the La Silla Quest survey less than one day after the explosion.

The host-galaxy of SN 2012fr, NGC 1365, is a member of the Fornax cluster and it has numerous published Tully–Fisher and Cepheid distance measurements. In the following, we adopt the Freedman et al. (2001) Cepheid distance modulus of μ = 31.27 ± 0.05 mag because of its internal consistency with the SNooPy analysis tool that we employ in this paper.25 The galactic reddening in the direction of NGC 1365 as computed from the Schlafly & Finkbeiner (2011) recalibration of the Schlegel et al. (1998) infrared dust maps is $E{(B-V)}_{\mathrm{MW}}=0.018\pm 0.003$ mag. When adopting a Fitzpatrick (1999) reddening law characterized by RV = 3.1, this corresponds to AV = 0.056 ± 0.009 mag. As discussed in Section 3.2, an analysis of the optical and NIR light curves indicates that SN 2012fr suffered little or no host-galaxy reddening.

The rest of this paper is organized as follows. Section 2 presents the observational data, consisting of the broadband ultraviolet (UV), optical, and NIR photometry. In Section 3, these observations are analyzed to derive the light-curve parameters, the host-galaxy reddening, an independent estimate of the distance based on the SN properties, a bolometric light curve, the rise time to maximum, and the mass of the 56Ni produced in the explosion. Section 4 discusses the early light curve and makes an assessment of the normality of SN 2012fr. Finally, the main conclusions of this study are summarized in Section 5.

2. Observations

2.1. Ultraviolet Photometry

The imaging of SN 2012fr was performed from space with the Swift UltraViolet Optical Telescope (UVOT). In this paper, we only present the uvw2, uvm2, and uvw1 passbands. According to Brown et al. (2016), the effective wavelengths of these filters (convolved with an SN Ia spectrum of SN 1992A, at +5 epoch) are 2010 Å, 2360 Å, and 2890 Å, respectively.

The photometry of SN 2012fr was computed from these images following the techniques described in detail by Brown et al. (2009, see their Section 2). In short, a 3'' to 5'' source aperture was used for each image , depending on the signal-to-noise ratio, to measure the counts at the position of the SN. The count rate from the underlying galaxy was then subtracted from these counts using images without the supernova present. The corrected counts were converted to the UVOT photometric system (Poole et al. 2008), while adopting the zero points provided by Breeveld et al. (2011). The measured instrumental magnitudes are corrected for a coincidence loss correction factor that is based on measurements made with a 5'' aperture. In addition, an aperture correction is applied based on an average point-spread function (PSF) in the Swift calibration database.

Covering ∼40 epochs ranging from −13 days to +120 days relative to the epoch of ${t}_{{B}_{\max }}$ (see Section 3), the high-cadence UV light curves of SN 2012fr are plotted in Figure 2. The corresponding photometry is listed in Table 1.

Figure 2.

Figure 2. UV, optical, and NIR light curves of SN 2012fr. For the Swift UV bands, the error bars are omitted for clarity. For the optical and NIR bands, the error bars are smaller in size than the symbols.

Standard image High-resolution image

Table 1.  UVOT Photometry of SN 2012fr in the Vega System

JD (days) uvw2 (mag) JD (days) uvm2 (mag) JD (days) uvw1 (mag) Phase (days)
2456229.32 19.139(186) 2456229.31 17.795(115) −13.78
2456231.99 16.798(077) 2456232.00 18.815(118) 2456231.99 15.090(058) −11.11
2456233.93 15.936(074) 2456233.94 17.676(101) 2456233.92 14.096(046) −9.17
2456235.68 15.625(074) 2456235.68 16.977(081) 2456235.71 13.559(044) −7.42
2456236.17 15.303(064) 2456236.16 13.502(042) −6.93
2456237.80 15.099(063) 2456237.80 16.603(070) 2456237.70 13.332(044) −5.30
2456239.17 16.372(072) −3.93
2456239.60 14.913(075) 2456239.60 16.327(108) 2456239.60 13.176(044) −3.50
2456241.64 14.878(068) 2456241.64 16.091(084) 2456241.64 13.165(043) −1.46
2456242.14 14.914(058) −0.96
2456243.13 14.980(063) 2456243.11 15.919(078) +0.03
2456244.90 14.958(054) +1.80
2456250.13 15.170(072) 2456250.14 15.950(073) 2456250.13 13.655(045) +7.03
2456251.70 15.351(075) 2456251.70 16.021(074) 2456251.69 13.839(047) +8.60
2456254.00 15.548(077) 2456254.00 16.242(077) 2456253.99 14.062(050) +10.90
2456262.36 16.540(088) 2456262.36 17.191(095) 2456262.35 15.089(067) +19.26
2456264.15 16.692(097) 2456264.15 17.296(111) 2456264.15 15.329(072) +21.05
2456268.23 17.077(103) 2456268.23 17.630(117) 2456268.22 15.673(075) +25.13
2456270.09 17.035(104) 2456270.10 17.685(124) 2456270.09 15.827(078) +26.99
2456280.18 17.578(110) 2456280.19 18.356(140) 2456280.18 16.822(091) +37.08
2456287.96 17.953(118) 2456287.96 18.434(137) 2456287.96 16.665(082) +44.86
2456296.02 18.132(131) 2456296.02 18.606(165) 2456296.01 16.899(090) +52.92
2456300.22 18.345(144) 2456300.23 18.726(166) 2456300.22 17.005(094) +57.12
2456303.96 18.207(139) 2456303.97 19.020(201) 2456303.96 17.109(099) +60.86
2456314.69 18.453(133) 2456314.69 19.082(177) 2456314.68 17.258(091) +71.59
2456332.85 18.727(166) 2456332.86 19.668(288) 2456332.85 17.690(116) +89.75
2456350.74 19.386(232) 2456350.75 19.815(306) 2456350.74 18.016(130) +107.64
2456356.82 19.010(178) 2456356.83 19.823(280) 2456356.82 17.963(121) +113.72

Note. The values in parentheses are 1σ measurement uncertainties and they are given in millimag. The phase is relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$.

Download table as:  ASCIITypeset image

2.2. Optical Photometry

2.2.1. CSP

Optical imaging of SN 2012fr was obtained with the 1 m Henrietta Swope telescope located at LCO, which was equipped with the same Johnson (BV) and Sloan (ugri) filter set and CCD detector that were used for the CSP-I (Contreras et al. 2010; Stritzinger et al. 2011; Krisciunas et al. 2017). All of the images were reduced in the manner described by Contreras et al. (2010) and Krisciunas et al. (2017), including the subtraction of the host-galaxy reference images that were obtained after the SN had fully disappeared. PSF photometry of the SN was computed with respect to a local sequence of stars that were calibrated to Landolt (1992) and Smith et al. (2002) standard fields, which were observed over the course of more than 60 photometric nights. The photometry of the local sequence in the standard system is given in Table 2.

Table 2.  Photometry of SN 2012fr Local Sequence Stars in the Standard System

Star ID R.A. Decl. u (mag) g (mag) r (mag) i (mag) B (mag) V (mag) Y (mag) J (mag) H (mag)
01 03:33:55.48 −36:07:19.0 16.836(010) 14.087(003) 12.846(002) 12.332(009) 14.748(003) 13.423(002)
02 03:33:33.67 −36:07:14.1 15.314(004) 13.879(002) 13.343(002) 13.160(008) 14.293(003) 13.560(002) 12.467(034) 12.167(023) 11.821(035)
03 03:33:49.55 −36:07:42.8 15.781(006) 14.774(003) 14.450(003) 14.370(008) 15.062(004) 14.569(003)
04 03:33:55.66 −36:08:38.2 16.881(007) 15.385(004) 14.833(003) 14.641(009) 15.803(006) 15.050(004)
05 03:33:40.67 −36:10:52.9 16.565(007) 15.368(004) 14.924(004) 14.789(009) 15.732(006) 15.091(005)
06 03:33:57.36 −36:09:24.0 18.750(074) 17.019(011) 16.429(007) 16.244(014) 17.463(014) 16.689(014)
07 03:33:54.96 −36:09:09.6 18.769(020) 17.017(015) 16.423(010) 16.237(014) 17.487(030) 16.659(013)
08 03:33:46.83 −36:03:04.1 15.247(007) 14.249(005) 13.916(005) 13.790(005)
12 03:33:16.73 −36:05:08.9 18.748(001) 16.033(012) 14.789(005) 14.267(010)
13 03:33:15.45 −36:03:47.4 18.102(007) 16.955(014) 16.555(007) 16.417(012)
102 03:33:31.07 −36:08:11.8 14.303(040) 13.933(173) 13.432(030)
103 03:33:45.26 −36:07:29.4 15.506(035) 15.011(033) 14.463(061)

Note. Values in parentheses are 1σ measurement uncertainties in millimag.

Download table as:  ASCIITypeset image

The definitive uBgVri-band photometry of SN 2012fr in the Swope natural system (Krisciunas et al. 2017) is given in Table 3 and the corresponding light curves are plotted in individual panels contained within Figure 2. The color curves are also plotted in Figure 3, showing the high quality of the data set. Comprising ∼120 epochs, the light curves track the flux evolution from −12 days to +160 days relative to ${t}_{{B}_{\max }}$, representing one of the most comprehensive data sets yet obtained of a Type Ia SN.

Figure 3.

Figure 3. Optical colors of SN 2012fr. The error bars are omitted for the sake of clarity.

Standard image High-resolution image

Table 3.  Optical Photometry of SN 2012fr in the Natural System of the Swope

JD (days) Phase (days) u (mag) g (mag) r (mag) i (mag) B (mag) V (mag)
2456230.70 −12.40 14.975(015) 13.952(009) 14.079(008) 14.576(009) 14.066(027) 13.990(015)
2456231.70 −11.40 14.220(017) 13.458(011) 13.609(010) 14.028(011) 13.524(025) 13.528(011)
2456232.70 −10.40 13.688(016) 13.098(010) 13.248(009) 13.646(009) 13.126(025) 13.178(009)
2456233.62 −9.48 13.351(018) 12.841(010) 13.022(009) 13.369(010) 12.846(026) 12.931(010)
2456234.66 −8.44 13.061(019) 12.778(012) 13.136(012) 12.621(026) 12.707(011)
2456235.60 −7.50 12.901(017) 12.450(011) 12.647(012) 13.026(013) 12.465(026) 12.569(013)
2456236.69 −6.41 12.758(017) 12.312(012) 12.515(012) 12.941(013) 12.333(028) 12.417(013)
2456238.72 −4.38 12.594(017) 12.298(010) 12.867(013) 12.113(030) 12.196(015)
2456239.67 −3.43 12.564(015) 12.058(009) 12.248(009) 12.843(010) 12.109(027) 12.128(011)
2456241.67 −1.43 12.559(015) 12.001(015) 12.149(009) 12.842(009) 12.021(025) 12.051(010)
2456242.64 −0.46 12.582(016) 11.994(010) 12.107(008) 12.851(010) 12.051(027) 12.016(011)
2456243.63 +0.53 12.621(019) 11.985(012) 12.084(012) 12.888(017) 12.018(030) 12.002(012)
2456244.65 +1.55 12.649(015) 11.988(011) 12.075(009) 12.913(011) 12.060(028) 11.994(010)
2456245.61 +2.51 12.687(019) 12.004(013) 12.061(015) 12.920(018) 12.057(030) 11.981(016)
2456246.63 +3.53 12.734(017) 12.010(011) 12.093(014) 12.946(014) 12.102(028) 12.000(013)
2456247.65 +4.55 12.776(018) 12.044(012) 12.096(011) 12.966(015) 12.140(029) 12.015(012)
2456248.62 +5.52 12.833(015) 12.072(011) 12.116(010) 13.043(013) 12.168(029) 12.035(010)
2456249.66 +6.56 12.890(016) 12.109(011) 12.148(009) 13.060(012) 12.237(029) 12.060(011)
2456250.67 +7.57 12.967(021) 12.153(014) 12.227(014) 13.094(016) 12.294(030) 12.102(017)
2456251.63 +8.53 13.053(019) 12.213(015) 12.253(014) 13.164(013) 12.342(033) 12.141(015)
2456252.68 +9.58 13.125(016) 12.269(009) 12.322(008) 13.214(010) 12.416(025) 12.193(010)
2456253.61 +10.51 13.222(016) 12.324(009) 12.387(008) 13.310(009) 12.480(025) 12.265(011)
2456254.64 +11.54 13.329(017) 12.370(011) 12.456(010) 13.378(011) 12.544(026) 12.311(010)
2456255.67 +12.57 13.434(018) 12.464(013) 12.554(013) 13.475(014) 12.646(026) 12.390(013)
2456256.60 +13.50 13.537(018) 12.528(013) 12.623(013) 13.571(014) 12.745(030) 12.433(015)
2456257.60 +14.50 13.655(018) 12.612(012) 12.682(013) 13.658(013) 12.835(028) 12.497(012)
2456258.62 +15.52 13.765(018) 12.692(011) 12.755(011) 13.712(012) 12.891(028) 12.582(013)
2456259.59 +16.49 13.884(018) 12.749(010) 12.808(009) 13.732(012) 13.000(030) 12.629(011)
2456261.62 +18.52 14.096(016) 12.890(009) 12.874(009) 13.747(010) 13.174(027) 12.716(010)
2456262.65 +19.55 14.240(016) 12.959(012) 12.901(012) 13.759(013) 13.268(027) 12.772(012)
2456263.65 +20.55 14.334(015) 13.022(009) 12.916(009) 13.714(011) 13.392(028) 12.791(011)
2456264.69 +21.59 14.445(016) 13.114(010) 12.920(009) 13.683(011) 13.464(027) 12.837(011)
2456265.63 +22.53 14.523(017) 13.171(011) 12.921(011) 13.648(013) 13.543(027) 12.885(012)
2456266.64 +23.54 14.630(016) 13.234(011) 12.922(010) 13.632(011) 13.642(028) 12.896(011)
2456267.65 +24.55 14.726(017) 13.308(012) 12.915(013) 13.579(013) 13.762(038) 12.935(013)
2456268.74 +25.64 14.816(015) 13.382(009) 12.926(008) 13.554(009) 13.817(027) 12.985(010)
2456269.78 +26.68 14.912(017) 13.451(010) 12.930(010) 13.518(013) 13.900(029) 13.028(013)
2456270.74 +27.64 14.985(015) 13.518(010) 12.946(010) 13.470(012) 13.982(027) 13.065(011)
2456271.74 +28.64 15.065(017) 13.590(012) 12.950(011) 13.429(013) 14.080(030) 13.116(013)
2456272.74 +29.64 15.125(015) 13.674(011) 12.983(010) 13.421(011) 14.105(029) 13.141(012)
2456273.71 +30.61 15.209(016) 13.731(011) 12.995(011) 13.411(012) 14.165(028) 13.186(012)
2456274.69 +31.59 15.284(016) 13.806(011) 13.034(010) 13.374(011) 14.258(028) 13.250(012)
2456275.72 +32.62 15.341(017) 13.880(011) 13.082(012) 13.375(014) 14.340(030) 13.305(014)
2456276.75 +33.65 15.409(017) 13.966(012) 13.141(010) 13.388(014) 14.401(030) 13.371(013)
2456278.61 +35.51 15.520(016) 14.104(011) 13.244(011) 13.477(012) 14.486(028) 13.515(011)
2456279.74 +36.64 15.581(017) 14.165(011) 13.334(010) 13.546(013) 14.545(036) 13.587(013)
2456280.73 +37.63 15.617(016) 14.229(010) 13.409(010) 13.633(010) 14.608(026) 13.667(012)
2456281.67 +38.57 15.648(015) 14.275(010) 13.456(009) 13.694(010) 14.646(025) 13.716(011)
2456282.76 +39.66 15.683(016) 14.330(010) 13.539(008) 13.789(009) 14.730(026) 13.784(009)
2456283.68 +40.58 15.736(016) 14.384(011) 13.595(010) 13.850(012) 14.747(026) 13.831(011)
2456284.72 +41.62 15.754(016) 14.423(009) 13.645(008) 13.910(010) 14.784(025) 13.897(009)
2456285.65 +42.55 15.806(016) 14.466(011) 13.698(012) 13.960(011) 14.828(027) 13.938(013)
2456286.65 +43.55 15.820(017) 14.492(011) 13.745(010) 14.009(012) 14.859(027) 13.973(012)
2456288.73 +45.63 15.860(017) 14.548(010) 13.816(009) 14.127(011) 14.908(026) 14.041(011)
2456289.60 +46.50 15.877(017) 14.574(012) 13.846(011) 14.158(013) 14.906(038) 14.079(012)
2456292.65 +49.55 15.950(016) 14.642(011) 13.960(010) 14.299(012) 14.963(027) 14.184(011)
2456293.69 +50.59 15.944(016) 14.656(012) 14.005(011) 14.348(015) 14.969(026) 14.205(011)
2456294.68 +51.58 15.958(015) 14.687(009) 14.040(009) 14.407(011) 14.999(026) 14.245(011)
2456295.66 +52.56 15.963(016) 14.692(009) 14.070(010) 14.435(012) 15.006(026) 14.266(011)
2456296.60 +53.50 15.988(016) 14.718(010) 14.088(011) 14.476(012) 15.016(025) 14.294(013)
2456297.67 +54.57 16.010(015) 14.738(009) 14.136(008) 14.518(008) 15.039(025) 14.319(009)
2456299.69 +56.59 16.041(016) 14.774(009) 14.207(009) 14.608(008) 15.066(024) 14.374(009)
2456305.60 +62.50 16.134(016) 14.881(010) 14.398(009) 14.840(010) 15.157(025) 14.547(011)
2456308.55 +65.45 16.188(017) 14.928(013) 14.496(014) 14.940(015) 15.205(027) 14.623(013)
2456311.61 +68.51 16.253(017) 14.982(011) 14.588(011) 15.059(014) 15.250(028) 14.700(012)
2456314.59 +71.49 16.289(016) 15.039(011) 14.682(010) 15.176(011) 15.303(027) 14.781(011)
2456315.63 +72.53 16.297(019) 15.057(010) 14.718(010) 15.239(021) 15.335(028) 14.815(012)
2456316.56 +73.46 16.351(017) 15.066(010) 14.743(010) 15.247(011) 15.329(028) 14.825(011)
2456317.58 +74.48 16.369(017) 15.085(012) 14.769(011) 15.265(014) 15.343(027) 14.850(012)
2456320.57 +77.47 16.420(018) 15.131(009) 14.866(009) 15.367(011) 15.408(028) 14.929(011)
2456321.56 +78.46 16.439(017) 15.152(010) 14.894(009) 15.418(011) 15.426(028) 14.947(011)
2456322.59 +79.49 16.470(016) 15.153(011) 14.920(011) 15.432(013) 15.450(028) 14.968(012)
2456323.62 +80.52 16.472(016) 15.178(010) 14.947(011) 15.480(014) 15.439(027) 14.986(013)
2456324.57 +81.47 16.492(016) 15.198(010) 14.997(009) 15.515(011) 15.474(027) 15.022(012)
2456325.56 +82.46 16.532(016) 15.214(011) 15.023(011) 15.531(014) 15.475(028) 15.033(012)
2456326.59 +83.49 16.540(016) 15.231(011) 15.057(010) 15.574(011) 15.512(028) 15.055(011)
2456327.56 +84.46 16.544(017) 15.255(010) 15.093(009) 15.626(011) 15.502(028) 15.083(011)
2456328.57 +85.47 16.588(017) 15.269(009) 15.117(009) 15.663(013) 15.534(027) 15.119(011)
2456329.56 +86.46 16.602(016) 15.279(011) 15.134(011) 15.681(014) 15.568(028) 15.120(013)
2456330.55 +87.45 16.620(017) 15.296(010) 15.188(010) 15.709(014) 15.616(030) 15.165(013)
2456331.53 +88.43 16.651(016) 15.313(010) 15.206(010) 15.596(027) 15.175(011)
2456332.52 +89.42 16.663(016) 15.324(010) 15.238(009) 15.779(013) 15.603(027) 15.192(011)
2456336.56 +93.46 16.758(016) 15.389(010) 15.353(010) 15.895(013) 15.663(028) 15.288(011)
2456337.57 +94.47 16.794(017) 15.420(009) 15.409(009) 15.941(013) 15.700(028) 15.320(011)
2456340.57 +97.47 16.840(016) 15.464(011) 15.466(012) 16.012(015) 15.725(028) 15.390(012)
2456343.57 +100.47 16.922(017) 15.507(009) 15.559(010) 16.099(013) 15.806(028) 15.432(011)
2456346.60 +103.50 17.012(027) 15.535(012) 15.657(011) 16.193(017) 15.839(029) 15.510(014)
2456351.58 +108.48 17.088(018) 15.627(010) 15.793(010) 16.349(013) 15.904(026) 15.593(011)
2456354.55 +111.45 17.163(016) 15.676(010) 15.897(011) 16.430(017) 15.975(029) 15.676(013)
2456357.55 +114.45 17.242(017) 15.736(009) 15.990(011) 16.543(014) 16.005(027) 15.731(011)
2456360.51 +117.41 17.300(016) 15.772(009) 16.058(011) 16.605(016) 16.066(027) 15.780(011)
2456364.53 +121.43 17.396(018) 15.838(010) 16.177(010) 16.740(015) 16.142(027) 15.850(011)
2456366.51 +123.41 17.449(017) 15.870(009) 16.227(011) 16.748(016) 16.166(027) 15.897(011)
2456368.51 +125.41 17.506(017) 15.894(009) 16.274(013) 16.799(017) 16.176(027) 15.932(011)
2456369.51 +126.41 17.500(017) 15.909(010) 16.305(013) 16.822(017)
2456370.55 +127.45 17.524(032) 15.919(009) 16.345(012) 16.881(015) 16.227(026) 15.978(010)
2456372.50 +129.40 17.551(018) 15.942(014) 16.395(012) 16.904(016) 16.237(026) 16.000(011)
2456376.50 +133.40 17.683(022) 16.010(009) 16.518(013) 17.021(016) 16.309(026) 16.081(010)
2456379.52 +136.42 17.724(025) 16.049(011) 16.575(013) 17.051(018) 16.354(027) 16.124(012)
2456382.49 +139.39 16.102(009) 16.666(012) 17.142(016) 16.412(026) 16.194(011)
2456385.49 +142.39 17.852(019) 16.141(009) 16.748(012) 17.206(016) 16.435(026) 16.237(011)
2456388.48 +145.38 17.933(021) 16.181(010) 16.830(013) 17.279(018) 16.480(030) 16.270(012)
2456398.48 +155.38 18.100(030) 16.336(015) 17.085(016) 17.474(020) 16.636(029) 16.484(018)
2456403.47 +160.37 16.392(012) 17.195(014) 17.573(019) 16.697(030) 16.542(013)
2456412.46 +169.36 16.524(013) 17.405(020)

Note. The phase is relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$. The values in parentheses are 1σ measurement uncertainties and they are given in millimag.

Download table as:  ASCIITypeset images: 1 2

As discussed in Contreras et al. (2010), photometry in the natural system is the purest form of the data and it provides the most transparent way to combine CSP photometry with data sets from other groups. Nevertheless, as requested by the referee, we provide S-corrections in Table 4 that, when added to the corresponding optical natural photometry magnitudes in Table 3, convert the photometry to the standard systems; i.e., Landolt (1992) for B and V , and Smith et al. (2002) for u, g, r, and i. The S-corrections were calculated using the Hsiao et al. (2007) spectral template and SN 2012fr spectra (Childress et al. 2013) when available (these values are given in parenthesis). The S-corrections are measured after the spectra are slightly altered to match the photometric colors of the corresponding phases.

Table 4.  S-corrections to Standard Photometric Systems

Filter B V u' g' r' i'
2456230.70 −0.040 −0.014 −0.014 −0.004 −0.000 −0.020
2456231.70 −0.036 −0.015 0.002 −0.004 −0.000 −0.016
2456232.70 −0.031 −0.022 0.015 −0.004 −0.000 −0.011
2456233.62 −0.031 −0.020 0.015 −0.005 −0.000 −0.012
2456234.66 −0.024 −0.014 0.030 −0.006 −0.000 −0.018
2456235.60 −0.026 −0.012 0.023 −0.006 −0.000 −0.017
2456236.69 −0.025 −0.013 0.011 −0.006 −0.001 −0.017
2456238.72 −0.020 −0.014(−0.008) −0.005 −0.005(−0.005) −0.000(−0.002) −0.017(−0.013)
2456239.67 −0.025 −0.015(−0.004) 0.007 −0.005(−0.005) −0.000(−0.002) −0.015(−0.013)
2456241.67 −0.026 −0.016 0.039 −0.005 0.000 −0.016
2456242.64 −0.025 −0.016(−0.016) 0.041 −0.006(−0.005) 0.000(−0.002) −0.018(−0.022)
2456243.63 −0.023 −0.017 0.051 −0.006 0.000 −0.020
2456244.65 −0.022 −0.018(−0.020) 0.053 −0.006(−0.006) 0.000(−0.002) −0.022(−0.025)
2456245.61 −0.020 −0.019(−0.011) 0.058 −0.006(−0.006) 0.000(−0.002) −0.024(−0.024)
2456246.63 −0.019 −0.020(−0.024) 0.066 −0.006(−0.006) 0.000(−0.002) −0.026(−0.032)
2456247.65 −0.018 −0.021(−0.014) 0.078 −0.006(−0.007) 0.001(−0.002) −0.026(−0.031)
2456248.62 −0.016 −0.022(−0.015) 0.089 −0.006(−0.007) 0.001(−0.002) −0.028(−0.032)
2456249.66 −0.014 −0.022(−0.015) 0.100 −0.006(−0.008) 0.001(−0.002) −0.029(−0.033)
2456250.67 −0.012 −0.023(−0.028) 0.112 −0.006(−0.007) 0.000(−0.002) −0.030(−0.038)
2456251.63 −0.011 −0.022(−0.029) 0.114 −0.006(−0.007) 0.000(−0.002) −0.030(−0.038)
2456252.68 −0.010 −0.021(−0.028) 0.128 −0.005(−0.008) 0.000(−0.002) −0.029(−0.038)
2456253.61 −0.008 −0.021 0.137 −0.005 0.000 −0.028
2456254.64 −0.005 −0.020 0.135 −0.004 0.000 −0.025
2456255.67 −0.005 −0.020 0.136 −0.002 0.000 −0.022
2456256.60 −0.008 −0.021 0.143 0.000 0.001 −0.019
2456257.60 −0.008 −0.021 0.149 0.002 0.001 −0.016
2456258.62 −0.003 −0.022 0.153 0.002 0.000 −0.014
2456259.59 −0.001 −0.023 0.157 0.002 0.000 −0.013
2456261.62 0.011(0.006) −0.025(−0.036) 0.165 0.004(−0.002) 0.001(−0.002) −0.012(−0.034)
2456262.65 0.013 −0.027 0.161 0.006 0.002 −0.013
2456263.65 0.014 −0.029 0.157 0.008 0.002 −0.013
2456264.69 0.015 −0.034 0.159 0.009 0.002 −0.011
2456265.63 0.014 −0.039 0.162 0.012 0.003 −0.008
2456266.64 0.013 −0.044 0.163 0.014 0.003 −0.005
2456267.65 0.010 −0.048 0.159 0.017 0.003 −0.000
2456268.74 0.007 −0.052 0.154 0.020 0.003 0.004
2456269.78 0.003 −0.055 0.148 0.023 0.003 0.008
2456270.74 0.002 −0.059 0.146 0.024 0.003 0.010
2456271.74 0.005 −0.061 0.144 0.025 0.003 0.012
2456272.74 0.010 −0.061 0.143 0.024 0.004 0.015
2456273.71 0.015 −0.060 0.142 0.022 0.004 0.017
2456274.69 0.018 −0.059 0.138 0.021 0.004 0.019
2456275.72 0.017 −0.060 0.129 0.022 0.004 0.021
2456276.75 0.017 −0.060 0.120 0.023 0.004 0.022
2456278.61 0.016(0.019) −0.060(−0.068) 0.109 0.024(0.021) 0.004(0.002) 0.024(0.002)
2456279.74 0.015 −0.058 0.100 0.024 0.004 0.025
2456280.73 0.015 −0.056 0.103 0.023 0.004 0.025
2456281.67 0.016 −0.054 0.107 0.022 0.004 0.025
2456282.76 0.017 −0.050 0.112 0.020 0.003 0.025
2456283.68 0.018 −0.049 0.117 0.018 0.003 0.025
2456284.72 0.016 −0.049 0.121 0.018 0.003 0.026
2456285.65 0.013 −0.050 0.125 0.018 0.003 0.028
2456286.65 0.010 −0.052 0.128 0.019 0.003 0.029
2456288.73 0.006 −0.055 0.131 0.021 0.003 0.030
2456289.60 0.006 −0.055 0.134 0.021 0.003 0.031
2456292.65 0.009 −0.052 0.140 0.020 0.003 0.033
2456293.69 0.010 −0.051 0.142 0.019 0.002 0.034
2456294.68 0.012 −0.050 0.144 0.018 0.002 0.035
2456295.66 0.012 −0.047 0.144 0.017 0.002 0.036
2456296.60 0.013 −0.048(−0.075) 0.144 0.018 0.002(−0.002) 0.036(0.015)
2456297.67 0.014 −0.047 0.144 0.018 0.002 0.037
2456299.69 0.016 −0.046 0.144 0.017 0.002 0.039
2456305.60 0.012 −0.040 0.145 0.013 0.001 0.039
2456308.55 0.010 −0.043 0.148 0.014 0.001 0.036
2456311.61 −0.003 −0.068 0.151 0.017 0.003 0.040
2456314.59 0.010 −0.039 0.153 0.014 −0.000 0.038
2456315.63 0.011 −0.043 0.151 0.014 0.000 0.038
2456316.56 0.012 −0.043 0.150 0.014 −0.000 0.037
2456317.58 0.013 −0.042 0.148 0.013 −0.000 0.037
2456320.57 0.017 −0.039 0.144 0.011 −0.001 0.035
2456321.56 0.018 −0.038 0.142 0.011 −0.001 0.034
2456322.59 0.018 −0.038 0.142 0.011 −0.001 0.034
2456323.62 0.019 −0.037 0.140 0.010 −0.001 0.033

Note. S-corrections, as given by ${M}_{\mathrm{std}}={M}_{\mathrm{nat}}+{S}_{M}$, are calculated for all phases using the Hsiao et al. (2007) spectral template. For those phases where a spectrum of SN 2012fr exists (Childress et al. 2013), the S-correction calculated from that spectrum is given in parentheses. For B and V, the filter functions and zero points are taken from Stritzinger et al. (2005). For u, g, r, and i, the BD+17°4708 spectrum that is given by Bohlin et al. (2014) is used, along with Fukugita et al.'s (1996) filter functions and Smith et al.'s (2002) standard magnitudes to obtain photometric zero points.

Download table as:  ASCIITypeset images: 1 2

2.2.2. TAROT

SN 2012fr was discovered in images obtained with the robotic telescope that is operated by the TAROT collaboration at the La Silla Observatory. The discovery image was taken on 2012 October 27.05 UT (JD 2456227.55) with an open filter (i.e., without a filter) and follow-up BVR images were obtained beginning two nights later. The telescope optics, filters, and detector are described in detail by Klotz et al. (2008). Throughput functions are presented in Appendix A for the open and V filters. The magnitudes for SN 2012fr and the local standards in the field are derived in the natural systems of these two filters.

2.2.3. La Silla-QUEST

In an attempt to better constrain the rise time of the SN, we examined the observing log of the La Silla-QUEST (LSQ) Low Redshift Supernova Survey (Baltay et al. 2013), which went into routine operations approximately one year before the discovery of SN 2012fr. The log indicated that the field of NGC 1365 was observed through a wide-band gr filter every night from 2012 October 23 to 27 UT. Subsequent examination of the images revealed that the SN was clearly visible on 2012 October 26.19 UT (see Figure 4) but was absent in the images obtained before this date. Unfortunately, the image of the SN was saturated in all of the images that were obtained after 2012 October 26 UT. The details of the telescope and detector system employed for the LSQ survey are given in Baltay et al. (2007, 2013). Appendix A derives the throughput function of the gr filter, and the natural system magnitudes for the SN and the local standards are presented.

Figure 4.

Figure 4. LSQ Survey images of NGC 1365 taken on 2012 October 25.34 UT and 2012 October 26.19 UT. The position of SN 2012fr is magnified in the upper left-hand corner of both images.

Standard image High-resolution image

2.2.4. Slooh

Klotz & Conseil (2012) reported that SN 2012fr was confirmed by E. Conseil from images obtained with the 0.5 m Slooh Space Camera robotic telescope at Mt. Teide, on the island of Tenerife. These observations were made with a set of Astrodon Tru-Balance LRGB E-Series filters using an FLI PL09000 CCD camera on 2012 October 27.13 UT; i.e., less than two hours after the TAROT discovery image. The Slooh observations are extremely important because they provide information on both the brightness of the SN and the color shortly after discovery. The throughput functions of the LRGB filters are estimated in Appendix A, and the natural system magnitudes are derived for the SN and the local standards in the field of SN 2012fr.

2.2.5. BOSS

Confirming images of SN 2012fr were taken by Stu Parker at Parkdale Observatory in New Zealand, starting 1.3 days after discovery by the TAROT collaboration (Parker 2012). Parker, who is a member of the Backyard Observatory Supernova Search (BOSS), used a 0.36 m Celestron telescope with an SBIG ST-10XME CCD camera and no filter for these observations. The throughput function of this combination is calculated in Appendix A, which we henceforth refer to as the BOSS open system, and the natural system magnitudes for the SN and the local standards are derived.

2.3. NIR Photometry

The majority of our NIR imaging was obtained on the LCO du Pont 2.5 m telescope with RetroCam, which is a Y JH-band imager employing a Hawaii-1 HgCdTe 1024 × 1024 pixel array. The field of view of RetroCam at the du Pont telescope is 3farcm× 3farcm5, with a pixel scale of 0farcs20. Additional NIR imaging of SN 2012fr was obtained in J1, J, and H filters using the FourStar camera attached to the 6.5 m Magellan Baade telescope (Persson et al. 2013). FourStar consists of a mosaic of four Hawaii-2RG HgCdTe detectors and each chip yields a field of view of 5' × 5'. The RetroCam and FourStar bandpasses are illustrated in Figure 5. It can be seen from this figure that the FourStar J1 filter covers ∼75% of the wavelength range of the RetroCam Y filter.

Figure 5.

Figure 5. Comparison between the transmission functions of the Swope + RetroCam YRC, JRC, and HRC filters and the Magellan Baade + FourStar J1FS, JFS, and HFS filters.

Standard image High-resolution image

The NIR images were reduced in the standard manner following the steps described by Contreras et al. (2010). In short, images were dark subtracted, flat fielded (and sky + fringing subtracted in the case of RetroCam images), and each dithered frame was aligned and combined. Host-galaxy reference images were subtracted from each combined image on a 15'' radius circle around the SN.

The NIR photometry of SN 2012fr was computed differentially with respect to a local sequence of stars, which were defined using RetroCam observations. The local sequence was calibrated in the Persson et al. (1998) JH photometric system using standard star fields that were observed during 10 photometric nights. For the Y-band, the local sequence stars were calibrated using the magnitudes for a subset of Persson et al. standards, as published by Krisciunas et al. (2017). The final NIR Y, J, and H magnitudes for the local sequence stars in these standard systems are listed in Table 2.

To compute magnitudes for the local sequence stars in the natural system of the three FourStar filters (hereafter, referred to as J1FS, JFS, and HFS), the following set of transformation relations were obtained via synthetic photometry of Castelli & Kurucz's (2003) model atmospheres:

Equation (1)

Equation (2)

Equation (3)

In these equations, Jstd and Hstd are the magnitudes in Persson et al.'s standard system; Ystd is the magnitude in the RetroCam Y-band standard system, as defined by Krisciunas et al. (2017); JFS, J1FS, and HFS are the magnitudes in natural system of FourStar; and ${(J-H)}_{\mathrm{std}}$ is the color index of the star in Persson et al.'s system. The Persson et al. and FourStar filter functions that are used to derive the synthetic magntiudes are made available on the CSP website.26

Making use of the local sequence, we now proceeded to compute the definitive NIR photometry of SN 2012fr in the natural system of the RetroCam Y JH filters (hereafter, referred to as YRC, JRC, and HRC, as in Krisciunas et al. 2017), which is listed in Table 5.

Table 5.  Near-infrared Photometry of SN 2012fr in the Natural System of du Pont

JD (days) Phase (days) Y (mag) J (mag) H (mag) Telescope
2456231.84 −11.26 13.675(027) 13.607(033) 13.712(031) BAA
2456232.87 −10.23 13.390(022) 13.302(020) 13.507(029) BAA
2456233.86 −9.24 13.177(022) 13.069(020) 13.269(029) BAA
2456234.71 −8.39 13.061(025) 12.953(020) 13.148(061) DUP
2456234.85 −8.25 13.028(022) 12.938(020) 13.207(029) BAA
2456235.73 −7.37 12.979(022) 12.848(020) 13.094(061) DUP
2456235.84 −7.26 12.974(022) 12.834(020) 13.072(031) BAA
2456236.67 −6.43 12.929(025) 12.799(020) 12.974(062) DUP
2456237.60 −5.50 12.894(025) 12.745(020) 12.990(061) DUP
2456241.86 −1.24 12.988(026) 12.776(025) 13.040(027) BAA
2456244.87 +1.77 13.224(022) 12.901(020) 13.134(027) BAA
2456246.86 +3.76 13.482(022) 13.048(020) 13.185(029) BAA
2456248.87 +5.77 13.624(022) 13.176(020) 13.235(028) BAA
2456249.87 +6.77 13.648(022) 13.286(020) 13.316(030) BAA
2456251.86 +8.76 13.780(022) 13.545(020) 13.398(040) BAA
2456252.84 +9.74 13.885(027) 13.639(020) 13.389(027) BAA
2456253.69 +10.59 13.861(022) 13.807(020) 13.383(027) DUP
2456254.68 +11.58 13.899(022) 13.984(020) 13.440(027) DUP
2456256.84 +13.74 13.949(022) 14.335(020) 13.450(028) BAA
2456258.59 +15.49 13.911(026) 14.521(020) 13.444(062) DUP
2456259.60 +16.50 13.875(025) 14.530(020) 13.398(061) DUP
2456264.59 +21.49 13.525(022) 14.459(020) 13.287(027) DUP
2456265.61 +22.51 13.459(026) 14.432(025) 13.254(030) DUP
2456266.61 +23.51 13.381(026) 14.393(025) 13.218(030) DUP
2456281.57 +38.47 12.660(022) 13.898(020) 13.291(027) BAA
2456284.60 +41.50 12.911(022) 14.163(020) 13.502(027) DUP
2456287.75 +44.65 13.094(022) 14.438(025) 13.698(028) DUP
2456293.61 +50.51 13.437(026) 14.871(025) 13.909(030) DUP
2456295.56 +52.46 13.522(022) 15.017(020) 13.972(027) DUP
2456296.54 +53.44 15.051(033) 14.039(027) BAA
2456296.62 +53.52 13.595(026) 15.092(025) 14.041(030) DUP
2456297.66 +54.56 13.655(026) 15.148(025) 14.085(030) DUP
2456298.53 +55.43 13.661(022) 15.131(020) 14.128(027) BAA
2456304.54 +61.44 13.967(022) 15.555(020) 14.350(029) BAA
2456318.59 +75.49 14.767(022) 16.388(021) 14.955(035) DUP
2456326.57 +83.47 15.170(022) 16.859(021) 15.280(029) DUP
2456327.56 +84.46 15.222(022) 16.900(021) 15.293(061) DUP
2456342.60 +99.50 15.861(025) 17.587(025) 15.926(063) DUP
2456343.59 +100.49 15.905(034) 17.565(028) 16.028(054) DUP
2456345.57 +102.47 16.003(026) 17.693(039) 16.092(033) DUP
2456346.55 +103.45 16.023(027) 17.709(033) 16.122(034) DUP
2456349.58 +106.48 16.138(026) 17.741(031) 16.134(064) DUP
2456352.58 +109.48 16.215(025) 17.871(022) 16.223(062) DUP
2456378.51 +135.41 16.911(022) 18.591(044) 16.897(045) DUP
2456382.53 +139.43 17.003(023) 18.723(036) 16.984(044) DUP

Note. The values in parentheses are 1σ measurement uncertainties and they are given in millimag. The phase is relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$.

Download table as:  ASCIITypeset image

This photometry is plotted in individual panels contained within Figure 2. Consisting of more than 40 epochs, the light curves follow the flux evolution from −11 to +140 days with respect to ${t}_{{B}_{\max }}$. The FourStar J1FS, JFS, and HFS photometry was S-corrected (Suntzeff 2000; Stritzinger et al. 2002) to match the RetroCam YRC, JRC, and HRC system using the following spectrophotometric prescription:

  • 1.  
    Spectral templates (Hsiao et al. 2007) matching the phases of each FourStar photometry epoch were employed. The template spectrum for each FourStar epoch was color matched separately to the photometry in the J1FS, JFS, and HFS bands, respectively.
  • 2.  
    For the J1FS band, the color-matching function was a second order polynomial calculated using CSP i and FourStar J1FS and JFS photometry. Likewise, for the JFS band, a second order polynomial was derived from the J1FS, JFS, and HFS photometry. For the HFS band, a linear polynomial was employed to color match the template to the JFS and HFS magnitudes.
  • 3.  
    The S-correction was then derived as the difference between the synthetic RetroCam and FourStar magnitudes derived from the mangled spectrum. Specifically: ${S}_{J{1}_{\mathrm{FS}}}\,={Y}_{\mathrm{RC}}-J{1}_{\mathrm{FS}};$ ${S}_{{J}_{\mathrm{FS}}}={J}_{\mathrm{RC}}-{J}_{\mathrm{FS}};$ and ${S}_{{H}_{\mathrm{FS}}}={H}_{\mathrm{RC}}-{H}_{\mathrm{FS}}$.

The resulting S-corrections are listed in Table 6. The excellent consistency between the S-corrected FourStar photometry and the RetroCam observations is illustrated in Figure 6. This is confirmed by looking at nights where nearly simultaneous measurements were made with both instruments. For the Y-band at phases of −8.3 and −7.3 days, we find differences between the S-corrected FourStar and RetroCam photometry of −0.022 ± 0.030 and −0.006 ± 0.031 mag, respectively. For the J band at phases of −8.3, −7.3, and +53.5 days, the differences are −0.016 ± 0.028, −0.017 ± 0.028, and +0.035 ±0.041 mag, respectively, while for the H band at the same phases, the differences are −0.061 ± 0.068, −0.020 ±0.060, and 0.000 ± 0.001 mag.

Figure 6.

Figure 6. Y JH light curves of SN 2012fr from −12 to +60 days relative to ${t}_{{B}_{\max }}$, acquired with the du Pont (+RetroCam) and the Magellan Baade (+FourStar) telescopes. The photometry obtained between the two facilities matches exceptionally well. In particular, note the excellent agreement between the FourStar S-corrected $J1\to Y$-band and RetroCam Y-band photometry.

Standard image High-resolution image

Table 6.  Baade FourStar S-corrections for SN 2012fr

Phase (days) ${S}_{J{1}_{\mathrm{FS}}}$ (mag) ${S}_{{J}_{\mathrm{FS}}}$ (mag) ${S}_{{H}_{\mathrm{FS}}}$ (mag)
−11 0.044 0.002 0.001
−10 0.032 0.001 0.002
−09 0.024 −0.001 0.002
−08 0.011 −0.001 0.002
−07 −0.001 −0.003 0.002
−01 −0.036 −0.003 0.000
+02 −0.012 −0.007 0.000
+04 −0.021 −0.010 0.001
+06 −0.016 −0.009 0.003
+07 −0.017 −0.007 0.005
+09 −0.024 −0.005 0.006
+10 −0.017 −0.007 0.006
+14 −0.027 −0.006 −0.013
+38 0.010 −0.004 0.001
+53 −0.006 −0.002
+55 0.063 −0.008 −0.001
+61 0.048 −0.013 −0.004

Note. The S-correction values were applied as: XRetrocam = XFourstar+S-corr. The phase was computed relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$ and rounded to match the phases of Hsiao et al. (2007) spectral templates.

Download table as:  ASCIITypeset image

3. Analysis

3.1. Light-Curve Parameters

The densely sampled light curves of SN 2012fr allow us to accurately measure the apparent magnitude at maximum and the decline-rate parameter, Δm15(X),27 in each filter, X, using a smooth Gaussian process's fitting curve. The values are summarized in Table 7. K-corrections were ignored as the small redshift (zhelio = 0.005457 according to NED28 ) of the host-galaxy, NGC 1365, has very little effect on these quantities. In addition, absolute magnitudes are provided as computed using the adopted average Cepheid distance discussed in Section 3.3 and assuming a host-galaxy dust reddening of $E{(B-V)}_{\mathrm{host}}\,=0.03\pm 0.03$ mag (see Section 3.2). The fit to the B-band light curve indicates that ${t}_{{B}_{\max }}$ occurred on JD 2456243.1 ± 0.3 and ${\rm{\Delta }}{m}_{15}(B)=0.82\pm 0.03$ mag. This rather slow decline-rate implies that SN 2012fr should be moderately over-luminous compared to a standard SN Ia with ${\rm{\Delta }}{m}_{15}(B)=1.1$ mag and an absolute B-band magnitude MB = −19.1 mag (Folatelli et al. 2010). Indeed, as shown in Table 7, SN 2012fr was ∼0.3 mag more luminous than this value.

Table 7.  SN 2012fr: Light-curve Parameters and Absolute Magnitudes

  ${t}_{\max }(\mathrm{JD})$ Δtpeak mXpeak ${m}_{\mathrm{Xpeak}}-{A}_{X}$ Δm15 ${M}_{X\mathrm{Ceph}}$
Band (days) (days) (mag) (mag) (mag) (mag)
u 2456240.5 ± 0.3 −2.8 12.59(02) 12.36(15) 0.87(03) −18.91(16)
B 2456243.1 ± 0.3 +0.0 12.04(02) 11.84(13) 0.82(03) −19.43(14)
g 2456243.5 ± 0.3 +0.6 11.96(02) 11.78(12) 0.70(03) −19.49(13)
V 2456245.1 ± 0.3 +2.0 11.99(02) 11.84(10) 0.67(03) −19.43(11)
r 2456245.7 ± 0.3 +2.6 12.06(02) 11.93(08) 0.79(03) −19.34(09)
i 2456240.1 ± 0.3 −3.0 12.84(02) 12.74(07) 0.59(03) −18.53(09)
Y 2456238.0 ± 1.0 −5.0 12.87(02) 12.81(04) 1.01(03) −18.46(06)
J 2456239.4 ± 1.0 −3.7 12.69(02) 12.65(03) 1.22(03) −18.62(06)
H 2456238.6 ± 1.0 −4.5 12.99(02) 12.96(02) 0.39(03) −18.31(05)
LUVOIR 2456242.7 ± 0.3 −0.4 ... ... ... ...
Δtpeak = tmax − tBmax
AX: Total reddening, i.e., Milky Way plus host-galaxy reddening.

Note. The values in parentheses are 1σ measurement uncertainties and they are given in hundredths of a magnitude.

Download table as:  ASCIITypeset image

As expected for a normal SN Ia, the NIR light curves reached primary maxima ∼4–5 days prior to ${t}_{{B}_{\max }}$. The absolute magnitudes in Y, J, and H were calculated using the Cepheid distance to NGC 1365 and they are given in Table 7. These values are also fully consistent with the average values for slow-to-mid-decliners given by Kattner et al. (2012) in Table 7.

3.2. Host-Galaxy Reddening

To determine if SN 2012fr suffered any reddening due to dust external to the Milky Way, we make use of three methods. In the first method, the (B − V) color curve of SN 2012fr is compared to the Lira relation (Lira 1995; Phillips et al. 1999). The Lira relation relies on the fact that the (B − V) color of normal unreddened SNe Ia follows a linear evolution with minimal scatter between 30 and 90 days past maximum (see Hoeflich et al. 2017 for a description of the physics underlying the Lira relation). Therefore, once corrected for galactic reddening, comparing the (B − V) color curve of any given SN Ia to the Lira relation provides an indication of the amount of dust reddening external to the Milky Way. Burns et al. (2014) provide the following fit to the Lira relation for a SN with ${\rm{\Delta }}{m}_{15}(B)$ = 0.82 mag based on the CSP-I data releases 1 and 2 (Contreras et al. 2010; Stritzinger et al. 2011):

Equation (4)

This calibration differs somewhat from an analysis done previously by Folatelli et al. (2010), who used a smaller sample of SNe Ia that was presumed to be unreddened. Meanwhile, Burns et al. (2014) used all of the SNe in the CSP-I sample with good photometric coverage extending beyond 40 days past ${t}_{{B}_{\max }}$. The late-time slopes of the $(B-V)$ light curves were measured, in addition to the value of $(B-V)$, at 45 days after ${t}_{{B}_{\max }}$ for each object separately. The median slope is used for the Lira law and the median absolute deviation is taken as its uncertainty. The observed distribution of $(B-V)$ colors at day +45 is modeled as the convolution of an intrinsic Gaussian distribution and exponential tail, which is similar to Jha et al. (2007) . The maximum and standard deviation of the resulting Gaussian distribution are taken as the Lira intercept and uncertainty, respectively.

Figure 7 illustrates a plot of the galactic reddening-corrected (B − V) color evolution of SN 2012fr from 25 to 95 days past ${t}_{{B}_{\max }}$. The Lira relation from Burns et al. (2014) is over-plotted as a solid line, as defined in Equation (4). Also shown for reference is the relation given by Folatelli et al. (2010). The comparison between the color evolution and the Lira relation from Burns et al. implies that SN 2012fr suffered minimal host-galaxy reddening, $E{(B-V)}_{\mathrm{host}}=0.03\pm 0.04$ mag. However, the excellent precision and sampling of the observations clearly reveal that the (B − V) temporal evolution of SN 2012fr was not linear at these epochs. From 35 to 60 days past ${t}_{{B}_{\max }}$, the slope is approximately −0.014 mag day−1, whereas from 60–95 days it is approximately −0.008 mag day−1. As shown by Burns et al. (2014) in their Figure 12, these values cover the range of slopes displayed by normal SNe Ia. Nevertheless, it is unusual to observe such a large change of slope in any single event.

Figure 7.

Figure 7. Galactic extinction-corrected (B − V) color evolution of SN 2012fr. The Lira relation as determined by Folatelli et al. (2010) is over-plotted as a dashed line while the solid line is the recalibration presented by Burns et al. (2014). In both cases, if the overall range of the Lira relation is considered, then the implication is that SN 2012fr suffered little or no host-galaxy reddening.

Standard image High-resolution image

The second method that we adopted to estimate the host-galaxy $E{(B-V)}_{\mathrm{host}}$ color excess relies on well-defined relations between maximum light intrinsic pseudo-colors6 29 and the decline-rate derived from a large sample of SNe Ia. This method is more often applicable to SNe Ia observations because it makes use of maximum light observations, which are normally more readily available than the post-maximum regime required for an accurate Lira relation analysis. This method is fully detailed in Burns et al. (2014) and we briefly describe it here.

The observed colors $(B-X)$, where X represents each filter except for B, are modeled as the sum of an intrinsic color that depends on the decline-rate parameter of the SN, the color excess $E{(B-X)}_{\mathrm{MW}}$ from the Milky Way dust, and the color excess $E{(B-X)}_{\mathrm{host}}$ from the host-galaxy ISM. The intrinsic colors are derived from an MCMC analysis of a training sample of SNe Ia from the CSP-I (Burns et al. 2014). The Milky Way component of the reddening is determined by assuming a value $E{(B-V)}_{\mathrm{MW}}$ from the Schlafly & Finkbeiner (2011) dust maps, a fixed value for the ratio of total-to-selective absorption ${R}_{V}^{\mathrm{MW}}=3.1$, and Fitzpatrick's (1999) reddening law. Finally, the host-galaxy component of the reddening is modeled with Fitzpatrick's (1999) reddening law and two free parameters: $E{(B-V)}_{\mathrm{host}}$ and ${R}_{V}^{\mathrm{host}}$, which are determined using MCMC methods. The resulting reddening $E{(B-V)}_{\mathrm{host}}\,=0.06\pm 0.02$ mag is quite low and, as a result, the posterior of ${R}_{V}^{\mathrm{host}}=4.2\pm 1.3$ more closely resembles the population distribution of RV from the training sample rather than the observed colors of SN 2012fr.

The Na i D interstellar absorption lines can be used to provide a third estimate of the host-galaxy reddening of SN 2012fr. For this purpose, we use the Keck HIRES echelle spectrum that was published by Childress et al. (2012). These authors measured weak absorption in the D1 and D2 lines due to gas in the Milky Way at a combined equivalent width (EW) of 118 ± 19 mÅ. Employing the fit to EW(Na i D) versus AV for Galactic stars using Munari & Zwitter's (1997) relation, as shown in Phillips et al. (2013) in their Figure 9, this value implies AV = 0.05 ± 0.03 mag. This is in excellent agreement with the Schlafly & Finkbeiner (2011) value given in Section 1. No Na iD absorption at the redshift of NGC 1365 is visible in the echelle spectrum of Childress et al. at a 3σ upper limit of EW(Na i D) = 82.8 mÅ, implying AV ≤ 0.04 ± 0.03 mag, or $E{(B-V)}_{\mathrm{host}}=0.01\pm 0.01$ mag if the interstellar gas has similar properties to the gas in the solar neighborhood.

Taken together, these three estimates are consistent with zero or negligible host-galaxy reddening of SN 2012fr. In the rest of this paper, we adopt a value of $E{(B-V)}_{\mathrm{host}}\,=0.03\pm 0.03$ mag, which is consistent with all three estimates. We also assume RV = 3.1, but the exact value is not critical since the reddening is low.

3.3. Distance to NGC 1365 as Derived from SN 2012fr

From the fits to the observed light curves of SN 2012fr using the SNooPy EBV method (Burns et al. 2011), the SN-based distance to NGC 1365 is found to be μ = 31.25 ±0.01stat ± 0.08syst mag. SNooPy also offers an alternative way to estimate the distance to the host-galaxy using the broadband light curves through the Tripp (1998) method. The functional form of the Tripp method relates the distance modulus of a SN Ia to its decline rate and color via:

Equation (5)

Here ${m}_{X}^{\max }$ is the observed K-corrected and galactic-reddening corrected magnitude at maximum, ${M}_{X}^{0}$ is the peak absolute magnitude of SNe Ia with ${\rm{\Delta }}{m}_{15}=1.1$ and zero dust extinction, ${b}_{X}$ is the slope of the luminosity versus decline-rate relation, ${\beta }_{X}^{{YZ}}$ is the slope of the luminosity-color relationship, and (Y − Z) is a pseudo-color at maximum. Note that the SNooPy parameter, ${\rm{\Delta }}{m}_{15}$, is the template-derived value of the decline-rate parameter, which correlates strongly with the directly measured value, ${\rm{\Delta }}{m}_{15}(B)$, but with some random and systematic deviations (see Figure 6 of Burns et al. 2011).

The Tripp method requires an accurate calibration between the relations of peak absolute magnitude versus ${\rm{\Delta }}{m}_{15}$ and pseudo-color. Here the calibrations presented by Folatelli et al. (2010) in Table 8, lines 2 and 6 are adopted, which are based on 26 well-observed SN Ia light curves published by Contreras et al. (2010) for the B-band, and 21 well-observed SN Ia light curves for J band of the same paper. The resulting estimates of distance modulus for NGC 1365 based on the Tripp method are μB = 31.14 ± 0.15 mag and μJ = 31.34 ± 0.14 mag.

Table 8.  grLSQ Light Curve of SN 2012fr

      Filter of
JD (days) Phasea (days) grLSQ (mag) Observation
2456225.72 −17.41 >20.38b grLSQ
2456225.84 −17.33 >20.34b grLSQ
2456226.69 −16.41 17.71 ± 0.03 grLSQ
2456226.77 −16.33 17.36 ± 0.02 grLSQ
2456227.55 −15.55 $16.03{\,}_{-0.10}^{+0.21}$ openTAROT
2456227.62 −15.48 16.24 ± 0.06 LSlooh
2456228.86 −14.24 15.21 ± 0.03 openBOSS
2456229.63 −13.47 14.66 ± 0.03 VTAROT
2456229.67 −13.43 14.56 ± 0.03 VTAROT
2456229.72 −13.38 14.56 ± 0.03 VTAROT
2456229.76 −13.34 14.50 ± 0.03 VTAROT
2456229.80 −13.30 14.52 ± 0.03 VTAROT
2456229.91 −13.20 14.41 ± 0.01 openBOSS
2456230.63 −12.47 13.99 ± 0.01 VTAROT
2456230.67 −12.43 13.94 ± 0.02 VTAROT
2456230.70 −12.40 13.95 ± 0.01 VCSP
2456230.72 −12.38 13.89 ± 0.01 VTAROT
2456230.76 −12.34 13.91 ± 0.03 VTAROT
2456230.80 −12.30 13.87 ± 0.01 VTAROT
2456230.94 −12.16 13.79 ± 0.02 openBOSS
2456231.75 −11.35 13.49 ± 0.01 VCSP
2456231.80 −11.30 13.45 ± 0.01 VTAROT
2456232.63 −10.47 13.18 ± 0.01 VTAROT
2456232.67 −10.43 13.13 ± 0.01 VTAROT
2456232.72 −10.38 13.12 ± 0.01 VTAROT
2456232.75 −10.35 13.14 ± 0.01 VCSP
2456232.76 −10.34 13.11 ± 0.01 VTAROT
2456233.70 −9.40 12.90 ± 0.01 VCSP
2456234.63 −8.47 12.69 ± 0.01 VTAROT
2456234.67 −8.43 12.66 ± 0.01 VTAROT
2456234.70 −8.40 12.68 ± 0.01 VCSP
2456234.72 −8.38 12.64 ± 0.01 VTAROT
2456234.76 −8.34 12.64 ± 0.01 VTAROT
2456235.67 −7.43 12.64 ± 0.01 VCSP
2456236.68 −6.42 12.39 ± 0.01 VCSP
2456238.72 −4.38 12.18 ± 0.01 VCSP
2456239.67 −3.43 12.12 ± 0.01 VCSP
2456241.73 −1.37 12.05 ± 0.01 VCSP
2456242.64 −0.46 12.02 ± 0.01 VCSP
2456243.62 +0.52 12.01 ± 0.01 VCSP
2456244.72 +1.62 12.00 ± 0.01 VCSP
2456245.61 +2.51 11.99 ± 0.01 VCSP
2456246.71 +3.61 12.01 ± 0.01 VCSP
2456247.65 +4.55 12.03 ± 0.01 VCSP
2456248.71 +5.61 12.06 ± 0.01 VCSP

Notes.

aThis phase was computed relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$. b3σ upper limit for nondetection.

Download table as:  ASCIITypeset image

These results are consistent with Freedman et al.'s (2001) Cepheid distance modulus of μ = 31.27 ± 0.05 mag, which was adopted in this work, indicating that SN 2012fr had a luminosity consistent with a normal SNe Ia with ${\rm{\Delta }}{m}_{15}(B)\,=0.82$ mag.

3.4. Bolometric Light Curve of SN 2012fr

The extended wavelength coverage of SN 2012fr afforded by observations spanning from UV through NIR wavelengths allows us to construct an essentially complete ultraviolet-optical-infrared bolometric light curve, which is usually termed as UVOIR bolometric curve in the literature. However, as pointed out by Brown et al. (2016), this does not clearly specifiy the boundaries of the wavelength domain over which the flux estimate is done. Consequently, in this paper we define the bolometric luminosity, LBol, as the luminosity between wavelengths of 1800 Å and infinity, corresponding to the sum of LUVOIR(3000–16600 Å) flux plus the contribution, Luvm2(1800–3000 Å), as deduced from the SWIFT photometry and the unobserved far-infrared ${L}_{\lambda \gt {\lambda }_{H}}(16600-\infty \,$Å).

Details of this calculation are given in Appendix B and the final bolometric light curve is plotted in Figure 8. The bolometric light curve derived by Pereira et al. (2013) for SN 2011fe is shown for comparison, which is one of the few other SNe Ia to have been well-observed in UV, optical, and NIR.

Figure 8.

Figure 8. Bolometric light curve of SN 2012fr compared to that of the normal SN Ia 2011fe (Pereira et al. 2013). The bolometric light curve shown for SN 2012fr is the average of the trapezoidal rule integration and spectral template fitting methods (see Appendix B). Top right: a Gaussian process smooth curve is fitted to the bolometric data to recover the time and amplitude of peak.

Standard image High-resolution image

At peak, SN 2012fr reached a maximum luminosity of LBol = (1.35 ± 0.14) × 1043 erg s−1, which is on the bright end of the normal SNe Ia distribution (cf. Figure 6 of Scalzo et al. 2014). The relative fractions of the NIR (λ > 10000 Å) integrated flux are plotted as a function of light-curve phase in Figure 9, showing that the NIR contribution to the bolometric light curve of SN 2012fr is nearly 15% at −12 days. It then falls to a minimum of ∼4% a few days after ${t}_{{B}_{\max }}$ and rises again steeply to 19% at +40 days, at which point it begins to slowly decrease again. The latter behavior is similar to that shown by Scalzo et al. (2014) in Figure 3 and it is responsible for the prominent shoulder in the light curve that can be seen between +20 and +45 days in Figure 8. Figure 9 shows that the UV contribution (λ < λu) to the bolometric light curve of SN 2012fr is generally of less importance than the NIR. The only exception is around the ${t}_{{B}_{\max }}$ epoch, when both contributions are similar, peak around 10% a few days before ${t}_{{B}_{\max }}$ and contribute >5% of the integrated flux only during the early epochs; i.e., before 20 days after ${t}_{{B}_{\max }}$. This behavior is consistent with previous attempts to quantify the UV contribution to the bolometric luminosity (e.g., see Suntzeff 2003).

Figure 9.

Figure 9. Estimated ratio of UV and NIR fluxes to the total flux for SN 2012fr. The dotted and dashed lines represent the ratios for estimated fluxes beyond the observed optical and NIR domain; i.e., fluxes for λ less than 3000 Å and for λ greater than the effective wavelength of the H band. For comparison purposes, the red solid line gives the flux beyond 10000 Å and the blue solid line shows the UV flux for wavelengths less than u-band effective wavelength.

Standard image High-resolution image

3.5. Rise Time

SN 2012fr was discovered by the TAROT collaboration on 2012 October 27.05 (UT), 15.6 days before ${t}_{{B}_{\max }}$. As reported by Klotz (2012), the SN was not visible to a limiting magnitude of R > 15.8 in an image of NGC 1365 that was taken three days earlier on 2012 October 24.05 (UT) with the same telescope. A more stringent non-detection of R > 19.3 on 2012 October 24.02 (UT) was obtained by J. Normand from stacked images taken with an 0.6 m telescope at the Observatoire des Makes (Klotz 2012). Thus, the rise time to ${t}_{{B}_{\max }}$ was constrained to somewhere between 18.6 and 15.6 days.

Fortunately, the LSQ images that are presented in this paper provide a much tighter constraint on the rise time because the SN was clearly visible on 2012 October 26.19 (UT) but was absent in an image of similar depth that was obtained on 2012 October 25.34 (see Figure 4). Thus, the SN was detected less than a day after the explosion, which occurred some time between 17.3 and 16.5 days before ${t}_{{B}_{\max }}$.

Our observations of SN 2012fr present a rare opportunity to study the early rising light curve of SN 2012fr. However, we first must S-correct the various measurements to the same filter bandpass. Given that the earliest detection and non-detection of SN 2012fr were made in the LSQ gr filter, we choose to convert the Slooh, TAROT, BOSS, and CSP photometry to LSQ gr magnitudes. The steps that are required are:

  • 1.  
    Slooh L filter. As discussed in Appendix A, the Slooh L-band photometry is essentially in the same natural system as the grLSQ observations. Hence, no S-correction is required for this measurement, which was obtained less than two hours after the TAROT discovery image.
  • 2.  
    TAROT open filter. The TAROT discovery image was obtained within approximately 1–2 days of explosion. This is more than a day before the first spectroscopic observation. Fortunately, the Slooh B and G filter observations, which were obtained less than two hours after the TAROT discovery image, provide color information that can be used to estimate the S-correction under the assumption that the spectrum at this epoch can be approximated by a blackbody. Figure 10 shows the magnitude difference $({{gr}}_{\mathrm{LSQ}}-{\mathrm{open}}_{\mathrm{TAROT}})$ as a function of ${(B-G)}_{\mathrm{Slooh}}$ as derived from synthetic photometry of main-sequence stars from the Pickles (1998) stellar library spectra and black bodies of varying temperature. From the photometry given in Appendix A, ${(B-G)}_{\mathrm{Slooh}}\,=0.28\pm 0.31$ mag. Using the blackbody curve, this implies $({{gr}}_{\mathrm{LSQ}}-{\mathrm{open}}_{\mathrm{TAROT}})={0.09}_{-0.08}^{+0.20}$ mag. As a sanity check on this result, we plot the $(B-V)$ color evolution for SN 2012fr in Figure 11. For comparison, observations of SN 2011fe are also shown. Combining the BSlooh and GSlooh magnitudes measured for the SN and the color–color plots in Appendix A gives $(B-V)=0.40\pm 0.31$ mag. This measurement is plotted in Figure 11 and it appears to be generally consistent with expectations if SN 2011fe is a valid comparison. Nevertheless, it should be kept in mind that the correction of the TAROT photometry to the LSQ system that was derived in the previous paragraph is strictly only valid for an object with a stellar or blackbody spectrum. Strong features such as the Ca ii and Si ii HVFs observed in the earliest spectra of SN 2012fr could affect the S-correction.
  • 3.  
    BOSS open filter. The first BOSS open filter observation was made 14.2 days before ${t}_{{B}_{\max }}$. This is only ∼0.3 days after the first spectrum of the SN obtained by Childress et al. (2012) and, therefore, we have used this spectrum (published by Childress et al. 2013) to calculate the S-correction. We find $({{gr}}_{\mathrm{LSQ}}-{\mathrm{open}}_{\mathrm{BOSS}})\,=-0.02\pm 0.01$ mag, where the error reflects the uncertainty in the spectrophotometric calibration of the spectrum. S-corrections for the BOSS open filter observations at −12.2 and −7.1 days were obtained from synthetic photometry of the Childress et al. (2013) spectra after color matching the spectra to the CSP photometry using first- or second- order polynomials (see Figure 12). Finally, the S-correction for the BOSS open filter observation at −13.2 days was obtained by interpolating the S-corrections for the −14.2 and −12.2 day spectra.
  • 4.  
    CSP and TAROT V filters. Although the V filter is narrower than the LSQ gr filter, they are well-matched in central wavelength. The CSP V observations began 12.4 days before ${t}_{{B}_{\max }}$ and the TAROT V filter imaging started one night earlier at −13.5 days. S-corrections for both filters derived from the color-matched Childress et al. (2013) spectra are plotted in Figure 12.

Figure 10.

Figure 10. Color–color plot for the conversion of TAROT open filter magnitudes to the natural system magnitudes for the LSQ gr filter using the ${(B-G)}_{\mathrm{Slooh}}$ measurement obtained less than two hours after the TAROT observation. The black curve show synthetic photometry carried out using the Pickles (1998) stellar atlas, and the red line corresponds to synthetic photometry of black bodies covering a range of temperature. The solid vertical line and the dashed lines on either side indicate the ${(B-G)}_{\mathrm{Slooh}}$ color measurement of 0.28 ± 0.31 mag.

Standard image High-resolution image
Figure 11.

Figure 11. Early-time (B − V) and $(g-r)$ color evolution of SN 2012fr. The color evolution of SN 2011fe is shown for comparison. The observations for SN 2012fr are taken from the CSP, Tarot, Swift, and LCOGT (Graham et al. 2017). The observations of SN 2011fe are from Swift (Brown et al. 2012) and synthetic photometry in the CSP filter bandpasses of the spectrophotometry of Pereira et al. (2013). The abscissa is corrected for time dilation.

Standard image High-resolution image
Figure 12.

Figure 12. Magnitude differences between synthetic photometry of the Childress et al. (2013) spectra of SN 2012fr in the CSP V, TAROT V, and BOSS open filters and synthetic photometry of the same spectra in the LSQ gr filter. The spectra were first color matched to the CSP BVgri photometry to improve their spectrophotometric precision. The single BOSS point at −14.5 days was derived from the first spectrum obtained of the SN at −14.5 days (Childress et al. 2012). Likewise, the CSP V and TAROT V points at −14.5 and −14.2 days were calculated from the first two spectra obtained by (Childress et al. 2013). The larger errors on these points reflect the fact that these spectra could not be color matched because the CSP photometry did not begin until −12.4 days.

Standard image High-resolution image

The resulting LSQ gr light curve is given in Table 8 and it is plotted in normalized flux units in Figure 13. In addition, Figure 13 indicates the LSQ non-detections and the epoch of the first spectrum. From the first detection of the SN at −16.4 days with respect to ${t}_{{B}_{\max }}$ to the first BOSS observation at −14.2 days, the light curve rises very close to linearly. After the first BOSS observation, the light curve rises more steeply to a second nearly linear phase that lasts from approximately −11.5 to −6.5 days. The non-detection at −17.33 days and the first detection at −16.41 days constrain the time of explosion to have occurred at JD 2456226.23±0.46, which is 16.87 ± 0.46 days before ${t}_{{B}_{\max }}$. The bolometric maximum was reached on JD 2456242.7±0.3 (see Table 7). Hence, the rise to bolometric maximum took a total of ∼16.47 ± 0.55 days.

Figure 13.

Figure 13. (Top) Relative fluxes (normalized to maximum) plotted as a function of the time with respect to ${t}_{{B}_{\max }}$ during the rising phase of the grLSQ-band light curve of SN 2012fr and extending to a few days past maximum. (Bottom) Enlargement of the same observations during the seven days following explosion. The error bars are not visible because they are the same size or smaller than the symbols that are used to plot the data. The broken power-law fits of the early light curves of SNe 2013dy and 2014J (Zheng et al. 2013, 2014) are plotted for comparison.

Standard image High-resolution image

3.6. 56Ni Mass

With a well-sampled bolometric light curve and a precise measurement of the bolometric rise time at hand, the amount of 56Ni synthesized during the explosion can be estimated using Arnett's rule (Arnett 1982), which relates the bolometric rise time, tr, and peak bolometric luminosity, Lpeak, to the energy deposition, ENi, within the expanding ejecta supplied by the radioactive decay chain ${}^{56}\mathrm{Ni}\,\to {}^{56}\mathrm{Co}\,\to {}^{56}\mathrm{Fe}$ (see, e.g., Stritzinger & Leibundgut 2005). This is formulated as follows:

Equation (6)

Arnett's rule is derived from semi-analytical solutions to the radiative transfer problem of the expanding SN Ia ejecta. It explicitly assumes equality between energy generation and luminosity; i.e., the factor α = 1. A number of years ago, Branch (1992) surveyed the explosion models that were then available in the literature and concluded that α = 1.2 ± 0.2 was a more appropriate value to use. While this value is still commonly employed for determining 56Ni masses (see Scalzo et al. 2014, and the references therein), Höflich & Khokhlov (1996) found that α ranged from 0.7 to 1.4 for a large variety of explosion models, with an average value of 1.0 ± 0.2. Meanwhile, Stritzinger & Leibundgut (2005) cited radiative transport calculations for two modern 3D deflagration models from the MPA group that were both consistent with α = 1.0.

An expression for ${E}_{\mathrm{Ni}}({t}_{r})$ is provided by Nadyozhin (1994, see his Equation (18)). After plugging in various constants, this yields the following relation for 1 M of 56Ni:

Equation (7)

By combining Equations (6) and (7), we obtain the following simple relation to estimate the 56Ni mass:

Equation (8)

Given the nonlinearity of Equation (8), we computed the error in MNi by simulating 105 computations using randomly drawn uncorrelated values for α, $E{(B-V)}_{\mathrm{host}}$, tr, the distance modulus, and the peak bolometric flux, Fpeak, assuming that these parameters are described by Gaussian distributions with mean and standard deviation vectors of μ = [1.0, 0.03, 16.47, 31.27, Fpeak] and σ = [0.2, 0.03, 0.60, 0.05, 0.05Fpeak], respectively. Arnett's parameter, α, the host-galaxy reddening, and the distance modulus are the dominant error sources, while the other parameters only exert a mild effect. The error in the host-galaxy reddening alone translates to an error of 10% in MNi. The uncertainty in α has an even larger effect and introduces a slight asymmetry to the marginalized distribution of the 56Ni mass. The final value derived from this analysis is ${M}_{\mathrm{Ni}}={0.60}_{-0.14}^{+0.16}\,{M}_{\odot }$.

Childress et al. (2015) gave an independent estimate of MNi for SN 2012fr based on measurements of the [Co iii] λ5893 emission in nebular-phase spectra. Adjusting their value for the distance modulus for NGC 1365, which is adopted in the present paper, and assuming the same reddening of $E{(B-V)}_{\mathrm{host}}\,=0.03\pm 0.03$ mag gives a value of MNi = 0.61 ± 0.07 M. This is fully consistent with our estimate from the bolometric light curve.

We note that Zhang et al. (2014) derived a much higher 56Ni mass from a bolometric light curve that was constructed from their own optical photometry, the same SWIFT ultraviolet observations used in the present paper, and NIR corrections taken from SN 2005cf (Wang et al. 2009). Adjusting their quoted value of 0.88 ± 0.08 M to the same distance modulus and host-galaxy reddening that we assume in this paper gives a value of 0.84 ± 0.08 M. From J. Zhang (private communication), the peak bolometric luminosity of Zhang et al. (2014) was corrected to Lpeak = 1.65 × 1043 erg s−1. The mismatch is mainly due to an overestimation of the NIR contribution. After this correction and putting both measurements at the same distance modulus, the difference amounts to 7%, half of which can be explained for our differences in the u-band domain flux.

4. Discussion

4.1. Early Light Curve

The earliest emission of a SN Ia probes the location in the ejecta of the 56Ni that powers the light curve and, therefore, is an important diagnostic of the explosion's physics (Piro & Nakar 2013). Early studies of the rise times of SNe Ia (Riess et al. 1999; Aldering et al. 2000; Goldhaber et al. 2001; Conley et al. 2006; Hayden et al. 2010; Ganeshalingam et al. 2011; González-Gaitán et al. 2012) assumed or concluded that at the earliest epochs the flux increase is proportional to t2. This so-called fireball model is physically motivated by the idea that the luminosity of a young SN Ia in a homologous expansion is most sensitive to its changing radius, and is less sensitive to changes in the temperature and photospheric velocity. More recently, Firth et al. (2015) also invoked a model with a single power law but they allowed the exponent to be a free parameter. In a variation of this, Shappee et al. (2016) also employed a power-law model but they allowed the exponent to be different for each filter.

Only recently have observations of individual SNe Ia allowed stringent tests of the validity of the power-law model. In two events, SNe 2011fe (Nugent et al. 2011) and 2012ht (Yamanaka et al. 2014), the fireball model (i.e., n = 2) seems to provide a good fit to the earliest measurements. The most convincing of these two cases is SN 2011fe, which was discovered in M101 by the Palomar Transit Factory (PTF) with a g-band magnitude of 17.35 (Nugent et al. 2011). PTF images that were obtained in the previous night showed no source at a limiting magnitude of g ≥ 21.5. Nugent et al. (2011) found that the fireball model provided a consistent fit to the first three nights of g-band measurements and they used it to infer a date of explosion only 11 hr before discovery. However, as shown by Piro & Nakar (2014), the bolometric light curve of SN 2011fe seems not to follow the fireball model.

Olling et al. (2015) presented observations of the light curves of three SN Ia that were followed nearly continuously by Kepler. The rising portion of the light curve of the brightest of these SNe, KSN 2011b, was found to be well-fitted by a single power-law but had an exponent of n = 2.44 ± 0.14.

Nevertheless, observations of SNe 2013dy (Zheng et al. 2013) and 2014J (Zheng et al. 2014; Goobar et al. 2015) have demonstrated quite clearly that the single power-law model does not apply to all SNe Ia. For both of these objects, high-cadence (essentially daily) observations both before and after explosion revealed that the early-time light curves were well-described by a varying power-law exponent. The flux rose nearly linearly during the first day and then transitioned over the next 2–4 days to a relation closer to the t2 law. Based on the results for these two events, Zheng et al. (2014) speculated that the varying power-law behavior may be common to SNe Ia and that previous results favoring the t2 law may have been due to a lack of high-cadence observations constraining the shape of the light curve at the earliest epochs. This conclusion is supported by recently published observations of the early light curve of iPTF16abc (Miller et al. 2018), which also display a nearly linear rise during the first three days following explosion.

Figure 13 shows a comparison between our observations of the rising light curve of SN 2012fr with the broken power-law model used by Zheng et al. (2013, 2014) to fit the light curves of SNe 2013dy and 2014J. The time of first light for these fits has been adjusted to coincide with the value that we have determined for SN 2012fr. Note that for both of the latter SNe, the fits were made to unfiltered photometry and, therefore, some caution should be taken in comparing these results with our grLSQ light curve of SN 2012fr. Nevertheless, the resemblance is remarkable and supports the view that this behavior of a nearly linear rise initially, followed by a steeper increase in luminosity, may be common among SNe Ia.

Piro & Morozova (2016) have investigated how the distribution of 56Ni in the ejecta of a SN Ia affects the earliest phases of the light curve using SNEC (SuperNova Explosion Code; Morozova et al. 2015), which is an open source Lagrangian radiation hydrodynamics code that allowed them to initiate a shock wave within a white dwarf model, explode the white dwarf, and follow the early light curve evolution. These authors found that models with more highly mixed 56Ni rise more quickly than do models with centrally concentrated 56Ni. Using a grid of 800 models generated with SNEC, we have attempted to match the early light curve of SN 2012fr. A comparison of one of our better fitting light curves along with a range of other light curves with varying 56Ni mixing is shown in Figure 14 (with the best-fit light curve in green). This demonstrates that the steepening of SN 2012fr's light curve over the first couple of days can naturally be accounted for by a moderately mixed 56Ni distribution. This corresponds to a 56Ni mass fraction of 0.05 at roughly 0.05M below the surface of the exploding white dwarf. If 56Ni were not mixed out to this shallow region, then the theoretical light curves would tend to rise too slowly in comparison to SN 2012fr. Whereas, if the 56Ni mixing were more strong, then we would not expect to see the change in slope in the early light curve.

Figure 14.

Figure 14. Observed flux in the LSQ gr-band as a function of time during the first seven days following explosion for a bare white dwarf model. The line colors indicate the level of 56Ni mixing, which correspond to the profiles shown in Figure 15.

Standard image High-resolution image

This fit suggests that the first data point for this event was ∼21 hr after explosion. Constraining the explosion time like this is important if we wish to be able to put limits on the interaction with a companion, as argued in Shappee et al. (2016). The corresponding 56Ni mass fraction distributions for the models from Figure 14 are shown in Figure 15.

Figure 15.

Figure 15. Profiles of the mass fraction 56Ni as a function of the mass coordinate in the white dwarf for the various levels of mixing corresponding to the light curves shown in Figure 14. These are constrained by the observations only out to the vertical dashed line.

Standard image High-resolution image

The observations only constrain the distribution of 56Ni for material interior to the vertical dashed line. Constraints on shallower material would require even earlier observations. Whether or not the 56Ni distribution we find for SN 2012fr is unique is unclear because degeneracies may arise between various physical parameters when only photometric data is considered (Noebauer et al. 2017). This point will be explored in more detail in a forthcoming work (A. L. Piro et al. 2017, in preparation). In the future, it will be useful to also have spectral information at these early epochs to uniquely determine the physical cause of the early light curve shape.

4.2. Spectroscopic Peculiarities

As more SNe Ia have been observed, several subclasses within the general phenomenon have been identified. Benetti et al. (2005) divided SNe Ia into three groups based on the light curve decline rate, ${\rm{\Delta }}{m}_{15}(B)$, and the post-maximum evolution of the velocity of the minimum of the Si ii λ6355 absorption. The FAINT group, as represented by the prototypical SN 1991bg (Filippenko et al. 1992; Leibundgut et al. 1993; Turatto et al. 1996), consists of fast-declining events (${\rm{\Delta }}{m}_{15}(B)\gt 1.5$ mag) with low Si ii velocities at maximum, which decrease rapidly with time. SNe Ia is then further split with normal decline rates (${\rm{\Delta }}{m}_{15}(B)\lt 1.5$ mag) into two further categories: the HVG group, which display a high temporal velocity gradient ($\dot{v}\gt 70$ $\mathrm{km}\,{{\rm{s}}}^{-1}$ day−1) in the days following ${t}_{{B}_{\max }}$; and the more common LVG events, which show a lower velocity gradient. The HVG events typically also have high Si ii velocities at maximum. This led Wang et al. (2009) to propose a parallel subtype classification that was solely based on a measurement of the Si ii λ6355 velocity within a week of maximum. Wang et al. termed those having velocities ≳11,800 $\mathrm{km}\,{{\rm{s}}}^{-1}$ at the epoch of ${t}_{{B}_{\max }}$ high-velocity (HV) SNe Ia. Lower-velocity events, not including the 1991bg-like and 1991T-like spectral types (Branch et al. 1993), were termed "normal" by Wang et al. As shown in Figure 2 of Foley et al. (2011) and Figure 6 of Silverman et al. (2012), ∼80%–90% of HV events also pertain to the HVG group. Conversely, most LVG objects have normal Si ii velocities in the Wang et al. classification system (Silverman et al. 2012).30

Branch et al. (2006, 2009) developed an independent classification system for SNe Ia that is based on a plot of the pseudo-equivalent widths of the Si ii λ5972 and λ6355 absorption features. Four groups were established: Core Normal (CN), Cool (CL), Broad Line (BL), and Shallow Silicon (SS). Figure 16 displays these groups in a branch diagram, using the data and classification criteria of Blondin et al. (2012). There is a rough correspondence between the BL and HV (or HVG) groups, with Blondin et al. (2012) and Folatelli et al. (2013) independently finding that two-thirds of SNe Ia that are classified as BL also belong to the HV subtype. The CL-type correlates well with the Benetti et al. FAINT and Wang et al. 1991bg groups, and the LVG (or Wang et al. normal) SNe generally encompass the CN and SS classes. Finally, the SS group includes not only 1991T-like SNe but also the peculiar 2002cx-like events (Li et al. 2003).

Figure 16.

Figure 16. Plot of the pseudo-equivalent width of the Si ii λ5972 absorption vs. that of the Si ii λ6355 absorption within five days from maximum light, illustrating the four subtypes of SNe Ia defined by Branch et al. (2006). Core Normal (CN; violet squares), Shallow Silicon (SS; green upward-looking triangles), Broad Line (BL; blue downward-looking triangles), and Cool (CL; red diamonds) subtypes taken from the data of Blondin et al. (2012) are plotted. The positions of the shallow velocity gradient SNe 2000cx, 2006is, 2009ig, and 2012fr are indicated by black circles.

Standard image High-resolution image

So, how does SN 2012fr fit into these classification schemes? Childress et al. (2013) found that SN 2012fr lies on the border separating the Branch et al. CN and SS subclasses (see Figure 16). The shallow SN 2012fr Si ii velocity gradient displayed by SN 2012fr also places it clearly in the Benetti et al. LVG group. Nevertheless, the high Si ii velocity at maximum of ∼12,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$ qualifies SN 2012fr as an HV event in the Wang et al. system. This is an unusual combination of classifications because 5% or less of SNe Ia in the CN + SS groups belong to the HV subtype and <10% of LVG SNe Ia are also classified as HV events (Blondin et al. 2012; Folatelli et al. 2013).

Figure 17 shows our measurements of the evolution of the expansion velocities of the Si ii λ6355, Si iii λλ4564,5740, S ii λλ5449,5622, Ca ii λ8662, and Fe ii λλ4924,5018, absorption minima from the Childress et al. (2013) spectra of SN 2012fr. The expansion velocities of the Mg ii λ10927 line measured from unpublished CSP-II NIR spectra of SN 2012fr are also included. For comparison, we plot the expansion velocities of the same features as measured from the spectra of the prototypical SN 2011fe published by Pereira et al. (2013) and Hsiao et al. (2013). The small amount of velocity stratification of the intermediate mass elements (IMEs) in SN 2012fr from −5 days onward is remarkable and is in stark contrast to that observed for SN 2011fe. Childress et al. (2013) concluded that there is either a shell-like density enhancement in the ejecta at a velocity of ∼12,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$ or a sharp cutoff in the radial distribution of the IMEs in the ejecta. As discussed by Quimby et al. (2007) in the context of the slow-declining SN 2005hj—which is a Wang et al. (2009) normal event that showed a nearly flat Si ii velocity gradient—a strong density enhancement is not predicted by standard delayed-detonation and deflagration models but is instead suggestive of rapidly expanding material interacting with overlying material. Examples of scenarios that produce a well-defined shell of IMEs include pulsating delayed-detonations (Khokhlov et al. 1993) or tamped detonations in a double-degenerate merger (Fryer et al. 2010). This possibility will be explored in more detail in a future paper (C. Cain et al. 2018, in preparation).

Figure 17.

Figure 17. Comparison of the velocity stratification of different ions for SNe 2011fe and 2012fr. The measurements correspond to the absorption minima of each line.

Standard image High-resolution image

In the case of SN 2012fr, Childress et al. (2013) argued that it is more likely that the small velocity evolution of the Si ii and Ca ii lines reflects the fact that these ions are physically confined to a narrow region in velocity space in the ejecta. This conclusion was based, in part, on the observation that the expansion velocity of the Fe ii lines in SN 2012fr began to slowly decrease ∼10 days after maximum while the Si ii and Ca ii lines remained at a nearly constant velocity of 12,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$ (see Figure 17). One interpretation of Figure 17 is that the inner edge of the IMEs in the ejecta of SN 2012fr was located at a velocity of ∼11,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$, which is in contrast to a typical SN Ia as represented by 2011fe for which an abundance tomography analysis indicates that the inner edge of the IMEs extended down to ∼5000 $\mathrm{km}\,{{\rm{s}}}^{-1}$ (Mazzali et al. 2015). This suggests a 56Ni distribution that extended out to higher velocities than normal. In other words, at the phase range where the photosphere of a typical SN Ia is still in the silicon-rich layer, the photosphere of SN 2012fr may already have receded into the 56Ni-rich layer due to its extended distribution. Höflich et al. (2002) presented one-dimensional, delayed-detonation models of a Chandrasekhar mass white dwarf that reproduce the luminosity–decline rate relation, from the sub-luminous to the most luminous SNe Ia, by varying the density at which the deflagration transitions to a detonation, ${\rho }_{\mathrm{tr}}$, from values of 8–27 ×  106 gm cm−3. The same model with a slightly higher transition density of ${\rho }_{\mathrm{tr}}=30\times {10}^{6}$ gm cm−3 produces a minimum velocity of Si/S of 12,983 $\mathrm{km}\,{{\rm{s}}}^{-1}$, 56Ni mass = 0.67 M, and MV = −19.35 mag (from private communication with Peter Höflich).

Maund et al. (2013), Childress et al. (2013), and Zhang et al. (2014) have pointed out that SN 2012fr shares certain properties with luminous events such as SN 1991T (Filippenko et al. 1992; Phillips et al. 1992; Ruiz-Lapuente et al. 1992), SN 1999aa (Li et al. 2001; Garavini et al. 2004), and the previously mentioned SN 2005hj: slow-declining light curves, weak Si ii λ6355 absorption, and shallow photospheric velocity gradients. From this coincidence and the strong presence of HVFs during the premaximum phases, Zhang et al. (2014) proposed that SN 2012fr may represent a subset of the 1991T-like SNe Ia viewed at an angle where the ejecta has a clumpy or shell-like structure. However, SN 2012fr differs from 1991T-like events in three important ways, as follows: the high photospheric velocity of the Si ii λ6355 absorption at maximum light, the strong Si ii absorption observed at early epochs, and the lack of strong features due to Fe iii at maximum.

4.3. A Distinct Subclass?

Figure 18 illustrates the unusual nature of the Si ii λ6355 photospheric velocity evolution of SN 2012fr that began around 15 days before ${t}_{{B}_{\max }}$. The 1σ dispersion about the average of the velocities for normal SNe Ia in the Wang et al. system (Folatelli et al. 2013) is shown for comparison, while the dashed lines correspond to a subset of the family of functions derived by Foley et al. (2011) to describe the velocity evolution of LVG and HVG SNe Ia. The extraordinarily shallow velocity evolution of SN 2012fr is unmatched, except for a handful of other HV events. The evolution of the Si ii velocity for the two best observed of these—SN 2006is (Folatelli et al. 2013) and SN 2009ig (Marion et al. 2013)—is shown for comparison in Figure 18. The positions of these two SNe in the branch diagram are also indicated in Figure 16.

Figure 18.

Figure 18. Evolution of the photospheric Si ii λ6355 velocity for the shallow velocity gradient SNe 2000cx, 2006is, 2009ig, and 2012fr. The shaded area shows the average and 1σ dispersion of Si ii velocities of normal SNe Ia reproduced from Folatelli et al. (2013), while the dashed lines represent a subset of the family of functions that describe the velocity evolution of normal and high-velocity gradient SNe Ia (Foley et al. 2011). The measurements for SNe 2000cx, 2009ig, and 2012fr are taken from Li et al. (2001), Marion et al. (2013), and Childress et al. (2013), while those for SN 2006is were made using the spectra of Folatelli et al. (2013).

Standard image High-resolution image

SN 2009ig was extensively observed by Foley et al. (2012) and Marion et al. (2013). This slow-declining SN (${\rm{\Delta }}{m}_{15}(B)=0.89$ mag) was caught very early and it displayed strong HVF absorption in Ca ii and Si ii in the first spectra obtained ∼14 days before ${t}_{{B}_{\max }}$. The earliest spectra were dominated by the HVFs. The first detection of photospheric features began 12 days before ${t}_{{B}_{\max }}$. Both the HVF and photospheric absorption remained visible until 6 days before ${t}_{{B}_{\max }}$, when the HVF absorption was no longer detectable. The optical spectral evolution of the SN 2009ig was remarkably similar to that of SN 2012fr (Childress et al. 2013), with the only significant differences being: (1) the higher velocity, broader HVF absorption in the earliest (day −14) spectrum of SN 2012fr; (2) the longer persistence of the HVF Ca ii and Si ii absorption in SN 2012fr; and (3) the narrower Ca ii, Fe ii, S ii, and Si ii photospheric absorption in SN 2012fr from maximum onward (seen most clearly in the Ca ii triplet). The photometric evolution of both SNe was also striking similar, with the most significant difference being in the I band, where the secondary maximum of SN 2012fr occurs several days later than that of SN 2009ig. Like SN 2012fr, SN 2009ig was discovered <1 day after explosion (Foley et al. 2012). The R-band light curve of SN 2009ig is compared with the grLSQ light curve of SN 2012fr in Figure 19, from which it is clear that the rise time to ${t}_{{B}_{\max }}$ of SN 2009ig was somewhat longer than that of SN 2012fr. The unfiltered observation corresponds to the discovery image of SN 2009ig. Unfortunately, the last non-detection was four days before discovery, making it difficult to determine with certainty that the early rise of SN 2009ig showed the same broken power-law morphology that was observed for SN 2012fr. However, we note that the shape of R-band light curve of SN 2009ig closely mimics that of the grLSQ light curve of SN 2012fr for much of the rise to maximum, which suggests a similar morphology at the earliest epochs following explosion.

Figure 19.

Figure 19. Relative fluxes (normalized to maximum) plotted as a function of the time with respect to ${t}_{{B}_{\max }}$ for the R-band light curve of SN 2009ig and the grLSQ light curve of SN 2012fr. The photometry for SN 2009ig is taken from (Foley et al. 2012). The unfiltered first observation was approximately S-corrected assuming a color (B − V) = 0.6 ± 0.1 and a sensitivity function similar to those of the Tarot and BOSS unfiltered photometry. The black arrow shows the last non-detection before discovery for SN 2009ig and the red arrow shows the same for SN 2012fr. The abscissa is corrected for time dilation.

Standard image High-resolution image

SN 2006is, has been discussed by Folatelli et al. (2013). Photometrically, it was also quite a slow decliner, with ${\rm{\Delta }}{m}_{15}(B)=0.80$ mag (Stritzinger et al. 2011). Unfortunately, spectroscopic observations did not begin until maximum, and so it is unknown whether this SN also displayed strong HVF absorption at early epochs. Nevertheless, as shown in Figure 14 of Folatelli et al. (2013), the maximum-light spectra of SNe 2006is and 2009ig were remarkably similar.

From this above discussion, we conclude that SNe 2006is, 2009ig, and 2012fr were similar events sharing the following characteristics:

  • 1.  
    All three were slow decliners (${\rm{\Delta }}{m}_{15}(B)=0.8\mbox{--}0.9$ mag).
  • 2.  
    All three occupied a similar region of the branch diagram near the edge of the CN and SS distributions (see Figure 16).
  • 3.  
    All three displayed unusually shallow Si ii velocity gradients that are consistent with LVG events in the Benetti et al. (2005) classification scheme but at velocities ≳12,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$, which place them in the Wang et al. (2009) HV class.

Moreover, SNe 2009ig and 2012fr displayed remarkably strong HVFs that persisted to maximum light. Unfortunately, SN 2006is was not discovered early enough to determine if it also shared this characteristic.

Very few other SNe Ia have displayed this particular combination of properties. The most notable exception is the peculiar SN 2000cx that was extensively observed by Li et al. (2001) and Candia et al. (2003). In addition to the high, nearly constant evolution of the Si ii velocities (see Figure 18), SN 2000cx displayed remarkably strong Ca ii HVFs at velocities of >20,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$ that persisted to maximum light (Branch et al. 2004; Thomas et al. 2004). Interestingly, the evolution of the expansion velocities of the IMEs as deduced from the Si ii λ6355 and S ii λλ5449,5622 lines (Li et al. 2001) closely resembles the evolution that was observed for SN 2012fr. Optical and NIR photometry of SN 2000cx also revealed certain anomalies in its photometric behavior, as follows: (1) an asymmetric B-band light curve, with a relatively fast rise from −10 days to maximum but then a slow post-maximum decline (${\rm{\Delta }}{m}_{15}(B)=0.93$ mag); (2) a weaker and earlier-occurring I-band secondary maximum than would be expected for such a slow decline rate; and (3) a $(B-V)$ color evolution displaying several peculiarities, including a brief plateau phase of nearly constant color that began ∼1 week after maximum and a strikingly bluer color than predicted by the Lira relation from ∼35 to 90 days after maximum.

It is tempting to postulate that SNe 2006is, 2009ig, and 2012fr were closely related to SN 2000cx and to its apparent twin, SN 2013bh (Silverman et al. 2013). The difference is that SN 2000cx was a fairly extreme SS event in the Branch et al. (2006, 2009) classification scheme (see Figure 16). Interestingly, SN 2012fr and, to a lesser extent, SN 2009ig also resembled SN 2000cx in displaying a B-band light curve that rose rather quickly to maximum but then declined more slowly than average after maximum. This is illustrated in Figure 20. Here, the normalized B light curves are plotted with respect to ${t}_{{B}_{\max }}$, with time dilation taken into account. The Goldhaber et al. (2001) B-band Parab-18 template is also plotted, although it has been stretched to match the observed decline rates of each SN. The unusually rapid rise to maximum of SN 2000cx stands out clearly in this figure and is mimicked to a large extent by SN 2012fr. SN 2009ig also initially appeared to rise somewhat more quickly than the Goldhaber et al. template, although this difference would not have been so obvious if this event had not been discovered so early. Li et al. (2001) also called attention to a peculiarity in the $(B-V)$ evolution of SN 2000cx, which showed a phase of nearly constant color at $(B-V)\sim 0.3$ mag between 6 and 15 days past ${t}_{{B}_{\max }}$. As seen in Figure 21, the $(B-V)$ evolution of SN 2012fr showed a similar, nearly constant color of $(B-V)\sim 0.2$ mag for a few days centered around day +10. Unfortunately, the photometry of SN 2009ig is insufficiently precise to discern whether it also displayed such a peculiarity. However, Figure 21 shows that a change in the slope of the $(B-V)$ evolution of this SN occurred around day +10, as it did for SN 2012fr.

Figure 20.

Figure 20. Comparison of B-band light curves for SN 2000cx (${\rm{\Delta }}{m}_{15}(B)=0.93$ mag), SN 2009ig (${\rm{\Delta }}{m}_{15}(B)=0.89$ mag), and SN 2012fr (${\rm{\Delta }}{m}_{15}(B)=0.82$ mag). The abscissa is corrected for time dilation. Goldhaber et al.'s (2001) B-band template is shown for comparison, stretched to the decline rates of the three SNe.

Standard image High-resolution image
Figure 21.

Figure 21. Comparison of the $(B-V)$ color evolution of SN 2000cx (green squares), SN 2009ig (black diamonds), and SN 2012fr (red circles). The abscissa is corrected for time dilation. The data for SN 2000cx are from Li et al. (2001), for SN 2009ig are from Foley et al. (2012), and for SN 2012fr are from this paper. No correction has been made for host-galaxy reddening for any of the SNes. The dashed line is the average Lira relation from Burns et al. (2014).

Standard image High-resolution image

5. Conclusions

We have presented densely sampled, high quality, six-band optical photometry of the Type Ia SN 2012fr in the Fornax Cluster member NGC 1365. The data span epochs from 13 days before to 140 days after the epoch of ${t}_{{B}_{\max }}$, with typical errors below 2%. We also present similarly high quality NIR and UV photometric data sets. Based on these observations, we conclude the following:

  • 1.  
    SN 2012fr was a slow declining (${\rm{\Delta }}{m}_{15}(B)\,=0.82\pm 0.03$ mag), luminous event. From the observed colors at maximum, the evolution of the $(B-V)$ color at later epochs, and the lack of detectable host-galaxy interstellar Na i D absorption in high-dispersion spectra, we estimate that SN 2012fr suffered little or no host-galaxy reddening, adopting a conservative value of $E{(B-V)}_{\mathrm{host}}\,=0.03\pm 0.03$ mag.
  • 2.  
    Analysis of the optical and NIR light curves shows that the luminosity of SN 2012fr was completely normal for its decline rate.
  • 3.  
    The images obtained by the LSQ survey tightly constrain the epoch of the explosion to 16.87 ± 0.46 days before ${t}_{{B}_{\max }}$ or 16.5 ± 0.5 days before the bolometric maximum. The luminosity of the SN increased nearly linearly at first, transitioning to a faster rising phase ∼2.5 days after explosion. This behavior is well-fitted by an explosion model with moderate mixing of 56Ni in the ejecta.
  • 4.  
    The densely sampled bolometric light curve derived from our UV, optical, and NIR photometry indicates that SN 2012fr peaked at a luminosity of LBol =1.35 ± 0.14 × 1043 erg s−1. By combining this value with the measured rise time, we estimate a 56Ni mass of 0.60 ± 0.15 M using Arnett's rule. This amount is consistent with an independent measurement of the 56Ni mass from nebular spectra.
  • 5.  
    Despite its normal luminosity, SN 2012fr displayed spectroscopic properties that set it apart from most other SNe Ia. The classification as a Branch et al. (2006) SS/CN event, with HV Si ii absorption in the Wang et al. (2009) system but with a shallow velocity gradient (LVG using the Benetti et al. 2005 criteria) is uncommon and reminiscent of the peculiar SN 2000cx. Like SN 2000cx, SN 2012fr also displayed a fast rise to maximum and a much slower post-maximum decline. Finally, the velocity of the inner edge of the IMEs in SN 2012fr appears to have been ∼11,000 $\mathrm{km}\,{{\rm{s}}}^{-1}$, or approximately double that of a typical SN Ia such as SN 2011fe. We call attention to two other SNe Ia, 2006is and 2009ig, that showed photometric and spectroscopic properties similar to those of SN 2012fr and we suggest that all three, along with SN 2000cx, may form a distinct subclass of SNe Ia.
  • 6.  
    The Cepheid-based distance to its host-galaxy is already available, which adds SN 2012fr to the still relatively small number of nearby SNe Ia that are suitable for measuring the value of the Hubble constant, H.

We acknowledge Wei Kang Zheng for his help in providing discovery images of SN 2009ig and for his useful discussion. The work of the CSP has been supported by the National Science Foundation under grants AST0306969, AST0607438, AST1008343, AST1613426, AST1613455, and AST1613472. M.S. and C.C. acknowledge the generous support from the Danish Agency for Science and Technology and Innovation through a Sapere Aude Level 2 grant. E.Y.H. is supported by Florida Space Grant Consortium. M.S. is also supported by a research grant (13261) from VILLUM FONDEN.

Appendix A: Throughput Functions and Zero Points for the LSQ, Tarot, Slooh, and BOSS Filters

A.1. La Silla-QUEST

The relative throughput of the LSQ gr filter is plotted in panel (a) of Figure 22. This was constructed by multiplying the filter and CCD quantum efficiency curves given by Baltay et al. (2007, 2013) with a reflectivity curve for aluminum and an atmospheric transmission spectrum appropriate for the La Silla Observatory. The transmission of the corrector of the ESO Schmidt telescope was assumed to be flat over the spectral region covered by the filter. LSQ images were processed by the LSQ pipeline before being given to us.

Figure 22.

Figure 22. Relative throughput functions of the LSQ, Tarot, Slooh, and BOSS filters.

Standard image High-resolution image

We measured the instrumental PSF photometry of SN 2012fr and field stars in the LSQ images using DAOPhot routines (Stetson 1987). The final magnitudes for the local sequence stars in the natural system of the LSQ gr filter were calculated using the following procedure:

  • 1.  
    Figure 23 shows color–color diagrams of $({V}_{\mathrm{CSP}}-{{gr}}_{\mathrm{LSQ}})$ and $({r}_{\mathrm{CSP}}-{{gr}}_{\mathrm{LSQ}})$ colors versus ${(g-r)}_{\mathrm{CSP}}$ derived from synthetic photometry of main-sequence stars,31 selected from the Pickles (1998) stellar library spectra. The zero point of the synthetic photometry using the LSQ gr filter throughput function was set by requiring grLSQ = 0.0 for Vega (α Lyr).32 Likewise, synthetic magnitudes in the CSP V, g, and r filters were calculated using the throughput functions and zero points given by Krisciunas et al. (2017).
  • 2.  
    The photometry of the local sequence stars in the field of SN 2012fr are plotted in Figure 23 as red circles. The zero point for the instrumental magnitudes of the local sequence stars was found by shifting the photometric measurements along the y-axis to match the sequence defined by the Pickles stars.

Figure 23.

Figure 23. Color–color plots to convert CSP V and r magnitudes to the natural system magnitudes for the LSQ gr filter. The black curves show synthetic photometry carried out using the Pickles (1998) stellar atlas. The red points correspond to observations of the local sequence stars in the field of SN 2012fr.

Standard image High-resolution image

The resulting magnitudes for the local sequence stars in the natural system of the LSQ gr filter are given in Table 9. Our measurement of the magnitude of SN 2012fr on 2012 October 26 UT in the natural system of the LSQ gr filter is given in Table 10.

Table 9.  LSQ, TAROT, BOSS, and Slooh Natural System Photometry of SN 2012fr Local Sequence Stars

      grLSQ openTAROT VTAROT openBOSS RSlooh GSlooh BSlooh LSlooh
Star ID R.A. Decl. (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag)
1 03:33:55.48 −36:07:19.0 12.748(007) 13.301(002) 12.844(034) 12.484(014) 14.403(028) 13.591(016) 13.242(006)
2 03:33:33.67 −36:07:14.1 13.293(007) 13.323(014) 13.134(023) 14.127(023) 13.659(020) 13.511(008)
3 03:33:49.55 −36:07:42.8 14.513(005) 14.376(011) 14.541(003) 14.386(026) 14.130(040) 14.975(042) 14.539(033) 14.534(010)
4 03:33:55.66 −36:08:38.2 14.994(005) 14.824(015) 15.006(005) 14.789(035) 14.399(048) 15.721(077) 15.222(074) 15.001(015)
5 03:33:40.67 −36:10:52.9 15.041(006) 14.822(015) 15.004(015) 14.863(025) 14.644(052) 15.507(060) 15.131(052) 15.074(017)
6 03:33:57.36 −36:09:24.0 16.654(012) 16.505(053) 16.608(028) 16.421(067) 16.620(288) 17.028(228) 16.909(365) 16.600(052)
7 03:33:54.96 −36:09:09.6 16.641(011) 16.501(058) 16.652(031) 16.407(059) 15.988(163) 17.751(516) 16.883(342) 16.685(054)
8 03:33:46.83 −36:03:04.1 13.848(018) 13.676(026) 14.380(137) 14.051(027)
9 03:33:30.90 −36:02:08.1 15.169(019) 15.805(008) 15.329(039) 14.731(050) 16.886(202) 16.311(179) 15.697(024)
10 03:33:39.87 −36:01:27.2 16.244(042) 16.564(015) 16.331(041) 15.903(179) 17.800(565) 16.677(278) 16.422(045)
11 03:33:10.58 −36:03:41.9 13.334(064) 13.072(031) 13.806(103) 13.487(017)
12 03:33:16.73 −36:05:08.9 15.198(005) 14.651(014) 15.238(005) 14.788(045) 14.274(051) 16.331(117) 15.437(078) 15.144(016)
13 03:33:15.45 −36:03:47.4 16.640(011) 16.428(048) 16.641(016) 16.477(095) 15.831(139) 17.407(319) 16.374(183) 16.524(054)
14 03:33:58.35 −36:02:18.5 16.447(056) 16.518(014) 16.401(067) 15.854(138) 16.936(242) 16.509(195) 16.478(049)
15 03:34:01.65 −36:08:45.9 15.095(015) 15.306(005) 15.100(046) 15.012(098) 15.741(073) 15.404(075) 15.286(017)
16 03:33:28.37 −36:14:38.0 15.573(022) 15.980(010) 15.670(037) 15.280(084) 16.740(166) 16.030(132)
17 03:33:15.26 −36:14:39.3 15.313(070) 14.890(052) 16.319(113) 15.967(124)
18 03:33:02.16 −36:13:28.2 14.588(013) 14.892(004) 14.661(110) 14.498(047) 15.590(055) 15.172(062) 14.929(012)
19 03:33:12.50 −36:11:10.3 16.137(010) 15.857(031) 16.077(010) 15.934(068) 15.662(112) 16.641(152) 16.563(222) 16.080(033)
20 03:33:52.46 −36:14:47.6 15.993(046) 15.674(114) 16.374(120) 16.358(188)
21 03:33:04.96 −36:02:30.6 16.184(039) 16.898(022) 16.367(133) 15.798(141) 18.654(331) 17.229(430) 16.801(066)
22 03:33:14.32 −36:12:13.6 16.394(048) 17.017(033) 16.505(052) 15.892(174) 18.372(820) 17.872(669) 16.976(073)
23 03:33:12.02 −36:11:49.2 15.905(035) 16.774(018) 16.128(056) 15.612(098) 17.460(446) 16.777(238) 16.623(051)

Note. The values in parentheses are 1σ measurement uncertainties and they are given in millimag.

Download table as:  ASCIITypeset image

Table 10.  LSQ, TAROT, BOSS, and Slooh Natural System Photometry of SN 2012fr

JD Phase grLSQ openTAROT VTAROT openBOSS RSlooh GSlooh BSlooh LSlooh
(days) (days) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag)
2456226.69 −16.41 17.705(029)
2456226.77 −16.33 17.364(021)
2456227.55 −15.55 15.940(061)
2456227.62 −15.48 15.663(139) 16.112(149) 16.426(274) 16.238(0.064)
2456228.86 −14.24 15.214(029)
2456229.63 −13.47 14.667(027)
2456229.67 −13.43 14.573(021)
2456229.72 −13.38 14.576(018)
2456229.76 −13.34 14.520(020)
2456229.80 −13.30 14.541(022)
2456229.91 −13.20 14.408 (011)
2456230.63 −12.47 14.059(014)
2456230.67 −12.43 14.004(015)
2456230.72 −12.38 13.953(012)
2456230.76 −12.34 13.976(026)
2456230.80 −12.30 13.937(012)
2456230.94 −12.16 13.791 (015)
2456231.80 −11.30 13.513(008)
2456232.63 −10.47 13.238(005)
2456232.67 −10.43 13.188(005)
2456232.72 −10.38 13.175(006)
2456232.76 −10.34 13.162(005)
2456234.63 −8.47 12.733(003)
2456234.67 −8.43 12.705(003)
2456234.72 −8.40 12.688(003)
2456234.76 −8.34 12.684(003)
2456234.80 −8.30 12.688(003)

Note. The phase was computed relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$. The values in parentheses are 1σ measurement uncertainties and they are given in millimag.

Download table as:  ASCIITypeset image

SN 2012fr was not detected in a pair of grLSQ images acquired on 2012 October 25.34 UT (see Figure 4). We performed aperture photometry at the location of SN 2012fr on these images using the IRAF apphot package and we calibrated the results using the local sequence stars. This procedure implies 3-σ limits of grLSQ > 20.38 mag and grLSQ > 20.34, respectively. This limit was verified by placing artificial sources with increasing magnitudes at the location of SN 2012fr in these images and then carrying out photometry on the resulting images. The artificial source was recovered when it was m = 20.34 mag, which is in excellent agreement with the previously given measured upper limits.

A.2. TAROT

The throughput of the TAROT open system was calculated from the detailed information given in Klotz et al. (2008). Panel (b) of Figure 22 shows the resulting sensitivity function, which includes atmospheric transmission typical of La Silla. Fully reduced TAROT images were provided to us by Alain Klotz.

We computed PSF photometry on the images with DAOPhot. Magnitudes for the local sequence stars in the natural system of the TAROT open filter were derived in the same manner as described previously for the LSQ gr filter. Figure 24 shows $({V}_{\mathrm{CSP}}-{\mathrm{open}}_{\mathrm{TAROT}})$ and $({r}_{\mathrm{CSP}}-{\mathrm{open}}_{\mathrm{TAROT}})$ colors for the Pickles main-sequence stars plotted versus ${(g-r)}_{\mathrm{CSP}}$, with the zero point for the openTAROT filter chosen to give a magnitude of 0.0 for Vega. The photometry of the local sequence stars in the field of SN 2012fr is plotted as red points after adjusting the zero points of the openTAROT instrumental magnitudes to provide the best fit to the Pickles stars. The final magnitudes for the local sequence stars in the natural system of the openTAROT filter are foundin Table 9, and the photometry of the SN is given in Table 10.

Figure 24.

Figure 24. Color–color plots for converting CSP V and r magnitudes to the natural system magnitudes for the TAROT open filter. The black curves show synthetic photometry carried out using the Pickles (1998) stellar atlas. The red points correspond to observations of the local sequence stars in the field of SN 2012fr.

Standard image High-resolution image

A.3. Slooh

Panel (c) of Figure 22 shows the throughput functions calculated for the Slooh LRGB filters using information provided by E. Conseil (private communication). We include the atmospheric transmission appropriate for the Teide Observatory. The processed Slooh images were provided by E. Conseil. These displayed a strong gradient in the sky background that was subtracted prior to computing instrumental PSF magnitudes.

Although SN 2012fr is clearly detected in the Slooh images, it has asymmetric, low signal-to-noise ratio PSFs that are unsuitable for measuring magnitudes with our standard photometry tools. Therefore, to analyze these data, it was necessary to develop a specialized tool that produces a three-dimensional model of the PSF using an isolated bright star in the images. The model is background subtracted and subsampled, and it is then fitted to the SN and local sequence stars in the image using an MCMC procedure where the amplitude and the center coordinates are fitted simultaneously. An example of the PSF signal subtraction for the B-band image of the supernova and one of the local standards is shown in Figure 25.

Figure 25.

Figure 25. Examples of PSF subtractions for the Slooh B-band image, with SN 2012fr in the top row, and a star in the bottom row. The left image in both rows shows the original data, the center image shows the scaled PSF model, and the right image shows the subtractions of the scaled PSF model from the original data.

Standard image High-resolution image

The BSlooh and GSlooh filters are moderately well-matched to the CSP B and V filters, as shown in Figures 26(a) and (b). Here, $({B}_{\mathrm{Slooh}}-{B}_{\mathrm{CSP}})$ and $({G}_{\mathrm{Slooh}}-{V}_{\mathrm{CSP}})$ are plotted versus ${(B-V)}_{\mathrm{CSP}}$ for the Pickles stars with the zero points of the BSlooh and GSlooh magnitudes set, assuming ${B}_{\mathrm{Slooh}}\,={V}_{\mathrm{Slooh}}=0.0$ for Vega. The observed colors for the sequence stars, plotted as red points, have been adjusted to fit the synthetic colors by varying the zero points of the instrumental magnitudes. The final magnitudes of the local sequence stars in the natural systems of the BSlooh and GSlooh filters derived from this procedure are given in Table 9. The photometry of the SN in these filters is listed in Table 10.

Figure 26.

Figure 26. Color–color plots for converting (a) BCSP magnitudes to natural system magnitudes in the BSlooh filter, (b) VCSP magnitudes to natural system magnitudes in the GSlooh filter, (c) rCSP magnitudes to natural system magnitudes in the RSlooh filter, and (d) grLSQ magnitudes to natural system magnitudes in the LSlooh filter. The black curves show synthetic photometry carried out using the Pickles (1998) stellar atlas. The red points correspond to observations of the local sequence stars in the field of SN 2012fr.

Standard image High-resolution image

We also attempted to match the RSlooh filter to the rCSP bandpass. The result is shown in the color–color diagram plotted in Figure 26(c), where $({R}_{\mathrm{Slooh}}-{r}_{\mathrm{CSP}})$ is plotted versus ${(g-r)}_{\mathrm{CSP}}$ for the Pickles stars and assuming ${R}_{\mathrm{Slooh}}=0.0$ for Vega. The observed trend that was measured from the photometry of the local sequence stars, as plotted by the red points, is seen to be consistent with the expectation from the synthetic photometry of the Pickles stars after adjusting for the zero point difference. The magnitudes of the local sequence stars in the natural system of the Slooh R filter are given in Table 9, and photometry of the SN in this filter is given in Table 10.

The LSlooh (luminance) filter throughput function closely resembles that of the LSQ gr filter. This is confirmed in the color–color diagram shown in Figure 26(d). Here, $({L}_{\mathrm{Slooh}}-{{gr}}_{\mathrm{LSQ}})$ is plotted versus ${(g-r)}_{\mathrm{CSP}}$ from synthetic photometry of the Pickles main-sequence stars, where the zero point for the LSlooh filter has been calculated assuming LSlooh = 0.0 for Vega. The slope is nearly flat over the whole color range of the Pickles stars, which is consistent with a color term that is essentially zero. The measured colors of the local sequence stars are plotted as red points after adjusting the zero point of the LSlooh instrumental magnitudes to provide the best fit to the Pickles stars. Therefore, we assume that the LSlooh filter is in the same natural system as the grLSQ filter. The photometry of the SN in the LSlooh filter is given in Table 10.

A.4. BOSS

The throughput of the BOSS open system was constructed from information provided by Stuart Parker, who observed SN 2012fr using a Celestron C14 f/10 telescope with a f/6.3 focal reducer attached before a Kodak KAF-3200ME CCD with coverglass. The transmission curves for the StarBright coatings on the telescope optics were found on the Celestron website, but they only covered the wavelength range of 4000–7500 Å. In response to our inquiry, Celestron advised us that the "transmission falloff in the IR and UV bands was pretty severe" (private communication). The transmission of the CCD coverglass is given in a specification sheet that was supplied by the manufacturer, but they only covered the wavelength range from 3500 to 8500 Å. No information on the transmission of the focal reducer could be found. Hence, constructing the BOSS open filter throughput required some guesswork. Panel (d) of Figure 22 shows our final approximation, including the transmission of the atmosphere.

The BOSS images were processed by Stuart Parker using CCDsoft routines. PSF photometry was then carried out with DAOPhot. Final magnitudes for the local sequence stars in the natural system of the BOSS open filter were derived in the same manner as described previously for the LSQ gr filter. Figure 27 shows the ($V-{\mathrm{open}}_{\mathrm{BOSS}}$) and ($r-{\mathrm{open}}_{\mathrm{BOSS}}$) colors for the Pickles main-sequence stars plotted versus ${(g-r)}_{\mathrm{CSP}}$, with the zero point for the synthetic magnitudes chosen to give ${\mathrm{open}}_{\mathrm{BOSS}}=0.0$ for Vega. The local sequence stars in the field of SN 2012fr are plotted as solid red circles in Figure 27 after adjusting the zero points of the openBOSS instrumental magnitudes to provide the best fit to the Pickles stars. Because of the uncertainties involved in constructing the throughput curve of the BOSS open system, these observations of the SN 2012fr local sequence are augmented by photometry (plotted with open red circles) of stars in the field of SN 2014do and a relatively low galactic latitude transient that was observed by both Parker and the CSP-II. In general, the observations are well-matched by the colors derived from the synthetic photometry. This provides confidence that the throughput function shown in panel (d) of Figure 22 is a reasonable representation of the BOSS open filter. The final magnitudes for the local sequence stars in the natural system of the BOSS open filter are given in Table 9 and the photometry of the SN in the natural system of the BOSS open filter is given in Table 10.

Figure 27.

Figure 27. Color–color plots for converting CSP V and r magnitudes to the natural system magnitudes for the BOSS open filter. The black curves show synthetic photometry carried out using the Pickles (1998) stellar atlas. The red points correspond to observations of the local sequence stars in the field of SN 2012fr.

Standard image High-resolution image

Appendix B: Bolometric Light Curve Construction

As discussed in detail by Brown et al. (2016), producing a bolometric light curve from photometry of an object such as a SN Ia with a spectral energy distribution (SED) that differs significantly from those of stars is, at best, an approximate procedure. Consequently, we employed two different methods, which are detailed in this appendix.

B.1. Photometric Trapezoidal Integration

This method estimates the flux density from the observed photometry by simple trapezoidal integration of the fluxes derived from each filter. In addition, estimates must be made for the UV (λ < λu) and infrared (λ > λH) contributions—the former from the Swift uvm2 photometry, and the latter from the Rayleigh–Jeans law. The following steps are involved:

  • 1.  
    The NIR Y JH magnitudes are fitted to match optical phase coverage.
  • 2.  
    The optical magnitudes are interpolated for nights where observations were not available.
  • 3.  
    Corrections are applied for Milky Way and host-galaxy dust reddening.
  • 4.  
    The resulting magnitudes in the natural system are transformed to AB magnitudes.
  • 5.  
    The AB magnitudes are then converted to monochromatic flux densities.
  • 6.  
    The bolometric flux at each phase is derived by integration of the flux densities via the trapezoidal rule.
  • 7.  
    Approximations are applied for the missing flux in the UV (λ < λu) and infrared (λ > λH) for each phase.
  • 8.  
    The resulting bolometric flux versus time is converted to luminosity assuming the Cepheid-based distance modulus for the host-galaxy, NGC 1365.

Note that the K-corrections can be neglected because the redshift of the host-galaxy, NGC 1365, is very small (zhelio = 0.005457).

We used the smooth curve Gaussian process fitting to match the photometry in the NIR Y JH bands to the phases of the optical photometry. The magnitudes were interpolated for the few nights when optical observations were not obtained. Next, the dust extinction from the Milky Way in the direction of SN 2012fr was corrected by subtracting the absorption values given in the third column of Table 11. For most of the filters, these were taken directly from NED. For the Y-band, these were calculated using Fitzpatrick's (1999) galactic-reddening law. Finally, the correction by host-galaxy reddening (see Section 3.2) was performed by subtracting the values given in the final column of Table 11.

Table 11.  Milky Way and Host-galaxy Extinction Corrections and AB Magnitude Offsets

  Effective Milky Way Host-galaxy AB Magnitude
Filter Wavelength (Å) Aλ (mag) Aλ (mag) Offset
uvm2 2360 0.140 0.266 +1.69
u 3550 0.087 0.147 −0.06
B 4300 0.074 0.126 −0.13
g 4700 0.068 0.115 −0.02
V 5400 0.056 0.095 −0.02
r 6200 0.047 0.080 −0.01
i 7500 0.035 0.064 +0.00
Y 10000 0.020 0.038 +0.64
J 12500 0.015 0.026 +0.90
H 16600 0.009 0.017 +1.34

Note. For the uvm2-filter, the AB shift was taken from https://swift.gsfc.nasa.gov/analysis/uvot_digest/zeropts.html. For the CSP-I uBgVriYJH filters, the AB offsets are given in Krisciunas et al. (2017) in Table 16.

Download table as:  ASCIITypeset image

The uvm2, u, B, g, V, r, i, Y, J, and H photometry were then converted to AB magnitudes using the offsets given in the final column of Table 11. Once these are applied, the flux in each band is approximated by:

Equation (9)

The u-to-H flux was then integrated using the trapezoidal rule. The flux beyond the H band was estimated by assuming that it follows the Rayleigh–Jeans law. This leads to the following expression for the integrated flux at wavelengths longward of λH:

We estimate the UV contribution between 1800 and 3000 Å assuming the the SED is flat there and we then estimate the contribution between 3000 Å and λu (3500 Å) by imposing the condition that the flux falls linearly to the flux estimated for uvm2-band from λu. The flux below 1800 Å is assumed to be zero. Hence,

This assumption fits well with the fall off of the flux at wavelengths below the u-band that are typically observed in UV spectra of SNe Ia (e.g., see Foley et al. 2016).

Finally, to calculate the bolometric luminosity in absolute terms, we assumed the Freedman et al. (2001) distance modulus of μ = 31.27 ± 0.05 mag as derived from Cepheid variable observations.

B.2. Spectral Template Fitting

Our second method to estimate the bolometric luminosity takes the spectral template for each phase from Hsiao et al. (2007) and then multiplies it by a function P(λ) such that the synthetic photometry measured on the template exactly matches the real photometry. This method is limited to $t-{t}_{{B}_{\max }}=+79$ because this is the last epoch for the spectral template. The steps can be summarized as follows:

  • 1.  
    The NIR Y JH magnitudes are interpolated to match the optical phase coverage.
  • 2.  
    The optical magnitudes are interpolated for nights where observations were not available.
  • 3.  
    Corrections are applied for Milky Way and host-galaxy dust reddening.
  • 4.  
    For each photometry epoch, the Hsiao et al. spectrum corresponding to the phase that is selected.
  • 5.  
    The template spectrum for each epoch is then matched in flux to the photometry via an iterative procedure:
    • (a)  
      First iteration. The spectrum is divided into wavelength bins corresponding to the non-overlapping filters uvm2, u, g, r, i, Y, J, and H which, in practice, cover almost all the UV-to-NIR wavelength domain, except for two gaps: one between i and Y, and other between J and H. Then, a step-wise function P(λ) is fitted such that:
      Equation (10)
      where X is the filter, mX is the reddening-corrected magnitude of the SN for that filter, SX is the transmission function of the filter, Fλ is the Hsiao et al. spectrum, and zpX is the zero point that was previously adjusted to match our system.
    • (b)  
      Second iteration. The step-wise function that was calculated in Step 5 is now converted to a piece-wise function (or polygonal line) with nodes at the blue and red limits of each filter, except for the u-band where we added an extra node at the effective wavelength to account for the rapid change of flux that occurs across this filter. The values of the nodes for the g-band bin are then imposed to ensure that the slope is determined by the step-wise values of u and r that were measured in the previous step. The nodes of Y and J also adopt the slope of the two corresponding step-wise values measured in the previous step. This smooths the final function P(λ). All of the other nodes are then determined and may be calculated via Equation (10). Figure 28 shows an example of this two-step procedure.
  • 6.  
    The modified Hsiao et al. spectra, $W(\lambda )=P(\lambda ){F}_{\lambda }$, are then integrated between the effective wavelengths λu and λH to calculate the bolometric flux at each phase.
  • 7.  
    For the infrared flux at wavelengths longward of λH, a Rayleigh–Jeans law is again assumed:
  • 8.  
    Finally, the flux is converted into luminosity by assuming the Cepheid-based distance modulus to NGC 1365.

Figure 28.

Figure 28. An example of the two-step process that was devised to match the Hsiao et al. (2007) spectral template to the observed CSP photometry for the same epoch. The dashed line in the lower panel shows the first iteration consisting of a step-wise function, P(λ), as calculated from Equation (10). In the second iteration, a continuous piece-wise function is derived from the step-wise function by forcing the value of the nodes of the g-band bin to have the slope determined by the measured step-wise values of u and r. The nodes for the remaining filters are then calculated from their step-wise values and the shaded regions correspond to the eight bins. The Hsiao et al. template spectrum is then multiplied by the final piece-wise function, as shown in the upper panel.

Standard image High-resolution image

The upper panel of Figure 29 displays the final bolometric light curves calculated using the photometric trapezoidal integration and spectral template fitting methods. In the lower panel of the figure, the ratio of each light curve to the average of the two is plotted versus phase with respect to the epoch of ${t}_{{B}_{\max }}$. It can be seen that the two methods produce light curves that are consistent at the ±5% level. The spectral template method yields a higher luminosity at phases before $t-{t}_{{B}_{\max }}\sim 30$ and a lower luminosity at later epochs. The final bolometric light curves for the two methods are given in Table 12.

Figure 29.

Figure 29. Bolometric light curves calculated with the photometric trapezoidal integration (Method 1) and the spectral template fitting (Method 2) techniques. The inset plot shows luminosity differences of both methods with respect to the average of the two methods.

Standard image High-resolution image

Table 12.  Bolometric Luminosity of SN 2012fr

JD Phase L1(Total) LUVOIR Luvm2 ${L}_{(\lambda \gt {\lambda }_{H})}$ L2(Total) LUVOIR Luvm2 ${L}_{(\lambda \gt {\lambda }_{H})}$
(days) (days) (1043 erg s−1) (%) (%) (%) (1043 erg s−1) (%) (%) (%)
2456230.70 −12.40 0.228(021) 98.15 0.30 1.85 0.229(022) 97.51 0.59 1.91
2456231.75 −11.35 0.362(035) 98.54 0.30 1.46 0.368(037) 97.95 0.57 1.48
2456232.75 −10.35 0.512(050) 98.77 0.32 1.23 0.519(053) 98.15 0.59 1.26
2456233.70 −09.40 0.655(066) 98.90 0.35 1.10 0.661(068) 98.21 0.65 1.14
2456234.70 −08.40 0.809(082) 99.00 0.39 1.00 0.812(085) 98.24 0.73 1.03
2456235.67 −07.43 0.927(095) 99.05 0.44 0.95 0.945(100) 98.22 0.82 0.96
2456236.69 −06.41 1.036(107) 99.09 0.52 0.91 1.063(114) 98.17 0.92 0.91
2456237.60 −05.50 1.117(115) 99.12 0.58 0.88 1.139(122) 98.11 1.00 0.88
2456238.72 −04.38 1.224(127) 99.18 0.66 0.82 1.240(134) 98.03 1.13 0.84
2456239.67 −03.43 1.276(133) 99.21 0.73 0.79 1.307(142) 97.91 1.29 0.81
2456241.73 −01.37 1.332(140) 99.26 0.89 0.74 1.368(151) 97.51 1.73 0.76
2456242.64 −00.46 1.327(139) 99.27 0.96 0.73 1.365(151) 97.32 1.94 0.74
2456243.63 +00.53 1.325(139) 99.30 1.02 0.70 1.350(149) 97.20 2.08 0.72
2456244.73 +01.63 1.307(137) 99.32 1.07 0.68 1.338(148) 97.13 2.18 0.69
2456245.61 +02.51 1.300(137) 99.34 1.09 0.66 1.322(147) 97.09 2.24 0.66
2456246.71 +03.61 1.256(132) 99.34 1.12 0.66 1.279(142) 97.05 2.31 0.64
2456247.65 +04.55 1.222(129) 99.35 1.13 0.65 1.245(138) 97.05 2.33 0.62
2456248.71 +05.61 1.178(124) 99.35 1.12 0.65 1.192(133) 97.01 2.38 0.61
2456249.73 +06.63 1.129(119) 99.35 1.11 0.65 1.148(128) 96.97 2.42 0.61
2456250.67 +07.57 1.070(112) 99.33 1.11 0.67 1.093(122) 96.89 2.48 0.63
2456251.71 +08.61 1.016(106) 99.32 1.09 0.68 1.029(115) 96.79 2.54 0.67
2456252.71 +09.61 0.961(100) 99.29 1.06 0.71 0.974(108) 96.83 2.50 0.67
2456253.64 +10.54 0.903(094) 99.26 1.04 0.74 0.918(102) 96.97 2.39 0.64
2456254.66 +11.56 0.845(088) 99.22 1.02 0.78 0.858(095) 97.05 2.29 0.66
2456255.67 +12.57 0.777(081) 99.15 1.01 0.85 0.783(086) 97.00 2.27 0.73
2456256.68 +13.58 0.726(075) 99.10 0.98 0.90 0.723(079) 96.98 2.21 0.82
2456257.60 +14.50 0.675(069) 99.02 0.96 0.98 0.671(073) 96.94 2.18 0.88
2456258.61 +15.51 0.629(064) 98.94 0.93 1.06 0.622(067) 97.09 1.98 0.93
2456259.59 +16.49 0.593(060) 98.87 0.89 1.13 0.595(063) 97.14 1.88 0.98
2456261.62 +18.52 0.534(053) 98.70 0.81 1.30 0.530(055) 97.32 1.54 1.15
2456262.65 +19.55 0.506(049) 98.60 0.78 1.40 0.503(051) 97.34 1.42 1.24
2456263.65 +20.55 0.488(047) 98.51 0.74 1.49 0.489(049) 97.37 1.32 1.30
2456264.66 +21.56 0.471(044) 98.43 0.70 1.57 0.468(046) 97.36 1.25 1.39
2456265.62 +22.52 0.459(043) 98.35 0.66 1.65 0.458(045) 97.38 1.16 1.46
2456266.63 +23.53 0.447(041) 98.27 0.63 1.73 0.442(042) 97.38 1.05 1.57
2456267.65 +24.55 0.439(039) 98.20 0.59 1.80 0.431(041) 97.39 0.94 1.67
2456268.74 +25.64 0.425(038) 98.11 0.57 1.89 0.416(038) 97.36 0.85 1.79
2456269.78 +26.68 0.416(036) 98.03 0.54 1.97 0.405(037) 97.32 0.80 1.89
2456270.74 +27.64 0.410(035) 97.99 0.52 2.01 0.399(036) 97.27 0.81 1.92
2456271.74 +28.64 0.403(034) 97.93 0.50 2.07 0.390(034) 97.19 0.84 1.97
2456272.74 +29.64 0.395(033) 97.88 0.48 2.12 0.376(033) 97.09 0.88 2.03
2456273.71 +30.61 0.388(032) 97.84 0.47 2.16 0.366(031) 97.00 0.91 2.08
2456274.69 +31.59 0.379(030) 97.80 0.46 2.20 0.359(030) 96.93 0.91 2.16
2456275.72 +32.62 0.368(029) 97.76 0.45 2.24 0.346(029) 96.81 0.92 2.27
2456276.75 +33.65 0.354(028) 97.71 0.45 2.29 0.332(027) 96.67 0.94 2.39
2456278.61 +35.51 0.328(025) 97.62 0.45 2.38 0.305(025) 96.41 0.99 2.61
2456279.74 +36.64 0.312(024) 97.56 0.46 2.44 0.293(024) 96.33 0.99 2.68
2456280.73 +37.63 0.295(022) 97.50 0.47 2.50 0.278(022) 96.20 1.01 2.79
2456281.64 +38.54 0.284(022) 97.48 0.48 2.52 0.270(022) 96.15 1.01 2.84
2456282.76 +39.66 0.268(020) 97.43 0.50 2.57 0.255(020) 96.06 1.02 2.92
2456283.68 +40.58 0.257(020) 97.42 0.51 2.58 0.244(020) 95.99 1.01 3.00
2456284.68 +41.58 0.246(019) 97.41 0.52 2.59 0.234(019) 95.95 1.00 3.05
2456285.65 +42.55 0.235(018) 97.41 0.53 2.59 0.224(018) 95.90 0.99 3.11
2456286.65 +43.55 0.227(017) 97.43 0.54 2.57 0.215(017) 95.88 0.98 3.15
2456287.75 +44.65 0.218(017) 97.45 0.55 2.55 0.206(017) 95.84 0.98 3.18
2456288.73 +45.63 0.210(016) 97.48 0.55 2.52 0.197(016) 95.79 1.02 3.19
2456289.60 +46.50 0.204(016) 97.51 0.56 2.49 0.192(016) 95.82 1.03 3.15
2456292.65 +49.55 0.182(014) 97.57 0.59 2.43 0.169(014) 95.62 1.21 3.17
2456293.66 +50.56 0.176(014) 97.60 0.60 2.40 0.163(014) 95.60 1.25 3.15
2456294.68 +51.58 0.169(014) 97.61 0.61 2.39 0.156(013) 95.52 1.31 3.16
2456295.63 +52.53 0.165(013) 97.65 0.61 2.35 0.153(013) 95.53 1.35 3.13
2456296.60 +53.50 0.160(013) 97.67 0.61 2.33 0.148(013) 95.54 1.37 3.09
2456297.67 +54.57 0.154(013) 97.69 0.62 2.31 0.144(012) 95.55 1.39 3.06
2456298.53 +55.43 0.150(012) 97.70 0.62 2.30 0.140(012) 95.57 1.39 3.04
2456299.69 +56.59 0.145(012) 97.72 0.62 2.28 0.136(012) 95.57 1.42 3.01
2456304.54 +61.44 0.126(011) 97.77 0.63 2.23 0.119(011) 95.62 1.49 2.90
2456305.60 +62.50 0.122(010) 97.78 0.63 2.22 0.116(010) 95.66 1.49 2.85
2456308.55 +65.45 0.112(010) 97.80 0.63 2.20 0.107(010) 95.76 1.49 2.75
2456311.61 +68.51 0.104(009) 97.85 0.63 2.15 0.099(009) 95.92 1.50 2.58
2456314.59 +71.49 0.096(009) 97.91 0.63 2.09 0.091(008) 96.07 1.50 2.42
2456315.63 +72.53 0.093(008) 97.91 0.63 2.09 0.088(008) 96.18 1.53 2.29
2456316.56 +73.46 0.091(008) 97.95 0.63 2.05 0.086(008) 96.23 1.52 2.25
2456317.58 +74.48 0.089(008) 97.99 0.63 2.01 0.084(008) 96.34 1.52 2.14
2456318.59 +75.49 0.087(008) 98.02 0.63 1.98 0.082(008) 96.38 1.52 2.10
2456320.57 +77.47 0.082(007) 98.08 0.63 1.92 0.079(007) 96.47 1.53 2.00
2456321.56 +78.46 0.080(007) 98.11 0.64 1.89 0.076(007) 96.50 1.54 1.96
2456322.59 +79.49 0.078(007) 98.16 0.63 1.84 0.075(007) 96.58 1.53 1.89
2456323.62 +80.52 0.076(007) 98.20 0.63 1.80 ... ... ... ...
2456324.57 +81.47 0.074(007) 98.22 0.64 1.78 ... ... ... ...
2456325.56 +82.46 0.072(007) 98.27 0.64 1.73 ... ... ... ...
2456326.58 +83.48 0.070(007) 98.30 0.65 1.70 ... ... ... ...
2456327.56 +84.46 0.068(006) 98.34 0.65 1.66 ... ... ... ...
2456328.57 +85.47 0.066(006) 98.36 0.66 1.64 ... ... ... ...
2456329.56 +86.46 0.065(006) 98.41 0.66 1.59 ... ... ... ...
2456330.55 +87.45 0.063(006) 98.43 0.67 1.57 ... ... ... ...
2456331.53 +88.43 0.062(006) 98.47 0.67 1.53 ... ... ... ...
2456332.52 +89.42 0.061(006) 98.51 0.67 1.49 ... ... ... ...
2456336.56 +93.46 0.056(005) 98.62 0.69 1.38 ... ... ... ...
2456337.57 +94.47 0.054(005) 98.62 0.71 1.38 ... ... ... ...
2456340.57 +97.47 0.051(005) 98.69 0.71 1.31 ... ... ... ...
2456342.60 +99.50 0.049(005) 98.72 0.73 1.28 ... ... ... ...
2456343.58 +100.48 0.048(005) 98.72 0.74 1.28 ... ... ... ...
2456345.57 +102.47 0.046(004) 98.74 0.75 1.26 ... ... ... ...
2456346.58 +103.48 0.045(004) 98.75 0.76 1.25 ... ... ... ...
2456349.58 +106.48 0.042(004) 98.76 0.79 1.24 ... ... ... ...
2456351.58 +108.48 0.041(004) 98.77 0.81 1.23 ... ... ... ...
2456352.58 +109.48 0.040(004) 98.76 0.82 1.24 ... ... ... ...
2456354.55 +111.45 0.038(004) 98.75 0.85 1.25 ... ... ... ...
 
L(1,2): Bolometric luminosity in the wavelength range λ > 1800 Å.
LUVOIR: Luminosity in the wavelength range 3000 < λ < 16600 Å.
Luvm2: Luminosity in the wavelength range 1800 < λ < 3000 Å.
${L}_{\lambda \gt {\lambda }_{H}}$ : Luminosity in the wavelength range λ > 16600 Å.

Note. The phase is relative to ${t}_{{B}_{\max }}=\mathrm{JD}\,2456243.1$. LUVOIR is the luminosity measured from the flux of 2012fr through uBgVriYJH-bands photometry.

Download table as:  ASCIITypeset images: 1 2

Footnotes

  • This paper includes data gathered with the 6.5 m Magellan Baade Telescope, located at Las Campanas Observatory, Chile.

  • 25 

    Freedman et al.'s (2001) distance modulus for NGC 1365 is consistent within the errors with those obtained more recently by Riess et al. (2016) from Cepheid variables and by Jang et al. (2018) from the tip of the red giant branch method.

  • 26 
  • 27 

    The decline-rate parameters, Δm15(X), are here defined as the magnitude difference of the light curve in filter X from peak to 15 days later. Historically, the value of ${\rm{\Delta }}{m}_{15}(B)$ was shown by Phillips (1993) to correlate with the absolute peak magnitude in such a way that more luminous objects show slower light curve decline rates.

  • 28 

    NED is the NASA/IPAC Extragalactic Database.

  • 29 

    A pseudo-color is defined as the difference between peak magnitudes of two passbands. In the case of SNe Ia, the time of peak brightness can vary by up to a few days from passband to passband.

  • 30 

    Note that the LVG subtype as defined by Benetti et al. includes 1991T-like objects.

  • 31 

    According to Finlator et al. (2000), 99% of the field stars observed in the Sloan Digital Sky Survey are on the main sequence.

  • 32 
Please wait… references are loading.
10.3847/1538-4357/aabaf8