Free Access
Issue
A&A
Volume 641, September 2020
Article Number A148
Number of page(s) 17
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202038118
Published online 23 September 2020

© ESO 2020

1. Introduction

A Type IIn supernova (SN) exhibits conspicuous, narrow Balmer lines in emission (Schlegel 1990) formed as a result of the interaction of rapidly expanding SN ejecta with circumstellar matter (CSM, e.g., Chugai & Danziger 1994). Progenitors of SNe IIn experience significant mass loss during their pre-SN evolution, leading to a local environment that contains dense CSM (e.g., Ofek et al. 2013, 2014a; Fraser et al. 2013). The direct detection of SNe IIn progenitors in pre-explosion imaging (Gal-Yam & Leonard 2009; Prieto et al. 2008) and environmental studies around SNe IIn (Anderson et al. 2012; Habergham et al. 2014; Taddia et al. 2015; Galbany et al. 2018) reveals that a range of masses, including both relatively low-mass progenitors such as red supergiants and high-mass progenitors such as luminous blue variables, serve as the progenitors of SNe IIn.

SNe IIn exhibit a heterogeneous light-curve (LC) evolution (e.g., Nyholm et al. 2020; Taddia et al. 2013, 2015; Stritzinger et al. 2012; Kiewe et al. 2012). Some SNe IIn have long-lasting emission and their declines are well-fit by a single power-law (cf. SN 1988Z; Stathakis & Sadler 1991; Turatto et al. 1993). Others exhibit rapid, exponentially declining LCs (cf. SN 1998S; Fassia et al. 2000). Some SN IIn have slowly evolving, Gaussian-shaped LC peaks (cf. SN 2006gy; Smith et al. 2007), while other SN IIn exhibit LCs with a plateau that is akin to SNe IIP (cf. SN 1994W; Sollerman et al. 1998). Several SNe IIn exhibit multiple luminosity peaks (Stritzinger et al. 2012; Nyholm et al. 2017) and some even have documented rise times in excess of 100 days (Miller et al. 2010; Moriya et al. 2019). This LC diversity largely originates from differences in the properties of the progenitors and CSM (Moriya et al. 2013a,b; Chatzopoulos et al. 2013; Ofek et al. 2014b; Dessart et al. 2016; Tsuna et al. 2019; Takei & Shigeyama 2020), compounded by viewing angle effects (e.g., Suzuki et al. 2019).

Another important characteristic of SNe IIn is that they often exhibit a number of observational properties linked to the presence of dust. The dust in SNe IIn is revealed by the infrared (IR) luminosity evolution (e.g., Fox et al. 2011, 2013; Szalai et al. 2019), as well as by the evolution of prominent emission features (e.g., Pozzo et al. 2004; Taddia et al. 2020). The dust in SNe IIn has been well-traced in SN 2010jl. SN 2010jl is a slowly-declining 1988Z-like object that exhibits a power-law luminosity decline at early times (Stoll et al. 2011; Zhang et al. 2012; Maeda et al. 2013; Fransson et al. 2014; Ofek et al. 2014c; Jencson et al. 2016). Although SN 2010jl was observed intensively over a range of wavelengths (Chandra et al. 2012, 2015; Williams & Fox 2015; Bevan et al. 2020), the origin of the dust traced by its IR observations is still a matter of debate. The dust might have formed within a cool dense shell located within the wake of the supernova shock wave which is positioned at the interface between the rapidly expanding supernova ejecta and the dense CSM (e.g., Smith et al. 2011; Gall et al. 2014; Maeda et al. 2013). Alternatively, the IR emission could be attributed to pre-existing dust located around the progenitor (Andrews et al. 2011; Fransson et al. 2014; Sarangi et al. 2018; Chugai 2018). The possibility that most IR emission observed in SNe IIn originates from pre-existing dust has also been proposed (e.g., Fox et al. 2011). To clarify the origin of dust in SNe IIn, more detailed observations of SNe IIn covering a wide wavelength range are required.

In this paper, we report on the observations of the long-lasting (1988Z-like) Type IIn SN 2014ab, initiated by the Carnegie Supernova Project-II (CSP-II; Phillips et al. 2019; Hsiao et al. 2019). Although this object spectroscopically resembles the well-observed SN 2010jl, it does exhibit different LC evolution and dust signatures. SN 2014ab offers a new and important clue on the nature of slowly-declining SNe IIn, including details about the dust in their nearby environment.

2. SN 2014ab: discovery, distance, and reddening

The Catalina Real-time Transient Survey (CRTS) reported the discovery of a transient with an apparent magnitude of 16.4 in MCG +01-35-037 (Howerton et al. 2014). First designated PSN J13480599+0723164 when it was posted to the web page of the Central Bureau for Astronomical Telegrams, the transient was identified by CRTS in an image taken by the Catalina Sky Survey (CSS) on March 9.43 UT (i.e., JD-2456725.93). With coordinates RA (J2000.0) = 13h48m05s.99 and Dec (J2000.0) = +07°23′16 . $ {{\overset{\prime\prime}{.}}} $40, the transient (aka SNhunt237) is located 1 . 5 $ 1{{\overset{\prime\prime}{.}}}5 $ north of the host galaxy (see Fig. 1).

thumbnail Fig. 1.

Finding chart constructed using a single r-band image obtained on 03 March 2014 with the Swope 1.0 m telescope. The position of SN 2014ab is indicated with a dotted circle and the approximate area contained within the circle is shown in the upper right inset Optical and NIR local sequence stars are indicated with yellow and red circles, respectively.

Within 24 h of the discovery (Mar. 10.25 UT), the CSP-II obtained a near-IR (NIR) spectrum with the Folded-port Infrared Echellette (FIRE; Simcoe et al. 2013) spectrograph attached to the Magellan Baade 6.5 m Telescope located at Las Campanas Observatory (LCO). This spectrum revealed an SN IIn and upon the reporting of its spectral classification, the transient obtained the International Astronomical Union (IAU) designation of SN 2014ab. The CSP-II classification was promptly confirmed by PESSTO (Public ESO Spectroscopic Survey for Transient Objects; Smartt et al. 2015) based on a visual-wavelength spectrum obtained with the New Technology Telescope (+ EFOSC; Buzzoni et al. 1984), located at the ESO La Silla Paranal Observatory (Fraser et al. 2014).

According to the NASA Extracgalactic Database (NED), the host-galaxy MCG +01-35-037 has a redshift z = 0.02352. Turning to the medium-dispersion spectra of SN 2014ab, obtained with the ESO La Silla Paranal Observatory’s Very Large Telescope (VLT) equipped with the X-shooter (see below), we obtain an averaged redshift of z = 0.02262 ± 0.00001 from the narrow Hα emission component in the three spectra. Adopting H0 = 73.2 ± 2.3 km s−1 Mpc−1 (Burns et al. 2018), ΩM = 0.27, and ΩΛ = 0.73, this corresponds to a luminosity distance of 94.3 ± 6.6 Mpc, or a distance modulus μ = 34.87 ± 0.15 mag.

According to Schlafly & Finkbeiner (2011) and reported in NED1, the Milky Way visual-band extinction in the direction of SN 2014ab is A V MW =0.083 $ A^{MW}_V = 0.083 $ mag. This value is inferred from infrared-based dust maps and assumes a Fitzpatrick (1999) reddening law with RV = 3.1.

An inspection of the visual-wavelength spectra of SN 2014ab presented below reveals no Na I D features at the redshift of the host galaxy and there is no evidence of the 5780 Å diffuse interstellar band. We assume, therefore, in the following that SN 2014ab suffers from minimal to no host-galaxy reddening.

3. Observations

3.1. WISE mid-infrared photometry

We include ten epochs of two-channel mid-IR (MIR) imaging of SN 2014ab serendipitously obtained by the Wide-field Infrared Survey Explorer (WISE) satellite (Wright et al. 2010; Mainzer et al. 2014), and downloaded from NASA/IPAC Infrared Science Archive2. Photometry retrieved from the archive includes W1-band measurements at 3.4 μm and W2-band measurements at 4.6 μm. This photometry was previously presented by Jiang et al. (2019). The first detection by WISE is at JD 2456671.59, which is 56 days before the optical discovery at JD 2456725.93 in the observer frame. The WISE photometry is listed in Table A.1.

3.2. Carnegie Supernova Project-II observations of SN 2014ab

The CSP-II obtained broad-band optical and NIR imaging with facilities located at LCO. In particular, twenty-two epochs of optical (ugriBV) imaging was taken with the Henrietta Swope 1.0 m telescope equipped with a direct optical camera containing an e2V CCD, while five epochs of NIR images were taken with the Irénée du Pont 2.5 m telescope equipped with RetroCam (Hamuy et al. 2006). The ugriBV-band photometry extends, relative to the first WISE detection, from +56 d to +485 d, with a significant gap between +140 d to +360 d when SN 2014ab was located behind the Sun. The NIR photometry is more limited in coverage, ranging between +56 d to +121 d. All dates presented here are in the observer frame.

