© ESO, 2017

1. Introduction

Our understanding of the late evolutionary stages and final fates of massive stars (with masses above ~8 M_⊙) remains incomplete, despite the decades of theoretical and observational work devoted to it. Generally, these stars are expected to end their lives as core-collapse supernovae (CCSNe) after evolving through a series of stages of nuclear burning until they have built up an iron core. The type of the resulting SN explosion is expected to depend on the initial mass of the progenitor star, along with an as yet unclear contribution from a number of other factors such as the metallicity, rotation and multiplicity of the star (e.g. Smartt 2009; Smith et al. 2011a; Georgy et al. 2012; Eldridge et al. 2013, and references therein).

CCSNe are divided into two main types, hydrogen-poor type Ib/c SNe (type Ic also being helium-poor) and hydrogen-rich type II SNe (Filippenko 1997). Type Ib and Ic SNe are believed to have progenitors that have lost their hydrogen and/or helium envelopes, while the progenitors of type II SNe still retain at least a part of their envelope, resulting in hydrogen lines in their spectra (e.g. Smartt 2009). Type II SNe are also divided into subtypes. The red supergiant (RSG) progenitors of type II-P (plateau) SNe still retain massive hydrogen envelopes that power the distinctive plateaus in their light curves (Grassberg et al. 1971), while type II-L (linear) progenitors are believed to have lost a significant part of their hydrogen envelope and type IIb progenitors almost all of it. Although originally considered two distinct populations of events (e.g. Barbon et al. 1979), types II-P and II-L are now more commonly thought to occupy a continuum of different levels of mass loss primarily influenced by the initial mass of the progenitor (e.g. Anderson et al. 2014; Sanders et al. 2015; González-Gaitán et al. 2015). Type IIn (narrow lines) SNe are characterized by the presence of a dense circumstellar medium (CSM) at the time of explosion, resulting in strong interaction between the SN ejecta and the CSM (e.g. Schlegel 1990; Stathakis & Sadler 1991). Non-terminal outbursts of massive stars in external galaxies, analogous to the Great Eruption of η Car, have sometimes been mistaken for true SNe, and are called SN impostors (e.g. van Dyk et al. 2000; Pastorello et al. 2013; Kankare et al. 2015).

For type Ib and Ic SNe, both very massive single Wolf-Rayet (WR) stars (e.g. Maeder & Lequeux 1982) and lower-mass stars stripped through interaction with a binary companion (e.g. Podsiadlowski et al. 1992) have been suggested as progenitors, and Eldridge et al. (2013) suggested a mix of both progenitor channels. Recently, Smith (2014) challenged the established picture of single-star mass loss; namely, that mass loss through conventional line-driven winds should not be sufficient to create a type Ib/c SN, and most of their progenitors should be stripped through binary interaction. Lyman et al. (2016) found the ejecta masses of type Ib/c and IIb SNe inconsistent with very massive stars, also indicating that interacting binaries are the dominant progenitor channel. The effects of multiplicity on massive star evolution, however, are not well understood either, despite its tremendous importance (e.g. Sana et al. 2012, 2013). Properly taking into account Roche-lobe overflow, mergers and other features of binaries is thus one of the most important challenges in understanding the evolution of massive stars.

Empirical studies of the progenitors of different types of SNe can help distinguish between different evolutionary schemes. Direct detections of type II-P SN progenitors in pre-explosion high-resolution images (e.g. Mattila et al. 2008; van Dyk et al. 2012; Maund et al. 2013; Fraser 2016) have already helped to establish them as RSGs with initial masses of ~8.5 to ~16.5 M_⊙ (Smartt 2009, 2015).

The detected type IIb progenitors are consistent with interacting binary systems (e.g. Aldering et al. 1994; Maund et al. 2004; van Dyk et al. 2011; Folatelli et al. 2014; 2015; Fox et al. 2014; Maeda et al. 2014), and a WR progenitor was ruled out for SN 2008ax (Folatelli et al. 2015) although the proposed binary companion of the SN 2011dh progenitor was disputed by Maund et al. (2015a). Maeda et al. (2015) found that the more extended progenitors may still be undergoing binary interaction, while for the less extended ones this phase would be over. The detected type IIb progenitors have been identified as yellow or blue supergiants (YSGs or BSGs, respectively) with initial masses between 13 and 19 M_⊙ (van Dyk et al. 2011; Maund et al. 2011; van Dyk et al. 2014; Folatelli et al. 2015).

Among other SN types, however, direct progenitor detections are scarce. A few type IIn events have been connected to a luminous blue variable (LBV) progenitor (e.g. Gal-Yam & Leonard 2009; Smith et al. 2011b; Fransson et al. 2014), along with events such as SN 2009ip and SN 1961V where the disappearance of the progenitor has not been firmly ascertained (e.g. van Dyk & Matheson 2012; Fraser et al. 2015). The (probably) only detected progenitor of a type II-L SN, that of SN 2009hd, was likely a high-mass RSG (Elias-Rosa et al. 2011). A possible progenitor for another type II-L, SN 2009kr (Fraser et al. 2010), was shown by Maund et al. (2015b) to probably be a small cluster instead of a single star.

No type Ic progenitors have been detected so far, possibly because of their faintness in the optical bands (Yoon et al. 2012). The only detected type Ib progenitor, that of iPTF13bvn (Cao et al. 2013), was suggested to be an initially 10−20 M_⊙ star in an interacting binary by Eldridge et al. (2015). More recently Eldridge & Maund (2016) suggested an initially 10−12 M_⊙ star that evolved into a helium giant; Folatelli et al. (2016) also favoured a binary progenitor, although none of the models they examined proved entirely satisfactory. van Dyk et al. (2016) excluded a ≳10 M_⊙ companion to the progenitor of the type Ic SN 1994I, which earlier had been considered the result of binary interaction (Nomoto et al. 1994). Heikkilä et al. (2016) used pre-explosion X-ray observations to investigate the possibility of high-mass X-ray binaries (in which one component is a neutron star or a stellar-mass black hole) being the progenitors of some type Ib, Ic or IIb SNe. Such progenitors were found to be rare.

In lieu of direct detections, SN progenitors can also be studied in various indirect ways. Hakobyan et al. (2009, 2016), Habergham, Anderson & James (2010) and Habergham, James & Anderson (2012), for example, examined the distances of different types of SNe from the nuclei of their host galaxies, and found type Ib/c SNe to be more centrally concentrated than type II. Kelly & Kirshner (2012) used the colours of SN environments to suggest a higher initial mass for type Ib/c progenitors than those of type II. Kuncarayakti et al. (2013a) used integral field spectroscopy to infer a higher environmental metallicity and initial mass for type Ic progenitors than type Ib, while Kuncarayakti et al. (2013b) found some type II progenitors to have initial masses comparable to type Ib/c progenitors. Leloudas et al. (2010) connected type Ib/c SNe and gamma-ray bursts to WR stars using the relative B- and K-band brightness of their environments. Spectral synthesis model fitting to spectra of SNe at the nebular phase has been used by Mazzali et al. (2010) and Jerkstrand et al. (2015) to suggest an initially 15 M_⊙ progenitor for type Ic SN 2007gr and a type IIb progenitor mass range of 12−16 M_⊙, respectively. Tomasella et al. (2013) used hydrodynamical modelling to constrain the ejecta mass of a type II-P event, SN 2012A, to 12.5 M_⊙. Dall’Ora et al. (2014) suggested an envelope mass of ~ 20M_⊙ for SN 2012aw, another type II-P event, using similar methods, challenging the ~8.5 to ~16.5 M_⊙ initial mass range.

Because ionizing radiation from young massive stars is responsible for creating H ii regions (Kennicutt 1998), correlations between different SN types and the Hα emission from their host galaxies can also be used to statistically study their progenitors. James & Anderson (2006, hereafter JA06), used a method called pixel statistics to study this correlation, and this method was used again with larger samples by Anderson & James (2008) and Anderson et al. (2012), hereafter AJ08 and A12, respectively. These results indicated a higher average initial progenitor mass for type Ic SNe than types Ib, II-P or IIn. Habergham et al. (2014, hereafter H14), applied the same method to interacting transients, that is type IIn SNe and SN impostors. Kangas et al. (2013, hereafter K13), used both the pixel statistics method and the distances to host galaxy nuclei to study SNe in strongly star-forming galaxies specifically, and also found the A12 mass sequence (with a stronger correlation between type Ic and Hα emission than in normal galaxies) and an increased centralization for type Ib/c. The results from Hα images were also supported with near ultraviolet (NUV) pixel statistics by both A12 and K13. However, the pixel statistics method has so far remained qualitative. Crowther (2013) argued that it can only provide weak constraints because the lifetimes of the giant H ii regions probed by A12 are longer than the lifetimes of some CCSN progenitors. The Crowther (2013) method of measuring distances from the SNe to the nearest H ii regions, however, corroborated the A12 result. Indirect methods are also affected by biases such as the Malmquist bias when it comes to constructing a sample of SNe.

Fig. 1 Original continuum-subtracted Hα intensity maps used in this study: the LMC on the left (from the SHASSA survey; Gaustad et al. 2001) and M 33 on the right (Hoopes & Walterbos 2000), cropped in the case of the LMC but otherwise unaltered. North is up and east is left. The black scale bar corresponds to 2 deg (~1.7 kpc) in the LMC and 20 arcmin (~4.9 kpc) in M 33.

