The Extraplanar Type II Supernova ASASSN-14jb in the Nearby Edge-on Galaxy ESO 467-G051

The Extraplanar Type II Supernova ASASSN-14jb in the Nearby Edge-on Galaxy ESO 467-G051

Nicolás Meza Instituto de Astrofísica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 782-0436 Macul, Santiago, Chile    J. L. Prieto Núcleo de Astronomía de la Facultad de Ingeniería y Ciencias, Universidad Diego Portales, Av. Ejército 441, Santiago, Chile Millennium Institute of Astrophysics, Santiago, Chile    A. Clocchiatti Instituto de Astrofísica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 782-0436 Macul, Santiago, Chile Millennium Institute of Astrophysics, Santiago, Chile    L. Galbany PITT PACC, Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260, USA    J.P. Anderson European Southern Observatory, Alonso de Córdova 3107, Casilla 19, Santiago, Chile    E. Falco Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA    C. S. Kochanek Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University, 191 W. Woodruff Avenue, Columbus, OH 43210, USA    H. Kuncarayakti Tuorla Observatory, Department of Physics and Astronomy, University of Turku, Väisäläntie 20, FI-21500 Piikkiö, Finland    J. Brimacombe Coral Towers Observatory, Cairns, Queensland 4870, Australia    T. W.-S. Holoien The Observatories of the Carnegie Institution for Science, 813 Santa Barbara Street, Pasadena, CA 91101, USA    B. J. Shappee Institute for Astronomy, University of Hawai’i, 2680 Woodlawn Drive, Honolulu, HI 96822, USA    K. Z. Stanek Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University, 191 W. Woodruff Avenue, Columbus, OH 43210, USA    T. A. Thompson Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University, 191 W. Woodruff Avenue, Columbus, OH 43210, USA
Received / Accepted
Key Words.:
supernovae: individual: ASASSN-14jb – supernovae: general


Aims:We present optical photometry and spectroscopy of the Type II supernova ASASSN-14jb, together with VLT MUSE IFU observations of its host galaxy and a nebular-phase spectrum.

Methods:This supernova, in the nearby galaxy ESO 467-G051 (), was discovered and followed-up by the All Sky Automated Survey for SuperNovae (ASAS-SN). We obtained well-sampled LCOGTN and optical and near-UV/optical light curves and several optical spectra in the early photospheric phases. ASASSN-14jb exploded  kpc above the star-forming disk of ESO 467-G051, an edge-on disk galaxy. The large projected distance from the disk and non-detection of any H II region in a 1.4 kpc radius in projection are in conflict with the standard environment of core-collapse supernova progenitors and suggests the possible scenario that the progenitor received a kick in a binary interaction.

Results:We present analysis of the optical light curves and spectra, from which we derive a distance of Mpc using state of the art empirical methods for Type II SNe, physical properties of the SN explosion (Ni mass, explosion energy, and ejected mass) and properties of the progenitor, namely the progenitor radius, mass and metallicity. Our analysis yields a Ni mass of M, an explosion energy of ergs and an ejected mass of M. We also constrain the progenitor radius to be R which seems to be consistent with the sub-Solar metallicity of Z derived from the supernova Fe II line. The nebular spectrum constrains strongly the progenitor mass to be in the range 10-12 M. From Spitzer data archive we detect ASASSN-14jb days past explosion and we derive a total dust mass of M from the 3.6 m and 4.5 m photometry. Using the FUV, NUV, ,, 3.6 m, and 4.5 m total magnitudes for the host galaxy, we fit stellar population synthesis models which gives an estimate of M , an age of 3.2 Gyr, and a SFR M/yr. We also discuss the low oxygen abundance of the host galaxy derived from the MUSE data, having an average of using the O3N2 diagnostic, with strong line methods and compare it with the supernova spectra, which is also consistent with a sub-Solar metallicity progenitor. Following recent observations of extraplanar H II regions in nearby edge-on galaxies, we derive the metallicity offset from the disk, being positive (but consistent with zero at 2 ), suggesting enrichment from disk outflows. We finally discuss the possible scenarios for the unusual environment for ASASSN-14jb and conclude that either the in-situ star formation or runaway scenario would imply a low mass progenitor, agreeing with our estimate from the supernova nebular spectrum. Regardless of the true origin of ASASSN-14jb we show in this work that the detailed study of the environment can roughly agree with the stronger constrains of the transient observations.


1 Introduction

Originally classified based on the absence (Type I) or presence (Type II) of hydrogen lines in their optical spectra (Minkowski, 1941), supernovae (SNe) represent the explosive ending of a star. Decades of research have added a considerable degree of complexity to the simple scheme of Minkowski. For more detail see Filippenko (1997) or Turatto et al. (2007). The great diversity of core-collapse supernovae (CCSNe) is understood as the result of a rich variety of parent systems. Initial differences in mass, radius, metallicity or rotation, and evolutionary differences in mass lost to stellar winds or interacting binary companions, would result in a wide distribution of envelope masses when the progenitor stars reach the time of core collapse (e.g., Heger et al., 2003; Kasen & Woosley, 2009; Dessart et al., 2013; Pejcha & Prieto, 2015). According to the chemical composition of the outer layers at the time of explosion the spectroscopic display will be of Type II, or Type IIb, Ib or Ic, with little or no presence of hydrogen. The latter are collectively called “stripped envelope SNe”. Finally, according to the total mass of hydrogen in the envelope, a bona fide Type II SN will show a slower or faster rate of decline after maximum and will be named “plateau” (IIP) or “linear” (IIL, Barbon et al., 1979). The convention has stuck although we know now that there is a continuous distribution of decline rates (e.g., Anderson et al., 2014; Sanders et al., 2015; Pejcha & Prieto, 2015; Galbany et al., 2016a) and that Type IIP and IIL have similar progenitors (Valenti et al., 2015). Extreme mass loss shortly before explosion may lead to the formation of a dense circumstellar medium (CSM). The interaction of SN ejecta with the CSM would produce narrow emission lines of hydrogen and the SN is named Type IIn in these cases (e.g., Dopita et al., 1984; Schlegel, 1990; Stathakis & Sadler, 1991; Chugai, 1994).

Progenitors of CCSNe have been identified in pre-explosion images (Smartt et al., 2009, 2015), and found to be red supergiants (RSG) with Zero Age Main Sequence (ZAMS) masses between and  M for Type II SNe or objects with masses between and  M at the time of explosion after a complex interacting binary evolution for Type IIb and Ib SNe (Folatelli et al., 2015; Eldridge et al., 2013; Folatelli et al., 2016).

As CCSNe progenitors are massive stars with relatively short lifetimes of  Myr, it is expected that the SNe are still associated with their birth places, namely spiral arms or H II regions (Bartunov et al., 1994; McMillan & Ciardullo, 1996), although CCSNe in early-type galaxies with residual star formation have been reported in the literature (Hakobyan et al., 2008).

The correlation between CCSNe and star forming regions is expected to decrease toward the lower-mass progenitors and get diluted by progenitors resulting from binary evolution which can result in longer timescales before explosion (Zapartas et al., 2017; Eldridge et al., 2017). Recent work has found that stripped envelope SNe are more closely associated with star forming regions than hydrogen rich SNe, as expected from the increasing progenitor mass of the sequence (e.g., Anderson et al., 2012; Galbany et al., 2014).

Observations of the SN host metallicity, either global (Prieto et al., 2008; Arcavi et al., 2010) or local (Modjaz et al., 2008; Anderson et al., 2010; Stoll et al., 2013; Galbany et al., 2014; Kuncarayakti et al., 2018), and the spatial distribution of CCSNe (e.g., Van Dyk et al., 1999; Petrosian et al., 2005; Mikhailova et al., 2007; Kangas et al., 2017), also provide constrains on the different progenitors. Chemical abundance studies generally associate stripped envelope SNe with higher metallicity hosts as expected from the probable increase of mass loss with metallicity (e.g., Henry & Worthey, 1999; Sánchez et al., 2014; Sánchez-Menguiano et al., 2016). Studies of the radial distribution of SNe find that Type Ib/c SNe are more centrally concentrated in their hosts (Anderson & James, 2009; Hakobyan et al., 2009), which is again consistent with the higher metal enrichment towards the center of galaxies. And studies of the height distribution of SNe in edge-on disk galaxies (Hakobyan et al., 2017), find that CCSNe are nearly twice as much concentrated towards the disk than thermonuclear SNe, a result consistent with the height scale of the stellar populations where their progenitors are expected to originate.

Some CCSNe, however, defy the common sense implicit in the previous description by appearing far from any identifiable birthplace. One striking example is SN 2009ip located  kpc from NGC 7259, the nearest spiral galaxy (e.g., Fraser et al., 2013; Mauerhan et al., 2013; Pastorello et al., 2013; Prieto et al., 2013). SN 2009ip was first identified as a SNe impostor in 2009 three years before exploding as a Type IIn supernova (Smith et al., 2014). Late time HST data rules out the presence of star forming regions comparable to Carinae or the Orion Nebula (Smith et al., 2016) at the explosion site. The possibility that the progenitor was a runaway star is also rejected. The peculiar velocity of the SN is smaller than 400 km/s and the high mass of the progenitor implies a lifetime shorter than the travel time from the nearest star forming site located at  kpc.

This paper introduces another example of this sort, the Type IIP like supernova ASASSN-14jb in the edge-on disk galaxy ESO 467-G051 (Brimacombe et al., 2014; Challis, 2014; Zhang & Wang, 2014). The SN exploded at 2.5 kpc from the center and 2.1 kpc above the galactic disk where no significant star forming region is detected. Also, it is interesting that ESO 467-G051 and NGC 7259 (the host of SN 2009ip) form an interacting pair. The SNe are separated by 2.4’ in the sky ( kpc).

We present here photometric and spectroscopic observations of ASASSN-14jb and Integral Field Spectroscopy (IFS) of its explosion site. We perform a thorough analysis to estimate physical parameters of the SN, the progenitor star, and the parent galaxy. The paper is organized as follows. In Section 2 we describe the photometric and spectroscopic observations. Section 3 contains our comparative analysis of the SN, including the photometric and spectroscopic evolution of the ASASSN-14jb. Section 4 contains the estimates of basic physical parameters of the SN, Section 5 includes the analysis of the host galaxy and its H II regions. Finally, Section 6 includes our discussion and Section 7 our conclusions.

2 Data

2.1 Discovery and Explosion Time

ASASSN-14jb was discovered on UT 2014-10-19.09 (MJD = 2456949.09,  mag, Brimacombe et al. 2014) at RA = 22:23:16.12, DEC = -28:58:30.78 (J2000.0) by the ongoing All Sky Automated Survey for SuperNovae (ASAS-SN; Shappee et al., 2014; Holoien et al., 2017a) from the “Cassius” 4-telescope unit at the Cerro Tololo Inter-American Observatory (CTIO), in Chile, hosted by the Las Cumbres Observatory (LCOGTN; Brown et al., 2013). It was spectroscopically classified as a young Type II on UT 2014-10-20 (Challis, 2014; Zhang & Wang, 2014). The last non-detection was reported to be on MJD with a 3 magnitude upper limit of  mag, 6 days before the time of discovery. We take the midpoint between this last non-detection and the discovery time, (MJD), as the “explosion” time. The uncertainty is estimated as half the interval between the last non-detection and the discovery epoch.

Figure 1: -band image of the field of ASASSN-14jb from LCOGTN. The circles show the position of ASASSN-14jb, the local photometric standards (numbered 1-10), and the nearby SN 2009ip in NGC 7259. The zoom around the position of ASASSN-14jb and its host galaxy, ESO 467-G051, shows a synthetic -band image obtained from the MUSE datacube. The MUSE image field of view is of 1.5x1 arcmin.

2.2 Photometry

Optical photometric observations of ASASSN-14jb were obtained with the ASAS-SN unit “Cassius” at CTIO in the -band and the LCOGTN 1-meter telescopes at CTIO, the Siding Spring Observatory (SSO), the South African Astronomical Observatory (SAAO), and the McDonald Observatory (MDO) in the and the Sloan Digital Sky Survey (SDSS) filters. All the ASAS-SN images are processed in an automated pipeline using the ISIS image subtraction package (Alard & Lupton, 1998; Alard, 2000), with further details given in Shappee et al. (2014). Using the IRAF111IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. apphot package, we performed aperture photometry on the subtracted images and then calibrated the results using the AAVSO Photometric All-Sky Survey (APASS; Henden et al., 2012). The ASAS-SN photometry is presented in Table 2.

The fully reduced LCOGTN 1m images (bias/overscan subtracted and flat-fielded) were retrieved from the LCOGTN data archive. We solved the astrometry of each CCD frame using the software package (Lang et al., 2010). We obtained calibrated magnitudes for local standard stars in the field from the APASS catalog. In Figure 1 we display one of the LCOGT images showing the position of the SN, the local standards, and the nearby SN 2009ip. The data of the local standard sequence is given in Table LABEL:tab:std. ASASSN-14jb was recorded with high signal-to-noise ratio (SNR) in the images of LCOGTN. This, together with the negligible background from the host galaxy at the SN site, prompted us to measure the SN flux using aperture photometry. The procedure is explained in detail in Appendix A.

As our observations extended only days after explosion, we searched the data in public domain for images of the field which could have recorded ASASSN-14jb at later times. We found some images in the ESO archive222 We retrieved ESO/NTT EFOSC images in filters that contained ASASSN-14jb’s explosion site at three different epochs during 2015, obtained by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO; Smartt et al., 2015, ESO program ID 191.D-0935). The CCD frames were overscan subtracted and flat-fielded with calibration frames obtained the same day using standard IRAF routines. Since the SN is fainter in these late-time images, we measured the SN brightness using PSF fitting photometry as implemented in IRAF DAOPHOT package. We calibrated the magnitudes using local standards from Pastorello et al. (2013). The early LCOGTN and late ESO/NTT photometry are presented in Table LABEL:tab:allphot.