The CSP-II data were reduced to photometry following standard prescriptions that are well documented in Krisciunas et al. (2017, and references therein). Prior to computing the photometry of the SN from each science image, deep optical and NIR templates of the host were obtained with the du Pont telescope equipped with a direct CCD camera and with the Magellan Baade telescope equipped with the NIR imager FourStar (Persson et al. 2013).

PSF photometry of SN 2014ab was computed from the template-subtracted science images and its instrumental magnitudes were calibrated relative to a local sequence of stars located in the field of the host. The optical local sequence consists of twenty-three stars and the NIR local sequence consists of six stars. The optical local sequence was calibrated relative to standard star fields observed over a minimum of three photometric nights, while the NIR sequence were calibrated relative to Persson et al. (1998) (JH band) and Krisciunas et al. (2017) (Y band) standard star fields observed over four and three photometric nights, respectively. Optical and NIR photometry of the local sequence stars in the “standard” photometric systems are listed in Table A.2 and Table A.3, respectively. Finally, the optical and NIR photometry of SN 2014ab in the natural system are listed in Table A.4.

The CSP-II obtained a single visual-wavelength spectrum of SN 2014ab with the du Pont telescope equipped with Wide Field Re-imaging CCD Camera (WFCCD) and four epochs of NIR spectra with the Magellan Baade telescope equipped with the FIRE spectrograph. In addition to these data, we obtained optical/NIR spectra from the VLT (+ X-shooter, Vernet et al. 2011; Program ID: 092.D-0645, PI Gall), along with a single NTT (+ EFOSC) spectrum reported by PESSTO. A journal of the spectroscopic observations is provided in Table 1. The X-shooter data were obtained in nodding (ABBA) mode and underwent basic reduction steps in “stare” mode for the UVB and VIS arms and “nodding” mode for the NIR arm using the EsoReflex workflow (Freudling et al. 2013). The UVB and VIS single-stare frames were combined using custom python scripts3 (Selsing et al. 2019) and optically extracted, slit loss corrected, and corrected for heliocentric velocity using a custom IDL programs. The telluric correction was performed by using Molecfit (Smette et al. 2015; Kausch et al. 2015). The other spectroscopic data were reduced following standard procedures as previously described by Hamuy et al. (2006) (optical) for and Hsiao et al. (2019) (NIR). Finally, the flux scale of each of the one-dimensional, visual-wavelength spectra were adjusted such that their synthetic colors match the more accurate broad-band colors computed from photometry (see, e.g., Boldt et al. 2014).

Table 1.

Journal of spectroscopic observations of SN 2014ab.

In addition to the spectra of SN 2014ab, we also present CSP-II spectra of SN 2010jl in Appendix B. As demonstrated below, SN 2014ab shares some common features with SN 2010jl. These spectra were reduced following the same manner as applied to the data of SN 2014ab.

4. Results

4.1. Photometry and B  −  V color-curve evolution

Plotted in Fig. 2 are the two-channel MIR LCs computed from WISE observations covering over 1600 days of evolution. The object reached a maximum between 200 and 400 days post-detection, followed by a slow decline (∼0.0015 mag day−1) over the next 1200 days. The object appears to have reached a W1-band peak brightness of ≤ − 22.5 mag and the W2-band reached a peak of ≤ − 23 mag.

thumbnail Fig. 2.

Mid-infrared light curves of SN 2014ab obtained with the WISE satellite. The open symbols show the raw data and the filled symbols show co-added magnitudes of the nearby epochs. We also present the photometry of KISS15s, which is another SN IIn with serendipitous WISE observations (Kokubo et al. 2019). The KISS15s photometry is plotted with the right and top axes. The time origin in KISS15s is its optical discovery.

The CSP-II optical and NIR LCs of SN 2014ab are plotted in Fig. 3. The optical photometry extends from +56 d to +485 d relative to the first WISE detection in the observer frame, while the NIR photometry extends between +56 d to +121 d. The photometry shown in the figure has been corrected for Galactic reddening, and placed on the absolute magnitude scaling adopting the distance computed in Sect. 2. The optical LCs at the time of discovery range from −17.5 mag to −19 mag, while the NIR absolute magnitudes reach around −20 mag. In the initial phases of our follow-up, the optical/NIR LCs slowly evolve until +123 d when follow-up ceased for ∼200 days. For example, the g-band LC declines 0.36 mag per 100 days. Follow-up then continued at around +350 d and extended until around +480 d. During this time the g-band LC decline rate was 1.4 mag per 100 days.

thumbnail Fig. 3.

Multi-color, optical (ugriBV) and NIR (YJH) light curves of SN 2014ab obtained by the CSP-II. The data have been corrected for Galactic extinction. Vertical lines indicate epochs that optical (green), NIR (red), or optical+NIR (blue) spectroscopic observations were obtained.

Figure 4 presents the B − V color evolution of SN 2014ab and other representative SNe IIn. SN 2014ab has little color evolution during our observations lasting for about 100 days. A lack of significant color evolution appears to be a hallmark of other slowly declining SNe IIn such as SN 1988Z and SN 2010jl. In particular, the color-curve evolution of SN 2010jl is very much like that of SN 2014ab. As discussed later, the spectra of SN 2014ab are also very similar to those of SN 2010jl.

thumbnail Fig. 4.

Observed B − V color curve of SN 2014ab compared with those of three different SNe IIn subtypes. The reddening is corrected based on the values reported in the literature. The explosion date, i.e., day zero, for SN 2014ab is arbitrarily set to be 100 days before the optical discovery, or 46 days before the WISE first detection in the rest frame. The comparison sample including objects of the slowly-declining 1988Z-like subtype (solid lines) including SN 2010jl (Fransson et al. 2014), SN 2006tf (Smith et al. 2008), and SN 1988Z (Turatto et al. 1993), as well as the 1994W-like plateau subtype (dashed lines) including SN 1994W (Sollerman et al. 1998), SN 2009kn (Kankare et al. 2012), and SN 2011ht (Mauerhan et al. 2013), and finally, the prototypical rapidly-declining SN 1998S (dot-dashed line, Fassia et al. 2000).

4.2. Spectral energy distributions and the UVOIR light curve

Here, we use the broad-band optical through MIR photometry of SN 2014ab to construct spectral energy distributions (SEDs) for each observed epoch. These are then fit with black-body (BB) functions in order to estimate the BB temperature (TBB), BB radius (RBB), and BB luminosity (LBB) profiles associated with the underlying emission region(s).

Plotted in Fig. 5, we present the SEDs of SN 2014ab, extending from +56 d to +485 d relative to the first WISE detection. To build these SEDs, first the photometry at each epoch was corrected for reddening and then each optical and NIR point was placed on the AB system using the AB offsets given in Krisciunas et al. (2017); see their Table 16. Then each photometric point was converted to flux following Eq. (F5) of Krisciunas et al. (2017). If there were gaps in the optical or NIR coverage at any particular epoch, they were estimated from neighboring photometry through the use of linear interpolation. The WISE photometry, which is given in the Vega system, was converted to flux using the flux zero points given in Wright et al. (2010).

thumbnail Fig. 5.

Spectral energy distributions (SEDs) of SN 2014ab. The first 16 quadrants from the top contain SEDs constructed from flux points covering optical, NIR, and MIR wavelengths and extend from +56 d to +140 d relative to the time of first WISE detection. Each of these SEDs are fit with a single black body (BB) function, where the optical flux is dominating. We excluded u and r band from the fit due to line blanketing and the strong Hα emission line, respectively. The SEDs plotted in the other quadrants are constructed from optical and MIR flux points and extend from +360 d to +485 d relative to first WISE detection. The late-phase SEDs clearly exhibit an excess of flux relative to a single BB fit to the optical flux points. We therefore fit the late SEDs with a two-component BB function. In the last quadrant the optical BB fit is not reliable as only g and i are available.

