Free Access
Issue
A&A
Volume 645, January 2021
Article Number A6
Number of page(s) 17
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202037744
Published online 21 December 2020

© ESO 2020

1. Introduction

Core-collapse supernovae (SNe) occur at the end of the evolution of massive stars (e.g., Baade & Zwicky 1934; Bethe et al. 1979; Woosley et al. 2002; Heger et al. 2003). They are essential for the chemical enrichment of galaxies and as sources of feedback (e.g., Larson 1974; Hopkins et al. 2014). Supernovae are classified observationally into two main groups (e.g., Filippenko 1997; Gal-Yam 2017): hydrogen-rich, Type II SNe (including types II-P, II-L and IIn) and stripped-envelope SNe, which either show the complete absence of hydrogen features in their spectra or which show evidence for hydrogen only at early times (including types IIb, Ib, Ic, and more exotic examples such as Type Ic-broadline and Ibn). The first group are thought to arise from massive stars which have retained most of their hydrogen-rich envelope until the explosion. This study focuses on the progenitors of these hydrogen-rich, Type II SNe (hereafter referred to as SN II).

The reason why some massive stars end their lives with most of their hydrogen-rich envelope still intact and others do not is not fully understood. The canonical picture considers the progenitor stars as single stars that evolve in isolation. They can lose their envelope through stellar-wind mass loss or eruptive mass loss episodes (Conti 1975). In this picture, SN II are proposed to originate from single stars with initial masses between ∼8 and ∼25 − 30 M at solar metallicity (e.g., Heger et al. 2003; Eldridge & Tout 2004; Georgy et al. 2012; Renzo et al. 2017).

Massive stars with initial masses in this range are very often found in close binary systems (Kobulnicky & Fryer 2007; Eggleton & Tokovinin 2008; Chini et al. 2012; Kiminki & Kobulnicky 2012; Sana et al. 2012; Dunstall et al. 2015; Moe & Di Stefano 2017; Almeida et al. 2017). A large fraction of them are expected to interact with their stellar companion during their lifetime. Such interactions can remove the envelope of the donor star and give rise to a stripped-envelope SNe. This has been the topic of various theoretical studies (e.g., Podsiadlowski et al. 1992; De Donder & Vanbeveren 1998; Yoon et al. 2010, 2017; Claeys et al. 2011; Eldridge et al. 2013; Zapartas et al. 2017a; Sravan et al. 2018; Laplace et al. 2020, and references therein) in order to help explain their high relative rates (e.g., Smartt et al. 2009; Smith et al. 2011; Li et al. 2011; Eldridge et al. 2013; Graur et al. 2017), their low ejecta masses (e.g., Drout et al. 2011; Taddia et al. 2015; Lyman et al. 2016; Prentice et al. 2019), and the difficulty of finding their progenitors (e.g., Eldridge et al. 2013; Cao et al. 2013; Van Dyk et al. 2018).

The possible role of binary interaction in the lives of hydrogen-rich SN II is less explored, but this topic is gaining interest (e.g., Podsiadlowski et al. 1992; Podsiadlowski 1992; Tutukov et al. 1992; Vanbeveren et al. 2013; Justham et al. 2014; Menon & Heger 2017; Menon et al. 2019; Eldridge et al. 2018, 2019). In a previous study (Zapartas et al. 2019, hereafter referred to as Z19), we used population synthesis simulations to investigate the rates of different binary evolutionary channels that can lead toward hydrogen-rich Type II SN progenitors. Our simulations presented in that work suggest that a significant fraction of all SN II progenitors, ∼30−50%, are expected to have interacted with their companion prior to explosion. This prediction is robust against several variations in the model assumptions. Almost all of these binary progenitors are either stars that gained mass through Roche-lobe overflow (RLOF) or stars that merged with their companion. The remainder of the SN II progenitors evolve without experiencing any significant type of interaction. They effectively evolve as single stars in isolation. These may be either true single stars or stars that have one or more companions orbiting too far away for mass transfer to occur.

If binary interaction indeed plays a significant role in the lives of the population of SN II progenitors, as our previous simulations suggest, then this would have various implications that deserve further investigation. For example, we would expect a diverse set of observable consequences, such as possible influence on their light curves and systematic differences in the offset between the location of the explosions and the nearest star-forming region (see Eldridge et al. 2018; Z19 for a discussion).

The aim of this work is to investigate, in a statistical manner, the expected impact of binary evolutionary channels on SN II progenitors. We focus on the distribution of final core masses. The final core mass is one of the primary parameters that characterize the final state of a massive star. It is thought to play a key role in determining its fate and is a primary parameter that affects the late pre-SN observables such as the luminosity of the progenitor (e.g., Giannone 1967; Sukhbold et al. 2016; Farrell et al. 2020). Other parameters, such as the details of the internal density, composition and angular momentum profile and magnetic field strength also play a role, but at present it is much harder, or even impossible, to make reliable statistical predictions for these quantities. Even the prediction of final core masses for stars that are the products of binary interaction is not free from caveats, but we consider this as a natural starting point.

We first present the initial – final core mass relation that holds for single stars and show how binarity breaks this one-to-one relation. We find that Type II SN progenitors in binary systems typically originate from lower initial masses than single stars that reach similar core masses, because they gained mass or merged with a companion before explosion. We find that binary star progenitors of SN II typically end their life with a larger core mass than they would have had if they had lived in isolation. We investigate the contribution of different single and binary evolutionary channels and show the implications for the distribution of final core masses expected for Type II SNe.

We then place our theoretical results in context of the observations, using the sample of SN II compiled by Smartt et al. (2009), Smartt (2015). This sample consists of 26 nearby SN II, for which high resolution, archival images are available of the explosion site. These images have been used to derive constraints on the properties of the progenitors (see also Van Dyk et al. 2003; Maund & Smartt 2005). We compare our predictions with these constraints, discuss the validity and limitations of this comparison and the impact of our model assumptions. We further discuss how binarity affects the so-called red supergiant problem, that is, the claimed lack of red supergiant progenitors with initial masses above ∼17 M (Smartt et al. 2009, but see also Davies & Beasor 2018, 2020). We also discuss the implications of binary interaction in other techniques that infer the initial mass of Type II progenitors and predict a systematic discrepancy among them.

The paper is structured as followed. We explain our method in Sect. 2, present our results in Sect. 3 and place them in context of SN progenitor detections (Sect. 4), after exploring the sensitivity of our findings to model assumptions. We discuss the implications of our results for other observational constraints of Type II SN progenitors in Sect. 5 and summarize our conclusions in Sect. 6. Our paper can be considered as a companion study to Z19, but can be read independently.

2. Method

In order to study the end fate of Type II SN progenitors we simulate the evolution of entire populations of single and binary stellar systems. For these simulations we employ the stellar population synthesis code binary_c, developed by Izzard et al. (2004, 2006, 2009) with updates described in de Mink et al. (2013) and Schneider et al. (2015). In this code, the evolution of the stars is based on analytical fitting formulae (Hurley et al. 2000) derived from the grid of detailed simulations of Pols et al. (1998). Hurley et al. (2002) presented the implementation of these evolutionary formulae in the context of population synthesis simulations, which consider also physical processes in binary systems.

This code allows us to explore the initial parameter space of our population, which is very extended due to the numerous possible configurations of binary systems. Using this code, we have created a grid of simulations to study the statistical properties of core-collapse SNe (Zapartas et al. 2017b, hereafter Z17). In this study we focus on models describing the evolution of Type II SN progenitors that retained their hydrogen-rich envelope until the explosion, as in the case of Z19. Our main grid consists of 104 single stars of various initial masses and 150 × 150 × 150 binary systems varying their initial mass, their mass ratio and their orbital period. We assign a probability of formation to each of these stellar systems according to our assumed initial distributions. For our standard assumptions we follow an initial mass function (IMF) for our single and primary stars from Kroupa (2001), and an initial mass ratio and period distribution following Sana et al. (2012) for initially O-type primaries. The latter becomes flatter for initially lower mass systems (Öpik 1924). The slope of the initial distributions above are referred to as αIMF, κ, and π, respectively. Our results are described in the context of populations that consist exclusively of single stars or alternatively of binary systems, but we present also our findings for more realistic populations that include a combination of these extreme cases, assuming a binary fraction of fbin = 0.3, 0.5, 0.7 and a mass dependent case fbin(M), as in Z17.

The speed of population synthesis codes allows us to investigate the sensitivity of our main findings to variations in our input assumptions, testing the robustness of our results. These assumptions concern the initial stellar population, as specified above, as well as the treatment of the physical processes that take place during the evolution of single and binary stars. We first present the results of our fiducial model, in which we choose a standard set of assumptions that we describe here. For a more detailed description of all the input model assumptions we refer to Z17. We also perform simulations in which we vary our main model assumptions one-by-one, as in Z17, Z19. An overview is given in Table A.1. For consistency, we link and refer to the model variations that are also presented in Z17 with the same model number.

To account for the effect of mass loss on the evolution of stars, we follow the standard rate prescriptions from Vink et al. (2000, 2001) and Nieuwenhuijzen & de Jager (1990) during different evolution phases of hydrogen-rich stars, as well as the rate of Hamann et al. (1995) reduced by a factor of 10 for stripped helium cores (Yoon & Langer 2005). We also consider the effect of varying wind mass loss efficiency, η, between the range of 0.1 to 3 (Models 25, 26 and 51−54). The chosen metallicity for our main results is solar-like (Asplund et al. 2009, Z = 0.014) although we also consider different, mainly lower metallicities (Models 40, 42, 44).

We treat stellar spin following Hurley et al. (2000) and account for the effects of tides in binary systems (Zahn 1977; Hut 1980, 1981; Hurley et al. 2002). We simulate stable Roche-lobe overflow removing mass from the donor so as to remain inside its Roche lobe but limiting its mass loss rate according to its thermal rate. The accretion of matter from the companion is limited to 10 times its thermal rate (Schneider et al. 2015) although we run also simulations in which we vary the efficiency of this process as a fixed parameter (β, Models 1 to 3). The mass lost from the system is assumed to carry the specific orbital angular momentum of the accretor, with variations also ranging from no angular momentum loss to losses as high as to form a Keplerian circumbinary disk with a radius described in Z17 (parameter γ, Models 4 and 5). We account for rejuvenation of the accreting star or merger according to Tout et al. (1997), de Mink et al. (2013) and Schneider et al. (2015), considering variations of the fraction of mass lost and mixed in case of main-sequence mergers (parameters μloss and μmix respectively, Models 6 to 9).

We consider the possibility of unstable mass transfer in cases of mass ratios more extreme than some critical value (Maccretor/Mdonor <  qcrit), depending on the evolutionary state of the donor. This accounts for the possibility of a runaway mass transfer process as well as the case of simultaneously filling of the Roche lobes by both stars. We follow Hurley et al. (2002) for the mass-ratio threshold of qcrit, apart from the case of a donor crossing the Hertzsprung gap where we follow de Mink et al. (2013). In the case of two main sequence stars, such an instability leads to contact and merging, while in the other cases it leads to common envelope evolution. We consider different values of the mass ratio stability criteria (see Table A.1). We focus mostly on variations during the early evolutionary phases of the donor and have not explored the uncertainty of stability for Red Supergiant donors in initially wide orbits (e.g., Ge et al. 2010; Shao & Li 2014; Pavlovskii & Ivanova 2015; Pavlovskii et al. 2017) as the number of type II progenitors originating from the latter is low in our fiducial simulation (Fig. 4 in Z19).