We supplemented our ground-based optical photometry with data from the Swift and Spitzer space-born observatories. Repeated optical and near-UV observations of ASASSN-14jb were obtained by Swift/UVOT in the passbands w2,m1 w1,u,b,v during the first days after explosion. We retrieved them from the Swift Optical/Ultraviolet Supernova Archive (SOUSA; Brown et al., 2014). Also, the field of ASASSN-14jb and its host have been observed nine times since 2010 in the 3.6 m and 4.5 m bands by the Spitzer IRAC instrument (Fazio et al., 2004). The observations obtained on 2015-09-13 (program ID 11053) and 2016-08-17 (program ID 12099) were taken when the SN was and days past explosion, respectively. We retrieved these images from the data archive333 and also images taken on 2014-09-05 (program ID 10139) to use as templates for image subtraction. We used HOTPANTS (Becker, 2015) for difference imaging and measured aperture photometry on the difference images. The SN is detected on 2015-09-13 and undetected on 2016-08-17.

Our near-UV and optical photometry is displayed in Figures 2 and 3. The first shows the near-UV and optical light curves for the first days after explosion, and the second the complete light curves in and bands. The Spitzer mid-infrared photometry is discussed in § 4.5.

Figure 2: ASAS-SN , LCOGTN , and Swift-UVOT , , , , light curves of ASASSN-14jb. Each band has been shifted to the value showed in the label for presentation purposes. The observed time is given in days since the estimated explosion epoch (MJD) in the bottom x-axis and in MJD in the top axis. The black squares are the non-detection derived from the ASAS-SN -band images.
Figure 3: Same as Figure 2 for the EFOSC and late time photometry. The decay slope is also showed.

2.3 Spectroscopy

A total of ten single-slit, low-resolution optical spectra of ASASSN-14jb were obtained with the FAST spectrograph (Fabricant et al., 1998) mounted on the Fred L. Whipple Observatory Tillinghast 1.5-m telescope ( Å, ) and with the Inamori-Magellan Areal Camera and Spectrograph (IMACS; Dressler et al., 2011) mounted on the Magellan Baade 6.5m telescope at Las Campanas Observatory ( Å, ). The CCD images were reduced and the 1D spectra were extracted and calibrated using standard routines in the IRAF packages ccdproc and onedspec, respectively. The epochs of the ten spectra obtained in the photospheric phase are presented in Table LABEL:tab:spectra. A montage of the spectra obtained in the photospheric phase is shown in Figure 4.

Figure 4: The full spectroscopic time series, in logarithmic scale and corrected by redshift (), for ASASSN-14jb in the photospheric phase (up to 80 days past explosion). The time next to each spectrum, in days, is measured since the estimated time of explosion (MJD). The vertical lines show the wavelengths of the some of the main features in the spectra of type II SNe. Lines with the symbol represent the position of strong telluric absorptions.

We also observed the field of ASASSN-14jb as part of the All-weather MUse Supernova Integral field Nearby Galaxies (AMUSING; Galbany et al. 2016b) with the Multi Unit Spectroscopic Explorer (MUSE, Bacon et al. 2010) on ESO’s Very Large Telescope UT4 (Yepun). MUSE is a state-of-the-art integral field spectrograph with a field of view of and spaxels. It covers the spectral range  Å with a resolving power . Our MUSE data were obtained on 2015-11-14 and consisted of two different pointings of four dithered exposures each with an integration time of 698 sec. The sky conditions were clear, and we measure a full-with-half-maximum for the stellar point-spread function of 108 at 6600 Å. We reduced the MUSE spectroscopy with version 1.2.1 of the EO(fpipeline (Weilbacher et al., 2014), but also checked and corrected the astrometric zeropoint using LCOGTN images (for further detail on MUSE data reduction see Prieto et al. 2016).

Although the MUSE observations were obtained  days after explosion (see Table LABEL:tab:spectra) when ASASSN-14jb was faint, the quality and depth of the exposures allow us to easily detect the SN and extract a nebular phase spectrum. We used a circular aperture with a radius of to extract the spectrum at the position of the SN with QFitsView444 To take care of any offset in the flux scale we estimated the magnitude of the SN at the time of the spectrum from the late time light curve () and compared it with the synthetic -band magnitude obtained from the nebular spectrum (). We then applied an offset of 0.35 mag (factor of 1.38 in flux) to obtain an approximate absolute flux calibration. The nebular spectrum of ASASSN-14jb is displayed in Figure 11.

3 Comparative Analysis

To put ASASSN-14jb in context we assembled a comparison sample representative of the wide range of properties of Type II SNe. These include the luminous Type IIP SN 2007pk (Inserra et al., 2013), normal Type IIP like SN 1999em (e.g., Hamuy et al., 2001), and the subluminous Type II SN 2005cs (e.g., Pastorello et al., 2009). We chose objects with data in the public domain, good coverage of the photospheric phase, ideally including photometry (8 SNe), and optical nebular spectra at days after explosion (6 SNe). Further details of the comparison set are given in Table 5.

Estimates of the distance and total foreground reddening are needed in order to compare ASASSN-14jb with the rest. We use a distance modulus of  mag, to be justified below.

We estimate extinction by dust in the Galaxy using the maps of Schlafly & Finkbeiner (2011). Taking the average in a circle with 5 arcminutes of diameter centered on the SN position we obtain  mag. Extinction by foreground dust in the host galaxy is expected to be low, as the SN exploded at a significant distance from the host galaxy disk. Consistent with this, the narrow Na I D doublet at the host galaxy redshift are not detected in our spectra, which is a clear sign of low reddening (Phillips et al., 2013). In addition, the MUSE data to be discussed below shows no significant star forming region near the explosion site. We conclude, hence, that ASASSN-14jb is affected only by Galactic extinction. Assuming a standard reddening law with a ratio of total to selective absorption , we estimate  mag.

3.1 Photometric Evolution

Figure 5 shows the absolute -band light curve of ASASSN-14jb together with those of the comparison set. The peak at  mag is within the usual range of Type IIP SNe, although lower than the average  mag in the sample of Type II SNe of Anderson et al. 2014, (from now onward A14). ASASSN-14jb is dimmer than the classic, well-studied Type II-SN 1999em by  mag, but it shows a similar behavior in its -band plateau, including the early double peak (Hamuy et al., 2001; Leonard et al., 2002). Up to days since explosion the photometry shows the behavior typical of Type IIP like SNe.

Figure 5: Absolute -band magnitude evolution for ASASSN-14jb (black filled stars) and the comparison sample of Type II SNe as a function of days after their estimated time of explosion. See Table 5 for references of each SNe, including distance and reddening adopted.
Figure 6: Extinction corrected colors for ASASSN-14jb (black filled stars) and the comparison sample of Type II SNe. Reddening values used are in Table 5.
Figure 7: Extinction corrected , , and colors for ASASSN-14jb (black filled stars) and the Type II SNe comparison sample. Reddening values used are in Table 5.

The measured slope of the plateau ( days after explosion) in the -band for ASASSN-14jb is  mag per 100 days, lower than the average -band plateau slope in the A14 sample ( mag per 100 days). The slow increase in brightness of the plateau is consistent with the correlation between plateau slope and absolute -band magnitude found by A14 and later work (e.g., Valenti et al., 2016). It is also an effect seen in bolometric light curves of CCSNe models resulting from the explosion of low mass progenitors (Sukhbold et al., 2016).

In Figures 6 and 7 we present the extinction corrected color evolution of ASASSN-14jb and the comparison set. The color of ASASSN-14jb is bluer than average at  days after explosion, being comparable only to LSQ 13fn (Polshaw et al., 2016) and SN 2007pk (Inserra et al., 2013) at 80 days after explosion. The change in color from the onset of the plateau up to 80 days is  mag, low when compared with the comparison sample. The color evolution of ASASSN-14jb in , and also indicates a comparatively bluer continuum, except in the color. This behavior may be due to a true difference in the temperature evolution or it may be an effect of lower line-blanketing in the blue part of the spectrum, resulting from low-metallicity of the progenitor. The similarity in the color evolution of ASASSN-14jb and LSQ 13fn, a Type II SN that probably comes from a low-metallicity progenitor (Polshaw et al., 2016), might suggest that a similar physical mechanism is at play.

3.2 Spectroscopic Evolution in the Plateau phase

As shown in Figure 4, the spectral evolution of ASASSN-14jb resembles the evolution of other hydrogen-rich Type II SNe. At very early times the ejecta are optically thick and the spectra show a hot blue continuum with relatively weak, but broad, Balmer lines without P-Cygni absorption troughs. In the first spectrum at  days after explosion, we detect broad Balmer lines in emission and the high ionization He II line. The FWHM of the Balmer lines and the He II lines in this earliest spectrum are  km/s and  km/s, respectively, and all the emission peaks are clearly blueshifted by  km/s. The strength of the He II line decreases substantially in the spectrum obtained a day later, at 5 days after explosion, and the Balmer lines start to develop P-Cygni profiles. We also detect the He I  line in the early spectra. At 6 days after explosion, the He II line has weakened further and it is undetected in the spectrum obtained at 7 days after explosion.

At phases later than  days, the ejecta cools and the hydrogen lines develop strong P-Cygni profiles. The blueshift of the peak of the Balmer lines decreases to  km/s. This is a common feature in Type II SNe during the photospheric epochs and it has been attributed to the shallow density profile of the ejecta compared to a typical stellar wind (Dessart & Hillier, 2005; Anderson et al., 2014). Lines from transitions in metals such as Fe II, Sc II, and Ti II also start to appear and significantly contribute to line blanketing. The Na I D doublet at appears now in the wavelength range where He I  was detected at early times. Comparison of the photospheric phase spectra of ASASSN-14jb with those of the comparison sample at days past explosion (see Figure 8) shows that the absorption features of the former in the range  Å are relatively weak.

In Figure 9 we present the time evolution of the H, H and Fe II lines. The last two spectra we obtained in the plateau phase, at 64 and 80 days after explosion, show clear absorption features on the blue side of both the H and H main absorption troughs. These high velocity (HV) features are at  km/s, while the main absorption features are at  km/s. The HV absorption lines cannot be explained by Si II  or Ba II absorption features and have been observed in of the Type II SNe studied by Gutiérrez et al. (2017).

Figure 8: Comparison of rest-frame photospheric phase spectra at days after explosion. The spectra are not corrected for extinction.
Figure 9: Evolution of the H (left), Fe II  (middle) and H (right) profiles, in velocity space, as a function of time after explosion (MJD), in the plateau phase for ASASSN-14jb (up to 80 days past explosion). In each panel, a vertical dashed line indicates zero velocity.
Figure 10: Expansion velocities of the H (top), Fe II  (middle), and H (bottom) lines as a function of time after explosion in the plateau phase, measured from the P-cygni profile minima of the Balmer lines and the Fe II line for ASASSN-14jb (black filled stars) and the Type II SNe comparison sample.

3.3 Expansion velocities

The ejecta of the Type II SNe achieve near homologous expansion after a few days (Bersten et al., 2011) and at early times the hydrogen atmosphere is opaque and the line forming region is in the outer layers at high velocities. As the temperature drops, the opacity decreases and the photosphere recedes in mass (Lagrangian) coordinate, and appears at lower velocities. The photosphere is tightly bound to the opacity drop caused by hydrogen recombination (Bersten et al., 2011). In practice one estimates the photospheric velocity by measuring the blueshift of maximum absorption of an specific P-Cygni profile. This can over or under estimate the true photospheric velocity (Dessart & Hillier, 2005, 2010; Takáts & Vinkó, 2012). Other methods used are the cross-correlation with library spectra, the comparison with detailed NLTE codes like CMFGEN or PHOENIX (Dessart & Hillier, 2005; Baron et al., 2005) or comparison with more simplified, parametrized LTE spectral modeling (e.g., SYNOW, Takáts & Vinkó 2012).

The expansion velocity measured from optical spectra is an important physical parameter directly related to the explosion energy and the plateau luminosity. These are also fundamental for distance measurements (Hamuy & Pinto, 2002). Strong Balmer lines and Fe II lines velocities were measured in our spectra. The spectra were first normalized to a global continuum fitting a black body or a local power law. Then, when possible, low degree polynomials were fitted to the P-Cygni profile absorptions and the minimum was taken as a proxy for the expansion velocity. We did not apply corrections by reddening or peculiar velocities. Details of the P-cygni profile fitting are in Appendix B. In Figure 10 we show the expansion velocities of ASASSN-14jb together with those of the comparison set. For LSQ 13fn we only present the velocities after 30 days due to the low signal to noise of the spectra. The velocities of SN 2013ab are directly taken from Bose et al. (2015). The measured velocities are also collected in Table 4 together with the average velocity of our comparison set and those of Gutiérrez et al. (2017).

To compare with well known correlations (Hamuy, 2003; Faran et al., 2014), we interpolated the velocities of our sample at 50 days after explosion using a Monte-Carlo approach. For each SN we generated a re-sampling of each velocity using the measured errors and following a non-linear least squares procedure, using the curve fit function from the NumPy library, we fit a power law model each time. For this procedure we select only velocities after 20 days, as early velocities do not follow the same decay rate as in the plateau. We then take the average and standard deviation of the re-samplings as the 50 days velocity and error, respectively. For all the lines measured the velocities for ASASSN-14jb are slightly under average. The Fe II velocity interpolated at 50 days is measured to be km/s, compared to the average of km/s of our comparison sample. Our sample average velocity is in agreement with the value obtained with a larger sample in Gutiérrez et al. (2017), of km/s at 53 days after explosion. The photospheric velocities inferred are in agreement with the velocity-luminosity correlations (see Section 3.2 below).

As is common in Type II SNe distance measurements (Hamuy & Pinto, 2002; de Jaeger et al., 2017), we fitted a power law to interpolate the velocities. The Fe II logarithmic decay presents an average of in our sample, compared to the value obtained for ASASSN-14jb, of . We observe a strong anti-correlation between the velocity decay slope and the -band absolute magnitude at 50 days. This is expected as there is an internal correlation between both parameters of the power law (de Jaeger et al., 2017) and the velocity scale correlates with the luminosity (Hamuy & Pinto, 2002).

In conclusion, ASASSN-14jb presents expansion velocities consistent, both from the scale and velocity decay, with below average luminosity Type II SNe.

Figure 11: Comparison of nebular spectra of observed Type II SNe at days past explosion.

3.4 Nebular Spectrum

The nebular spectrum of ASASSN-14jb is shown in Figure 11 together with a sample of Type II SNe. It shows the typical emission lines of Type IIP like SNe at nebular phases: H, [Ca II] , [O I] , Ca II , [Fe II] , with a boxy profile, Na I and [O I] , in order of observed intensity (Silverman et al., 2017). Also, a strong [C I] line is present, which has been detected in other Type IIP like SNe such as SN 2004dj (Silverman et al., 2017). A weak He I line at is observed, but with no He I counterpart at . The clear separation observed in different multiplets like [OI], [CaII] and Ca II + [CI] is a sign of a low explosion energy for ASASSN-14jb.

The blueshift with respect to rest-frame of  km/s observed in the peaks of the strongest lines is particularly interesting. It cannot be caused by the same physical effect as in the plateau because in this optically-thin phase the outer density profile is less relevant to the emission (Anderson et al., 2014). Other hypotheses to explain this blueshift may be: 1) a shift in velocity relative to the host due to dynamical effects or an intrinsically high-velocity of the progenitor star; 2) dust production; 3) Some lines are still optically thick at 400 days past explosion; or 4) an asymmetry in the explosion itself. Notable is the double peak observed in which is not observed in other lines, although this is also seen in the nebular spectra of other Type IIP like SNe (Elmhamdi et al., 2003; Chugai et al., 2005; Utrobin & Chugai, 2009).