The results of A12 showed a weaker correlation with Hα emission for type IIn than other CCSNe, indicating a lower initial progenitor mass. A12 argued that while the diversity within the subtype allows for some progenitors to be LBVs, most IIn progenitors are instead lower-mass stars, such as asymptotic giant branch (AGB) stars exploding as electron-capture SNe. To reconcile this with the identified massive LBV progenitors of type IIn SNe, Smith & Tombleson (2015) examined the LBVs in the Large Magellanic Cloud (LMC) and found them to be relatively isolated from O-type stars, which was argued to at least partially explain the weak correlation with Hα emission for the SNe. However, Humphreys et al. (2016) pointed out that this LBV sample included both classical and lower-luminosity post-RSG LBVs along with some unrelated stars, and that the classical LBVs are in fact closely associated with the O stars. A subset of ambiguous type IIn/Ia-CSM SNe also exists. It is unclear whether these events are thermonuclear (Fox et al. 2015) or core-collapse SNe (Inserra et al. 2016), which further complicates the already diverse type IIn. A fraction of type Ia SNe disguised as type IIn may partially explain the weak correlation with Hα emission.

In this paper, we used Hα intensity maps of two nearby galaxies, the LMC and Messier 33 (M 33), to study the correlation between massive stars and star-forming regions in a way similar to A12 and K13. Nearby galaxies like these provide a unique opportunity for this study, as they are close enough for individual bright stars to be resolved and, unlike the Milky Way, we can see practically all of the Hα-emitting regions in them from the outside. Two very different galaxies were included in order to check the consistency of the results. We then compared these results to those of different CCSN types in other galaxies in order to test the validity of the pixel statistics method and attempt to derive more quantitative constraints for CCSN progenitors. We describe our Hα images and the catalogues we used for stellar coordinates in Sect. 2. In Sect. 3 we describe our analysis methods. The results of the analysis are presented in Sect. 4, with an investigation of possible systematic effects in Sect. 5. Our findings are discussed in Sect. 6. Finally, in Sect. 7 we present our conclusions.

2. Data

Here we describe the Hα images, catalogues of stellar coordinates and parameters, and SN samples used in the pixel statistics analysis.

2.1. Hα images

We used two local galaxies to study the correlations between stars and H ii regions: the LMC, an SB(s)m-type dwarf galaxy at a well-established distance of ~50 kpc, and M 33, an SA(s)cd-type spiral galaxy at a distance between ~700 and ~900 kpc (we adopted the most recent Cepheid distance, 839 kpc, from Gieren et al. 2013). The Hα images of these galaxies are shown in Fig. 1, with no alterations except cropping in the case of the LMC.

The continuum-subtracted Hα image of the LMC was obtained from the Southern Hemisphere All-Sky Survey Atlas¹ (SHASSA; Gaustad et al. 2001). The observations were performed using the Swarthmore robotic camera at the Cerro Tololo Inter-American Observatory (CTIO) in Chile. Details on the observations and reduction procedures are in the survey paper by Gaustad et al. (2001). The uncropped image covers an area of 11.8 deg × 12.4 deg, which was cropped for our analysis into a size of 7.9 deg × 8.2 deg to only show the areas with Hα emission. The pixel scale in the image is 47.64 arcsec pix^-1, corresponding to 11.6 pc pix^-1 at the distance of 50 kpc. The image has a spatial resolution (full width half maximum, FWHM) of 4.0 arcmin because of a median filter smoothing algorithm used in the survey for easier removal of foreground stars. At 50 kpc this corresponds to a resolution of 58.2 pc, which is sufficient for this study.

The continuum-subtracted Hα image of M 33 was obtained from the NASA/IPAC Extragalactic Database (NED)². This mosaic was originally made by Hoopes & Walterbos (2000) using images from the 0.6-m Burrell-Schmidt telescope at Kitt Peak National Observatory, taken in 1995. The image on NED covers an area of 1.15 deg × 1.15 deg. The image has a 2.03 arcsec pix^-1 scale and a spatial resolution (FWHM) of 3.0 arcsec, corresponding to 8.3 pc pix^-1 and 12.2 pc, respectively, at the distance of 839 kpc. For details on the observations and reductions, see Hoopes & Walterbos (2000).

2.2. Stellar samples

Fig. 2 Hα images of the LMC (left) and M 33 (right), with our RSG, YSG, WN and WC samples overplotted in order to illustrate the spatial coverage of their respective catalogues. The black scale bar corresponds to 2 deg (~1.7 kpc) in the LMC and 20 arcmin (~4.9 kpc) in M 33.

Fig. 3 Hα images of the LMC, with our main-sequence B and O star samples as well as the SG B[e] stars and LBVs overplotted. A concentration of B stars is visible in a rectangular area (marked in red) with no strong Hα sources, indicating a spatial bias in this subsample (see also Sect. 6). The black scale bar corresponds to 2 deg (~1.7 kpc).

The coordinates and parameters (such as temperature, spectral type and/or luminosity) of the stellar samples in the two galaxies were obtained from various sources, described below. Because of the nature of our analysis (Sect. 3), the aim in each case was to obtain a list of coordinates with good spatial coverage and to avoid biases against regions of high or low Hα surface brightness. The individual catalogue papers contain details on the observations, reductions and the selection and confirmation of candidates in each catalogue; these will be briefly summarized below. The spatial coverage of the supergiant and WR catalogues in both galaxies is illustrated in Fig. 2, while the coverage of the main-sequence B and O stars, as well as SG B[e] stars and LBVs in the LMC, is shown in Fig. 3.

2.2.1. Supergiants

A catalogue of spectroscopically classified RSGs and YSGs in M 33 was compiled by Drout et al. (2012). The observations were performed using the Hectospec multi-fibre spectrograph on the 6.5-m MMT telescope at the Fred Lawrence Whipple Observatory in Arizona, US. The RSG and YSG candidates in M 33 were photometrically selected by Drout et al. from the Local Group Galaxy Survey (LGGS; Massey et al. 2006), applying the following criteria: the YSG candidates had V < 18.7 mag, U−B > −0.4 mag and 0 mag ≤ B−V ≤ 1.4 mag, while for the RSG candidates V < 20.0 mag and V−R > 0.6 mag. An additional criterion of B−V > −1.599(V−R)² + 4.18(V−R)−0.83 mag (Massey et al. 2009) was adopted to separate likely foreground dwarfs from the supergiants. The rest of the foreground stars were excluded spectroscopically using the radial velocities of each star. Using the upper and lower temperature limits of YSGs set by Drout et al. (2012) at 7500 K (the line between the A and F spectral classes) and 4800 K (the line between the G and K spectral classes), respectively, the catalogue contains 188 RSGs, 74 YSGs and 47 BSGs. Spatially, the catalogue covers the entire disk of M 33, but only two RSGs in the central square kpc region are included because of crowding effects. As this region is quite small and does not in any way dominate the Hα emission of M 33 (Fig. 1), we consider any biases in the spatial coverage with respect to Hα surface brightness insignificant (this assumption is tested and found valid in Sect. 5). Luminosity-wise, the catalogue was argued by its authors to be complete down to a luminosity of log L/L_⊙ ~ 4.8, and contains 154 stars below this limit as well.

A similar catalogue for the LMC was compiled by Neugent et al. (2012) by observing 64 fields within the galaxy with the 138-fibre multi-object spectrometer Hydra on the Cerro Tololo 4-m telescope. YSG and RSG candidates were selected by Neugent et al. from the USNO CCD Astrograph Catalogue Part 3 (UCAC3) using magnitudes from the Two Micron All Sky Survey³ (2MASS; Skrutskie 2006) and spectroscopically classified. For RSG candidates (J−K)_2MASS < 0.9 mag and K_2MASS ≤ 10.2 mag. For YSG candidates the K brightness limit depends on (J−K) as described in Neugent et al. (2010) – at (J−K)_2MASS = 0 mag the limit is K_2MASS = 12.6 mag and at (J−K)_2MASS = 0.9 mag it is K_2MASS = 10.2 mag. Foreground red and yellow dwarf stars were excluded using radial velocities. Eight stars were classified merely as “possible LMC supergiants” (as opposed to “probable”) and were excluded from the present study. Using the same limiting temperatures as for M 33, the catalogue contains 543 RSGs, 109 YSGs and 163 BSGs. The spatial coverage of this catalogue, as demonstrated in Neugent et al. (2012), is not perfect over the entire galaxy but the 64 fields are spread roughly evenly over the galaxy and thus the catalogue should not be biased toward regions of any particular Hα surface brightness. Such coverage is sufficient for this study (this assumption is also tested in Sect. 5). The catalogue was deemed by its authors to be missing no more than a few YSGs, while the RSG completeness is estimated to be roughly 30 per cent (out of some ~1800 RSGs).

As the BSGs in the catalogues do not include stars with temperatures higher than ~ 11 000 K, we did not include the BSGs in our analysis.

2.2.2. WR stars