Next, each SED was fit with a BB function. When we compute the best-fit BB functions to the SEDs, the u- and r-band flux points were excluded because they are affected by significant line blanketing and the strong Hα emission line, respectively (Sect. 4.3). We find that the obtained BB fit matches well to the continuum of the spectra when available. The BB temperature at the early epochs is around 8000 K. An inspection of Fig. 5 reveals that from the +360 d SED and onwards, an excess of MIR flux is present relative to expectations of a single BB function fit (blue dashed line). We therefore fit the SEDs extending from +360 d to +484 d with an additional BB-function component, which is plotted as red dashed line, while the combined BB functions are overplotted in magenta. We refer to the first hot BB component as BBhot and the second warm BB component as BBwarm. BBhot remains around 8000 K, although the uncertainty is large at the late epochs. We cannot tell if BBwarm is affected by the presence of emission lines in the MIR bands because we lack MIR spectra. We therefore turn to the MIR spectra of the intermediate-luminosity-red-transient NGC 300-2008-OT, which were obtained around 400 days and 600 days (Ohsawa et al. 2010). This gap transient is thought to have been powered by circumstellar interaction (Bond et al. 2009; Berger et al. 2009). These spectra reveal no prevalent emission features. The BBwarm component fit to the MIR flux points of SN 2014ab suggests T BB warm 2000 $ T_{BB}^{\mathrm{warm}} \sim 2000 $ K, which is similar to the value inferred for SN 2010jl (1500 − 2000 K, Fransson et al. 2014). We note that T BB warm $ T_{\mathrm{BB}}^{\mathrm{warm}} $ is mostly constrained by the two MIR photometric bands. We also tried BB fits by forcing T BB warm = 1000 K $ T_{\mathrm{BB}}^{\mathrm{warm}}=1000\,\mathrm{K} $ and found that the BB fit significantly overestimates the W2 photometry. Therefore, T BB warm $ T_{\mathrm{BB}}^{\mathrm{warm}} $ would not be far from 2000 K (≃1500 − 2000 K) and is close to the dust evaporation temperature.

Next, by integrating the BB function fits we obtain the UVOIR LC of SN 2014ab as plotted in the top panel of Fig. 6. Specifically, we plot L BB hot $ L^{\mathrm{hot}}_{\mathrm{BB}} $ (blue), L BB warm $ L^{\mathrm{warm}}_{\mathrm{BB}} $ (red) and the summation of the two (black). We find that on +56 d, SN 2014ab reached L BB tot 1.0 × 10 43 $ L^{\mathrm{tot}}_{\mathrm{BB}} \approx 1.0\times10^{43} $ erg s−1, and by +360 d this dropped by an order of magnitude. Furthermore, we find that at +360 d that the L BB warm $ L^{\mathrm{warm}}_{\mathrm{BB}} $ dominates over the emission of L BB hot $ L^{\mathrm{hot}}_{\mathrm{BB}} $ by a factor of ∼3 and this increases to a factor of ∼5 by +485 d. Comparison of the T BB warm $ T^{\mathrm{warm}}_{\mathrm{BB}} $ and T BB hot $ T^{\mathrm{hot}}_{\mathrm{BB}} $ values reveals that the former characterized by T BB warm 2000 $ T^{\mathrm{warm}}_{\mathrm{BB}} \sim 2000 $ K is around a factor of ∼4 lower than the latter over the duration covered by our late phase data. Over the same time, R BB warm $ R^{\mathrm{warm}}_{\mathrm{BB}} $ is found to be around a factor of ∼6 higher than R BB hot $ R^{\mathrm{hot}}_{\mathrm{BB}} $. These attributes suggest that the emission of the BBwarm component is related to dust reprocessing of shorter wavelength photons to longer wavelength photons.

thumbnail Fig. 6.

Bolometric LC (top) and the BB-temperature (middle) and BB-radius (bottom) profiles of SN 2014ab, computed from fitting single and double BB functions to the SEDs of SN 2014ab.

The radiated energy from +56 d to +140 d in Fig. 6 is 6.5 × 1049 erg and that from +360 d to +455 d is 6.2 × 1048 erg. The total radiated energy during our observation is therefore 7 × 1049 erg. Because of the limited LC coverage of our observation, the actual radiated energy should be much larger than 7 × 1049 erg.

The bolometric luminosity we presented is obtained by fitting BB functions to the SEDs from photometry. When we have both spectra and photometry, we compare the results of the BB fits from spectra and photometry. We found that the uncertainty from our use of the photometry to fit the BB function is up to 10%. The systematic uncertainty caused by the use of the two BB functions to fit the SEDs is, however, uncertain.

4.3. Spectroscopy

Plotted in Fig. 7, we present our visual-wavelength and NIR spectra of SN 2014ab extending from +54 d to +183 d, relative to the first WISE detection in the rest frame. The spectra are dominated by narrow emission lines, justifying a SN IIn classification.

thumbnail Fig. 7.

Optical through NIR wavelength spectroscopy of SN 2014ab, extending from +54 d to 183 d (restframe days) past the first WISE detection. The numbers in the parentheses indicate the flux offset adopted for the spectra in the unit of 10−16 erg s−1 cm−2 Å−1. The shaded wavelength regions indicate areas of strong telluric absorption. The last three spectra are not calibrated by photometry because of the lack in the photometric information.

Figure 8 shows the temporal evolution of prevalent hydrogen emission-line profiles. The hydrogen emission-line profiles exhibit little evolution over the first 120 days of spectroscopic observations. The broad component matches the Lorentzian function with a full-width half maximum (FWHM) velocity of 3500 km s−1. This Lorentzian line profile is formed by the multiple electron scattering in the optically thick CSM (e.g., Chugai 2001, but see also Huang & Chevalier 2018). In superposing the broad component, we find a narrow emission component as shown in the top panel of Fig. 8. The peak of the narrow Hα component in the X-shooter spectra match well and are used to determine the redshift of SN 2014ab (z = 0.02262 ± 0.00001). The difference between the host galaxy redshift (z = 0.02352) and the SN redshift (z = 0.02262) corresponds to the line-of-sight velocity difference of 270 km s−1. The difference may originate from the host galaxy’s rotation and is within the normal distribution of radial velocity of SN host galaxies (Galbany et al. 2018). We can also see a tentative absorption component in the narrow Hα lines with a minimum at −80 km s−1. This feature could be produced from an unshocked wind with a characteristic velocity of 80 km s−1. Such an absorption with the minimum velocity of arouund 100 km s−1 is often observed in SNe IIn (e.g., Taddia et al. 2013). The peak of the broad component is shifted by −500 km s−1.

thumbnail Fig. 8.

Time-evolution of the narrow- (top panel) and broad-component (middle panel) Hα profiles, and the broad-component (bottom panel) of the Pa-β profile, plotted in velocity space. Over-plotted as dashed cyan lines are Lorentzian functions with the FWHM velocity of 3500 km s−1 and the peak shifted by −500 km s−1. The redshift is determined by the average of the narrow emission peak. Both broad components are overall symmetric and do not exhibit significant evolution over the duration of our followup observations.

The broad emission lines show little evolution during our spectroscopic observations (Fig. 8). The emission-line luminosities of the representative Balmer lines are ≃7 × 1041 erg s−1 (Hα) and ≃2 × 1041 erg s−1 (Hβ). The flux ratio of Hα to Hβ is ≃4, and is similar to that of SN 2010jl at around 100 days after the estimated explosion date (Fransson et al. 2014). Similarly, the emission line luminosities of the representative Paschen lines are ≃2 × 1041 erg s−1 (Paα) and ≃9 × 1040 erg s−1 (Paβ). The luminosity of Paα is difficult to measure because of the strong telluric absorption and we therefore measure it by using the +97 d spectrum which has a better telluric correction.

5. Discussion

5.1. Comparison of SN 2014ab to other SNe IIn

We now compare SN 2014ab with other SNe IIn. The explosion date of SN 2014ab is not well-constrained. Here, we arbitrarily assume that the explosion date is 100 days before the optical discovery in the rest frame which is 46 days before the first WISE detection. With this assumption on the explosion date, the LC of SN 2014ab matches that of SN 2010jl well, which is shown to be very similar to SN 2014ab later. Figure 9 presents the LC comparison. In the first 200 days, the LC of SN 2014ab declines slowly and it resembles slowly declining SNe IIn such as SN 1988Z and SN 2010jl. SN 2014ab is significantly brighter than SN 1988Z and it is slightly fainter than SN 2010jl. The optical LC of SN 2014ab is found to be similar to that of KISS15s in the early phases (Kokubo et al. 2019). However, their MIR LCs have different properties (Fig. 2). Specifically, the MIR luminosity of KISS15s rises gradually for 700 days, while in the case of SN 2014ab, the MIR luminosity increases suddenly at around +300 d and subsequently gradually declines.

thumbnail Fig. 9.

Absolute V-band LC of SN 2014ab compared with well-observed examples of the three subtypes of SN IIn discussed by Taddia et al. (2013). This includes the long-lasting 1988Z-like objects (circles) SN 1988Z (Turatto et al. 1993), SN 2006tf (Smith et al. 2008), SN 2010jl (Stoll et al. 2011; Zhang et al. 2012; Maeda et al. 2013; Fransson et al. 2014), and KISS15s (Kokubo et al. 2019), the 1994W-like objects (squares) including SN 1994W (Chugai et al. 2004), SN 2009kn (Kankare et al. 2012), and SN 2011ht (Mauerhan et al. 2013), which show a LC pleateau, and finally, the prototypical, rapidly-declining SN 1998S (triangles, Fassia et al. 2000).