4 Distance and Physical parameters of ASSASN-14jb

Progenitors of Type IIP like SNe have been constrained by pre-explosion images to be RSGs with ZAMS masses of M (Smartt et al., 2009, 2015). The discrepancy between these initial mass constraints and the predictions from evolutionary codes (Heger et al., 2003; Ekström et al., 2012) has been called the “RSG problem”. The data set we have compiled allow us to constrain the basic physical parameters of the progenitor of ASSASN-14jb, compare them with the expected values for SN Type IIP like progenitors, and contribute to this discussion. But doing so requires a good estimate of the distance to the SN.

4.1 Distance

We measure the distance towards ASSASN-14jb using some of the most recent empirical methods: The Photospheric Magnitude Method (Rodríguez et al., 2014), the Standard Candle Method (Hamuy & Pinto, 2002), and the Photometric Color Method, (de Jaeger et al., 2017). Further details of our hypothesis and calculations are given in Appendix C.

The three different estimates of distance are shown in Figure 12. They are consistent with each other within each method’s intrinsic dispersion. They are also consistent with the distance obtained from the Virgo-inflow corrected recession velocity of the host galaxy and the Hubble law, assuming . Given this, we simply take as the final distance the weighted average, considering the statistical errors only, of  mag ( Mpc).

Figure 12: Distance modulus estimates for ASASSN-14jb using the PMM, SCM and PCM methods (green dots). The black horizontal line is the distance modulus from the host galaxy redshift and a Hubble law with km/s/Mpc, and the grey area shows a dispersion of km/s as a reference error for the CMB frame redshift. The thick green line is the average from the Tully-Fisher measurements obtained from the NED database and the green area represent the standard deviation of the measurements, all normalized to km/s/Mpc. The thick blue line is the adopted distance modulus in this work, , obtained as the weighted average of the PMM, SCM and PCM methods and the standard deviation is showed as a blue fill. We also show as a reference the distance measured to SN 2009ip in NGC 7259 (Potashov et al., 2013).

4.2 Temperature Evolution and the Progenitor Radius

The early emission, a few days after explosion, for a Type II SN comes from the adiabatic cooling of the outer envelope and depends mainly on the pre-supernova radius (Rabinak & Waxman, 2011; Sapir & Waxman, 2017). Observationally, the pre-supernova radius has been constrained from early UV and optical photometry, using semi-analytic models (Rabinak & Waxman, 2011; Nakar & Sari, 2010), on individual Type II SNe (Gezari et al., 2010; Valenti et al., 2014; Bose et al., 2015; Huang et al., 2016, 2018; Tartaglia et al., 2018), and on larger samples to make statistic studies (González-Gaitán et al., 2015; Gall et al., 2015; Rubin et al., 2016).

We build the spectral energy distribution (SED) for ASASSN-14jb at early times, from 4 to 26 days after explosion to derive an estimate of the progenitor radius. Using the Swift/UVOT near-UV and optical photometry, we build the SED of ASASSN-14jb at each epoch, after correcting for extinction assuming an standard Galactic extinction law with total to selective extinction ratio (Cardelli et al., 1989). We fit a blackbody function to the SED and obtain the photospheric temperature and radius evolution, assuming the distance to the host galaxy obtained in Section 3.2.

The results are presented in Figure 13. The observed photospheric temperatures are in the range 5-21 and the radii are in the range 2-18 R. We observe a change in the evolution  days after explosion. The observed temperature logarithmic decay increases from to and the photospheric radius growth rate increases from to .

Rabinak & Waxman 2011, hereafter RW11, developed analytic expressions for the early evolution of the photospheric radius and temperature , expected to be valid before hydrogen recombination is significant. The expressions depend on progenitor properties at the time of explosion such as radius , opacity , and density profile. They also depend on the explosion energy and ejecta mass .


where for RSGs, depends on the progenitor density profile, is the ejected mass in solar masses, is the explosion energy in units of  erg, is the opacity normalized by 0.34, is the progenitor radius in units of  cm and is the time since shock breakout in units of seconds.

It is well known that the emission from Type II SNe differs from a pure black body (Eastman et al., 1996; Dessart & Hillier, 2005). When the atmosphere is highly ionized, the electron scattering opacity dominates and the photon thermalization layer differs from the photosphere, the former being underneath at higher temperatures. Because of this the observed temperature (or color temperature , as it is usually called) can differ significantly from the photospheric temperature . This effect, together with the increasing line blanketing and departures from the plane parallel atmosphere approximation, can make the observed SED significantly different from a black body. This effect can be corrected using a theoretically estimated dilution factor , so that the emergent flux is defined as,


where is the luminosity distance and is the Planck function (Eastman et al., 1996; Dessart & Hillier, 2005). The dilution factor is in general wavelength dependent and directly affects the angular diameter distance estimate to a supernova. However, as we take the luminosity distance from empirical methods independent of 3, the dilution factor only affects the photospheric radius estimation and not the temperature. Note, however, that in the RW11 scheme , which translates in a factor of increase in the photospheric radius (see equation 38 in RW11). In our case we choose a representative value of for this high temperature regime, as obtained in Dessart & Hillier (2005); Pejcha & Prieto (2015).

To test the effect of the deviations from a black body, we used a semi-empirical approach. We first selected the 15  progenitor mass supernova explosion model from Dessart et al. (2013) that best fits the observed near-UV colors (named m15z8m8). Then, at each epoch we fit the SED using,


a model built from the weighted geometric mean of a black body and the normalized explosion model. The normalization ensures that the model’s intrinsic brightness does not affect the fitted parameters. In principle, the weighting parameter is unconstrained but to ensure a smooth correction we fixed its value to be between and selected the model with the minimum residuals. Figure 13 show our results for both approaches. We see that both temperature and radius evolution become smoother when including the models, probably because the hybrid approach accounts better for line blanketing. Since the available models start from  days past explosion, all the SEDs before that time were fitted with a single black body. The derived slope in the temperature changes to after 10 days. The photospheric radius follows roughly a linear expansion with a logarithmic slope of during the same epochs.

Now we apply the RW11 models using the modified black body. We fitted the temperature evolution assuming the fiducial values of , , , and use only temperatures greater than  eV because RW11 inform that their models are more reliable in this range. The radius we obtain is  R. Noting the impact of the time since explosion in the result, we varied the shock-breakout (SBO) epoch and look for the one that minimizes the of the fit. We then obtain  R, for a SBO occurring  days earlier than our initial estimate (reduced ).

Inverting Equation 2 is another form of understanding possible variations of :


where is the ratio of the color temperature and true photospheric temperature. Following RW11 we take as a good approximation at early times 555We note that this equation is a corrected version of Equation 33 in RW11. The original one is , which is inconsistent with Equation 13 of the same paper (our Equation 1)..

The values obtained using Equation 5 with early time observations are consistent with  R. As expected, the radius decreases as the actual decay of the photospheric temperatures is faster than the model, and the varies due to the changing opacity of the ejecta.

We choose as our best estimate of the progenitor radius, the RW11 fit with an SBO offset of  days, which gives an estimate of  R. This value is consistent with the color evolution of the mildly sub-solar (Z = 0.4 Z) and compact ( R) RSG progenitor (m15z8m3) used in Dessart et al. (2013). As the photospheric radius growth strongly depends on the explosion energy to ejecta mass ratio (equation 1) we can estimate the ejecta mass as a function of the explosion energy. Fixing , as before we obtain .

Figure 13: Black-body temperature (blue dots) and radius evolution (red dots) as inferred from the SED of data. The SED was corrected by a reddening value of  mag. The horizontal line indicates the 1 eV level, where recombination is expected to become relevant afterwards. The vertical black solid line indicates the assumed explosion epoch at MJD and the vertical dashed lines is at the assumed offset () in the explosion epoch, for the optimal cooling model from RW11. The red dashed lines shows a power-law fit to the photospheric radii (only using the data points with ) while the blues dashed lines shows the same for the Temperatures plus a line with a factor of 1.2 (i.e., the lines shows a fit to ). The black curved dashed lines are the fits assuming the RW11 for and .

4.3 Bolometric Light Curve and the Ni Mass

We constructed a bolometric light curve of ASASSN-14jb, by applying the semi-analytic bolometric corrections from Pejcha & Prieto (2015) for the extinction corrected colors. Figure 14 shows the result, together with the bolometric light curves of the comparison sample. The bolometric luminosity of ASSASN-14jb  days after explosion is , which is under-luminous compared to SN 1999em, SN 2012aw or SN 2013ab.

The late-time bolometric light curve can be used to estimate the Ni and total masses ejected in the explosion. We assume a model of partial trapping of gamma rays (Clocchiatti & Wheeler, 1997), where the energy released by Cobalt decay is scaled by the transparency factor , and the characteristic time is given by . We assumed a fixed density slope and constant gamma-ray opacity (which involves assuming an electron fraction per nucleon of ). The late-time light curve and fitted model are plotted in Figure 15. The best fit provides  days, which would correspond to an ejecta mass of for an explosion energy of  erg (1 foe). The Ni mass results  M, which is slightly lower than the median of 0.031 M corresponding to the Type II SNe sample studied by Müller et al. (2017).

Late-time light curves that decay faster than Cobalt are expected in SNe with low mass ejecta, as in the case of stripped enveloped SNe and some fast declining Type II SNe (e.g., SN 2014G, SN 2013hj and ASASSN-15nx, Bose et al., 2016; Terreran et al., 2016; Bose et al., 2018). Given the low energetics shown in, for example, the expansion velocities, the ejected mass estimated above is consistent with the low mass range of Type IIP SNe progenitors.

The ejected mass estimated with the late-time light curve can be compared with that obtained in the previous section using the early light curve. The scaling in this case is  M which implies an “agreement” pair () at (6 M,0.25 foe).

Anderson et al. 2014 discusses the observations of higher deviation from the expected full trapping decay for the brightest Type II SNe in their sample (in the -band) and they suggest that this would be explained by the expected more diluted (i.e., less massive ejecta or more extended) hydrogen associated with faster declining SNe. In our case, the ejecta mass is slightly below what would be needed for a relatively low energy explosion of 0.25 foe to have a high gamma ray opacity () at days.

We can also estimate the ejected Ni mass using its correlation with the FWHM of H (Maguire et al., 2012). A value of  Å gives us a Ni mass of  M. This is lower than, but consistent with, the estimate obtained from the late-time bolometric light curve.

The low nickel mass agrees with the explosion parameters according to the classic explosion models for Type II SNe (Popov, 1993; Litvinova & Nadezhin, 1985) and observed correlations (Hamuy, 2003).

Figure 14: Comparison of bolometric light curves using the bolometric corrections of Pejcha & Prieto (2015). Each SN was corrected for extinction. See Table 5 for references of the comparison sample.
Figure 15: Bolometric light curves in the radioactive tail derived from the and -band photometry from EFOSC imaging, using the bolometric corrections of Pejcha & Prieto (2015) (PP15, blue dots) and Bersten & Hamuy (2009) (MH09, green squares). The decay is faster than the expected from full trapping of Co decay (showed as a dashed black line). The non-full trapping model from Clocchiatti & Wheeler (1997) gives a Ni Mass of for the Pejcha & Prieto (2015) bolometric corrections. The time scale for this fit is days, which corresponds to a ejecta mass of 12 for an explosion energy of 1 foe.

4.4 Nebular Spectra and the Progenitor Mass

In the nebular phase, spectra of Type IIP like SNe show prominent H, O I, Ca II and Fe II emission lines. The [O I] is used to estimate the Oxygen core mass and therefore the main sequence mass of the progenitor star (e.g., Uomoto, 1986; Fransson & Chevalier, 1989; Elmhamdi et al., 2003; Jerkstrand et al., 2012; Kuncarayakti et al., 2015). We compare our spectra with models from Dessart et al. (2013), Jerkstrand et al. (2014) and Jerkstrand et al. (2018). (Dessart et al., 2013, hereafter D13), in particular, presented radiative transfer models for Type II SNe using a grid of 15 M stellar models from MESA STAR, with varying mixing length, metallicity, and rotation. We choose the nebular spectra in their “m15” model, which reach up to 400 days past explosion (see Table 1 in D13 for details). (Jerkstrand et al., 2014, hereafter J14), calculated models for 12, 15, 19 and 25 M stars evolved using the KEPLER code, assuming Solar metallicity and no rotation. We also included the recent model of 9 M from (Jerkstrand et al., 2018, hereby J18).

In Figure 16 and 17 we compare the data and the models for the [O I] and H lines. In the case of D13 models m15z2m2 and m15z8m3 provide the best fit. They correspond to a progenitor metallicities of and , oxygen mass of and  M, and nickel masses of and  M, respectively (for details see Table 1 in D13). In the case of the J14 and J18 models the 9 M and 12 M are closest to the data, although all the models over-predict the H emission line. A problem with this direct comparison is that the explosion models presented in D13, J14 and J18 produce higher Ni masses than measured for ASASSN-14jb. If the spectra are normalized by the H emission to approximately take care of this difference, the J14 models of 15 M  and 9 M become closest to the data. On the other hand, the D13 models reveal a high degeneracy given their differences in stellar evolution.

Other approach we can use to alleviate the impact in the spectra of the Ni mass discrepancy between ASASSN-14jb and the models, is to normalize the flux in a given line by the Ni decay power. Doing this for the [O I] line, we obtain a ratio of . This leads to a surprising better consistency with the higher mass range of J14 models ( M), but again to a better fit of the 9 M. According to Jerkstrand et al. (2018) this results both from the oxygen shell in the 9 M  model lying closer to the Ni and a significant contribution of Fe I to the [O I] doublet.