We also obtained catalogues of spectroscopically classified WR stars in both galaxies. WR stars with CNO cycle products, especially nitrogen, on their surfaces, are classified WN, while those with triple-α products, especially carbon, on the surface are classified WC. WC stars are considered a later stage of evolution than WN, and stars that reach the WC stage before their deaths have a higher minimum initial mass (~25 and ~40 M_⊙ for WN and WC stars, respectively; Crowther 2007). There is, however, also evidence (Georgy et al. 2012; McClelland & Eldridge 2016) that some WC stars require a different evolutionary channel than WN stars and may have initial masses lower than what Crowther (2007) concluded. WN stars with strong hydrogen signatures are classified WNh (Smith et al. 1996). WNh stars are still considered to be core-hydrogen-burning, and their initial masses are believed to be higher than those of WN or WC stars (≳60 M_⊙; Crowther 2007; Smith & Conti 2008).

For the LMC, we used the catalogue of WN stars by Hainich et al. (2014), considered by the catalogue authors to be complete (or very close to complete) over the entire galaxy. This catalogue is mostly based on the earlier BAT99 catalogue (Breysacher et al. 1999), with some more recent re-classifications. This catalogue contains a total of 94 WN stars. For the WC stars, we used the massive star catalogue of Bonanos et al. (2009); the 24 WC stars, 22 of type WC4 and two labeled early-type WC, in this catalogue are also largely based on BAT99 and thought to represent a complete sample. For M 33, we used the catalogue of WR stars by Neugent & Massey (2011). This catalogue is considered by its authors to be ~ 95 per cent complete, and contains 139 WN stars, 52 WC stars and ten Ofpe/WN9 stars. Stars with an uncertain WR status (classification for example “WN?”) were excluded – there is one of these in the LMC and three in M 33.

2.2.3. Other evolved stars

We also used the LBV stars listed by Smith & Tombleson (2015). These include all confirmed LBVs in the LMC and likely, well-studied LBV candidates with shells. We excluded the five unrelated stars and the duplicate entry pointed out by Humphreys et al. (2016), bringing the total number to ten LBVs – three classical LBVs and seven lower-luminosity LBVs. In M 33, there are only four confirmed LBVs (Hubble-Sandage variables B, C and 2 and Var 83; Humphreys et al. 2014), but we included them as well (the actual number of LBV stars in M 33 was estimated to be in the hundreds by Massey et al. 2007). As Smith & Tombleson (2015) considered supergiant B[e] (SG B[e]) stars to be the lower-mass analogs of LBVs, these stars in the LMC were included for comparison. The coordinates of 12 spectroscopically classified SG B[e] stars were obtained from Zickgraf (2006) and Bonanos et al. (2009); this includes every such star discovered in the LMC.

2.2.4. Main-sequence stars

As well as these evolved stars, we analyzed massive main-sequence stars in the LMC with spectroscopic classification. Their coordinates were obtained from the SIMBAD Astronomical Database⁴ (Wenger et al. 2000), with a search radius of five degrees from the centre of the LMC as reported in NED. With such a search we obtained the coordinates of main-sequence stars between the B2V and O3V spectral types as reported in April 2015. Very few stars classified as main-sequence spectral types later than B2V could be found in the LMC (for example, only five stars of type B3V as opposed to 92 of type B2V). The main-sequence subsamples are most likely spatially biased, and were only used for qualitative purposes. The Bonanos et al. (2009) catalogue contains hundreds of main-sequence stars, but being a compilation of catalogues that mostly target specific regions in the LMC, as a main-sequence catalogue it is also spatially biased. As the SIMBAD search results simply contain higher numbers of stars, we used them instead.

2.3. SN samples

To compare the pixel statistics (Sect. 3) of the stars to previous statistics of CCSNe, we used the distances and results of the sample of SNe reported in A12 which include the results of AJ08 as well. The host galaxies of the SNe in AJ08 and A12 were observed using the 1-m Jacobus Kapteyn Telescope (JKT), the 2-m Liverpool Telescope (LT) and the 2.5-m Isaac Newton Telescope (INT) at the Observatorio del Roque de los Muchachos on La Palma, Canary Islands, and with the European Southern Observatory (ESO) 2.2-m telescope at La Silla, Chile. In some cases we also made use of the results of K13. The Hα images in K13 were observed with the Andalucia Faint Object Spectrograph and Camera (ALFOSC) on the 2.5-m Nordic Optical Telescope (NOT) on La Palma. Details of these observations can be found in the aforementioned papers.

Samples of SNe are potentially affected by issues in precise SN classification. Type Ib, Ic and IIb SNe can be difficult to distinguish spectroscopically (Filippenko 1997). Although type II-L SNe generally exhibit a weaker P Cygni absorption profile at the Hα line than type II-P SNe (Schlegel 1996), separating these types can be difficult without light curves which are not always available – not to mention that, for example, Anderson et al. (2014) argue that there is no unambiguous borderline between types II-P and II-L at all. The A12 SN subsamples may therefore be affected by mixing between the subtypes (for example, one type Ib/IIb SN was included as half a type Ib and half a type IIb). Type IIn SNe, furthermore, are a very diverse class with several possible progenitor channels, with narrow hydrogen emission lines being the common feature. Therefore the sample of type IIn SNe may not include representative fractions of different SNe within this type. These effects may result in, for example, larger scatter in the correlations between SNe and Hα emission.

3. Analysis

Here we describe the methods we used to connect CCSNe to possible progenitor populations.

3.1. The pixel statistics method

The pixel statistics analysis method was introduced by Fruchter et al. (2006), who applied it to blue light to study the long gamma-ray burst distribution, and JA06. It is described in detail by JA06, and has since been used to study correlations between CCSNe and emission at different wavelengths, mainly Hα, by AJ08, A12, K13 and H14.

Briefly put, each pixel in the Hα image of a galaxy is assigned a number, called its normalized cumulative rank (NCR), that tells us what fraction of the Hα emission originates in pixels with lower surface brightness. When the coordinates of a SN or, in the case of this study, a star, are connected to a pixel and thus an NCR, this acts as an indicator of any correlation between the SN or star and the underlying Hα emission. Each “background” pixel with no emission has an NCR of zero, while the brightest pixel of the galaxy has an NCR of one. This requires that the bias level and sky emission have been properly subtracted. When NCRs are calculated for a sample of objects, the mean NCR and distribution of the NCRs can be used to statistically compare the samples, for example, with an Anderson-Darling (AD) test. A distribution that exactly follows the underlying emission has a mean value of 0.5 and a uniform distribution. As an example, 70 per cent of the sample have NCR ≤ 0.7 since 70 per cent of the emission comes from pixels with NCR ≤ 0.7.

In earlier studies, starting with JA06, a higher NCR(Hα) for a set of SNe has been interpreted as a shorter lifespan of the SN progenitor and thus an initially more massive progenitor star. The studies of, for example, A12 and K13 probed structures at a projected scale of roughly 100 or 200 pc depending on the distance of the galaxy. Any correlation with Hα emission thus does not indicate a correlation with small H ii regions created by individual massive stars but with larger H ii complexes. Crowther (2013) argued that the ages of these complexes are measured in tens of Myr, as opposed to smaller H ii regions which only last as long as their central star continues to emit ionizing radiation. The lifespan of a large H ii complex is thus comparable to that of a star with initial mass around 10 M_⊙ and the same complex can host more than one generation of more massive stars.

Nonetheless, as a star grows older, a few things happen to lower its NCR. The star gradually drifts away from its native H ii complex, and the complex may grow fainter or dissipate entirely as massive stars within it explode as SNe, sweeping the gas away, and stop emitting ionizing radiation. Therefore even though a relatively low-mass SN progenitor may still be inside its H ii complex when it explodes, it is still less likely to be there. Some such stars have a high NCR, being born early in the life of the complex or having drifted into another H ii region. Others, born at later times, have a lower (possibly zero) NCR as the complex disappears around them – the timing depends on the lifespan of the star. On the other hand, the higher-mass stars have the shortest lifetimes and the highest fluxes of ionizing radiation, and stay closer to the centre of the H ii region.

Some stars are born into smaller H ii regions than others. Because complexes of different size and brightness have different NCRs, and some smaller and fainter regions are obscured by noise in the Hα image, even high-mass stars occasionally have low NCRs. The smallest H ii regions also have the shortest lifespans, which may further lower the NCRs of stars born in them although, on the spatial scales probed by A12 and K13, the small H ii regions either blend into the large ones or do not have a high NCR to begin with. Conversely, projection effects may result in higher NCRs for stars lying on the line of sight to an H ii complex. However, statistically, the longer the star lives the less likely it is to have a high NCR. Therefore, samples of stars or SNe should still form a sequence of descending initial mass with descending NCR. Testing this assumption is one of the aims of this paper.

Some [N ii] λλ6548,6584 emission is included in the Hα filters used in the SHASSA survey (Gaustad et al. 2001) and by Hoopes & Walterbos (2000). The level of contamination differs with the metallicity of the emitting region, and metallicity gradients in galaxies could bias NCR distributions. However, using Eq. (6) in Helmboldt et al. (2004), we estimated that the [N ii]/Hα flux ratios in the LMC and M 33 should be about 0.14 and 0.2, respectively. Furthermore, because of the transmission curve of the filter used by Gaustad et al. (2001), only a third of the [N ii] emission would be observed, while the peak transmission (at the Hα line) was 78 per cent. For M 33, Hoopes & Walterbos (2000) estimated the [N ii] contribution in their study to be five per cent or less of the total flux. Therefore, metallicity differences between the two galaxies or between regions inside one galaxy should not significantly affect observed fluxes in the Hα filters or the resulting NCRs. Previous NCR studies of SNe, which we used for comparison, were similarly contaminated. The metallicity difference between galaxies is not a concern because NCRs are by definition normalized to the brightest region in a galaxy. Metallicity-related NCR biases should be mainly due to internal metallicity gradients (or their absence) inside galaxies. Assuming that the [N ii] lines mainly fall on the wings of the Hα filter, we conclude that the effect on our NCR comparisons should be minimal. Line strengths in H ii regions are also affected by other factors, such as differences in dust content, ionization, temperature and density in the regions, so that Hα emission does not perfectly trace star formation. However, estimating the biases these effects may cause in the NCRs is outside the scope of this paper. For now, we neglected them for simplicity, consistently with other NCR studies.

