Free Access
Issue
A&A
Volume 630, October 2019
Article Number A76
Number of page(s) 32
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201935537
Published online 23 September 2019

© ESO 2019

1. Introduction

Type Ia supernovae (SNe Ia) are among the brightest transient objects in the Universe, typically reaching peak absolute magnitudes of ≲ − 19 mag. They are also very homogeneous and have peak brightnesses varying ∼1−2 mag (e.g., Hamuy et al. 1996). However, through empirical calibrations, such as the relations between absolute peak magnitude and the post-maximum brightness decline (e.g., Δm15(B), the parameter accounting for magnitude decay after 15 days post-maximum in the B band; Phillips 1993) and peak color (Riess et al. 1996; Tripp 1998), this scatter around peak brightness is reduced to ∼0.1 mag. Owing to their standardizability and because they can be discovered from the ground up to high redshifts (z ≲ 1.0), SNe Ia have been widely used in recent decades to put constraints on cosmological parameters and the Hubble constant (e.g., Howell 2011; Burns et al. 2018), they were the key to the discovery of the accelerated expansion of the Universe (Riess et al. 1998; Perlmutter et al. 1999), and remain one of the key probes of dark energy (e.g., Betoule et al. 2014; DES Collaboration 2019).

However, their success as cosmological probes is not matched by our understanding of their progenitor systems, the details on how they explode, and the underlying physical interpretation of the known empirical relations that enable these objects to provide precise cosmological distances. Further improvement requires a much better understanding of the physical model of the explosion, observational constraints on which kind of progenitors can produce SNe Ia, a better control of reddening, and a complete characterization of peculiar (non-standardizable) types.

The fact that SNe Ia are discovered in all types of galaxies, including passive galaxies with no signs of star formation, points to old stars being the most probable progenitors (> 100s Myr; Mannucci et al. 2006). At the same time, their uniform brightness suggests that their progenitor systems form a homogeneous family. The consensus is that a SN Ia results from a degenerate (held by electron degeneracy pressure) carbon-oxygen (C-O) white dwarf (WD) undergoing a thermonuclear runaway (Hoyle & Fowler 1960), and that SNe Ia originate from close binary systems. However, the nature of the system and the explosion mechanism is still under debate. Potential progenitor systems may either consist of two WDs, called a double degenerate (DD) system, or a single WD with a non-degenerate companion, called a single degenerate (SD) system. Candidate progenitor systems have been inferred for both SD and DD systems. Although the lack of donor stars surviving the explosion may favor DD systems (Tucker et al. 2019), other examples would favor SD progenitors, such as the potential He donor star identified for the well-observed 2002cx-like type Ia SN 2012Z (McCully et al. 2014a; Stritzinger et al. 2015) and the massive asymptotic giant branch star inferred for SN 2002ic (Hamuy et al. 2003).

Within these progenitor scenarios four classes of explosion scenarios are most favored: (1) the dynamical merging of two C-O WDs in a binary system after losing angular momentum via gravitational radiation, in which the thermonuclear explosion is triggered by the heat of the merging process (García-Berro & Lorén-Aguilar 2017); (2) the explosion of a C-O WD with a mass close to MCh triggered by compressional heating near the WD center, caused by accretion from the donor star, which may be a red giant, a main sequence star of mass less than 7−8 M, a He star, or the tidally disrupted WD from a DD system (Piersanti et al. 2004); (3) the explosion of a sub-MCh C-O WD triggered by detonating a thin surface He layer on the WD (double detonations; Kromer et al. 2010; Shen & Moore 2014); and (4) direct collision of two WDs in tertiary systems, where the third star induces oscillations in the eccentricity of the other two increasing gravitational-wave losses and finishing in a head-on collision (e.g., Thompson 2011; Mazzali et al. 2018). The debate about dynamical mergers, MCh explosions, and double detonations is ongoing, but all channels may contribute to the SNe Ia population.

The connection of these progenitor scenarios and explosion mechanisms with SN Ia subtypes is still not well established. More extreme over- or subluminous SNe Ia may require progenitors deviating from those discussed above. It is of utmost importance to clarify which progenitor configuration and explosion mechanism lead to SNe Ia with particular properties and at what rates for a number of fields including stellar evolution and SN cosmology. The presence of these peculiar over- and subluminous SNe Ia, and other heterogeneous objects that do not follow the empirical relations (Taubenberger 2017), may challenge the picture of all SNe Ia coming from the same family of progenitors.

In particular, subluminous SN 1991bg-like objects seem to break such relations. Photometrically, these objects are characterized as fainter at peak and having faster light curve rise and decline rates, putting them at the faint and fast extreme of the light curve width-luminosity relation (Δm15 ≳ 1.7 mag). The maxima in the near-infrared (NIR) bands occur after the maximum in optical bands, but one of the main peculiarities is the lack of either a shoulder or a secondary maximum in the NIR light curves. The latter can be explained by the recombination of Fe III to Fe II occurring earlier: the Fe II line blanketing absorbs flux in the blue that is re-emitted at longer wavelengths, causing the first and second NIR peaks to merge to form a single, slightly delayed maximum (Kasen & Woosley 2007; Höflich et al. 2017). As a result, their colors are redder at maximum (B − V)Bmax ≳ 0.5 mag and their B − V color curves peak earlier (∼15 days from peak maximum in B-band) compared to normal SNe Ia (∼30 days). Spectroscopically, their main characteristics are the presence of prominent Ti II and strong O Iλ7744 lines, indicating low ionization and temperature, deeper Si IIλ5972 than in normal SNe Ia, and relatively low ejecta velocities with fast decline and no high-velocity features. To first order, the luminosity of a supernova is dependent on the mass of radioactive 56Ni synthesized in the explosion, where less luminous objects produce smaller amounts of 56Ni (Arnett 1982; Pinto & Eastman 2000; Mazzali et al. 2001). There is some evidence, however, that a variation of the total ejecta mass may also be required (Stritzinger et al. 2006; Scalzo et al. 2014; Blondin et al. 2017). Given their low luminosity, 1991bg-like SNe synthesize rather low amounts (≲0.1 M) of radioactive 56Ni (Stritzinger et al. 2006). In addition, most 1991bg-like SNe Ia are found in elliptical galaxies with rather old stellar populations and no ongoing or recent star formation, although the well-observed SN 1999by occurred in the spiral Sa-type galaxy NGC 2841 (Höflich et al. 2002; Garnavich et al. 2004).

In this paper, we present a comprehensive dataset collected for the peculiar SN 2016hnk (see also Sell et al. 2018), which exploded in the Sb-type galaxy MCG -01-06-070. SN 2016hnk showed many similarities with 1991bg-like SNe Ia, including lack of hydrogen features, typical deep Si IIλ6355 absorption, other typical SN Ia features with absorption velocities within 6000 < v < 10 000 km s−1, and very strong calcium lines in the optical. With an ATLAS-c absolute magnitude at discovery of −16.1, it was also reported to be similar to the unusual Ca-rich transient PTF09dav (Pan et al. 2016). Our goal with this complete dataset is to understand the nature of this peculiar object. The paper is structured as follows. In Sect. 2 we describe the discovery and the photometric and spectroscopic follow-up campaign. In Sect. 3, integral field spectroscopy (IFS) of the SN host galaxy from the PMAS/PPak Integral field Supernova hosts COmpilation (PISCO) survey (Galbany et al. 2018) is used to study the SN local environment. Later we present the observational analysis in two sections: the photometric view in Sect. 4 and the spectroscopic view in Sect. 5, including an observational characterization within the 1991bg-like group. We make a comparison with theoretical models from the abundance stratification technique and detailed non-local thermodynamical equilibrium (non-LTE) modified models of SN 1999by in Sect. 6. Finally, we discuss the origin of this object in Sect. 7. Throughout the paper we employed a standard flat ΛCDM cosmology with ΩM = 0.2865 and H0 = 69.32 km s−1 Mpc−1 (Hinshaw et al. 2013).

2. SN 2016hnk observations and reduction

The discovery of SN 2016hnk (aka ATLAS16dpc) was reported by the Asteroid Terrestrial-impact Last Alert System (ATLAS; Tonry et al. 2018) using the ACAM1 instrument in the ATLAS Haleakala telescope on 2016 October 27.47 (JD = 2457688.97) at α = 02h13m1663, δ = −07°39′4080 (J2000; see Fig. 1) and with a magnitude of 17.91 in the cyan-ATLAS filter. The object exploded in the type SB galaxy MCG -01-06-070 of the Morphological Catalogue of Galaxies (Vorontsov-Vel’Yaminov & Arkhipova 1962), 11.5 arcsec east of the center, which at the redshift of the galaxy z = 0.016 corresponds to 3.71 kpc (see inner panel of Fig. 1). The main parameters of the supernova and host galaxy are listed in Table 1.

thumbnail Fig. 1.

Grayscale B-band image obtained with the Las Cumbres 1 m telescope at Cerro Tololo Inter-American Observatory on 2016 November 3. North is up; east is left. SN 2016hnk is denoted with two blue lines at the bottom inner panel, next to its host galaxy MCG -01-06-070. Red circles represent all 35 standard stars used for photometric calibration in the optical and NIR.

Table 1.

Properties of SN 2016hnk and its host galaxy.

The last reported non-detection of the object is from 2016 October 19.53 (JD = 2457681.03) in the orange-ATLAS filter (> 18.66 mag) ∼8 days before the first detection. Forced photometry of prediscovery images, which were differenced with the all-sky reference frame (see Stalder et al. 2017), confirmed this as an actual detection. This pushed the last non-detection to the previous image of the region on 2016 October 10.56 (JD = 2457672.06) in the orange-ATLAS filter, which provides a constraint on the explosion epoch to within nine days between the last non-detection and first detection. Throughout the paper, we assume that the explosion epoch is the mid-point between the last non-detection and the first detection with an uncertainty of half of this time, on JD = 2457676.5 ± 4.5, however, we base SN epochs on the date of peak brightness (see Sect. 4).

SN 2016hnk was initially classified by the NOT (Nordic Optical Telescope) Un-biased Transient Survey (NUTS; Cannizzaro et al. 2016a,b) on 2016 Nov 01.08 UT as a type I SN. However, they noted similarities to SN 1991bg-like spectra. The spectrum was relatively red and had strong lines dominated by the Ca NIR triplet, O Iλ7774, and Si IIλ6355, and a Si II velocity of ∼10 000 km s−1. A few hours later, on November 2 at 20:16 UT, the Public ESO Spectroscopic Survey for Transient Objects (PESSTO; Smartt et al. 2015) confirmed the redshift and provided a definite classification as a SN Ia spectroscopically similar to SN 1991bg (Dimitriadis et al. 2016). A third report a day later, on November 3 at 04:24 UT, provided a more detailed classification (Pan et al. 2016). Although the report was cautious, it pointed out similarities between SN 2016hnk and the peculiar Ca-rich object PTF09dav (Sullivan et al. 2011; Kasliwal et al. 2012), from their peak magnitudes, similar colors, and red spectra with overall similar features. In particular, two deep absorptions at ∼5350 and ∼5540 Å and another absorption at ∼7120 Å (in the rest frame) were identified as Sc II and Ti II, respectively, by Sullivan et al. (2011). However, in Sect. 5 we show that these features are formed by different elements removing any resemblance between both objects, other than the extremely strong calcium features and low luminosity. Finally, another spectrum was obtained with Hiltner MDM-OSMOS on October 29.29, however observations were not reported.

Given its peculiar nature, i.e., it was unusually faint and showed spectral features that clearly distinguish it from other 1991bg-like SNe, SN 2016hnk was recognized as a rare and potentially interesting object. A follow-up campaign was started immediately the day after the PESSTO Telegram with the New Technology Telescope (NTT) and the Las Cumbres Observatory network, which ran for more than two months. All the spectra described above and those obtained during the PESSTO and Las Cumbres campaigns, together with extensive photometry from different sources, are presented in this work and described below.

2.1. Photometry

Most of the optical photometric coverage was obtained by the Las Cumbres Observatory 1.0 m telescope network (Brown et al. 2013) through PESSTO allocated time and a number of dedicated SN follow-up programs. In addition, supplementary data were obtained from the 2.0 m Liverpool Telescope (LT; Steele et al. 2004), the twin 0.5 m ATLAS telescope system (Tonry et al. 2018), the 0.5 m DEdicated MONitor of EXotransits and Transients (DEMONEXT; Villanueva et al. 2018), the NUTS survey using the ALhambra Faint Object Spectrograph and Camera (ALFOSC) at the NOT 2.56 m telescope, and prediscovery data from the All-Sky Automated Survey for SuperNovae (ASAS-SN; Shappee et al. 2014). All data have been reduced following standard procedures with dedicated pipelines for each instrument, which consist of bias subtraction, trimming of the over scan region, and flat fielding.

Near-infrared follow-up broadband imaging was obtained by both PESSTO and NUTS. PESSTO obtained NIR imaging with the 3.58 m NTT equipped with the Son of ISAAC (SOFI; Moorwood et al. 1998), while NUTS obtained data with NOTCAM on the NOT 2.56 m telescope. From these observations we were able to obtain six epochs of JHK-band photometry extending from +5 to +72 d from peak. The NTT images were reduced using the PESSTO pipeline described in Smartt et al. (2015), while the NUTS images were reduced using a slightly modified version of the NOTCam Quicklook v2.5 reduction package1. In both cases, the reduction procedures consisted of flat-field, distortion, cross-talk, and illumination corrections, trimming, sky subtraction, and stacking of the dithered images After proper sky subtraction, obtained through dithering, the individual science exposures were aligned and combined to produce a definitive stacked science image, to which we then applied accurate astrometry.

Given the location of the SN within the galaxy, we also obtained deep observations of the field in all bands and in all telescopes after more than a year from the discovery of the SN, to be used as a template for galaxy subtraction.

We performed photometry on the SN and up to 35 reference stars located within the field covered by these various instruments using the SuperNOva PhotometrY (SNOOPY2) package. First, each individual exposure is aligned to the United States Naval Observatory (USNO) B catalog (Monet et al. 2003), and a point spread function (PSF) with a width of the seeing is fitted to all reference stars shown in Fig. 1 to calculate the zero points. We then subtracted the host galaxy templates using hotPants (Becker 2015), and performed PSF fitting on the SN in the subtracted image (see Fig. 2 for an example). Finally, we calibrated all SN optical photometry against the Sloan Digital Sky Survey Data Release 15 (SDSS; ugriz; Aguado et al. 2019), the AAVSO Photometric All-Sky Survey (APASS; BV; Henden et al. 2015), and NIR photometry to the Two Micron All Sky Survey (2MASS; JHK; Skrutskie et al. 2006), accounting for corresponding airmass and color corrections in a standard manner.

thumbnail Fig. 2.

g-band imaging of SN 2016hnk and its host galaxy, before (left) and after (middle) subtracting the host galaxy template using hotpants integration in SNOOPY. This example shows how the subtracted image is clean enough to perform PSF fitting, as shown in the right panel where the SN PSF has been removed. A small residual of the galaxy core is shown on the middle and right panels (white blob), however it does not affect SN PSF fitting.

We present ATLAS photometry in c and o bands, and ASASS-SN V-band photometry in Tables C.1 and C.2, respectively. All other photometric observations and SN magnitudes are given in Tables C.3 (optical) and C.4 (NIR) in the AB system.

2.2. Spectroscopy

After the reported peculiarity of the object, we began an intensive optical spectroscopic follow-up campaign using several facilities. Most of the data (12 epochs, including 1 NIR spectrum) were obtained by the PESSTO collaboration with the NTT using the ESO Faint Object Spectrograph and Camera (EFOSC2; Buzzoni et al. 1984) instrument and SOFI (for the NIR spectrum). A description of the reduction pipeline can be found in Smartt et al. (2015).

These observations are complemented with data from the following sources: the LT using the Spectrograph for the Rapid Acquisition of Transients (SPRAT), and reduced and calibrated in wavelength using the pipeline described in Piascik et al. (2014), and the NOT using ALFOSC through the NUTS survey. These data were reduced in a standard manner with ALFOSCGUI; the Yunnan Faint Object Spectrograph and Camera (YFOSC) mounted on the Lijiang 2.4 m telescope (LJT; Zhang et al. 2014); the MDM 2.4 m Hiltner telescope at Kitt Peak National Observatory using the Ohio State Multi-Object Spectrogragh (OSMOS; Martini et al. 2011). The reductions follow normal procedures (Villanueva et al. 2016, 2018); the 6.5 m MMT telescope with the Blue Channel instrument (Schmidt et al. 1989). Standard procedures to bias-correct, flat-field, and flux calibrate the data using IRAF were followed; the 8.1 m Gemini South Telescope with the FLAMINGOS-2 spectrograph reduced following the prescriptions in Hsiao et al. (2019); and the two discovery spectra reported in the Transient Name Server and Astronomical Telegrams (ATELs) from the 4.1 m Southern Astrophysical Research Telescope (SOAR) and the 11 m Southern African Large Telescope (SALT) are also included.

Before performing any analysis, we scaled the spectra to match the photometry following a procedure known as color-matching (Hsiao et al. 2007). We first fit a third order spline to both the interpolated observed photometry with Gaussian processes at the epoch of each spectrum, and to the same synthetic bands from the observed spectrum. Next, the observed spectrum was multiplied by a wavelength-dependent factor that matches the two splines. This procedure ensures that colors from the corrected spectra are the same as those from the photometric measurements, which we assume to be more reliable.

Finally, we obtained six spectra of the nebular phase at four different epochs from +346 to +377 d with X-shooter (Vernet et al. 2011) under ESO program 0100.D-0738(A). X-shooter provides coverage of the full near-UV (UVB), optical (VIS), and NIR wavelength range (300−2500 nm) in a single exposure in three different channels. We used slits of 13 and 12 for the UVB and VIS arms, respectively. For the NIR arm, we opted for the 09JH slit to optimize the signal-to-noise (S/N) in the J and H bands while blocking the K-band region. With this setup, we achieved a relatively homogeneous slit transmission of ∼85% throughout the wavelength range. Exposure times were fixed to 3060 s in the UVB channel, 3000 s in the VIS, and 3000 s in the NIR. The six sky-subtracted (but not telluric-removed) observations were combined in a single spectrum to enhance the S/N. The full spectral log is presented in Table C.5.

3. Host galaxy and local environment of SN 2016hnk

Integral field spectroscopy of the type SB galaxy MCG -01-06-070 was obtained on 2018 October 11 as part of PISCO (Galbany et al. 2018).