As a final estimate we calculate an upper limit to the emitting oxygen mass using the analytic formula provided by J14,


where is the luminosity of the [O I] doublet in erg/s, is its Sobolev escape probability, and the equilibrium temperature in K. A temperature in the range 3900–4300 K was estimated from the ratio between the lines [O I] and [O I] 666Assuming as in J14 and estimating the fluxes by fitting a skewed Gaussian to the line profiles.. The resulting values of fall in the range 0.09–0.18 M, which is above that of the 9 M model of J18 and 1.5–3 times lower than that of the 14 M model of J14.

We conclude that the nebular spectra confidently points at a low mass progenitor for CCSNe standards, in the range  M. This is supported by surveys of theoretical explosions spanning a wide grid of progenitor properties (e.g., Sukhbold et al., 2016), which show that models above M are much more efficient at producing oxygen than below that ZAMS mass.

Figure 16: Comparison of the nebular spectra of ASASSN-14jb with the NLTE models of Dessart et al. (2013) at days past explosion. The vertical dashed lines shows the rest wavelength of H and the [O I] 6300, 6364 doublet.
Figure 17: Comparison of the nebular spectra of ASASSN-14jb with the NLTE models of Jerkstrand et al. (2014, 2018) at days past explosion. The vertical dashed lines shows the rest wavelength of H and the [O I] 6300, 6364 doublet.

4.5 Warm Dust in the Ejecta

ASASSN-14jb is clearly detected in the 2015-09-13 Spitzer images, with  mag and  mag (magnitudes in the AB system), but undetected on 2016-08-17. The implied absolute magnitudes (, ) and red color ( mag) at days after explosion are consistent with warm dust formation in other Type II SNe (e.g., Prieto et al., 2012; Szalai et al., 2018). We used a black body function with dust emissivity to fit the mid-infrared SED and found  L,  K, and  AU. Using Equation 1 in Prieto et al. (2009) (from Dwek et al. 1983), we obtain a total dust mass of  M. These estimates are consistent with models of newly formed dust in the SN ejecta and observations for other Type II SNe at a similar epoch after explosion (e.g., Sarangi & Cherchneff, 2015).

4.6 Progenitor Metallicity

The pseudo-equivalent widths (pEWs) of the Fe II  line during the photospheric epoch of Type IIP SNe are a proxy for the progenitor metallicity. This has been shown both by models (D14) and correlation of observed pWEs with metallicities of H II regions in the host galaxies (Anderson et al., 2016). We measured the pEW evolution during the plateau phase of ASSASN-14jb and compare it with those of D13 models. The result is shown in Figure 18. The models span four progenitor metallicities (2 , 1 , 0.4 and 0.1 ) and four different mixing-length scales for the Solar composition model, which greatly change the progenitor radius. We see that the observed pEWs of the Fe II line are most consistent with the curve for 0.4 , but the Solar metallicity model for the progenitor with largest radius is also close enough. If we linearly interpolate the models and, assume that the dispersion given by different radius at Solar metallicity is a fair estimate of the uncertainty, we obtain Z = for ASASSN-14jb.

Figure 18: Measurements of the pseudo-equivalent width of the Fe II line for ASASSN-14jb (black stars), with the correspondent epoch since explosion (MJD). The solid lines show Dessart et al. (2013) models for progenitors with different metallicities and the dashed line model progenitors with Solar metallicities and different initial radii.

5 The Host Galaxy of ASASSN-14jb

Statistical constrains on properties of the progenitor stars can be also establish through the study of the galactic environment of the SNe (e.g., Stoll et al., 2013; Anderson et al., 2015; Galbany et al., 2016c). Relevant to our discussion are recent integral field spectroscopy observations (IFS) which have provided metallicities and ages of nearby stellar populations (Kuncarayakti et al., 2013; Galbany et al., 2016b, 2018; Kuncarayakti et al., 2018).

Of particular interest are the studies of H II regions, which rapidly fade away after the massive stars which provide the ionizing photons explode as SNe (e.g., Pols et al., 1998). Since typical lifetimes are up to a few tens of Myr they trace well the ongoing star formation. Typical sizes of giant H II regions are in the hundreds hundreds parsecs (Sánchez et al., 2012). Our MUSE observations of ASASSN-14jb were done with a seeing of 108 in the -band, providing a resolution of pc at the distance of the host.

We built 2D emission maps of ESO 467-G051 from the MUSE datacube following the procedures of Galbany et al. 2016b. In brief, we fitted the spectrum of each spaxel using a combination of the stellar continuum modeled with the stellar population synthesis code STARLIGHT (Cid Fernandes et al., 2009), with Gaussian profiles at the central wavelengths of the emission lines detected with SNR greater than 3, and we corrected the line ratios using the Galactic reddening map of Schlafly & Finkbeiner (2011) and an estimate of the internal extinction in ESO 467-G051 using the Balmer decrement. In addition, we used the H II-Explorer code (Sánchez et al., 2012) to detect and analyze H II regions in the datacube.

5.1 Overall Properties of ESO 467-G051

The host galaxy of ASASSN-14jb, ESO 467-G051, is an edge-on, Scd galaxy (de Vaucouleurs et al., 1991). Our observations, together with extensive data in the literature, make it possible to measure integrated fluxes from the UV up to the IR.

To do so, we (1) used co-added LCGOTN images with good seeing and measured the total magnitudes using a large elliptical aperture in ds9; (2) used a near-infrared -band image of the field from the ESO data archive777Image obtained with SOFI at the NTT by the PESSTO program (ID 191.D-0935) with SOFI mounted on the ESO/NTT to measure the magnitude of the host within an elliptical aperture calibrated to 2MASS (Cutri et al., 2003) stars in the field; (3) obtained archival ultraviolet magnitudes from GALEX (Morrissey et al., 2007) and mid-infrared magnitudes from Spitzer’s S4G survey (Sheth et al., 2010). The integrated apparent magnitudes are presented in Table 6. Using our distance, the absolute magnitudes of ESO 467-G051, corrected by Galactic extinction, are  mag,  mag, and  mag.

Figure 19 shows a histogram with the distribution of host galaxy magnitudes of Type II SNe discovered by ASAS-SN in 2013-2017 from Holoien et al. (2017a, b, c, 2018). The host galaxy of ASASSN-14jb is in the lower absolute magnitude (lower stellar mass) range of the Type II SNe discovered by ASAS-SN.

The stellar mass, age and SFR of ESO 467-G051 can be estimated from these data. We fitted stellar population synthesis (SPS) models to the , , , , 3.6 m, and 4.5 m magnitudes using FAST (Kriek et al., 2009). We assumed a Galactic extinction law, a Salpeter IMF, an exponentially declining star formation law (, with  Gyr), and the Bruzual & Charlot (2003) SPS models. The best-fit model provides the following parameters:  M,  Gyr, and  M/yr, where the uncertainties are . The SFR can be also estimated from the H extinction-corrected fluxes (Kennicutt, 1998, see below,), which provides  M/yr. This is fully consistent with the SED SPS fit.

Neutral hydrogen (HI 21 cm) observations suggest that the host galaxy of ASASSN-14jb is undergoing a direct encounter with 467-ESOG050 or NGC 7259, the host galaxy of SN 2009ip (km/s, Nordgren et al., 1997). The total HI flux is  Jy km/s and the derived HI mass is 888 M, D = 23.5 Mpc. Using our distance estimate ( Mpc) the HI mass scales to . This gives a gas fraction of .

The maximal HI rotational velocity is km/s (Hyperleda database999 In our H velocity map we observe this velocity limit more clearly on the East side of the galaxy.

Figure 19: Histogram with the distribution of absolute -band magnitudes of the host galaxies of 94 Type II SNe discovered by ASAS-SN in 2013-2017 (Holoien et al. 2017a, b, c, 2018). The dashed vertical line marks the absolute magnitude of the host galaxy of ASASSN-14jb, ESO 467-G051, indicating a relatively low stellar mass for a Type II SN in this sample.

5.2 H and Oxygen Abundance Maps

Our map of the H dereddened emission line fluxes is given in Figure 20. It includes H II regions, diffuse emission in the host, and the nebular emission from ASASSN-14jb. The position of the H II regions detected with the H II-Explorer code, and that of the SN, are marked. ASASSN-14jb is at a vertical distance of  kpc from the disk major axis. The H II region nearest to the SN is at of  kpc.

The spectra of the H II regions can be used to estimate the gas-phase oxygen abundance. We did so using the O3N2 and N2 strong emission line methods (Marino et al., 2013). The median obtained, with confidence intervals, were and , respectively. A new strong line calibration by Dopita et al. 2016, hereafter D16, which uses [S II], [N II], and H lines, provides significantly lower values, with a median of . The D16 and O3N2 maps are presented in Figure 21. Median values of the maps are and , for the spaxels with S/N .

The distributions of oxygen abundances obtained with the N2, O3N2 and D16 nebular emission line diagnostics is given in Figure 22. All the mean values are sub-Solar (using as reference , from Asplund et al. 2009), and significantly below (Galbany et al., 2018), the average for H II regions near Type II SNe of the the PMAS/PPAK Integral-Field Supernova Host Compilation (PISCO). Our host galaxy -band absolute magnitude, of  mag, compares well with the mean -band magnitude of  mag in the recent sample of Type II SNe in low-metallicity, dwarf galaxies (Gutiérrez et al., 2018). Also, as showed in Gutiérrez et al. (2018), lower luminosity galaxies hosts lower metallicity Type II SNe, filling the space between the 0.1 and 0.4 Z metallicity models. This is expected given the strong correlation between stellar mass (or galaxy luminosity) and metallicity. It is worthwhile to mention that the O3N2 map shows lower oxygen abundances in the cores of the H II regions. This could be a consequence of the unaccounted internal gradient of the ionization parameter of each H II region in the empirical strong line methods (Krühler et al., 2017).

Figure 20: H MUSE map of the host galaxy and the supernova explosion site. Black circles marks the position of the H II regions detected with H II-Explorer. Coordinates are the right ascension and declination offsets from the center of ESO 467-G051, in arcminutes. The 1 kpc scale of the image is marked in the upper left.
Figure 21: D16 (top) and O3N2 (bottom) oxygen abundance maps of the host galaxy ESO467-G051 as seen by MUSE, shifting the D16 map by +0.4 dex to provide an easier comparison. Coordinates are the right ascension and declination offsets from the center of ESO 467-G051, in arc-minutes. The 1 kpc scale of the image is marked in the upper left.

5.3 D16 and the N/O ratio

The significant discrepancy of the D16 with other strong line oxygen abundance calibrations calls our attention.This abundance scale relies on a specific relation between the N/O and O/H ratios (Dopita et al., 2016) for which it is useful to consider a diagnostic independent on any a priori relation. Pérez-Montero (2014) developed a semi-empirical code, H II-CHI-MISTRY, that derives O/H, N/O and the ionization parameter U, using the reddening-corrected fluxes of [O II] 3727 , [O III] 4363, 5007 , [N II] 6584 and [S II] 6717,6731, relative-to- H. When the O auroral lines are absent, as is our case, the code uses a limited grid of empirically constrained models to provide abundances that are consistent with the “direct” electron temperature measurement method. We applied H II-CHI-MISTRY v3.0 with the accessible lines [O I] 5007 , [N II] 6584 and [S II] 6717,6731. The output N/O ratio is shown as a function of the oxygen abundances in Figure 23. We observe a clear offset from the D16 calibration. The mean of the oxygen abundances is while the N/O ratio has a mean of . There is a trend with the ionization parameter so that higher oxygen abundances show lower N/O and U values (same with H luminosity instead of U) and, vice versa, the lower abundances correspond to higher N/O ratios and also they are closer to the D16 prescription. The negative correlation between N/O and O/H, together with the trend with H luminosity, has been observed in blue compact dwarf galaxies (Kumari et al., 2018).

Figure 22: Distribution of the gas-phase oxygen abundances used for ESO467-G051 H II regions, based on the D16, O3N2 and N2 calibrations. The dashed lines represent the median for each distribution and the black vertical line indicates the average abundance for Type II SNe hosts from the PISCO compilation using the D16 calibration (Galbany et al., 2018).
Figure 23: Oxygen abundances versus N/O ratios obtained with the code developed in Pérez-Montero 2014. Each point has the errorbars delivered by the code and each point is colored according to the ionization parameter U delivered. The solid orange line indicates the relation as assumed in the D16 calibration and the blue dots correspond to the observed values for our H II regions abundances using D16. The dashed lines indicate the asymptotes for orange line, at the very high and very low O/H ratios regimes.

6 Discussion

ASASSN-14jb is a normal Type IIP SN located at a projected distance of  kpc from the mid-plane and  kpc from the closest, diffuse H II region of ESO 467-G051. Since Type II SNe come from stars with short life times, the unusual location begs for an explanation. Several hypothesis have been advanced for SNe with no clear parent galaxy. One is that the host is an undetected low-surface brightness (LSB) galaxy, or even a small fragment of a disrupted satellite galaxy (Abadi et al., 2009). Other interesting possibility is that the progenitor of ASASSN-14jb traveled  kpc from its birthplace driven by gravitational interaction.

6.1 Type II SNe from Runaway Stars

Different scenarios have been proposed for SNe coming from runaway stars. One is the Binary Supernova Scenario (BSS) where they are pushed by the interaction with a massive binary companion that exploded as a supernova. Other is the Dynamical Ejection Scenario (DES), where the impulse comes from close encounters in clusters, including those with binaries, very massive stars, or intermediate or massive black-holes (Gvaramadze et al., 2009; Hansen & Milosavljević, 2003; Gualandris & Portegies Zwart, 2007). Yet another is the hyper-velocity star (HVS) scenario, resulting from the so-called Hills Mechanism, where velocities higher than 100 km/s are provided by interaction with a supermassive black-hole at the center of the galaxy (see Brown, 2015, for a review on HVS).

The approximate lifetime of a star with ZAMS mass of in single stellar evolution is  Myr. Covering the  kpc from the nearest star forming region in this time would require a projected space velocity of  km/s. Considering the lifetime of the least massive stars that undergo core collapse in the models of Zapartas et al. (2017),  Myr for  M progenitors, reduces the requirement to  km/s. Both are plausible velocities given the observed population of runaway stars in the Galaxy and the Magellanic Clouds (e.g., Schnurr et al., 2008). The velocity distribution of early type runaway stars in the Milky Way also supports these estimates (Silva & Napiwotzki, 2011). OB stars at high Galactic latitudes exceed the minimum velocity needed to reach 1 kpc off the Galactic disk. Scaling the result to ESO 467-G051 to account for its smaller mass, the maximum height above the disk that would be reached for a given velocity would be comparable to that of ASASSN-14jb. Let us now consider an illustrative case of the BSS scenario, with a primary star of 25 at ZAMS. The primary would explode Myr after formation (e.g., Pols et al., 1998, assuming half solar metallicity in their grid), leaving  Myr as the flight time for a secondary of 10 at ZAMS. Traveling the  kpc in this time requires a velocity of  km/s. Finally, the HVS scenario is also feasible. The projected distance from the position of ASASSN-14jb to the host galaxy center is  kpc. Taking the lifetime of the longest living CCSN progenitors of Zapartas et al. (2017)(48) the lower limit for the space velocity is  km/s. Comparing this with the case of the Type II SN 2006bx (Zinn et al., 2011), the progenitor of which would have been an HVS with  km/s, an HVS progenitor for ASASSN-14jb seems unnecessary.

On the other hand, theory is not as supportive of the runaway progenitor scenario. Eldridge et al. (2011) studied the distribution of spatial velocities and distances traveled for progenitors of CCSNe in the BSS scenario, ignoring the galaxy potential. They obtained a total average of  pc (non-Gaussian distribution) with a  km/s space velocity for the secondary of the binary system exploding as a Type IIP like SN at Solar metallicity (), while gives mildly higher velocities. A more recent study by Renzo et al. (2018) confirms these results, showing that of binary systems have a disruption velocity  km/s after the explosion of the primary. Hence, the bulk of the ejected populations does not reach the velocities needed by the progenitor of ASASSN-14jb and taking the host galaxy potential into account would make the scenario less probable. There are, however, tails in the ejected velocity distributions where the velocities are consistent with the requirements of the ASSASN-14bj progenitor.

6.2 Disk Thickness and Extraplanar Star Formation

Hakobyan et al. (2017) studied the vertical distribution of both SNe Ia and CCSNe in edge-on, late-type disk galaxies using a sample of 102 historical SNe. For Sb-Sc type galaxies the scale height, in units of the isophotal radius at a surface brightness in -band of 25 mag/arcsec or , for CCSN. ASASSN-14jb has , more than above this mean. We estimated the height scale of for ESO 467-G051 by fitting a a profile to the H profile along the perpendicular between the disk and the SN. Fits to the or the -band images gave similar results.

The scale height of each supernova type is expected to match that of the stellar population supplying their progenitors (e.g., Kangas et al., 2017). For Type II SNe we expect that the vertical distribution follows the distribution of recent star formation (OB stars). In particular, as the thick disk population is generally significantly older (Dalcanton & Bernstein, 2002; Yoachim & Dalcanton, 2008; Elmegreen et al., 2017), we are less likely to observe a Type II SN in the thick disk. Nevertheless, Howk et al. (2018b) compiled a sample of 6 extraplanar H II regions detected in nearby, edge-on disk galaxies. The projected vertical distances of these regions above their disk cover the range  kpc, up to in units of , which is within a factor of two of the height of ASASSN-14jb in the disk of ESO 467-G051. Given this, the scenario of in-situ star formation for the progenitor of ASASSN-14jb is possible (see Howk et al. 2018b for a discussion on the necessary conditions for extraplanar star formation).

The MUSE nebular spectrum can be used to put constrains on the in-situ star formation at the position of ASASSN-14jb. We fit the spectrum at the wavelength of the [O III] line and obtain a 3 upper limit to the emission line luminosity of  erg/s. We did not use H because the strong and broad SN nebular emission line dominates at those wavelengths. Comparing the [O III] luminosity upper limit with the ionization models of Xiao et al. (2018), we conclude that any underlying H II region would have to be at least  Myr old, implying M M for a single star progenitor. All massive O stars would have been gone and exploded as CCSNe in such conditions. This indicates the need for a relatively low-mass progenitor of ASASSN-14jb in the in-situ star formation scenario, as well.

6.3 Abundance Offset

Another relevant piece of information is the expected abundance offset, , of a former, or undetected, star forming region at the position of the SN (Howk et al., 2018a, b). To do so we need estimates of the oxygen abundance at the supernova position and the disk, given in the same scale. Howk et al. (2018a) show that relative abundances in the O3N2 scale are reliable so we chose this one.

For the reference disk () we take the average of the H II regions within 1 of their vertical distribution, which corresponds to a vertical distance of kpc. The reult is  dex. For ASASSN-14jb we obtain the oxygen abundance from our estimate of the iron abundance derived in Section 4.6 of . We adopt the linear relation between [O/H] and [Fe/H] given by Stoll et al. (2013), and obtain  dex. The oxygen abundance offset implied results  dex.

Recent data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey (Majewski et al 2018), provides for an independent estimate. Rojas-Arriagada et al. (2018, private comm., submitted) study abundance gradients using red giant branch stars, assumed to be likely bulge members by their spatial location. Adopting the same strategy on APOGEE data of our host we obtain relation between and similar to that of Stoll et al. (2013). A linear fit of the form , using a re-sampling technique to estimate the errors, provides coefficients and was obtained. This finally translates to an offset of  dex.

Another possible estimate of the the abundance at the supernova height is the extrapolation of the vertical gradient derived from the H II regions in the disk. Doing a linear fit as a function , we obtain a gradient of dex per height normalized by , which for ASASSN-14jb height corresponds to an oxygen abundance offset of  dex.

In Figure 24 we show the vertical oxygen abundance gradient derived for the extraplanar H II region compilation of Howk et al. (2018a) and ASASSN-14jb, as a function of the height in kiloparsec in units of . There are three interesting points to note from these plots: 1) There is a 2 offset between our abundance offset estimates (direct EW measurements from the SN versus vertical gradient of the disk abundance profile); 2) The metallicity estimate from ASASSN-14jb corresponds to a positive offset relative to the mean reference measurement of H II regions in the galaxy disk; and 3) Including the normalization, the hypothetical H II region where the progenitor of ASASSN-14jb could have formed (in an in-situ star-formation scenario) is farther than the Howk et al. (2018a) sample.

