Free Access
Issue
A&A
Volume 613, May 2018
Article Number A35
Number of page(s) 18
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201731923
Published online 28 May 2018

© ESO 2018

1 Introduction

A core-collapse (CC) supernova (SN) is produced when the stellar core of a massive star collapses onto itself at the end of its lifetime. Supernovae are among the brightest phenomena in the Universe. At peak brightness, an SN may rival or even outshine the entire host galaxy where it resides. The explosion distributes heavy elements forged inside the stellar interior (Hoyle & Fowler 1960; Nomoto et al. 2006), driving chemical enrichment of the interstellar medium (Matteucci & Greggio 1986; Timmes et al. 1995), and may trigger new waves of star formation (McCray & Kafatos 1987; Thornton et al. 1998). Thus, SNe are important players in cosmic evolution.

A star needs to be massive enough to experience core collapse at the end of its life. It has been established that in general stars with initial masses of around 8 M and above finish their lifetimes with a CC (see reviews by Langer 2012; Smartt 2009, from both theoretical and observational standpoints). One lingering question is how the observed diversity in SNe relates to the distribution of massive star parameters in terms of initial mass and pre-SN evolutionary state, or in other words, what kind of massive star gives rise to what type of CCSN?

Some of the massive star progenitors of SNe have been directly identified in archival images taken before the explosion. This method currently provides the most reliable constraints on the physical parameters of SN progenitors. The hydrogen-rich type II SNe, which are the most abundant of the CCSN subtypes (Li et al. 2011), are found to be predominantly produced by red supergiant (RSG) progenitors (Smartt et al. 2009). These SN II progenitors are born with initial masses between ~8 and 20 M and spend their main sequence lifetime as late-O/early-B spectral type stars. By the end of stellar evolution, they still retain most of their H envelope and appear as RSGs, before eventually exploding as type II SNe.

The H-poor CCSNe1, which encompass types Ib (H spectral lines absent, He present), Ic (both H and He absent), and IIb (He present, little H) (Filippenko 1997; Gal-Yam 2016), are thought to arise from massive stars that have lost most of their outer envelope due to some mechanism. One way for a massive star to lose its outer envelope is through metallicity-driven stellar winds. With higher metallicity, the stellar wind becomes stronger due to increased line opacity, thus the rate of mass removed from the star would be higher (Z0.6–0.7 Vink et al. 2001). The progenitors of SESNe are thought to be classical Wolf–Rayet (WR) stars that are without an envelope and exhibit strong stellar winds. At solar metallicity, the lower limit of the initial zero age main sequence (ZAMS) mass for a star to eventually evolve into a WR star is around 25 M (Crowther 2007).

An alternative to the metallicity-driven wind mechanism is mass loss via close binary interaction (Podsiadlowski et al. 1992). In this scenario, the SN progenitor star is in a close binary system, and as stellar evolution proceeds the progenitor expands and overfills its Roche lobe, resulting in mass transfer and significant mass loss. With such a mechanism, the star does not need to be initially as massive as WR star progenitors (≳25 M) in order to lose its outer envelope and eventually explode as an SESN. The detections of SESN progenitors in pre-explosion images in recent years seem to support this binary scenario.

Five SN IIb progenitor candidates have been identified in pre-explosion images: SNe 1993J, 2008ax, 2011dh, 2013df, and 2016gkg. It is very interesting that for all of these cases, the observations seem to be consistent with the binary progenitor scenario (Maund et al. 2004; Folatelli et al. 2014, 2015; Maeda et al. 2015; Kilpatrick et al. 2017; although in some cases a single-progenitor scenario cannot be ruled out, Maund et al. 2015; Sravan et al. 2018). For these cases, the derived progenitor mass is typically higher than most SN II progenitors (~15 M).

In contrast to SN II and IIb progenitors, efforts to detect the progenitors of SNe Ib and Ic have proved unfruitful (see e.g. Eldridge et al. 2013). Currently, there is only one case of SN Ib/Ic progenitor detection, that of iPTF13bvn, an SN Ib (Cao et al. 2013). The progenitor star was constrained to be a sub-WR-mass star (MZAMS < 25 M) in a massive binary system through a number of independent methods (Bersten et al. 2014; Kuncarayakti et al. 2015; Folatelli et al. 2016; Eldridge & Maund 2016; Hirai 2017). Along with the observational constraints on SN IIb progenitors, this questions the importance of massive, single WR stars in the production of SESNe.

An independent way to constrain SN progenitors is by studying their environments. Unlike direct detection methods that are strictly limited to the availability of useful archival images, and analyses based on SN light curves and spectra that are restricted to the time window when the SN Is still visible, environment studies do not depend on these limitations and still can be performed much later after the SN has faded. Therefore, it is possible to build statistically large samples of SN environments and derive strong progenitor constraints from there (see Anderson et al. 2015, for a review). From the environments, the estimate for metallicity and age of the stellar population that gave rise to the SN can be derived. This subsequently reflects the birth metallicity and mass of the SN progenitor, which are the two fundamental factors driving massive star evolution up to the terminal SN explosion (see e.g. Heger et al. 2003; Georgy et al. 2009).

The birth mass and metallicity of SN progenitors are thought to strongly affect the mass loss, and thus the degree of envelope stripping and the eventual SN type. Several studies have used proxies such as host central metallicity and SN radial distance from host centres for metallicity estimates of SN progenitors (e.g. Prieto et al. 2008; Anderson & James 2009; Taddia et al. 2016). A more direct approach to measure metallicity is by obtaining the spectrum of the explosion site, and using strong line diagnostics to derive metallicity. Anderson et al. (2010) and Leloudas et al. (2011) showed that the gas-phase metallicities measured at SN explosion sites do not show statistically significant differences for SNe of different types, while exhibiting an overall trend of higher metallicity for SESN explosion sites. Modjaz et al. (2008) showed that the environments of broad-lined type Ic SNe (IcBL) that are associated with long gamma-ray bursts (GRB) are more metal-poor compared to their GRB-less counterparts. Taddia et al. (2015) studied the explosion site metallicities of interacting SNe and concluded that the metallicity of SNe IIn is different compared to that of SN impostors. Combining direct and indirect metallicity measurements, Taddia et al. (2013) showed that SNe similar to SN 1987A prefer low metallicity environments.

If SESNe are produced by classical WR stars that are single and massive (≳25 M) at birth, there may be a correlation between birth mass and degree of envelope stripping, such that the most massive stars would lose H and He more easily compared to their lower mass counterparts. This scenario implies that SESNe are, on average, more massive than SNe II. Anderson et al. (2012) investigated this issue using a pixel statistics technique (Fruchter et al. 2006; James & Anderson 2006) that essentially constrains the association of SNe with Hα emission from a nearby H II region. As stars age, they drift further from the star-forming region birthplace and show less association with the Hα emission. Thus, more massive stars (including SN progenitors) should show higher association with star formation than the lower mass ones. It was found that indeed SNe Ic show the highest association with star formation, thus highest progenitor mass, followed by SNe Ib, which somewhat resembles the SNe II distribution. Kangas et al. (2017) went further with this technique and compared the pixel statistics of SNe with that of main sequence and evolved stars. They confirmed the finding of Anderson et al. (2012), and suggested the following pairs of SN types and progenitors: SNe Ic are consistent with WN stars ≳ 20 M, while SNe Ib and II occupy the lower mass range between around 9 and 15 M.

With the recent deployment of integral field spectroscopy (IFS/IFU spectroscopy), the field of SN environments has received a significant boost. Integral field spectroscopy allows the collection of both spatial and spectral information of the environments in an efficient manner, while at the same time returning a large amount of information. With IFS, the spectra of the exact explosion site, the immediate local surroundings, and – depending on the size of the IFU field of view(FoV) – the host galaxy can be obtained simultaneously. Thus, in contrast to slit spectroscopy, IFS offers the possibility of identifying and extracting the parent stellar population of the SN rather than simply integrating over the slit aperture (Kuncarayakti et al. 2013a; Kuncarayakti et al. 2013b), and also comparing its properties to the wider environment and the rest of the host galaxy (Galbany et al. 2014, 2016a,b, 2017; Krühler et al. 2017).

In this paper we aim to refine the work presented in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b) in constraining the mass and metallicity of SN progenitors by way of spectral analysis of the parent stellar population. Similarly, IFU spectroscopy of the SN explosion site was utilized to identify the parent population and extract the spectrum, from which the SN progenitor properties were derived assuming coevality. In total, 24 SN sites were used in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b), and with the new dataset presented in this study, this number has now increased more than threefold (83 SN sites). Furthermore, the current study involves other CCSN subtypes not covered in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b), namely SNe IIb, IIn, and IcBL, and thus offers a more complete understanding of massive star evolution and the endpoint SNe.

Following the introduction, the observations and data reductions are presented in Sect. 2. Analyses of progenitor metallicity, initial mass, and star formation history at the explosion site are presented in Sect. 3, followed by discussions in Sect. 4. Section 5 concludes the paper.