The outcome of the CE is the spiral-in of the components resulting either in a merger or the ejection of the common envelope, according to the αCE-prescription (Webbink 1984). In our fiducial model we assume αCE = 1 and a gravitational binding energy parameter λCE according to Dewi & Tauris (2000, 2001) and Tauris & Dewi (2001). We consider other values for both in our variations (Table A.1).

We account for supernova kicks, both due to instantaneous mass loss (Blaauw 1961) as well as a possilbe natal kick to the remnant due to asymmetries in the explosion. These kicks may disrupt a binary system in most of the cases (Eldridge et al. 2011; Renzo et al. 2019). In our fiducial model we select random natal kicks from a Maxwellian distribution of root-mean square σ = 265 km s−1, based on the work of (Hobbs et al. 2005) on observed radio pulars. Natal SN kicks are uncertain, with evidence for lower kicks in some cases (e.g., Verbunt & Cator 2017). To test the sensitivity of our results to this assumption we also consider extreme variations of no and or very high natal kicks, affecting the rate of disrupted binary systems and thus of mass gainer type II SN progenitors (Z19). For a more in depth study of the impact of kick to star ejected from binary systems, we refer the reader to Renzo et al. (2019). If the system remains bound, we recalculate its orbital parameters. In our model variations we also consider the possibility of core-collapse events resulting in direct fallback onto a black hole with no observable transient (e.g., O’Connor & Ott 2011; Sukhbold et al. 2016; Adams et al. 2017a) for cases in which the final core of the progenitor is more massive than the core mass corresponding to a single star of Mmax, ccSN (Models 29, 52, 53, 54). We also run simulations where we artificially change the core mass criterion for a collapse, essentially slightly varying the minimum mass threshold for a core-collapse SN, Mmin, ccSN (Models 30 and 31).

3. Results

In this section we discuss the impact of binary interaction on the distribution of the final core masses of Type II SNe progenitors. Throughout this section and those that follow, we use the term “final core mass” to refer to the final mass of the helium core or, more accurately, the hydrogen-depleted core, MHe, preSN, of the supernova progenitor immediately prior to the explosion as predicted in our simulations. We use the term “initial mass” to indicate the zero-age main sequence mass in the case of single stars, MZAMS. In the case of binary systems, we define MZAMS as the mass of the initially more massive primary star, unless the SN originates from the initially less massive secondary star.

We first discuss the relation between the final core mass and the initial mass for single stars and how binaries affect this (Sect. 3.1). We then show the impact of different binary evolutionary channels (Sect. 3.2). We end with our predictions for the distribution of final core masses as we would expect from an idealized survey of observed Type II SNe events (Sect. 3.3).

3.1. Binary progenitors of Type II SNe typically originate from lower mass stars than single star progenitors

For single stars there is a tight relation between the initial mass of a star and its final core mass. This is because the evolution and final fate of a single star is primarily a function of its initial mass. Other parameters, such as for example, the metallicity and the initial rotation rate play a role, but are less important (e.g., Heger et al. 2003; Brott et al. 2011a). This statement breaks when considering extreme values, but such cases are very rare, certainly for events in the local universe. This means that for a given set of assumptions (for various uncertain internal physical processes, such as convective overshooting, rotational mixing, etc.), the initial-final core mass relation is effectively a one-to-one mapping.

The left panel of Fig. 1 shows the initial – final core mass relation for our fiducial simulation of single-star progenitors of Type II SNe. Under the assumptions followed in these simulations, Type II SNe result from stars with initial masses in the range between Mmin, ccSN ≈ 7.6 M and Mmin, WR ≈ 27 M. Stars that are initially less massive than Mmin, ccSN end their life as white dwarfs, whereas the ones more massive than Mmin, WR lose their hydrogen-rich envelope due to winds prior to explosion and become Wolf-Rayet stars. If they explode successfully, they would result in stripped-envelope SNe instead of the Type II SNe that we focus on here. The shape of the relation and the mass boundaries depend on the physical assumptions, in particular the treatment of internal mixing processes and stellar wind mass loss. However, for fixed assumptions, we always expect this relation between the initial mass and the final core mass to be tight for single-star progenitors.

thumbnail Fig. 1.

Initial mass to final core mass relation of SN II progenitors. We present the normalized 2D density distribution of the final core mass, MHe, preSN, with the initial mass, MZAMS, for all SN II progenitors in case they have an evolution in isolation (left panel) or in case they experience mass exchange with a binary companion (right panel). The color mapping is in logarithmic scale. We also depict the fitting formula for the one-to-one relation of single stars (Eq. (1); white dashed line) and for the median of the initial mass distribution for the case of binary progenitors (Eq. (2); red dashed line).

Mass exchange in binary systems can radically affect the pre-SN structure of a massive star (e.g., Podsiadlowski et al. 1992) and therefore breaks the one-to-one mapping between the initial stellar mass and the final core mass of a SN progenitor. This is illustrated in the right panel of Fig. 1, which shows a very broad distribution for binary progenitors.

We find that binary star progenitors of Type II SNe typically end their life with a larger core mass than they would have had if they had lived in isolation. This can be seen more clearly by comparing the distribution for binary progenitors with the white dashed line, which shows an analytical fit to the relation we obtained for single stars. This can occur, for example, after gaining mass from their companion through stable mass transfer or after a merger (as we show in more depth in Sect. 3.2). If such mass gain happens early in the evolution of the progenitor, we expect that its core can still grow as the star adapts its interior structure to its new mass (Hellings 1983, 1984; Braun & Langer 1995; Dray & Tout 2007). The star, therefore, ends its life with a core mass that is larger than if it would have been a single star.

The red dashed line shows a similar analytical fit to the median of the initial mass distribution for binary progenitors. The offset between the white and red dashed line emphasizes the systematic trend that binary progenitors tend to end their lives with larger core masses than expected for single stars with the same initial mass. Stating this more quantitatively: for a given final core mass, we find that binary progenitors of Type II SNe typically start their lives with an initial mass that is 15−25% lower than single stars. This implies that if one would not account for the possibility that a Type II SN progenitor experienced binary interaction, then one would typically overestimate its true initial mass.

The fact that almost no progenitors end their lives with core masses that are lower than expected according to the single star relation is a consequence of several effects. First of all, the mass of the core of a star is largely determined early in its evolution, when the star has completed its main sequence phase. After this a typical progenitor from a binary route seems unlikely to significantly reduce the mass of its core compared to a single star progenitor. Although stars may be spun-up by binary interactions, the transition between helium core and hydrogen-rich envelope is marked by a sharp chemical gradient, which largely inhibits rotational mixing that might otherwise reduce the mass of the core (e.g., Mestel 1957; Mestel & Moss 1986). However, in some cases other mixing processes can be effective in eroding the core of a merged star (e.g., Justham et al. 2014), which we neglect in this study. In principle, mass loss from the star could be also thought of as a mechanism to reduce the core size, but mass loss primarily affects the decoupled envelope. The core is not significantly affected unless the envelope is removed almost entirely and in this study we focus on Type II SNe, which still have at least part of their hydrogen envelopes left until the moment of explosion. Only for progenitors with low final core masses, with MHe, preSN ≲ 4 M, we find binary progenitors to lie below the single star relation as can be seen in Fig. 1. Most of these systems originate from mergers of two evolved stars, called reverse mergers, that we briefly discuss in the next subsection.

For convenience we provide the analytical fit to the one-to-one relation between the initial mass, MZAMS, and the final core mass, MHe, preSN for single stars, depicted by the white dashed lines in Fig. 1. We find that the single-star relation is well described (with a relative error of < 0.5%) by the second-order polynomial,

M He , preSN , singles M = 0.33 M ZAMS M + 0.002 ( M ZAMS M ) 2 0.59 , $$ \begin{aligned} \frac{M_{\rm He,preSN, singles}}{M_{\odot }} = 0.33 \frac{M_{\rm ZAMS}}{M_{\odot }}+ 0.002 \left(\frac{M_{\rm ZAMS}}{M_{\odot }}\right)^2-0.59, \end{aligned} $$(1)

in the range of 8 <  MZAMS/M <  25. Later in this study, we use the inverse of Eq. (1) to attribute to each pre-SN core mass an “equivalent single-star initial mass” that would have formed such a massive core if it was in isolation, even in cases that the progenitor in reality exchanged mass with its companion.

For binary progenitors the distribution is very broad, but if we want to describe the median of the distribution with a similar polynomial, we get:

M He , preSN , binaries M = 0.56 M ZAMS M + 0.002 ( M ZAMS M ) 2 1.67 , $$ \begin{aligned} \frac{M_{\rm He,preSN, binaries}}{M_{\odot }} = 0.56 \frac{M_{\rm ZAMS}}{M_{\odot }}+ 0.002 \left(\frac{M_{\rm ZAMS}}{M_{\odot }}\right)^2 - 1.67, \end{aligned} $$(2)

in the range of 7 <  MZAMS/M <  18, with a relative error of < 9% between the fitting formula and the median of the distribution. Equation (2) is shown as the red dashed line in the right panel of Fig. 1.

3.2. Role of different binary channels

The progenitors of Type II SNe originate from a variety of binary-evolution channels, which include stars that gain mass through Roche-lobe overflow and also various types of mergers. We presented an extensive discussion of these channels in Z19. Here, we briefly summarize and discuss how they contribute to the spread in the relation between the initial and final core mass.

To discuss the trends with mass, we choose to subdivide the progenitors into groups that each have similar final core masses. For each group we show the distribution of the initial masses of their binary progenitors. The group with the lowest final core masses, between 2 and 3 M, is shown in the bottom panel of Fig. 2. The panels above show the results for progenitors that end their life with higher final core masses, increasing in steps of 1 M, up to progenitors with final core masses between 9 and 10 M, which are shown in the top panel. The contributions of the different dominant binary channels are shown in different colors. The range of initial masses of single stars that would have reached the same core mass range is marked with a black dashed line. Its height has no physical meaning.

thumbnail Fig. 2.

Distribution of initial masses, MZAMS, of SN II progenitors experiencing binary interaction (red line) for increasing ranges of final core masses, MHe, preSN (on the right part). The contribution of each of the main binary scenarios is also shown in color below the distribution of binary-interacted progenitors. For reference we include the initial mass range of a single star corresponding to each MHe, preSN bin (black dashed line), which is shown in an arbitrary height. It becomes clear that the binary progenitors of SN II can originate from a wider range of lower on average MZAMS compared to single stars to form cores of similar mass, due to mass accretion or merging during their evolution. The results are for our standard assumptions although similar trends are found in all our model variations. The values on the y-axis show the number of SN II per bin in MZAMS per 106M of stars formed.

The figure highlights once more that binary progenitors with very similar final core masses originate from a wide range of initial masses. The main diagonal trend that is visible across all diagrams is caused by the fact that progenitors with higher final core masses (shown in the top panels) originate on average from initially more massive stars, both in the case of single stars and stars in binary systems. It is also important to note the decrease of the total number of SNe for higher core masses, mostly due to the initial mass function. The width of the distribution increases for systems with higher final core masses. This is because they typically result from more massive binary systems, where the amount of mass that a progenitor can gain is larger. The figure also shows that progenitors with very similar final core masses originate from very different evolutionary histories (see legend of Fig. 2).