One important characteristic that makes SN 2014ab distinct from other slowly declining SNe IIn is the optical luminosity drop that occurred between +200 d and +400 d (Fig. 9). The optical LCs of slowly evolving SNe IIn sometimes start to decline rapidly from several hundred days after the explosion. The optical LC drop displayed by SN 2014ab is the most abrupt among slowly declining SNe IIn. The change in the optical LC decline rate in slowly-declining SNe IIn is often accompanied by an IR luminosity increase. Indeed, the sudden optical luminosity decline in SN 2014ab follows after an increase in luminosity at MIR wavelengths.

Figure 10 shows the spectral comparison of SN 2014ab with different SNe IIn subtypes. This includes the slowly declining 1988Z-like Type IIn SN 2010jl, the rapidly declining Type IIn SN 1998S and, finally, the 1994W-like Type IIn SN 2009kn with its LC akin to normal SNe IIP. The spectrum of SN 2010jl closely resembles that of SN 2014ab, with both SNe having a very similar shape of their continua and the same spectral features. Furthermore, the hydrogen-line flux ratios measured from the spectra of both objects are nearly identical. Comparison of the spectra of these two objects with that of SN 2009kn reveals that the latter has much narrower emission-line features. Finally, the spectrum of SN 1998S is quite distinct from that of SN 2014ab. The spectral comparison also confirms that SN 2014ab is a member of the slowly declining, 1988Z-like subclass of SNe IIn.

thumbnail Fig. 10.

Comparison of our +67 d visual-wavelength (top) and NIR-wavelength (bottom) spectrum of SN 2014ab with other SN IIn subtypes. Comparison objects are the slowly declining SN 2010jl (Borish et al. 2015), the 1994W-like SN 2009kn (Kankare et al. 2012), and the rapid-declining SN 1998S (Fassia et al. 2001). Prominent ions are indicated and labeled. Regions suffering from strong telluric lines are shaded. Note that the NIR spectrum of SN 1998S is composed of data taken on two different but close in time epochs as indicated. Note that the previously unpublished visual-wavelength spectrum of SN 2010jl was obtained by the CSP-II (see Appendix B for details).

Inspection of the emission-lines associated with Hα and the other hydrogen lines in SN 2014ab have the centroid of the broad emission line blue-shifted by 500 km s−1 as shown in Fig. 8. A similar blue-shift of the order of 100 km s−1 in the emission peak is also found in SN 2010jl, although the shift velocity is time-dependent (Gall et al. 2014; Fransson et al. 2014). In the case of SN 2014ab, we do not find the time dependent shift of the peak of the broad component as in SN 2010jl. In addition, we do not find the bluer emission peaks in bluer emission lines as observed in SN 2010jl. Indeed, the shape of the hydrogen emission lines in the spectra of SN 2014ab do not strongly depend on either emission wavelengths or time (Fig. 11).

thumbnail Fig. 11.

Peak of the hydrogen emission lines. We use the optical spectrum at +67 d and the NIR spectrum at +69 d. The peak of the broad component does not show a wavelength-dependent shift expected if dust is formed at the emission line forming region, i.e., the emission lines at the shorter wavelengths to have more blue-shifted peaks. The flux is normalized to match the emission peak.

An interesting difference found between the spectra of SN 2014ab and SN 2010jl is in the blended emission line of the O Iλ8446 and Ca II NIR triplet λλ8498, 8542, 8662 features. These lines appear to be brighter in SN 2014ab than in SN 2010jl. The Ca II NIR triplet is likely produced within the cold dense shell where the hydrogen emission lines also originate (e.g., Dessart et al. 2016). Since the strength of the hydrogen emission lines are similar in SN 2014ab and SN 2010jl, the difference in the emission strength found here is presumably due to the O I line rather than the Ca II NIR triplet. Indeed, the emission peak is bluer in SN 2014ab. Since the emission from the inner SN ejecta is hidden by the electron scattering and no other O I lines are observed, the O Iλ8446 emission is likely to be excited by Lyβ fluorescence (Bowen 1947; Bhatia & Kastner 1995; Kastner & Bhatia 1995).

5.2. Luminosity evolution

As we have discussed so far, SN 2014ab is similar to SN 2010jl in many respects. Figure 12 shows a comparison between the bolometric LCs of SN 2014ab and SN 2010jl. The bolometric LC of SN 2014ab is obtained by integrating the two blackbody components shown in Sect. 4.2 and the bolometric LC of SN 2010jl is obtained by integrating the photometry in optical (3400 − 9000 Å) and NIR (9000 − 24 000 Å) as given in Fransson et al. (2014). We assumed that the explosion date of SN 2014ab is 100 days before the optical discovery, which provides a good match to the bolometric LC of SN 2010jl. The explosion date, however, is uncertain. Under this assumption, the bolometric LC can be fit with the same power-law function as in SN 2010jl with a different scaling factor (Fransson et al. 2014):

L = 1.2 × 10 44 ( t 1 day ) 0.54 e r g s 1 , $$ \begin{aligned} L=1.2\times 10^{44}\left(\frac{t}{1~\mathrm{day} }\right)^{-0.54}~\mathrm erg~s^{-1} , \end{aligned} $$(1)

thumbnail Fig. 12.

Bolometric LC of SN 2014ab compared to that of SN 2010jl. Here, the explosion epoch of SN 2014ab is assumed to occur one hundred days prior to the optical discovery. The “hot + warm” luminosity is obtained by adding the two blackbody components and the “hot” luminosity is obtained only with the optical blackbody component in SN 2014ab. The LCs of SN 2010jl are from Fransson et al. (2014). The “optical + NIR” luminosity is the luminosity in the wavelength range from 3600 Å to 24 000 Å. The “optical” luminosity wavelength range is from 3600 Å to 9000 Å. We show the power-law bolometric LC fit that is analytically expected for interacting SNe and the total nuclear deposition energy from the decay of 56Ni and 56Co.

which is plotted in Fig. 12.

Assuming that the mass loss from the progenitor is steady, the CSM density structure can be expressed as

ρ CSM ( r ) = M ˙ 4 π V wind r 2 , $$ \begin{aligned} \rho _\mathrm{CSM} (r) = \frac{\dot{M}}{4\pi V_\mathrm{wind} }r^{-2}, \end{aligned} $$(2)

where is the mass-loss rate and Vwind is the wind velocity. The bolometric luminosity from the shock interaction with the shock velocity of Vshock can be expressed as

L = ϵ 2 M ˙ V wind V shock 3 , $$ \begin{aligned} L = \frac{\epsilon }{2}\frac{\dot{M}}{V_\mathrm{wind} }V_\mathrm{shock} ^3, \end{aligned} $$(3)

where ϵ is the conversion efficiency from kinetic energy to radiation energy (e.g., Moriya et al. 2013a). The mass-loss rate, therefore, can be estimated as

M ˙ = 0.3 ( ϵ 0.3 ) 1 ( L 10 43 e r g s 1 ) ( V shock 3000 k m s 1 ) 3 V wind 80 k m s 1 M y r 1 . $$ \begin{aligned} \dot{M}=0.3 \left(\frac{\epsilon }{0.3}\right)^{-1}\left(\frac{L}{10^{43}~\mathrm erg~s^{-1} }\right)\left(\frac{V_\mathrm{shock} }{3000\,\mathrm km\,s^{-1} } \right)^{-3}\frac{V_\mathrm{wind} }{80~\mathrm km~s^{-1} }\ {M}_\odot ~\mathrm yr^{-1} . \end{aligned} $$(4)

The bolometric luminosity of SN 2014ab is around 1043 erg s−1 at around 100 − 200 days (Fig. 12). Although Vshock is not constrained in SN 2014ab, well-observed SNe IIn including SN 2010jl (Fransson et al. 2014) and SN 2013L (Taddia et al. 2020) typically show Vshock ≃ 3000 km s−1. The conversion efficiency, ϵ, is another uncertain parameter. Numerical simulations often indicate that this parameter is on the order of 0.1 (van Marle et al. 2010; Moriya et al. 2013b). It has also been suggested to be time dependent (Tsuna et al. 2019; Takei & Shigeyama 2020). For simplicity, we assume that ϵ is constant and adopt the value ϵ ≃ 0.3, which is on the order of 0.1. These assumptions are reasonable enough for our order-of-magnitude estimate. Assuming the wind velocity estimated obtained from the narrow Hα P Cygni profile (≃80 km s−1, Sect. 4.3), Eq. (4) then implies a mass-loss rate of SN 2014ab’s progenitor of 0.3 M yr−1. This estimated mass-loss rate is similar to that found by Fransson et al. (2014) for SN 2010jl (∼0.1 M yr−1). We note that the mass-loss rate estimated here should be regarded as a very rough estimate, given the uncertain shock velocity information and ϵ. It is an order-of-magnitude estimate.