An important caveat to consider in the previous analysis is that, because the galaxy is edge-on, we are probably observing mostly the outer H II regions of the host galaxy that generally will not be representative of the whole disk. In general star-forming galaxies have higher metallicity gas in their inner regions, having a common radial gradient of dex/ (or dex) up to 2 (e.g., Sánchez et al., 2012; Pilyugin et al., 2014; Sánchez-Menguiano et al., 2016). As we do not actually know the radial coordinate of ASASSN-14jb and the observed H II regions, there might be a systematic error on what we consider the reference abundance for the disk, being higher if the observed H II regions are in average metal poor at the external disk. If we consider up to dex in the metallicity offset, given by the possible bias in our H II regions, ASASSN-14jb metallicity is still consistent with the disk at the level. This is consistent with the possibility that the progenitor star actually came from the disk as a runaway or that the gas from which it formed at this height has significant contribution of enriched material from the disk.

Figure 24: Metalicity gradient obtained from the extraplanar H II-regions from Howk et al. (2018b) and from ASASSN-14jb, as a function of the vertical distance to the disk normalized by the isophotal radius . ASASSN-14jb metallicity offset was obtained from the metallicity measurements from the FeII line pEW and the calculated mean oxygen abundance of the host, ESO467-G051, from the O3N2 diagnostic. A metallicity of is used for ASASSN-14jb as obtained interpolating the models pEWs from Dessart et al. (2013) and this was converted to an oxygen abundance using a linear scaling between [O/H] and [Fe/H] as done in Stoll et al. (2013) (marked with the unfilled square) while the calibration from Rojas, (Submitted) was used in the value marked as a black and filled square. We also estimated the oxygen abundance at the supernova position extrapolating the vertical gradient obtained from our H II-region abundances using the regions at any azimuthal position, dubbed “linear”, and from a two dimensional model, dubbed “planar”, which are marked as a black triangle and square, respectively.

6.3.1 SN-Host connection

We have constrained the properties of ASASSN-14jb both from the observations of the transient and from the observations of its host. Both sets of results provide a consistent picture.

From the early and late light curves we estimated an ejected mass ( M), given a low explosion energy ( foe) consistent with the transient energetics, and the low estimate of the Ni mass. Also, our analysis of the nebular spectra points to a low mass progenitor of  M.

The analysis of the environment of ASASSN-14jb, led us to two possible explanations for the unusual location of the transient, which are the extraplanar in-situ star formation and the runaway scenarios. In both cases, a low mass progenitor is required because their longer lifetime allows for a longer flight time in the runaway case and a longer time delay between the explosion of the more massive ionizing stars versus the less massive ones, given the strong upper limit of the possible underlying H II region in the in-situ scenario.

We have also obtained a constraint on the progenitor metallicity from a direct measurement of the Fe II pseudo equivalent width. The value of is consistent, or mildly higher, than the values inferred for the gas-phase metallicity of the disk. This, again we argue, is consistent with both scenarios as the progenitor metallicity could reflect that the star formed in the disk and was ejected or that the enriched gas itself was ejected and contributed to the abundance of the in-situ extraplanar region where the progenitor could have been born. The efficiency of extraplanar enrichment depends on the depth of the galaxy potential well, which is a decreasing function of the stellar mass (Peeples & Shankar, 2011; Brook et al., 2014; Christensen et al., 2016; Chisholm et al., 2018). Hence, the extraplanar relative enrichment is expected to be larger in low mass galaxies, as the ASASSN-14jb host, than, for example, Milky Way-like galaxies. The same reasoning applies for runaways stars as the maximum height that they achieve would be a decreasing function of the potential well depth, as well.

7 Conclusions

ASASSN-14jb shows light curves and spectral evolution consistent with typical Type II SNe. The SN had a peak absolute -band magnitude of  mag and an estimated bolometric luminosity at 50 days of , which put it in the low luminosity tail of the brightness distribution. Both the decline slope in the plateau phase ( mag per 100 days) and the expansion velocity at 50 days ( km/s) follow the observed correlations with absolute magnitudes.

The color evolution of ASASSN-14jb points to a bluer continuum than other Type IIP like SNe. This behavior may be due to a true difference in the temperature evolution or it may be an effect of lower line-blanketing in the blue part of the spectrum due to a lower metallicity of the progenitor star.

The spectral evolution presents a couple of interesting peculiarities. All our plateau observations show weaker lines (lower pseudo-EW) compared with a set of prototypical Type II SNe from the literature, which also points to a relatively bluer continuum or lower metallicity progenitor. Also, in photospheric spectra obtained at 64 and 80 days after explosion, we clearly detect high velocity absorption components at km/s in the H and H lines. These absorption components cannot be Si II  or Ba absorption features, and signal a moderate interaction with circumstellar material.

We used three empirical methods to measure Type II SNe distances and found internal consistency in the distance estimates, considering their uncertainties. The weighted average distance modulus of ASASSN-14jb was found to be  mag ( Mpc).

From the late-time, nebular-phase photometry we estimate a Ni mass of , slighly lower than the median for normal Type II SNe. Based on the near-UV photometry from Swift/UVOT and analytic models from RW11 we estimate a progenitor radius of , which is consistent with the color evolution of the mildly sub-Solar and compact RSG progenitor model (m15z8m3) presented in D13. Comparing our nebular-phase spectrum with models from D13 and J14 we conclude ASSASN14jb had a low to moderate mass progenitor (). The early light curve analysis yields  M, while at nebular times it implies  M. Both constraints intersect at the pair () (6 M, 0.25 foe). The estimated pseudo-equivalent width of the Fe II line and its time evolution are more consistent with the theoretical expectations for a metallicity progenitor. Given that the progenitor seems to be relatively compact ( ), the metallicity could be lower. More models would be needed to constrain this further.

We used the MUSE data to constrain the host galaxy gas-phase oxygen abundances. From the O3N2 and N2 strong line methods (Marino et al., 2013) the H II regions have a median abundance of and , respectively. Using the D16 (Dopita et al., 2016) calibration we obtained a median of . These values are significantly below the average of the PISCO sample mean for Type II SNe nearby H II regions, of  dex. This is independent of the diagnostics used.

Finally, we discussed scenarios for the unusual explosion site of ASASSN-14jb. We conclude that, although the probablity of ejection from the disk is low from theory (Renzo et al., 2018) and the specific mechanism to initiate the star formation at this heights is uncertain, both scenarios require a low mass progenitor with a ZAMS mass of . This estimate is consistent with the physical parameters derived for the explosion.