2 Observations and data reductions

2.1 Sample description

The SNe used in this work were selected from the Asiago Supernova Catalog (Barbon et al. 1999). The following selection criteria were used:

  • non-thermonuclear (type Ia) SNe;

  • δ ≤ 30°, to enable observations from the southern hemisphere;

  • host galaxy redshift cz ≤ 3000 km s−1, to limit the survey volume;

  • host galaxy inclination ≤ 65° (following Crowther 2013), to minimize chance superposition and other effects that might be attributed to a line of sight through the host galaxy disk;

  • only relatively recent SNe, discovered not earlier than 1970 and not after January 2013, were selected. This is to ensure reliable typing/position and minimal contamination from the late-time SN emissions2.

New IFU data of 62 SNe following these criteria were obtained between 2014 and 2015. One object, SN 1992bd (type II), was not considered in the analysis because its is located very close to the host galaxy centre which contains an active nucleus; this makes contamination-free extraction of the parent stellar population light difficult. Twenty-two3 SNe from Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b) were added to this sample to construct the total sample of 83 SNe (see Table A.1). In the subsequent analyses, all non-interacting H-rich type II SNe are collectively put under one group of SN II without considering the traditional subtypes of IIP and IIL as these hydrogen-rich explosions seem to comprise a continuum without obvious dichotomy in the light curve shape (Anderson et al. 2014; Sanders et al. 2015; Galbany et al. 2016c; Rubin et al. 2016).

Figure 1 shows the pie chart of relative frequency of each SN subtype in the sample. For comparison, the relative frequencies of SNe from the LOSS survey (Smith et al. 2011a) are also shown. The relative frequency of SNe in the two samples are consistent when the sample is divided into hydrogen-poor and hydrogen-rich. Within the hydrogen-poor SN subset itself there are differences, but they are not statistically significant. The relative frequencies of Smith et al. (2011a) is recently updated by Shivvers et al. (2017), however, the hydrogen-poor to hydrogen-rich SN ratio is effectively the same.

Figure 2 shows the histogram of the SN host distances. All of the objects are closer than 40 Mpc; 90% are within ~30 Mpc and the median distance is 17.6 Mpc. At these distances, the typical spatial resolution in the optical data as constrained by natural seeing is better than 100 pc. To date, the current study has the shortest median distance and thus the highest spatial resolution compared to previously published studies focusing on SN environments (see e.g. Anderson et al. 2010, 2015; Leloudas et al. 2011; Modjaz et al. 2011; Sanders et al. 2012; Galbany et al. 2014, 2016a).

thumbnail Fig. 1

Relative frequency of each SN type in the current sample (large pie chart), compared to that of LOSS (Smith et al. 2011a, small pie chart in pale colours).

thumbnail Fig. 2

Histogram of SN host distance. The upper abscissa indicates the projected linear scale in pc/′′. The verticaldashed line indicates the median value.

2.2 Optical observations and data reductions

All optical observations for the new dataset were conducted with the Very Large Telescope (VLT) at the Cerro Paranal Observatory, Chile, in classical visitor mode. The instruments VIMOS (Le Fèvre et al. 2003) and MUSE (Bacon et al. 2014) were used for this purpose.

VIMOS was used in IFU mode and with a medium-resolution (MR) grating, while MUSE was used in the Wide Field Mode. VIMOS gives 0.33′′ /spaxel scale within a 13′′ × 13′′ FoV, while MUSE gives 0.2′′/spaxel within a 1′× 1′ FoV. Table 1 lists the instrument configurations and the resulting spatial and spectral characteristics. The observations were performed over approximately one year (2014–2015) in five different runs: 2014 April, 2014 July, 2014 November, 2015 February (VIMOS4), and 2015 May (MUSE5).

Table A.1 lists the IFU instruments used for each object, and the typical seeing size during the exposure. The sky conditions during the observations were generally clear, although some were conducted in slightly cloudy conditions. Spectrophotometric standard stars were observed for the purpose of flux calibration.

VIMOS observations were done in 2 × 1800 s on-source exposures, targeting the SN sites. In MUSE observations, the whole SN host galaxy was observed, taking advantage of the wide field of view. In several instances multiple pointings were taken due to the large extent of the galaxy. For each MUSE pointing, four dithered positions of 450 s exposure time were integrated on-source resulting in a total exposure time of 1800 s per pointing. A separate sky pointing was also observed for each object for the purpose of background sky subtraction.

Raw data reduction was carried out using VIMOS and MUSE data reduction pipelines, run using the Reflex interface (Freudling et al. 2013). For both VIMOS and MUSE raw data, the reduction procedures of bias subtraction, flat-fielding, wavelength, and flux calibration were applied; and finally cube reconstruction was performed. MUSE cubes were sky subtracted using the blank sky pointings and further corrected for atmospheric effects using the Zurich Atmospheric Package (ZAP; Soto et al. 2016). Sky subtraction in VIMOS data was done in the subsequent spectral extraction step, where the sky background was defined with an annulus around the object extraction aperture.

From the IFU datacube, the parent stellar population of the SN was identified. Physically, these are young stellar clusters/H II regions lying within one seeing radius from the SN positions. The SN coordinates were obtained from the Asiago SN database and were used to localize the SN position in the datacube. The World Coordinate System (WCS) of the datacube was generated by the reduction pipeline and reflects the telescope pointing accuracy. The accuracy of the WCS coordinates was around 1′′ –1.5′′. This was determined using foreground stars and several SNe IIn that were still visible. Whenever possible, the SN position in the datacube was further checked with available broad-band images of the SN or reported offsets towards the centre of the host galaxy. A one-dimensional spectrum of the stellar population was extracted from the datacube by using an aperture, whose radius was set to the size of seeing FWHM. Visualization and extraction were done using the QFitsView6 (Ott 2012) software. Analysis of the one-dimensional spectrum was carried out using IRAF7, as described in the next section.

Table 1

Instrument configurations used in this study.

2.3 Infrared observations and data reductions

The optical IFU data are supplemented with near-infrared (NIR) IFU observations, primarily for the purpose of constraining star formation history (see Sect. 3.3).

VLT/SINFONI (Eisenhauer et al. 2003; Bonnet et al. 2004) was used in the K band with the 0.1′′ /spaxel scale, giving a field of view of 3′′× 3′′ (Table 1). Adaptive optics were used with either natural guide stars or laser guide stars when applicable, resulting in near diffraction-limited observations. The observations were conducted between 2012 and 2015 in several runs in both visitor and service modes8. The typical on-source total integration time per object is 6900 s, taken in 300 s dithers for both object and sky frames. Ten SN sites from the optical sample were observed.

The raw data were reduced using the SINFONI instrument pipeline within the GASGANO interface, which includes the standard procedures of bias subtraction, flatfielding, distortion correction, wavelength and flux calibration, and sky subtraction. The resulting datacubes were subsequently analysed using QFitsView and IRAF. The uncertainty of the SN position in the SINFONI datacube was estimated to be similar to the 1′′ –1.5′′ uncertainty in VIMOS/MUSE datacubes, as the instruments use the same VLT telescope. Due to the narrow SINFONI FoV, no foreground stars could be used to check the WCS accuracy. In some cases, the general appearance of the SN field in the K-band image generated from the datacube were compared to the optical image (see Sect. 3.3), which was used to confirm the above-mentioned estimate.

3 Analysis and results

3.1 Metallicity

Gas-phase abundance of oxygen is used as the proxy for metallicity. The widely used N2 index was used to calculate 12 + log(O/H), according to the calibration by Marino et al. (2013). Due to the unavailability of the Hβ line in the data, the O3N2 index was only available for a part of the sample, thus the N2 index was used for the whole sample. The N2 index uses the ratio of emission lines Hα and [N II] λ6584, which are closely spaced in wavelength and thus robust against uncertainties introduced by reddening or flux calibration. As the 12 + log(O/H) values from Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b) are originally in the Pettini & Pagel (2004) scale, the corresponding 12 + log(O/H) values on Marino et al. (2013) scale were calculated and used. The solar oxygen abundance was taken as 12 + log(O/H) = 8.69 (Asplund et al. 2009).

The Hα and [N II] λ6584 lines were measured in the one-dimensional spectrum by fitting a Gaussian curve after subtracting the stellar continuum using a polynomial function. This approach does not make assumptions regarding the underlying stellar population behind the gas; therefore, it does not take into account the underlying stellar absorption. The resulting Hα and [N II] line strengths are then directly used for the calculations of the N2 index and subsequently metallicity. In several cases, no emission was detected at the explosion sites, or the SN was still bright at the time of the observation9. For such cases, the nearest H II region spectrum is extracted and the metallicity estimate was adopted for the metallicity value at the explosion site. This procedure is only applicable for SN sites observed with MUSE, due to the large field of view allowing the observation of the whole host galaxy (SNe 1978G, 1983K, 1988E, and 2011fh). In these cases, the average projected offset from the SN positions to the nearest H II region is 540 parsec.