With this in mind, we calculated the NCRs of each star in each of our samples. For comparison, we also generated a sample of 250 random positions within each of the two galaxies, with a uniform spatial distribution inside the visible extent of the galaxy. For the LMC this corresponds to a circle with a radius of 3.5 deg, and for M 33 an ellipse with a = 33 arcmin, b = 20 arcmin and position angle 22.7 deg. The distribution of the random positions is shown in Fig. 4. These steps were repeated for each of the images we analyzed in this study, as described below.

Fig. 4 Hα images of the LMC (left) and M 33 (right), with the coordinates of our randomly generated uniformly distributed samples of positions overplotted.

3.2. Comparison to CCSN pixel statistics

In order to make our NCR results comparable to those of the SNe in A12, we simulated observations of the two galaxies at a distance, spatial resolution and signal-to-noise level comparable to the galaxies in A12. All the galaxy images in this sample were re-binned to a pixel scale of ~1 arcsec pix^-1. The median distances of the type II-P and II-L SN host galaxies in A12 are close to 20 Mpc and those of type Ib, Ic, IIn and IIb SNe are all close to 35 Mpc. Therefore we simulated distances of 20 and 35 Mpc, and using these we obtained our main results. The images were first convolved to have a spatial resolution (FWHM) close to what 1 arcsec would correspond to at distances of 20 and 35 Mpc, that is 97 and 170 pc, respectively. This resolution is comparable to the sizes of giant H ii regions (Crowther 2013). These are equivalent to 8.4 and 14.7 pix in the original LMC image and 11.8 and 20.6 pix in the M 33 image. The images were then re-binned to have a pixel scale as close as possible to 1 arcsec pix^-1 at this simulated distance. This was done using the iraf tasks gauss and blkavg. In addition, using the task mknoise in the package artdata, Poisson noise was added to the images in order to match the signal-to-noise ratio (S/N) per pixel with the images used by A12. The S/N per pixel of the images used by A12 was on the order of 50 to 100 in H ii regions of average brightness and on the order of 10 to 20 in faint ones.

To make the spatial scales of the Hα images comparable to the images in A12, the LMC image was re-binned by eight to simulate a distance of 20 Mpc and by 15 for 35 Mpc, while the M 33 image was re-binned by 12 and by 21 for these distances. In the original image, this re-binning would correspond to a pixel scale of over 6 arcmin pix^-1 for the LMC and over 24 arcsec pix^-1 for M 33. The possible uncertainties in the coordinates of the stars in these galaxies, or in the coordinate systems of the images themselves (on the order of arcseconds or less depending on the galaxy), are insignificant compared to this pixel size. Therefore the positions and NCRs of the stars in our samples are sufficiently accurate. However, the uncertainties in SN positions and in the coordinate systems of the images of their host galaxies are up to ~1 arcsec (as was conservatively estimated by K13), which is comparable to the pixel size of the images of A12. Therefore the NCRs of the SNe and local galaxy stars should not be directly compared. To compensate for this, we simulated how an accurate NCR distribution would change with the addition of positional uncertainties by performing the following Monte Carlo analysis. After calculating the NCRs of each stellar sample using the original, accurate coordinates, an offset with a random direction and distance was added to each coordinate set. The offset distances were drawn from a Gaussian probability distribution with a standard deviation of σ = 0.5 arcsec. This accounts for the errors in the world coordinate system fitting typical for the images in K13 (0.2 arcsec) and uncertainties in the coordinates reported for SNe. Thus 95 per cent of the positional shifts are 2σ = 1 arcsec or less. This was repeated 1000 times for each stellar sample and a new mean NCR was then calculated. The errors of the NCRs with accurate coordinates are reported as standard errors of the mean (SEM), while the errors reported for the NCRs with positional errors are , where σ_MC is the standard deviation of the Monte Carlo NCRs for the sample and N the original size of the stellar sample. This is analogous to the SEM in order to be comparable to the SEMs reported by A12 and K13. The NCR distributions with offsets could then finally be compared to those of the SNe. Both the accurate NCRs and those calculated with offsets are listed in Sect. 4.

The distributions of NCRs that include positional errors should not be compared to others with an AD test using the Monte Carlo distribution as it is – despite the Monte Carlo distribution having 1000 times more individual values, the analysis does not increase information. Therefore, before any AD tests with these distributions were made, we re-binned them by sorting them in order of increasing NCR and taking the average of every 1000 values.

4. Results

We divided our samples from catalogues and SIMBAD (described in Sect. 2) into sub-samples as follows. Each of the spectral types of the LMC main-sequence stars is its own sub-sample. The Wolf-Rayet samples were split into WN and WC stars in both galaxies. The WN stars were split into early- (WN2-WN5) and late-type (WN6+) stars. In the LMC we analyzed samples of WN stars including WNh stars, simply labeled “WN”, and without WNh stars, labeled “WN (no H)”. While WN and WNh stars are expected to be systematically different (Crowther 2007; Smith & Conti 2008), this approach eases the comparison between the two galaxies. In M 33, the Neugent & Massey (2011) catalogue only contains one WNh star. However, Hamann, Gräfener & Liermann (2006) listed many WNh stars in the Milky Way, providing circumstantial evidence that they should also exist in larger numbers in M 33; one star in the Neugent & Massey (2011) catalogue is labeled H-rich but not classified WNh; one of the classification sources (Abbott et al. 2004) intentionally omitted the Smith et al. (1996) WNh criteria; and one WN7h star from another source (Drissen et al. 2008) is simply reported as a WN7 in Neugent & Massey (2011). Therefore, in M 33, WN and WNh stars seem to be mixed and cannot be reliably separated. WC stars were also split into early- and late-type; however, in the LMC, the WC stars are all early-type. In M 33, we also analyzed the ten stars listed as Ofpe/WN9 as their own sub-sample. We set a cutoff in the RSG sample at log L/L_⊙ ≥ 4.6, close to the luminosity of RSGs with initial mass 9 M_⊙ (Smartt et al. 2009), which is consistent with the lower limits of type II-P progenitors estimated by Smartt et al. (2009) and Smartt (2015). Another cutoff was set at log L/L_⊙ ≥ 4.8, which is the completeness limit of the M 33 supergiant catalogue by Drout et al. (2012), for the purpose of comparison between the two galaxies. We also split the YSG samples at log L/L_⊙ = 4.8 for the same purpose.

The mean NCRs of the different stellar samples, with both accurate coordinates and simulated positional errors, are presented in Table 1 for the LMC and Table 2 for M 33. The errors reported for the mean NCRs are the SEM in the case of accurate coordinates and, as explained in Sect. 3, the analogous from the Monte Carlo analysis. The mean NCR values for different types of transient events from A12 and H14 are presented in Table 3.

Fig. 5 Top panel: cumulative NCR distributions of the main sequence stars of different spectral type in the LMC at a simulated distance of 35 Mpc, illustrating the trend of increasing NCR with earlier spectral type. Bottom panel: distributions of different types of evolved stars in the LMC at the distance of 35 Mpc, similarly illustrating the dependence of NCR on stellar type and hence initial mass. In both panels, shorter-lived stars tend toward the lower right and longer-lived stars toward the upper left.

Within the main-sequence sample in the LMC, it is evident that, in most cases, stars with an earlier spectral type – and thus higher initial mass – have a higher mean NCR, apart from the O8V and O9V stars, whose NCRs are consistent within the uncertainties. The mass-NCR correspondence is also visible when comparing evolved stars. Wolf-Rayet stars, with a lower initial mass limit around 20 or 25 M_⊙ (Crowther 2007; Hainich et al. 2014), have a substantially higher mean NCR than the supergiants which mostly occupy a lower initial mass range. For example, at 20 Mpc in the LMC, the WN sample has a mean NCR of 0.561 ± 0.031 while the RSG sample has 0.182 ± 0.010. Setting a higher luminosity cutoff (corresponding to higher minimum initial mass) to the supergiant samples also results in a higher mean NCR. For example, at 20 Mpc in the M 33, the RSGs with log L/L_⊙ ≥ 4.6 have a mean NCR of 0.261 ± 0.024 while those with log L/L_⊙ ≥ 4.8 have 0.330 ± 0.033. Figure 5 shows the cumulative distributions of the main-sequence and evolved stars (with accurate coordinates) in the LMC at a simulated distance of 35 Mpc to illustrate the mass-NCR sequence.

4.1. Differences between galaxies