NM acknowledges the Insitute of Astrophysics at PUC, the Millennium Institute of Astrophysics (MAS), and the Astronomy Nucleus for University Diego Portales (UDP) for supporting this research. We thank Luc Dessart, Enrique Perez-Montero, Ondrej Pejcha, Franz Bauer, Melina Bersten, Lin Xiao, and Antonia Bevan for valuable discussions. Support for NM and JLP was provided in part by FONDECYT through the grant 1151445. Support for NM, JLP, and AC was also provided by the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. The authors thank Las Cumbres Observatory and its staff for their continued support of ASAS-SN. ASAS-SN is supported by the Gordon and Betty Moore Foundation through grant GBMF5490 to the Ohio State University and NSF grant AST-1814440. Development of ASAS-SN has been supported by NSF grant AST-0908816, the Center for Cosmology and AstroParticle Physics at the Ohio State University, the Mt. Cuba Astronomical Foundation, the Chinese Academy of Sciences South America Center for Astronomy (CASSACA), and by George Skestos. This research has made use of data from the Public ESO Spectroscopic Survey of Transient Objects (PESSTO; Smartt et al., 2015, ESO program ID 191.D-0935) and Spitzer IRAC data from 2015-09-13 (program ID 11053), 2016-08-17 (program ID 12099) and 2014-09-05 (program ID 10139). This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, Chile. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics. We acknowledge the use of the HyperLeda database ( This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Observations made with the NASA Galaxy Evolution Explorer (GALEX) were used in the analyses presented in this manuscript. Some of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST).


  • Abadi et al. (2009) Abadi, M. G., Navarro, J. F., & Steinmetz, M. 2009, ApJ, 691, L63
  • Alard & Lupton (1998) Alard, C., & Lupton, R. H. 1998, ApJ, 503, 325
  • Alard (2000) Alard, C. 2000, A&AS, 144, 363
  • Alatalo et al. (2016) Alatalo, K., Cales, S. L., Rich, J. A., et al. 2016, ApJS, 224, 38
  • Anderson & James (2009) Anderson, J. P., & James, P. A. 2009, MNRAS, 399, 559
  • Anderson et al. (2010) Anderson, J. P., Covarrubias, R. A., James, P. A., Hamuy, M., & Habergham, S. M. 2010, MNRAS, 407, 2660
  • Anderson et al. (2012) Anderson, J. P., Habergham, S. M., James, P. A., & Hamuy, M. 2012, MNRAS, 424, 1372
  • Anderson et al. (2014) Anderson, J. P., Dessart, L., Gutierrez, C. P., et al. 2014, MNRAS, 441, 671
  • Anderson et al. (2014) Anderson, J. P., González-Gaitán, S., Hamuy, M., et al. 2014, ApJ, 786, 67
  • Anderson et al. (2015) Anderson, J.P., James, P.A., Habergham, S.M., Galbany, L., & Kuncarayakti, H. 2015, PASA, 32, e019
  • Anderson et al. (2016) Anderson, J. P., Gutiérrez, C. P., Dessart, L., et al. 2016, A&A, 589, A110
  • Arcavi et al. (2010) Arcavi, I., Gal-Yam, A., Kasliwal, M. M., et al. 2010, ApJ, 721, 777
  • Arnett (1996) Arnett, D. 1996, Supernovae and Nucleosynthesis: An Investigation of the History of Matter, from the Big Bang to the Present, by D. Arnett. Princeton: Princeton University Press, 1996.,
  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
  • Bacon et al. (2010) Bacon, R., Accardo, M., Adjali, L., et al. 2010, Proc. SPIE, 7735, 773508
  • Baldwin et al. (1981) Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5
  • Barbon et al. (1979) Barbon, R., Ciatti, F., & Rosino, L. 1979, A&A, 72, 287
  • Baron et al. (2004) Baron, E., Nugent, P. E., Branch, D., & Hauschildt, P. H. 2004, ApJ, 616, L91
  • Baron et al. (2005) Baron, E., Nugent, P. E., Branch, D., & Hauschildt, P. H. 2005, 1604-2004: Supernovae as Cosmological Lighthouses, 342, 351
  • Bartunov et al. (1994) Bartunov, O. S., Tsvetkov, D. Y., &
  • Becker (2015) Becker, A. 2015, Astrophysics Source Code Library, ascl:1504.004 Filimonova, I. V. 1994, PASP, 106, 1276
  • Bersten & Hamuy (2009) Bersten, M. C., & Hamuy, M. 2009, ApJ, 701, 200
  • Bersten et al. (2011) Bersten, M. C., Benvenuto, O., & Hamuy, M. 2011, ApJ, 729, 61
  • Bertin & Arnouts (1996) Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393
  • Blaauw (1993) Blaauw, A. 1993, Massive Stars: Their Lives in the Interstellar Medium, 35, 207
  • Blinnikov & Bartunov (1993) Blinnikov, S. I., & Bartunov, O. S. 1993, A&A, 273, 106
  • Bose et al. (2013) Bose, S., Kumar, B., Sutaria, F., et al. 2013, MNRAS, 433, 1871
  • Bose et al. (2015) Bose, S., Valenti, S., Misra, K., et al. 2015, MNRAS, 450, 2373
  • Bose et al. (2016) Bose, S., Kumar, B., Misra, K., et al. 2016, MNRAS, 455, 2712
  • Bose et al. (2018) Bose, S., Dong, S., Kochanek, C. S., et al. 2018, arXiv:1804.00025
  • Brimacombe et al. (2014) Brimacombe, J., Kiyota, S., Holoien, T. W.-S., et al. 2014, The Astronomer’s Telegram, 6592,
  • Brook et al. (2014) Brook, C. B., Stinson, G., Gibson, B. K., et al. 2014, MNRAS, 443, 3809
  • Brown et al. (2013) Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031
  • Branch & Wheeler (2017) Branch, D., & Wheeler, J. C. 2017, Supernova Explosions: Astronomy and Astrophysics Library, ISBN 978-3-662-55052-6. Springer-Verlag GmbH Germany, Cap. 1.3
  • Brown et al. (2014) Brown, P. J., Breeveld, A. A., Holland, S., Kuin, P., & Pritchard, T. 2014, Ap&SS, 354, 89
  • Brown (2015) Brown, W. R. 2015, ARA&A, 53, 15
  • Bruzual & Charlot (2003) Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000
  • Cardelli et al. (1989) Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245
  • Challis (2014) Challis, P. 2014, The Astronomer’s Telegram, 6600,
  • Christensen et al. (2016) Christensen, C. R., Davé, R., Governato, F., et al. 2016, ApJ, 824, 57
  • Chevalier & Soderberg (2010) Chevalier, R. A., & Soderberg, A. M. 2010, ApJ, 711, L40
  • Chisholm et al. (2018) Chisholm, J., Tremonti, C., & Leitherer, C. 2018, MNRAS,
  • Chugai (1994) Chugai, N. N. 1994, Circumstellar Media in Late Stages of Stellar Evolution, 148
  • Chugai et al. (2005) Chugai, N. N., Fabrika, S. N., Sholukhova, O. N., et al. 2005, Astronomy Letters, 31, 792
  • Cid Fernandes et al. (2009) Cid Fernandes, R., Schoenell, W., Gomes, J. M., et al. 2009, Revista Mexicana de Astronomia y Astrofisica Conference Series, 35, 127
  • Clocchiatti & Wheeler (1997) Clocchiatti, A., & Wheeler, J. C. 1997, ApJ, 491, 375
  • Cutri et al. (2003) Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog, 2246
  • Dalcanton & Bernstein (2002) Dalcanton, J. J., & Bernstein, R. A. 2002, AJ, 124, 1328
  • Dall’Ora et al. (2014) Dall’Ora, M., Botticella, M. T., Pumo, M. L., et al. 2014, ApJ, 787, 139
  • Dessart & Hillier (2005) Dessart, L., & Hillier, D. J. 2005, A&A, 439, 671
  • Dessart & Hillier (2005) Dessart, L., & Hillier, D. J. 2005, A&A, 437, 667
  • Dessart & Hillier (2005) Dessart, L., & Hillier, D. J. 2005, The Fate of the Most Massive Stars, 332, 427
  • Dessart & Hillier (2010) Dessart, L., & Hillier, D. J. 2010, MNRAS, 405, 2141
  • Dessart et al. (2013) Dessart, L., Hillier, D. J., Waldman, R., & Livne, E. 2013, MNRAS, 433, 1745 (D13)
  • de Vaucouleurs et al. (1991) de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., et al. 1991, Third Reference Catalogue of Bright Galaxies (New York: Springer)
  • Dopita et al. (1984) Dopita, M. A., Evans, R., Cohen, M., & Schwartz, R. D. 1984, ApJ, 287, L69
  • Dopita et al. (2016) Dopita, M. A., Kewley, L. J., Sutherland, R. S., & Nicholls, D. C. 2016, Ap&SS, 361, 61
  • Dressler et al. (2011) Dressler, A., Bigelow, B., Hare, T., et al. 2011, PASP, 123, 288
  • Dwek et al. (1983) Dwek, E., A’Hearn, M. F., Becklin, E. E., et al. 1983, ApJ, 274, 168
  • Eastman et al. (1996) Eastman, R. G., Schmidt, B. P., & Kirshner, R. 1996, ApJ, 466, 911
  • Ekström et al. (2012) Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146
  • Eldridge et al. (2011) Eldridge, J. J., Langer, N., & Tout, C. A. 2011, MNRAS, 414, 3501
  • Eldridge et al. (2013) Eldridge, J. J., Fraser, M., Smartt, S. J., Maund, J. R., & Crockett, R. M. 2013, MNRAS, 436, 774
  • Eldridge et al. (2017) Eldridge, J. J., Stanway, E. R., Xiao, L., et al. 2017, PASA, 34, e058
  • Elmegreen et al. (2017) Elmegreen, B. G., Elmegreen, D. M., Tompkins, B., & Jenks, L. G. 2017, ApJ, 847, 14
  • Elmhamdi et al. (2003) Elmhamdi, A., Danziger, I. J., Chugai, N., et al. 2003, MNRAS, 338, 939
  • Fabricant et al. (1998) Fabricant, D., Cheimets, P., Caldwell, N., & Geary, J. 1998, PASP, 110, 79
  • Faran et al. (2014) Faran, T., Poznanski, D., Filippenko, A. V., et al. 2014, MNRAS, 442, 844
  • Fazio et al. (2004) Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10
  • Filippenko (1988) Filippenko, A. V. 1988, AJ, 96, 1941
  • Filippenko et al. (1993) Filippenko, A. V., Matheson, T., & Ho, L. C. 1993, ApJ, 415, L103
  • Filippenko (1997) Filippenko, A. V. 1997, ARA&A, 35, 309
  • Folatelli et al. (2015) Folatelli, G., Bersten, M. C., Kuncarayakti, H., et al. 2015, ApJ, 811, 147
  • Folatelli et al. (2016) Folatelli, G., Van Dyk, S. D., Kuncarayakti, H., et al. 2016, ApJ, 825, L22
  • Fransson & Chevalier (1989) Fransson, C., & Chevalier, R. A. 1989, ApJ, 343, 323
  • Fraser et al. (2013) Fraser, M., Inserra, C., Jerkstrand, A., et al. 2013, MNRAS, 433, 1312
  • Galbany et al. (2014) Galbany, L., Stanishev, V., Mourão, A. M., et al. 2014, A&A, 572, A38
  • Galbany et al. (2016c) Galbany, L., Stanishev, V., Mourão, A. M., et al. 2016, A&A, 591, A48
  • Galbany et al. (2016b) Galbany, L., Anderson, J. P., Rosales-Ortega, F. F., et al. 2016, MNRAS, 455, 4087
  • Galbany et al. (2016a) Galbany, L., Hamuy, M., Phillips, M. M., et al. 2016, AJ, 151, 33
  • Galbany et al. (2018) Galbany, L., Anderson, J. P., Sánchez, S. F., et al. 2018, ApJ, 855, 107
  • Gall et al. (2015) Gall, E. E. E., Polshaw, J., Kotak, R., et al. 2015, A&A, 582, A3
  • Gezari et al. (2010) Gezari, S., Rest, A., Huber, M. E., et al. 2010, ApJ, 720, L77
  • González-Gaitán et al. (2015) González-Gaitán, S., Tominaga, N., Molina, J., et al. 2015, MNRAS, 451, 2212
  • Gualandris & Portegies Zwart (2007) Gualandris, A., & Portegies Zwart, S. 2007, MNRAS, 376, L29
  • Gutiérrez et al. (2014) Gutiérrez, C. P., Anderson, J. P., Hamuy, M., et al. 2014, ApJ, 786, L15
  • Gutiérrez et al. (2017) Gutiérrez, C. P., Anderson, J. P., Hamuy, M., et al. 2017, ApJ, 850, 89
  • Gutiérrez et al. (2017) Gutiérrez, C. P., Anderson, J. P., Hamuy, M., et al. 2017, ApJ, 850, 90
  • Gutiérrez et al. (2018) Gutiérrez, C. P., Anderson, J. P., Sullivan, M., et al. 2018, MNRAS, 479, 3232
  • Green et al. (2014) Green, A. W., Glazebrook, K., McGregor, P. J., et al. 2014, MNRAS, 437, 1070
  • Gvaramadze et al. (2009) Gvaramadze, V. V., Gualandris, A., & Portegies Zwart, S. 2009, MNRAS, 396, 570
  • Hakobyan et al. (2008) Hakobyan, A. A., Petrosian, A. R., McLean, B., et al. 2008, A&A, 488, 523
  • Hakobyan et al. (2009) Hakobyan, A. A., Mamon, G. A., Petrosian, A. R., Kunth, D., & Turatto, M. 2009, A&A, 508, 1259
  • Hakobyan et al. (2014) Hakobyan, A. A., Nazaryan, T. A., Adibekyan, V. Z., et al. 2014, Multiwavelength AGN Surveys and Studies, 304, 339
  • Hakobyan et al. (2014) Hakobyan, A. A., Nazaryan, T. A., Adibekyan, V. Z., et al. 2014, MNRAS, 444, 2428
  • Hakobyan et al. (2016) Hakobyan, A. A., Karapetyan, A. G., Barkhudaryan, L. V., et al. 2016, MNRAS, 456, 2848
  • Hakobyan et al. (2017) Hakobyan, A. A., Barkhudaryan, L. V., Karapetyan, A. G., et al. 2017, MNRAS, 471, 1390
  • Hamuy & Pinto (2002) Hamuy, M., & Pinto, P. A. 2002, ApJ, 566, L63
  • Hamuy et al. (2001) Hamuy, M., Pinto, P. A., Maza, J., et al. 2001, ApJ, 558, 615
  • Hamuy (2003) Hamuy, M. 2003, ApJ, 582, 905
  • Hansen & Milosavljević (2003) Hansen, B. M. S., & Milosavljević, M. 2003, ApJ, 593, L77
  • Heger et al. (2003) Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003, ApJ, 591, 288
  • Henden et al. (2012) Henden, A. A., Levine, S. E., Terrell, D., Smith, T. C., & Welch, D. 2012, Journal of the American Association of Variable Star Observers (JAAVSO), 40, 430
  • Henry & Worthey (1999) Henry, R. B. C., & Worthey, G. 1999, PASP, 111, 919
  • Hillebrandt & Niemeyer (2000) Hillebrandt, W., & Niemeyer, J. C. 2000, ARA&A, 38, 191
  • Holoien et al. (2017a) Holoien, T. W.-S., Stanek, K. Z., Kochanek, C. S., et al. 2017a, MNRAS, 464, 2672
  • Holoien et al. (2017b) Holoien, T. W.-S., Brown, J. S., Stanek, K. Z., et al. 2017b, MNRAS, 467, 1098
  • Holoien et al. (2017c) Holoien, T. W.-S., Brown, J. S., Stanek, K. Z., et al. 2017c, MNRAS, 471, 4966
  • Holoien et al. (2018) Holoien, T. W.-S., Brown, J. S., Vallely, P. J., et al. 2018, MNRAS, submitted, arXiv:1811.08904
  • Hoogerwerf et al. (2001) Hoogerwerf, R., de Bruijne, J. H. J., & de Zeeuw, P. T. 2001, A&A, 365, 49
  • Howk et al. (2018b) Howk, J. C., Rueff, K. M., Lehner, N., et al. 2018, ApJ, 856, 167
  • Howk et al. (2018a) Howk, J. C., Rueff, K. M., Lehner, N., et al. 2018, ApJ, 856, 166
  • Huang et al. (2016) Huang, F., Wang, X., Zampieri, L., et al. 2016, ApJ, 832, 139
  • Huang et al. (2018) Huang, F., Wang, X.-F., Hosseinzadeh, G., et al. 2018, MNRAS, 475, 3959
  • Inserra et al. (2013) Inserra, C., Pastorello, A., Turatto, M., et al. 2013, A&A, 555, A142
  • Jerkstrand et al. (2012) Jerkstrand, A., Fransson, C., Maguire, K., et al. 2012, A&A, 546, A28
  • Jerkstrand et al. (2014) Jerkstrand, A., Smartt, S. J., Fraser, M., et al. 2014, MNRAS, 439, 3694 (J14)
  • Jerkstrand (2017) Jerkstrand, A. 2017, arXiv:1702.06702
  • Jerkstrand et al. (2018) Jerkstrand, A., Ertl, T., Janka, H.-T., et al. 2018, MNRAS, 475, 277 (J18)
  • José (2016) José, J. 2016, Stellar Explosions: Hydrodynamics and Nucleosynthesis, CRC/Taylor and Francis
  • Kangas et al. (2017) Kangas, T., Portinari, L., Mattila, S., et al. 2017, A&A, 597, A92
  • Kasen & Woosley (2009) Kasen, D., & Woosley, S. E. 2009, ApJ, 703, 2205
  • Kennicutt (1998) Kennicutt, R. C., Jr. 1998, ARA&A, 36, 189
  • Kewley et al. (2001) Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121
  • Khazov et al. (2016) Khazov, D., Yaron, O., Gal-Yam, A., et al. 2016, ApJ, 818, 3
  • Kirshner & Kwan (1974) Kirshner, R. P., & Kwan, J. 1974, ApJ, 193, 27
  • Kriek et al. (2009) Kriek, M., van Dokkum, P. G., Labbé, I., et al. 2009, ApJ, 700, 221
  • Krühler et al. (2017) Krühler, T., Kuncarayakti, H., Schady, P., et al. 2017, A&A, 602, A85
  • Kuncarayakti et al. (2013) Kuncarayakti, H., Doi, M., Aldering, G., et al. 2013, AJ, 146, 31
  • Kuncarayakti et al. (2015) Kuncarayakti, H., Maeda, K., Bersten, M. C., et al. 2015, A&A, 579, A95
  • Kuncarayakti et al. (2018) Kuncarayakti, H., Anderson, J. P., Galbany, L., et al. 2018, A&A, 613, A35
  • Kumari et al. (2018) Kumari, N., James, B. L., Irwin, M. J., Amorín, R., & Pérez-Montero, E. 2018, MNRAS, 476, 3793
  • Lang et al. (2010) Lang, D., Hogg, D. W., Mierle, K., Blanton, M., & Roweis, S. 2010, AJ, 139, 1782
  • Leonard et al. (2002) Leonard, D. C., Filippenko, A. V., Li, W., et al. 2002, AJ, 124, 2490
  • Leonard et al. (2002) Leonard, D. C., Filippenko, A. V., Gates, E. L., et al. 2002, PASP, 114, 35
  • Litvinova & Nadezhin (1985) Litvinova, I. Y., & Nadezhin, D. K. 1985, Soviet Astronomy Letters, 11, 145
  • McMillan & Ciardullo (1996) McMillan, R. J., & Ciardullo, R. 1996, ApJ, 473, 707
  • Maguire et al. (2012) Maguire, K., Jerkstrand, A., Smartt, S. J., et al. 2012, MNRAS, 420, 3451
  • Marino et al. (2013) Marino, R. A., Rosales-Ortega, F. F., Sánchez, S. F., et al. 2013, A&A, 559, A114
  • Mauerhan et al. (2013) Mauerhan, J. C., Smith, N., Filippenko, A. V., et al. 2013, MNRAS, 430, 1801
  • McEvoy et al. (2017) McEvoy, C. M., Dufton, P. L., Smoker, J. V., et al. 2017, ApJ, 842, 32
  • McSwain et al. (2007) McSwain, M. V., Boyajian, T. S., Grundstrom, E. D., & Gies, D. R. 2007, ApJ, 655, 473
  • Mikhailova et al. (2007) Mikhailova, G. A., Bartunov, O. S., & Tsvetkov, D. Y. 2007, Astronomy Letters, 33, 715
  • Minkowski (1941) Minkowski, R. 1941, PASP, 53, 224
  • Modjaz et al. (2008) Modjaz, M., Kewley, L., Kirshner, R. P., et al. 2008, AJ, 135, 1136
  • Moriya et al. (2016) Moriya, T. J., Pruzhinskaya, M. V., Ergon, M., & Blinnikov, S. I. 2016, MNRAS, 455, 423
  • Morrissey et al. (2007) Morrissey, P., Conrow, T., Barlow, T. A., et al. 2007, ApJS, 173, 682
  • Mould et al. (2000) Mould, J. R., Huchra, J. P., Freedman, W. L., et al. 2000, ApJ, 529, 786
  • Müller et al. (2017) Müller, T., Prieto, J. L., Pejcha, O., & Clocchiatti, A. 2017, ApJ, 841, 127
  • Nakar & Sari (2010) Nakar, E., & Sari, R. 2010, ApJ, 725, 904
  • Nomoto et al. (1993) Nomoto, K., Suzuki, T., Shigeyama, T., et al. 1993, Nature, 364, 507
  • Nordgren et al. (1997) Nordgren, T. E., Chengalur, J. N., Salpeter, E. E., & Terzian, Y. 1997, AJ, 114, 913
  • Olivares E. et al. (2010) Olivares E., F., Hamuy, M., Pignata, G., et al. 2010, ApJ, 715, 833
  • Pastorello et al. (2009) Pastorello, A., Valenti, S., Zampieri, L., et al. 2009, MNRAS, 394, 2266
  • Pastorello et al. (2013) Pastorello, A., Cappellaro, E., Inserra, C., et al. 2013, ApJ, 767, 1
  • Peeples & Shankar (2011) Peeples, M. S., & Shankar, F. 2011, MNRAS, 417, 2962
  • Pejcha & Prieto (2015) Pejcha, O., & Prieto, J. L. 2015, ApJ, 806, 225
  • Pejcha & Prieto (2015) Pejcha, O., & Prieto, J. L. 2015, ApJ, 799, 215
  • Pérez-Montero (2014) Pérez-Montero, E. 2014, MNRAS, 441, 2663
  • Petrosian et al. (2005) Petrosian, A., Navasardyan, H., Cappellaro, E., et al. 2005, AJ, 129, 1369
  • Phillips et al. (2013) Phillips, M. M., Simon, J. D., Morrell, N., et al. 2013, ApJ, 779, 38
  • Pilyugin et al. (2014) Pilyugin, L. S., Grebel, E. K., & Kniazev, A. Y. 2014, AJ, 147, 131
  • Pols et al. (1998) Pols, O. R., Schröder, K.-P., Hurley, J. R., Tout, C. A., & Eggleton, P. P. 1998, MNRAS, 298, 525
  • Polshaw et al. (2016) Polshaw, J., Kotak, R., Dessart, L., et al. 2016, A&A, 588, A1
  • Popov (1993) Popov, D. V. 1993, ApJ, 414, 712
  • Potashov et al. (2013) Potashov, M., Blinnikov, S., Baklanov, P., & Dolgov, A. 2013, MNRAS, 431, L98
  • Prieto et al. (2008) Prieto, J. L., Stanek, K. Z., & Beacom, J. F. 2008, ApJ, 673, 999
  • Prieto et al. (2009) Prieto, J. L., Sellgren, K., Thompson, T. A., & Kochanek, C. S. 2009, ApJ, 705, 1425
  • Prieto et al. (2012) Prieto, J. L., Lee, J. C., Drake, A. J., et al. 2012, ApJ, 745, 70
  • Prieto et al. (2013) Prieto, J. L., Brimacombe, J., Drake, A. J., & Howerton, S. 2013, ApJ, 763, L27
  • Prieto et al. (2016) Prieto, J. L., Krühler, T., Anderson, J. P., et al. 2016, ApJ, 830, L32
  • Rabinak & Waxman (2011) Rabinak, I., & Waxman, E. 2011, ApJ, 728, 63
  • Renzo et al. (2018) Renzo, M., Zapartas, E., de Mink, S. E., et al. 2018, arXiv:1804.09164
  • Rodríguez et al. (2014) Rodríguez, Ó., Clocchiatti, A., & Hamuy, M. 2014, AJ, 148, 107
  • Rodríguez et al. (2014) Rodríguez, Ó., Pignata, A., et al., in preparation (R18)
  • Rubin et al. (2016) Rubin, A., Gal-Yam, A., De Cia, A., et al. 2016, ApJ, 820, 33
  • Sana et al. (2012) Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444
  • Sana (2017) Sana, H. 2017, The Lives and Death-Throes of Massive Stars, 329, 110
  • Sánchez et al. (2012) Sánchez, S. F., Rosales-Ortega, F. F., Marino, R. A., et al. 2012, A&A, 546, A2
  • Sánchez et al. (2014) Sánchez, S. F., Rosales-Ortega, F. F., Iglesias-Páramo, J., et al. 2014, A&A, 563, A49
  • Sánchez-Menguiano et al. (2016) Sánchez-Menguiano, L., Sánchez, S. F., Pérez, I., et al. 2016, A&A, 587, A70
  • Sanders et al. (2015) Sanders, N. E., Soderberg, A. M., Gezari, S., et al. 2015, ApJ, 799, 208
  • Sapir & Waxman (2017) Sapir, N., & Waxman, E. 2017, ApJ, 838, 130
  • Sarangi & Cherchneff (2015) Sarangi, A., & Cherchneff, I. 2015, A&A, 575, A95
  • Schlafly & Finkbeiner (2011) Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103
  • Schlegel (1990) Schlegel, E. M. 1990, MNRAS, 244, 269
  • Schmidt et al. (1992) Schmidt, B. P., Kirshner, R. P., & Eastman, R. G. 1992, ApJ, 395, 366
  • Schnurr et al. (2008) Schnurr, O., Moffat, A. F. J., St-Louis, N., Morrell, N. I., & Guerrero, M. A. 2008, MNRAS, 389, 806
  • Shappee et al. (2014) Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48
  • Sheth et al. (2010) Sheth, K., Regan, M., Hinz, J. L., et al. 2010, PASP, 122, 1397
  • Silva & Napiwotzki (2011) Silva, M. D. V., & Napiwotzki, R. 2011, MNRAS, 411, 2596
  • Silverman et al. (2017) Silverman, J. M., Pickett, S., Wheeler, J. C., et al. 2017, MNRAS, 467, 369
  • Smartt et al. (2009) Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 395, 1409
  • Smartt (2015) Smartt, S. J. 2015, PASA, 32, e016
  • Smartt et al. (2015) Smartt, S. J., Valenti, S., Fraser, M., et al. 2015, A&A, 579, A40
  • Smith et al. (2014) Smith, N., Mauerhan, J. C., & Prieto, J. L. 2014, MNRAS, 438, 1191
  • Smith et al. (2016) Smith, N., Andrews, J. E., & Mauerhan, J. C. 2016, MNRAS, 463, 2904
  • Stathakis & Sadler (1991) Stathakis, R. A., & Sadler, E. M. 1991, MNRAS, 250, 786
  • Stoll et al. (2013) Stoll, R., Prieto, J. L., Stanek, K. Z., & Pogge, R. W. 2013, ApJ, 773, 12
  • Sukhbold et al. (2016) Sukhbold, T., Ertl, T., Woosley, S. E., Brown, J. M., & Janka, H.-T. 2016, ApJ, 821, 38
  • Szalai et al. (2018) Szalai, T., Zsíros, S., Fox, O. D., Pejcha, O., & Müller, T. 2018, arXiv:1803.02571
  • Takáts & Vinkó (2012) Takáts, K., & Vinkó, J. 2012, MNRAS, 419, 2783
  • Takáts et al. (2014) Takáts, K., Pumo, M. L., Elias-Rosa, N., et al. 2014, MNRAS, 438, 368
  • Tartaglia et al. (2018) Tartaglia, L., Sand, D. J., Valenti, S., et al. 2018, ApJ, 853, 62
  • Terreran et al. (2016) Terreran, G., Jerkstrand, A., Benetti, S., et al. 2016, MNRAS, 462, 137
  • Turatto et al. (2007) Turatto, M., Benetti, S., & Pastorello, A. 2007, Supernova 1987A: 20 Years After: Supernovae and Gamma-Ray Bursters, 937, 187
  • Uomoto (1986) Uomoto, A. 1986, ApJ, 310, L35
  • Utrobin & Chugai (2009) Utrobin, V. P., & Chugai, N. N. 2009, A&A, 506, 829
  • Valenti et al. (2014) Valenti, S., Sand, D., Pastorello, A., et al. 2014, MNRAS, 438, L101
  • Valenti et al. (2015) Valenti, S., Sand, D., Stritzinger, M., et al. 2015, MNRAS, 448, 2608
  • Valenti et al. (2016) Valenti, S., Howell, D. A., Stritzinger, M. D., et al. 2016, MNRAS, 459, 3939
  • Van Dyk et al. (1999) Van Dyk, S. D., Peng, C. Y., Barth, A. J., & Filippenko, A. V. 1999, AJ, 118, 2331
  • van Zee & Haynes (2006) van Zee, L., & Haynes, M. P. 2006, ApJ, 636, 214
  • Weilbacher et al. (2014) Weilbacher, P. M., Streicher, O., Urrutia, T., et al. 2014, Astronomical Data Analysis Software and Systems XXIII, 485, 451
  • Xiao et al. (2018) Xiao, L., Stanway, E. R., & Eldridge, J. J. 2018, MNRAS, 477, 904
  • Yaron et al. (2017) Yaron, O., Perley, D. A., Gal-Yam, A., et al. 2017, Nature Physics, 13, 510
  • Yoachim & Dalcanton (2008) Yoachim, P., & Dalcanton, J. J. 2008, ApJ, 683, 707
  • Zapartas et al. (2017) Zapartas, E., de Mink, S. E., Izzard, R. G., et al. 2017, A&A, 601, A29
  • Zhang & Wang (2014) Zhang, J., & Wang, X. 2014, The Astronomer’s Telegram, 6601,
  • Zinn et al. (2011) Zinn, P.-C., Grunden, P., & Bomans, D. J. 2011, A&A, 536, A103
  • de Jaeger et al. (2015) de Jaeger, T., González-Gaitán, S., Anderson, J. P., et al. 2015, ApJ, 815, 121
  • de Jaeger et al. (2017) de Jaeger, T., González-Gaitán, S., Hamuy, M., et al. 2017, ApJ, 835, 166