In addition, metallicity from the calibration of Dopita et al. (2016) was also calculated. Their new calibration similarly uses Hα and [N II] λ658410 lines, and additionally [S II] λλ6717, 6731. The results of metallicity analysis are unchanged when this calibration is applied.

The left panel in Fig. 3 shows the cumulative distribution function of metallicity of different SN types. The SNe are colour-coded according to their subtypes. The distributions do not show any particular pattern and are practically superposed one on top of another. The lowest metallicity bin of the sample is occupied by a number of SNe II and Ic only, between 12 + log(O/H) = 8.0 and 8.2 dex. Metallicity measurements for extreme SNe such as superluminous (SL) SNe (Leloudas et al. 2015) and SNe IcBL (Modjaz et al. 2011), both converted into the Marino et al. (2013) scale, are also shown for comparison. It is readily evident that the metallicity of normal CCSN events is typically higher than those. The two SNe IcBL in our sample, SNe 1998bw and 2009bb, have 12 + log(O/H) metallicities of 8.30 and 8.49 dex, respectively.

Two-sample Kolmogorov–Smirnov (K–S) statistical tests were performed to analyse whether two different SN types arise from the same population. The right panel in Fig. 3 shows the matrix for the K–S test result. None of the compared SN subtype pairs shows significant difference, with K–S p < 10%.

thumbnail Fig. 3

Left: observed cumulative distributions of 12 + log(O/H) for different SN subtypes. Colour-coding is as follows: SNe Ic blue,Ib green, IIb orange, IIn purple, II red, and SESNe (Ic + Ib + IIb combined) dashed cyan. For comparison, the literature distributions of SN IcBL (dashed line), SLSN-I (dash-dotted line), and SLSN-II (dotted line) are also shown. Right: matrix for the K–S test result between different SN subtypes. The colour bar indicates the probability that samples were drawn from the same metallicity distribution, where the black end indicates significantly different distributions.

3.2 Progenitor age and initial mass

Stellar evolution is constrained by the amount of nuclear fuel available for burning in the stellar core. More massive stars burn fuel more rapidly than their less massive counterparts. Within a coeval stellar population, stars that were born with the highest mass have the shortest lifespans and die out first, leaving behind stars with lower masses. Such a simple stellar population is formed out of a single, homogeneous molecular cloud in an instantaneous burst of star formation, hence all the stars share the same age and metallicity. Therefore, by determining the age of the stellar population it is possible to estimate the age, and thus lifetime, of the last star that died out. This star was observed as the SN, and all the current remaining stars in the stellar population must have initial mass lower than that of the SN progenitor star.

In this work, as in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b), the stellar population that is present at the SN explosion site is assumed to be the parent stellar population from which the SN progenitor emerged. The age of the stellar population is taken as the SN progenitor lifetime span, and this is subsequently converted into an initial mass estimate. Due to a number of possible factors this initial mass estimate is effectively an upper limit, i.e. the actual SN progenitor masses are more likely to be lower than those derived from stellar population age. For example, including continuous star formation or a binary population may extend the lifetime of Hα emission used in age determination, thus the stellar population age derived with such assumptions will be considerably older than the instantaneous-burst or single-star populations for a given Hα equivalent width (EW; see e.g. Leitherer et al. 1999; Crowther 2013; Eldridge et al. 2017).

In this study, the age indicator HαEW is used to constrain the stellar population age. The indicator is measured as the ratio between the emission line and continuum fluxes at Hα. It is effectively a measure of the number ratio between the ionizing stars of OB spectral class responsible for the Hα emission line and the other non-ionizing (lower mass) stars in the stellar population, assuming a fixed IMF. As the stellar population ages, these ionizing stars will be reduced in number while the number of lower mass stars remains constant, therefore resulting in the decline of HαEW.

Measurements of HαEW were performed on the one-dimensional spectra of the stellar population, using IRAF/splot. The spectral continuum was normalized using a polynomial function and the emission line was fit using a Gaussian function. As in the metallicity estimate, this measurement involves a very small window in wavelength space, thus is robust against uncertainties introduced by reddening and flux calibration.

Simple stellar population (SSP) models from Starburst99 (Leitherer et al. 1999) were used to compare the observed age indicators and infer the age of the stellar population. These models have been shown to be reliable for analysing young stellar populations (Kuncarayakti et al. 2016b). While HαEW is the primary age indicator used in this study, a number of other spectral lines are also useful as age indicators. These include the EWs of Brγ and CO ν = 2–0 band (the 2.3 μm overtone; hereafter CO 2.3 μ) in the K band (see Sect. 3.3), and also to a limited extent the near-infrared Ca-triplet lines (see Kuncarayakti et al. 2013a; Kuncarayakti et al. 2013b). The observed EW is compared to the SSP model at the corresponding metallicity to estimate the stellar population age. Figure 4 shows the evolution of these age indicators.

The observed HαEW corresponds to the age of the stellar population, through SSP models. Starburst99 SSP models in the corresponding metallicities (Z = 0.004, 0.008, and 0.02) were used, assuming single stars born in an instantaneous star formation and distributed in mass according to the Salpeter IMF. The cumulative distribution of the measured HαEW of the SN parent populations is given in the left panel in Fig. 5. This is the empirical measure of the SN progenitor lifetime in its rawest form. In the figure it is apparent that the SN IIn distribution prefers low HαEW, followed by SN IIb, then SN Ib and SN II which are quite close to each other, and finally the SN Ic distribution that shows a preference for high HαEW. The majority of SN IIn explosion sites do not show any Hα emission, while nearly all SN Ic sites have significant Hα emission. Even if the next nearest H II regions to the SNe IIn are used for the HαEW measurement (see Sect. 3.1), it is still apparent that SNe IIn have the lowest HαEW among the SN types. The result of the K–S test in the right panel in Fig. 5 confirms that SNe IIn are statistically different from most of the other SN types. On the other hand, the differences between the other SN types in terms of parent population HαEW are not statistically significant.

The parent population age, in turn, translates into the age of the previous most massive star in the population which exploded as an SN. This is essentially the lifetime of the progenitor star, which is governed by its initial mass. Stellar evolution models by Bressan et al. (1993, for Z = 0.02) and Fagotto et al. (1994, for Z = 0.008 and 0.004) were used to derive the star initial mass from the lifetime. Figure 6 shows the cumulative distribution of the derived age and initial mass for different SN subtypes. A similar trend to that in Fig. 5 is observed: SNe Ic occupy the young age and high-mass end, while SNe IIn are characterized by old age/low mass, and in between the SNe II, IIb, and Ib together show similar distributions in age and initial mass. Again, the differences between the SN subtypes are not found to be statistically significant, except in the case of SNe IIn.

3.3 Star formation history

As discussed in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b), one of the most important caveats in the determination of stellar population age is the uncertainty in the star formation history. The age-dating of the stellar population is based on the assumption that star formation was instantaneous, while it is uncertain whether this is actually the case. What is observed is the combined light of many different stars and interstellar gas clouds, and it is not possible to ascertain the actual history of star formation in the region. However, using a combination of age indicators it is possible to probe distinct stellar populations of different ages, thus giving clues to whether the star formation was instantaneous or not.

Using near-infrared IFU data obtained with SINFONI, here we attempt to alleviate this problem. The age indicators Brγ emission EW and CO 2.3 μ absorption EW probe entirely different SSP ages. Brγ exclusively probes stellar populations younger than ~8 Myr, while CO probes the population older than that (Fig. 4). Moreover, CO EW values between 11 and 14 Å indicate stellar population ages corresponding to the typical SN II progenitor mass of 8–16 M obtained from direct progenitor detection (Smartt et al. 2009). If both indicators are observed to be present at the same location, it may suggest thatthe star formation in the area did not proceed in one single burst (i.e. there are two distinct populations with average ages of <8 Myr and >8 Myr). This immediately separates the regions with single or multiple bursts of star formation. Although the actual star formation history is not recovered, being able to assess whether the star formation was a single burst or not is an important point for the current work. In addition, this dividing line at ~8 Myr corresponds to the lifetime of ~25 M stars, which is thought to be the lower limit of the ZAMS mass of single WR stars (Crowther 2007). Thus, the detection of strong Brγ emission with EW corresponding to ≲8 Myr may be used to constrain the presence of such massive stars.