Observations were performed with the Potsdam Multi Aperture Spectrograph (PMAS) in PPak mode (Roth et al. 2005) mounted to the 3.5 m telescope of the Centro Astronómico Hispano Alemán (CAHA) located at Calar Alto Observatory. PPak is a bundle of 331 fibers of 27 diameter ordered in a hexagonal shape with a filing factor of the field of view (FoV) of 55% (Verheijen et al. 2004; Kelz et al. 2006). Three 900 s exposures were performed with the V500 grating of 500 lines mm−1, which provides a spectral resolution of ∼6 Å and covers the whole optical range from 3750 to 7300 Å. The second and third exposures were shifted by a few arcsec with respect to the first pointing to provide spectroscopic coverage of the full hexagonal ∼1.3 arcmin2 FoV with 1″ × 1″ pixels, which corresponds to ∼4000 spectra. Reduction was done with the pipeline used for the Calar Alto Legacy Integral Field Area survey (CALIFA) Data Release 3 (v2.2; all details can be found in Sánchez et al. 2016).

The analysis was performed in a similar way as that presented in Galbany et al. (2014, 2016a). We used a modified version of STARLIGHT (Cid Fernandes et al. 2005; López Fernández et al. 2016, priv. comm.) to estimate the fractional contribution of different simple stellar populations (SSP) from the Granada-Miles base (González Delgado et al. 2015) with different ages and metallicities to the stellar continuum in the spectra; we added dust effects as a foreground screen with a Fitzpatrick (1999) reddening law and RV = 3.1. We then obtained pure gas emission spectra by subtracting the best SSP fit from the observed spectra and estimated the flux of the most prominent emission lines after correcting for dust attenuation from the Balmer decrement (assuming case B recombination; Osterbrock & Ferland 2006; a Fitzpatrick 1999 extinction law; and RV = 3.1). Line fluxes were used to estimate a number of galactic parameters, such as the ongoing star formation rate (SFR; from the Hα emission line), the weight of young-to-old populations (from the Hα equivalent width; HαEW), and oxygen abundances.

In this work we are mainly interested in studying the local properties of the SN 2016hnk environment and the integrated properties of its host. For the latter, we repeated the above procedures on an integrated spectrum of the host galaxy, summing up the spectra in all spaxels, and accounting for spatial covariance in the error budget. Wide-field IFS allows SN explosion site parameters to be compared not only to those of the overall host, but to all other stellar populations found within hosts, therefore exploiting the full capabilities of the data. Following Galbany et al. (2016b) we developed a method to characterize the SN environment comparing its properties to all other nebular clusters in the galaxy. Using our extinction-corrected Hα map we selected star-forming HII regions across the galaxy with a modified version of HIIEXPLORER3 (Sánchez et al. 2012), which is a package that detects clumps of higher intensity in a map by aggregating adjacent pixels. This procedure selected 56 HII clumps with an average radius of ∼700 pc. Once the HII regions were identified, the same analysis described above was performed on the extracted spectra.

Figure 3 shows two-dimensional maps of the continuum-subtracted Hα emission with the HII region segmentation performed with HIIEXPLORER on top (left), the Hα equivalent width (middle), and the rotational velocity field (right). In all panels the red star represents the SN position. SN 2016hnk occurred between two HII clumps, within the small disc created by the bar present in the host galaxy. One of these clumps has both high Hα brightness and HαEW, pointing to a strong component of young stellar populations and ongoing star formation. However the SN is closer to the clump showing bright Hα emission but lower HαEW, which means ongoing star formation in a location dominated by older nonionizing stellar populations, as it is seen in the star formation history of Fig. 5. At the SN location the rotational velocity of the galaxy is 82 ± 6 km s−1, which we subtracted in all SN velocities measured in Sect. 5.

thumbnail Fig. 3.

Continuum-subtracted Hα emission (left), HαEW (middle), and Hα velocity (right) maps of MCG -01-06-070. HII regions extracted with HIIEXPLORER are indicated in brown polygons in the left panel. The red star and the black/white × represent the locations of SN 2016hnk and MCG -01-06-070 nucleus, respectively.

The resulting integrated spectrum of MCG -01-06-070 is presented in Fig. 4 with the best STARLIGHT SSP fit. The recovered star formation history is also shown in the inner plot. The stellar mass of the galaxy is 4.84 × 1010 M, and the average stellar age of the galaxy is 690 Myr, although we were able to recover three main peaks of formation: a very old component at around 6 Gyr, an intermediate burst at around 1 Gyr, and a young component of a few Myr. The best fit also provides an oversolar stellar metallicity of Z = 0.03 and a nonsignificant optical stellar extinction of AV = 0.02 ± 0.06 mag.

thumbnail Fig. 4.

Top: MCG -01-06-070 integrated spectrum (blue) and best SSP fit from STARLIGHT (red). Vertical strips (in beige) represent masked regions in the fit due to known emission lines (such as Balmer lines, oxygen, and WR bumps) and regions with skylines (in brown). Inner panel: star formation history of the spectra (in purple), where brown ticks on top represent the ages of the models used in the fit. The S/N of the spectrum, optical extinction, average stellar age, and average metallicity are reported in the top left corners. Bottom rows: distributions of the Hα equivalent width, SFR, oxygen abundance, and average stellar age of all 180 SN Ia host galaxies in PISCO (Galbany et al. 2018). The large colored dots represent the position of MCG -01-06-070, and the three small dots represent other 1991bg-like SNIa hosts in PISCO.

The bottom panels of Fig. 4 aim to compare the properties of SN 2016hnk and the other three 1991bg-like SNIa host galaxies with those of all other 180 SN Ia host galaxies in PISCO (Galbany et al. 2018; published sample plus observations up to March 2019 to be presented elsewhere). We present distributions of the HαEW, SFR, oxygen abundance (12 + log10 O/H) in the Dopita et al. (2016) scale, and the average light-weighted stellar age. For MCG -01-06-070 we found an HαEW of 7.04 ± 0.18 Å, a SFR of 0.649 ± 0.234 M yr−1, and an oxygen abundance of 8.64 ± 0.11 dex. The big colored dot represents the position of MCG -01-06-070 parameters in such distributions, and the vertical and horizontal strips show the actual parameter value of this galaxy and the rank within PISCO. The other three 1991bg-like SN Ia hosts are also represented with smaller colored dots. Overall, MCG -01-06-070 has no extreme properties and lies within the 35−50% rank in all measured parameters within the PISCO sample. However, it has a younger average stellar population age and, equivalently, a higher HαEW than the other 1991bg-like SN Ia hosts, which show less evidence of ongoing or recent star formation.

Regarding the local environment of the SN, in Fig. 5 we present the spectrum of the closest H II region from the segmentation described above. We also show two panels at the top right centered at Balmer Hα and Hβ wavelengths, where the best stellar continuum fit is subtracted from the observations leaving the pure emission spectra. In those, the [O III], [N II], and [S II] emission lines (in pale pink) and Hα and Hβ (in orange) are highlighted. The extinction-corrected fluxes of these seven lines were used to calculate the HαEW, the SFR, and the oxygen abundance, which are shown in the bottom panels together with the distributions of the same parameters measured from the spectra of all 56 HII regions found in our segmentation of the host galaxy. We included the distribution of the parameters as measured from the spectra, and once weighted by the contribution of the HII region to the total SFR (Lyman et al. 2018). For the HII region closest to SN 2016hnk we get a HαEW of 4.39 ± 0.04 Å, a SFR of −3.44 ± 0.34 M yr−1, an oxygen abundance of 8.65 ± 0.11 dex, and an average stellar age of 1.2 Gyr; there is a dominant peak (∼65%) at the average age and only a ∼4% contribution from young (< 15 Myr) populations. Interestingly, in this case we find that these parameter values deviate from the average properties of all HII regions. The local HαEW is at the lower end of the distribution indicating a residual contribution of young populations compared to other regions of the galaxy and can also be seen in the stellar age distribution, which agrees with the expected old progenitors for 1991bg-like SNe Ia. Conversely, this SN is in an average star-forming region (25% to 50% rank once weighted for the SFR) and has an average oxygen abundance (∼50th percentiles).

thumbnail Fig. 5.

Top: similar to Fig. 4, the spectrum of the SN 2016hnk parent HII region (blue) and best SSP fit from STARLIGHT (red) are shown. Right panels: stellar subtracted spectra around the Hα and Hβ wavelengths. Bottom rows: distributions of the gas-based parameters of all 56 HII regions in MCG -01-06-070, weighted by the contribution of the region to the total SFR of the galaxy (dark lines) or unweighted (light lines). In this case, colored horizontal bands represent the band where the SN 2016hnk parent HII region lies within.

The Hα to Hβ flux ratio also provides an estimate of the color excess E(B − V) under certain assumptions (case B recombination; Osterbrock & Ferland 2006). We get a ratio of 4.36 ± 0.35, which compared to the theoretical value of 2.86 for a non-extincted ratio and using a Fitzpatrick (1999) extinction law, gives an E(B − V) of 0.363 ± 0.049 mag.

4. Observational photometric properties

Final ocBVugrizJHK light curves are presented in Fig. 6. SN 2016hnk photometry has rather sparse coverage before peak, since it was only observed in ATLAS oc bands and ASAS-SN V band during the rise. However, from the shape of the light curves in all bands shown in Fig. 6, we can infer that the follow-up campaign started very close to the epoch of peak brightness (tmax).

thumbnail Fig. 6.

Optical and NIR light curves. A gray solid vertical line represents the epoch of peak brightness in the B band obtained from a SNOOPY fit to the cyan and orange ATLAS light curves, which are shown in orange and cyan solid lines. The ticks at the bottom represent epochs with optical spectra available. The dotted vertical lines delineate the last non-detection (yellow) and first detection (orange) epochs. The open symbols with downward arrows correspond to ASAS-SN (circles), ATLAS (crosses/diamonds), and u-band (triangles) non-detections.

To determine tmax precisely, we used the SuperNovae in object-oriented Python (SNOOPY; Burns et al. 2011, 2014) fitter, which allows us to consider different extinction laws and amounts of reddening in our light curve fits. Also, SNOOPY makes studies involving color curves very flexible; such include the study of color laws, the construction of bolometric light curves, and estimates of the dimensionless sBV parameter, which gives the stretch factor of the B − V color curve with respect to the typical behavior of a SN Ia having a of 30 days. The sBV parameter provides a more physically significant interpretation of the light curve width than similar parameters, as it measures the time at which iron recombines from Fe III to Fe II (e.g., Burns et al. 2014; Wygoda et al. 2019). Therefore, sBV may be providing more information when it comes to fast decliners and peculiar types, such as SN 2016hnk. We fit the light curves of the three bands that have pre-maximum observations, as well as the B band, with the 1991bg-like built-in model, and obtain a tmax of MJD = 56789.5 ± 3.3 (ATLAS-co fits are included in Fig. 6). This estimate constrains the SN rise time to 17.4 ± 5.5 days, hence to the range from 11.9 to 22.9 days. This value is consistent with most SNe Ia and 1991bg-like SNe, but larger than in Ca-rich transients, which have faster rise times on the order 9−15 days (Taubenberger 2017).

A SNOOPY light curve fit using the Markov Chain Monte Carlo fitter and the 91bg model was simultaneously performed in all bands and to the whole temporal extent, fixing the tmax and also inputting the E(B − V) obtained from the analysis of the most nearby environment as a Gaussian prior. Our initial attempt provided a reasonable fit up to 15−20 days past maximum (See Fig. 7), but the most remarkable mismatch was found at late times (> +20 d) where SN 2016hnk showed an excess in brightness in all bands. Best-fit parameters were RV = 2.1 ± 0.4, E(B − V) = 0.547 ± 0.064 mag, and sBV = 0.438 ± 0.030.

thumbnail Fig. 7.

Multiband light curves of SN 2016hnk (filled circles) compared with a 1991bg-like template (lines). Best-fit SNOOPY parameters are RV = 2.1 and a host-galaxy color excess of E(B − V)host = 0.55 mag. We note that beyond +20 d the template fit does not succeed in capturing the SN brightness.

Figure 8 shows SN 2016hnk B band compared to a variety of SNe Ia: the normal SN Ia 2011fe (Richmond & Smith 2012); three 2002es-like SNe Ia, 2002es (Ganeshalingam et al. 2012), 2010lp (Kromer et al. 2013), iPTF14atg (Cao et al. 2015); a 2002cx-like SN Ia, 2008ae (Krisciunas et al. 2017); the transitional SN 1986G (Phillips et al. 1987); the highly extincted SN Ia 2006X (Wang et al. 2008; Krisciunas et al. 2017); and a sample of 1991bg-like SNe Ia including 1991bg (Filippenko et al. 1992; Leibundgut et al. 1993; Turatto et al. 1996; Krisciunas et al. 2004), 1999by (Ganeshalingam et al. 2010; Garnavich et al. 2004), 2005bl (Taubenberger et al. 2008; Friedman et al. 2015; Hicken et al. 2009; Contreras et al. 2010), 2005ke (Friedman et al. 2015; Hicken et al. 2009; Contreras et al. 2010; Krisciunas et al. 2017), and 2006mr (Contreras et al. 2010; Stritzinger et al. 2010), the fastest SN Ia ever found. These 1991bg-like SNe are used in the following sections as a base for comparison. All light curves have been rest-framed and shifted vertically to match the SN 2016hnk at peak brightness. For reference, two vertical lines correspond to epochs at peak and at +15 days, and three horizontal lines cross the +15 d line at the position where light curves with Δm15(B) = 1.3 (similar to SN 2016hnk), 1.5 (intermediate), and 1.7 mag (approximately the minimum for 1991bg-like SN Ia) should cross. SN 2016hnk has a Δm15(B) most similar to 2002es-like objects, but with two peculiarities: this SN shows a much narrower light curve before +15 d, and has a much shallower decay at later epoch. In fact, up to 12−13 days past maximum, SN 2016hnk follows very closely the low-luminosity transitional SN Ia 1986G light curve, being only slightly wider than all other 1991bg-like SN Ia light curves. Focusing on the late (> +15 days) light curve, we clearly see that SN 2016hnk has the brightest tail of all objects in comparison. This figure clearly shows that Δm15(B) is not capable of capturing all the diversity in SNe Ia and may mask intrinsic differences in light curve shapes, as it is seen in SNe 2002es and 2016hnk. We expand on the explanation of this light excess in Sect. 4.1.

thumbnail Fig. 8.

B-band light curves of SN 2016hnk and a variety of SNe Ia of different subtypes: normal SN Ia (2011fe), 2002es-like (2002es, 2010lp, iPTF14atg), 2002cx-like (2008ae), transitional (1986G), and 1991bg-like (1991bg, 1999by, 2005bl, 2005ke, 2006mr). All data were rest-framed and scaled to match the SN 2016hnk peak. Vertical lines represent epochs of the peak and +15 days, and horizontal lines indicate the crossing point at 15 days of a SN with Δm15(B) = 1.3, 1.5, and 1.7 mag, respectively. SN 2016hnk has one of the lowest Δm15(B) but is the object with the brightest tail in both bands.