An important binary channel that leads to Type II SN progenitors (Z19) consists of secondary stars that have gained mass through stable mass transfer (marked as mass gainers, shown in dark yellow). These are systems in which the binary is disrupted due to the prior explosion of the primary. This leaves the secondary as a single “walkaway” or “runaway” star (e.g., Renzo et al. 2019). Since it is not longer bound to a binary companion, it has the space to swell to giant dimensions at the end of its life without initiating a new phase of reverse mass transfer. This channel is important for progenitors that end their lives with high final core masses, as can be seen in the top panels. Typically they originate from stars that are initially half as massive as single stars. The contribution of this channel is smaller for progenitors with lower core masses, as can be seen in the bottom panels. This is because the accretor needs to become massive enough to explode as Type II SN and also due to the requirement that the primary star needs to be massive enough to explode in the first place and unbind the binary system. We also see in the bottom panels that the mass gainers of lower mass originate from a narrow initial mass range. This is because there is less mass available to accrete in these systems.

We also find a significant contribution of progenitors that result from channels that involve the merger of two stars. We distinguish whether it concerns the merger of two main-sequence stars (listed as MS+MS mergers, shown in orange), a post-main-sequence star with a main-sequence one (listed as postMS+MS mergers, shown in violet) or mergers where both stars are evolved and the merger is initiated by the evolutionary expansion of the secondary star (listed as reverse mergers, shown in brown). We caution that the evolution and final fate of merger products is particularly uncertain. For a discussion and overview of the limitations of the model predictions we refer to Z17 and Z19.

Mergers of the first type, MS+MS, occur early in the evolution of both stars. The resulting merger product is a star whose evolution may resemble that of a more massive single star (Schneider et al. 2019). It ends its life with a final core that is more massive than the individual stars would have had if they had evolved in isolation. These mergers play a more important role in the production of SN progenitors with more massive final cores, as can be seen in the top panels. This is because more massive stars are more often found in very close binary systems, where the stars start to interact already in the main sequence (e.g., Sana et al. 2012) and because massive stars expand more than lower mass ones during their main sequence, triggering the RLOF and the eventual merging.

Type II SN progenitors that are products of mergers of an evolved star with its main sequence companion (postMS+MS) originate from binary systems with the smallest difference in initial mass to single stars, compared to the other binary channels. The evolved primary star in this path has already formed a distinct helium core at the moment of merging and so we assume that the steep chemical gradient at the core-envelope boundary suppresses rotational mixing across the boundary (Mestel 1957; Mestel & Moss 1986). Further, the rate at which a post-merger helium core grows in mass in our models is mainly determined by the helium core mass itself. So, in our simulations, these progenitors lead to a pre-SN helium core mass similar to that of a single star only a few solar masses more massive than that of the initial primary-star mass. However, this neglects the possibility that the helium core mass may be reduced as a consequence of the merger, for example, by convective dredge-up during the post-merger relaxation phase (Justham et al. 2014). Indeed, models of the blue supergiant progenitor of the Type II-peculiar SN 1987A, which is thought to have been produced by this type of binary merger (Podsiadlowski et al. 1990; Podsiadlowski 1992), need to invoke extra mixing and so reduction of the core mass (Podsiadlowski 1992; Menon & Heger 2017). However, the core erosion found in the Justham et al. (2014) merger models is less significant for lower amounts of mass accretion during the merger process, and these are the mergers that eventually lead to red supergiant progenitors such as the ones in the observational sample to which we compare our models in Sect. 4.

Mergers of the third type, reverse mergers involving two evolved stars that have already formed a core, are favored among progenitors with the lowest final core masses as can be seen in the bottom panels. We also see in each panel that these mergers originate from stars with the lowest initial masses. This is because they typically result from binary systems where the primary is a star of intermediate-mass, too low to undergo core collapse in case it had an isolated evolution. If the primary had been more massive it would have been more likely to explode and create a compact object, possibly disrupting the system and anyway unable to produce another core-collapse supernova in case of a merger. The helium core or even the white dwarf formed from the primary needs to be engulfed by the evolving secondary that has gained mass beforehand. Because of their complex evolutionary history, with many instances of mass gain and loss for the two components, this is the only channel for which we find a non-negligible contribution that leads to lower MHe, preSN than if they were in isolation.

3.3. Initial mass and final core mass distribution for Type II SN progenitors

Getting accurate parameters from the observations of individual SNe is challenging, but the sample of well-studied events is increasing. This means that we are starting to get constraints on the distribution of the characteristics of SN progenitors, including estimates of their initial and final core masses. We therefore analyze our simulations to obtain predictions for the distributions of initial masses and of final core masses for Type II SNe.

We show the initial mass distribution of SN II in Fig. 3. For a population consisting only of single stars (black dashed line) the distribution is the initial mass function itself. All stars ranging between Mmin, ccSN ≈ 7.6 M and Mmin, WR ≈ 27 M lead to a hydrogen-rich event. The distribution of initial masses of progenitors of SN II that experienced binary interaction (red, combining all the possible binary channels) is more spread and it is shifted toward lower initial masses. A fraction of the progenitors originate from intermediate-mass stars, with MZAMS <  Mmin, ccSN, which form a core massive enough to collapse only due to merging with their companion (Z17). The distribution for binary progenitors of SN II for masses above Mmin, ccSN roughly follows the slope of the IMF but is lower in number of events. This is partly because the distributions are normalized to the total mass of stars formed and the total mass of a binary system is by definition higher in each mass bin of MZAMS (which depicts the initial mass of the one of the two stellar components that finally explodes). Another reason is that stars with MZAMS <  Mmin, WR can lead to stripped-envelope SNe instead of a Type II ones if they are in binary systems (shown in blue, stacked on top of the red, in Fig. 3). This is mainly due to mass stripping of the hydrogen-rich envelope of the primary star by its companion. Alternatively, for high initial masses close to Mmin, WR, the same can occur due to merging. If the merger product has eventually a strong enough wind mass loss to eject its hydrogen-rich envelope prior to explosion, it can avoid a Type II event. In summary, many stars in binaries in the mass range of ∼8 − 27 M produce a hydrogen-poor SN. At the same time the possibility of progenitors with initially masses lower than Mmin, ccSN partially compensate to the number of Type II events. The overall effect is the shift of the distribution of initial mass of Type II SNe to lower values.

thumbnail Fig. 3.

Distribution of initial masses, MZAMS, of SN progenitors with MZAMS <  Mmin, WR ≈ 27 M. For a population of single stars (dashed black line) we expect the distribution to follow the IMF, with all progenitors between Mmin, ccSN ≈ 7.6 M and Mmin, WR leading to Type II SNe. The shape of the distribution changes for progenitors that experienced binary interaction (solid black line). The progenitors of SN II arising from a population of binary interacting systems (red region) originate from lower on average MZAMS, even below Mmin, ccSN. In addition, many stars with MZAMS <  Mmin, WR lead to stripped-envelope SNe due to binary interaction (blue region, stacked on top of the red region), avoiding a hydrogen-rich event. The results are for our standard assumptions and show the number of SNe per bin in MZAMS per solar mass of star formation.

We also investigate the subsequent effect of binaries on the pre-SN core mass distribution, MHe, preSN. We show the results for our fiducial simulations for single and binary stars in Fig. 4 in black and red, respectively. On the top x-axis we also show the initial mass of a single star that would form a core as massive as shown in the bottom x-axis. To calculate this “equivalent single-star initial mass” we use the inverse of Eq. (1), independent of whether the actual evolutionary path of the progenitor was influenced by binary interaction.

thumbnail Fig. 4.

Distribution of the pre-SN helium core mass of SN II progenitors in case they experienced isolated evolution (black dashed line) or have followed a channel of binary interaction (red solid line). The contribution of each different binary scenario is also shown (dashed and/or dotted lines). The equivalent single-star initial masses corresponding to each MHe, preSN (inverse of Eq. (1)) is shown on the top x-axis. The results are for our standard assumptions and show the number of SN II per logarithmic bin in MHe, preSN per solar mass of star formation.

The core mass distribution that we obtain for binary progenitors is somewhat shallower than for single stars. We fit the distributions with a power-law dN/d M He,preSN M He,preSN α He $ {\rm d}N/{\rm d}M_{\rm He,preSN} \propto M_{\rm He,preSN}^{\alpha_{\rm He}} $, which is represents a line of slope (αHe + 1) in Fig. 4. We find a best fitting exponent αHe of around −2.0 for single stars, while we find a lower value of around −1.9 for binary progenitors. Both are shallower than the distribution of initial masses, which follows a standard IMF with a power-law slope of −2.3 (Salpeter 1955; Kroupa 2001). For single stars the difference is the effect of the relation between the initial and the final core mass (Eq. (1)). For binaries the resulting distribution is the combined effect of the four dominant channels of binary mass gainers and mergers that contribute. The individual contributions are shown with dashed and dotted lines in Fig. 4 and are discussed below.

The distributions for mass gainers and for MS+MS mergers stand out because they are very shallow, nearly flat, and even rising at the high mass end in case of the MS+MS mergers. These are also the channels in which the progenitor on average gains the largest amount of mass. In our fiducial simulations, these are the two channels that contribute most significantly to the production of Type II SN progenitors with high final core masses, as we have already shown in Sect. 3.2, see the top panels in Fig. 2. These two channels are thus the most effective in enabling stars with lower initial masses to end their lives appearing as if they had been initially much more massive. They are the main reason that the final core mass distribution of SN II is flatter for binaries, even though they originate on average from lower initial masses (red in Fig. 3).

We see a different behavior for mergers of the second type (postMS+MS). They have a slope that is roughly similar to the slope we obtain from our single stellar evolutionary models. The reason for this is that the merger occurs after the evolved star has formed a well-defined core which afterwards does not change significantly. Thus, in our models this channel of merger has little effect on the final core mass distribution.

Finally, reverse mergers show the steepest core distribution, peaking at low final core masses. This is because they primarily originate from intermediate-mass binary systems where the primary star is close to or below Mmin, ccSN ∼ 7.6 M and does not explode prior to merging. Thus, they tend to form SN progenitors with final core mass equivalent to single stars with initial mass not much above the minimum mass for core collapse.

The slightly shallower helium core mass distribution when accounting for all binary channels implies a larger relative contribution of stars with massive final cores, as compared to a population of single stars with our assumed IMF. One way to visualize this is by how flat the IMF would need to be for a purely single-star population to produce such a core-mass slope. In Table A.1 we compute the IMF slope, αsingle_IMF, that would have led a population of single stars to generate a final helium core-mass distribution that best fits the outcomes of each of our model binary populations. We further discuss these model variations in the next section, when comparing with observations.

4. Implications for interpreting direct detections of Type II SN progenitors and uncertainties

4.1. The sample of Type II SN progenitors

In this section we place our results in the context of observational constraints from direct progenitor detections. We discuss the implications for different observational techniques in Sect. 5.

Various groups have devoted large efforts to these searches, mostly using Hubble Space Telescope archival images (e.g., Van Dyk et al. 1999, 2003; Smartt et al. 2001, 2009; Maund & Smartt 2005; Kochanek et al. 2012; Davies & Beasor 2018). Smartt (2015) made a compilation of SN progenitor observations. After briefly describing this sample here, we compare with our predictions in Sect. 4.2, and discuss the effects of model variations and potential physical implications in Sects. 4.3 and 4.4.