Brγ emission and CO absorption are observed in some of the SN sites in the SINFONI sample. In the SN sites where Brγ and CO appear together, no spatial correlation between the two lines is observed (e.g. Fig. 7, panels (iii) and (iv)). A more quantitative description of this observation is obtained by analysing the pixel statistics (Fruchter et al. 2006; James & Anderson 2006) of Brγ and CO. In the pixel statistics technique, the IFU spaxels are ranked in a normalized cumulative rank (NCR) in which the brightest pixel has NCR = 1, while the faintest has NCR = 0. This is done for both Brγ and CO, and additionally for the Hα and [N II] emission lines in the optical data. In panel (v) in Fig. 7, the pixel statistics for the two NIR age indicators are plotted as contours. No correlation is observed between the two, while for comparison the pixels containing Hα and [N II] show high correlation. This means that the brightest Hα pixels are also the brightest [N II] pixels, while there is no such behaviour for Brγ emission and CO absorption observed.

In general there is no spatial correlation between Brγ and CO within an area of a few hundred pc at the explosion sites (typical SINFONI FoV size for the objects in the sample); however, it has been observed that these two indicators can still appear together at the same position. Comparing the spectra of regions A (strong K-band source, indicative of high stellar mass) and B (strong, localized Brγ emission) in Fig. 7, panel (vi), it is apparent that both regions contain the Brγ and CO lines simultaneously.

Table 2 lists the SN explosion sites for which we have collected NIR SINFONI data, together with the values of Brγ and CO EWs and the corresponding ages. To measure the Brγ and CO EWs at the SN position, a one-dimensional spectrum was extracted from the SINFONI datacube within one optical seeing radius. Then, Brγ and CO lines in the spectrum were measured in the same fashion as the optical lines, and similarly the EW values were compared to Starburst99 SSP models to obtain the age estimate. In the last column of Table 2, the stellar population age derived from HαEW is presented for comparison. The ages derived from Brγ and Hα are consistent with each other within ~15%, as expected if these line emissions come from the same component in the stellar population.

In the SN sites where both Brγ and CO are detected, the corresponding ages show that these two age indicators indeed do not probe the same age range. It is interesting to note that the derived CO ages, within the errors, correspond to the age range that give rise to RSG progenitors. On the other hand, the presence of Brγ emission suggests that more massive stars are also present at the same location. Therefore, these instances illustrate the case where there could be multiple star formation bursts in an SN explosion site. While the youngest starbursts generally occur in compact, localized regions (Brγ, Fig. 7 panel (iii)), the older populations are more spread throughout the region (CO, Fig. 7 panel (iv)), forming a background. Therefore, the old population may contaminate the young population in the same line of sight. As the SN exploded in this kind of region, the non-instantaneous star formation history may confuse the age determination. The Hα/BrγEW method that is used in this work is sensitive towards the youngest stellar populations at the massive end of the IMF; in thiscontext, contamination from older stellar populations could be considerable and as a result the Hα/BrγEW age is to be treated effectively as a lower limit, i.e. an upper limit for the progenitor mass.

Clues on the star formation history can also be obtained from the host H II region itself. Larger H II regions tend to contain multiple stellar populations resulting from several past starbursts, while smaller more compact H II regions tend to have a nearly instantaneous starburst (Crowther 2013). These giant H II regions typically have luminosities of the order of log L(Hα) ~ 39 erg s−1 and sizes of a few hundred pc. We examine the H II regions at the SN sites in our sample to identify the type of H II region (hence star formation history) that dominates our sample. Figure 8 shows the cumulative distribution of the observed flux and luminosity of the H II regions. Most of the H II regions (90%) have a luminosity lower than log L(Hα) = 39 erg s−1. This suggests that our SN sites are dominated by relatively small H II regions that are expected to be characterized by nearly instantaneous star formation history. Nevertheless, even though efforts have been made to ensure that we are probing small regions, we note that the effects of chance superposition between the SN and H II region may always be present in this kind of environment study. Monte Carlo simulations predict that this could affect up to ~50% of the observed associations (Kuncarayakti et al. 2013b). As seen in the NIR IFU data, in three out of five SN sites where a young population is detected via Brγ emission, a potentially older population is also detected via CO absorption. The risk of chance superposition increases with lower mass progenitors, which are older (longer-lived) hence have more time to drift away from the actual star formation birthplace. In addition, studies of resolved stellar populations around type IIP SNe show that there are cases where relatively low-mass progenitor candidates (~8 M) coexist with higher mass stars (15–60 M) within 100 pc of the SN position (Maund 2017). Only when an appropriate age component can be identified from the mixed stellar population can the SN progenitor initial mass be tightly constrained. While this caveat introduces more scatter in the statistics when comparing SNe of different types, general trends of differences between different SN types are still expected to be observed.

thumbnail Fig. 4

Evolution of Starburst99 age indicators, shown for Z = 0.02. The upper abscissa represents the stellar lifetime at the corresponding SSP age. Hα is indicated in red, Brγ green, CO 2.3 μ purple, and CaT blue. The CO EW values where typical SN II progenitors (MZAMS ≈ 8–16 M) are expected to be found are encompassed within the purple shaded area.

thumbnail Fig. 5

Left: observed cumulative distribution of log(HαEW) for different SN subtypes. SN sites which show no Hα emission are plotted as having log(HαEW) = 0. The dashed purple line indicates SN IIn distribution using the nearest H II regions (see text for description). Right: matrix for K–S test result. The colour-coding is the same as in Fig. 3.

thumbnail Fig. 6

Left and centre: observed cumulative distribution of parent stellar population age and progenitor initial mass for different SN subtypes. SN sites which show no Hα emission are plotted as having an age of 40 Myr, corresponding to the single-star SN progenitor lower mass limit of 8 M. Right: matrix for K–S test result. The colour-coding is the same as in Fig. 3.

thumbnail Fig. 7

Maps of the explosion site of SN 2009dq generated from the datacubes, in continuum-subtracted Hα (i), K-band (ii), Brγ (iii), andCO (iv). As the maps are color-coded with white for high-count pixels and black for low-count pixels, the pixels with strong CO absorption appear black in panel (iv). The white 3′′   × 3′′ square in panel (i) represents the SINFONI FoV, and the white circle denotes the SN position within 1′′ radius. North is up and east is left in all maps. Cyan circles in panels (iii) and (iv) indicated with A and B have identical radii of 0.25′′. Panel (v) shows the SINFONI spaxels plotted according to the Brγ and CO NCRs; these are shown as contours at 10%, 50%, and 90% levels, while for comparison, the VIMOS spaxels are plotted in red plus marks according to the Hα and [N II] NCR. Panel (vi) shows the extracted spectra of regions A and B, in the Brγ and CO spectral regions. The two lines are indicated with green-shaded regions. Wavelength is in the observer frame.

Table 2

SN sites observed with SINFONI.

thumbnail Fig. 8

Cumulative distribution of host H II region flux and luminosity.

4 Discussion

Massive stars are thought to end their lives as SNe after stellar evolutionary processes have guided them from their ZAMSphase through to the pre-SN stage. The initial birth mass and metallicity are considered to be the most fundamental physical parameters driving stellar evolution (see e.g. Heger et al. 2003). As the strength of mass loss via stellar wind scales with stellar mass and metallicity, one would expect the most highly stripped SNe to be associated with progenitor stars with the highest mass and metallicity, in comparison to the other SN subtypes.

As shown in Fig. 3, the differences in metallicity between different SN subclasses are not significant. This is in contradiction with what is expected from single-star evolution theory, where metallicity-driven winds are crucial: type Ic SNe, which are the most highly stripped, would show the highest metallicity, followed by type Ib, and finally the H-rich type II SNe. The observations, on the other hand, reveals that this is not the case. Some SNe Ic are even located in the low-metallicity part of the distribution in the current sample. This result strengthens the notion that metallicity may not play an important role in deciding theresulting SN type, in agreement with other works based on SN environments (Anderson et al. 2010, 2015; Leloudas et al. 2011; Galbany et al. 2016a). The environments of broad-lined SNe IcBL are found to be relatively metal poor compared to the normal CCSNe, in agreement with previous studies (Modjaz et al. 2011; Galbany et al. 2016a). However, we note that there are only two such SNe in the current sample. The explosion site of SN 1998bw (the first SN to be associated with a GRB: 980425; Galama et al. 1998; Krühler et al. 2017) in this study shows a lower metallicity of 12 + log(O/H) = 8.30 dex compared to the GRB-less SN 2009bb (Pignata et al. 2011), 12 + log(O/H) = 8.49 dex. Levesque et al. (2010a), using slit spectroscopy of the explosion site, concluded that the high metallicity of SN 2009bb site is consistent with typical GRB-less SNe IcBL and not with GRB hosts. Their metallicity value recalculated on the Marino et al. (2013) N2 scale is 12 + log(O/H) = 8.52 dex. These two different cases illustrate the importance of metallicity in deciding whether an SN IcBL progenitor would also produce GRB or not (Modjaz et al. 2008; Levesque et al. 2010b). Progenitors with higher metallicity are not able to spin fast enough and thus produce high angular momentum essential for GRB jet production, eventually producing a GRB-less SN IcBL (Woosley & Bloom 2006).