Appendix A Photometry

Here we describe the photometric pipeline used to derive the magnitudes of ASASSN-14jb in the standard system. For each frame previously reduced we do the following steps to obtain the aperture photometry for the supernova:

  • We run SExtractor (Bertin & Arnouts, 1996) to obtain an initial estimate for the FWHM of the stars in the image, using bright stars far from the edges of the frame, and global background level.

  • We re-run SExtractor using the initial estimates and fixing the aperture to 1.6 times the median FWHM of the stars. After cross-matching this catalog with our local standard star catalog we reject the standard objects that have photometric uncertainties .

  • With the zero-point and uncertainty for each local standard star (no color term) we calculate the mean zero-point of the image with its final error including both statistical photometric uncertainties and the rms of the zero-point. The typical value of the light curves deviations were  mag.


    So with this the final photometry of the supernova and its variance are:

  • Regarding color terms, as it is showed in Valenti et al. (in prep)101010See, the only significant color term for LCOGT cameras in observations obtained in the filters are on the -band. We apply this correction using the color but with the -band magnitudes only corrected trough a mean zero-point term. That is, we follow the equations:


    with the corresponding error propagation. The -band color term correction obtained was important only at early times when the SN is bluest.

To verify our results from aperture photometry, we ran PSF fitting photometry through our pipeline using the softwares DOPHOT and DAOPHOT in IRAF for selected frames, and in both cases we obtained consistent results. We chose the SExtractor aperture photometry because it was the most robust method. The PSF photometry of DOPHOT was limited because of the fixed gaussian profile shape and DAOPHOT required too much tuning of the parameters to obtain good results on frames that had significant sky variations or were noiser than average.