Generally speaking, the NCR results in Tables 1 and 2 are quite similar between the two galaxies used in this study, despite the differences between the galaxies. M 33 is a late-type spiral galaxy with a metallicity gradient between roughly solar metallicity (Z_⊙) in the nucleus and 0.5 Z_⊙ in the outer disk (Magrini et al. 2007). Its baryonic mass is about 10¹⁰M_⊙ (Corbelli 2003). The LMC is a dwarf spiral galaxy with a baryonic mass about 3 ×10⁹M_⊙ (Kim et al. 1998) and exhibits signs of a bar structure. The LMC has a metallicity generally around 0.4 Z_⊙ without a clear gradient (Piatti & Geisler 2013).

The mean NCRs of the RSG samples with log L/L_⊙ ≥ 4.6 and log L/L_⊙ ≥ 4.8 are inconsistent within a 1σ limit but consistent within 2σ. They are also formally consistent according to the AD test: the probabilities for having identical distributions are 41 and 29 per cent respectively using the 20 Mpc distance. Because of completeness issues, the RSGs with log L/L_⊙ ≤ 4.6 in M 33 are biased toward lower values than in the LMC (0.095 ± 0.022 and 0.155 ± 0.011, respectively, at 20 Mpc). This is probably caused by a larger fraction of missing low-luminosity RSGs against bright regions in M 33. The YSGs with log L/L_⊙ ≥ 4.8 have both consistent mean NCRs and a 72 per cent AD test probability. The mean NCRs of WN-star samples (including WNh stars) are also self-consistent, with AD test probabilities of 23 and 26 per cent for early and late WN respectively. There is thus no apparent metallicity effect on the WN progenitors, in agreement with Hainich et al. (2014), who find a lower mass limit for WN stars in the LMC (~20 M_⊙) comparable to the limit in the Milky Way (~25 M_⊙). They suggest that this may be because the reduced mass loss caused by low metallicity implies less angular momentum loss and thus faster rotation, leading to lower minimum WR star mass.

The only significant difference between the galaxies concerns the early WC stars, which have a lower NCR in M 33 than in the LMC (0.551 ± 0.049 and 0.656 ± 0.045, respectively, with a 27 per cent AD test probability) – the late WC stars in M 33 more closely resemble the latter in terms of mean NCR, with a 34 per cent AD test probability. The simplest explanation could be that, at a lower metallicity, a higher initial mass would be required to make a WC star through the wind mass loss channel (Crowther 2007), which would result in higher NCRs for WC stars in the LMC at the same relative evolutionary stage as those in M 33. However, removing the central regions of M 33 with the highest metallicity (an ellipse with a ~ 6 kpc, b ~ 3.5 kpc and position angle 22.7 deg), and the 13 early WC stars inside this area, results in a mean NCR of 0.499 ± 0.073 at 20 Mpc for the remaining 15 early WC stars. Because of the large uncertainty, this is still consistent with the mean NCR of 0.551 ± 0.049 obtained using the whole galaxy, but does not seem to support the metallicity explanation for the ~2σ difference between the two galaxies. We also note that, while in the LMC the NCRs of WC stars are higher than those of WN stars, as expected from progenitor mass differences and as found by Leloudas et al. (2010) in M83 and NGC 1313, this is not the case in M 33. The results in M 33 thus remain peculiar.

4.2. Comparison with SN types

All of the stellar samples studied here, except the sample of RSGs with log L/L_⊙< 4.6 in M 33, show a mean NCR higher than that of the completely random distribution. This indicates that all these samples are correlated with the Hα emission at least weakly. The mean NCR of SN impostors from H14 (0.133 ± 0.086) is consistent with a random distribution, although the errors are quite large and the median distance of the impostor sample, 8.5 Mpc, somewhat different from our simulated distances.

In order to use our analysis to set constraints on CCSN progenitors, we compared the distributions of CCSN NCRs reported in A12, K13 and H14, and their mean NCRs listed in Table 3, to those of our evolved star samples at the appropriate distance. We obtained median distances for the SN samples using the redshifts (radial velocities) of the host galaxies reported in A12 and H14 and converting them into distances using H₀ = 70 km s^-1 Mpc^-1, Ω_m = 0.3, Ω_Λ = 0.7. The radial velocities of the K13 sample galaxies were obtained from the NED. As the median distances of type Ic SNe in K13 and the type II-P and II-L SNe in the A12 sample are closer to 20 Mpc, we compared them to our results using the 20 Mpc images. Similarly, the median distances of types Ib, Ic (A12), IIb and IIn are closer to 35 Mpc and therefore those types were compared to our results from the 35 Mpc images. The A12 S/N was used for all these comparisons, and the star NCR distributions calculated including positional errors. We used the AD test to compare the distributions of samples containing at least ten objects. The results from these comparisons are presented in Tables 4 and 5, with comparisons where the mean NCRs are also consistent highlighted. In previous NCR studies such as A12 and K13, the Kolmogorov-Smirnov test was used. The AD test, however, is more sensitive to differences between the distributions, especially near the endpoints of the cumulative distribution. Comparisons between consistent distributions of stars and SNe (with consistent mean NCRs and/or with an AD probability ≥ten per cent) are presented in Fig. 6, along with LBVs versus type IIn SNe to demonstrate their difference.

The best matches for each type can be found using Tables 1, 2, 4 and 5, defined here as both being formally consistent with the SN type according to the AD test (an AD probability ≥ten per cent) and having a consistent mean NCR. These are as follows for each SN type:

• Type II-P SNe are best matched by the RSG samples with logL/L_⊙ ≥ 4.6 or log L/L_⊙ ≥ 4.8 and SG B[e] stars. The YSGs in the LMC are consistent with type II-P according to the AD test but do not have a consistent mean NCR.

• Type II-L SNe are best matched by RSGs with log L/L_⊙ ≥ 4.8 in M 33, YSGs (with or without the cutoff) in both galaxies, as well as early WN stars without hydrogen and SG B[e] stars in the LMC. The samples of RSGs with log L/L_⊙ ≥ 4.6 or log L/L_⊙ ≥ 4.8, LBVs and early-type WN stars are consistent with type II-L according to the AD test, but inconsistent based on the mean NCR. The small sample size of type II-L affects these results.

• Type IIb SNe are best matched by RSGs with log L/L_⊙ ≥ 4.8 and YSGs (with or without the cutoff) in both galaxies, as well as early WN stars without hydrogen and SG B[e] stars in the LMC. LBVs and early WN stars in the LMC, Ofpe/WN9 stars, RSGs with log L/L_⊙ ≥ 4.6 and early WC stars in M 33 are consistent with type IIb according to the AD test but inconsistent based on the mean NCR. This is again at least partly due to the small sample size of type IIb.

• Type Ib SNe are best matched by the RSG samples with log L/L_⊙ ≥ 4.6 or log L/L_⊙ ≥ 4.8. In the LMC the AD test probabilities are lowered by having a different fraction of zero-NCR stars – otherwise, as seen in Fig. 6, the distributions are quite similar. SG B[e] stars also match type Ib SNe.

• Type Ic SNe (A12; random sample of spiral galaxies) are best matched by early WN stars (both in general and those without hydrogen) and LBVs. Several other distributions are consistent with type Ic according to the AD test but inconsistent based on the mean NCR, namely YSGs with log L/L_⊙ ≥ 4.8 in both galaxies, SG B[e] stars in the LMC as well as Ofpe/WN9 stars and early WC stars in M 33.

• Type Ic SNe (K13; strongly star-forming galaxies) are best matched by Ofpe/WN9 stars in M 33 and WR stars (both WN and WC) in both galaxies. The LBVs are also formally consistent with this SN sample based on the AD test but inconsistent based on the mean NCR.

• Type IIn SNe are only matched by the full sample of RSGs in M 33. The RSGs with log L/L_⊙ ≤ 4.6 are consistent with type IIn according to the AD test but inconsistent based on the mean NCR. Even these distributions seem quite different from type IIn by eye in Fig. 6. Furthermore, the LBVs show a distribution (and a mean NCR) very different from type IIn.

Fig. 6 Comparisons between the cumulative NCR distributions of various SN types and stellar samples in the LMC and/or M 33. Panel a) type II-P SNe; b) type II-L; c) type IIb; d) type Ib; e) type Ic from A12; f) type Ic from K13 and g) type IIn. The stellar samples have been selected based on formal consistency according to the AD test (≥ten per cent probability of the distributions being the same) and similar mean NCRs. In the case of the type IIn SNe, we also include the LBVs to demonstrate the difference of these distributions.

5. Systematic effects

A general effect immediately seen in Tables 1 and 2 is that including positional errors in the analysis results in slightly lower mean NCRs. For low initial NCRs the effect is small, and the value may even increase slightly. This is unsurprising because a small error in the position of a star in a high-NCR pixel mostly lowers its NCR, meaning that the neighbouring pixels are likely to be fainter, while a star in a low-NCR pixel can move to a fainter or brighter pixel. With σ = 0.5 arcsec and using the A12 pixel scale, ~70 per cent of the stars stay in the same pixel after the error is applied, and as a result the effect of this error is not very large. There is a small increase in the effect of the positional error with distance; this is simply because the same positional error in angular coordinates corresponds to a larger error in physical coordinates as the distance of the galaxy is increased.