A significant difference between SN 2014ab and SN 2010jl appears at the late phases beginning around 400 days after explosion. The bolometric luminosity (hot + warm in Fig. 12) of SN 2014ab decreases by a factor of ≃7 during the gap in observations between 200 and 400 days, while over the same period, the bolometric luminosity (optical + NIR in Fig. 12) of SN 2010jl decreases by a factor of ≃1.5. One possible reason for this discrepancy could be due to differences in the structure of their CSM. Because the interaction luminosity is proportional to (Eq. (4)), the luminosity decrease by a factor of 0.15 corresponds to the decrease by a factor of 0.15. Thus, the mass-loss rate of the progenitor of SN 2014ab may have been around 10−2M yr−1 and increased to the order of 0.1 M yr−1 shortly before the explosion, while the progenitor of SN 2010jl maintained the large mass-loss rate for a longer time. The difference in the mass-loss timescale may originate from the difference in the progenitor mass. The duration of extensive mass loss induced by pulsational pair-instability SNe, for example, strongly depends on the progenitor mass (e.g., Woosley 2017).

We note that the late-phase bolometric luminosity decline is consistent with the nuclear decay rate of 56Co (Fig. 12). However, the required initial 56Ni mass to account for the late-phase luminosity is 3 M. If the ejecta have such a large amount of 56Ni, it is likely to have an impact on the LC in the early phase because the total luminosity from the nuclear decay is larger than the observed luminosity over the entire early phase (Fig. 12). We also note that we cannot rule out the possibility that the majority of the radioactive energy in the ejecta is thermalized and diffused with a long diffusion time to avoid the early effect on the luminosity. In such a case, however, the diffusion time in the ejecta needs to be more than 200 days following “Arnett’s rule” (Arnett 1982), which would require an extremely large ejecta mass. It may be argued that the NIR luminosity originates from light echoes and only the optical (hot) component should be regarded as the true bolometric luminosity. If this is the case, the initial 56Ni mass required to account for the luminosity becomes 0.8 M (Fig. 12) and it does not affect the early luminosity much. However, this amount of 56Ni is still much larger than would normally be synthesized in SNe (Anderson 2019; Meza & Anderson 2020).

Finally, we assumed that the dense CSM is spherically symmetric and comes from steady wind when analyzing the bolometric LC and deriving the mass-loss rate in this section. While we do not have clues on the CSM geometry in SN 2014ab, the dense CSM around SN 2010jl is likely to be aspherical (e.g., Patat et al. 2011). We do not know whether or not if the mass loss is eruptive (cf. Moriya et al. 2014). Despite of the possible asphericity and non-steady mass loss, Fransson et al. (2014) estimate the mass-loss rate for SN 2010jl by assuming the spherical symmetry and steady mass loss for simplicity as we have done in this section. Thus, the estimated mass-loss rate of ∼0.3 M yr−1 should be, again, regarded as a rough estimate.

5.3. Origin of dust

Although SN 2014ab and SN 2010jl share many common characteristics, one big difference is in the signatures of dust. The IR excess indicates the existence of warm dust in both cases. However, while the peak of the broad emission lines observed in SN 2010jl shifts towards shorter wavelengths over time (Smith et al. 2011; Gall et al. 2014; Fransson et al. 2014), this is not the case in SN 2014ab (Fig. 8). A shift of the broad emission peak to shorter (bluer) wavelengths is a signature of ongoing dust formation (cf. Kozasa et al. 2009). The peak of the broad hydrogen emission lines in SN 2010jl also shows the wavelength dependence expected from the dust absorption, that is, the peak of the emission lines with shorter wavelengths shifts more to the blue (Maeda et al. 2013; Gall et al. 2014). As demonstrated in Fig. 11, this phenomenon is not observed in SN 2014ab. The overall blue-shift of the broad emission lines can be explained by the acceleration of the unshocked CSM by precursor photons as suggested by Fransson et al. (2014) for SN 2010jl (see also Moriya et al. 2011; Chevalier & Irwin 2011). It has also been suggested that blue-shift is affected by the radiation transfer effect in addition to the pre-shock acceleration (Dessart et al. 2015). SN 2014ab can be a good object to investigate the blue-shift caused by the hydrodynamic and radiation transfer effects because it may be less affected by dust as we discuss below.

The fact that we do not see clear observational signatures of ongoing dust formation in SN 2014ab may indicate that the dust found to be associated with it is not formed within the cool, dense shell where the emission lines are formed during our early-phase observations. The BBwarm component does not exist at early phases. At late phases with T warm BB $ T^{BB}_{\mathrm{warm}} $ estimated to be ∼2000 K, we find that this temperature is too high to originate from the inner SN ejecta located below the cold dense shell (e.g., Wooden et al. 1993; Nozawa et al. 2003). Indeed, T warm BB 2000 K $ T^{\mathrm{BB}}_{\mathrm{warm}}\sim 2000\,\mathrm{K} $ is close to the dust evaporation temperature. Therefore, the pre-existing dust in the dense CSM may be heated by the SN shock or radiation, which drives an increase in the MIR luminosity (Fox et al. 2010). Interestingly, the MIR luminosity keeps decreasing after the sudden MIR luminosity increase that occurs between 200 days and 300 days (Fig. 2). In the case of KISS15s, shown in the same figure, the MIR luminosity keeps increasing with the increasing mass of the presumably newly formed dust. Although we cannot rule out the possibility that dust is newly formed in the late phases, the MIR luminosity decline in SN 2014ab is probably indicative of little to no ongoing dust formation. The decline in the MIR luminosity could be explained by the destruction of the pre-existing dust by the SN shock wave as the ejecta expands over time.

In summary, dust formation in SN 2014ab appears to be inefficient and the observed emission seems to be associated with pre-existing dust formed in the dense CSM prior to the final demise of the progenitor star. SN 2014ab indicates that the origin of dust emission observed in SNe IIn may not always be due to newly formed dust but rather sometimes could be due to pre-existing dust. The exact reason why dust formation in SN 2014ab is inefficient is not clear and warrants further investigation. Finally, we note that we cannot exclude the possibility that newly-formed dust may contribute to the late-phase MIR radiation. Indeed, both newly formed and pre-existing dust can appear at the same time, as has been suggested for SN 2006jd (Stritzinger et al. 2012) and other SNe IIn (e.g., Gerardy et al. 2002; Fox et al. 2009).

6. Conclusions

We have reported the optical and NIR observations of Type IIn SN 2014ab by CSP-II, as well as the serendipitous MIR observations by WISE. Although the CSP-II observations started immediately after the optical discovery, the MIR observation showed that SN 2014ab had already been bright at 56 days before the optical discovery and the explosion date of SN 2014ab is uncertain. Nonetheless, we find that SN 2014ab shows similarities with SN 2010jl, one of the best-observed slowly declining SNe IIn. SN 2014ab is fainter than SN 2010jl by a factor of ∼2, but its LC evolution is similar to that of SN 2010jl. The spectroscopic properties are almost identical to those of SN 2010jl. Based on the bolometric LC and assuming a spherically symmetric dense CSM and steady mass loss, the mass-loss rate of the progenitor immediately before the explosion is estimated to be on the order of 0.1 M yr−1.

The MIR observation of SN 2014ab revealed that dust emission with the blackbody temperature of around 2000 K (close to the dust evaporation temperature) exists from around 350 days after the WISE first detection. The MIR luminosity increases up to 350 days and then declines afterwards. Although we observe the warm dust emission at late phases, no clear observational signatures of dust formation within the cool dense shell such as the blue-shift of the broad line emission peaks with time and the wavelength-dependent shift of the emission peak are observed at early phases. These features were observed in SN 2010jl. In addition, KISS15s, which is another SN 2010jl-like SN IIn with serendipitous MIR observations by WISE, showed a gradual increase in the MIR luminosity, rather than the continuous decline found in SN 2014ab. Based on these differences, we presume that the dust emission observed in SN 2014ab originates from pre-existing dust located within the dense circumstellar matter heated by the SN shock or radiation, and the dust formation at the cool dense shell is much less efficient in SN 2014ab as compared with SN 2010jl; however, the possibility of some contribution from newly-formed dust in the late-phase radiation is not excluded. The continuous MIR luminosity decline of SN 2014ab may also indicate that the heated dust is continuously destroyed by the shock. Our study of SN 2014ab reveals a diversity of plausible origins for dust emission and dust formation in SNe IIn.


Acknowledgments