Appendix B P-Cygni profile fitting

Expansion velocities were estimated using the velocity shift at maximum abpsorption in the, continuum normalized, P-Cygni profiles. For a given supernova we took each individual spectra and, after masking the regions for the strong lines H and H and NaI a black body is fitted. If the spectra is after 20 days past explosion the range  Å is also removed to avoid the contamination from line blanketing. Also a simple power law is fitted and if the residuals in this case are better this fit is adopted as the continuum. after estimating the continuum the spectra and wavelength is normalized as , where the 0 subscript refers to the line position (e.g.  Å for H).
To find a minima of a profile throughout a polinomial fitting, a function written in Python code was created. If no initial guess is given for the maximum abpsorption velocity (or equivalently ) first the data is smoothed with a high degree polinomial to find a guess. Then a new fitting window is defined, centered on the this guess, and a new low degree polinomial (from 2 up to 8) is fitted to find a new minima. In this final fit a weighting is applied to the data to avoid problems on the defined edges of the fitting window and the adopted guess. This weighting goes on the residuals as :


The e-fold scale is used to only consider the points more closed to the assumed profile minima. For the sample of this work a systematic error was considered using different values for of 10,50 and 100 Å, at least. This procedure effectively measures how biased is the chosen minima .

We estimate a global statistical error, for each , using a Jacknife procedure, where each data point is masked before fitting and a mean and standard deviation are obtained from the distribution.

Figure 25: Example of a P-Cygni profile fit. The x-axis is in units of 1000 km/s. The blue dots are the data points for the normalized spectra and the red points are the points actually used for the fitting procedure. The vertical line marks the position of the final velocity adopted and the cyan vertical indicates zero velocity at the line position.

Appendix C Distance Measurement

The methods for Type II SNe distance measurement can be categorized in theoretically based, such as the Expanding Photosphere Method (EPM, Kirshner & Kwan, 1974; Schmidt et al., 1992) or the spectra-fitting expanding atmosphere method (SEAM, Baron et al., 2004), and purely empirical such as the Standard Candle Method (SCM , Hamuy & Pinto, 2002), the Photospheric Magnitude Method (PMM Rodríguez et al., 2014, 2014, henceforth R14; R18) or the Photometric Color Method (DJ17 de Jaeger et al., 2017).

We will compute distances to ASSASN-14bj using the three empirical ones. We will not apply -corrections since the low redshift of the source makes them very small (e.g., less than 0.01 mag for the -band).

c.1 Photospheric Magnitude Method

The PMM is a method to measure distances to Type II SNe that relies on a time-variable standardization of the photospheric magnitude evolution. This method can be viewed as a generalization for various epochs of the SCM (Hamuy & Pinto 2002) in the sense that both rely on the physical connection between the expansion velocity (a measure of the explosion energy) and the plateau luminosity, but the PMM has a different approach in the calibration. For a distance to be measured we need a good constraint on the shock-breakout epoch, magnitudes corrected from both host and local extinction and corrected by the -correction, and expansion velocities inferred from the minima of the P-Cygni profiles (usually from the Fe II lines) in the photospheric phase. Explicitly we have,


where is the rest-frame and extinction corrected magnitude at the wavelength , is the photospheric magnitude, is the factor that considers the surface area of the photosphere and depends on the photospheric radius, which for homologous expansion is simply the product of photospheric velocity and the time since the beginning of the expansion (see R14 for details). Although we have named the ‘explosion epoch‘, this reference time is actually associated with the shock breakout. As our constraint on is of the order of the shock propagation time to reach a RSG surface after core-collapse (Arnett, 1996) we do not make an explicit distinction. We have measured the velocities of ASASSN-14jb from the Fe II line (see Figure 10) and the extinction corrected magnitudes in the -band without applying a -correction. Averaging the distances obtained at the three epochs with Fe II velocities, we obtain an average distance modulus of mag.

To make a direct comparison of the distance measurements of other methods we need consider that the PMM was calibrated using local distance ladder measurements, so the method carries an implicit value of the Hubble constant. As the value used by R14 is for the -band calibration, we will take this value as reference in the following analysis. We also use now the CMB frame redshift of for the ASASSN-14jb host. (Mould et al., 2000).

c.2 Standard Candle Method

The most common method to measure Type IIP like distances is the SCM which relies on the tight correlation between brightness and the expansion velocity at a given epoch in the plateau phase. Similar to the Type Ia SNe calibrations, the SCM parameters are fitted directly in the Hubble diagram following the model (Hamuy & Pinto, 2002),


where is the absolute magnitude in the -band, is the photospheric velocity, usually at 50 days, and () are the parameters to fit using a statistically significant sample. The distance modulus in the DJ17 prescription takes the form,


. We use the parameters obtained it DJ17 for the SDSS color combination111111. With the interpolated velocity in at 40 days, using a power law model, of and the color at 45 days mag, the distance modulus is mag.

c.3 Photometric Colour Method

The PCM, presented in (de Jaeger et al., 2015), exploits the correlation between the light curve decline rate at the plateau, , and the luminosity (Anderson et al., 2014). In this method the distance modulus is obtained via,


where the values of the parameters taken from DJ17 are . Measuring the slope in the -band light curve we find . With the corresponding color we obtain mag.

Star B V g r i
1 () () () () ()
2 () () () () ()
3 () () () () ()
4 () () () () ()
5 () () () () ()
6 () () () () ()
7 () () () () ()
8 () () () () ()
9 () () () () ()
10 () () () () ()
Table 1: magnitudes of the local APASS standard star sequence used.
UT Date MJD V ( error) Telescope
(days) (days) (mag)
2014-10-13 56943.098 -3.485 ASAS-SN
2014-10-19 56949.093 2.510 16.861 (0.129) ASAS-SN
2014-10-20 56950.060 3.477 16.284 (0.069) ASAS-SN
2014-10-22 56952.058 5.475 16.154 (0.065) ASAS-SN
2014-10-23 56953.051 6.468 16.072 (0.065) ASAS-SN
2014-10-29 56959.084 12.501 16.139 (0.069) ASAS-SN
2014-11-04 56965.117 18.534 16.040 (0.109) ASAS-SN
2014-11-07 56968.064 21.481 16.026 (0.112) ASAS-SN
2014-11-11 56972.059 25.476 16.071 (0.063) ASAS-SN
2014-11-12 56973.144 26.561 15.907 (0.063) ASAS-SN
2014-11-16 56977.056 30.473 16.118 (0.078) ASAS-SN
2014-11-21 56982.058 35.475 16.089 (0.068) ASAS-SN
2014-11-24 56985.046 38.463 16.106 (0.068) ASAS-SN
2014-12-01 56992.053 45.470 16.125 (0.088) ASAS-SN
2014-12-04 56995.031 48.448 16.230 (0.103) ASAS-SN
2014-12-08 56999.035 52.452 16.213 (0.084) ASAS-SN
2014-12-09 57000.034 53.451 15.995 (0.081) ASAS-SN
2014-12-14 57005.026 58.443 15.984 (0.108) ASAS-SN
2014-12-19 57010.044 63.460 16.030 (0.062) ASAS-SN
2015-04-10 57122.421 175.838 ASAS-SN
Table 2: ASAS-SN -band photometry of ASASSN-14jb, including pre-discovery and post-plateau upper limits. (*) : Days after the estimated explosion time, (MJD).

LCOGT and NTT/EFOSC photometry of ASASSN-14jb. (*) : Days after the estimated explosion time, .

UT Date MJD Telescope
2014-10-19 56949.785 3.202 16.077 (0.023) 16.402 (0.019) LCOGT 1m/SAAO
2014-10-19 56949.901 3.318 16.072 (0.020) 16.376 (0.014) LCOGT 1m/SAAO
2014-10-20 56950.393 3.810 15.986 (0.022) 16.237 (0.015) LCOGT 1m/SSO
2014-10-20 56950.802 4.219 15.960 (0.021) 16.159 (0.014) LCOGT 1m/SAAO
2014-10-21 56951.536 4.953 16.084 (0.036) 16.016 (0.020) 15.973 (0.021) 16.080 (0.021) 16.231 (0.033) LCOGT 1m/SSO
2014-10-22 56952.532 5.949 16.055 (0.025) 16.022 (0.016) 15.948 (0.021) 16.055 (0.013) 16.217 (0.021) LCOGT 1m/SSO
2014-10-23 56953.527 6.944 15.992 (0.025) 15.970 (0.017) 15.944 (0.018) 15.984 (0.012) 16.116 (0.023) LCOGT 1m/SSO
2014-10-24 56954.141 7.558 15.995 (0.027) 15.941 (0.016) 15.913 (0.022) 15.982 (0.014) 16.195 (0.027) LCOGT 1m/MDO
2014-10-24 56954.543 7.960 16.021 (0.033) 15.960 (0.016) 15.927 (0.020) 15.945 (0.011) 16.091 (0.022) LCOGT 1m/SSO
2014-10-26 56956.830 10.247 16.026 (0.023) 15.974 (0.016) 15.948 (0.012) 16.128 (0.020) LCOGT 1m/SAAO
2014-10-29 56959.546 12.963 16.097 (0.023) 15.980 (0.015) 15.967 (0.021) 15.911 (0.011) 16.084 (0.020) LCOGT 1m/SSO
2014-11-01 56962.518 15.935 16.166 (0.036) 16.020 (0.026) 16.064 (0.026) 15.910 (0.017) 16.096 (0.042) LCOGT 1m/SSO
2014-11-04 56965.576 18.993 16.211 (0.028) 15.955 (0.019) 16.023 (0.030) 15.861 (0.014) 15.991 (0.023) LCOGT 1m/SSO
2014-11-07 56968.505 21.921 16.309 (0.032) 15.984 (0.024) 16.046 (0.029) 15.860 (0.023) 15.934 (0.037) LCOGT 1m/SSO
2014-11-10 56971.505 24.921 16.348 (0.029) 15.949 (0.020) 16.127 (0.024) 15.833 (0.016) 15.948 (0.027) LCOGT 1m/SSO
2014-11-13 56974.763 28.180 16.440 (0.028) 16.018 (0.017) 16.193 (0.022) 15.866 (0.012) 15.980 (0.021) LCOGT 1m/SAAO
2014-11-17 56978.766 32.183 16.510 (0.032) 16.023 (0.019) 16.212 (0.023) 15.867 (0.014) 15.939 (0.023) LCOGT 1m/SAAO
2014-11-21 56982.777 36.194 16.566 (0.034) 15.998 (0.019) 16.223 (0.024) 15.853 (0.014) 15.861 (0.029) LCOGT 1m/SAAO
2014-11-25 56986.781 40.198 16.582 (0.026) 15.980 (0.017) 16.274 (0.022) 15.848 (0.013) 15.902 (0.021) LCOGT 1m/SAAO
2014-11-28 56989.854 43.271 16.633 (0.031) 15.972 (0.021) 16.259 (0.024) 15.839 (0.014) 15.830 (0.022) LCOGT 1m/SAAO
2014-12-01 56992.776 46.193 16.679 (0.050) 15.979 (0.025) 16.255 (0.025) 15.847 (0.023) 15.805 (0.033) LCOGT 1m/SAAO
2014-12-04 56995.836 49.253 16.660 (0.059) 15.961 (0.026) 16.259 (0.030) 15.784 (0.021) 15.751 (0.028) LCOGT 1m/SAAO
2014-12-08 56999.102 52.519 16.681 (0.029) 15.962 (0.016) 16.262 (0.021) 15.802 (0.012) 15.797 (0.019) LCOGT 1m/CTIO
2014-12-13 57004.784 58.201 16.669 (0.027) 15.955 (0.015) 16.239 (0.020) 15.785 (0.011) 15.788 (0.019) LCOGT 1m/SAAO
2014-12-16 57007.794 61.211 16.678 (0.026) 15.954 (0.016) 16.273 (0.020) 15.815 (0.011) 15.802 (0.020) LCOGT 1m/SAAO
2014-12-26 57017.040 70.457 16.750 (0.032)