The strong dependence of stellar wind with metallicity in massive stars should also be reflected in the number ratio of H-poor and H-rich SNe within different metallicity bins. With high metallicity, wind and mass loss become stronger, resulting in more H-poor massive stars and eventually SESNe. Figure 9 (upper panel) shows the observed number ratio between H-poor SESNe and H-rich SNe II. In general, the trend over metallicity is flat11. This suggests that other factors more dominant than metallicity could be at play in producing SESNe. The picture changes when one examines the SN Ic/SN Ib number ratio (Fig. 9, lower panel). As metallicity increases, a rise in the Ic/Ib number ratio is seen. This suggests that within the SESN group, metallicity is affecting the production of SNe Ic more than it does Ib. Indeed, in our sample SNe Ic are found with higher median metallicity than are SNe Ib. A similar trend is also found by Galbany et al. (2016a), where the association of SN Ic with high metallicity is even more pronounced in non-targeted surveys. More SNe Ic are produced in high-metallicity environments (cf. Arcavi et al. 2010; Arcavi 2011; Kelly & Kirshner 2012; Graur et al. 2017b), indicating that metallicity-dependent stellar winds are indeed affecting the production of this particular subclass. This indirect evidence leaves open the possibility that massive single stars with strong winds are contributing to the production of SN Ic, at least more strongly than SN Ib and possibly also IIb. However, we note that binary progenitor models also predict that to some extent metallicity affects the SESN production in a way that higher metallicity would still produce more highly stripped progenitors (Eldridge et al. 2017; Yoon et al. 2017). A standard binary model typically leaves a thin H layer after the Roche-lobe overflow, and as it is difficult to remove the He layer this way, a final push by metallicity-driven wind is thought to be important in expelling the remaining layer and in producing an SN Ic progenitor.

A more stringent constraint for progenitor mass comes from the age analysis. The observed HαEW and the derived age/initial mass points toward higher mass progenitors for SNe Ic. SNe Ib, and IIb, while similarly deprived of the outer hydrogen envelope, correspond more closely to the lower mass SN II progenitors (see also Kangas et al. 2017). Again, this suggests that these two He-rich subclasses are probably dominated by sub-WR mass progenitors that have lost the outer envelope through binary interactions. However, the small differences means that these results are not statistically significant.

Evidence of binary interaction affecting the production of SESNe is mounting. Smith et al. (2011a) argued that the observed relative fractions of CCSN subtypes cannot be reconciled with IMF calculations if only single WR star progenitors are expected for SESNe. There are simply not enough massive stars above 25 M to account for all the observed SESNe. Lyman et al. (2016) showed that the explosion parameters of most SESNe do not correspond to very massive progenitors, but to moderately massive stars that were likely the product of binary interaction. While it has been recently established that the majority of massive stars experience interaction with a binary companion during their lifetimes (Sana et al. 2012), population synthesis studies further suggest that interaction may bring stars of intermediate mass (4–8 M) to core collapse as a result of mass transfer or merger (Zapartas et al. 2017). This accounts for 14 14+15$^{&#x002B;15}_{-14}$% of all CCSNe, and results in a distribution tail of “delayed” CCSNe that occur in stellar populations older than 50 Myr where all the massive stars above MZAMS ~ 8 M have already disappeared. In this case, these SNe would appear in regions without Hα emission. Within our sample, 20% of the SNe occur in such non-star-forming regions. Zapartas et al. (2017) speculated that the various SNe from this delayed population may not form a homogeneous class, but also noted that some type II and IIn SNe can be producedthrough this scenario.

It is clear that in order to explode as SESN the progenitor star should experience significant mass loss, whose mechanism may take the form of strong binary interaction or metallicity-driven winds. Type IIn SNe, on the other hand, are evidently interacting with nearby CSM that is thought to be the result of mass loss activity of the progenitor. At least some of these SNe IIn are associated with stars of extremely high initial mass that resemble luminous blue variable (LBV) stars (see e.g. Gal-Yam et al. 2007; Gal-Yam & Leonard 2009; Smith et al. 2011b). However, in many instances the environments of SNe IIn do not suggest recent star formation (Anderson et al. 2015). In comparison with other SN types, it has also been shown that SNe IIn exhibit the least association with ongoing star formation compared to the other CCSNe (Habergham et al. 2014). Kangas et al. (2017) showed that the SN IIn population does not share a similar spatial distribution in host galaxy Hα light with LBV stars. They are instead best matched with the RSG stars, whose distribution suggests lower mass compared to LBV stars. An opposing view against LBVs as highly massive stars comes from Smith & Tombleson (2015), who argue that most LBV stars are isolated, being a product of binary evolution, and therefore consistent with SNe IIn exploding inpassive regions. All the evidence points to the majority of SN IIn coming from the low-mass end of massive stars. The current study confirms this view. SNe IIn are found to be the least associated with recent star formation among CCSNe. While the environments suggest a low-mass origin, the fact that some SNe IIn are strongly associated with massive LBV stars would indicate that the SN IIn population is a mix of different kinds of progenitor systems. The majority of SNe IIn are produced by relatively low-mass stars and not massive progenitors. It is interesting to note that while SNe Ic and IIn must both suffer significant progenitor mass loss prior to the SN, they represent the opposite ends in CCSN progenitor initial mass distribution. Nevertheless, the mass loss mechanism for SN Ic and IIn progenitors may be different. SN IIn progenitors must experience significant mass loss at least immediatelybefore the core collapse, within years to centuries (Smith 2017). This suggests that the high mass-loss rate is timed with the core collapse, and may be caused by late-stage burning instabilities (Smith & Arnett 2014). On the other hand, the WR star progenitors of SNe Ic lose mass through relatively steady, strong stellar winds (Crowther 2007; Smith 2014). If SN IIn progenitors only experience severe mass loss immediately before the core-collapse, this could be consistent with a scenario in which low-mass stars with average mass loss rates undergo outbursts before the core-collapse.

A direct link between SNe Ic and IIn has recently been discovered. SN 2017dio showed SN Ic spectral characteristics at the early phase, and signs of an associated H-rich CSM (Kuncarayakti et al. 2018). At later phases, the spectral characteristics became more similar to SNe IIn as the CSM interaction became more intense. This nearby H-rich CSM cannot come from the SN progenitor star itself as it has been stripped of the H and He layers (hence the SN Ic spectral appearance). It has been suggested that the SN progenitor had a secondary companion that was H-rich due to mass transfer from the primary, and that underwent an unstable LBV phase or transfer of mass back to the primary at the time of the explosion. This might explain the presence of the H-rich CSM. Another event, SN 2014C, initially showed a type Ib spectrum around maximum light before transforming into showing type IIn spectra at late times (Milisavljevic et al. 2015; Margutti et al. 2017). This event was interpreted as an explosion of a H-free progenitor star inside a CSM cavity, and as the SN ejecta traverses outward, it encountered the H-rich CSM and produced IIn-like spectral signatures. The CSM must have been produced by mass loss episodes of the progenitor preceding the terminal SN. Such events are rare and highlights our incomplete understanding of pre-SN massive star evolution.

In the context of the whole CCSN subclass, the current study confirms that a strictly single-star progenitor scenario does not likely hold. Figure 10 shows single-star theoretical predictions of massive star evolution on the initial mass-metallicity plane (Georgy et al. 2009). In such a scenario, progenitors of different SN subtypes should not overlap on the mass-metallicity (MZ) plane. However, observational data obtained in this study suggest that they are actually significantly overlapping. Most notable are the H-poor SNe that are less massive than 25 M and fall into the SN II progenitor domain. These SNe are presumably produced by progenitor stars in a binary system, since as single stars they are not massive enough to evolve into WR stars. Theoretical predictions that take into account binary progenitors (e.g. Zapartas et al. 2017; Eldridge et al. 2017) do predict that stars lower than 25 M can remove the outer envelope through binary interaction. These progenitors would appear practically anywhere in the MZ diagram (Eldridge & Stanway, priv. comm.).

However, there could still be some influence of single progenitors in the production of CCSNe. According to the degree of envelope stripping, SESNe can be sorted as type IIb, then Ib, and finally Ic. Although not apparent in metallicity distributions, this pattern is visible in Figs. 5 and 6, albeit without statistical significance. The median initial mass estimates for SNe IIb, Ib, and Ic progenitors are, respectively, 29.7, 31.3, and 32.6 M. The two broad-line IcBL SNe are similarly of 36 M, higher than the median of the other SN subtypes. For comparison, SNe II progenitors have a median initial mass of 29.7 M, similar to SNe IIb. These SN II progenitors have their hydrogen envelope largely intact at the time of the SN and are traditionally thought to occupy the lowest mass range of CCSNe at ≲ 20 M, verified through various independent methods (e.g. direct detections, Smartt et al. 2009, and nebular phase analysis, Jerkstrand et al. 2012). With different prescriptions for mass loss and reddening, however, it is still possible to have RSG SN II progenitors at around 25 M (Beasor & Davies 2016). It is to be noted that the initial mass estimates presented in the current work do not necessarily reflect the true value, but are more likely to be upper limits (see Sects. 3.2 and 3.3). The nature of the environment studies does not allow stringent characterization of each individual event based on the underlying stellar population; nevertheless, the statistics are useful in order to compare and constrain different SN subclasses.

On average, SNe with a higher degree of stripping are associated with younger and more massive progenitors. This mass sequence is expected in single SN progenitor populations (Heger et al. 2003; Georgy et al. 2009), and this trend has been previously shown via SN environment observations (Anderson et al. 2012). Considering previous studies in the literature and the results presented in this work, the view that there is a mix between single and binary SN progenitor populations appears to be consistent with the observed signatures. Massive stellar evolution does not seem to work as a straight pipeline in which a certain fate of SN type awaits a star of a particular initial mass, but mass loss (whether induced by binary interaction or not) does greatly affect the outcome of the evolution (cf. Smith 2014).

thumbnail Fig. 9

Number ratios of different SN types in metallicity bins: N(SESN)/N(SN II) (upper panel), and N(Ic)/N(Ib) (lower panel). Horizontal error bars indicate the bin size and vertical error bars are Poisson errors.

thumbnail Fig. 10

Diagram showing the plane of progenitor star initial mass and metallicity. Shaded areas are stellar evolution predictions from rotating single-star models of Georgy et al. (2009), for SN II (pink), Ib (light green), and Ic (blue-green). Data points from this work are colour-coded: blue for SN Ic, bright green for Ib, red for II, purple for IIn, orange for IIb, and cyan forIcBL. Uncertainties in metallicity reflect the 0.16 dex error in the 12 + log(O/H) N2 calibration of Marino et al. (2013).

5 Summary and conclusions

In this work, we present progenitor initial mass and metallicity constraints for distinct subclasses of CCSNe, namely SNe types II, IIn, IIb, Ib, Ic, and IcBL. As in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b), the parent stellar population of the SN was identified and analysed using IFU spectroscopy, where the integrated spectrum was used to derive the metallicity and age. Assuming coevality between the SN progenitor and stellar population, these two parameters were adopted for the SN progenitor and used forthe derivation of its initial mass estimate. Unlike in Kuncarayakti et al. (2013a); Kuncarayakti et al. (2013b), where the SN sites were preferentially selected for bright H II regions, the current work uses a relatively unbiased, distance-limited approach while still keeping the sample distance short. Ninety per cent of the sample falls within ~30 Mpc distance, with median sample distance of 18 Mpc and typical projected linear size per spatial resolution better than 100 pc.