As an independent check to increase the robustness of our findings, we used the Color-Decline-Magnitude intercept Calibration (CDMagic; Höflich et al. 2017) to estimate the host-galaxy reddening and preferred dust extinction law. First, following Goldhaber et al. (2001), we measured the stretch factor s that multiplies the temporal axis of the light curve providing the best match of the Hsiao et al. (2007) SN Ia template to the observations. In our case, we restricted the fit to the first 15 days, to avoid the dependence of the light curve on secondary parameters and the light excess detected for SN 2016hnk at later times, and obtained an s = 0.59 ± 0.04. The modified decline rate Δm15, s(B) is then defined as the magnitude difference between the peak and the reference epoch, 15 days × s (measured in the template), divided by the s factor [Δm15, s = , which in our case is = = 1.803 ± 0.200 mag. The measured value of Δm15, s(B) determines the theoretical model from Höflich et al. (2017) that is used as a reference to correct the B- and V-band light curves, the B versus (B − V) and V versus (B − V) diagrams, and the (B − V) evolution of SN 2016hnk. In the CDMagic fit we fixed the to the value previously obtained with SNOOPY and the distance modulus measured from the assumed cosmology (μ = 34.195 mag), and only left as free parameters E(B − V) and RV. While the color excess is estimated from the shift between the uncorrected (as observed and Milky Way dust-corrected) and the corrected points that best match the template, the RV is measured from the direction of this shift in the color-magnitude diagrams. In Fig. 9 we show, as an example, the absolute V and vs. B − V diagram, including a number of SNe Ia from the Carnegie Supernova Project I (CSP-I; Krisciunas et al. 2017) and the 1991bg-like SN Ia 2005ke (empty triangles) for reference, models from Höflich et al. (2017) with different stretch factor (green lines), and SN 2016hnk data before (blue) and after corrections (red). We obtained an E(B − V) of 0.45 ± 0.08 mag, consistent with previous methods, and an RV = 2.1 ± 0.6, which we use hereafter; however, we note that our results are not significantly dependent on the particular choice of RV.

thumbnail Fig. 9.

Representative panel for the CDMagic method, where a simultaneous fit in the B- and V-band light curves, the B vs. (B − V) and V vs. (B − V) diagrams, and the (B − V) evolution, is performed. More details are given in, for example, Fig. 13 in Höflich et al. (2017). We show the absolute V-band vs. (B − V) diagram populated with eight normal SNe Ia from the CSP-I (Krisciunas et al. 2017; SNe 2004ef, 2005eo, 2005el, 2005iq, 2005ki, 2005M, 2006bh, 2006D; in differnet symbols), the SN Ia 1991bg-like 2005ke (empty triangles), and the uncorrected (blue) and corrected (red) SN 2016hnk data points. Models from Höflich et al. (2017) are shown in green for reference; the two arrows are the reddening vectors representing the direction to which the data points are shifted when affected by RV = 3.2 and RV = 1.0. SN 2016hnk corrected data points lie on top of the other subluminous SN Ia 2005ke, and from the shift we estimate an RV = 2.1 ± 0.6.

The observed apparent B-band peak magnitude is 19.003 ± 0.088 (at epoch +0.7 d). Using tmax as a reference, we estimated Δm15(B) by interpolating the flux at 15 days with a third order spline and obtained a value of 1.279 ± 0.071 mag. We then applied the correction provided by Phillips et al. (1999) to get the true Δm15(B)true = Δm15(B)obs + 0.1 × E(B − V), and found 1.324 ± 0.096 mag. In order to find the absolute B-band peak magnitude, we corrected the observed value for Milky Way reddening of E(B − V) = 0.0224 ± 0.0008 mag from extinction maps of Schlafly & Finkbeiner (2011) and the average Milky Way value for RV of 3.1, and the SN host galaxy dust extinction with the color excess of E(B − V) = 0.45 mag and an RV of 2.1 found above. We obtained an = 0.069 mag and = 0.945 mag, which correspond to = 0.092 mag and = 1.395 mag. Considering both corrections we obtained a corrected absolute B band magnitude of −16.656 ± 0.162. This is at the low end of the luminosity distribution for SNe Ia (Ashall et al. 2016a), and confirms SN 2016hnk as a subluminous SN Ia.

This is clearly seen in Fig. 10, where we show the absolute magnitude versus light curve width relation for SN 2016hnk and SNe Ia from the CSP-I (Krisciunas et al. 2017). In the top panel we used the Δm15(B) parameter. While this parameter correlates strongly with the SN peak brightness, it is known to become degenerate for the fainter SNe Ia showing a branch that follows the same relation as brighter SNe Ia, and a second branch that breaks the relation. In this plot, SN 2016hnk is an outlier, and the reason for this is discussed above and can be seen in Fig. 8. It has a very low Δm15(B) when compared to other subluminous SNe Ia, which usually have larger values (> 1.7 mag). In the bottom panel we used the sBV parameter instead. Burns et al. (2014) showed that this parameter keeps more information and breaks the degeneracy in the peak brightness versus light curve width relation for the fainter objects. In this case, SN 2016hnk perfectly follows the relation.

thumbnail Fig. 10.

Luminosity width relation as a function of Δm15(B) (top) and sBV (bottom) for a selection of CSP-I SNe Ia from Krisciunas et al. (2017) with the addition of SN 2016hnk (black circle). We note that SN 2016hnk sits off the Δm15(B)−MB relationship because of the degeneracy in Δm15(B) for fast declining SNe Ia.

In Fig. 11 we present several SN 2016hnk color curves (top panel), and the B − V color curve compared to the sample of 1991bg-like SNe Ia presented above, SNe 2011fe, and 2002es (bottom panel). The slope of the B − V color during the Lira law regime (between 35 and 80 days after maximum) provides an independent test on the nature of the reddening process in play (Lira et al. 1998; Phillips et al. 1999). Extinction by interstellar material (ISM) produces a shift in the B − V color curve to higher (red) values, without significant changes on the decline rate, whereas sufficiently close dust (≲10 pc, either circumstellar material −CSM− or ISM) can also affect the decline rate due to an evolving E(B − V) (Förster et al. 2013; Bulla et al. 2018a).

thumbnail Fig. 11.

SN 2016hnk color curves (top panel), the B − V color curve of SN 2016hnk, the 1991bg-like comparison sample, and SNe 2002es and 2011fe for reference. The shaded region projects the SN 2016hnk color curve down by 0.5 mag, on top of the bulk of other 1991bg-like SNe Ia. In pink we show the linear fit the to B − V curve from 30 to 70 days past B-band maximum.

It can clearly be seen that the observed B − V curve is well above (redder) those other subluminous SNe at all epochs. While other 1991bg-like objects have B − V colors at maximum light around 0.4−0.8 mag, SN 2016hnk has 1.5 mag. Up to +20 days, the color curve is consistently half a magnitude larger when compared to other 1991bg-like SNe Ia. This is in agreement with the E(B − V) values found previously in the CDMagic and SNOOPY analyses, and confirms a large amount of reddening. On the other hand, the post maximum slope is clearly shallower, which is opposite of what is expected for objects affected by large amounts of extinction (Förster et al. 2013). We fit a straight line between epochs (+30, +70) and found a slope of −0.005 ± 0.001 mag day−1, which is in the low regime than is expected for objects with insignificant reddening and is in conflict with the expectations for CSM scattering. However, this may also be caused by either difference in the central density of the progenitor at the explosion (Gall et al. 2018a) or by differences in the transmission filters, which may cause differences in the color curve slope (Höflich et al. 2017).

4.1. Late light curve excess

Figure 8 shows that SN 2016hnk presents the shallowest brightness decline at later epochs (≳12 days) compared with a variety of thermonuclear events. To explain this light excess we performed a number of tests confirming, first that SN 2016hnk is similar to SN 1991bg when similar host-galaxy reddening with low extinction law is artificially applied to SN 1991bg, although it is still redder and no Na I D lines are detected at any epoch in its spectra; and second when a simple light-echo (LE) model is removed from the observed SN 2016hnk light curve, a SN 1991bg-like model fits reasonably well to the observations at late epochs, however no LE (spectrum at maximum) features are detected in the late-time spectra. Details of those analysis are left for Appendix A, in this section we concentrate on our most plausible explanation.

At late epochs, the slow decline in the SN 2016hnk light curves suggests an additional energy source beyond 56Co to 56Fe radioactive decay. The most likely mechanism for this particular light curve shape is multiple scattering with either a circumstellar or interstellar dust cloud intervening in the line of sight, which would in turn causes the total-to-selective extinction ratio to be steeper (lower RV; Goobar 2008). Amanullah & Goobar (2011) presented predictions of how scattering by dust affects the SN light curve shape depending on the distance to the interfering cloud and its thickness (opacity). Qualitatively, the light curve suffers some broadening before 15 days and produces a light excess after 15 days (as observed in, e.g., SN 2006X), which is more pronounced with thicker dust clouds (prolonged scattering of SN photons). We see similarities between simulated light curves shown in Fig. 3 of Amanullah & Goobar (2011) for the 1018 cm (∼0.1 pc) case, with a broadening and a plateau at latter times, and the SN 2016hnk B-band light curve when compared to all other 1991bg-like objects. Importantly, the net effect of this multiple scattering would be a reduction of the Δm15(B) parameter (of up to ∼0.4 mag) for similar E(B − V) as in SN 2016hnk, as we have also shown in Fig. 8.

We used the Monte Carlo simulations by Bulla et al. (2018b) to estimate the distance to the intervening dust cloud and, at the same time, try to fit multiband light curves of SN 2016hnk with this model. In Table 2 we provide the dust properties of the model for each band, which are based on a LMC-type dust composition from Weingartner & Draine (2001). We use the 1991bg-like SN 2005ke (Krisciunas et al. 2017) light curves as a reference in our simulation. In Fig. 12, we show the BVri reference light curves and the resulting simulated light curves affected by nearby dust with E(B − V) = 0.5 mag at a distance of 1.0 ± 0.5 pc from the SN (blue filled strip), and at larger distance of 100 pc (blue dashed line). While dust at larger distances would affect the light curve by dimming the light and keeping the observed light curve shape, nearby dust would broaden it at early times and provide an extra energy source at later times. SN 2016hnk observations are included in all panels to show the good agreement between data and model. Based on this simulation, SN 2016hnk would have the same Δm15(B) than SN 2005ke, if it was not affected by the dust cloud.

thumbnail Fig. 12.

BVri observed light curves (in circles), compared with the multi-scattering dust model of Bulla et al. (2018b) using 1991bg-like SN 2005ke as a reference (black dotted line). Simulated light curves affected by nearby dust with E(B − V) = 0.5 mag at a distance of 1.0 ± 0.5 pc from the SN (blue filled strip) provide very reasonable fits to the observations, while a model in which dust is located at larger distance (100 pc; blue dashed line) just dims the light curves by keeping the same shape.

Table 2.

Scattering parameters (albedo, average of the cosine of the scattering angle, and absorption cross section divided by dust mass) for optical photons from Weingartner & Draine (2001) corresponding to interstellar extinction in the LMC.

5. Observational spectral properties

In this section, we describe our spectral data in terms of their observed properties, and then we present the results of a tomographic analysis of SN 2016hnk in Sect. 6.1 and a comparison to explosion models in Sect. 6.2.

The seq21. 13, covering epochs from +0.8 to +68.5 days. In Fig. 14 we show spectra of SN 2016hnk and the subluminous SN comparison sample at three different epochs: around maximum light, at one week past maximum, and around three weeks past maximum. We also include spectra of PTF09dav (Sullivan et al. 2011) because SN 2016hnk was claimed to be similar to this Ca-rich SN, and of SN 2002es (Blondin et al. 2012) to enable comparisons to this different SN Ia subclass. In all three panels we can see the resemblance of SN 2016hnk to other 1991bg-like SNe Ia, especially with respect to the Ti II feature at ∼4400 Å, O Iλ7774, and other intermediate mass elements (IME) characteristic of SNe Ia. However, there are a number of features that make SN 2016hnk special:

thumbnail Fig. 13.

Optical spectral sequence for SN 2016hnk including 22 observations listed in Table C.5 (shifted by 1.5e−16 erg s−1 cm−2 Å−1). In gray we present the color-corrected, Milky Way reddening-corrected, and rest-framed spectra, and in black a smoothed version using a 10 pixel box. Rest-frame epochs for each spectrum are attached in blue.

thumbnail Fig. 14.

Spectra of SN 2016hnk at three different epochs (around maximum, after a week, and around three weeks past maximum light) compared to other subluminous 1991bg-like SNe Ia, PTF09dav (Ca-rich), and SN 2002es (2002es-like SN Ia). The main typical features of 1991bg-like SNe Ia are identified in the top panel. Vertical lines in the top panel are located at the minima of a few features and serve as a comparison of the velocities with respect other objects. All spectra have been Milky Way extinction-corrected and shifted to the rest frame.

– First, this object has by far the broadest NIR Ca II triplet absorption. This absorption remains wide and deep compared to other 1991bg-like objects until at least three weeks past maximum (see bottom panel of Fig. 14).

– A second peculiarity are the three features between 6400 and 7300 Å, especially in the spectrum around maximum light, which fall in a wavelength region that is very flat in other objects. We identified these features to be caused by Fe II, Co I, and Co II lines (see discussion in Sect. 6). Some shallower features are also present in SN 2005bl, and Hachinger et al. (2009) identified these as being mostly caused by Ti II, however their Ti abundance is much larger than that found through explosion models and they do not reproduce all of the features.

– As Folatelli et al. (2013) previously noted, some extremely cool SNe show a double structure in the profile of the Si IIλ5972 absorption with a weak component on the blue side of the main absorption. This weaker absorption was attributed to the Na I D doublet that grows stronger at lower temperature, or to Fe II lines at low velocity. SN 2016hnk is the only object in the sample that has this bluer absorption detached from the main Si II feature. Related to this, in the +3 week spectra all other 1991bg-like SNe show a deep Na I D absorption, however SN 2016hnk does not show such an absorption (see also Appendix A).

– All SN 2016hnk spectra are clearly redder than those of the comparison sample. All spectra in Fig. 14 have been corrected for Milky Way reddening and shifted to the rest frame, so that reddening excess has to be caused by host-galaxy reddening or intrinsic color. As shown in the previous section, we estimated a color excess of 0.4−0.5 mag, which is corroborated in this work when comparing SN 2016hnk to the spectra of other 1991bg-like SNe Ia.

– Another interesting thing about SN 2016hnk is the presence of two narrow lines between 5200 and 5600 Å, which we identify as the typical S II “W” feature seen in all SNe Ia, and that were used to claim similarities to the peculiar Ca-rich SN PTF09dav (Pan et al. 2016). Sullivan et al. (2011) were able to reproduce the observed shape of these features in PTF09dav spectra by fitting synthetic spectra simulated with SYNAPPS (Thomas et al. 2011) and inputting Sc II instead of S II. We show in Sect. 6 that we are able to reproduce reliably these and other features with typical elements found in other 1991bg-like SNe but assuming lower ionization/excitation.

In the top panel of Fig. 14 we highlight a few features in SN 2016hnk that we characterized in the spectral sequence in terms of the evolution of their pseudo-equivalent width (pEW) and velocities. We show our measurements in Fig. 15 and list them in Tables C.6 and C.7. The most remarkable characteristics of SN 2016hnk compared to all other subluminous Ia are as follows: (i) This object shows stronger Ca II H&K and NIR features and larger expansion velocities. However, we note that there are indications of a shallow feature on the blue side of the Ca NIR in the spectra of SNe 1991bg and 2005bl, which may simply be deeper in SN 2016hnk. (ii) While SN 20016hnk has larger S IIpEW than other SNe, the expansion velocities are clearly lower. This can be seen in the top panel of Fig. 14 where we plotted vertical lines corresponding to the velocities of SN 2016hnk to better illustrate the difference from other objects. (iii) While the Si II velocity and pEW measured from the 6355Å line are consistent with other subluminous SNe, the Si IIλ5972 velocity is significantly lower and the pEW decreases very rapidly. This may be caused by the detachment of the bluer feature. (iv) We see indications of two different behaviors regarding the O Iλ7774 line within the group of subluminous SNe Ia: SNe 2016hnk, 1991bg, and 2006mr show larger pEW and velocities than SNe 1999by, 2005bl, and 2005ke. However, as we show in Sect. 6, this line is caused by a combination of O I, Mg II, and Si II, and this excess may be caused by the large concentration of Mg II in SN 2016hnk.

thumbnail Fig. 15.

Evolution of pseudo-equivalent widths pEW and velocities as a function of time for the most prominent features highlighted in the top panel of Fig. 14. Shaded regions represent the 1σ dispersion of COOL (brown) and BROAD LINE (orange) evolution in the CSP-I sample from Folatelli et al. (2013).

The pEW and velocities of the two silicon features have been widely used to identify subgroups within SNe Ia faling in different regions on a number of diagrams. For instance, Branch et al. (2006) proposed mapping SN properties using the pEWs of the two Si II features at 5972 and 6355 Å at the epoch of maximum brightness, and Wang et al. (2009) proposed a diagram of the velocity and the pEW of the same Si IIλ6355 feature. While the latter uses the same line, the former uses two different lines, taking into account different levels of excitation. Less luminous SNe Ia have stronger Si II 5972 features at maximum light. This is caused by the recombination from Si III to Si II happening earlier in less luminous objects, and the 5972 Å feature having a higher excitation level compared to the 6355 Å feature, which is saturated at maximum for all SNe Ia (Hachinger et al. 2008). Furthermore, if Si lines are at low velocities the 6355 Å feature can also be contaminated by Fe IIλ6516 (see Sect. 6.1). Another difference between the 5972 and 6355 Å Si II features is that while the range of pEW for Si 5972 Å is due to differences of depth, the range of pEW in the 6355 Å feature is because of width (e.g., a larger or small Si region in velocity space), as can be seen in the relation between velocity and pEW of the Wang et al. (2009) diagram.

We present both diagrams in Fig. 16 with the background filled in with measurements from the CfA SN Ia sample in Blondin et al. (2012) and the CSP-I sample in Folatelli et al. (2013), and with measurements for SN 2016hnk and the comparison sample shown on top. As expected, in the Branch et al. (2006) diagram (left panel), all 1991bg-like SNe Ia fall in the COOL region, but SN 2016hnk falls in the BROAD LINE region because, although it has consistent Si IIλ6355 pEW, has a shallower Si IIλ5972 feature compared to all other 1991bg-like SNe Ia.

thumbnail Fig. 16.

Branch et al. (2006) and Wang et al. (2009) diagrams. The left panel is populated with measurements from the CfA sample by Blondin et al. (2012) (open symbols) and from the CSP-I sample by Folatelli et al. (2013) (filled symbols). Right panel: only the measurements from Folatelli et al. (2013) are included. In both panels we added our own measurements of SN 2016hnk and the sample of subluminous 1991bg-like SNe Ia. SN 2002es was added only in the left panel because with a velocity of ∼5800 km s−1 it falls outside the range shown in the right panel.

Groups are a bit mixed in the Wang et al. (2009) diagram (Fig. 16, right panel), with overlapping regions between COOL and NORMAL and COOL and BROAD LINE SNe Ia. Accordingly, some 1991bg-like SNe Ia fall on top of the other two groups of CSP-I SNe. In contrast to the left panel of Fig. 16 (Branch et al. diagram) the location of SN 2016hnk in this diagram is more consistent with the COOL objects, having the largest Si IIpEW of all objects with similar expansion velocities. This diagram is efficient in separating high-velocity SNe, which all fall on the right side of the plot, and overluminous 1991T-like shallow silicon SNe Ia, which fall on the bottom side of the plot. However this figure is not able to clearly separate the core normal and the cool groups, as the Branch et al. diagram does, mostly because it does not capture excitation differences.

In Fig. 17 we present two additional historical diagrams from Benetti et al. (2005) (in the left column), which combine spectroscopic parameters with the light curve width parameter Δm15(B). In the top panels we show the R(Si II) parameter, the ratio of depth of the Si II 5972 over 6355 Å features, and in the bottom panels we show the Si II 6355 Å velocity gradient as measured in the range from peak up to the epoch when the feature disappears. We also populated these diagrams with literature measurements from Benetti et al. (open symbols) and the CSP-I (filled symbols). Corresponding parameters for our sample of 1991bg-like SNe Ia are listed in Table 3. As before, all 1991bg-like SNe Ia fall in the region of the diagram where all other subluminous objects from the literature are located, however SN 2016hnk falls in an unpopulated region closer to NORMAL SNe Ia. Although its R(Si II) would be lower but consistent with COOL SNe, the parameter that is shifting SN 2016hnk far from the COOL group is its low Δm15.

thumbnail Fig. 17.

Benetti et al. (2005) diagrams of R(Si II) and Si II velocity gradient as a function of Δm15 (left column). We replicated these diagrams using Δm15, s (center) and sBV (right column) instead as a SN light curve width parameter. In all panels, open symbols represent measurements from the original SN sample in Benetti et al. (2005), while filled symbols are measurements from the CSP-I (Folatelli et al. 2013).

Table 3.

Parameters measured in this work and used in the diagrams presented in Fig. 17.

As a comparison we repeat those same two diagrams using the Δm15, s and sBV parameters, under the hypothesis that they provide a more physically significant interpretation of the light curve width because they measure the time at which iron recombines from Fe III to Fe II and breaks the degeneracy in the absolute magnitude versus Δm15 (Phillips) relation. We filled the diagrams with measurements of Δm15, s for the CSP-I from Höflich et al. (2017) and the original Benetti et al. (2005) sample measured in this work. Similarly, we used sBV for the CSP-I sample from Burns et al. (2014) and the Benetti et al. (2005) sample measured in this work. In both cases, while 1991bg-like SNe Ia are kept in the region defined by the COOL objects, SN 2016hnk is also shifted toward the COOL group. This shows how the Δm15 parameter is ambiguous when the SN evolves too fast (but Δm15, s and sBV do not have this issue).

5.1. Near-infrared spectra

The two NIR spectra available are presented in Fig. 18 together with comparison spectra of a normal SN Ia (SN 2011fe; Mazzali et al. 2014), a transitional object (iPTF13ebh; Hsiao et al. 2015), and a SN 1991bg-like (SN 1999by; Höflich et al. 2002). SN 2016hnk spectra show a few features between 9500 and 13 000 Å, up to the point where a telluric band starts to affect the data. According to Marion et al. (2003, 2009), the two deeper fetaures at ∼10 000 and 11 500 Å could be caused by Fe II but they are deeper in SN 2016hnk than in other SNe Ia. Only SN 1999by show similar shallower features at those wavelengths. However, at a week past maximum SN 1999by did not transition from the Si/S/Ca to the Fe/Ni-dominated regime (see Fig. 3 in Höflich et al. 2002). SN 2016hnk shows lower velocities and the transition is faster (around 1day), so that the Fe features in the NIR appear earlier. A more detailed line identification is left for Sect. 6 in which the NIR spectra are compared in more detail to SN 1999by models from Höflich et al. (2002).

thumbnail Fig. 18.

Near-infrared spectra of SN 2016hnk at around 6.1 and 8.5 (±3.1) days and, as a comparison, spectra of SN 2011fe (normal SN Ia), iPTF13ebh (transitional SN Ia), and SN 1999by (SN Ia 1991bg-like). We also include delayed-detonation models for SN 1999by some 24−27 days after the explosion, which corresponds to about 8−11 days past-maximum light in the B band (see Sect. 6.2 for more details). The SN 2016hnk spectra are dominated in the J band by blends of singly-ionized iron group elements, Ca II, and Si II, while the H band appears very flat with the H-band break indicated in red. The two deep features at around 10 000 and 11 500 Å are only seen shallower in SN 1999by spectra.

Although in the middle of the telluric feature, the Si IIλ13650 seems to also be present at ∼13 500 Å. Finally, at ∼15 500 Å, we detect the H-band break, a complex formed by Fe II/Co II/Ni II. The strength of the H-band break correlates with Δm15(B), with brighter SNe Ia having the strongest H-band break (Hsiao et al. 2013), and the outer edge of the velocity of this region is correlated with light curve shape (Ashall et al. 2019a). The outer-edge of this feature has a velocity of vedge = 5380 ± 490 km s−1 as measured in the SN 2016hnk spectrum at +6.5 days. Sub-Chandrasekhar mass explosions are expected to have vedge values of around 7000 km s−1, while MCh models of subluminous SN Ia have lower values, around ∼5000 km s−1 (Ashall et al. 2019b). Therefore, in this picture, SN 2016hnk is compatible with high mass, near MCh models.

5.2. Nebular spectrum

The resulting stack of the six individual exposures results in a low S/N spectrum, at an average epoch of +360 days post-explosion. Carefully scanning all regions from 3000 to 25 000 Å, we find that the only visible feature is a double emission at ∼7300 Å, as reproduced in Fig. 19, which we associate to the 3 → 1 and 2 → 1 transitions of the Ca II doublet at 7291.5 and 7323.9 Å, respectively. The peaks of the emission are redshifted by ∼740 km s−1. In the top panel we show a zoom into the 2D spectra where the trace of the SN light is clearly seen at those wavelengths. A poor telluric correction is responsible for the feature around 7500 Å, as the 2D spectrum unequivocally demonstrates. Although similar narrow [Ca II] features have been previously seen in 2002cx-like SN Ia spectra (e.g., SN 2005hk, SN 2008A, and SN 2010ae; McCully et al. 2014b; Foley et al. 2016; Stritzinger et al. 2014) and are possibly also present in the 2002es-like SN Ia spectrum of SN 2010lp (Taubenberger et al. 2013; see Fig. 20), the SN 2016hnk spectrum does not show either the narrow [Fe II] λ7155 and [Ni II] λ7378, nor the broad [Fe II] λ7155 accompanying features, respectively. In addition, the lack of other features at bluer wavelengths can be attributed to the low S/N and high reddening that is affecting the SN. We provide more details in Sect. 6.3.

thumbnail Fig. 19.

Top: two-dimensional frame of the combined X-shooter nebular phase (+360 d) spectrum of SN 2016hnk after sky subtraction around the region where we detected the only feature. The trace can be clearly seen around 730 nm. The other feature around 750 nm is due to a bad telluric subtraction, since it is present along the spatial direction (vertical). Middle: combined 1D X-shooter optical spectrum showing the only two features detected in the whole UV+Opt+NIR wavelength range. Bottom: zoom-in to the feature we identified as being due to calcium emission. The blue and orange shaded regions represent the shift of ∼700 km s−1 between the peak of the emission and the rest-frame wavelength.

thumbnail Fig. 20.

Comparison of SN 2016hnk nebular spectra (in black) to other 1991bg-like, 2002es-like, 2002cx-like, and Ca-rich objects. While all these objects show Ca emission at late times, only SN 2002cx-like (and possibly 2002es-like; e.g., SN 2010lp) objects show the narrow emission lines similar to SN 2016hnk. The yellow shaded region is zoomed in on the right panel, where the vertical dotted lines represent the rest-frame wavelength of the two [Ca II] emission lines. We note that while SN 2005hk lines have negative velocities, SN 2008A lines have positive velocities similar to SN 2016hnk.

6. Spectral modeling

To get insights into the possible progenitor system and explosion mechanism that explain our observations we turn to modeling of the observed spectra. In Sect. 6.1, we want to evaluate spectral line identifications and spectral fits using the method of abundance stratification. This approach is rather independent from the explosion scenario because the overall envelope structures of most explosions scenarios are similar, and nuclear physics dominates the specific explosion energy.

In the subsequent sections, we use detailed time-dependent non-LTE models. Based on results of the abundance stratification method and guided by the observed light curves and spectra during the photospheric phase, we identify an explosion model from the literature that matches the observations of SN 2016hnk. This turns out to be the non-LTE model for SN 1999by from Höflich et al. (2002). Abundance stratification does not allow us to fit NIR spectra because of the assumption that the photosphere can be approximated as a blackbody. Therefore, for verification, we check that the non-LTE model mentioned above can reproduce the NIR spectrum of SN 2016hnk without further tuning of the parameters; we also identity the features produced in the NIR (CI, Mg II, and other IME lines). The non-LTE models are then used to revisit some of the line identifications of the previous section and compared to the abundance stratification results.

In Sect. 6.3, we modify our explosion model with the goal of reproducing the narrow [Ca II] feature during the nebular phase. We change the initial conditions of the WD and calculate the initial phase of the expansion up to about ten days (Höflich et al. 2017). We also construct detailed non-LTE nebular spectra, which include gamma- and positron transport, but we neglect radiative and recombination timescales. We discuss why energy input and transport by radioactive 56Ni decay are important in the central region.

6.1. Abundance stratification

The abundance stratification technique utilizes the fact that a SN Ia is in homologous expansion (vph(r)∝r) from ∼10 s after explosion (Röpke 2005). With a given density profile, time from explosion (texp), photospheric velocity (vph), bolometric luminosity (Lbol), and abundance structure, optimally fitting synthetic spectra are produced. This approach has been used for many SNe including subluminous and transitional SNe Ia such as SNe 2005bl, 1986G, 2007on, and 2011iv (e.g., Hachinger et al. 2009; Ashall et al. 2016b, 2018), normal SNe Ia (e.g., Tanaka et al. 2011; Ashall et al. 2014; Mazzali et al. 2014), and stripped envelope SNe (e.g., Ashall et al. 2019c; Prentice et al. 2018). A brief description of the code can be found in Appendix B.

The code requires an input density profile, in the case of SN 2016hnk we used the W7 density profile (Nomoto et al. 1984). The W7 density profile looks similar to delayed-detonation models of subluminous SNe Ia from Höflich et al. (2017). This density profile comes from a fast deflagration explosion of a Ch-mass WD. Spectral models were produced between +2.5 to +7.5 days relative to B-band maximum, to obtain the broad structure and chemical composition of the outer layers of the ejecta. We used a rise time of 16.5 days, which gives the best fit and is consistent with the value estimated for SN 2016hnk from non-detections in Sect. 2. Figure 21 shows the synthetic spectra and models at three epochs (corresponding to +19.0, +21.8, +24.0 days from explosion). All observed spectra were first corrected for host galaxy and foreground galactic extinction. The strongest lines contributing to each feature are labeled in the top panel of the plot. The bluest part of the spectra (< 3900 Å) is dominated by blends of many Fe group and metal lines, which cause a large amount of line blanketing and reprocesses the flux into the optical. However, the overall red shape of the spectrum is intrinsic and because SN 2016hnk is less luminous and cooler than a standard object.

thumbnail Fig. 21.

Time series of data (black) and spectral models (blue) for SN 2016hnk during the photospheric phase. The data were corrected for foreground and host galaxy extinction. The observer-frame time from B-band maximum and the rest-frame time from explosion is provided for each model. Line identifications from the model are shown in the top plot.

The strongest lines in the ∼3900 Å feature are the Ca II H&K lines λλ3934,3969 with some contribution from Ti IIλλ3901,3913. The feature at ∼4200 Å is dominated by Ti IIλλ4395,4418,4444,4501, Fe IIλλ4178,4233,4297,4303,4352,4417,4491, and Mg IIλ4481. Moving redward the feature at ∼5000 Å consists of Fe IIλλ4924,5018,5169,5198,5235,5276,5317, Si IIλλ5041,5056, and Mg Iλ5184. The notch at ∼5300 Å is produced by S IIλλ5432,5454 and Si IIλλ5466,5467. The absorption at ∼5450 Å consists of S IIλλ5606,5639 and Si IIλλ5670,5707, although we note that in the data there could be an emission line in this area of the spectra. The feature at ∼5800 Å is produced by Si IIλλ5958,5979 and the small absorption between the two Si II features is made by Fe IIλλ6147,6150,6238. The absorption at ∼6200 Å is produced by Si IIλλ6347,6371; the unusual red side of the feature is dominated by Fe IIλλ6456,6516, which gives it a broader asymmetrical shape. This is seen in most subluminous SNe Ia (see Fig. 14). However, for brighter objects, where the Si region is further out in velocity space, this Fe IIλ6516 line is on the emission component of the Si IIλ6355 feature, as the velocity of the Fe and Si-rich regions may not scale in by the same factor. In the latter case it may be easy to misidentify this feature as C II.

Between 6200 and 7700 Å there are a number of small absorption features not typically seen in brighter SNe Ia. Many of these are low excitation Fe II lines that are produced at lower temperatures, which is consistent with the low luminosity of SN 2016hnk. The strongest lines in the models between the Si II 6355 Å feature and the O Iλ7774 features are Fe IIλλ6456,6516,7308,7462,7711, Si IIλλ6818,6830, and Al IIλλ7042,70574. The feature at ∼7770 Å is dominated by O Iλλ7772,7774,7775 with strong contributions from Si IIλλ7849,7850 and Mg IIλλ7877,7896,7896. The deep feature at 8200 Å, often attributed to only Ca II, is produced by Ca IIλλ8498,8542,8662, O Iλλ8446,8446,8447, and Mg IIλλ8213,8235. It is important to note that the blue side of this feature consists of Mg II, Ca II, and O I, however Mg II is not one of the main sources of opacity in the feature and its contribution is minor. Finally the absorption at 8900 Å is produced by Mg IIλλ9218,9244. Throughout the spectra modeled in Fig. 21 the line identification does not change significantly; only the strength of each line varies.

Table 4 shows the abundances as a function of velocity that were used to compute the synthetic spectra. We also provided the abundances in the shells outside the formation region of the first spectrum. These shells are required to produce the stratified abundances and to match the observations. For example, if we assume the abundances in the 6500 km s−1 shell are the same as the very outer layer, the models would contain high-velocity Si absorption. This would be caused by the ejecta at high velocities containing 68% Si.

Table 4.

Input parameters for the early-time spectral models.

In the first epoch the modeled vph is 6500 km s−1, the abundances near the photosphere are dominated by Si (∼68%), S (∼24%), and C+O (7%), where Ca, Ti, Cr, V, and Fe make up the remaining 1%. The next epoch has a vph of 5000 km s−1, and the Si and S abundances have increased slightly to ∼70% and ∼29%, respectively. This model has no C or O at the photopshere; Ca, Ti, Cr, V, and Fe make up the rest of the mass. By the third epoch the vph is 4700 km s−1 and the Si abundance is ∼72%, S is ∼27%, and the remaining mass consists of Ca, Ti, Cr, V, and Fe. For a summary of errors using this method, see Mazzali et al. (2008).

The derived abundance distribution obtained through spectral modeling is consistent with a low luminosity delayed-detonation SN Ia. Our models are sensitive out to velocities of ∼15 000 km s−1. Intermediate mass elements dominate the spectra at the epochs observed, however the presence of significant O is required above 9000 km s−1. It can be assumed that the bulk of the 56Ni is located below our innermost shell. There could be a small abundance of 56Ni above the innermost photosphere, but it is not large enough to be optically thick. There is less effective burning in this SN explosion compared to a more luminous event.

In these models at +24 d past explosion there are 0.8 M above the photosphere. As there is a large amount of material below the photosphere, it is not possible for us to constrain the total mass of the explosion. However the velocity of the H-band break (see above) and the nebular phase spectrum helps; see Sect. 6.2.

A comparison of Fig. 22 to the abundance distributions of other transitional and subluminous SNe Ia show similar results. However, for SN 2016hnk we find no evidence of 56Ni in our photospheric models, whereas in the subluminous SN 2005bl Hachinger et al. (2009) found 56Ni out to ∼7000 km s−1. A comparison to the abundance distribution to the transitional SN 1986G (Ashall et al. 2016b) shows similar results to SN 2016hnk; the only discrepancy is that SN 1986G has 56Ni located out to 5500 km s−1. As discussed below, we note that our abundance distribution is similar to the delayed-detonation model of SN 1999by of Höflich et al. (2002).

thumbnail Fig. 22.

Abundance distribution as a function of velocity for SN 2016hnk. The vertical dashed lines denote the photosphere velocity of the +2.5 and +7.5 d models. We note that because of the small time step between our first and last model our results are only sensitive between 4500 and 10 000 km s−1.

6.2. Explosion models photospheric phase

We calculated the final phase of the thermonuclear runaway, the explosion, light curves, and spectra with our code HYDRA Höflich et al. (2017), and references therein. We used delayed-detonation models Khokhlov (1989) in which the explosion starts as a slow laminar deflagration. The free parameters are the main sequence mass MMS and metallicity Z of the progenitor, the central density of the explosion ρc. The flame starts as laminar deflagration and transitions to a detonation at a density ρtr, which in spherical geometry is equivalent to the amount of burning during the deflagration phase. To set up the explosion, we calculated the final phase of the accretion toward the thermonuclear runaway with our code but neglected the second derivatives to allow for large time steps similar to Nomoto (1982). In this case, we adjusted the constant accretion rate so that the thermonuclear runaway occurs at a given central density ρc (Höflich et al. 2017). To first order, the explosion is triggered when compressional heating due to increasing mass by accretion from the companion star cannot be balanced by heat conduction and neutrino cooling (including the convective URCA process). Close to the Accretion Induced Collapse (AIC), the detailed parameters such as the accretion rate depend sensitively on the equation of state and neutrino cooling, and the history of the binary evolution (see final discussion and conclusions).

As mentioned in Sect. 6.1, SN 2016hnk has many observational characteristics similar to the subluminous SN 1999by. To produce subluminous SNe Ia in the framework of DDT models, the pre-expansion of the WD needs to be sufficiently large during the deflagration phase that the density has dropped enough so that little 56Ni is produced during the subsequent detonation phase, and a partial suppression of Rayleigh-Taylor (RT) mixing is needed (Höflich et al. 2002, 2017; Hristov et al. 2018). The model for SN1999by had ρc  =  2 × 109 g cm−3, ρtr  =  8 × 106 g cm−3, and MS = 5 M⊙, Z = Z. Our tomography approach for the early spectra resulted in an abundance structure that closely resembled that of the delayed-detonation model for SN1999by, which shows that the Nuclear Statistical Equilibrium (NSE) elements are within ≈4500 km s−1, followed by a Si/S-rich layer and a transition to layers of explosive carbon burning (Mg/O) at ≈9000 km s−1.

Near-infrared observations provide a useful constraint on the explosion because its sensitivity to whether the photosphere is formed in a mostly C/O rich mixture, IME elements, or NSE elements. As such, these observations provide a useful probe for the time of the explosion and the explosion model. In Fig. 18 a NIR spectrum of SN 2016hnk is compared to the non-LTE delayed-detonation model of SN 1999by at 26 days past explosion, corresponding to about 11 days past maximum. Overall, the observed and theoretical spectra agree. The theoretical spectrum is dominated by blends of singly ionized Si, Ca, Fe, and Co. The strongest features are Fe II (≈1.02 μm), Si II (≈1.20 μm), Fe II (≈1.06 μm), Ca II, Fe II and Co II (1.16 μm), Co II and Fe II (1.15 and 1.38 μm), and a multiplet of Fe II/Co II/NiII, which is weak but appears in the H-band region. We note that the feature at 0.98 μm is a Fe II dominated blend of Ca II/Fe II/Co II. In the model it is more defined or in the observation the suppression of flux is stronger. Because it is formed in a steeply declining wing, it is ambiguous whether the metallicity, Z, of SN 2016hnk is significantly larger than Z in the model leading to a suppression of flux and a widening of the adjoining feature to form a shoulder, or if Z is significantly smaller than Z resulting in a narrower feature. This strongly places SN 2016hnk into the realm of the subluminous 1991bg-like SNe Ia and emphasizes the NIR as a chronometer for the time since maximum because of the rapid spectral changes. The time of this model is consistent with our light curve-based estimates for maximum light and with abundance stratification models.

The features between 6000 and 7500 Å in the photospheric phase spectra are not normally seen in bright SNe Ia. With the abundance stratification models we produce many of these lines but at different flux levels. Therefore, to see if our line identifications are correct, we compare SN 2016hnk to the full non-LTE models of SN 1999by from Höflich et al. (2002).

To conclude the photospheric section, the element identification should be revisited in light of the explosion based non-LTE models. Figure 23 shows the +2.5 ± 3.3 d spectra of SN 2016hnk and the non-LTE model of SN 1999by at +7 d in the wavelength region 6200−7800 Å. We note that we limited our analysis to the published spectra from Höflich et al. (2002). We note that in the model at +2 d these lines have not yet formed, but given the error on the time of maximum from the observations and the quick evolution of SN 2016hnk, we speculate that these lines begin to emerge in the models of SN 1999by between +2 and +7 d, which is consistent with the observations and NIR models. The model and data have similar spectral features, shown with the vertical dashed lines. The line identification is mainly consistent with the abundance stratification models, where Fe II lines dominate the spectra in the 6000−7500 Å region.

thumbnail Fig. 23.

Comparison of SN 2016hnk (black) +2.5 d past maximum light and the non-LTE model of SN 1999by at +7 days past maximum light (blue) from Höflich et al. (2002). The telluric features are denoted in the observed spectra.

Some of the other features in the 6200−7800 Å region are Fe IIλ6664 and Co Iλ7084. The small absorption at 6550 Å is one feature that cannot be explained by the models. However, the Höflich et al. (2002) models assume no mixing. With a small amount of mixing it could be expected that this feature is caused by Co IIλλ6571,6576, as mixing enhances the Co abundance in these layers. In the wavelength region 6000−7500 Å we favor the line identifications from the non-LTE models, as these solve for atomic level populations explicitly and are produced through full radiation-hydrodynamical simulations and not a parameterized to fit the data.

6.3. Explosion models for the nebular phase

Owing to the similarity of SN 2016hnk to SN 1999by, we adapted the previous MCh mass spectral simulations of SN 1999by (Höflich et al. 2002, 2017) to produce models for SN 2016hnk. In SN 2016hnk an interesting observation in the nebular spectrum is the appearance of a single, narrow dominating spectral feature at ≈7200 Å that has a Doppler width of ≈1000 km s−1. This is very different from late-time spectra of normal bright SNe Ia that are dominated by iron-group elements. Because the delayed-detonation model for SN 1999by does not produce narrow features on less than a 2000 − 3000 km s−1 scale, the inner region of the original model needs to be modified.

We note that, in principle, multidimensional effects such as RT instabilities produce plumes with an individual velocity spread on the scale of ≈1000 km s−1. However they produce multiple plumes, which would make the total velocity spread too large to explain the observation of SN 2016hnk. Moreover, various observations and analyses strongly suggest that the RT instabilities are present but need to be partially suppressed (Höflich et al. 2006, 2017; Fesen et al. 2015; Diamond et al. 2018; Hristov et al. 2018). Therefore, in this work spherical delayed-detonation models are used, and we focus on nuclear physics based effects in combination with progenitor parameters to explain the observations of SN 2016hnk. The hydrodynamical simulations are in the comoving (mass) frame to avoid numerical advection of abundances, which may affect both nuclear burning and small-scale structures. The simulations utilize a nuclear network of 315 isotopes during the early phase of the explosion based on the nuclear reaction network using the library REACLIB (Cyburt et al. 2010) and weak rates from Martínez-Pinedo et al. (2014), Langanke & Martínez-Pinedo (2003), and Fuller et al. (1985). We use detailed, time-dependent non-LTE models for atomic level populations, including γ-ray and positron transport and radiation-hydrodynamics (RH; Höflich 1995, 2003; Penney & Höflich 2014; Höflich et al. 2017). After an initial deflagration phase, a detonation is triggered when the density at the flame front drops to a transition density, ρtr, which is a parameterization of the amount of deflagration burning.

The model presented originates from a progenitor with a main sequence mass (MMS) of 7 M with solar metallicity. The models of SN 2016hnk require a higher central density (ρc  =  6 × 109 g cm−3) than that of SN 1999by, which implies a larger binding energy, and to compensate for the longer phase to lift the WD, we used a slightly higher ρtr  =  107 g cm−3. The rate of subluminous SNe Ia diminish at high redshifts indicating a long evolutionary timescale (González-Gaitán et al. 2011). Thus, we assume an increased 22Ne/20Ne ratio (≈ factor of 2) in the center of the WD because of gravitational settling over billions of years (Deloye & Bildsten 2002). We note that most of the Ne is inherited from the formation of the progenitor star in the molecular cloud. In effect, our model starts at a slightly lower Ye compared to models without gravitational settling. The front starts with a laminar speed. To increase the electron capture, we assume a slow deflagration front (C1  =  0.1 of Eq. (1) in Domínguez & Höflich 2000).

The density and velocity structure of the model during the homologous expansion, after ≈10 days, is shown in Figs. 24 and 25. The density structure is very similar to the model used for SN 1999by in most respects. As in Höflich et al. (2002) models and also in W7 (Nomoto et al. 1984), a local maximum in the outer density structure, M ≈ 0.98 MWD, marks the transition from burned to unburned material. The feature at M ≈ 0.15 MWD marks the imprint of the refraction wave during the deflagration to detonation transition, which is inherent to spherical models. The increased density in the central region of the explosion5 is not present in the Höflich et al. (2002) models, as it is caused by radiation hydrodynamical effects. After the explosion, the envelope enters the phase of “almost” free expansion on timescales of seconds to minutes. Adiabatic cooling results in a rapid drop of temperature in the envelope and an overall free expansion. This explains the similarity in the outer layers between the new model and the old SN 1999by model. However, the energy production by radioactive 56Ni decay corresponds to a change of the specific energy equivalent to a velocity of ≈3000 km s−1. This effect changes the density structure of low-velocity layers because the “dynamics” is dominated by pressure equilibrium. In models with lower ρc, the central region becomes partially optically thin for γ-rays and low energy photons. Thus, the entire inner region is affected by the heating. However, in our high-density model with ρc = 6 × 109 g cm−3, the average γ-ray optical depth τg stays larger than 2−3, even at twice the decay time of 56Ni, leading to a compression of layers with ≲1500 km s−1 (Fig. 24). We note that these layers are well inside the 56Ni layers as discussed below. However, this increase has severe effects on the γ-ray heating because the central τg is an order of magnitude larger than the rest of the envelope even at +350 d. In part, this is responsible for the late-time spectral formation discussed below.

thumbnail Fig. 24.

Structure during the phase of homologous expansion as a function of the mass coordinate in WD-masses for our subluminous model with high central density. We give the normalized density as log(ρ/ρc) (left scale), velocity (right scale), and the mean optical depth log(τγ) (left scale) in γ-rays at day 350.

6.3.1. Ca line in the late-time spectrum

The main effect on Ye is the density of burning. A central density of 6 × 109 g cm−3 is required to reduce Ye to values that produce centrally located Ca as observed. With regard to the electron capture, the settling of Ne and slow deflagration phase has the same effect of an increased ρc by ≈10%. The transition to a detonation is triggered at ρtr = 107 g cm−3 to compensate for the larger binding energy of the initial WD. The WD has a final explosion energy Ekin = 8.83 × 1050 erg and produces M(56Ni) = 0.108 M. The overall density and chemical structure is given in Fig. 25, and the element production is shown in Table 5. The overall chemical structure resembles the Höflich et al. (2002) model for SN 1999by and, more importantly, the results of the abundance tomography analysis (Table 4): (a) the unburned C-O layers reach down to about 15 000 km s−1; (b) products of explosive carbon burning (O-Ne-Mg) are seen down to ≈9500 km s−1; and (c) followed by incomplete burning (Si-S) down to about 3500 km s−1.

thumbnail Fig. 25.

Same as Fig. 24 but showing the overall abundance structure at large times as a function of mass coordinate (top) and expansion velocity (bottom). It is very similar to the subluminous model for SN 1999by and consistent with the chemical structure based on our abundance tomography given in Table 3. In addition, we show the initial 56Ni → 56Co → 56Fe at the time of the explosion (purple dotted line).

Table 5.

Asymptotic isotopic yields at large times for the delay-detonation model for SN 2016hnk discussed in Sect. 6.2.

In the inner region, electron capture shifts the nuclear statistic equilibrium from 56Ni to neutron rich isotopes, as the Ye decreases toward the central region; see Fig. 26. The decrease is from a Ye of 0.497, which produces 56 − 58Ni, 55Co, 54Fe, to a Ye of 0.43 − 0.44, which produces 58 − 60Fe, 54 − 56Mn, 52 − 56Cr and most importantly 47 − 48Ca.

thumbnail Fig. 26.

Distribution of most abundant isotopes in log(Xi) and Ye as a function of velocity for the central region are dominated by neutron-rich isotopes at several minutes after the explosion. Subsequently, radioactive isotopes decay 55Cr (3.5 m) → 55Mn, 55Fe(2.7 y) → 55Mn, 59Fe (44.5 d) → 59Co, 56Ni (6 d) → 56Co (77 d) → 56Fe, 57Co (271 d) → 57Fe.

In Fig. 27, the comparison between the theoretical and observed late-time spectrum is given. The observation shows narrow emission from the forbidden [Ca II] doublet at λλ7291,7324 whose lower level is the ground-state of Ca II. This narrow [Ca II] doublet dominates with respect to the other features because most of the energy emitted in Ca is confined to a very narrow region of central Ca powered by radioactive decay. The strength of these other features is only up to the 5% of the Ca doublet and is lost below the noise level. The central energy deposition is large because that region is still optically thick to γ-rays, whereas γ-rays from the outer layers mostly escape (Fig. 24). The underlying continuum is formed by many overlapping allowed and forbidden transitions of mostly Cr, Mn, and Fe group elements, which distribute the energy over a wide wavelength range from the UV to mid-infrared and produce many lines commonly seen in SNe Ia, but on much lower levels of about ∼3% relative to the strength of the dominant narrow Ca feature. We note that the density at the inner region is still high at ≈3 × 107 particles cm−3. Thus, forbidden lines dominate for Ca for transitions to the ground state but not for the more complex ions. In the model, the central feature is formed by 48Ca with expansion velocities vexp up to about 700 km s−1, whereas a shoulder is the result of 40Ca of the region expanding at v ≈ 3000 − 9000 km s−1. Within the scatter in the observed spectrum (black, Fig. 27), the synthetic spectra agree with the observations.

thumbnail Fig. 27.

Comparison of the nebular spectrum of SN 2016hnk with our high-density delayed-detonation model in the rest frame assuming a redshift of z = 0.01610. The flux is in 10−18 erg s−1 cm−2 on the same scale as Fig. 19. The telluric region in the original spectrum is plotted in gray, around 7500 Å. The theoretical spectrum is dominated by a narrow, forbidden [Ca II] doublet at 7291 and 7324 Å above a quasi-continuum at about the 3% level formed by many transitions of mostly heavier elements including Fe, Co, Ni, and Cr.

The comparison with the smoothed spectrum (green, Fig. 27) shows some discrepancies, namely the hint of the 40Ca shoulder in the models as just discussed and, in particular, a Doppler shift of the central Ca feature by about 700 km s−1. Possible factors that could contribute to this offset and that can be discriminated by future observations of high-resolution spectra include a peculiar motion of the SN progenitor within the galaxy or a passive drag of preexisting turbulent fields during the smoldering phase (Höflich & Stein 2002; Yoon & Langer 2004; Zingale & Dursi 2007; Zingale et al. 2011) during the 2−3 s of a typical deflagration phase. The former would show variations but would be correlated with the rotation of the host galaxy, while the latter would show variations between objects. We note that the smoothed spectrum shows structure on the 10% level that is not predicted, but is also not at a significant level. The strong jump at 7500 Å is telluric.

7. Discussion and conclusions

We present photometric and spectroscopic observations of SN 2016hnk and IFS of its type SB host galaxy MCG -01-06-070. For the first time, these data are combined with a unique late-time spectrum at ≈360 d and with IFS data providing a unique link to the progenitor and environment of subluminous SNe Ia and, possibly, their nature. The IFS places SN 2016hnk close to a region with significant Hα emission and high dust content near the nucleus of MCG -01-06-070. For the ISM, we consistently find E(B − V)∼0.4 mag with three different methods and a low RV = 2.1 ± 0.6 typical of highly extinguished SNe Ia (Burns et al. 2014).

The SN 2016hnk light curve is most similar to 2002es-like SNe Ia up to +15 days past maximum, however, optical and NIR spectra of SN 2016hnk during the photospheric phase show that the SN closely resembles the prototypes for subluminous SNe Ia, namely SNe 1991bg and 1999by. The Δm15, s(B) of 1.80 ± 0.20 mag and sBV of 0.438 ± 0.030 confirms the similarity between SN 2016hnk and subluminous SNe Ia. Based on detailed analyses of the optical and NIR light curves, we detect a flux excess after ≈15 d caused by scattering processes with an interstellar medium at around 1.0 pc from the SN, suggesting that the subluminous SN 2016hnk is embedded in a dense environment. We note that the extinction by dust is significantly larger in our object than in SNe 1991bg, 1999by and 2005ke (Lucy 1999; Höflich et al. 2002; Patat et al. 2012; Penney & Höflich 2014). The underrepresentation of subluminous SNe Ia in high-redshift samples is usually attributed to long progenitor evolution, but may also be magnified by subluminous SNe Ia occurring in dusty regions or at low galactocentric distances (Galbany et al. 2012), or simply be attributed to selection effects (González-Gaitán et al. 2011).

The early-time spectra of SN 2016hnk have unique absorption features between 6000 and 7000 Å, which we identified as low excitation Fe II, Co II, and Co I lines. Furthermore, it was shown that the line on the redward side of the Si II 6355 Å feature is not produced by C II but rather by Fe IIλ6516. In the instance of subluminous SNe Ia the ionization state of carbon at the epochs studied would be neutral and carbon would not be present in the optical spectrum, furthermore it could be expected to be further out in velocity space and not visible at these epochs. This Fe II line has also been seen by a few groups for normal bright SNe Ia (see Höflich 1995) and with the PHOENIX code (Baron, priv. comm.).

We performed a detailed comparison of the optical spectra to spectral abundance stratification and tomography. The results are consistent with a MCh mass WD with characteristics found in a previous analysis of SN 1999by based on delayed-detonation models (Höflich et al. 2002). These characteristics include the following: (a) low maximum-light luminosity and ionization state; (b) low photospheric velocities; (c) low-velocity iron-group elements; and (d) an extended region (in velocity space) of products of explosive C-burning and incomplete O- and Si-burning. The model spectrum places our early-time NIR spectrum at +11 days after maximum light, which is ≈2 days later than inferred from maximum-light curve fits. Given the fast evolution in the NIR and the large uncertainty in tmax, we show that NIR spectra can be used efficiently to refine the time of maximum light and place the first optical spectrum at about two to three days after maximum. Both the tomography and the previous analysis of SN 1999by provide consistent line identifications for spectra dominated by neutral and single-ionized Fe, S/Si II, and O I.

Moreover, and unlike other SN Ia subtypes, there are indications of asphericity for subluminous SNe Ia. The early continuum polarization of SN 1999by shows a good agreement with the same models from Höflich et al. (2002) used in this work when an asphericity of about ∼20% is included (Howell et al. 2001). However, the resulting anisotropy in the luminosity is about ∼5% and hardly effects the spectra (Patat et al. 2012).

The late-time nebular spectrum obtained by X-shooter shows only one prominent, narrow feature with a width of ≈700 km s−1. This feature has not been predicted by the delayed-detonation model for SN 1999by even though most of the optical and NIR features at early phases are rather similar between SN 2016hnk and SN 1999by, as discussed above. This narrow late-time feature can be attributed to the [Ca II] doublet at λλ7291,7324. It can be understood as a result of a nucleosynthesis effect, namely electron capture. High central densities of a WD close to the limit for an accretion induced collapse, with Ye less than 0.44, results in a shift of the nuclear statistical equilibrium toward 48Ca. In our model, a separate Ca layer is produced at expansion velocities less than 750 km s−1. As discussed in Sect. 6.2, the exact density depends on the initial conditions. For example, the gravitational settling of heavy elements (i.e., 22Ne and 20Ne) over long timescales may increase the central density by about 10% to achieve lower Ye. There is a high ratio of 48Ca/40Ca in solar composition (Thielemann et al. 2002), and the result may be suggestive of some subluminous SNe Ia as a main production sites for 48Ca. However, we have late-time, high-resolution observations only for SN 2016hnk and for none of the other subluminous SNe Ia, and we lack high S/N observations of the underlying iron-group lines that are predicted on the few percent level compared to the [Ca II] feature. We note that sub-MCh models seem to fall well short of producing 48Ca (Brachwitz et al. 2000), despite having a compressional detonation wave, which increases the peak density in a WD (Höflich et al. 2018).

We note our interpretation of the narrow nebular [Ca II] feature may not be unique. However, RT instabilities are unlikely to cause this feature. As discussed in Sect. 6.2, RT instabilities would produce many plumes rather than a single narrow feature and, even if one plume were to dominate from the region of incomplete burning, we would expect a mixture of Ca/Fe and, potentially, S/Si. Conversely, a narrow feature similar to the observations could be produced within the ejecta in a narrow jet-like Ca-rich structure. In this scenario, the single strong feature does not consist of iron-group elements and must be viewed almost equatorially. If viewed pole on we would not see a narrow feature, but rather a high-velocity feature.

Another explanation could be in the framework of interacting binaries with material stripped from the donor star. This can occur within both the MCh and sub-MCh scenarios. In principle, the low-velocity Ca may be attributed to material stripped from a companion star. However, in this case, we may expect very strong forbidden [O I] lines at 6300, 6364 (and 6393) Å regardless of whether the accreted material is H, He, or C-O rich. To rule out or confirm this possibility we would require a time series of spectra of more objects at those late epochs.

Finally, we want to address the possible link to progenitor systems. Potential progenitor systems may either consist of systems of two WDs (called double degenerate), a single WD with a donor which may be main sequence, red giant, or He star (called a single degenerate), or explosions within an asymptotic red giant branch star (Branch et al. 1995; Han et al. 2010; Piersanti et al. 2014; Denissenkov et al. 2015; Nomoto & Leung 2018; Yoon et al. 2007).

However, the link of SN2016hnk is highly uncertain. As discussed, SN2016hnk is in a galaxy of low SFR. This may point toward a progenitor scenario with a long evolutionary timescale, but there is still some ongoing star formation. This problem cannot be solved by direct analysis of the ejecta due to the lack of “classical” observations such as high-precision light curves, as well as flux and polarization spectra well before maximum light (Patat et al. 2012; Hoeflich et al. 2013; Gall et al. 2018b). The lack of early time light curves, optical, or IR spectra do not allow us to probe the outermost layers to measure Z in the unburned layers, or the amount of stable iron-group elements.

This leaves the high central density of the WD and the nearby material in the environment as a hint toward the progenitor. Within our models, we need a low accretion rate during the final phase of accretion to achieve a high central density, which results in a low Ye. The nearby material may indicate a high mass loss rate prior to the explosion, for example, during a common envelope phase. For many of the common progenitor systems, a drop in the final accretion rate may be expected because binary systems tend to become wider due to mass loss. However, the low accretion rate depends on the treatment of the convective URCA process (Paczyński 1972; Barkat & Wheeler 1990; Stein & Wheeler 2006; Denissenkov et al. 2015; Piersanti et al. 2017), and a faster accretion rate may automatically lead to an explosion close to the AIC.


2

SNOOPY is a package for SN photometry using PSF fitting and/or template subtraction developpd by E. Cappellaro. A package description can be found at http://sngroup.oapd.inaf.it/snoopy.html

4

We note that we doubt the Al II line identification in this case as it requires an abundance that is larger than what is predicted by explosion models.

5

We note this is different from the central density of the exploding WD.

Acknowledgments

Endavant les atxes. LG would like to thank the organizers of the Transients in New Surveys: the Undiscovered Country workshop held in July 2018 in Leiden, where this work started as a result of a huge collaborative effort, and the Lorentz Center staff for providing such a great work environment. CA and EYH acknowledge the support provided by the National Science Foundation under grant No. AST-1613472. PH acknowledges support by the National Science Foundation (NSF) grant 1715133, and thanks F.K. Thielemann and G. Martinez-Pinedo for helpful discussions and advice on the current state of electron-capture rates. SGG acknowledges support by FCT under Project CRISP PTDC/FIS-AST- 31546. MS is supported in part by a generous grant (13261) from VILLUM FONDEN and a project grant from the Independent Research Fund Denmark. EB acknowledges partial support from NASA grant NNX16AB25G, and thanks the Aarhus University Research Fund (AUFF) for a guest researcher grant. SB is partially supported by China postdoctoral science foundation grant No.2018T110006. MB acknowledges support from the Swedish Research Council (Vetenskapsrådet), the Swedish National Space Board, and the research environment grant “Gravitational Radiation and Electromagnetic Astrophysical Transients (GREAT)”. CPG and MS acknowledge support from EU/FP7-ERC grant number 615929. JH acknowledges financial support from the Finnish Cultural Foundation. DAH, CM, and GH were supported by NSF grant AST-1313484. Support for TdJ have been provided by US NSF grant AST-1211916, the TABASGO Foundation, Gary and Cynthia Bengier. SWJ acknowledges support from NSF award AST-1615455. KM acknowledges support from STFC (ST/M005348/1) and from H2020 through an ERC Starting Grant (758638). Research by DJS is supported by NSF grants AST-1821987 and 1821967. XW is supported by the National Natural Science Foundation of China (NSFC grants 11325313 and 11633002) JZ is supported by the National Natural Science Foundation of China (NSFC, grants 11773067, 11403096), the Youth Innovation Promotion Association of the CAS (grants 2018081), and the Western Light Youth Project of CAS. Based in part on data obtained with the LCO network under programs CLN2016A-005 (PI: Prieto), NAOC2017AB-001 (PI: Dong), CON2016A-008 (PI: Stanek), and ARI2017AB-001 (PI: Bersier). This work has made use of data from the Asteroid Terrestrial-impact Last Alert System (ATLAS) project. ATLAS is primarily funded to search for near earth asteroids through NASA grants NN12AR55G, 80NSSC18K0284, and 80NSSC18K1575; byproducts of the NEO search include images and catalogs from the survey area. The ATLAS science products have been made possible through the contributions of the University of Hawaii Institute for Astronomy, the Queen’s University Belfast, and the Space Telescope Science Institute. ASAS-SN is supported by the Gordon and Betty Moore Foundation through grant GBMF5490 to the Ohio State University and NSF grant AST-1515927. Development of ASAS-SN has been supported by NSF grant AST-0908816, the Mt. Cuba Astronomical Foundation, the Center for Cosmology and AstroParticle Physics at the Ohio State University, the Chinese Academy of Sciences South America Center for Astronomy (CASSACA), the Villum Foundation, and George Skestos. We thank the Las Cumbres Observatory (LCO) and its staff for its continuing support of the ASAS-SN project. NUTS is supported in part by the Instrument Center for Danish Astrophysics (IDA). This research uses data obtained through the Telescope Access Program (TAP), which has been funded by “the Strategic Priority Research Program- The Emergence of Cosmological Structures” of the Chinese Academy of Sciences (Grant No.11 XDB09000000), the Special Fund for Astronomy from the Ministry of Finance, and the National Astronomical Observatories of China. Funding for the LJT has been provided by the Chinese Academy of Sciences (CAS) and the People’s Government of Yunnan Province. Based on observations obtained at the Gemini Observatory under program GS-2016B-Q-22 (PI: Sand). Gemini is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the NSF (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). The data were processed using the Gemini IRAF package. We thank the queue service observers and technical support staff at Gemini Observatory for their assistance.

References

  1. Abolfathi, B., Aguado, D. S., Aguilar, G., et al. 2018, ApJS, 235, 42 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aguado, D. S., Ahumada, R., Almeida, A., et al. 2019, ApJS, 240, 23 [NASA ADS] [CrossRef] [Google Scholar]
  3. Amanullah, R., & Goobar, A. 2011, ApJ, 735, 20 [NASA ADS] [CrossRef] [Google Scholar]
  4. Arnett, W. D. 1982, ApJ, 253, 785 [NASA ADS] [CrossRef] [Google Scholar]
  5. Ashall, C., Mazzali, P., Bersier, D., et al. 2014, MNRAS, 445, 4427 [NASA ADS] [CrossRef] [Google Scholar]
  6. Ashall, C., Mazzali, P., Sasdelli, M., & Prentice, S. J. 2016a, MNRAS, 460, 3529 [NASA ADS] [CrossRef] [Google Scholar]
  7. Ashall, C., Mazzali, P. A., Pian, E., & James, P. A. 2016b, MNRAS, 463, 1891 [NASA ADS] [CrossRef] [Google Scholar]
  8. Ashall, C., Mazzali, P. A., Stritzinger, M. D., et al. 2018, MNRAS, 477, 153 [NASA ADS] [CrossRef] [Google Scholar]
  9. Ashall, C., Hsiao, E. Y., Hoeflich, P., et al. 2019a, ApJ, 875, L14 [NASA ADS] [CrossRef] [Google Scholar]
  10. Ashall, C., Hoeflich, P., Hsiao, E. Y., et al. 2019b, ApJ, 878, 86 [NASA ADS] [CrossRef] [Google Scholar]
  11. Ashall, C., Pian, E., Mazzali, P. A., et al. 2019c, MNRAS, 487, 5824 [NASA ADS] [CrossRef] [Google Scholar]
  12. Barbon, R., Benetti, S., Cappellaro, E., Rosino, L., & Turatto, M. 1990, A&A, 237, 79 [NASA ADS] [Google Scholar]
  13. Barkat, Z., & Wheeler, J. C. 1990, ApJ, 355, 602 [NASA ADS] [CrossRef] [Google Scholar]
  14. Becker, A. 2015, Astrophysics Source Code Library [record ascl:1504.004] [Google Scholar]
  15. Benetti, S., Cappellaro, E., Mazzali, P. A., et al. 2005, ApJ, 623, 1011 [NASA ADS] [CrossRef] [Google Scholar]
  16. Betoule, M., Kessler, R., Guy, J., et al. 2014, A&A, 568, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Blondin, S., Matheson, T., Kirshner, R. P., et al. 2012, AJ, 143, 126 [NASA ADS] [CrossRef] [Google Scholar]
  18. Blondin, S., Dessart, L., Hillier, D. J., & Khokhlov, A. M. 2017, MNRAS, 470, 157 [NASA ADS] [CrossRef] [Google Scholar]
  19. Brachwitz, F., Dean, D. J., Hix, W. R., et al. 2000, ApJ, 536, 934 [NASA ADS] [CrossRef] [Google Scholar]
  20. Branch, D., Livio, M., Yungelson, L. R., Boffi, F. R., & Baron, E. 1995, PASP, 107, 1019 [NASA ADS] [CrossRef] [Google Scholar]
  21. Branch, D., Dang, L. C., Hall, N., et al. 2006, PASP, 118, 560 [NASA ADS] [CrossRef] [Google Scholar]
  22. Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031 [NASA ADS] [CrossRef] [Google Scholar]
  23. Bulla, M., Goobar, A., & Dhawan, S. 2018a, MNRAS, 479, 3663 [NASA ADS] [CrossRef] [Google Scholar]
  24. Bulla, M., Goobar, A., Amanullah, R., Feindt, U., & Ferretti, R. 2018b, MNRAS, 473, 1918 [NASA ADS] [CrossRef] [Google Scholar]
  25. Burns, C. R., Stritzinger, M., Phillips, M. M., et al. 2011, AJ, 141, 19 [NASA ADS] [CrossRef] [Google Scholar]
  26. Burns, C. R., Stritzinger, M., Phillips, M. M., et al. 2014, ApJ, 789, 32 [NASA ADS] [CrossRef] [Google Scholar]
  27. Burns, C. R., Parent, E., Phillips, M. M., et al. 2018, ApJ, 869, 56 [NASA ADS] [CrossRef] [Google Scholar]
  28. Buzzoni, B., Delabre, B., Dekker, H., et al. 1984, The Messenger, 38, 9 [NASA ADS] [Google Scholar]
  29. Cannizzaro, G., Fraser, M., Benetti, S., et al. 2016a, NUTS Transient Classification Report for 2016-11-02, Transient Name Server Classification Report, No. 2016-872 [Google Scholar]
  30. Cannizzaro, G., Fraser, M., Benetti, S., et al. 2016b, ATel, 9703 [Google Scholar]
  31. Cao, Y., Kulkarni, S. R., Howell, D. A., et al. 2015, Nature, 521, 328 [NASA ADS] [CrossRef] [Google Scholar]
  32. Cid Fernandes, R., Mateus, A., Sodré, L., Stasińska, G., & Gomes, J. M. 2005, MNRAS, 358, 363 [NASA ADS] [CrossRef] [Google Scholar]
  33. Contreras, C., Hamuy, M., Phillips, M. M., et al. 2010, AJ, 139, 519 [NASA ADS] [CrossRef] [Google Scholar]
  34. Cyburt, R. H., Amthor, A. M., Ferguson, R., et al. 2010, ApJS, 189, 240 [NASA ADS] [CrossRef] [Google Scholar]
  35. Deloye, C. J., & Bildsten, L. 2002, ApJ, 580, 1077 [NASA ADS] [CrossRef] [Google Scholar]
  36. Denissenkov, P. A., Truran, J. W., Herwig, F., et al. 2015, MNRAS, 447, 2696 [NASA ADS] [CrossRef] [Google Scholar]
  37. DES Collaboration (Abbott, T. M. C., et al.) 2019, ApJ, 872, L30 [NASA ADS] [CrossRef] [Google Scholar]
  38. Diamond, T. R., Höflich, P., Hsiao, E. Y., et al. 2018, ApJ, 861, 119 [NASA ADS] [CrossRef] [Google Scholar]
  39. Dimitriadis, G., Pursiainen, M., Smith, M., et al. 2016, ATel, 9704 [Google Scholar]
  40. Domínguez, I., & Höflich, P. 2000, ApJ, 528, 854 [NASA ADS] [CrossRef] [Google Scholar]
  41. Dopita, M. A., Kewley, L. J., Sutherland, R. S., & Nicholls, D. C. 2016, Ap&SS, 361, 61 [NASA ADS] [CrossRef] [Google Scholar]
  42. Fesen, R. A., Höflich, P. A., & Hamilton, A. J. S. 2015, ApJ, 804, 140 [NASA ADS] [CrossRef] [Google Scholar]
  43. Filippenko, A. V., Richmond, M. W., Branch, D., et al. 1992, AJ, 104, 1543 [NASA ADS] [CrossRef] [Google Scholar]
  44. Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
  45. Folatelli, G., Morrell, N., Phillips, M. M., et al. 2013, ApJ, 773, 53 [NASA ADS] [CrossRef] [Google Scholar]
  46. Foley, R. J., Jha, S. W., Pan, Y.-C., et al. 2016, MNRAS, 461, 433 [NASA ADS] [CrossRef] [Google Scholar]
  47. Förster, F., González-Gaitán, S., Folatelli, G., & Morrell, N. 2013, ApJ, 772, 19 [NASA ADS] [CrossRef] [Google Scholar]
  48. Friedman, A. S., Wood-Vasey, W. M., Marion, G. H., et al. 2015, ApJS, 220, 9 [NASA ADS] [CrossRef] [Google Scholar]
  49. Fuller, G. M., Fowler, W. A., & Newman, M. J. 1985, ApJ, 293, 1 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  50. Galbany, L., Miquel, R., Östman, L., et al. 2012, ApJ, 755, 125 [NASA ADS] [CrossRef] [Google Scholar]
  51. Galbany, L., Stanishev, V., Mourão, A. M., et al. 2014, A&A, 572, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Galbany, L., Stanishev, V., Mourão, A. M., et al. 2016a, A&A, 591, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Galbany, L., Anderson, J. P., Rosales-Ortega, F. F., et al. 2016b, MNRAS, 455, 4087 [NASA ADS] [CrossRef] [Google Scholar]
  54. Galbany, L., Anderson, J. P., Sánchez, S. F., et al. 2018, ApJ, 855, 107 [NASA ADS] [CrossRef] [Google Scholar]
  55. Gall, C., Stritzinger, M. D., Ashall, C., et al. 2018a, A&A, 611, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Gall, C., Stritzinger, M. D., Ashall, C., et al. 2018b, A&A, 611, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Ganeshalingam, M., Li, W., Filippenko, A. V., et al. 2010, ApJS, 190, 418 [NASA ADS] [CrossRef] [Google Scholar]
  58. Ganeshalingam, M., Li, W., Filippenko, A. V., et al. 2012, ApJ, 751, 142 [NASA ADS] [CrossRef] [Google Scholar]
  59. García-Berro, E., & Lorén-Aguilar, P. 2017, in Dynamical Mergers, eds. A. W. Alsabti, & P. Murdin (Springer International Publishing AG), 1237 [Google Scholar]
  60. Garnavich, P. M., Bonanos, A. Z., Krisciunas, K., et al. 2004, ApJ, 613, 1120 [NASA ADS] [CrossRef] [Google Scholar]
  61. Goldhaber, G., Groom, D. E., Kim, A., et al. 2001, ApJ, 558, 359 [NASA ADS] [CrossRef] [Google Scholar]
  62. González Delgado, R. M., García-Benito, R., Pérez, E., et al. 2015, A&A, 581, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. González-Gaitán, S., Perrett, K., Sullivan, M., et al. 2011, ApJ, 727, 107 [NASA ADS] [CrossRef] [Google Scholar]
  64. Goobar, A. 2008, ApJ, 686, L103 [NASA ADS] [CrossRef] [Google Scholar]
  65. Hachinger, S., Mazzali, P. A., Tanaka, M., Hillebrandt, W., & Benetti, S. 2008, MNRAS, 389, 1087 [NASA ADS] [CrossRef] [Google Scholar]
  66. Hachinger, S., Mazzali, P. A., Taubenberger, S., Pakmor, R., & Hillebrandt, W. 2009, MNRAS, 399, 1238 [NASA ADS] [CrossRef] [Google Scholar]
  67. Hamuy, M., Phillips, M. M., Suntzeff, N. B., et al. 1996, AJ, 112, 2391 [NASA ADS] [CrossRef] [Google Scholar]
  68. Hamuy, M., Phillips, M. M., Suntzeff, N. B., et al. 2003, Nature, 424, 651 [NASA ADS] [CrossRef] [Google Scholar]
  69. Han, D.-H., Park, C., Choi, Y.-Y., & Park, M.-G. 2010, ApJ, 724, 502 [NASA ADS] [CrossRef] [Google Scholar]
  70. Henden, A. A., Levine, S., Terrell, D., & Welch, D. L. 2015, Am. Astron. Soc. Meet. Abstr., 225, 336.16 [Google Scholar]
  71. Hicken, M., Challis, P., Jha, S., et al. 2009, ApJ, 700, 331 [NASA ADS] [CrossRef] [Google Scholar]
  72. Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [NASA ADS] [CrossRef] [Google Scholar]
  73. Höflich, P. 1995, ApJ, 443, 89 [NASA ADS] [CrossRef] [Google Scholar]
  74. Höflich, P. 2003, in Stellar Atmosphere Modeling, eds. I. Hubeny, D. Mihalas, & K. Werner, ASP Conf. Ser., 288, 185 [NASA ADS] [Google Scholar]
  75. Höflich, P., & Stein, J. 2002, ApJ, 568, 779 [NASA ADS] [CrossRef] [Google Scholar]
  76. Höflich, P., Gerardy, C. L., Fesen, R. A., & Sakai, S. 2002, ApJ, 568, 791 [NASA ADS] [CrossRef] [Google Scholar]
  77. Höflich, P., Gerardy, C. L., Marion, H., & Quimby, R. 2006, New Astron. Rev., 50, 470 [NASA ADS] [CrossRef] [Google Scholar]
  78. Hoeflich, P., Dragulin, P., Mitchell, J., et al. 2013, Front. Phys., 8, 144 [NASA ADS] [CrossRef] [Google Scholar]
  79. Höflich, P., Hsiao, E. Y., Ashall, C., et al. 2017, ApJ, 846, 58 [NASA ADS] [CrossRef] [Google Scholar]
  80. Höflich, P., Ashall, C., Fisher, A., et al. 2018, ArXiv e-prints [arXiv:1810.05910] [Google Scholar]
  81. Howell, D. A. 2011, Nat. Commun., 2, 350 [NASA ADS] [CrossRef] [Google Scholar]
  82. Howell, D. A., Höflich, P., Wang, L., & Wheeler, J. C. 2001, ApJ, 556, 302 [NASA ADS] [CrossRef] [Google Scholar]
  83. Hoyle, F., & Fowler, W. A. 1960, ApJ, 132, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Hristov, B., Collins, D. C., Höflich, P., Weatherford, C. A., & Diamond, T. R. 2018, ApJ, 858, 13 [NASA ADS] [CrossRef] [Google Scholar]
  85. Hsiao, E. Y., Conley, A., Howell, D. A., et al. 2007, ApJ, 663, 1187 [NASA ADS] [CrossRef] [Google Scholar]
  86. Hsiao, E. Y., Marion, G. H., Phillips, M. M., et al. 2013, ApJ, 766, 72 [NASA ADS] [CrossRef] [Google Scholar]
  87. Hsiao, E. Y., Burns, C. R., Contreras, C., et al. 2015, A&A, 578, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Hsiao, E. Y., Phillips, M. M., Marion, G. H., et al. 2019, PASP, 131, 014002 [NASA ADS] [CrossRef] [Google Scholar]
  89. Jester, S., Schneider, D. P., Richards, G. T., et al. 2005, AJ, 130, 873 [NASA ADS] [CrossRef] [Google Scholar]
  90. Kasen, D., & Woosley, S. E. 2007, ApJ, 656, 661 [NASA ADS] [CrossRef] [Google Scholar]
  91. Kasliwal, M. M., Kulkarni, S. R., Gal-Yam, A., et al. 2012, ApJ, 755, 161 [NASA ADS] [CrossRef] [Google Scholar]
  92. Kelz, A., Verheijen, M. A. W., Roth, M. M., et al. 2006, PASP, 118, 129 [NASA ADS] [CrossRef] [Google Scholar]
  93. Khokhlov, A. M. 1989, MNRAS, 239, 785 [NASA ADS] [Google Scholar]
  94. Krisciunas, K., Suntzeff, N. B., Phillips, M. M., et al. 2004, AJ, 128, 3034 [NASA ADS] [CrossRef] [Google Scholar]
  95. Krisciunas, K., Contreras, C., Burns, C. R., et al. 2017, AJ, 154, 211 [NASA ADS] [CrossRef] [Google Scholar]
  96. Kromer, M., Sim, S. A., Fink, M., et al. 2010, ApJ, 719, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  97. Kromer, M., Pakmor, R., Taubenberger, S., et al. 2013, ApJ, 778, L18 [NASA ADS] [CrossRef] [Google Scholar]
  98. Langanke, K., & Martínez-Pinedo, G. 2003, Rev. Mod. Phys., 75, 819 [NASA ADS] [CrossRef] [Google Scholar]
  99. Leibundgut, B., Kirshner, R. P., Phillips, M. M., et al. 1993, AJ, 105, 301 [NASA ADS] [CrossRef] [Google Scholar]
  100. Lira, P., Suntzeff, N. B., Phillips, M. M., et al. 1998, AJ, 115, 234 [NASA ADS] [CrossRef] [Google Scholar]
  101. López Fernández, R., Cid Fernandes, R., González Delgado, R. M., et al. 2016, MNRAS, 458, 184 [NASA ADS] [CrossRef] [Google Scholar]
  102. Lucy, L. B. 1999, A&A, 345, 211 [NASA ADS] [Google Scholar]
  103. Lyman, J. D., Taddia, F., Stritzinger, M. D., et al. 2018, MNRAS, 473, 1359 [NASA ADS] [CrossRef] [Google Scholar]
  104. Mannucci, F., Della Valle, M., & Panagia, N. 2006, MNRAS, 370, 773 [NASA ADS] [CrossRef] [Google Scholar]
  105. Marion, G. H., Höflich, P., Vacca, W. D., & Wheeler, J. C. 2003, ApJ, 591, 316 [NASA ADS] [CrossRef] [Google Scholar]
  106. Marion, G. H., Höflich, P., Gerardy, C. L., et al. 2009, AJ, 138, 727 [NASA ADS] [CrossRef] [Google Scholar]
  107. Marino, S., González-Gaitán, S., Förster, F., et al. 2015, ApJ, 806, 134 [NASA ADS] [CrossRef] [Google Scholar]
  108. Martínez-Pinedo, G., Lam, Y. H., Langanke, K., Zegers, R. G. T., & Sullivan, C. 2014, Phys. Rev. C, 89, 045806 [NASA ADS] [CrossRef] [Google Scholar]
  109. Martini, P., Stoll, R., Derwent, M. A., et al. 2011, PASP, 123, 187 [NASA ADS] [CrossRef] [Google Scholar]
  110. Mazzali, P. A. 2000, A&A, 363, 705 [NASA ADS] [Google Scholar]
  111. Mazzali, P. A., & Lucy, L. B. 1993, A&A, 279, 447 [NASA ADS] [Google Scholar]
  112. Mazzali, P. A., Nomoto, K., Cappellaro, E., et al. 2001, ApJ, 547, 988 [NASA ADS] [CrossRef] [Google Scholar]
  113. Mazzali, P. A., Sauer, D. N., Pastorello, A., Benetti, S., & Hillebrandt, W. 2008, MNRAS, 386, 1897 [NASA ADS] [CrossRef] [Google Scholar]
  114. Mazzali, P. A., Sullivan, M., Hachinger, S., et al. 2014, MNRAS, 439, 1959 [NASA ADS] [CrossRef] [Google Scholar]
  115. Mazzali, P. A., Ashall, C., Pian, E., et al. 2018, MNRAS, 476, 2905 [NASA ADS] [CrossRef] [Google Scholar]
  116. McCully, C., Jha, S. W., Foley, R. J., et al. 2014a, Nature, 512, 54 [NASA ADS] [CrossRef] [Google Scholar]
  117. McCully, C., Jha, S. W., Foley, R. J., et al. 2014b, ApJ, 786, 134 [NASA ADS] [CrossRef] [Google Scholar]
  118. Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 [NASA ADS] [CrossRef] [Google Scholar]
  119. Moorwood, A., Cuby, J.-G., & Lidman, C. 1998, The Messenger, 91, 9 [NASA ADS] [Google Scholar]
  120. Nomoto, K. 1982, ApJ, 253, 798 [NASA ADS] [CrossRef] [Google Scholar]
  121. Nomoto, K., & Leung, S.-C. 2018, Space Sci. Rev., 214, 67 [NASA ADS] [CrossRef] [Google Scholar]
  122. Nomoto, K., Thielemann, F.-K., & Yokoi, K. 1984, ApJ, 286, 644 [Google Scholar]
  123. Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (University Science Books) [Google Scholar]
  124. Paczyński, B. 1972, Astrophys. Lett., 11, 53 [NASA ADS] [Google Scholar]
  125. Pan, Y.-C., Duarte, A. S., Foley, R. J., et al. 2016, ATel, 9705 [Google Scholar]
  126. Patat, F., Höflich, P., Baade, D., et al. 2012, A&A, 545, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  127. Penney, R., & Höflich, P. 2014, ApJ, 795, 84 [NASA ADS] [CrossRef] [Google Scholar]
  128. Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565 [NASA ADS] [CrossRef] [Google Scholar]
  129. Phillips, M. M. 1993, ApJ, 413, L105 [NASA ADS] [CrossRef] [Google Scholar]
  130. Phillips, M. M., Phillips, A. C., Heathcote, S. R., et al. 1987, PASP, 99, 592 [NASA ADS] [CrossRef] [Google Scholar]
  131. Phillips, M. M., Lira, P., Suntzeff, N. B., et al. 1999, AJ, 118, 1766 [NASA ADS] [CrossRef] [Google Scholar]
  132. Phillips, M. M., Simon, J. D., Morrell, N., et al. 2013, ApJ, 779, 38 [NASA ADS] [CrossRef] [Google Scholar]
  133. Piascik, A. S., Steele, I. A., Bates, S. D., et al. 2014, in Ground-based and Airborne Instrumentation for Astronomy V, Proc. SPIE, 9147, 91478H [Google Scholar]
  134. Piersanti, L., Tornambé, A., & Castellani, V. 2004, MNRAS, 353, 243 [NASA ADS] [CrossRef] [Google Scholar]
  135. Piersanti, L., Tornambé, A., & Yungelson, L. R. 2014, MNRAS, 445, 3239 [NASA ADS] [CrossRef] [Google Scholar]
  136. Piersanti, L., Bravo, E., Cristallo, S., et al. 2017, ApJ, 836, L9 [NASA ADS] [CrossRef] [Google Scholar]
  137. Pinto, P. A., & Eastman, R. G. 2000, ApJ, 530, 744 [NASA ADS] [CrossRef] [Google Scholar]
  138. Poznanski, D., Ganeshalingam, M., Silverman, J. M., & Filippenko, A. V. 2011, MNRAS, 415, L81 [NASA ADS] [CrossRef] [Google Scholar]
  139. Poznanski, D., Prochaska, J. X., & Bloom, J. S. 2012, MNRAS, 426, 1465 [NASA ADS] [CrossRef] [Google Scholar]
  140. Prentice, S. J., Ashall, C., Mazzali, P. A., et al. 2018, MNRAS, 478, 4162 [NASA ADS] [CrossRef] [Google Scholar]
  141. Richmond, M. W., & Smith, H. A. 2012, J. Am. Assoc. Var. Star Obs., 40, 872 [Google Scholar]
  142. Riess, A. G., Press, W. H., & Kirshner, R. P. 1996, ApJ, 473, 88 [NASA ADS] [CrossRef] [Google Scholar]
  143. Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009 [NASA ADS] [CrossRef] [Google Scholar]
  144. Röpke, F. K. 2005, A&A, 432, 969 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  145. Roth, M. M., Kelz, A., Fechner, T., et al. 2005, PASP, 117, 620 [NASA ADS] [CrossRef] [Google Scholar]
  146. Sánchez, S. F., Rosales-Ortega, F. F., Marino, R. A., et al. 2012, A&A, 546, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  147. Sánchez, S. F., García-Benito, R., Zibetti, S., et al. 2016, A&A, 594, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  148. Scalzo, R., Aldering, G., Antilogus, P., et al. 2014, MNRAS, 440, 1498 [NASA ADS] [CrossRef] [Google Scholar]
  149. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  150. Schmidt, G. D., Weymann, R. J., & Foltz, C. B. 1989, PASP, 101, 713 [NASA ADS] [CrossRef] [Google Scholar]
  151. Sell, P. H., Arur, K., Maccarone, T. J., et al. 2018, MNRAS, 475, L111 [NASA ADS] [CrossRef] [Google Scholar]
  152. Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48 [NASA ADS] [CrossRef] [Google Scholar]
  153. Shen, K. J., & Moore, K. 2014, ApJ, 797, 46 [NASA ADS] [CrossRef] [Google Scholar]
  154. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  155. Smartt, S. J., Valenti, S., Fraser, M., et al. 2015, A&A, 579, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  156. Stalder, B., Tonry, J., Smartt, S. J., et al. 2017, ApJ, 850, 149 [NASA ADS] [CrossRef] [Google Scholar]
  157. Steele, I. A., Smith, R. J., Rees, P. C., et al. 2004, in Ground-based Telescopes, ed. J. M. Oschmann, Jr., Proc. SPIE, 5489, 679 [CrossRef] [Google Scholar]
  158. Stehle, M., Mazzali, P. A., Benetti, S., & Hillebrandt, W. 2005, MNRAS, 360, 1231 [NASA ADS] [CrossRef] [Google Scholar]
  159. Stein, J., & Wheeler, J. C. 2006, ApJ, 643, 1190 [NASA ADS] [CrossRef] [Google Scholar]
  160. Stritzinger, M., Mazzali, P. A., Sollerman, J., & Benetti, S. 2006, A&A, 460, 793 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  161. Stritzinger, M., Burns, C. R., Phillips, M. M., et al. 2010, AJ, 140, 2036 [NASA ADS] [CrossRef] [Google Scholar]
  162. Stritzinger, M. D., Hsiao, E., Valenti, S., et al. 2014, A&A, 561, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  163. Stritzinger, M. D., Valenti, S., Hoeflich, P., et al. 2015, A&A, 573, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  164. Sullivan, M., Kasliwal, M. M., Nugent, P. E., et al. 2011, ApJ, 732, 118 [NASA ADS] [CrossRef] [Google Scholar]
  165. Tanaka, M., Mazzali, P. A., Stanishev, V., et al. 2011, MNRAS, 410, 1725 [NASA ADS] [CrossRef] [Google Scholar]
  166. Taubenberger, S. 2017, Handbook of Supernovae (Springer International Publishing AG), 317 [CrossRef] [Google Scholar]
  167. Taubenberger, S., Hachinger, S., Pignata, G., et al. 2008, MNRAS, 385, 75 [NASA ADS] [CrossRef] [Google Scholar]
  168. Taubenberger, S., Kromer, M., Pakmor, R., et al. 2013, ApJ, 775, L43 [NASA ADS] [CrossRef] [Google Scholar]
  169. Thielemann, F.-K., Argast, D., Brachwitz, F., et al. 2002, Ap&SS, 281, 25 [NASA ADS] [CrossRef] [Google Scholar]
  170. Thomas, R. C., Nugent, P. E., & Meza, J. C. 2011, PASP, 123, 237 [NASA ADS] [CrossRef] [Google Scholar]
  171. Thompson, T. A. 2011, ApJ, 741, 82 [NASA ADS] [CrossRef] [Google Scholar]
  172. Tonry, J. L., Denneau, L., Heinze, A. N., et al. 2018, PASP, 130, 064505 [NASA ADS] [CrossRef] [Google Scholar]
  173. Tripp, R. 1998, A&A, 331, 815 [NASA ADS] [Google Scholar]
  174. Tucker, M. A., Shappee, B. J., Vallely, P. J., et al. 2019, MNRAS, submitted [arXiv:1903.05115] [Google Scholar]
  175. Turatto, M., Benetti, S., Cappellaro, E., et al. 1996, MNRAS, 283, 1 [NASA ADS] [CrossRef] [Google Scholar]
  176. Verheijen, M. A. W., Bershady, M. A., Andersen, D. R., et al. 2004, Astron. Nachr., 325, 151 [NASA ADS] [CrossRef] [Google Scholar]
  177. Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  178. Villanueva, S., Eastman, J. D., Gaudi, B. S., et al. 2016, in Ground-based and Airborne Telescopes VI, Proc. SPIE, 9906, 99062L [Google Scholar]
  179. Villanueva, Jr., S., Gaudi, B. S., Pogge, R. W., et al. 2018, PASP, 130, 015001 [NASA ADS] [CrossRef] [Google Scholar]
  180. Vorontsov-Vel’Yaminov, B. A., & Arkhipova, V. P. 1962, Morphological Catalogue of Galaxies. Part 1 [Google Scholar]
  181. Wang, X., Li, W., Filippenko, A. V., et al. 2008, ApJ, 675, 626 [NASA ADS] [CrossRef] [Google Scholar]
  182. Wang, X., Filippenko, A. V., Ganeshalingam, M., et al. 2009, ApJ, 699, L139 [NASA ADS] [CrossRef] [Google Scholar]
  183. Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 [NASA ADS] [CrossRef] [Google Scholar]
  184. Wygoda, N., Katz, B., & Elbaz, Y. 2019, MNRAS, 484, 3951 [NASA ADS] [CrossRef] [Google Scholar]
  185. Yoon, S.-C., & Langer, N. 2004, A&A, 419, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  186. Yoon, S.-C., Podsiadlowski, P., & Rosswog, S. 2007, MNRAS, 380, 933 [NASA ADS] [CrossRef] [Google Scholar]
  187. Zhang, J.-J., Wang, X.-F., Bai, J.-M., et al. 2014, AJ, 148, 1 [NASA ADS] [CrossRef] [Google Scholar]
  188. Zingale, M., & Dursi, L. J. 2007, ApJ, 656, 333 [NASA ADS] [CrossRef] [Google Scholar]
  189. Zingale, M., Nonaka, A., Almgren, A. S., et al. 2011, ApJ, 740, 8 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Light curve excess tests

We performed a number of tests to try to explain the late light curve excess observed in SN 2016hnk. In the main text we provided the most plausible explanation, multi-scattering by dust clouds at 1.0 ± 0.5 pc from the SN, while we give more details of other tests in this appendix.

A.1. Reddening SN 1991bg spectrum and Na I D absorption

From Figs. 8 and 14, it is obvious that SN 2016hnk light and color curves show some peculiarities. The first question to address is whether SN 2016hnk is reddened or is its red color intrinsic? We took SN 1991bg spectrum at maximum light and, after correcting for Milky Way extinction and shifting to the rest frame, we artificially applied the same reddening we found for SN 2016hnk (E(B − V) = 0.4−0.5 mag with an RV = 2.1). Figure A.1 shows the SN 1991bg spectrum at maximum (in black), together with the corresponding SN 2016hnk spectrum also at maximum (in blue; also corrected for Milky Way reddening and shifted to the rest frame), and three SN 1991bg reddened spectra with E(B − V) = 0.40, 0.45, and 0.50 mag (in red, purple, and orange, respectively). Besides the differences already pointed out in Sect. 2.2 (broader NIR calcium absorption, and features in the range 6400−7300 Å), we can see an almost perfect match between SN 2016hnk and all reddened SN 1991bg spectra redward of 5200 Å. However, the blue part of the spectrum is where most of the discrepancy comes from. Both the aggregate Fe II, Si II, and Mg II features at ∼5000 Å, and the TI II at ∼4200 Å, are deeper in SN 2016hnk than in the reddened SN 1991bg spectra, which makes SN 2016hnk color curves even redder than expected by simple dust extinction. Thus, we conclude that there should be an intrinsic effect, in addition to reddening responsible for the large red colors of SN 2016hnk.

thumbnail Fig. A.1.

SN 1991bg (black) and SN 2016hnk (blue) spectra at maximum, together with three artificially reddened SN 1991bg spectra using an RV = 2.1 and E(B − V) = 0.40 (red), 0.45 (purple), and 0.50 (orange) mag. All spectra were first corrected for Milky Way extinction and rest-framed.

Another peculiarity of SN 2016hnk is the lack of Na I D absorption at any epoch, even when we obtain large E(B − V) (> 0.4 mag) with different methods. Although the relation between extinction and Na I D pEW has been put in question when measured in low-resolution spectra (Poznanski et al. 2011), it is usually accepted that the presence of Na I D pEW absorption is an indication of extinction by ISM (Phillips et al. 2013). The relations from Barbon et al. (1990) and Poznanski et al. (2012) predict both Na I D pEW of at least 1.3 Å, from the E(B − V)∼0.4 mag measured in SN 2016hnk. We artificially input Na I D absorption with pEW of 1.3 Å in several SN 2016hnk spectra at different epochs, as shown in Fig. A.2, taking into account the spectral resolution of each spectrum. Even for the latter spectrum at +59.7 d that has the lowest S/N, we can significantly discard any indication of Na I D absorption. Therefore, in the case of SN 2016hnk, the large extinction does not come together with Na I D absorption.

thumbnail Fig. A.2.

Zoom into SN 2016hnk spectra around the Na I D doublet (rest-frame 5890 and 5896 Å; in vertical black dotted lines) absorption, at seven epochs from +2.5 to +59.7 days after maximum. Thick lines represent Na I D absorption with pEW = 1.3 Å at the spectral resolution of each spectrum. Red vertical lines indicate the position where the Na I D doublet from the Milky Way should be found. In all cases, we do not detect the Na I D lines in SN 2016hnk spectra.

A.2. Simple light echo model

Another possible explanation for the light excess in the SN 2016hnk light curve after +15 days past maximum would be an intervening dust cloud close to the line of sight that is producing a LE (Marino et al. 2015). To test this possibility, we constructed a simple model accounting for different contributions of the LE to the total observed light in the form,

(A.1)

where is the actual emitted flux from the SN in the X band; is the flux we observed in the same band; is the re-emitted flux from the cloud; α is a multiplicative factor that corresponds to the contribution of the LE, with respect to the actual emission from the SN (and which we assumed to be band independent), in the sense that α = 0 means no contribution at all and α = 1 corresponds to the extreme case of perfect reflection, resulting in fSN = fLE; and S is a step function that has a value of one for tdelay <  t <  tend and zero at other epochs.

To estimate we assumed that early epochs are free of LE contribution and performed a SNOOPY fit (fixing tmax and RV to the same values found on Sect. 2.1) but only up to +15 days past maximum. Then, the peak flux of the best-fit model was taken as for each band independently, and used to generate a series of LE-free light curves with different combinations of α parameters from 0.02 to 0.14, in steps of 0.02; tdelay from 5 to 30 days, in steps of 5 days; and tend from 40 days on, in steps of 5 days. See Fig. A.3 as an example. Those sets of light curves were fitted with SNOOPY keeping the assumptions on the tmax and RV, but extending the fit to the whole time extent. The best fit, shown in Fig. A.4, was for a combination of parameters (α, tdelay, tend) = (0.10, 15, 80). We note that the best tend value is right after the last photometric point, meaning that the LE is contributing flux at least until the epoch when we stopped our follow up campaign. The model also provided the best estimates of sBV = 0.380 ± 0.032 and E(B − V) = 0.458 ± 0.064 mag. sBV is a bit lower than the reported in the main text and E(B − V) is consistent with the value found using CDMagic.

thumbnail Fig. A.3.

Effect of subtracting different fractional contribution of the LE to the SN 2016hnk B band observed light curve.

thumbnail Fig. A.4.

Multiband SN 2016hnk SNOOPY fits using a 91bg template and with a 10% contribution of a LE starting at +15 days until the end of the light curve. Fits in all band are reasonably better than those without a LE shown in Fig. 7.

Although light curve fits improve significantly when this simple model of the LE is removed from observations, if this excess is due to a LE, we should further see a mixture of spectral features in the spectra. At latter epochs when the LE is affecting the observations (> +15 d in our model), we should see spectral features from earlier epochs that are being reflected and that are contaminating the light directly received by the observer from the SN. This contamination is strongly dependent on the configuration of the ISM around the SN but, again, in the simplest scenario (single cloud at a certain distance from the SN) this would correspond to a mixture of the spectrum at the actual epoch and a spectrum at an early epoch eq27.

To test if these features are present at late epochs we took the spectrum of a normal 1991bg-like SNIa, SN 2005ke at +23.6 days and added up incremental contributions of SN 2016hnk spectrum at maximum light, which represents the contamination from the LE. The resulting mixed spectrum is compared to the observed spectrum of SN 2016hk at +20.3 days to check if the mixed features from those two spectra reproduce those seen at late epochs in SN 2016hnk. Starting from no contribution from the LE, and increasing in steps of 1% up to 50% contribution to the total light, we find the most similar spectrum for a 39% LE contribution (based on χ2 minimization). In this best case, shown in Fig. A.5, we do not find a convincing resemblance between the two spectra (observed, in black, and simulated, in purple) that confirms the hypothesis of a LE causing those features. For this reason, we do not consider the LE explanation for the late-time light excess and favor the multi-scattering approach detailed in the main text.

thumbnail Fig. A.5.

Simulation of a LE contaminated spectrum at three week past maximum, based on a fiducial subluminous spectrum at +3 weeks past maximum from SN 2005ke (contributing 61% to the total), and adding 39% of the SN 2016hnk spectrum at maximum light, which represents the contamination from the LE. Comparing this LE-affected spectrum (in purple) to the observed SN 2016hnk (in black) we do not see any feature that definitely confirms that the latter is affected by any feature reflected from earlier epochs.

Appendix B: Spectral modeling technique

The code is a 1D Monte Carlo radiation transport code that follows the propagation of photon packets through a SN atmosphere. The photon packets can undergo Thomson scattering and line absorption. Photon packets have two fates: they either re-enter the photosphere in a process known as backscattering or escape the SN atmosphere after possibly undergoing scattering or line absorption. A modified nebular approximation is used to account for non-LTE effects. The code makes use of the Schuster-Schwarzschild approximation, which assumes that radiative energy is emitted from an inner blackbody. This has the advantage that the radiation transport below the photosphere is not required to be known, but can cause excess flux in the red (> 6500 Å) past maximum light. We note that the values we obtain for vph are, in most cases, lower than those measured from the individual ions in the spectra in Sect. 5 because the lines are produced in the ejecta above the photosphere. For further details on the abundance stratification code and method used in this work see Mazzali & Lucy (1993), Lucy (1999), Mazzali (2000), and Stehle et al. (2005).

Appendix C: Photometry, local sequence, and spectral log

All photometric and spectroscopic observations are publicly available at https://github.com/lgalbany/SN2016hnk. The spectra are also available at the Weizmann Interactive Supernova Data Repository (WISeREP; https://wiserep.weizmann.ac.il).

Table C.1.

ASASSN V-band photometry and upper limits.

Table C.2.

ATLAS c- and o-band photometry and upper limits.

Table C.3.

Photometry for SN 2016hnk. Epochs are given in the rest frame.

Table C.4.

Near-infrared photometry for SN 2016hnk. Epochs from peak are given in the rest frame.

Table C.5.

Journal of spectroscopic observations for SN 2016hnk.

Table C.6.

Measured velocities (in units of 103 km s−1). Epochs are measured from the B-band maximum brightness.

Table C.7.

Measured pseudo-equivalent widths (in Å). Epochs are measured from the B-band maximum brightness.

All Tables

Table 1.

Properties of SN 2016hnk and its host galaxy.

Table 2.

Scattering parameters (albedo, average of the cosine of the scattering angle, and absorption cross section divided by dust mass) for optical photons from Weingartner & Draine (2001) corresponding to interstellar extinction in the LMC.

Table 3.

Parameters measured in this work and used in the diagrams presented in Fig. 17.

Table 4.

Input parameters for the early-time spectral models.

Table 5.

Asymptotic isotopic yields at large times for the delay-detonation model for SN 2016hnk discussed in Sect. 6.2.

Table C.1.

ASASSN V-band photometry and upper limits.

Table C.2.

ATLAS c- and o-band photometry and upper limits.

Table C.3.

Photometry for SN 2016hnk. Epochs are given in the rest frame.

Table C.4.

Near-infrared photometry for SN 2016hnk. Epochs from peak are given in the rest frame.

Table C.5.

Journal of spectroscopic observations for SN 2016hnk.

Table C.6.

Measured velocities (in units of 103 km s−1). Epochs are measured from the B-band maximum brightness.

Table C.7.

Measured pseudo-equivalent widths (in Å). Epochs are measured from the B-band maximum brightness.

All Figures

thumbnail Fig. 1.

Grayscale B-band image obtained with the Las Cumbres 1 m telescope at Cerro Tololo Inter-American Observatory on 2016 November 3. North is up; east is left. SN 2016hnk is denoted with two blue lines at the bottom inner panel, next to its host galaxy MCG -01-06-070. Red circles represent all 35 standard stars used for photometric calibration in the optical and NIR.

In the text
thumbnail Fig. 2.

g-band imaging of SN 2016hnk and its host galaxy, before (left) and after (middle) subtracting the host galaxy template using hotpants integration in SNOOPY. This example shows how the subtracted image is clean enough to perform PSF fitting, as shown in the right panel where the SN PSF has been removed. A small residual of the galaxy core is shown on the middle and right panels (white blob), however it does not affect SN PSF fitting.

In the text
thumbnail Fig. 3.

Continuum-subtracted Hα emission (left), HαEW (middle), and Hα velocity (right) maps of MCG -01-06-070. HII regions extracted with HIIEXPLORER are indicated in brown polygons in the left panel. The red star and the black/white × represent the locations of SN 2016hnk and MCG -01-06-070 nucleus, respectively.

In the text
thumbnail Fig. 4.

Top: MCG -01-06-070 integrated spectrum (blue) and best SSP fit from STARLIGHT (red). Vertical strips (in beige) represent masked regions in the fit due to known emission lines (such as Balmer lines, oxygen, and WR bumps) and regions with skylines (in brown). Inner panel: star formation history of the spectra (in purple), where brown ticks on top represent the ages of the models used in the fit. The S/N of the spectrum, optical extinction, average stellar age, and average metallicity are reported in the top left corners. Bottom rows: distributions of the Hα equivalent width, SFR, oxygen abundance, and average stellar age of all 180 SN Ia host galaxies in PISCO (Galbany et al. 2018). The large colored dots represent the position of MCG -01-06-070, and the three small dots represent other 1991bg-like SNIa hosts in PISCO.

In the text
thumbnail Fig. 5.

Top: similar to Fig. 4, the spectrum of the SN 2016hnk parent HII region (blue) and best SSP fit from STARLIGHT (red) are shown. Right panels: stellar subtracted spectra around the Hα and Hβ wavelengths. Bottom rows: distributions of the gas-based parameters of all 56 HII regions in MCG -01-06-070, weighted by the contribution of the region to the total SFR of the galaxy (dark lines) or unweighted (light lines). In this case, colored horizontal bands represent the band where the SN 2016hnk parent HII region lies within.

In the text
thumbnail Fig. 6.

Optical and NIR light curves. A gray solid vertical line represents the epoch of peak brightness in the B band obtained from a SNOOPY fit to the cyan and orange ATLAS light curves, which are shown in orange and cyan solid lines. The ticks at the bottom represent epochs with optical spectra available. The dotted vertical lines delineate the last non-detection (yellow) and first detection (orange) epochs. The open symbols with downward arrows correspond to ASAS-SN (circles), ATLAS (crosses/diamonds), and u-band (triangles) non-detections.

In the text
thumbnail Fig. 7.

Multiband light curves of SN 2016hnk (filled circles) compared with a 1991bg-like template (lines). Best-fit SNOOPY parameters are RV = 2.1 and a host-galaxy color excess of E(B − V)host = 0.55 mag. We note that beyond +20 d the template fit does not succeed in capturing the SN brightness.

In the text
thumbnail Fig. 8.

B-band light curves of SN 2016hnk and a variety of SNe Ia of different subtypes: normal SN Ia (2011fe), 2002es-like (2002es, 2010lp, iPTF14atg), 2002cx-like (2008ae), transitional (1986G), and 1991bg-like (1991bg, 1999by, 2005bl, 2005ke, 2006mr). All data were rest-framed and scaled to match the SN 2016hnk peak. Vertical lines represent epochs of the peak and +15 days, and horizontal lines indicate the crossing point at 15 days of a SN with Δm15(B) = 1.3, 1.5, and 1.7 mag, respectively. SN 2016hnk has one of the lowest Δm15(B) but is the object with the brightest tail in both bands.

In the text
thumbnail Fig. 9.

Representative panel for the CDMagic method, where a simultaneous fit in the B- and V-band light curves, the B vs. (B − V) and V vs. (B − V) diagrams, and the (B − V) evolution, is performed. More details are given in, for example, Fig. 13 in Höflich et al. (2017). We show the absolute V-band vs. (B − V) diagram populated with eight normal SNe Ia from the CSP-I (Krisciunas et al. 2017; SNe 2004ef, 2005eo, 2005el, 2005iq, 2005ki, 2005M, 2006bh, 2006D; in differnet symbols), the SN Ia 1991bg-like 2005ke (empty triangles), and the uncorrected (blue) and corrected (red) SN 2016hnk data points. Models from Höflich et al. (2017) are shown in green for reference; the two arrows are the reddening vectors representing the direction to which the data points are shifted when affected by RV = 3.2 and RV = 1.0. SN 2016hnk corrected data points lie on top of the other subluminous SN Ia 2005ke, and from the shift we estimate an RV = 2.1 ± 0.6.

In the text
thumbnail Fig. 10.

Luminosity width relation as a function of Δm15(B) (top) and sBV (bottom) for a selection of CSP-I SNe Ia from Krisciunas et al. (2017) with the addition of SN 2016hnk (black circle). We note that SN 2016hnk sits off the Δm15(B)−MB relationship because of the degeneracy in Δm15(B) for fast declining SNe Ia.

In the text
thumbnail Fig. 11.

SN 2016hnk color curves (top panel), the B − V color curve of SN 2016hnk, the 1991bg-like comparison sample, and SNe 2002es and 2011fe for reference. The shaded region projects the SN 2016hnk color curve down by 0.5 mag, on top of the bulk of other 1991bg-like SNe Ia. In pink we show the linear fit the to B − V curve from 30 to 70 days past B-band maximum.

In the text
thumbnail Fig. 12.

BVri observed light curves (in circles), compared with the multi-scattering dust model of Bulla et al. (2018b) using 1991bg-like SN 2005ke as a reference (black dotted line). Simulated light curves affected by nearby dust with E(B − V) = 0.5 mag at a distance of 1.0 ± 0.5 pc from the SN (blue filled strip) provide very reasonable fits to the observations, while a model in which dust is located at larger distance (100 pc; blue dashed line) just dims the light curves by keeping the same shape.

In the text
thumbnail Fig. 13.

Optical spectral sequence for SN 2016hnk including 22 observations listed in Table C.5 (shifted by 1.5e−16 erg s−1 cm−2 Å−1). In gray we present the color-corrected, Milky Way reddening-corrected, and rest-framed spectra, and in black a smoothed version using a 10 pixel box. Rest-frame epochs for each spectrum are attached in blue.

In the text
thumbnail Fig. 14.

Spectra of SN 2016hnk at three different epochs (around maximum, after a week, and around three weeks past maximum light) compared to other subluminous 1991bg-like SNe Ia, PTF09dav (Ca-rich), and SN 2002es (2002es-like SN Ia). The main typical features of 1991bg-like SNe Ia are identified in the top panel. Vertical lines in the top panel are located at the minima of a few features and serve as a comparison of the velocities with respect other objects. All spectra have been Milky Way extinction-corrected and shifted to the rest frame.

In the text
thumbnail Fig. 15.

Evolution of pseudo-equivalent widths pEW and velocities as a function of time for the most prominent features highlighted in the top panel of Fig. 14. Shaded regions represent the 1σ dispersion of COOL (brown) and BROAD LINE (orange) evolution in the CSP-I sample from Folatelli et al. (2013).

In the text
thumbnail Fig. 16.

Branch et al. (2006) and Wang et al. (2009) diagrams. The left panel is populated with measurements from the CfA sample by Blondin et al. (2012) (open symbols) and from the CSP-I sample by Folatelli et al. (2013) (filled symbols). Right panel: only the measurements from Folatelli et al. (2013) are included. In both panels we added our own measurements of SN 2016hnk and the sample of subluminous 1991bg-like SNe Ia. SN 2002es was added only in the left panel because with a velocity of ∼5800 km s−1 it falls outside the range shown in the right panel.

In the text
thumbnail Fig. 17.

Benetti et al. (2005) diagrams of R(Si II) and Si II velocity gradient as a function of Δm15 (left column). We replicated these diagrams using Δm15, s (center) and sBV (right column) instead as a SN light curve width parameter. In all panels, open symbols represent measurements from the original SN sample in Benetti et al. (2005), while filled symbols are measurements from the CSP-I (Folatelli et al. 2013).

In the text
thumbnail Fig. 18.

Near-infrared spectra of SN 2016hnk at around 6.1 and 8.5 (±3.1) days and, as a comparison, spectra of SN 2011fe (normal SN Ia), iPTF13ebh (transitional SN Ia), and SN 1999by (SN Ia 1991bg-like). We also include delayed-detonation models for SN 1999by some 24−27 days after the explosion, which corresponds to about 8−11 days past-maximum light in the B band (see Sect. 6.2 for more details). The SN 2016hnk spectra are dominated in the J band by blends of singly-ionized iron group elements, Ca II, and Si II, while the H band appears very flat with the H-band break indicated in red. The two deep features at around 10 000 and 11 500 Å are only seen shallower in SN 1999by spectra.

In the text
thumbnail Fig. 19.

Top: two-dimensional frame of the combined X-shooter nebular phase (+360 d) spectrum of SN 2016hnk after sky subtraction around the region where we detected the only feature. The trace can be clearly seen around 730 nm. The other feature around 750 nm is due to a bad telluric subtraction, since it is present along the spatial direction (vertical). Middle: combined 1D X-shooter optical spectrum showing the only two features detected in the whole UV+Opt+NIR wavelength range. Bottom: zoom-in to the feature we identified as being due to calcium emission. The blue and orange shaded regions represent the shift of ∼700 km s−1 between the peak of the emission and the rest-frame wavelength.

In the text
thumbnail Fig. 20.

Comparison of SN 2016hnk nebular spectra (in black) to other 1991bg-like, 2002es-like, 2002cx-like, and Ca-rich objects. While all these objects show Ca emission at late times, only SN 2002cx-like (and possibly 2002es-like; e.g., SN 2010lp) objects show the narrow emission lines similar to SN 2016hnk. The yellow shaded region is zoomed in on the right panel, where the vertical dotted lines represent the rest-frame wavelength of the two [Ca II] emission lines. We note that while SN 2005hk lines have negative velocities, SN 2008A lines have positive velocities similar to SN 2016hnk.

In the text
thumbnail Fig. 21.

Time series of data (black) and spectral models (blue) for SN 2016hnk during the photospheric phase. The data were corrected for foreground and host galaxy extinction. The observer-frame time from B-band maximum and the rest-frame time from explosion is provided for each model. Line identifications from the model are shown in the top plot.

In the text
thumbnail Fig. 22.

Abundance distribution as a function of velocity for SN 2016hnk. The vertical dashed lines denote the photosphere velocity of the +2.5 and +7.5 d models. We note that because of the small time step between our first and last model our results are only sensitive between 4500 and 10 000 km s−1.

In the text
thumbnail Fig. 23.

Comparison of SN 2016hnk (black) +2.5 d past maximum light and the non-LTE model of SN 1999by at +7 days past maximum light (blue) from Höflich et al. (2002). The telluric features are denoted in the observed spectra.

In the text
thumbnail Fig. 24.

Structure during the phase of homologous expansion as a function of the mass coordinate in WD-masses for our subluminous model with high central density. We give the normalized density as log(ρ/ρc) (left scale), velocity (right scale), and the mean optical depth log(τγ) (left scale) in γ-rays at day 350.

In the text
thumbnail Fig. 25.

Same as Fig. 24 but showing the overall abundance structure at large times as a function of mass coordinate (top) and expansion velocity (bottom). It is very similar to the subluminous model for SN 1999by and consistent with the chemical structure based on our abundance tomography given in Table 3. In addition, we show the initial 56Ni → 56Co → 56Fe at the time of the explosion (purple dotted line).

In the text
thumbnail Fig. 26.

Distribution of most abundant isotopes in log(Xi) and Ye as a function of velocity for the central region are dominated by neutron-rich isotopes at several minutes after the explosion. Subsequently, radioactive isotopes decay 55Cr (3.5 m) → 55Mn, 55Fe(2.7 y) → 55Mn, 59Fe (44.5 d) → 59Co, 56Ni (6 d) → 56Co (77 d) → 56Fe, 57Co (271 d) → 57Fe.

In the text
thumbnail Fig. 27.

Comparison of the nebular spectrum of SN 2016hnk with our high-density delayed-detonation model in the rest frame assuming a redshift of z = 0.01610. The flux is in 10−18 erg s−1 cm−2 on the same scale as Fig. 19. The telluric region in the original spectrum is plotted in gray, around 7500 Å. The theoretical spectrum is dominated by a narrow, forbidden [Ca II] doublet at 7291 and 7324 Å above a quasi-continuum at about the 3% level formed by many transitions of mostly heavier elements including Fe, Co, Ni, and Cr.

In the text
thumbnail Fig. A.1.

SN 1991bg (black) and SN 2016hnk (blue) spectra at maximum, together with three artificially reddened SN 1991bg spectra using an RV = 2.1 and E(B − V) = 0.40 (red), 0.45 (purple), and 0.50 (orange) mag. All spectra were first corrected for Milky Way extinction and rest-framed.

In the text
thumbnail Fig. A.2.

Zoom into SN 2016hnk spectra around the Na I D doublet (rest-frame 5890 and 5896 Å; in vertical black dotted lines) absorption, at seven epochs from +2.5 to +59.7 days after maximum. Thick lines represent Na I D absorption with pEW = 1.3 Å at the spectral resolution of each spectrum. Red vertical lines indicate the position where the Na I D doublet from the Milky Way should be found. In all cases, we do not detect the Na I D lines in SN 2016hnk spectra.

In the text
thumbnail Fig. A.3.

Effect of subtracting different fractional contribution of the LE to the SN 2016hnk B band observed light curve.

In the text
thumbnail Fig. A.4.

Multiband SN 2016hnk SNOOPY fits using a 91bg template and with a 10% contribution of a LE starting at +15 days until the end of the light curve. Fits in all band are reasonably better than those without a LE shown in Fig. 7.

In the text
thumbnail Fig. A.5.

Simulation of a LE contaminated spectrum at three week past maximum, based on a fiducial subluminous spectrum at +3 weeks past maximum from SN 2005ke (contributing 61% to the total), and adding 39% of the SN 2016hnk spectrum at maximum light, which represents the contamination from the LE. Comparing this LE-affected spectrum (in purple) to the observed SN 2016hnk (in black) we do not see any feature that definitely confirms that the latter is affected by any feature reflected from earlier epochs.

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.