We use the sample of nearby Type II SNe compiled by Smartt et al. (2009) and Smartt (2015). This is a time- and volume-limited sample of SNe for which searches for progenitors were feasible. It contains 26 Type II SN events. For each event high-resolution images taken before the explosion allowed direct detection of the SN progenitor, or could be used to calculate upper limits on the progenitor’s luminosity. Smartt et al. (2009) compared the photometric measurements (or upper limits) with predictions for the pre-explosion luminosities from single stellar models by Eldridge & Tout (2004).

The final core masses for their best fitting progenitor models are not quoted in the paper by Smartt et al. (2009) directly, but they can be readily reconstructed using their Eq. (1). The result of this is shown as square symbols in Fig. 5, where the horizontal error bars indicate confidence intervals. This diagram is the analog to Fig. 6 in Smartt (2015), except that we chose to replace the initial mass with the final core mass on the bottom horizontal axis. The latter has the important advantage that it is more directly linked to what is inferred from pre-SN observations, circumventing uncertainties in the initial-final core mass relation. This is particularly important when considering binary progenitors, for which there is no unique mapping between the initial and final core mass, as we showed in Sect. 3.1.

thumbnail Fig. 5.

Cumulative distribution of final core masses, MHe, preSN, of Type II SN progenitors. White squares are the inferred core mass or upper limits of observed SN II progenitors from Smartt (2015). For comparison we show our predicted cumulative distribution for progenitors that evolved as single stars (black dashed line) or including possible binary interactions (red solid line), following our standard assumptions. The equivalent single-star initial mass, MZAMS, corresponding to each MHe, preSN is shown on the top x-axis. The equivalent pre-SN progenitor luminosity, corresponding to each MHe, preSN value according to single STARS models which were also used in Smartt et al. (2009), is shown on the bottom x-axis. We also show the range of uncertainty (red-shaded region with white dashed lines) of the expected distributions found for our simulated populations in which we vary the treatment of binary physics. These include among others the efficiency of mass transfer and of common envelope evolution, the initial period and mass ratio distribution and the binary fraction.

In the remainder of this section we discuss and overplot the helium core mass distributions derived from our simulations. When comparing our predictions with the observations, it should be kept in mind that the core masses of SN progenitors are inferred from the progenitor luminosity using single-star models. We are mostly confident in this approach for the binary progenitors which experience interaction in an early evolutionary stage, which are expected to continue their further evolution in a fashion that is similar to that of a single star. At present it is unknown exactly how accurate the estimates are for progenitors that result from more complex binary channels; our fiducial simulations suggest a contribution of 11% for reverse mergers, and 14% for mergers involving one postMS star, with both paths probably leading to more complex evolution than stars in isolation. A subset of these may lead to pre-explosion structures that differ from those that can be achieved by stars that evolve in isolation, for example those that end their lives as blue supergiants (e.g., Podsiadlowski 1992; Justham et al. 2014; Menon et al. 2019). Blue supergiant SN progenitors do not follow even an approximate core-mass-luminosity relation, unlike red supergiants (see, e.g., Fig. 3 of Justham et al. 2014; for red supergiants see also, e.g., Farrell et al. 2020). Overall we expect that this method should at least be reasonably reliable for stars which die as red supergiants, as in this observational sample.

This complexity of possible binary interaction of the progenitor is added on top of other sources of uncertainties. In part there are caveats that are inherent to the observations and the possible role of selection effects, such as dust extinction of detected red supergiant progenitors (e.g., Walmswell & Eldridge 2012; Beasor & Davies 2016) or limited accuracy of the exact position of the progenitor in the Hertzsprung-Russel diagram (Davies & Beasor 2018). There are also uncertainties in the link between the inferred pre-SN mass of the progenitor star and its ZAMS mass that exist anyway in single stellar evolution, including the mass-loss rate through winds (e.g., Smith 2014; Renzo et al. 2017), the overshooting of the convective core (e.g., Ribas et al. 2000; Claret 2007; Brott et al. 2011b; Straniero et al. 2019) and the theoretical preSN luminosities of the progenitors (e.g., Farrell et al. 2020).

4.2. Comparison with our fiducial models

In Fig. 5 we overplot the cumulative final core mass distributions for hydrogen-rich SN progenitors from our simulations. In black we show the results for our simulations for single stars and in red we show the distribution for our fiducial simulation for binary star progenitors. Comparing the black and red line, we see that accounting for binary interaction does not have a large effect on the distribution. The distribution for binary progenitors is slightly shallower than that for single stars, as we showed in Sect. 3.3.

The data roughly follows the theoretical distribution for the lower mass end, perhaps depending on the interpretation of the upper limits, but one could argue that there is a deviation between the model predictions and the data at the high mass end. The observed distribution appears to be steeper (or alternatively to have a lower maximum mass) than predicted by our models. In other words, our models over-predict the number of hydrogen-rich stars which reach core collapse with high helium core masses, when compared to the observed sample of Type II SN progenitors. This is the same discrepancy as noted by Smartt et al. (2009), who coined the term “red supergiant problem” for the apparent lack of detections of Type II SNe resulting from progenitors with an initial mass higher than ≈16.5 M. This corresponds to a helium core mass of ∼5.75 M (via Eq. (1) of Smartt et al. 2009).

Our simulations show that the inclusion of the effects of binary interaction does not help to solve or alleviate the claimed red supergiant problem. If anything, for our standard set of assumptions it makes it slightly worse. To quantify this effect (still focusing only on stars which retain their hydrogen-rich envelope), we calculate the fraction of potential Type II SN progenitors, Rmassive, with a core more massive than 5.75 M. Our models including binary systems predict a slightly larger Rmassive, about 0.2, compared to 0.18 for a single-star population (Table A.1). If stars with core masses above that value result in “failed” events directly collapsing onto a black hole (as suggested in Smartt et al. 2009 as a possibility and further discussed in Sect. 4.4), Rmassive would represent the fraction of hydrogen-rich stars at core collapse that would not produce a SN event.

An alternative way to compare our models to the observational sample is to calculate the probability for each model population, P0, massive, that in an observed sample of 26 SN II progenitors, none would have been found with a helium core more massive than 5.75 M. That is the same as the chance of picking randomly 26 SN progenitors from each model population (which is weighted according to initial conditions and taking into account stellar and binary evolution) and finding zero stars with equivalent single-star initial mass higher than 16.5 M. These are presented in Table A.1 for all our simulations. Even given the caveats above, we consider this to be a helpful indication. The probability of zero progenitor detections with helium core masses above 5.75 M, in a sample of 26 trials, is roughly a factor of two lower for our fiducial binary model than for purely single stars.

It is important to note however that this probability calculation is still based on the mass plane of the SN progenitors (of their cores in particular). Davies & Beasor (2020) recently, by redoing the analysis of the claimed lack of high-mass progenitors from the point of view of observed red supergiant luminosities, argue that it may be an effect of the small sample size and the steep luminosity distribution of the SN progenitors, and thus of lower statistical significance. In our work we refrain from directly predicting pre-SN luminosity distributions due to all the extra theoretical and observational uncertainties that our results would be subjected to, as discussed in Sect. 4.1.

4.3. Impact of assumptions concerning the treatment of binary physics

Our predictions are subject to several uncertain assumptions. Here we discuss how variations in these assumptions affect the final core-mass distribution. For a full list of the model variations that we considered, see Table A.1. For further discussion of the model variations we refer to Z17 and Z19.

Our results are not very sensitive to assumptions concerning the treatment of binary physics or the initial distribution of mass ratios and orbital periods. This can be seen from the red-dashed region in Fig. 5, which shows the collective area spanned by model variations (these are listed in Table A.1 as the first group). Each variation in binary physics assumptions is shown with white dashed lines inside the region. Most of the simulations result in shallower distributions for the final core mass than our results for single stars. None of the variations considered here solves the discrepancy between observations and the models at the high mass end. For these variations, the probability of not detecting a more luminous progenitor than in the Smartt (2015) sample remains below P0, massive ≲ 10−2, typically close to P0, massive ≈ 10−3, as for our fiducial model and purely single-star populations (also shown in Table A.1).

Simulations that predict a steeper cumulative core mass distribution are variations where we adopt extreme assumptions. These include a model variation where we assume that no mass is accreted during Roche-lobe overflow (Model 01), which prevents secondary stars from gaining mass, and a model variation in which we ignore the natal kick of compact objects (Model 10), which drastically reduces the number of binary systems that are disrupted upon the first SNe. Such disruption is needed to prevent the secondary from being stripped by the primary star and therefore reduces the number of secondaries that end their lives as a Type II SN. Finally, considering low efficiencies for common-envelope ejection (Models 12–14) leads to an increase in the contribution of reverse mergers, which dominate our predicted population of Type II SNe events with low final core masses. This also results in a steeper distribution in final core masses. The model variations that lead to the shallowest distributions include those where we assumed a very high accretion efficiency (Model 03) or high efficiency for common envelope ejection preventing the contribution of mergers (Models 15 to 17).

4.4. Impact of failed SNe, minimum mass for core-collapse SNe, stellar winds and the initial mass function

Since the inclusion of binary interaction does not significantly alter the final core mass distribution, we investigate which other model assumptions can potentially alleviate the apparent discrepancy between the predicted distribution for SN II and the observational data, known as the red supergiant problem.

Failed SNe. In the last decade, there has been a debate about whether some stars do not produce a transient event because of full collapse onto a black hole (O’Connor & Ott 2011; Piro 2013; Ertl et al. 2016). This is linked with the details of the explosion mechanism (e.g., Janka 2012), subsequently determining the compact remnant that they leave behind (e.g., Sukhbold et al. 2016). Furthermore, Gerke et al. (2015), Adams et al. (2017a) reported the disappearance of a red supergiant and proposed that this is an example of such an event. The possibility of failed SNe was suggested by Smartt et al. (2009) as a potential solution to the apparent red supergiant problem.

We thus consider the simple possibility that cores above a mass threshold lead to a collapse with no visible transient. In the upper-left panel of Fig. 6 we show the effect of variations in which we exclude all SNe that result from a progenitor with a helium core more massive than MHe, preSN = 4.8, 5.5, 6 and 6.8 M. These core masses correspond to a single-star initial mass of Mmax, ccSN ≈ 15, 16.5, 18 and 20 M, respectively, where we used Eq. (1) in this work (Models 52–54 and 29). The small differences in the final core masses compared to Smartt et al. (2009) are due to the slightly different initial – final core mass relation for single stars in that study from our Eq. (1). Decreasing Mmax, ccSN shifts the cumulative distribution toward lower final core masses, changing the shape of the predicted curve to be closer to the values inferred from the observational sample. This can lead to probabilities P0, massive of order unity for not finding SN progenitors with more massive cores than the observed sample, as we exclude high mass core stars from our progenitor population (see Table A.1).

thumbnail Fig. 6.

Same as Fig. 5 (with logarithmic x-axis here), but showing the results for other varying assumptions (blue lines): (a) the mass threshold for failed SN collapsing onto black holes Mmax, ccSN (Models 52, 53, 54 and 29), (b) the wind mass loss efficiency η (Models 25, 26 and 51), (c) the minimum mass for SN, Mmax, ccSN (Models 30 and 31) and (d) the slope of IMF αIMF (Models 32 to 34). These physical parameters influence the whole population, both single stars and in binary systems. The list of values depicted in the upper left corner shows a sequence of the varying distributions (blue lines) spanning from top left to bottom right of each panel.