The following conclusions were derived from the analysis:

  • Metallicity differences between SN types are not statistically significant. The implication is that metallicity does not play a critical role in pre-SN mass loss and deciding the outcome SN. Nonetheless, there is a subtle effect of metallicity on the production of SNe Ib and Ic. The Ic/Ib number ratio tends to increase as metallicity increases. In comparison, the number ratio between SESN and SN II across metallicity is flat.

  • SNe Ic appear to be the most associated with the youngest stellar populations and most massive progenitors. They are followed consecutively by SNe Ib, then IIb and II. SN IIn is at the opposite end of the spectrum, being associated with older populations and less massive progenitors. The differences in progenitor mass estimates are not significant, except when comparing SN IIn with other SN types. On average, the two SNe IcBL appear to have more massive progenitors than the other CCSN progenitors.

  • Near-IR IFS shows that there could be multiple stellar populations with different ages present at the SN site. As Hα probes the youngest stellar populations, mass constraints derived from it should be treated as an upper limit.

  • On the initial mass-metallicity plane, SN progenitors do not conform to single-star stellar evolution predictions. Combined with evidence derived from other methods, this confirms that the progenitors of CCSNe cannot only come from single stars. There must be a significant contribution from massive interacting binaries in CCSN production.

  • Along with initial mass, mass loss is one of the most important parameters in driving massive star evolution and decidingthe endpoint SN. This is apparent when comparing SNe Ic and IIn, where the progenitors of either type must similarly undergo significant mass loss, and yet they populate the two opposing ends of the mass spectrum.

Acknowledgements

We thank the referee for the useful suggestions that improved the paper, and Thomas Krühler for his indispensable help on working with MUSE. HK acknowledges support provided by CONICYT through FONDECYT grant 3140563, and the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. KM acknowledges support provided by Japan Society for the Promotion of Science (JSPS) through KAKENHI Grant 17H02864 and through JSPS Open Partnership Bilateral Joint Research Project between Japan and Chile. LG and KM acknowledge support from FINCA visitor programme. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programmes 089.D-0367, 091.D-0482, 093.D-0318, 094.D-0290, and 095.D-0172. It is a pleasure to thank K. Nomoto, M. Tanaka, S. Mattila, and T. Kangas for fruitful discussions; J. Lyman for alerting us to the discovery of SN 2017ahn; and G. Pignata for sharing the FITS image of SN 2009bb. HK is indebted to the late D.H. Nugroho for his suggestions in working with datacubes.

Appendix A: Long tables

Table A.1

Objects.

Table A.1

continued.

Table A.2

Results.