In addition to the steps described in Sect. 3.2 to make our images comparable to those of A12, we also investigated the effects of the removal of the central region in M 33, as well as changing the distance, binning and S/N. We performed these tests on multiple subsamples with different mean NCR to test for the biases as a function of NCR. We additionally tested for any spatial biases in the RSG sample in the LMC caused by the patchy coverage of the catalogue used.

In M 33, we tested the effect of removing the central square kiloparsec region from the Hα image and excluding the stars located in this region from the stellar samples. As mentioned above, the Drout et al. (2012) RSG sample is biased against this region because of crowding effects; therefore removing the region and the two RSGs located in it should make the NCR result more accurate. We also tested what effect this has on some other NCRs. Briefly put, there is no significant effect on the NCR results; the re-calculated values are listed in Table 6. In particular, the RSGs, which are the primary subjects of this test, show slightly increased mean NCRs, but are still well inside the errors of the original values. The YSGs, which might also have been affected by the crowding bias, show practically no change either. This confirms our earlier assertion that the spatial coverage of the M 33 RSG sample is in practice unaffected by the crowding bias reported by Drout et al. (2012), assuming that the RSG population in the central region does not significantly differ from the population in the outer regions.

We also calculated some test NCRs at the distance of 75 Mpc, which is at the upper end of the galaxy distances in the A12 and K13 samples, to test how distance affects the NCR determination. This was again done using the appropriate convolution – 1 arcsec seeing at this distance corresponding to a resolution of about 360 pc, equivalent to 31.5 pix for the LMC and 44.1 pix for M 33 – and by re-binning by 31 for the LMC and by 44 for M 33. To illustrate the effect on the images of the LMC and M 33, the convolved and re-binned Hα images at different simulated distances are presented in Fig. 7. Furthermore, we calculated the test NCRs for different noise levels. For each distance (20, 35 and 75 Mpc), in addition to the Poisson noise required to achieve the S/Ns per pixel mentioned above, we made images with five times more and five times less noise. Tables 7 and 8 show the effects of varying distance and noise on some NCRs in the LMC and M 33, respectively, selected to cover a wide range of mean NCRs.

Fig. 7 Hα images of the LMC and M 33 at different simulated distances, using a ~1 arcsec pix-1 binning with a simulated distance and S/Ns similar to those in A12 and K13. North is up and east is left. The black scale bar corresponds to 2 kpc in the LMC and 5 kpc in M 33; the linear scale is the same in all images.

As can be seen from Tables 7 and 8, increasing the distance generally shifts NCRs closer to ~0.6 in the LMC and ~0.5 in M 33, that is, mean values under ~0.6 or ~0.5 are increased with increasing distance and values over these are decreased. The effect is significant for several samples, and demonstrates the need to simulate the appropriate distance before comparing our results to those of SNe. As the distance increases, the spatial resolution worsens, and sharp features are smeared over a larger apparent area. After the convolution, an initially low-NCR star may end up in a pixel with higher flux, while an initially very high NCR does not change much, as the brightest pixels remain so even with the convolution. Adding noise (and thus lowering the S/N) generally has a very different effect: initially low NCRs are lowered further. High values are also lowered, but they change much less because as the noise level increases, more pixels become indistinguishable from the background and their NCR becomes zero. This affects faint pixels first, resulting in more zero NCRs in initially low-NCR samples. Figure 8 illustrates this effect on the cumulative distributions of the NCRs of some stellar types. The NCRs of pixels that are relatively bright but not the very brightest also decrease: less and less flux originates in the pixels fainter than them, while the bright pixels themselves are comparatively unaffected by the noise. The difference between “high” and “medium” (A12) S/N is small above NCR ~ 0.15, indicating that not much Hα emission is lost at the A12 noise levels and that the S/N in the A12 study is adequate. With further degradation of the S/N, however, the effect on the NCR becomes significant up to NCR ~ 0.4.

Furthermore, we calculated some test NCRs with different binning. The NCRs of K13 were calculated using images with 0.38 arcsec pix^-1 and, in the case of type Ic SNe, are somewhat different from those in A12. To investigate what bias, if any, the different binning may introduce, we re-binned the 35 Mpc image (close to the median distance of both samples) to have this pixel scale. This meant re-binning the original LMC image by six and that of M 33 by eight. We also compared these to images with the pixel scale of the unbinned ALFOSC Hα images of K13, 0.19 arcsec pix^-1 (binning by three for the LMC and by four for M 33). The effect of increasingly coarse binning without convolution, as shown in Tables 9 and 10, is qualitatively the same as that of increased distance, which includes both convolution and re-binning. More and more pixels are averaged into one, and an initially low-NCR star can end up in a higher-NCR pixel. Concomitantly, high NCRs can be diluted by the inclusion of faint pixels in the same bin, but the highest values generally remain so. This effect is, however, not significant for most samples.

Fig. 8 Illustration of the effect of noise in the Hα image on the NCR distribution. The distributions with the dashed lines were obtained using an S/N approximately five times better than that of the A12 images; the solid lines represent the A12 S/N; the dotted lines represent a five times worse S/N. The main effect is an increase in the number of zero-NCR stars, and this effect is most significant on initially low NCRs.

Finally, we tested whether the spatial coverage of the Neugent et al. (2012), with 64 fields spread seemingly randomly over the visible extent of the LMC and with gaps in-between, introduces a spatial bias in the NCR analysis. Conceivably, the regions between the observed fields might host preferentially low or high numbers of RSGs. Using the 2MASS catalogue, the Naval Observatory Merged Astrometric Dataset⁵ (NOMAD) and the Wide-field Infrared Survey Explorer catalogue⁶ (WISE), we searched for RSG candidate stars that fulfil the following criteria: the proper motions of the stars, obtained from NOMAD, are compatible with the Magellanic Clouds (González-Fernández et al. 2015); the magnitudes of the stars are K_2MASS ≤ 10.2 mag, adhering to the magnitude limit for RSGs in the Neugent et al. (2012) catalogue and corresponding to an absolute magnitude of M_K ≤ −8.3 mag; and the stars are visible in the WISE 3.4 and 4.6 μm bands. Roughly three times more 2MASS stars in the LMC (1670) match these criteria than are included in the Neugent et al. (2012) RSGs (543), which agrees with the completeness estimate by the catalogue authors. A visual comparison between the spatial distributions of these candidates and the RSG sample shows no obvious bias toward any particular surface brightness. The mean NCR for the candidates, using the image of the LMC at 20 Mpc with no positional errors, is 0.188 ± 0.006, which is consistent with the mean NCR of the RSG sample at that distance (0.182 ± 0.010). Furthermore, an AD test between the two distributions reveals that they are formally consistent (18 per cent probability of a common population). Thus we conclude that no significant biases are introduced by the coverage of the Neugent et al. (2012) catalogue.

6. Discussion

Looking at the results of main-sequence stars in Tables 1 and 2, with or without including positional errors, one can easily see a trend of higher NCR with earlier spectral type and thus higher initial mass. The trend is illustrated in Fig. 5. This sequence of mean NCRs provides evidence that NCR is a useful indicator of progenitor lifetime and mass. Unfortunately, a direct comparison between the NCRs of main-sequence stars and SNe is not feasible. SNe are by definition at the terminal stage, while the main-sequence stars in a galaxy are, on average, at the midpoint of their life assuming a constant star formation rate in the galaxy over the past few tens of Myr. Therefore an average OB main-sequence star still has time to move during its lifetime and change its NCR. Furthermore, the coordinates and spectral types of these stars are from SIMBAD and thus from multiple independent studies that may have concentrated on specific locations such as 30 Doradus. Therefore there may be severe completeness biases in the samples of main-sequence stars reported in SIMBAD. A sign of a bias is visible by eye in Fig. 3, where many main-sequence B stars are concentrated in a (marked) rectangular area with no strong Hα emission sources. The NCRs of main-sequence stars should only be regarded as qualitative evidence of a lifetime-NCR relation.

More fruitful and quantitative results can be obtained using the spatially unbiased, or even complete, samples of evolved stars. This is not only because of the biases mentioned, but because the time a star spends as a supergiant, LBV or WR star after the main sequence is short. The longest of these phases, the RSG phase, lasts up to about 1 Myr. Thus at a (projected) velocity of ~10 km s^-1, normal for an RSG in the LMC (Olsen & Massey 2007), a star can move about 10 pc – roughly half a pixel in our simulated 20 Mpc images – during this phase. The YSG and LBV phases last on the order of 10⁴ or 10⁵ yr in single star evolution models (Neugent et al. 2012; Smith 2014), while the WR phase lasts a few times 10⁵ yr (Crowther 2007). Therefore, unlike during the main-sequence phase, a star does not move much during these phases and the H ii region similarly does not change much, and hence its NCR should be roughly the same as that of the resulting SN. If the YSG phase does not end in a SN (doing so requires rotation and enhanced mass loss according to Georgy 2012), it will be followed by another evolutionary stage, but this stage should also be short.