Earlier we referred to the calculated fraction of hydrogen-rich stars at core collapse which produce failed SNe, Rmassive, assuming that the helium core mass governs this outcome (Table A.1). For most models we use a threshold helium core mass of ∼5.75 M which corresponds to that of an initial mass of 16.5 M for a single star, following Smartt et al. (2009). This fraction of failed hydrogen-rich core-collapse events is around 0.2 for both our population of single stars and our fiducial simulation, as well as for most of our model variations. This ratio is consistent with up-to-date results of surveys that monitor for observational evidence of this phenomenon (Kochanek et al. 2008; Gerke et al. 2015; Adams et al. 2017a,b). The value is only moderately influenced even for extreme assumptions in the IMF slope, wind efficiency and metallicity. Models 29 and 52–54 consider different values of this parameter, Mmax, ccSN, which naturally cause the predicted ratio to change.

Stellar-wind mass loss. Altering the efficiency parameter of wind mass loss, η, results in significant variation in our predictions (see the upper-right panel of Fig. 6; Models 25, 26, 5, 52, 53 and 54). For lower wind efficiencies, the cumulative distribution is becoming even shallower, increasing the discrepancy with observations. This is because stars are less affected by winds and thus even progenitors with MZAMS ≳ 30 − 35 M are able to produce Type II SNe. The reason that these distributions have very small differences is that stars above MZAMS ≳ 38 M reach the Humphreys-Davidson limit (Humphreys & Davidson 1979) and get stripped of their hydrogen-rich envelope anyway, even for low values of η.

Increasing the mass-loss rate by a factor of 1.5 to 3 (Model 26, 53, 54) leads to much steeper distributions due to the high wind mass-loss rates that cause progenitors above single-star masses of 20−13 M, respectively, to lose their hydrogen-rich envelope. They end their lives as a stripped-envelope SN instead of a Type II SN. This also severely affects the relative rates of those subclasses (Z19). Although the distribution for η = 1.5 seems visually quite consistent with observations, it should be primarily seen as a test for extreme assumptions of high mass loss, because, if anything, wind mass-loss efficiencies may be lower than in our standard assumptions (Smith 2014). Changing the metallicity (Models 39–44, not shown in Fig. 6) has a similar effect, since this also modifies the wind mass loss rates (e.g., Vink et al. 2001; Mokiem et al. 2007). Our model predictions for a super-solar metallicity of Z = 0.03 appear to provide a good match to the data. However, we have no indications, for example from the host galaxies, that the majority of the sample consists of SNe originating from metal-rich progenitors.

IMF slope. Varying the slope of the IMF (Models 32–34; bottom-right panel of Fig. 6) does not strongly affect the shape of the final core-mass distribution. Even an extremely steep IMF (e.g., αIMF = −3.0 of Model 34) only increases the probability of missing a higher-mass SN II progenitor by chance to P0, massive ∼ 2 × 10−2.

Minimum mass threshold for core-collapse SN. In models 30 and 31 (bottom-left panel of Fig. 6), we change the minimum core mass for core-collapse SNe. A lower minimum threshold leads to a distribution which rises earlier in mass and thus sightly decreases the expected contribution of high mass progenitors above the observationally established mass limit.

The parameters described in this subsection govern the initial conditions and the evolution of the whole population, both single stars and binary systems. The sensitivity of our results to these assumptions is almost independent of whether we include binary systems or not in our simulations. The impact of these parameters on the final core mass distribution requires further careful inspection but this is not the aim of this paper in which we focus on the consequences of binary interactions.

5. Predicted discrepancy of inferred initial masses of Type II SN progenitors between different observational methods

It is still challenging to establish a direct link between a SN event and its progenitor star. To this end, significant work has been invested in alternative methods to probe the characteristics of the progenitors of Type II SNe, in addition to searching for the progenitor in pre-SN images.

A way to directly constrain the SN progenitor is to study the late-time, nebular-phase spectra of the SN when the ejecta become optically thin (Jerkstrand et al. 2012). This provides an estimate of the oxygen core mass of the stellar progenitor. A few studies have implemented this method, finding again a lack of high-mass progenitors (e.g., Jerkstrand et al. 2014; Tomasella et al. 2013). This is consistent with the red supergiant problem inferred from direct progenitor detections, discussed in Sect. 4 (although see Anderson et al. 2018 for a reportedly high-mass progenitor case in a low-metallicity environment). Studies of the SN light curve evolution, that primarily constrain the ejecta masses of the event, are in general agreement with the above studies (e.g., Valenti et al. 2016; Martinez & Bersten 2019; Eldridge et al. 2019).

All the techniques mentioned above (including the pre-SN imaging of the progenitors) directly probe the state of the star very close in time or at the moment of explosion. These studies put important constraints on the final progenitor state, but we should note that they do not directly probe the initial mass of the SN progenitor. In fact, we suggest avoiding translating the SN progenitor properties on the initial mass plane. Firstly, because it propagates the uncertainties of stellar evolution from formation up to explosion (e.g., Davies & Beasor 2020; Farrell et al. 2020, see also Sect. 4.1). And especially because, in the case of binary interactions, the initial mass loses its predictive strength about the end result, due to the diversity of the different binary paths. We argue that the methods that estimate the pre-SN core mass and the stellar luminosity of the SN progenitor are more closely related to what is defined as the final, pre-SN core mass in this study. Thus, if the possibility of binary mass gain or merging is not taken into account, we expect an overestimation of the inferred initial masses of the SN II progenitor stars from these techniques (Sect. 3.1).

An independent way to probe the initial mass of the progenitor system is to infer the delay time between the formation of a star and its SN explosion, by studying the environment of the progenitor. This is achieved by effectively “age-dating” the parent population of the progenitor and thus inferring its initial mass. This method can mainly be applied to nearby SNe (Murphy et al. 2011; Williams et al. 2014, 2018). As the local rate of SNe is low, the same technique has been used for SN remnants too (Maoz & Badenes 2010; Jennings et al. 2014; Maund 2017; Díaz-Rodríguez et al. 2018; Auchettl et al. 2019). Several other studies investigate the association of SN II with star-forming regions, through observations of nebular emission lines of the SN environment (e.g., Anderson et al. 2012; Habergham et al. 2014; Xiao et al. 2018; Kuncarayakti et al. 2018; Schady et al. 2019), thus also providing an estimate of the initial mass of the progenitor star. This technique has been also used to other stellar objects, like Luminous Blue Variables that may be progenitors of type IIn SNe. Smith & Tombleson (2015) and Aghakhanloo et al. (2017) argue toward a binary interaction history of these objects to explain their relative isolation.

The lifetime of a star is not strictly determined by its initial mass as it may be affected by mass exchange with a companion (Crawford 1955; De Donder & Vanbeveren 2003; Z17; Xiao et al. 2019). However, the initial masses of the stars in a binary determine the approximate timescale until binary interaction, which is usually the dominant component of the full time between formation and explosion. Thus, within uncertainties, the “age-dating” techniques estimate the actual initial masses of the SN progenitors (even the ones that experienced binary interaction). As we show in Fig. 3, the expected initial mass distribution of SN II progenitors is shifted to lower initial masses (and longer delay-times, Z17) than a canonical Salpeter IMF. Indeed, Jennings et al. (2014) find a discrepancy in the inferred initial mass distribution of SN II progenitors compared to a Salpeter IMF. They report that there may be a bias against the highest mass progenitors in their sample, but we argue that this result may be partially an outcome of the possible binary history of many of the progenitors.

Putting all these together, we argue that the two groups of observational methods discussed (the ones that study directly the SN progenitor and the others that examine its surrounding environment) probe different characteristics of the system that exploded and therefore potentially infer different initial masses. In Fig. 7 we depict the theoretically expected 2D histogram of the inferred initial masses from the two techniques, for a realistic population of 50% single stars and 50% binary systems (e.g., Sana et al. 2012; Duchêne & Kraus 2013; Moe & Di Stefano 2017). We calculate the MZAMS that age-dating methods would derive (x-axis) by computing the delay-time of the SN progenitors in our simulations and then transforming it to an initial mass of a single star that would have the same lifetime (for details see Fig. 1 of Z17). At the same time, to compute the MZAMS that direct observations of the SN progenitors would infer, we transform the simulated final helium core to an initial mass, using the inverse of Eq. (1) in this study.

thumbnail Fig. 7.

Discrepancy in the inferred initial masses between methods that probe the pre-SN core mass or luminosity of the progenitor vs. the ones that age-date its surrounding environment. Our results (2D probability density distribution in green) show the transformation of the simulated delay-time of the SN progenitor to MZAMS (x-axis), assuming a lifetime-mass relation for single stars (Fig. 1 of Z17). The y-axis is the transformation of the simulated MHe, preSN to MZAMS using the inverse of Eq. (1). The points are nearby Type II SNe with observed estimates of MZAMS with both groups of techniques, shown in Table 1.

All single stars have the same value of MZAMS from the two methods (and are found on top of the x = y line in Fig. 7). However, we predict that for a significant fraction of SN II, the initial mass of a progenitor inferred from age-dating the surrounding population is systematically lower than inferred from direct observational methods that study pre-SN properties. This is because a star may take long to explode due to its initial low mass but it has a more massive core at the moment of its explosion than corresponding to a single-star scenario of the same initial mass (due to merging or mass accretion during its lifetime). This discrepancy is potentially observable both for specific events and in statistical samples.

We compare our results with a sample of Type II SNe that have estimates of their initial masses with both methods. The sample can be found in Table 1. These events have estimates of their MZAMS based on pre-SN imaging of their progenitor (Smartt 2015; Davies & Beasor 2018; Van Dyk et al. 2019) or modeling of the SN nebular spectrum (Silverman et al. 2017). We made the argument that these techniques probe more closely the pre-SN core (corresponding to the y-axis in Fig. 7). At the same time, these SNe have independent estimates of their MZAMS, based on studies of their surrounding environment that constrain the age of the progenitors (corresponding to the x-axis in Fig. 7).

Table 1.

Sample of Type II SN progenitors with inferred initial masses from both directly probing the SN progenitor and age-dating the surrounding environment.

Part of the sample has been presented in Schady et al. (2019), where they used Hα emission lines in the surrounding region as an age tracer. They show that this method is highly sensitive to the modeling of the stellar environment, because binary products make a region look younger than its actual age. This can lead to a possible underestimation of the delay-time and thus an overestimation of the initial mass of a SN II progenitor. Schady et al. (2019) address the reported high initial masses derived in Kuncarayakti et al. (2018) where single-stellar populations have been assumed for the age-dating (2nd column in Table 1, in parenthesis), arguing toward lower initial masses instead if binarity is taken into account (3rd column in Table 1). In Fig. 7 we show their derived values where the full effect of binarity on the time evolution of Hα emission lines is taken into account (blue squares in Fig. 7). They argue that these values of MZAMS should be considered lower limits because possible photon leakage decreases the effect of binarity in the age-dating via Hα emission lines (Xiao et al. 2018). The observed values of MZAMS, where the effect of binarity in the age-dating is considered, are roughly in-line with our prediction of the discrepancy between the two groups of observational methods. In fact, for some events the initial masses from studying the environment are much lower than our theoretical predictions, although they should be considered lower limits anyway.