References

  1. Anderson, J. P., & James, P. A. 2009, MNRAS, 399, 559 [NASA ADS] [CrossRef] [Google Scholar]
  2. Anderson, J. P., Covarrubias, R. A., James, P. A., Hamuy, M., & Habergham, S. M. 2010, MNRAS, 407, 2660 [NASA ADS] [CrossRef] [Google Scholar]
  3. Anderson, J. P., Habergham, S. M., James, P. A., & Hamuy, M. 2012, MNRAS, 424, 1372 [NASA ADS] [CrossRef] [Google Scholar]
  4. Anderson, J. P., González-Gaitán, S., Hamuy, M., et al. 2014, ApJ, 786, 67 [NASA ADS] [CrossRef] [Google Scholar]
  5. Anderson, J. P., James, P. A., Habergham, S. M., Galbany, L., & Kuncarayakti, H. 2015, PASA, 32, e019 [NASA ADS] [CrossRef] [Google Scholar]
  6. Arcavi, I. 2011, in Death of Massive Stars: Supernovae and Gamma-Ray Bursts, (Cambridge University Press) Proc. IAU 734 [Google Scholar]
  7. Arcavi, I., Gal-Yam, A., Kasliwal, M. M., et al. 2010, ApJ, 721, 777 [NASA ADS] [CrossRef] [Google Scholar]
  8. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bacon, R., Vernet, J., Borisova, E., et al. 2014, The Messenger, 157, 13 [NASA ADS] [Google Scholar]
  10. Barbon, R., Buondí, V., Cappellaro, E., & Turatto, M. 1999, A&AS, 139, 531 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Beasor, E. R., & Davies, B. 2016, MNRAS, 463, 1269 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bersten, M. C., Benvenuto, O. G., Folatelli, G., et al. 2014, AJ, 148, 68 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bonnet, H., Abuter, R., Baker, A., et al. 2004, The Messenger, 117, 17 [NASA ADS] [Google Scholar]
  14. Bressan, A., Fagotto, F., Bertelli, G., & Chiosi, C. 1993, A&AS, 100, 647 [NASA ADS] [Google Scholar]
  15. Cao, Y., Kasliwal, M. M., Arcavi, I., et al. 2013, ApJ, 775, L7 [NASA ADS] [CrossRef] [Google Scholar]
  16. Crowther, P. A. 2007, ARA&A, 45, 177 [NASA ADS] [CrossRef] [Google Scholar]
  17. Crowther, P. A. 2013, MNRAS, 428, 1927 [NASA ADS] [CrossRef] [Google Scholar]
  18. Dopita, M. A., Kewley, L. J., Sutherland, R. S., & Nicholls, D. C. 2016, Ap&SS, 361, 61 [NASA ADS] [CrossRef] [Google Scholar]
  19. Eisenhauer, F., Abuter, R., Bickert, K., et al. 2003, Proc. SPIE, 4841, 1548 [Google Scholar]
  20. Eldridge, J. J., & Maund, J. R. 2016, MNRAS, 461, L117 [NASA ADS] [CrossRef] [Google Scholar]
  21. Eldridge, J. J., Fraser, M., Smartt, S. J., Maund, J. R., & Crockett, R. M. 2013, MNRAS, 436, 774 [NASA ADS] [CrossRef] [Google Scholar]
  22. Eldridge, J. J., Stanway, E. R., Xiao, L., et al. 2017, PASA, 34, e058 [NASA ADS] [CrossRef] [Google Scholar]
  23. Fagotto, F., Bressan, A., Bertelli, G., & Chiosi, C. 1994, A&AS, 105, 29 [NASA ADS] [Google Scholar]
  24. Filippenko, A. V. 1997, ARA&A, 35, 309 [NASA ADS] [CrossRef] [Google Scholar]
  25. Folatelli, G., Bersten, M. C., Benvenuto, O. G., et al. 2014, ApJ, 793, L22 [NASA ADS] [CrossRef] [Google Scholar]
  26. Folatelli, G., Bersten, M. C., Kuncarayakti, H., et al. 2015, ApJ, 811, 147 [NASA ADS] [CrossRef] [Google Scholar]
  27. Folatelli, G., Van Dyk, S. D., Kuncarayakti, H., et al. 2016, ApJ, 825, L22 [NASA ADS] [CrossRef] [Google Scholar]
  28. Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Fruchter, A. S., Levan, A. J., Strolger, L., et al. 2006, Nature, 441, 463 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  30. Galama, T. J., Vreeswijk, P. M., van Paradijs, J., et al. 1998, Nature, 395, 670 [NASA ADS] [CrossRef] [Google Scholar]
  31. Galbany, L., Stanishev, V., Mourão, A. M., et al. 2014, A&A, 572, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Galbany, L., Stanishev, V., Mourão, A. M., et al. 2016a, A&A, 591, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Galbany, L., Anderson, J. P., Rosales-Ortega, F. F., et al. 2016b, MNRAS, 455, 4087 [NASA ADS] [CrossRef] [Google Scholar]
  34. Galbany, L., Hamuy, M., Phillips, M. M., et al. 2016c, AJ, 151, 33 [NASA ADS] [CrossRef] [Google Scholar]
  35. Galbany, L., Mora, L., González-Gaitán, S., et al. 2017, MNRAS, 468, 628 [NASA ADS] [CrossRef] [Google Scholar]
  36. Gal-Yam, A. 2016, in Handbook of Supernovae, eds. A. W. Alsabti, & P. Murdin (Berlin: Springer) [Google Scholar]
  37. Gal-Yam, A., & Leonard, D. C. 2009, Nature, 458, 865 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  38. Gal-Yam, A., Leonard, D. C., Fox, D. B., et al. 2007, ApJ, 656, 372 [NASA ADS] [CrossRef] [Google Scholar]
  39. Georgy, C., Meynet, G., Walder, R., Folini, D., & Maeder, A. 2009, A&A, 502, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Graur, O., Bianco, F. B., Huang, S., et al. 2017a, ApJ, 837, 120 [NASA ADS] [CrossRef] [Google Scholar]
  41. Graur, O., Bianco, F. B., Modjaz, M., et al. 2017b, ApJ, 837, 121 [NASA ADS] [CrossRef] [Google Scholar]
  42. Habergham, S. M., Anderson, J. P., James, P. A., & Lyman, J. D. 2014, MNRAS, 441, 2230 [NASA ADS] [CrossRef] [Google Scholar]
  43. Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003, ApJ, 591, 288 [NASA ADS] [CrossRef] [Google Scholar]
  44. Hirai, R. 2017, MNRAS, 466, 3775 [NASA ADS] [CrossRef] [Google Scholar]
  45. Hoyle, F., & Fowler, W. A. 1960, ApJ, 132, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. James, P. A., & Anderson, J. P. 2006, A&A, 453, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Jerkstrand, A., Fransson, C., Maguire, K., et al. 2012, A&A, 546, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Kangas, T., Portinari, L., Mattila, S., et al. 2017, A&A, 597, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Kelly, P. L., & Kirshner, R. P. 2012, ApJ, 759, 107 [NASA ADS] [CrossRef] [Google Scholar]
  50. Kilpatrick, C. D., Foley, R. J., Abramson, L. E., et al. 2017, MNRAS, 465, 4650 [NASA ADS] [CrossRef] [Google Scholar]
  51. Krühler, T., Kuncarayakti, H., Schady, P., et al. 2017, A&A, 602, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Kuncarayakti, H., Doi, M., Aldering, G., et al. 2013a, AJ, 146, 30 [NASA ADS] [CrossRef] [Google Scholar]
  53. Kuncarayakti, H., Doi, M., Aldering, G., et al. 2013b, AJ, 146, 31 [NASA ADS] [CrossRef] [Google Scholar]
  54. Kuncarayakti, H., Maeda, K., Bersten, M. C., et al. 2015, A&A, 579, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Kuncarayakti, H., Maeda, K., Anderson, J. P., et al. 2016a, MNRAS, 458, 2063 [NASA ADS] [CrossRef] [Google Scholar]
  56. Kuncarayakti, H., Galbany, L., Anderson, J. P., Krühler, T., & Hamuy, M. 2016b, A&A, 593, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Kuncarayakti, H., Maeda, K., Ashall, C. J., et al. 2018, ApJ, 854, L14 [NASA ADS] [CrossRef] [Google Scholar]
  58. Langer, N. 2012, ARA&A, 50, 107 [NASA ADS] [CrossRef] [Google Scholar]
  59. Le Fèvre, O., Saisse, M., Mancini, D., et al. 2003, Proc. SPIE, 4841, 1670 [NASA ADS] [CrossRef] [Google Scholar]
  60. Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 [NASA ADS] [CrossRef] [Google Scholar]
  61. Leloudas, G., Gallazzi, A., Sollerman, J., et al. 2011, A&A, 530, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Leloudas, G., Schulze, S., Krühler, T., et al. 2015, MNRAS, 449, 917 [NASA ADS] [CrossRef] [Google Scholar]
  63. Levesque, E. M., Soderberg, A. M., Foley, R. J., et al. 2010a, ApJ, 709, L26 [NASA ADS] [CrossRef] [Google Scholar]
  64. Levesque, E. M., Berger, E., Kewley, L. J., & Bagley, M. M. 2010b, AJ, 139, 694 [NASA ADS] [CrossRef] [Google Scholar]
  65. Li, W., Leaman, J., Chornock, R., et al. 2011, MNRAS, 412, 1441 [NASA ADS] [CrossRef] [Google Scholar]
  66. Lyman, J. D., Bersier, D., James, P. A., et al. 2016, MNRAS, 457, 328 [NASA ADS] [CrossRef] [Google Scholar]
  67. Maeda, K., Hattori, T., Milisavljevic, D., et al. 2015, ApJ, 807, 35 [NASA ADS] [CrossRef] [Google Scholar]
  68. Margutti, R., Kamble, A., Milisavljevic, D., et al. 2017, ApJ, 835, 140 [NASA ADS] [CrossRef] [Google Scholar]
  69. Marino, R. A., Rosales-Ortega, F. F., Sánchez, S. F., et al. 2013, A&A, 559, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Matteucci, F., & Greggio, L. 1986, A&A, 154, 279 [NASA ADS] [Google Scholar]
  71. Maund, J. R. 2017, MNRAS, 469, 2202 [NASA ADS] [CrossRef] [Google Scholar]
  72. Maund, J. R., Smartt, S. J., Kudritzki, R. P., Podsiadlowski, P., & Gilmore, G. F. 2004, Nature, 427, 129 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  73. Maund, J. R., Arcavi, I., Ergon, M., et al. 2015, MNRAS, 454, 2580 [NASA ADS] [CrossRef] [Google Scholar]
  74. McCray, R., & Kafatos, M. 1987, ApJ, 317, 190 [NASA ADS] [CrossRef] [Google Scholar]
  75. Michałowski, M. J, Gentile, G., Krühler, T., et al. 2018, A&A, submitted [Google Scholar]
  76. Milisavljevic, D., Margutti, R., Kamble, A., et al. 2015, ApJ, 815, 120 [NASA ADS] [CrossRef] [Google Scholar]
  77. Modjaz, M., Kewley, L., Kirshner, R. P., et al. 2008, AJ, 135, 1136 [NASA ADS] [CrossRef] [Google Scholar]
  78. Modjaz, M., Kewley, L., Bloom, J. S., et al. 2011, ApJ, 731, L4 [NASA ADS] [CrossRef] [Google Scholar]
  79. Mould, J. R., Huchra, J. P., Freedman, W. L., et al. 2000, ApJ, 529, 786 [NASA ADS] [CrossRef] [Google Scholar]
  80. Nomoto, K., Tominaga, N., Umeda, H., Kobayashi, C., & Maeda, K. 2006, Nucl. Phys. A, 777, 424 [NASA ADS] [CrossRef] [Google Scholar]
  81. Ott, T. 2012, Astrophysics Source Code Library [record ascl:1210.019] [Google Scholar]
  82. Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 [NASA ADS] [CrossRef] [Google Scholar]
  83. Pignata, G., Stritzinger, M., Soderberg, A., et al. 2011, ApJ, 728, 14 [NASA ADS] [CrossRef] [Google Scholar]
  84. Podsiadlowski, P., Joss, P. C., & Hsu, J. J. L. 1992, ApJ, 391, 246 [NASA ADS] [CrossRef] [Google Scholar]
  85. Prieto, J. L., Stanek, K. Z., & Beacom, J. F. 2008, ApJ, 673, 999 [NASA ADS] [CrossRef] [Google Scholar]
  86. Rubin, A., Gal-Yam, A., De Cia, A., et al. 2016, ApJ, 820, 33 [NASA ADS] [CrossRef] [Google Scholar]
  87. Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444 [Google Scholar]
  88. Sanders, N. E., Soderberg, A. M., Levesque, E. M., et al. 2012, ApJ, 758, 132 [NASA ADS] [CrossRef] [Google Scholar]
  89. Sanders, N. E., Soderberg, A. M., Gezari, S., et al. 2015, ApJ, 799, 208 [NASA ADS] [CrossRef] [Google Scholar]
  90. Shivvers, I., Modjaz, M., Zheng, W., et al. 2017, PASP, 129, 054201 [NASA ADS] [CrossRef] [Google Scholar]
  91. Smartt, S. J. 2009, ARA&A, 47, 63 [NASA ADS] [CrossRef] [Google Scholar]
  92. Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 395, 1409 [NASA ADS] [CrossRef] [Google Scholar]
  93. Smith, N. 2014, ARA&A, 52, 487 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  94. Smith, N. 2017, in Handbook of Supermovae, eds. A. W. Alsabti, & P. Murdin (Cham: Springer International Publishing AG) [Google Scholar]
  95. Smith, N., & Arnett, W. D. 2014, ApJ, 785, 82 [NASA ADS] [CrossRef] [Google Scholar]
  96. Smith, N., & Tombleson, R. 2015, MNRAS, 447, 598 [NASA ADS] [CrossRef] [Google Scholar]
  97. Smith, N., Li, W., Filippenko, A. V., & Chornock, R. 2011a, MNRAS, 412, 1522 [NASA ADS] [CrossRef] [Google Scholar]
  98. Smith, N., Li, W., Miller, A. A., et al. 2011b, ApJ, 732, 63 [NASA ADS] [CrossRef] [Google Scholar]
  99. Soto, K. T., Lilly, S. J., Bacon, R., Richard, J., & Conseil, S. 2016, MNRAS, 458, 3210 [NASA ADS] [CrossRef] [Google Scholar]
  100. Sravan, N., Marchant, P., Kalogera, V., & Margutti, R. 2018, ApJ 852, L17, [NASA ADS] [CrossRef] [Google Scholar]
  101. Taddia, F., Sollerman, J., Razza, A., et al. 2013, A&A, 558, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Taddia, F., Sollerman, J., Fremling, C., et al. 2015, A&A, 580, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Taddia, F., Moquist, P., Sollerman, J., et al. 2016, A&A, 587, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Tartaglia, L., Sand, D., Valenti, S., et al. 2017, ATel, 10058, 58 [Google Scholar]
  105. Thornton, K., Gaudlitz, M., Janka, H.-T., & Steinmetz, M. 1998, ApJ, 500, 95 [NASA ADS] [CrossRef] [Google Scholar]
  106. Timmes, F. X., Woosley, S. E., & Weaver, T. A. 1995, ApJS, 98, 617 [NASA ADS] [CrossRef] [Google Scholar]
  107. Tully, R. B., Rizzi, L., Shaya, E. J., et al. 2009, AJ, 138, 323 [NASA ADS] [CrossRef] [Google Scholar]
  108. Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  109. Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 [NASA ADS] [CrossRef] [Google Scholar]
  110. Yoon, S.-C., Dessart, L., & Clocchiatti, A. 2017, ApJ, 840, 10 [NASA ADS] [CrossRef] [Google Scholar]
  111. Zapartas, E., de Mink, S. E., Izzard, R. G., et al. 2017, A&A, 601, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