We thank the anonymous referee for the constructive comments that improved the quality of the manuscript. We thank Morgan Fraser and Takaya Nozawa for important discussions, and Ryan Foley, Howie Marion and Bob Kirshner for obtaining FIRE spectra of SN 2010jl. T.J.M. is supported by the Grants-in-Aid for Scientific Research of the Japan Society for the Promotion of Science (JP17H02864, JP18K13585, JP20H00174). M.S. and S.H. are supported in part by generous grants (project numbers 13261 and 28021) from VILLUM FONDEN. M.S. and F.T. are also funded by a project grant (8021-00170B) from the Independent Research Fund Denmark. N.B.S. acknowledges support from the NSF through grant AST-1613455, and through the Texas A&M University Mitchell/Heep/Munnerlyn Chair in Observational Astronomy. C.G. was supported by a VILLUM FONDEN Young Investigator grant (project number 25501). This work was supported by a VILLUM FONDEN Investigator grant to J.H. (project number 16599). L.G. was funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 839090. This work has been partially supported by the Spanish grant PGC2018-095317-B-C21 within the European Funds for Regional Development (FEDER). This work is supported by the Japan Society for the Promotion of Science Open Partnership Bilateral Joint Research Project between Japan and Chile (JPJSBP120209937). The CSP-II has been funded by the USA’s NSF under grants AST-0306969, AST-0607438, AST-1008343, AST-1613426, AST-1613455, and AST-1613472, and in part by a Sapere Aude Level 2 grant funded by the Danish Agency for Science and Technology and Innovation (PI M.S.). NTT (+ ELFOSC) spectrum was obtained by ESO program 191.D-0935. VLT (+ X-shooter) spectra were obtained under ESO program 092.D-0645, and the data were retrieved from the ESO archive. This work made use of the Open Supernova Catalog (Guillochon et al. 2017; https://sne.space/) and WISeREP (Yaron & Gal-Yam 2012; https://wiserep.weizmann.ac.il). This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

References

  1. Anderson, J. P. 2019, A&A, 628, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Anderson, J. P., Habergham, S. M., James, P. A., & Hamuy, M. 2012, MNRAS, 424, 1372 [NASA ADS] [CrossRef] [Google Scholar]
  3. Andrews, J. E., Clayton, G. C., Wesson, R., et al. 2011, AJ, 142, 45 [NASA ADS] [CrossRef] [Google Scholar]
  4. Arnett, W. D. 1982, ApJ, 253, 785 [NASA ADS] [CrossRef] [Google Scholar]
  5. Berger, E., Soderberg, A. M., Chevalier, R. A., et al. 2009, ApJ, 699, 1850 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bevan, A. M., Krafton, K., Wesson, R., et al. 2020, ApJ, 894, 111 [CrossRef] [Google Scholar]
  7. Bhatia, A. K., & Kastner, S. O. 1995, ApJS, 96, 325 [NASA ADS] [CrossRef] [Google Scholar]
  8. Boldt, L. N., Stritzinger, M. D., Burns, C., et al. 2014, PASP, 126, 324 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bond, H. E., Bedin, L. R., Bonanos, A. Z., et al. 2009, ApJ, 695, L154 [NASA ADS] [CrossRef] [Google Scholar]
  10. Borish, H. J., Huang, C., Chevalier, R. A., et al. 2015, ApJ, 801, 7 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bowen, I. S. 1947, PASP, 59, 196 [NASA ADS] [CrossRef] [Google Scholar]
  12. Burns, C. R., Parent, E., Phillips, M. M., et al. 2018, ApJ, 869, 56 [NASA ADS] [CrossRef] [Google Scholar]
  13. Buzzoni, B., Delabre, B., Dekker, H., et al. 1984, The Messenger, 38, 9 [NASA ADS] [Google Scholar]
  14. Chandra, P., Chevalier, R. A., Irwin, C. M., et al. 2012, ApJ, 750, L2 [NASA ADS] [CrossRef] [Google Scholar]
  15. Chandra, P., Chevalier, R. A., Chugai, N., Fransson, C., & Soderberg, A. M. 2015, ApJ, 810, 32 [NASA ADS] [CrossRef] [Google Scholar]
  16. Chatzopoulos, E., Wheeler, J. C., Vinko, J., Horvath, Z. L., & Nagy, A. 2013, ApJ, 773, 76 [NASA ADS] [CrossRef] [Google Scholar]
  17. Chevalier, R. A., & Irwin, C. M. 2011, ApJ, 729, L6 [Google Scholar]
  18. Chugai, N. N. 2001, MNRAS, 326, 1448 [NASA ADS] [CrossRef] [Google Scholar]
  19. Chugai, N. N. 2018, MNRAS, 481, 3643 [CrossRef] [Google Scholar]
  20. Chugai, N. N., & Danziger, I. J. 1994, MNRAS, 268, 173 [NASA ADS] [Google Scholar]
  21. Chugai, N. N., Blinnikov, S. I., Cumming, R. J., et al. 2004, MNRAS, 352, 1213 [NASA ADS] [CrossRef] [Google Scholar]
  22. Dessart, L., Audit, E., & Hillier, D. J. 2015, MNRAS, 449, 4304 [NASA ADS] [CrossRef] [Google Scholar]
  23. Dessart, L., Hillier, D. J., Audit, E., Livne, E., & Waldman, R. 2016, MNRAS, 458, 2094 [NASA ADS] [CrossRef] [Google Scholar]
  24. Fassia, A., Meikle, W. P. S., Vacca, W. D., et al. 2000, MNRAS, 318, 1093 [NASA ADS] [CrossRef] [Google Scholar]
  25. Fassia, A., Meikle, W. P. S., Chugai, N., et al. 2001, MNRAS, 325, 907 [NASA ADS] [CrossRef] [Google Scholar]
  26. Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
  27. Fox, O., Skrutskie, M. F., Chevalier, R. A., et al. 2009, ApJ, 691, 650 [NASA ADS] [CrossRef] [Google Scholar]
  28. Fox, O. D., Chevalier, R. A., Dwek, E., et al. 2010, ApJ, 725, 1768 [NASA ADS] [CrossRef] [Google Scholar]
  29. Fox, O. D., Chevalier, R. A., Skrutskie, M. F., et al. 2011, ApJ, 741, 7 [NASA ADS] [CrossRef] [Google Scholar]
  30. Fox, O. D., Filippenko, A. V., Skrutskie, M. F., et al. 2013, AJ, 146, 2 [NASA ADS] [CrossRef] [Google Scholar]
  31. Fransson, C., Ergon, M., Challis, P. J., et al. 2014, ApJ, 797, 118 [NASA ADS] [CrossRef] [Google Scholar]
  32. Fraser, M., Magee, M., Kotak, R., et al. 2013, ApJ, 779, L8 [NASA ADS] [CrossRef] [Google Scholar]
  33. Fraser, M., Blagorodnova, N., Walton, N., et al. 2014, ATel, 5968, 1 [Google Scholar]
  34. Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gal-Yam, A., & Leonard, D. C. 2009, Nature, 458, 865 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  36. Galbany, L., Anderson, J. P., Sánchez, S. F., et al. 2018, ApJ, 855, 107 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gall, C., Hjorth, J., Watson, D., et al. 2014, Nature, 511, 326 [NASA ADS] [CrossRef] [Google Scholar]
  38. Gerardy, C. L., Fesen, R. A., Nomoto, K., et al. 2002, ApJ, 575, 1007 [NASA ADS] [CrossRef] [Google Scholar]
  39. Guillochon, J., Parrent, J., Kelley, L. Z., & Margutti, R. 2017, ApJ, 835, 64 [NASA ADS] [CrossRef] [Google Scholar]
  40. Habergham, S. M., Anderson, J. P., James, P. A., & Lyman, J. D. 2014, MNRAS, 441, 2230 [NASA ADS] [CrossRef] [Google Scholar]
  41. Hamuy, M., Folatelli, G., Morrell, N. I., et al. 2006, PASP, 118, 2 [NASA ADS] [CrossRef] [Google Scholar]
  42. Howerton, S., Drake, A. J., Djorgovski, S. G., et al. 2014, Cent. Bureau Electron. Tel., 3826, 1 [Google Scholar]
  43. Hsiao, E. Y., Phillips, M. M., Marion, G. H., et al. 2019, PASP, 131, 014002 [NASA ADS] [CrossRef] [Google Scholar]
  44. Huang, C., & Chevalier, R. A. 2018, MNRAS, 475, 1261 [NASA ADS] [CrossRef] [Google Scholar]
  45. Jencson, J. E., Prieto, J. L., Kochanek, C. S., et al. 2016, MNRAS, 456, 2622 [NASA ADS] [CrossRef] [Google Scholar]
  46. Jiang, N., Wang, T., Mou, G., et al. 2019, ApJ, 871, 15 [CrossRef] [Google Scholar]
  47. Kankare, E., Ergon, M., Bufano, F., et al. 2012, MNRAS, 424, 855 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kastner, S. O., & Bhatia, A. K. 1995, ApJ, 439, 346 [NASA ADS] [CrossRef] [Google Scholar]
  49. Kausch, W., Noll, S., Smette, A., et al. 2015, A&A, 576, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Kiewe, M., Gal-Yam, A., Arcavi, I., et al. 2012, ApJ, 744, 10 [NASA ADS] [CrossRef] [Google Scholar]
  51. Kokubo, M., Mitsuda, K., Morokuma, T., et al. 2019, ApJ, 872, 135 [CrossRef] [Google Scholar]
  52. Kozasa, T., Nozawa, T., Tominaga, N., et al. 2009, in Dust in Supernovae: Formation and Evolution, eds. T. Henning, E. Grün, J. Steinacker, et al., ASP Conf. Ser., 414, 43 [Google Scholar]
  53. Krisciunas, K., Contreras, C., Burns, C. R., et al. 2017, AJ, 154, 211 [NASA ADS] [CrossRef] [Google Scholar]
  54. Maeda, K., Nozawa, T., Sahu, D. K., et al. 2013, ApJ, 776, 5 [CrossRef] [Google Scholar]
  55. Mainzer, A., Bauer, J., Cutri, R. M., et al. 2014, ApJ, 792, 30 [NASA ADS] [CrossRef] [Google Scholar]
  56. Mauerhan, J. C., Smith, N., Silverman, J. M., et al. 2013, MNRAS, 431, 2599 [NASA ADS] [CrossRef] [Google Scholar]
  57. Meza, N., & Anderson, J. P. 2020, A&A, in press, https://doi.org/10.1051/0004-6361/201937113 [Google Scholar]
  58. Miller, A. A., Silverman, J. M., Butler, N. R., et al. 2010, MNRAS, 404, 305 [NASA ADS] [Google Scholar]
  59. Moriya, T., Tominaga, N., Blinnikov, S. I., Baklanov, P. V., & Sorokina, E. I. 2011, MNRAS, 415, 199 [NASA ADS] [CrossRef] [Google Scholar]
  60. Moriya, T. J., Maeda, K., Taddia, F., et al. 2013a, MNRAS, 435, 1520 [NASA ADS] [CrossRef] [Google Scholar]
  61. Moriya, T. J., Blinnikov, S. I., Tominaga, N., et al. 2013b, MNRAS, 428, 1020 [NASA ADS] [CrossRef] [Google Scholar]
  62. Moriya, T. J., Maeda, K., Taddia, F., et al. 2014, MNRAS, 439, 2917 [NASA ADS] [CrossRef] [Google Scholar]
  63. Moriya, T. J., Tanaka, M., Morokuma, T., et al. 2019, ApJ, 882, 70 [CrossRef] [Google Scholar]
  64. Nozawa, T., Kozasa, T., Umeda, H., Maeda, K., & Nomoto, K. 2003, ApJ, 598, 785 [NASA ADS] [CrossRef] [Google Scholar]
  65. Nyholm, A., Sollerman, J., Taddia, F., et al. 2017, A&A, 605, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Nyholm, A., Sollerman, J., Tartaglia, L., et al. 2020, A&A, 637, A73 [CrossRef] [EDP Sciences] [Google Scholar]
  67. Ofek, E. O., Sullivan, M., Cenko, S. B., et al. 2013, Nature, 494, 65 [NASA ADS] [CrossRef] [Google Scholar]
  68. Ofek, E. O., Sullivan, M., Shaviv, N. J., et al. 2014a, ApJ, 789, 104 [NASA ADS] [CrossRef] [Google Scholar]
  69. Ofek, E. O., Arcavi, I., Tal, D., et al. 2014b, ApJ, 788, 154 [NASA ADS] [CrossRef] [Google Scholar]
  70. Ofek, E. O., Zoglauer, A., Boggs, S. E., et al. 2014c, ApJ, 781, 42 [NASA ADS] [CrossRef] [Google Scholar]
  71. Ohsawa, R., Sakon, I., Onaka, T., et al. 2010, ApJ, 718, 1456 [CrossRef] [Google Scholar]
  72. Patat, F., Taubenberger, S., Benetti, S., Pastorello, A., & Harutyunyan, A. 2011, A&A, 527, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Persson, S. E., Murphy, D. C., Krzeminski, W., Roth, M., & Rieke, M. J. 1998, AJ, 116, 2475 [NASA ADS] [CrossRef] [Google Scholar]
  74. Persson, S. E., Murphy, D. C., Smee, S., et al. 2013, PASP, 125, 654 [NASA ADS] [CrossRef] [Google Scholar]
  75. Phillips, M. M., Contreras, C., Hsiao, E. Y., et al. 2019, PASP, 131, 014001 [NASA ADS] [CrossRef] [Google Scholar]
  76. Pozzo, M., Meikle, W. P. S., Fassia, A., et al. 2004, MNRAS, 352, 457 [NASA ADS] [CrossRef] [Google Scholar]
  77. Prieto, J. L., Kistler, M. D., Thompson, T. A., et al. 2008, ApJ, 681, L9 [NASA ADS] [CrossRef] [Google Scholar]
  78. Sarangi, A., Dwek, E., & Arendt, R. G. 2018, ApJ, 859, 66 [CrossRef] [Google Scholar]
  79. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  80. Schlegel, E. M. 1990, MNRAS, 244, 269 [NASA ADS] [Google Scholar]
  81. Selsing, J., Malesani, D., Goldoni, P., et al. 2019, A&A, 623, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Simcoe, R. A., Burgasser, A. J., Schechter, P. L., et al. 2013, PASP, 125, 270 [NASA ADS] [CrossRef] [Google Scholar]
  83. Smartt, S. J., Valenti, S., Fraser, M., et al. 2015, A&A, 579, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Smette, A., Sana, H., Noll, S., et al. 2015, A&A, 576, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Smith, N., Li, W., Foley, R. J., et al. 2007, ApJ, 666, 1116 [NASA ADS] [CrossRef] [Google Scholar]
  86. Smith, N., Chornock, R., Li, W., et al. 2008, ApJ, 686, 467 [NASA ADS] [CrossRef] [Google Scholar]
  87. Smith, N., Li, W., Miller, A. A., et al. 2011, ApJ, 732, 63 [NASA ADS] [CrossRef] [Google Scholar]
  88. Sollerman, J., Cumming, R. J., & Lundqvist, P. 1998, ApJ, 493, 933 [NASA ADS] [CrossRef] [Google Scholar]
  89. Stathakis, R. A., & Sadler, E. M. 1991, MNRAS, 250, 786 [NASA ADS] [Google Scholar]
  90. Stoll, R., Prieto, J. L., Stanek, K. Z., et al. 2011, ApJ, 730, 34 [NASA ADS] [CrossRef] [Google Scholar]
  91. Stritzinger, M., Taddia, F., Fransson, C., et al. 2012, ApJ, 756, 173 [NASA ADS] [CrossRef] [Google Scholar]
  92. Suzuki, A., Moriya, T. J., & Takiwaki, T. 2019, ApJ, 887, 249 [CrossRef] [Google Scholar]
  93. Szalai, T., Zsíros, S., Fox, O. D., Pejcha, O., & Müller, T. 2019, ApJS, 241, 38 [NASA ADS] [CrossRef] [Google Scholar]
  94. Taddia, F., Stritzinger, M. D., Sollerman, J., et al. 2013, A&A, 555, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Taddia, F., Sollerman, J., Fremling, C., et al. 2015, A&A, 580, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Taddia, F., Stritzinger, M. D., Fransson, C., et al. 2020, A&A, 638, A92 [CrossRef] [EDP Sciences] [Google Scholar]
  97. Takei, Y., & Shigeyama, T. 2020, PASJ, 72, 67 [CrossRef] [Google Scholar]
  98. Tsuna, D., Kashiyama, K., & Shigeyama, T. 2019, ApJ, 884, 87 [NASA ADS] [CrossRef] [Google Scholar]
  99. Turatto, M., Cappellaro, E., Danziger, I. J., et al. 1993, MNRAS, 262, 128 [NASA ADS] [CrossRef] [Google Scholar]
  100. van Marle, A. J., Smith, N., Owocki, S. P., & van Veelen, B. 2010, MNRAS, 407, 2305 [NASA ADS] [CrossRef] [Google Scholar]
  101. Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Williams, B. J., & Fox, O. D. 2015, ApJ, 808, L22 [CrossRef] [Google Scholar]
  103. Wooden, D. H., Rank, D. M., Bregman, J. D., et al. 1993, ApJS, 88, 477 [NASA ADS] [CrossRef] [Google Scholar]
  104. Woosley, S. E. 2017, ApJ, 836, 244 [NASA ADS] [CrossRef] [Google Scholar]
  105. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
  106. Yaron, O., & Gal-Yam, A. 2012, PASP, 124, 668 [NASA ADS] [CrossRef] [Google Scholar]
  107. Zhang, T., Wang, X., Wu, C., et al. 2012, AJ, 144, 131 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Photometric data

Photometric data of SN 2014ab used in our analysis are provided in this appendix. This includes two-channel WISE photometry of SN 2014ab listed in Table A.1, optical and NIR local sequences of SN 2014ab in Tables A.2 and A.3, respectively, and optical/NIR photometry of SN 2014ab in Table A.4.

Table A.1.

WISE photometry of SN 2014ab in the Vega system.

Table A.2.

Optical photometry of the local sequence for SN 2014ab in the “standard” system.

Table A.3.

NIR photometry of the local sequences for SN 2014ab in the “standard” system.

Table A.4.

Photometry of SN 2014ab in the “natural” system.

Appendix B: CSP-II spectra of SN 2010jl

In Fig. B.1, we present unpublished NIR spectra of SN 2010jl, which is shown in this paper to have common features with SN 2014ab. A journal summarizing the spectra obtained by the CSP-II of SN 2010jl is provided in Table B.1. The spectroscopic data was reduced following standard techniques. The visual-wavelength spectrum of SN 2010jl is shown in Fig. 10, while five NIR spectra obtained with the Baade (+ FIRE) telescope are plotted in Fig. B.1, and the rest-frame axes assumes the host-galaxy redshift of 0.0107. Optical spectra of SN 2010jl at similar epochs have already been presented and analysed in the literature, but the NIR spectra cover late epochs than those published in the literature to our knowledge. These data are presented here for the benefit of the community.

thumbnail Fig. B.1.

NIR spectra of SN 2010jl obtained by the CSP-II. The epochs are in the rest frame relative to the first detection which occurreed on JD = 2455479.0. The number in the parenthesis next to the epochs indicate the constant value added to each spectrum for presentation purposes and are in the unit of 10−14 erg s−1 cm−2 Å−1.

Table B.1.

Journal of spectroscopic observations of SN 2010jl.

All Tables

Table 1.

Journal of spectroscopic observations of SN 2014ab.

Table A.1.

WISE photometry of SN 2014ab in the Vega system.

Table A.2.

Optical photometry of the local sequence for SN 2014ab in the “standard” system.

Table A.3.

NIR photometry of the local sequences for SN 2014ab in the “standard” system.

Table A.4.

Photometry of SN 2014ab in the “natural” system.

Table B.1.

Journal of spectroscopic observations of SN 2010jl.

All Figures

thumbnail Fig. 1.

Finding chart constructed using a single r-band image obtained on 03 March 2014 with the Swope 1.0 m telescope. The position of SN 2014ab is indicated with a dotted circle and the approximate area contained within the circle is shown in the upper right inset Optical and NIR local sequence stars are indicated with yellow and red circles, respectively.

In the text
thumbnail Fig. 2.

Mid-infrared light curves of SN 2014ab obtained with the WISE satellite. The open symbols show the raw data and the filled symbols show co-added magnitudes of the nearby epochs. We also present the photometry of KISS15s, which is another SN IIn with serendipitous WISE observations (Kokubo et al. 2019). The KISS15s photometry is plotted with the right and top axes. The time origin in KISS15s is its optical discovery.

In the text
thumbnail Fig. 3.

Multi-color, optical (ugriBV) and NIR (YJH) light curves of SN 2014ab obtained by the CSP-II. The data have been corrected for Galactic extinction. Vertical lines indicate epochs that optical (green), NIR (red), or optical+NIR (blue) spectroscopic observations were obtained.

In the text
thumbnail Fig. 4.

Observed B − V color curve of SN 2014ab compared with those of three different SNe IIn subtypes. The reddening is corrected based on the values reported in the literature. The explosion date, i.e., day zero, for SN 2014ab is arbitrarily set to be 100 days before the optical discovery, or 46 days before the WISE first detection in the rest frame. The comparison sample including objects of the slowly-declining 1988Z-like subtype (solid lines) including SN 2010jl (Fransson et al. 2014), SN 2006tf (Smith et al. 2008), and SN 1988Z (Turatto et al. 1993), as well as the 1994W-like plateau subtype (dashed lines) including SN 1994W (Sollerman et al. 1998), SN 2009kn (Kankare et al. 2012), and SN 2011ht (Mauerhan et al. 2013), and finally, the prototypical rapidly-declining SN 1998S (dot-dashed line, Fassia et al. 2000).

In the text
thumbnail Fig. 5.

Spectral energy distributions (SEDs) of SN 2014ab. The first 16 quadrants from the top contain SEDs constructed from flux points covering optical, NIR, and MIR wavelengths and extend from +56 d to +140 d relative to the time of first WISE detection. Each of these SEDs are fit with a single black body (BB) function, where the optical flux is dominating. We excluded u and r band from the fit due to line blanketing and the strong Hα emission line, respectively. The SEDs plotted in the other quadrants are constructed from optical and MIR flux points and extend from +360 d to +485 d relative to first WISE detection. The late-phase SEDs clearly exhibit an excess of flux relative to a single BB fit to the optical flux points. We therefore fit the late SEDs with a two-component BB function. In the last quadrant the optical BB fit is not reliable as only g and i are available.

In the text
thumbnail Fig. 6.

Bolometric LC (top) and the BB-temperature (middle) and BB-radius (bottom) profiles of SN 2014ab, computed from fitting single and double BB functions to the SEDs of SN 2014ab.

In the text
thumbnail Fig. 7.

Optical through NIR wavelength spectroscopy of SN 2014ab, extending from +54 d to 183 d (restframe days) past the first WISE detection. The numbers in the parentheses indicate the flux offset adopted for the spectra in the unit of 10−16 erg s−1 cm−2 Å−1. The shaded wavelength regions indicate areas of strong telluric absorption. The last three spectra are not calibrated by photometry because of the lack in the photometric information.

In the text
thumbnail Fig. 8.

Time-evolution of the narrow- (top panel) and broad-component (middle panel) Hα profiles, and the broad-component (bottom panel) of the Pa-β profile, plotted in velocity space. Over-plotted as dashed cyan lines are Lorentzian functions with the FWHM velocity of 3500 km s−1 and the peak shifted by −500 km s−1. The redshift is determined by the average of the narrow emission peak. Both broad components are overall symmetric and do not exhibit significant evolution over the duration of our followup observations.

In the text
thumbnail Fig. 9.

Absolute V-band LC of SN 2014ab compared with well-observed examples of the three subtypes of SN IIn discussed by Taddia et al. (2013). This includes the long-lasting 1988Z-like objects (circles) SN 1988Z (Turatto et al. 1993), SN 2006tf (Smith et al. 2008), SN 2010jl (Stoll et al. 2011; Zhang et al. 2012; Maeda et al. 2013; Fransson et al. 2014), and KISS15s (Kokubo et al. 2019), the 1994W-like objects (squares) including SN 1994W (Chugai et al. 2004), SN 2009kn (Kankare et al. 2012), and SN 2011ht (Mauerhan et al. 2013), which show a LC pleateau, and finally, the prototypical, rapidly-declining SN 1998S (triangles, Fassia et al. 2000).

In the text
thumbnail Fig. 10.

Comparison of our +67 d visual-wavelength (top) and NIR-wavelength (bottom) spectrum of SN 2014ab with other SN IIn subtypes. Comparison objects are the slowly declining SN 2010jl (Borish et al. 2015), the 1994W-like SN 2009kn (Kankare et al. 2012), and the rapid-declining SN 1998S (Fassia et al. 2001). Prominent ions are indicated and labeled. Regions suffering from strong telluric lines are shaded. Note that the NIR spectrum of SN 1998S is composed of data taken on two different but close in time epochs as indicated. Note that the previously unpublished visual-wavelength spectrum of SN 2010jl was obtained by the CSP-II (see Appendix B for details).

In the text
thumbnail Fig. 11.

Peak of the hydrogen emission lines. We use the optical spectrum at +67 d and the NIR spectrum at +69 d. The peak of the broad component does not show a wavelength-dependent shift expected if dust is formed at the emission line forming region, i.e., the emission lines at the shorter wavelengths to have more blue-shifted peaks. The flux is normalized to match the emission peak.

In the text
thumbnail Fig. 12.

Bolometric LC of SN 2014ab compared to that of SN 2010jl. Here, the explosion epoch of SN 2014ab is assumed to occur one hundred days prior to the optical discovery. The “hot + warm” luminosity is obtained by adding the two blackbody components and the “hot” luminosity is obtained only with the optical blackbody component in SN 2014ab. The LCs of SN 2010jl are from Fransson et al. (2014). The “optical + NIR” luminosity is the luminosity in the wavelength range from 3600 Å to 24 000 Å. The “optical” luminosity wavelength range is from 3600 Å to 9000 Å. We show the power-law bolometric LC fit that is analytically expected for interacting SNe and the total nuclear deposition energy from the decay of 56Ni and 56Co.

In the text
thumbnail Fig. B.1.

NIR spectra of SN 2010jl obtained by the CSP-II. The epochs are in the rest frame relative to the first detection which occurreed on JD = 2455479.0. The number in the parenthesis next to the epochs indicate the constant value added to each spectrum for presentation purposes and are in the unit of 10−14 erg s−1 cm−2 Å−1.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.