There is immediate evidence for the validity of the NCR method in the results using evolved stars as well. The initial masses of both LBVs and WN stars are in the ≥ 25 M_⊙ range (Smith et al. 2004; Crowther 2007, respectively, including the low-luminosity LBVs), and the mean NCRs of LBVs and WN stars in general are consistent within the errors. WC stars have a higher mean NCR because of their higher minimum initial mass of ~ 40 M_⊙ (Crowther 2007). Even though WC stars are at a relatively later stage of evolution than WN stars, the higher mass range means they are still younger on average. The NCRs of these massive evolved stars are significantly higher than those of the mostly lower-mass YSGs and RSGs. Inside the RSG and YSG samples, a higher luminosity cutoff (and hence higher minimum mass) results in a higher NCR as well. Furthermore, SG B[e] stars have a lower NCR than LBVs, which is compatible with the Smith & Tombleson (2015) argument that they are lower-mass analogs of LBVs. The pixel statistics method was criticized by Crowther (2013) and Smith & Tombleson (2015), who argued that, as the brightest giant H ii complexes are the longest-lived ones, an association with them should not be a measure of lifetime, and that the most massive stars should be located in small individual H ii regions. However, as discussed in Sect. 3, we would still expect a lower likelihood for a longer-lived star to remain in its native giant H ii region (into which the smaller regions would blend at the spatial scales probed by A12), and our NCR results indicate that this is true.

In addition to connecting SNe to their progenitor stars, the mean NCRs of the subsamples of evolved and main-sequence stars can, in principle, constrain stellar evolutionary channels. If a sample of evolved stars has a higher mean NCR than one of main-sequence stars, stars in the latter sample cannot predominantly be the precursors of the former – even though some exceptions to this rule would be possible. However, there are problems with this approach. As the NCR of a main-sequence star can change significantly before the following stages of its evolution, a direct comparison between the NCR distributions is no more fruitful than between the distributions of main-sequence stars and SNe. Furthermore, an unbiased catalogue of main-sequence stars would be required for meaningful constraints. One can, nonetheless, use this approach as another consistency check on the NCR method. Some mean NCRs of main-sequence stars in the LMC do seem incompatible with some of the evolved stars: types B1V and B2V (1-σ upper limit of 0.554 at 20 Mpc without positional errors) vs. WC stars (≥ 0.611) or late WN stars (≥ 0.618); types B0V and later (≤ 0.649) vs. classical LBVs (≥ 0.659). With other evolved types, such constraints cannot be made, as their mean NCRs are consistent with or lower than those of any of our main-sequence subsamples. As the precursors of WR stars and LVBs are believed to be much more massive than B-type main-sequence stars, these results are consistent with our present view of stellar evolution.

6.1. SN progenitor masses

As the lifetime-NCR connection now seems secure, we can try to use it to constrain the SN progenitors. Our first star-SN comparison is between type II-P SNe and RSGs, as the latter are already firmly established as the progenitors of the former. A luminosity of log L/L_⊙ = 4.6 corresponds to an initial mass of ~9 M_⊙ for an RSG with a weak dependence on metallicity according to Smartt et al. (2009), who use the Cambridge University stars code (Eldridge & Tout 2004). This is close to the minimum initial mass required for a type II-P progenitor (e.g. Smartt 2009, 2015). RSGs with this luminosity cutoff provide a good match to type II-P SNe in both galaxies (Sect. 4.2). The comparison is affected by lifetime effects. While the sample of type II-P SNe, with detected SNe covering the entire sky, probes the entire initial mass function (IMF) of type II-P progenitors, the probability of a star to be observed as an RSG depends on the corresponding RSG lifetime. Therefore, the RSG samples include fewer high-mass RSGs with high NCRs than expected from the IMF, and the RSGs with log L/L_⊙ ≥ 4.8 are also consistent with type II-P SNe in the LMC. In M 33, the more top-heavy RSG distribution in the catalogue (because of completeness issues) serves to offset this. As such, we conclude that the type II-P and RSG (log L/L_⊙ ≥ 4.6) distributions are consistent. Thus the NCR method succeeds in reproducing a progenitor range for type II-P SNe similar to that predicted by stellar evolution models and confirmed by direct progenitor detections.

The A12 IIb sample size is small and it is difficult to narrow down the possible progenitors. However, the similarity between our YSG samples and the type IIb SNe is consistent with the detected supergiant progenitors of type IIb SNe and with the 12−16 M_⊙ initial mass range suggested by Jerkstrand et al. (2015). The initial mass range of YSGs in the catalogues (Drout et al. 2012; Neugent et al. 2012) is between about 10 and 25 M_⊙ according to Geneva stellar evolution models, depending on rotation. The average luminosity of the YSGs in the LMC catalogue is log L/L_⊙ = 4.7, corresponding to about 15 M_⊙, while in M 33 the higher completeness limit results in an average of log L/L_⊙ = 5.0, corresponding to between 15 and 20 M_⊙. Type II-L SNe are consistent with similar stellar samples and have a mean NCR consistent with type IIb SNe, indicating a similarity in progenitor mass as well. Other results (e.g. Elias-Rosa et al. 2011; Faran 2014; Terreran et al. 2016) are also consistent with masses around 15 M_⊙ for type II-L progenitors. The small sample size of type II-L, again, makes it difficult to draw robust conclusions. However, the presence of a binary companion may result in stronger mass loss in type IIb progenitors despite the similar mass, leaving behind less hydrogen than would be expected from a single star.

As for type Ib and Ic SNe, A12 and K13 found their NCR distributions to differ significantly. The indication of this is that while type Ic progenitors may be relatively massive (whether single or binary), type Ib progenitors are not, on average, much more massive than the RSG progenitors of type II-P and therefore are likely to be in interacting binaries. As Fig. 6 shows, and as is apparent from Tables 1, 2, 4 and 5, type Ib SNe are consistent with RSGs with a log L/L_⊙ ≥ 4.6 or log L/L_⊙ ≥ 4.8 cutoff. While RSGs themselves cannot be the immediate progenitors of stripped-envelope SNe such as type Ib, this does provide an indication that the initial masses of type Ib progenitors are in a range similar to that of type II-P progenitors, located in interacting binaries. This is consistent with the suggested binary progenitor of iPTF13bvn, which had an initial mass of 10−12 M_⊙ (Eldridge & Maund 2016), similar to those of the RSG progenitors of several type II-P SNe (Smartt 2015). This result is also corroborated by, for example, Kuncarayakti et al. (2013a) who find type Ib progenitors to be older than those of type Ic and likely to be located in interacting binaries. The slightly higher mean NCR compared to type II-P is likely a result of the higher median distance of the Ib sample (Sect. 5).

Type Ic SNe, on the other hand, have an NCR distribution, and a mean NCR, similar to the early-type WN stars. The initial mass range of single WN stars is ≳ 25 M_⊙ (Crowther 2007) at Z_⊙, and Hainich et al. (2014) find a limit of ≳ 20 M_⊙ in the LMC; the similarity of this stellar type with type Ic SNe lends support to type Ic progenitors sharing this initial mass limit. Whether type Ic SNe are actually the end-points of single WN evolution (possibly through a subsequent short WC phase) or originate in more massive interacting binaries than type Ib cannot be said based on our results. A mix of WC stars and lower-mass interacting binaries is also possible. The fact that no type Ic SN progenitor has yet been detected is sometimes used to exclude WR stars as a significant type Ic SN progenitor channel (e.g. Smith 2014). However, Yoon et al. (2012) and Groh et al. (2013) argued that these stars would be optically faint shortly before the SN despite their high bolometric luminosity, and the lack of detected progenitors is thus not surprising. The K13 sample of type Ic SNe has an even higher mean NCR and a distribution more similar to WC or Ofpe/WN9 stars or a mixture of early- and late-type WN stars. The difference between the NCR distributions in A12 and K13 is likely a combination of factors: a different median distance and binning in the two studies contributes to the higher NCRs in K13, but a selection effect or a physical difference between the samples of host galaxies is probably needed as well. The K13 galaxies were selected for a high far-infrared luminosity and thus a high star formation rate. It is possible that a recent starburst episode in these galaxies results in type Ic progenitors being younger (and thus more massive, perhaps WC stars) on average – or still being preferentially located inside their native H ii region – thus artificially increasing their NCR. Differences between the NCR distributions of other CCSN types between A12 and K13 are not significant; longer-lived progenitors may not be similarly affected by a very recent starburst episode.

LBVs have been suggested by various studies as the immediate progenitors of at least some type IIn SNe (Trundle et al. 2008; Gal-Yam & Leonard 2009; Smith et al. 2011b; van Dyk & Matheson 2012; Fransson et al. 2014). The initial masses of these progenitors are estimated to be above 30 or even 50 M_⊙ (Gal-Yam & Leonard 2009; Smith et al. 2011b). It is, however, clear from Fig. 6 that the NCR distributions of type IIn SNe and LBVs do not match. We also note that some lower-luminosity LBVs are actually post-RSG stars according to the models of Meynet & Maeder (2005) and Ekström et al. (2012) and should be treated separately. The three classical LBVs in the LMC are consistent with stars of types O8V through O4V, a result similar to that of Humphreys et al. (2016) (although the numbers of these objects in both studies are very small). Lower-luminosity LBVs, on the other hand, are consistent with YSGs with log L/L_⊙ ≥ 4.8. Both of these types have a mean NCR significantly higher than that of type IIn SNe; therefore only a small fraction of IIn progenitors can be LBVs. Type IIn is a very diverse class of transients, likely to include SNe from multiple, very different channels such as super-AGB stars exploding as electron capture SNe (as suggested for 1994W-like SNe by e.g. Kankare et al. 2012), extreme RSGs (as suggested for 1998S-like SNe by e.g. Taddia et al. 2015) or type Ia events disguised as IIn (e.g. Fox et al. 2015). A combination of the latter three types, with a few LBVs added into the mix, could account for the observed NCR distribution. For example, RSGs with log L/L_⊙< 4.5, corresponding to M ≲ 8 M_⊙ (Smartt et al. 2009), at 35 Mpc in the LMC have a mean NCR of 0.143 ± 0.012, which is also close to the NCR of type Ia SNe (A12). A population of 70 per cent such RSGs, which probably include a significant fraction of AGB stars, and 30 per cent LBVs drawn from both LBV subclasses, results in a mean NCR of 0.258 ± 0.026. This is still consistent with the type IIn SNe within the errors.