1

Also known as the stripped-envelope (SE) SNe.

2

One SN (2017ahn, type II) was discovered much later after the IFU observations (Tartaglia et al. 2017). The explosion site of this SN was serendipitously covered in the IFU data, and therefore it was included in the sample.

3

SNe 1961I and 1964L were omitted as they do not satisfy the ≥1970 criterion.

4

ESO Programme ID 093.D-0318 and 094.D-0290.

5

ESO Programme ID 095.D-0172. The same dataset is also used in Krühler et al. (2017) and Michałowski et al. (2018).

7

Iraf is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation.

8

ESO Programme ID 089.D-0367, 091.D-0482, 093.D-0318, and 094.D-0290.

9

This late-time emission occurs for a number of interacting type IIn SNe where circumstellar material (CSM) interaction is still in progress and which results in an enduring light display. SN 1978K was found to be still bright 36 years after the explosion, as was discussed in a separate paper (Kuncarayakti et al. 2016a).

10

The [N II] line wavelength was incorrectly written as 6484 Å in the paper (Dopita et al., priv. comm.).

11

Anderson et al. (2015) found a similar result using direct explosion site metallicity measurement. However, see e.g. Graur et al. (2017a) and references therein, who include indirect metallicity estimates and argue for a lower number ratio at low metallicity.

All Tables

Table 1

Instrument configurations used in this study.

Table 2

SN sites observed with SINFONI.

All Figures

thumbnail Fig. 1

Relative frequency of each SN type in the current sample (large pie chart), compared to that of LOSS (Smith et al. 2011a, small pie chart in pale colours).

In the text
thumbnail Fig. 2

Histogram of SN host distance. The upper abscissa indicates the projected linear scale in pc/′′. The verticaldashed line indicates the median value.

In the text
thumbnail Fig. 3

Left: observed cumulative distributions of 12 + log(O/H) for different SN subtypes. Colour-coding is as follows: SNe Ic blue,Ib green, IIb orange, IIn purple, II red, and SESNe (Ic + Ib + IIb combined) dashed cyan. For comparison, the literature distributions of SN IcBL (dashed line), SLSN-I (dash-dotted line), and SLSN-II (dotted line) are also shown. Right: matrix for the K–S test result between different SN subtypes. The colour bar indicates the probability that samples were drawn from the same metallicity distribution, where the black end indicates significantly different distributions.

In the text
thumbnail Fig. 4

Evolution of Starburst99 age indicators, shown for Z = 0.02. The upper abscissa represents the stellar lifetime at the corresponding SSP age. Hα is indicated in red, Brγ green, CO 2.3 μ purple, and CaT blue. The CO EW values where typical SN II progenitors (MZAMS ≈ 8–16 M) are expected to be found are encompassed within the purple shaded area.

In the text
thumbnail Fig. 5

Left: observed cumulative distribution of log(HαEW) for different SN subtypes. SN sites which show no Hα emission are plotted as having log(HαEW) = 0. The dashed purple line indicates SN IIn distribution using the nearest H II regions (see text for description). Right: matrix for K–S test result. The colour-coding is the same as in Fig. 3.

In the text
thumbnail Fig. 6

Left and centre: observed cumulative distribution of parent stellar population age and progenitor initial mass for different SN subtypes. SN sites which show no Hα emission are plotted as having an age of 40 Myr, corresponding to the single-star SN progenitor lower mass limit of 8 M. Right: matrix for K–S test result. The colour-coding is the same as in Fig. 3.

In the text
thumbnail Fig. 7

Maps of the explosion site of SN 2009dq generated from the datacubes, in continuum-subtracted Hα (i), K-band (ii), Brγ (iii), andCO (iv). As the maps are color-coded with white for high-count pixels and black for low-count pixels, the pixels with strong CO absorption appear black in panel (iv). The white 3′′   × 3′′ square in panel (i) represents the SINFONI FoV, and the white circle denotes the SN position within 1′′ radius. North is up and east is left in all maps. Cyan circles in panels (iii) and (iv) indicated with A and B have identical radii of 0.25′′. Panel (v) shows the SINFONI spaxels plotted according to the Brγ and CO NCRs; these are shown as contours at 10%, 50%, and 90% levels, while for comparison, the VIMOS spaxels are plotted in red plus marks according to the Hα and [N II] NCR. Panel (vi) shows the extracted spectra of regions A and B, in the Brγ and CO spectral regions. The two lines are indicated with green-shaded regions. Wavelength is in the observer frame.

In the text
thumbnail Fig. 8

Cumulative distribution of host H II region flux and luminosity.

In the text
thumbnail Fig. 9

Number ratios of different SN types in metallicity bins: N(SESN)/N(SN II) (upper panel), and N(Ic)/N(Ib) (lower panel). Horizontal error bars indicate the bin size and vertical error bars are Poisson errors.

In the text
thumbnail Fig. 10

Diagram showing the plane of progenitor star initial mass and metallicity. Shaded areas are stellar evolution predictions from rotating single-star models of Georgy et al. (2009), for SN II (pink), Ib (light green), and Ic (blue-green). Data points from this work are colour-coded: blue for SN Ic, bright green for Ib, red for II, purple for IIn, orange for IIb, and cyan forIcBL. Uncertainties in metallicity reflect the 0.16 dex error in the 12 + log(O/H) N2 calibration of Marino et al. (2013).

In the text

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

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

Initial download of the metrics may take a while.