Williams et al. (2018) have also compared initial mass estimates from the two different method groups for a sample of SN events. We include these events in our sample too, with updated MZAMS estimates from pre-SN images by Davies & Beasor (2018) and Van Dyk et al. (2019) for most of them. Williams et al. (2018, adding to previous work from Williams et al. 2014) derive an estimate of the age of the region by modeling the resolved stellar population around the SN (2nd column in Table 1, not in parenthesis). They report a broad agreement between their initial mass estimates and the ones from pre-SN images (red stars in Fig. 7). So, there is no obvious trend for SNe that show a discrepancy between these two methods, although SN2017eaw is a clear exception. Williams et al. (2014, 2018) used single stellar modeling, not taking into account the effect of binarity on the observed color-magnitude diagram of a stellar population. Binaries would lead to hotter and more luminous stars, potentially leading to underestimations of a population’s age (e.g., Eldridge & Stanway 2009; Schneider et al. 2015). At the same time, they considered only the recent star-formation history of a region up to 50 Myr before the SN, based on the maximum lifetime of a single star progenitor. De Donder & Vanbeveren (2003) and Z17 show that we should expect even higher delay-time for SN progenitors originating from intermediate-mass stars. Thus, allowing for the possibility of longer delay-times than 50 Myr may result in lower derived initial progenitor masses. It would be interesting to see the results if those effects are considered, as both of them qualitatively are expected to lower the MZAMS estimates of these studies.

Although the theoretical and observational uncertainties in both groups of methods are significant, it is interesting that some observational points are found above the x = y axis. As a case study, we can discuss the significant discrepancy of the SN 2017eaw initial mass values. Here we used MZAMS  ∼  15 M derived by Van Dyk et al. (2019) from pre-SN images, but Kilpatrick & Foley (2018) and Rui et al. (2019) also infer an initial mass around 12 − 13 M. Interestingly, all of them are higher than the estimate of 8 . 8 0.2 + 2.0 $ 8.8^{+2.0}_{-0.2} $ from age-dating of Williams et al. (2014). Conclusive statements are difficult to make without strong constraints from both of the groups of techniques, but a high discrepancy in the values derived from these methods can potentially be an indicator of a binary history of a Type II SN progenitor. This indication will be even stronger in case it is combined with an exceptionally long delay-time of the progenitor (> 50 Myr) which would correspond to a single star initial mass below the mass threshold for a collapse Mmin, ccSN ∼ 7.5 M. These late SNe can be a signature of binary interaction of the progenitor too (Z17). A bigger sample of Type II SNe with accurate estimates of the pre-SN core or luminosity and, at the same time, with careful study of the age of their environment will shed more light on possible binary history of their progenitors.

On the other hand, it is interesting to discuss a few events that are found well below the x = y line, which is not expected in our findings. For example, the 9.5+0.7M upper limit of SN 2004dg in pre-SN images seems to be lower than the one derived by age-dating, even when a surrounding population of binaries is considered. Given our theoretical expectations, a physical reason for this to occur would be that the progenitor was heavily obscured and therefore its inferred mass from the nondetection of it is an underestimate. Further studies of these group of unexpected events would be needed to shed more light on their evolutionary history.

In our analysis we did not include the SNe reported in Williams et al. (2014, 2018) and Schady et al. (2019) where the progenitor was fully or partially stripped of its envelope (Type Ib iPTF13bvn, and Type IIb’s SN1993J and SN2011dh) because our predictions apply for hydrogen-rich, Type II SNe only. We also excluded SN1987A because single star models struggle to simulate its blue supergiant progenitor, thus they cannot confidently infer what its initial mass would be if it had somehow evolved from a single star.

We note that a potential complication for the methods of age-dating is the possibility that the SN progenitor was ejected from its initial binary system prior to its own explosion. This is exactly the case for SN II progenitors that gained mass, according to our study. These stars can travel tens to hundreds of pc away from their birthplace until explosion (Eldridge et al. 2011; Renzo et al. 2019). We expect a statistically weaker correlation of SNe from mass gainer progenitors to young, star-forming regions compared to the bulk population of SNe II. The relative spatial isolation of the progenitor of SN2017eaw (Van Dyk et al. 2019) also fits to this possibility.

6. Summary and conclusions

This work investigates the effect of binary interactions on the initial and final core mass distribution of hydrogen-rich Type II supernova progenitors. We expand on Zapartas et al. (2019), where we find that around 1/3 to 1/2 of Type II SN progenitors have experienced mass accretion or merged with a companion prior to explosion.

  • We demonstrate that binary interactions lead to a broad relation between initial mass and final core mass because of the variety of possible evolutionary paths. Consequently, this breaks the tight relation between pre-explosion properties and initial stellar mass.

  • We find that progenitors that have experienced mass accretion or merging with a companion, for a given helium core mass, typically originate from lower initial masses than would be expected from single-star evolution (Figs. 1 and 2). We argue that pre-SN images of Type II SN progenitors give a more direct insight into the final helium core masses than total masses, since the final stellar luminosity is more closely related to the former (e.g., Giannone 1967). Thus, if one does not consider the possibility of a binary history for a Type II SN progenitor, it may often lead to overestimation of its inferred initial mass.

  • We find that the initial mass distribution of SN II progenitors is broadly shifted to lower initial masses than expected from a canonical Salpeter IMF of single stars (Fig. 3). Our fiducial binary population simulation produces a marginally shallower distribution of final helium core masses than from a population of purely single stars. This is best fit by a power-law slope of approximately −1.9, compared to about −2.0 for a single-star population (Fig. 4).

  • Binary interactions do not appear to significantly affect the red supergiant problem, that is, the reported lack of SN progenitors with final inferred helium cores that correspond to initial single-star masses above ∼16.5 M (Smartt et al. 2009; Smartt 2015). This seems to hold for all our realistic variations in binary physics assumptions (Fig. 5).

  • We also show the effect of other physical processes and conditions that affect the final core distribution (Fig. 6). The inclusion of failed SNe or the increase of the typically used wind mass-loss rates result in a steeper distribution skewed toward lower core masses.

  • We argue that the initial masses inferred from the surrounding environments of Type II SN or SN remnants are typically closely associated with the actual initial primary-star mass of the progenitor, even for binary progenitors, as long as binarity is properly taken into account when age-dating the region. On the other hand, the luminosity of a red supergiant progenitor would correlate most strongly with its final core mass, and inferences based on the pre-SN state would tend to overestimate the initial stellar mass of the progenitor if they ignore the possibility of mass gain or merging (Sect. 5). Thus, we predict the inferred initial masses from studies of the surrounding environment to be systematically lower than inferred in the techniques that estimate the pre-SN core or luminosity (Fig. 7), consistent with the new observational analysis of Schady et al. (2019). An accurate and significant discrepancy in the values derived from these methods for a particular Type II SN may even suggest a binary history for the progenitor of that SN.

We suggest that future studies of Type II supernova progenitor populations and explosion diversity should not assume a one-to-one mapping between final progenitor properties and initial stellar mass. This suggestion applies even where individual cases do not seem to require an explanation involving binary interactions.

Acknowledgments

We are grateful to Rob Izzard for providing the population synthesis code, binary_c, used in this study and for his technical support. We are also grateful for the very useful discussions with Katie Auchettl, Ben Davies, Maria Drout, Tassos Fragos, Avishai Gal-Yam, Ylva Götberg, Avishai Gilkis, George Meynet, and Silvia Toonen. This project was funded in part by the European Union’s Horizon 2020 research and innovation program from the European Research Council (ERC, Grant agreement No. 715063), and by the Netherlands Organization for Scientific Research (NWO) as part of the Vidi research program BinWaves with project number 639.042.728. EZ acknowledges support from the Swiss National Science Foundation Professorship grant (project number PP00P2 176868) and from the Federal Commission for Scholarships for Foreign Students for the Swiss Government Excellence Scholarship (ESKAS No. 2019.0091).