SN impostor events have been connected to LBVs (e.g. van Dyk et al. 2002; Maund et al. 2006; Tartaglia et al. 2015); however, there is evidence (Smith & Tombleson 2015, and references therein) that lower-mass stars could also be responsible for similar events (e.g. a 12−17 M_⊙ star suggested to be the progenitor of NGC300 OT2008-1 by Gogarten et al. 2009). These events were investigated along with type IIn SNe by H14, but their mean NCR is only 0.133 ± 0.086 – although the median distance of their sample is only 8.5 Mpc and, as we have seen, increasing the distance will also increase low NCRs. One also needs to keep in mind that as SN impostors are significantly fainter than bona fide SNe, their detection may be heavily biased against bright backgrounds, and such regions often have high NCRs (this also causes the low median distance of the sample). As such, the NCR of the impostors may be heavily biased toward low values.

6.2. Effects of binary evolution

The initial mass ranges of WR stars reported by Crowther (2007) and Hainich et al. (2014), the evolutionary tracks presented in Neugent & Massey (2011) and Hainich et al. (2014) and the initial masses of RSG progenitors of CCSNe reported by Smartt et al. (2009) and others are the results of single star models. However, binary interaction dominates the evolution of massive stars: over 70 per cent of all O-type stars are expected to exchange mass with a companion, with a binary merger in a third of the cases (Sana et al. 2012). Only about half of apparently single massive stars, those with an amplitude of radial velocity variation less than 10 km s^-1, are truly single (de Mink et al. 2014), with the other half being wide binaries before mass transfer or products of mass transfer or a merger. Therefore it is worthwhile to discuss possible effects that the multiplicity of a star may have on the NCR-mass connection.

The time a star spends as an RSG in a binary system is shortened by the mass transfer by a factor of two or three (Eldridge et al. 2008). Thus the binary RSGs that are included in the sample make the average duration of the RSG phase shorter. Therefore the NCRs of RSGs in binary systems should be closer to those of the resulting SNe than expected from single stars. However, as the phase is already short compared to the total lifetime of the star, the change in the mean NCR should not exceed our 1σ uncertainty.

The apparent binary fraction in the Hainich et al. (2014) catalogue is about 15 per cent, with an additional 20 per cent labelled binary suspects. However, based on the statistics reported by de Mink et al. (2014), the total fraction of WN stars affected by binary evolution may be as high as about 60 per cent. Eldridge et al. (2008) found that while the change in the lifetime of a WR star of a particular mass is minor between single- and binary-star models, the minimum initial mass for a WR star in a multiple system decreases to around 15 M_⊙. In addition, the mean NCR may be changed by effects of mergers and mass transfer, which in some cases lengthen the lifespan of the star beyond what is expected from a single star of the same initial mass (de Mink et al. 2014). Because of the uncertainties involved in the mass-lifetime connection in binary systems, it is difficult to say what the effect on the NCRs is. This convolution is, however, already included in the NCRs of a particular SN or stellar sample. The average initial masses of the WR stars may in any case be lowered by the effect of binaries, which would also impact the inferred progenitor mass range of type Ic SNe.

7. Conclusions

We have used the pixel statistics method to calculate the spatial correlations of massive stars and Hα emission in two nearby galaxies, the LMC and M 33. We have simulated the typical distances, noise levels and binning of previous pixel statistics studies of more distant SN host galaxies and, for the first time, directly compared the NCR distributions of different types of SNe to different stellar samples. We find the following:

• The pixel statistics method shows a correlation between a shorterlifetime, and hence higher initial mass, and higher NCR. Thisapplies to both the main sequence OB stars in the LMC andevolved stars in both galaxies. In addition, the NCRs of the stellarsamples are generally consistent between the two galaxiesdespite different metallicities, structures and star formationhistories. Thus we find NCR to be a valuable tool in studying theprogenitor stars of CCSNe in the local Universe indirectly.

• The results from the pixel statistics method do not change significantly between an S/N typical to A12 and one five times higher. Whether the pixel scale is typical to A12 or five times smaller likewise does not have a significant effect. Therefore the aforementioned studies do not contain obvious biases resulting from these factors.

• The NCR distribution of RSGs with initial masses ≳9 M_⊙ matches that of type II-P SNe, consistently with the II-P progenitor masses derived from direct progenitor detections. Similarly, the NCRs of YSGs with an average initial mass of ~15 M_⊙ are consistent with those of type IIb SNe. Type II-L SNe are consistent with sharing a progenitor mass range with type IIb, suggesting another factor such as the presence or absence of a binary companion as the main difference between these types.

• The NCR distribution of type Ib SNe is consistent with stars with initial masses ≳9 M_⊙, indicating that the dominant progenitors for these SNe are not single WR stars but lower-mass interacting binaries. The NCR distribution of type Ic SNe, on the other hand, is consistent with early-type WN stars, typically with initial masses ≳20 M_⊙ according to single-star evolution models, with a contribution from binaries down to 15 M_⊙.

• Both classical (initially ≳50 M_⊙) and low-luminosity (initially ≳25 M_⊙) LBV stars have an NCR distribution inconsistent with type IIn SNe, with a significantly higher mean NCR. Therefore, the progenitors of type IIn SNe are probably mostly stars with a lower initial mass.

Acknowledgments

We thank the anonymous referee for many good suggestions that helped improve the paper, and Philip Massey and Brian Skiff for advice on the stellar catalogues in the LMC. This research has made use of the Southern H-Alpha Sky Survey Atlas (SHASSA), which is supported by the National Science Foundation; the SIMBAD database, operated at CDS, Strasbourg, France; the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration; and data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This work was partly supported by the European Union FP7 programme through ERC grant number 320360. R.G.I. thanks the STFC for his Rutherford Grant number ST/L003910/1.

References

All Tables

All Figures

	Fig. 1 Original continuum-subtracted Hα intensity maps used in this study: the LMC on the left (from the SHASSA survey; Gaustad et al. 2001) and M 33 on the right (Hoopes & Walterbos 2000), cropped in the case of the LMC but otherwise unaltered. North is up and east is left. The black scale bar corresponds to 2 deg (~1.7 kpc) in the LMC and 20 arcmin (~4.9 kpc) in M 33.
In the text

	Fig. 2 Hα images of the LMC (left) and M 33 (right), with our RSG, YSG, WN and WC samples overplotted in order to illustrate the spatial coverage of their respective catalogues. The black scale bar corresponds to 2 deg (~1.7 kpc) in the LMC and 20 arcmin (~4.9 kpc) in M 33.
In the text

	Fig. 3 Hα images of the LMC, with our main-sequence B and O star samples as well as the SG B[e] stars and LBVs overplotted. A concentration of B stars is visible in a rectangular area (marked in red) with no strong Hα sources, indicating a spatial bias in this subsample (see also Sect. 6). The black scale bar corresponds to 2 deg (~1.7 kpc).
In the text

	Fig. 4 Hα images of the LMC (left) and M 33 (right), with the coordinates of our randomly generated uniformly distributed samples of positions overplotted.
In the text

Fig. 5 Top panel: cumulative NCR distributions of the main sequence stars of different spectral type in the LMC at a simulated distance of 35 Mpc, illustrating the trend of increasing NCR with earlier spectral type. Bottom panel: distributions of different types of evolved stars in the LMC at the distance of 35 Mpc, similarly illustrating the dependence of NCR on stellar type and hence initial mass. In both panels, shorter-lived stars tend toward the lower right and longer-lived stars toward the upper left.

In the text

Fig. 6 Comparisons between the cumulative NCR distributions of various SN types and stellar samples in the LMC and/or M 33. Panel a) type II-P SNe; b) type II-L; c) type IIb; d) type Ib; e) type Ic from A12; f) type Ic from K13 and g) type IIn. The stellar samples have been selected based on formal consistency according to the AD test (≥ten per cent probability of the distributions being the same) and similar mean NCRs. In the case of the type IIn SNe, we also include the LBVs to demonstrate the difference of these distributions.

In the text

	Fig. 7 Hα images of the LMC and M 33 at different simulated distances, using a ~1 arcsec pix-1 binning with a simulated distance and S/Ns similar to those in A12 and K13. North is up and east is left. The black scale bar corresponds to 2 kpc in the LMC and 5 kpc in M 33; the linear scale is the same in all images.
In the text

Fig. 8 Illustration of the effect of noise in the Hα image on the NCR distribution. The distributions with the dashed lines were obtained using an S/N approximately five times better than that of the A12 images; the solid lines represent the A12 S/N; the dotted lines represent a five times worse S/N. The main effect is an increase in the number of zero-NCR stars, and this effect is most significant on initially low NCRs.

In the text