References

  1. Adams, S. M., Kochanek, C. S., Gerke, J. R., Stanek, K. Z., & Dai, X. 2017a, MNRAS, 468, 4968 [NASA ADS] [CrossRef] [Google Scholar]
  2. Adams, S. M., Kochanek, C. S., Gerke, J. R., & Stanek, K. Z. 2017b, MNRAS, 469, 1445 [NASA ADS] [CrossRef] [Google Scholar]
  3. Aghakhanloo, M., Murphy, J. W., Smith, N., & Hložek, R. 2017, MNRAS, 472, 591 [NASA ADS] [CrossRef] [Google Scholar]
  4. Almeida, L. A., Sana, H., Taylor, W., et al. 2017, A&A, 598, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Anderson, J. P., Habergham, S. M., James, P. A., & Hamuy, M. 2012, MNRAS, 424, 1372 [NASA ADS] [CrossRef] [Google Scholar]
  6. Anderson, J. P., Dessart, L., Gutiérrez, C. P., et al. 2018, Nat. Astron., 2, 574 [NASA ADS] [CrossRef] [Google Scholar]
  7. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  8. Auchettl, K., Lopez, L. A., Badenes, C., et al. 2019, ApJ, 871, 64 [NASA ADS] [CrossRef] [Google Scholar]
  9. Baade, W., & Zwicky, F. 1934, Proc. Natl. Acad. Sci., 20, 254 [NASA ADS] [CrossRef] [Google Scholar]
  10. Beasor, E. R., & Davies, B. 2016, MNRAS, 463, 1269 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bethe, H. A., Brown, G. E., Applegate, J., & Lattimer, J. M. 1979, Nucl. Phys. A, 324, 487 [NASA ADS] [CrossRef] [Google Scholar]
  12. Blaauw, A. 1961, Bull. Astron. Inst. Neth., 15, 265 [NASA ADS] [Google Scholar]
  13. Braun, H., & Langer, N. 1995, A&A, 297, 483 [NASA ADS] [Google Scholar]
  14. Brott, I., Evans, C. J., Hunter, I., et al. 2011a, A&A, 530, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Brott, I., de Mink, S. E., Cantiello, M., et al. 2011b, A&A, 530, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cao, Y., Kasliwal, M. M., Arcavi, I., et al. 2013, ApJ, 775, L7 [NASA ADS] [CrossRef] [Google Scholar]
  17. Chini, R., Hoffmeister, V. H., Nasseri, A., Stahl, O., & Zinnecker, H. 2012, MNRAS, 424, 1925 [Google Scholar]
  18. Claeys, J. S. W., de Mink, S. E., Pols, O. R., Eldridge, J. J., & Baes, M. 2011, A&A, 528, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Claret, A. 2007, A&A, 475, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Conti, P. S. 1975, Mem. Soc. Roy. Sci. Liege, 9, 193 [Google Scholar]
  21. Crawford, J. A. 1955, ApJ, 121, 71 [NASA ADS] [CrossRef] [Google Scholar]
  22. Davies, B., & Beasor, E. R. 2018, MNRAS, 474, 2116 [Google Scholar]
  23. Davies, B., & Beasor, E. R. 2020, MNRAS, 493, 468 [NASA ADS] [CrossRef] [Google Scholar]
  24. De Donder, E., & Vanbeveren, D. 1998, A&A, 333, 557 [NASA ADS] [Google Scholar]
  25. De Donder, E., & Vanbeveren, D. 2003, New Astron., 8, 415 [NASA ADS] [CrossRef] [Google Scholar]
  26. de Mink, S. E., Langer, N., Izzard, R. G., Sana, H., & de Koter, A. 2013, ApJ, 764, 166 [NASA ADS] [CrossRef] [Google Scholar]
  27. Dewi, J. D. M., & Tauris, T. M. 2000, A&A, 360, 1043 [NASA ADS] [Google Scholar]
  28. Dewi, J. D. M., & Tauris, T. M. 2001, ASP Conf. Ser., 229, 255 [Google Scholar]
  29. Díaz-Rodríguez, M., Murphy, J. W., Rubin, D. A., et al. 2018, ApJ, 861, 92 [NASA ADS] [CrossRef] [Google Scholar]
  30. Dray, L. M., & Tout, C. A. 2007, MNRAS, 376, 61 [NASA ADS] [CrossRef] [Google Scholar]
  31. Drout, M. R., Soderberg, A. M., Gal-Yam, A., et al. 2011, ApJ, 741, 97 [NASA ADS] [CrossRef] [Google Scholar]
  32. Duchêne, G., & Kraus, A. 2013, ARA&A, 51, 269 [Google Scholar]
  33. Dunstall, P. R., Dufton, P. L., Sana, H., et al. 2015, A&A, 580, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Eggleton, P. P., & Tokovinin, A. A. 2008, MNRAS, 389, 869 [Google Scholar]
  35. Eldridge, J. J., & Stanway, E. R. 2009, MNRAS, 400, 1019 [NASA ADS] [CrossRef] [Google Scholar]
  36. Eldridge, J. J., & Tout, C. A. 2004, MNRAS, 353, 87 [NASA ADS] [CrossRef] [Google Scholar]
  37. Eldridge, J. J., Langer, N., & Tout, C. A. 2011, MNRAS, 414, 3501 [NASA ADS] [CrossRef] [Google Scholar]
  38. Eldridge, J. J., Fraser, M., Smartt, S. J., Maund, J. R., & Crockett, R. M. 2013, MNRAS, 436, 774 [NASA ADS] [CrossRef] [Google Scholar]
  39. Eldridge, J. J., Xiao, L., Stanway, E. R., Rodrigues, N., & Guo, N. Y. 2018, PASA, 35, 49 [NASA ADS] [CrossRef] [Google Scholar]
  40. Eldridge, J. J., Guo, N. Y., Rodrigues, N., Stanway, E. R., & Xiao, L. 2019, PASA, 36, e041 [CrossRef] [Google Scholar]
  41. Ertl, T., Janka, H.-T., Woosley, S. E., Sukhbold, T., & Ugliano, M. 2016, ApJ, 818, 124 [NASA ADS] [CrossRef] [Google Scholar]
  42. Farrell, E., Groh, J., Meynet, G., & Eldridge, J. 2020, MNRAS, 494, L53 [CrossRef] [Google Scholar]
  43. Filippenko, A. V. 1997, ARA&A, 35, 309 [NASA ADS] [CrossRef] [Google Scholar]
  44. Gal-Yam, A. 2017, in Observational and Physical Classification of Supernovae, eds. A. W. Alsabti, & P. Murdin, 195 [Google Scholar]
  45. Ge, H., Hjellming, M. S., Webbink, R. F., Chen, X., & Han, Z. 2010, ApJ, 717, 724 [NASA ADS] [CrossRef] [Google Scholar]
  46. Georgy, C., Ekström, S., Meynet, G., et al. 2012, A&A, 542, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Gerke, J. R., Kochanek, C. S., & Stanek, K. Z. 2015, MNRAS, 450, 3289 [NASA ADS] [CrossRef] [Google Scholar]
  48. Giannone, P. 1967, ZAp, 65, 226 [Google Scholar]
  49. Graur, O., Bianco, F. B., Modjaz, M., et al. 2017, ApJ, 837, 121 [NASA ADS] [CrossRef] [Google Scholar]
  50. Habergham, S. M., Anderson, J. P., James, P. A., & Lyman, J. D. 2014, MNRAS, 441, 2230 [NASA ADS] [CrossRef] [Google Scholar]
  51. Hamann, W.-R., Koesterke, L., & Wessolowski, U. 1995, A&A, 299, 151 [NASA ADS] [Google Scholar]
  52. Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003, ApJ, 591, 288 [NASA ADS] [CrossRef] [Google Scholar]
  53. Hellings, P. 1983, Ap&SS, 96, 37 [NASA ADS] [CrossRef] [Google Scholar]
  54. Hellings, P. 1984, Ap&SS, 104, 83 [NASA ADS] [CrossRef] [Google Scholar]
  55. Hobbs, G., Lorimer, D. R., Lyne, A. G., & Kramer, M. 2005, MNRAS, 360, 974 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hopkins, P. F., Kereš, D., Oñorbe, J., et al. 2014, MNRAS, 445, 581 [NASA ADS] [CrossRef] [Google Scholar]
  57. Humphreys, R. M., & Davidson, K. 1979, ApJ, 232, 409 [NASA ADS] [CrossRef] [Google Scholar]
  58. Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 [NASA ADS] [CrossRef] [Google Scholar]
  59. Hurley, J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 [NASA ADS] [CrossRef] [Google Scholar]
  60. Hut, P. 1980, A&A, 92, 167 [NASA ADS] [Google Scholar]
  61. Hut, P. 1981, A&A, 99, 126 [NASA ADS] [Google Scholar]
  62. Izzard, R. G., Tout, C. A., Karakas, A. I., & Pols, O. R. 2004, MNRAS, 350, 407 [NASA ADS] [CrossRef] [Google Scholar]
  63. Izzard, R. G., Dray, L. M., Karakas, A. I., Lugaro, M., & Tout, C. A. 2006, A&A, 460, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Izzard, R. G., Glebbeek, E., Stancliffe, R. J., & Pols, O. R. 2009, A&A, 508, 1359 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Janka, H.-T. 2012, Annu. Rev. Nucl. Part. Sci., 62, 407 [NASA ADS] [CrossRef] [Google Scholar]
  66. Jennings, Z. G., Williams, B. F., Murphy, J. W., et al. 2014, ApJ, 795, 170 [NASA ADS] [CrossRef] [Google Scholar]
  67. Jerkstrand, A., Fransson, C., Maguire, K., et al. 2012, A&A, 546, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Jerkstrand, A., Smartt, S. J., Fraser, M., et al. 2014, MNRAS, 439, 3694 [NASA ADS] [CrossRef] [Google Scholar]
  69. Justham, S., Podsiadlowski, P., & Vink, J. S. 2014, ApJ, 796, 121 [NASA ADS] [CrossRef] [Google Scholar]
  70. Kilpatrick, C. D., & Foley, R. J. 2018, MNRAS, 481, 2536 [NASA ADS] [CrossRef] [Google Scholar]
  71. Kiminki, D. C., & Kobulnicky, H. A. 2012, ApJ, 751, 4 [NASA ADS] [CrossRef] [Google Scholar]
  72. Kobulnicky, H. A., & Fryer, C. L. 2007, ApJ, 670, 747 [NASA ADS] [CrossRef] [Google Scholar]
  73. Kochanek, C. S., Beacom, J. F., Kistler, M. D., et al. 2008, ApJ, 684, 1336 [NASA ADS] [CrossRef] [Google Scholar]
  74. Kochanek, C. S., Khan, R., & Dai, X. 2012, ApJ, 759, 20 [NASA ADS] [CrossRef] [Google Scholar]
  75. Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
  76. Kuncarayakti, H., Anderson, J. P., Galbany, L., et al. 2018, A&A, 613, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Laplace, E., Götberg, Y., de Mink, S. E., Justham, S., & Farmer, R. 2020, A&A, 637, A6 [CrossRef] [EDP Sciences] [Google Scholar]
  78. Larson, R. B. 1974, MNRAS, 169, 229 [NASA ADS] [CrossRef] [Google Scholar]
  79. Li, W., Leaman, J., Chornock, R., et al. 2011, MNRAS, 412, 1441 [NASA ADS] [CrossRef] [Google Scholar]
  80. Lyman, J. D., Bersier, D., James, P. A., et al. 2016, MNRAS, 457, 328 [NASA ADS] [CrossRef] [Google Scholar]
  81. Maoz, D., & Badenes, C. 2010, MNRAS, 407, 1314 [NASA ADS] [CrossRef] [Google Scholar]
  82. Martinez, L., & Bersten, M. C. 2019, A&A, 629, A124 [CrossRef] [EDP Sciences] [Google Scholar]
  83. Maund, J. R. 2017, MNRAS, 469, 2202 [NASA ADS] [CrossRef] [Google Scholar]
  84. Maund, J. R., & Smartt, S. J. 2005, MNRAS, 360, 288 [NASA ADS] [CrossRef] [Google Scholar]
  85. Menon, A., & Heger, A. 2017, MNRAS, 469, 4649 [NASA ADS] [Google Scholar]
  86. Menon, A., Utrobin, V., & Heger, A. 2019, MNRAS, 482, 438 [NASA ADS] [CrossRef] [Google Scholar]
  87. Mestel, L. 1957, ApJ, 126, 550 [NASA ADS] [CrossRef] [Google Scholar]
  88. Mestel, L., & Moss, D. L. 1986, MNRAS, 221, 25 [NASA ADS] [CrossRef] [Google Scholar]
  89. Moe, M., & Di Stefano, R. 2017, ApJS, 230, 15 [Google Scholar]
  90. Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2007, A&A, 465, 1003 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Murphy, J. W., Jennings, Z. G., Williams, B., Dalcanton, J. J., & Dolphin, A. E. 2011, ApJ, 742, L4 [NASA ADS] [CrossRef] [Google Scholar]
  92. Nieuwenhuijzen, H., & de Jager, C. 1990, A&A, 231, 134 [NASA ADS] [Google Scholar]
  93. O’Connor, E., & Ott, C. D. 2011, ApJ, 730, 70 [NASA ADS] [CrossRef] [Google Scholar]
  94. Öpik, E. 1924, Tartu Obs. Publ., 25, 6 [Google Scholar]
  95. Pavlovskii, K., & Ivanova, N. 2015, MNRAS, 449, 4415 [NASA ADS] [CrossRef] [Google Scholar]
  96. Pavlovskii, K., Ivanova, N., Belczynski, K., & Van, K. X. 2017, MNRAS, 465, 2092 [NASA ADS] [CrossRef] [Google Scholar]
  97. Piro, A. L. 2013, ApJ, 768, L14 [NASA ADS] [CrossRef] [Google Scholar]
  98. Podsiadlowski, P. 1992, PASP, 104, 717 [NASA ADS] [CrossRef] [Google Scholar]
  99. Podsiadlowski, P., Joss, P. C., & Rappaport, S. 1990, A&A, 227, L9 [NASA ADS] [Google Scholar]
  100. Podsiadlowski, P., Joss, P. C., & Hsu, J. J. L. 1992, ApJ, 391, 246 [NASA ADS] [CrossRef] [Google Scholar]
  101. Pols, O. R., Schröder, K.-P., Hurley, J. R., Tout, C. A., & Eggleton, P. P. 1998, MNRAS, 298, 525 [NASA ADS] [CrossRef] [Google Scholar]
  102. Prentice, S. J., Ashall, C., James, P. A., et al. 2019, MNRAS, 485, 1559 [NASA ADS] [CrossRef] [Google Scholar]
  103. Renzo, M., Ott, C. D., Shore, S. N., & de Mink, S. E. 2017, A&A, 603, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Renzo, M., Zapartas, E., de Mink, S. E., et al. 2019, A&A, 624, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Ribas, I., Jordi, C., & Giménez, Á. 2000, MNRAS, 318, L55 [Google Scholar]
  106. Rui, L., Wang, X., Mo, J., et al. 2019, MNRAS, 485, 1990 [CrossRef] [Google Scholar]
  107. Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
  108. Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444 [Google Scholar]
  109. Schady, P., Eldridge, J. J., Anderson, J., et al. 2019, MNRAS, 490, 4515 [CrossRef] [Google Scholar]
  110. Schneider, F. R. N., Izzard, R. G., Langer, N., & de Mink, S. E. 2015, ApJ, 805, 20 [Google Scholar]
  111. Schneider, F. R. N., Ohlmann, S. T., Podsiadlowski, P., et al. 2019, Nature, 574, 211 [NASA ADS] [CrossRef] [Google Scholar]
  112. Shao, Y., & Li, X.-D. 2014, ApJ, 796, 37 [Google Scholar]
  113. Silverman, J. M., Pickett, S., Wheeler, J. C., et al. 2017, MNRAS, 467, 369 [NASA ADS] [Google Scholar]
  114. Smartt, S. J. 2015, PASA, 32, e016 [NASA ADS] [CrossRef] [Google Scholar]
  115. Smartt, S. J., Gilmore, G. F., Trentham, N., Tout, C. A., & Frayn, C. M. 2001, ApJ, 556, L29 [NASA ADS] [CrossRef] [Google Scholar]
  116. Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 395, 1409 [NASA ADS] [CrossRef] [Google Scholar]
  117. Smith, N. 2014, ARA&A, 52, 487 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  118. Smith, N., & Tombleson, R. 2015, MNRAS, 447, 598 [NASA ADS] [CrossRef] [Google Scholar]
  119. Smith, N., Li, W., Filippenko, A. V., & Chornock, R. 2011, MNRAS, 412, 1522 [NASA ADS] [CrossRef] [Google Scholar]
  120. Sravan, N., Marchant, P., Kalogera, V., & Margutti, R. 2018, ApJ, 852, L17 [NASA ADS] [CrossRef] [Google Scholar]
  121. Straniero, O., Dominguez, I., Piersanti, L., Giannotti, M., & Mirizzi, A. 2019, ApJ, 881, 158 [CrossRef] [Google Scholar]
  122. Sukhbold, T., Ertl, T., Woosley, S. E., Brown, J. M., & Janka, H.-T. 2016, ApJ, 821, 38 [Google Scholar]
  123. Taddia, F., Sollerman, J., Leloudas, G., et al. 2015, A&A, 574, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Tauris, T. M., & Dewi, J. D. M. 2001, A&A, 369, 170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  125. Tomasella, L., Cappellaro, E., Fraser, M., et al. 2013, MNRAS, 434, 1636 [NASA ADS] [CrossRef] [Google Scholar]
  126. Tout, C. A., Aarseth, S. J., Pols, O. R., & Eggleton, P. P. 1997, MNRAS, 291, 732 [NASA ADS] [Google Scholar]
  127. Tutukov, A. V., Yungelson, L. R., & Iben, I., Jr. 1992, ApJ, 386, 197 [NASA ADS] [CrossRef] [Google Scholar]
  128. Valenti, S., Howell, D. A., Stritzinger, M. D., et al. 2016, MNRAS, 459, 3939 [NASA ADS] [CrossRef] [Google Scholar]
  129. Van Dyk, S. D., Peng, C. Y., Barth, A. J., & Filippenko, A. V. 1999, AJ, 118, 2331 [NASA ADS] [CrossRef] [Google Scholar]
  130. Van Dyk, S. D., Li, W., & Filippenko, A. V. 2003, PASP, 115, 1 [NASA ADS] [CrossRef] [Google Scholar]
  131. Van Dyk, S. D., Zheng, W., Brink, T. G., et al. 2018, ApJ, 860, 90 [NASA ADS] [CrossRef] [Google Scholar]
  132. Van Dyk, S. D., Zheng, W., Maund, J. R., et al. 2019, ApJ, 875, 136 [CrossRef] [Google Scholar]
  133. Vanbeveren, D., Mennekens, N., Van Rensbergen, W., & De Loore, C. 2013, A&A, 552, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  134. Verbunt, F., & Cator, E. 2017, JApA, 38, 40 [Google Scholar]
  135. Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, A&A, 362, 295 [Google Scholar]
  136. Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  137. Walmswell, J. J., & Eldridge, J. J. 2012, MNRAS, 419, 2054 [NASA ADS] [CrossRef] [Google Scholar]
  138. Webbink, R. F. 1984, ApJ, 277, 355 [NASA ADS] [CrossRef] [Google Scholar]
  139. Williams, B. F., Peterson, S., Murphy, J., et al. 2014, ApJ, 791, 105 [NASA ADS] [CrossRef] [Google Scholar]
  140. Williams, B. F., Hillis, T. J., Murphy, J. W., et al. 2018, ApJ, 860, 39 [CrossRef] [Google Scholar]
  141. Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015 [NASA ADS] [CrossRef] [Google Scholar]
  142. Xiao, L., Stanway, E. R., & Eldridge, J. J. 2018, MNRAS, 477, 904 [NASA ADS] [CrossRef] [Google Scholar]
  143. Xiao, L., Galbany, L., Eldridge, J. J., & Stanway, E. R. 2019, MNRAS, 482, 384 [NASA ADS] [CrossRef] [Google Scholar]
  144. Yoon, S.-C., & Langer, N. 2005, A&A, 443, 643 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  145. Yoon, S.-C., Woosley, S. E., & Langer, N. 2010, ApJ, 725, 940 [NASA ADS] [CrossRef] [Google Scholar]
  146. Yoon, S.-C., Dessart, L., & Clocchiatti, A. 2017, ApJ, 840, 10 [NASA ADS] [CrossRef] [Google Scholar]
  147. Zahn, J.-P. 1977, A&A, 57, 383 [NASA ADS] [Google Scholar]
  148. Zapartas, E., de Mink, S. E., Van Dyk, S. D., et al. 2017a, ApJ, 842, 125 [NASA ADS] [CrossRef] [Google Scholar]
  149. Zapartas, E., de Mink, S. E., Izzard, R. G., et al. 2017b, A&A, 601, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  150. Zapartas, E., de Mink, S. E., Justham, S., et al. 2019, A&A, 631, A5 [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Additional table

Table A.1.

Model variations of the assumptions that affect the binary population specifically as well as our whole stellar population.

All Tables

Table 1.

Sample of Type II SN progenitors with inferred initial masses from both directly probing the SN progenitor and age-dating the surrounding environment.

Table A.1.

Model variations of the assumptions that affect the binary population specifically as well as our whole stellar population.

All Figures

thumbnail Fig. 1.

Initial mass to final core mass relation of SN II progenitors. We present the normalized 2D density distribution of the final core mass, MHe, preSN, with the initial mass, MZAMS, for all SN II progenitors in case they have an evolution in isolation (left panel) or in case they experience mass exchange with a binary companion (right panel). The color mapping is in logarithmic scale. We also depict the fitting formula for the one-to-one relation of single stars (Eq. (1); white dashed line) and for the median of the initial mass distribution for the case of binary progenitors (Eq. (2); red dashed line).

In the text
thumbnail Fig. 2.

Distribution of initial masses, MZAMS, of SN II progenitors experiencing binary interaction (red line) for increasing ranges of final core masses, MHe, preSN (on the right part). The contribution of each of the main binary scenarios is also shown in color below the distribution of binary-interacted progenitors. For reference we include the initial mass range of a single star corresponding to each MHe, preSN bin (black dashed line), which is shown in an arbitrary height. It becomes clear that the binary progenitors of SN II can originate from a wider range of lower on average MZAMS compared to single stars to form cores of similar mass, due to mass accretion or merging during their evolution. The results are for our standard assumptions although similar trends are found in all our model variations. The values on the y-axis show the number of SN II per bin in MZAMS per 106M of stars formed.

In the text
thumbnail Fig. 3.

Distribution of initial masses, MZAMS, of SN progenitors with MZAMS <  Mmin, WR ≈ 27 M. For a population of single stars (dashed black line) we expect the distribution to follow the IMF, with all progenitors between Mmin, ccSN ≈ 7.6 M and Mmin, WR leading to Type II SNe. The shape of the distribution changes for progenitors that experienced binary interaction (solid black line). The progenitors of SN II arising from a population of binary interacting systems (red region) originate from lower on average MZAMS, even below Mmin, ccSN. In addition, many stars with MZAMS <  Mmin, WR lead to stripped-envelope SNe due to binary interaction (blue region, stacked on top of the red region), avoiding a hydrogen-rich event. The results are for our standard assumptions and show the number of SNe per bin in MZAMS per solar mass of star formation.

In the text
thumbnail Fig. 4.

Distribution of the pre-SN helium core mass of SN II progenitors in case they experienced isolated evolution (black dashed line) or have followed a channel of binary interaction (red solid line). The contribution of each different binary scenario is also shown (dashed and/or dotted lines). The equivalent single-star initial masses corresponding to each MHe, preSN (inverse of Eq. (1)) is shown on the top x-axis. The results are for our standard assumptions and show the number of SN II per logarithmic bin in MHe, preSN per solar mass of star formation.

In the text
thumbnail Fig. 5.

Cumulative distribution of final core masses, MHe, preSN, of Type II SN progenitors. White squares are the inferred core mass or upper limits of observed SN II progenitors from Smartt (2015). For comparison we show our predicted cumulative distribution for progenitors that evolved as single stars (black dashed line) or including possible binary interactions (red solid line), following our standard assumptions. The equivalent single-star initial mass, MZAMS, corresponding to each MHe, preSN is shown on the top x-axis. The equivalent pre-SN progenitor luminosity, corresponding to each MHe, preSN value according to single STARS models which were also used in Smartt et al. (2009), is shown on the bottom x-axis. We also show the range of uncertainty (red-shaded region with white dashed lines) of the expected distributions found for our simulated populations in which we vary the treatment of binary physics. These include among others the efficiency of mass transfer and of common envelope evolution, the initial period and mass ratio distribution and the binary fraction.

In the text
thumbnail Fig. 6.

Same as Fig. 5 (with logarithmic x-axis here), but showing the results for other varying assumptions (blue lines): (a) the mass threshold for failed SN collapsing onto black holes Mmax, ccSN (Models 52, 53, 54 and 29), (b) the wind mass loss efficiency η (Models 25, 26 and 51), (c) the minimum mass for SN, Mmax, ccSN (Models 30 and 31) and (d) the slope of IMF αIMF (Models 32 to 34). These physical parameters influence the whole population, both single stars and in binary systems. The list of values depicted in the upper left corner shows a sequence of the varying distributions (blue lines) spanning from top left to bottom right of each panel.

In the text
thumbnail Fig. 7.

Discrepancy in the inferred initial masses between methods that probe the pre-SN core mass or luminosity of the progenitor vs. the ones that age-date its surrounding environment. Our results (2D probability density distribution in green) show the transformation of the simulated delay-time of the SN progenitor to MZAMS (x-axis), assuming a lifetime-mass relation for single stars (Fig. 1 of Z17). The y-axis is the transformation of the simulated MHe, preSN to MZAMS using the inverse of Eq. (1). The points are nearby Type II SNe with observed estimates of MZAMS with both groups of techniques, shown in Table 1.

In the text

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

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

Initial download of the metrics may take a while.