EDP Sciences
Home All issues Volume 502 / No 2 (August I 2009) A&A, 502 2 (2009) 611-622 Full HTML
Free Access
Issue
A&A
Volume 502, Number 2, August I 2009
Page(s) 611 - 622
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/200811339
Published online 15 June 2009

The different progenitors of type Ib, Ic SNe, and of GRB

C. Georgy1 - G. Meynet1 - R. Walder1,2 - D. Folini2 - A. Maeder1

1 - Geneva Observatory, Geneva University, 1290 Sauverny, Switzerland
2 - École Normale Supérieure, Lyon, CRAL, UMR CNRS 5574, Université de Lyon, France

Received 13 November 2008 / Accepted 15 May 2009

Abstract
Aims. We discuss the properties of the progenitors of core collapse supernovae of type Ib and Ic and of long soft gamma ray bursts, as they can be deduced from rotating stellar models of single stars at various metallicities.
Methods. The type of the supernova progenitor was determined from the surface abundances at the pre-supernova stage. The type of the supernova event was obtained from the masses of hydrogen and helium ejected at the time of the core-collapse supernova event.
Results. We find that the minimum amount of helium ejected by a core-collapse supernova (of whatever type) is around 0.3 ${M}_\odot $. There is no difference between the WC and WO stars in the ejected masses of helium, CNO elements, and heavy elements. Also no difference is expected between the chemical composition of a WC star resulting from a normal or a homogeneous evolution. The progenitors of type Ib supernovae are WNL, WNE, or less massive WC stars. Those of type Ic are WC and WO stars. WO stars are produced in a limited mass range (around 60 ${M}_\odot $) and only at low metallicity (for $Z\le\sim 0.010$) as already found. The WO stars are the progenitors of only a small fraction of type Ic. Present stellar models indicate that, at solar metallicity, there is about 1 type Ib supernova for 1 type Ic, and this ratio rises to 3 type Ic for 1 type Ib SN at twice solar metallicity. At this metallicity, type Ic's are more frequent than type Ib's because most massive stars that go through a WNE stage evolve further into a WC/WO phase. Current models can account for the observed number ratios SN Ib/SN II and SN Ic/SN II and for their observed variation with the metallicity. In case no supernova occurs when a black hole is formed, single-star models can still account for more than half of the observed (SN Ib+SN Ic)/SN II ratio for $Z \ge Z_\odot$. For the gamma ray burst rate, our models produce too large a number for such an event, even if we restrict the progenitor to the WO stars. This confirms that only a fraction of the WC / WO stars evolve toward gamma ray burst event, most likely those arising from stars that were initially very rapid rotators.

Key words: stars: Wolf-Rayet - stars: supernovae: general - gamma rays: bursts

1 Introduction

Type Ib supernovae are core-collapse supernovae whose spectrum shows no hydrogen lines. The spectra of type Ic show no hydrogen and helium lines (see e.g. Wheeler et al. 1987; Nomoto et al. 1994). The progenitors of these core collapse supernovae are thus believed to be stars stripped of their original H-rich envelope for type Ib's and also of their He-rich envelope for type Ic's. Progenitors are therefore Wolf-Rayet stars. The observed frequency at solar metallicity of type Ibc supernovae is about 20% the frequency of type II supernovae (see e.g. Cappellaro et al. 1999), which represents a significant fraction of all core collapse supernovae. In four cases, the typical spectrum of a type Ic supernova has been observed associated with a long soft gamma ray burst (GRB) event (Woosley & Bloom 2006), indicating a privileged link between type Ic's and the most powerful supernova explosions observed in the Universe. That type Ibc represents a significant proportion of core collapse supernovae and the link between type Ic's and the long soft GRBs justify the study of the evolution leading to such events.

In this context, a very interesting feature observed by Prantzos & Boissier (2003) is the increase with the metallicity in the number fraction of type Ibc to type II supernovae. Such a trend is reproduced well by rotating single-star models accounting for metallicity-dependent stellar winds (Meynet & Maeder 2003). On the other hand, a well known scenario explaining this trend comes from close binary-star evolution (see for example Vanbeveren et al. 2007; Eldridge et al. 2008; Podsiadlowski et al. 1992). In this scenario, the hydrogen-rich envelope is removed either through a Roche lobe overflow process or during a common envelope phase, producing a WR star.

In the present work, we take the opportunity of recent observations by Prieto et al. (2008) from which separated frequencies for type Ib and Ic can be deduced to get a step further. The questions we want to address are: how many type Ib and Ic supernovae normalized to the number of type II supernovae can we expect from single-star models at different metallicities, how the numbers derived from single-star models compare with the observed values, and how the results change no supernova event occurs if every time a black hole forms. Answers to those questions will provide some hints as to the importance of the single star scenario in explaining the evolution leading to the type Ibc supernovae events and will provide useful quantities for more comparisons with similar outputs obtained from close binary evolution scenarios.

Other questions will be addressed, for instance, the absence of H and He lines does not necessarily imply that the associated elements are completely absent from the ejecta. It only tells us that the physical conditions (the chemical composition being only one of them) are such that no lines of these elements are formed at the time of the supernova event. Actually, the quantities of hydrogen and helium ejected by these two types of supernovae are not known. In this paper, guided by stellar models for single stars, we estimate the minimum quantities of hydrogen and of helium, we may expect in type Ib and type Ic supernovae. We also indicate the relations of filiation between type Ib, Ic supernovae and the WNL, WNE, WC, and WO stars.

The second section of this paper briefly recalls the main physical ingredients of the stellar models used in this study. Based on these models, we discuss in Sect. 3 the nature of the supernova progenitors, and in Sect. 4 the chemical composition of the supernova ejecta. In Sect. 5, we discuss the link between the chemical composition of the ejecta and the nature of the supernova. The nature of the stellar remnant is discussed in Sect. 6. Theoretical rates for different supernova types are estimated in Sect. 7 and compared with observations. Section 8 discusses the implications for the long soft GRB progenitors and finally conclusions are presented in Sect. 9.

2 Stellar models

All the rotating stellar models we use here come from Meynet & Maeder (2005,2003) (hereafter Papers X and XI). A complete description of the physics used can be found in these two papers. We just recall here some important characteristics:

  • all the models considered in this work are computed with an initial equatorial velocity of $300 ~{\rm km}~{\rm s^{-1}}$. This initial velocity implies time-averaged equatorial velocities between 160 and 250 km s-1 depending on the initial mass and metallicity (see Tables 1 in Papers X and XI). These values are in the range of observed values for OB stars; for instance, Huang & Gies (2006) present projected rotational velocities for 496 OB stars belonging to 19 young galactic clusters with estimated ages between 6 and 73 Myr. Mean $v\sin i$ values of 139, 154, and 151 km s-1 were obtained for groups of O9.5-B1.5, B1.5-B5.0, and B5.0-B9.0 type stars. These authors derived the underlying probability distribution for the equatorial velocities v and obtained a peak at 200 km s-1. Dufton et al. (2006) obtain a peak of v at 250 km s-1 with a full width half maximum of approximately 180 km s-1 for the unevolved targets in the galactic clusters NGC 3293 and 4755;
  • mass loss rates are those of Vink et al. (2001,2000), which take wind bistability into account. The rates of de Jager et al. (1988) are used outside the domain of application of the rates of Vink et al. (2001,2000). During the WR phase, the mass loss rates are those of Nugis & Lamers (2000). These mass loss rates, which account for the clumping effects in the winds, are lower by a factor of 2-3 than the ones used in our previous non-rotating stellar grids (Meynet et al. 1994);
  • the dependence of mass loss on rotation (Maeder & Meynet 2000) was taken into account, as is the dependence on the metallicity. During the non-WR phases of the present models, we assumed that the mass loss rates depend on the initial metallicity as $\dot M(Z)=(Z/Z_\odot)^{1/2} \dot M(Z_\odot)$ (Vink et al. 2001; Kudritzki & Puls 2000).
Models were computed for four metallicities: Z = 0.004, Z = 0.008, Z = 0.020 (standard) and Z = 0.040. The initial masses considered are indicated in Table 1. All the models were computed up to the end of the core helium-burning phase. The very short durations of the advanced phase imply that the star at the end of the core He-burning phase has already reached its final mass. Moreover, the decoupling between the rapid evolution of the core governed by the neutrino emission and the slow evolution of the envelope does not allow the envelope to change a lot during the advanced phases. Thus we consider that the properties of the outer layers of our models obtained at the end of the core He-burning phase are very near the one that would have obtained had we pursued the computation until the presupernova stage[*].

In previous papers (see e.g. Maeder et al. 2008; Meynet & Maeder 2000; Maeder & Meynet 2001, Paper X; Paper XI), we compared the outputs of the rotating models to many observed features of massive stars and showed that the models accounting for the effects of rotation gave much better agreement for fitting observed features such as:

  • surface enrichments;
  • the blue to red supergiant ratio in the SMC;
  • the variation with metallicity in the number ratio of WR to O-type stars;
  • the observed number ratio of WN to WC stars in the LMC and SMC;
  • the variations with metallicity in the number ratio of type Ibc to type II supernovae;
  • the existence of Wolf-Rayet stars showing both at their surface H- and He-burning products.
In Papers X and XI we discussed the last four points above in details. Here, as indicated in the introduction, we focus on the nature of the progenitors for the type Ib and type Ic SNe, taken separately.

3 Nature of supernova progenitors

The nature of the supernova progenitors was determined by the surface composition of the star at the presupernova stage. The following classification scheme was adopted:

  • all the stars not ending their lifetime as a Wolf-Rayet star end their nuclear life as a blue or a red supergiant. They are hereafter called supergiant (SG);
  • the star is considered to be a WR star when $\log T_{\rm eff} > 4.0$ and the mass fraction of hydrogen at the surface $X_{{\rm S}}$ is inferior to 0.4 (same criterion as in Meynet & Maeder 2005);
  • Wolf-Rayet stars with hydrogen at their surface ( $X_{\rm S} > 10^{-5}$) are WNL stars;
  • Wolf-Rayet stars without hydrogen on their surface ( $X_{\rm S} < 10^{-5}$) and with surface carbon abundance inferior to nitrogen abundance are WNE stars;
  • Wolf-Rayet stars with carbon abundance on the surface superior to nitrogen abundance are WC or WO stars. To distinguish between WC and WO, we follow Smith & Maeder (1991): if the ratio of $\frac{{\rm C}+{\rm O}}{{\rm He}}$ (in number) is less than 1, we have a WC, otherwise a WO.
A star enters the WR phase as a WNL, then may evolve through the sequence WNL  $\rightarrow$ WNE  $\rightarrow$ WC  $\rightarrow$ WO. A star can end its evolution at any of these stages. We do not consider the WN/WC phase here, which can occur as a very short transition phase between the WNE and WC phases.

Table 1:   For all our models, various quantities are given here (see text for details).

Table 1 gives for the 4 metallicities considered here the initial mass of the star $M_{\rm ini}$, the mass at the end of the helium-burning phase $M_{\rm end~He} $, the mass of the He-core, $M_{\rm He}$, defined as the mass interior to which the mass fraction of helium (or of its products after fusion) is higher than 0.75, the mass of the carbon-oxygen core, $M_{\rm CO} $, defined as the mass interior to which the mass fraction of carbon plus oxygen is higher than 0.75, the mass of the remnant $M_{\rm rem} $ (given by Hirschi et al. 2005), the type of the progenitor, the SN type produced (see Sect. 5), and the remnant type (see Sect. 6).

In Fig. 1, we show how the initial mass ranges of various supernova progenitors vary as a function of the metallicity. To build this figure we had to adopt a method to determine the mass limits between the different zones. To illustrate, let us consider one specific case. At Z=0.02, from Table 1 we know that the 25 ${M}_\odot $ stellar model ends its lifetime as a WNL star and the 40 ${M}_\odot $ as a WC star, but we know neither at which mass the transition from WNL to WC star occurs nor whether some stars will end their life in the WNE stage. To find the transition between the WNL and WNE stage, we deduce from the chemical structure of our 25 ${M}_\odot $ at the end of its evolution, the mass that should be removed in order for the star to become a WNE star (in that case 0.35 ${M}_\odot $ or 1.4% of its total initial mass). From our 40 ${M}_\odot $ stellar model at the end of its evolution, we deduce the mass it has lost since its entry into the WNE phase (1.53 ${M}_\odot $ or 6.1% of its initial mass). Assigning, by convention, a minus sign to this last quantity (mass lost) and a positive sign to the former one (mass which has to be lost), plotting these two quantities in a plane ``mass-to be/or- lost'' versus ``initial mass'', the abscissa of the intersection of the line connecting these two points with the horizontal line at ordinate 0 gives the initial mass for which the endpoint of the evolution just corresponds to the stage where ``no mass'' has to be lost to enter into the WNE phase and where no mass has been lost since the entry into that phase. Thus this is the minimum initial mass to end its evolution as a WNE star. Using similar procedures for all the limits, we obtain the regions shown in Fig. 1. Note also that we used the results of Ekström et al. (2008) for complementing the figure with Pop III stellar models.

Looking at that figure, we can make the following remarks.

  • as expected, the lower mass limit for having a WR progenitor decreases when the metallicity increases. We see that the dependence on the metallicity of this limit is much stronger at low Z than at high Z;
  • the mass range of stars ending their life as a WN star is relatively narrow. Below this mass range, the stars do not succeed in entering the WR phase; above it, mass loss rates are efficient enough for allowing the star to evolve further into the WC or WO phase;
  • we note that at metallicities higher than about 0.02, the WN progenitors are more or less equally distributed among WNL and WNE stars, while for metallicities below 0.01 all WN progenitors are WNL stars. Actually present stellar models predicts such a narrow range of masses for WNE progenitors below Z=0.008 that we neglected it. This can be explained in the following way: first let us point out that WNE stars can be produced only when the star has evolved beyond the core H-burning phase. Indeed, regions without hydrogen are present in the star only after the main-sequence phase has ended. At high metallicity, because of strong stellar winds, the outer part of the H-burning core is uncovered at an early stage when the size of the He-burning core is still small and has not yet transformed a significant part of the former H-burning core into carbon and oxygen. This allows the intermediate region between the He-burning core and the H-burning shell to have a large extension in mass. Thus, the WNE phase has a significant duration (see Fig. 9 of Meynet & Maeder 2005), which allows a narrow range of initial masses for which the endpoint of the evolution occurs when the star is in the WNE phase. At low metallicity, the outer part of the H-burning core is uncovered at later evolutionary stages when the He-burning core has grown in mass, reducing the extent in mass of the region with chemical composition typical of the WNE phase. This reduces the duration of the WNE phase and also the initial mass range of stars ending their evolution in that stage;
  • above 30 ${M}_\odot $ and for metallicities higher than 0.008, all stars end their life as WC/WO stars. For metallicities below 0.008, the minimum initial mass for stars ending their lifetime as WC stars rapidly increases when the metallicity becomes lower;
  • the WO progenitors region appears as an ``island'' in the ``WC'' sea, located at low metallicity. As first explained by Smith & Maeder (1991), the low metallicity position of this island comes from the fact that, in order to form a WO star, the He-burning core must be uncovered at a very advanced stage of its evolution when a lot of helium has been converted into carbon and oxygen. This is more easily realized at low metallicity where the mass loss rates are lower, hence the core is revealed at a later stage. We also note that the WO progenitors do arise neither from the most massive nor from the less massive stars having their He-burning core uncovered. This reflects the delicate interplay between two counteracting effects. On one hand, in order for a star to become a WO star, strong enough stellar winds are necessary; on the other hand, as explained above, the stellar winds should not be too strong, otherwise no WO stars but instead WC stars are formed!

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg1.eps} \end{figure} Figure 1:

Ranges of masses of different types of supernova progenitors at different metallicities. For each area, the type of progenitor is indicated in the figure (see text). The shading on the right indicates the area where formation of a black hole (BH) is expected; elsewhere, the remnant is a neutron star (NS) (color figure available online).

Table 2:   Initial mass $M_{\rm ini}$, total released mass $M_{\rm rel}$, mass and fraction of ejected H ($M_{\rm H}$), He ( $M_{\rm He}$), C ($M_{\rm C}$), N ($M_{\rm N}$), O ($M_{\rm O}$) and of heavy elements ( $M_{\rm heavy}$), with all masses in ${M}_\odot $.

Table 3:   Mass limits between type II and type Ib SNe ( $M_{\rm lim}^{\rm II/Ib}$), and between type Ib SNe and Ic SNe ( $M_{\rm lim}^{\rm Ib/Ic}$) in ${M}_\odot $ (see text).

Table 4:   Mass ranges for various supernova progenitors, supernova types, and compact remnants at different metallicities (in ${M}_\odot $).

4 Chemical composition of the ejecta

The masses of hydrogen and helium ejected at the time of the explosive event are the key quantities for determining the type of the supernova event. These quantities are given for each of our stellar models in Table 2, which indicates for each metallicity and initial mass considered the total mass ejected at the time of the supernova explosion (Col. 2), the mass of hydrogen, of helium, of carbon, nitrogen and oxygen, and of heavy elements (Cols. 3 to 8). Heavy elements are the sum of the mass ejected under the form of elements heavier than hydrogen and helium. For carbon and oxygen, we used relations given by Maeder (1992) between the mass of the CO core and the ejected mass of these two elements at the time of the supernova event.

From the results shown in Table 2, we can deduce the following trends.

  • most of supergiant progenitors eject more than 70% of their ejecta in the form of hydrogen and helium. At any given metallicity, the proportion of the ejecta under the form of hydrogen and helium is greater in lower initial mass stars;
  • the 20 ${M}_\odot $ stellar model at Z=0.04 is a special case of supergiant progenitor in the sense that it only ejects small amounts of hydrogen and helium. This model is a blue supergiant ( $\log T_{\rm eff}=4.338$) at the end of its stellar lifetime and is actually on the verge of becoming a Wolf-Rayet star. It should lose only 0.09 ${M}_\odot $ to uncover layers with a hydrogen mass fraction inferior to 0.4. Probably a slightly more massive star would show an evolutionary path of the type O-red SG-blue SG-WR. Let us mention here that the single-star evolutionary scenario predicts a strong anti-correlation in the number of RSG and of WR stars in single-aged populations (Leitherer 1999). Close binary evolution may produce WR stars more easily from lower initial mass stars (Podsiadlowski et al. 1992) and thus explain the simultaneous presence of RSG and WR stars even in stellar populations with age spread inferior to a few Myr. It would be interesting to obtain from well-observed, single-aged stellar populations quantitative constraints on the occurrence of RSG and WR stars and to compare them with population synthesis models accounting for single and close binary evolution;
  • our computed WNL progenitors eject less than 1 ${M}_\odot $ of hydrogen, but more than 1 ${M}_\odot $ of helium (more generally every time some hydrogen is present in the ejecta, the mass of ejected helium is above $\sim $2 ${M}_\odot $). They eject greater quantities of carbon and oxygen than supergiant progenitors. About 3/4 of the mass of the ejecta is in the form of heavy elements;
  • WC and WO progenitors eject more than 90% of their ejecta in the form of heavy elements. In these cases, the mass of ejected helium is between 0.31 and 0.54 ${M}_\odot $;
  • we see that the minimum helium content in the ejecta of a core collapse supernova is about $0.3~{M}_\odot$. This means in particular that type Ic SNe should at least contain this minimum amount of helium in their ejecta. This agrees with the results of non rotating models presented in Eldridge & Tout (2004), and with binary models of Eldridge (private communication). There is to date no clear measurement of the mass of helium contained in the ejecta of type Ic SNe (although see Elmhamdi et al. 2006). It would be very interesting to find a correlation or an anticorrelation between the mass of ejected helium and the mass of another element produced in the same region of the star that would be more easily detectable. Typically, such correlations exist between the ejected mass of He and of N, as can be seen in Fig. 3. We see that the correlation depends on the metallicity so that to use such theoretical guideline for the He-content, some measurement of the metallicity is needed (such correlations would allow an idea of the metallicity in case the amounts of helium and of nitrogen in the ejecta can be measured!). To advance on the question of the He-content in type Ic supernova ejecta, we would need similar correlations but of course with another element, since in type Ic's the amount of nitrogen is zero. Probably 22Ne would be an interesting candidate, since it traces the He-rich regions of the star. The difficulty would then be to distinguish the 22Ne abundances from that of 20Ne in the spectrum analysis. In a future work, we shall study this question;
  • There is no difference in the chemical composition of the ejecta of the WC and the WO stars.

5 The supernova types

As already emphasized in the introduction, the conditions for producing type Ib and type Ic supernovae are based on the presence/absence of H, He lines in the supernova spectrum. Here we adopt a criterion based on the quantities of hydrogen and of helium ejected at the time of the supernova event to link a given chemical structure of the supernova progenitor to a given supernova type. Let us call $m^{\rm min}_{\rm H}$ the minimum quantity of hydrogen that should be present in the ejecta for the spectrum of the supernova to present H lines and $m^{\rm min}_{\rm He}$, the minimum quantity of helium that should be present in the ejecta for the spectrum of the supernova to present He lines. In that case a given initial mass star M presenting at the end of its evolution a quantity $m_{\rm H}(M)$ of hydrogen in its envelope, respectively $m_{\rm He}(M)$, of helium will produce

  • a type II supernova if $m_{\rm H}(M) > m^{\rm min}_{\rm H}$;
  • a type Ib supernova if $m_{\rm H}(M) < m^{\rm min}_{\rm H}$ and $m_{\rm He}(M) > m^{\rm min}_{\rm He}$;
  • a type Ic supernova if $m_{\rm H}(M) < m^{\rm min}_{\rm H}$ and $m_{\rm He}(M) < m^{\rm min}_{\rm He}$.
The quantities $m^{\rm min}_{\rm H}$ and $m^{\rm min}_{\rm He}$ are not known a priori. In the case of $m^{\rm min}_{\rm H}=0$, we obtain the limiting masses, $M_{\rm lim}^{\rm II/Ib}$, indicated in the second column of Table 4, where stars with initial masses below $M_{\rm lim}^{\rm II/Ib}$ explode as type II supernovae and stars with masses above $M_{\rm lim}^{\rm II/Ib}$ explode as type Ib or type Ic supernovae. We see that, when the metallicity becomes higher, this mass limit decreases as expected. Interestingly also, we see that adopting $m^{\rm min}_{\rm H} = 0.25$, or 0.5 ${M}_\odot $, pushes $M_{\rm lim}^{\rm II/Ib}$ towards lower values, thereby decreasing the mass range for type II supernovae and increasing it for type Ibc supernovae. Therefore $m^{\rm min}_{\rm H}=0$ leads to a lower limit for the number fraction of type Ibc to type II supernovae. Except at the metallicity Z=0.004, the changes remain quite modest. In the following we adopt the value $m^{\rm min}_{\rm H}=0$ as our reference value keeping in mind that it may underestimate the number fraction of type Ibc to type II supernovae.
 \begin{figure} \par\includegraphics[angle = 0, width=9cm]{11339Fg2.eps} \end{figure} Figure 2:

Same as Fig. 1 for various types of SNe. The SN type is given for $m^{\rm min}_{\rm H}=0$, $m^{\rm min}_{\rm He}=0.6$ and supposing that an SN occurs even when a BH is formed (see text).

The choice of $m^{\rm min}_{\rm He}$ is more delicate. We first note that there are no models without any helium in the ejecta. For all the metallicities and masses considered here, the minimum value of the mass of helium in the ejecta is 0.31 ${M}_\odot $. Second, we see that passing from the value 0.4 to 0.5 changes the mass limit a lot between the type Ib and the type Ic. Still, some changes are brought when one passes from 0.5 to 0.6, but in general are much more modest except at solar metallicity. Choosing values for $m^{\rm min}_{\rm He}$ between 0.6 and 1-1.5 ${M}_\odot $ would not change the results significantly when adopting 0.6. From these considerations it appears that 0.6 is a kind of limiting value, because below that value, the mass limit between type Ib and type Ic presents an high sensitivity to the exact value adopted for $m^{\rm min}_{\rm He}$ (for instance, choosing 0.4 would make type Ic appearing only in a narrow interval of metallicities). Above it the sensitivity of the mass is much weaker and the mass limits are much less dependent on the exact value chosen. In that respect the value of 0.6 does appear to us as the most reasonable choice. Also this choice gives a monotonic decrease in the mass limit between type Ib and type Ic when the metallicity increases and allows all WC and WO stars to be type Ic progenitors. Thus, in the following we adopt the limiting values corresponding to $m^{\rm min}_{\rm He}=0.6$.

We are now able to determine the type of the supernova it will give for each model (see Col. 7 of Table 1). Table 4 gives for each metallicity the initial mass ranges leading to different progenitors, SN events, and stellar remnants. A comparison between the mass ranges for the different types of progenitors and of SN types shows that all SG stars produce a type II SN. The WNL stars generally lead to type II SNe except may be those WNL stars on the verge of becoming WNE stars. They present nearly no hydrogen on their surface (see the case of 25 ${M}_\odot $ model at Z=0.02). WNE stars contain no hydrogen, but enough helium to become a type Ib SN. WC stars end their life as type Ic except the less massive ones that contain too much helium to be classified as a type Ic SN and thus would explode as a type Ib. Finally, all the WO stars lead to a type Ic SN event.

Using Table 4, we can build Fig. 2, showing the different supernova types expected for various initial masses and metallicities. As for Fig. 1, we use the models of Ekström et al. (2008) to complete the figure for Z = 0.0. We see that type Ib SNe arise from a mass band recovering the WNL / WNE band (see Fig. 1) and the lower range of WC stars. Interestingly, the estimated initial mass of the progenitor of SN 2008ak, which exhibits a transition spectrum between type II and type Ib, is in the range of $25{-}30~{M}_\odot$ for a metallicity around solar (Pastorello et al. 2008; Crockett et al. 2008). This agrees with the present stellar models. Type Ic SNe cover a much broader range of initial masses than type Ib. This trend reflects that most stars which, at a given stage, in their evolution have peeled off their H-rich envelope do not stop at that stage. Their evolution drives them beyond up to the stage where most of their He-rich layers have been peeled off.

One of the main uncertainties of our models is the mass loss rate. Some recent studies suggest that this rate could be overestimated by a factor of 2 (Vanbeveren et al. 2007). We briefly discuss here qualitatively how the results presented here would be modified if we had considered lower mass loss rates. The main differences are:

  • the entrance in the WR phase is delayed, so this phase is shorter;
  • at a given metallicity, the minimum initial mass to obtain a WR star is increased;
  • according to this last point, the minimum initial mass providing a type Ibc SN at a given metallicity is also increased.
Globally, the decrease in the mass loss rate in our models corresponds to a shift in the limits shown in in Figs. 1 and 2 up along the metallicity (vertical) axis.
 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg3.eps} \end{figure} Figure 3:

Correlation between the mass of nitrogen in the supernova ejecta and the mass of helium. The SN II region (red) and SN Ib region (blue) are represented, as well as the SN Ic region (on the left of the short dashed line). Masses are given in ${M}_\odot $ (color figure available online).

6 Nature of the stellar remnants

Superposed to Figs. 1 and 2, we have indicated, as a light grey zone, the regions where a black hole (BH) might be formed instead of a neutron star (NS) (see also Table 4). For drawing this zone, we assumed that $2.7\;{M}_\odot$ is the maximum mass of a NS. This value is compatible with the one given by Shapiro & Teukolsky (1983) and also with the recent discovery of a massive NS of $2.1~{M}_\odot$ (Freire et al. 2008). By adopting this mass limit and the relation $M_{\rm NS}$ versus $M_{\rm CO} $ deduced from the models of Hirschi et al. (2005), we can estimate the mass ranges of stars producing a NS for each metallicity, or a BH. For $Z\le 0.01$, all stars with an initial mass higher than $30~{M}_\odot$ finish their life as a BH, that is, all WR stars and the most massive red- or blue-supergiants. The minimum initial mass to become a BH is higher at higher metallicities. There is also a maximum initial mass to become a BH at high metallicity. This comes from the intense stellar winds that peel off the most massive stars and prevent them from having sufficient massive cores to produce a BH. As a result, the mass range of BH progenitors decreases when more and more metal-rich environments are considered, and for $Z \ge \sim 0.040$, no BH is expected from the single-star scenario (see Fig. 1).

We note that, except at low metallicity ($Z\le 0.01$), the majority of WNL and WNE stars produce a NS. WC stars may produce either NS or BH for metallicities higher than about 0.01. Below that metallicity they would only produce BH. WO stars produce only BH. In the same manner looking at Fig. 2, we see that most of type Ib's produce a NS. Type Ic's produce either a BH or a NS for metallicities above 0.02. For lower metallicities, only BH are expected to result from a type Ic event.

7 Frequency of type Ib, Ic SNe, and comparison with observations

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg4.eps} \end{figure} Figure 4:

For each metallicity, the left column represents the ratio type Ib SNe/type II SNe, and the right one the ratio type Ic/type II. The contribution of each kind of progenitor is also represented in each column. We used our reference case to built this figure (color figure available online).

In Fig. 4 we show the variation with metallicity in the number ratio of type Ib, Ic SNe for each type of progenitor with respect to type II SNe. To obtain these results, we used a Salpeter initial mass function (IMF), and apply the relation (in case of the Ib/II ratio):

\begin{displaymath}\frac{{\rm SN~I}_{\rm b}}{{\rm SN~II}} = \frac{\int_{M_{\rm m... ...{-2.35}{\rm d}M}{\int_8 ^{M_{\rm max,~ II}} M^{-2.35}{\rm d}M} \end{displaymath} (1)

where $M_{\rm min,~ b}$ is the minimum initial mass giving a type Ib SN, $M_{\rm max,~ b}$ the maximum initial mass giving a type Ib SN, and $M_{\rm max,~ II}$ the maximum initial mass giving a type II SN . We used the mass limits given in Table 4, and we supposed that an SN event occurs, even when a BH forms. Similar expressions are used to compute the other ratios with respect to type II SNe.

The first striking point when looking at Fig. 4 is that, for most of the metallicities studied here, type Ic are more frequent than type Ib. This was already apparent in Fig. 2, but here it is confirmed after the mass ranges involved have been weighted by the IMF. Only at solar metallicities would present models predict a similar frequency between Ib and Ic. A second point deserving to be mention is that the vast majority of the type Ibc supernovae come from WC or WO stars. WN stars contribute only to a relatively small proportion of the type Ib's and only at metallicities equal to or above solar.

These frequencies are compared with recent observed ratios in Fig. 5, and our models with the observed rates given by Prantzos & Boissier (2003); Smartt et al. (2009); Cappellaro et al. (1999) and Prieto et al. (2008)[*]. We also add the prediction of binary star models of Eldridge et al. (2008) and of Fryer et al. (2007). The gray areas indicate the variation of our computed ratio taking the mass limits $M_{\rm lim}^{\rm II/Ib}$ obtained varying $m^{\rm min}_{\rm H}$ between 0 and 0.5.

Our models reproduce the observations up to $Z = {Z}_\odot$[*]. At metallicities higher than about 0.02, present models predict too few type Ibc supernovae, although in view of the big error bars, it is difficult to ascertain that this is a real deficit. The increase in the ratio of type Ibc/type II SNe with the metallicity naturally arises from the metallicity dependence of the stellar winds. At high metallicity stellar winds are stronger, making the formation of stars without an H-rich envelope that explode as a types Ib or Ic supernova easier.

The results obtained including close binary evolutionary scenarios are also able to reproduce the general increase in the SN Ibc/SN II ratio with respect to the metallicity. In these models most of the type Ibc supernovae occur as a result of mass transfer in close binary systems. Thus we face here the situation where two very different models (single stars with rotation/close binary evolution with mass transfer) are both able to give a reasonable fit to the data. Actually both scenarios probably contribute to the observed populations of type Ibc supernovae. However, it would be interesting to know their relative importance and how their relative importance changes with the metallicity. It might be that both scenarios predict different behaviors for the way the frequencies of the type Ib and type Ic SNe vary as a function of the metallicity. Below we discuss the predictions of single-star rotating models and compare them with the available observations.

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg5.eps} \end{figure} Figure 5:

Rate of SN Ibc/SN II if all models produce an SN (solid line) or if models producing a BH do not explode in an SN (dashed line). Grey areas are the corresponding estimated errors from our models. Pentagons are observational data from Prieto et al. (2008), triangles are data from Prantzos & Boissier (2003), circle, and square, are measurements at solar metallicity from Cappellaro et al. (1999) and Smartt et al. (2009) respectively. The dotted line represents the binary models of Eldridge et al. (2008), and the dotted-dashed line the rate obtained with the binary models of Fryer et al. (2007).

From the data of Prieto et al. (2008), we extract separate observed rates for type Ib and type Ic SNe. They are represented for the SN Ib/SN II rate, and for SN Ic/SN II in Fig. 6. We note that, in the high metallicity range, the observed rate of type Ib SNe is lower than the rate of type Ic by about a factor two. The trend deduced from the present theoretical models is in good agreement with the observed one.

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg6.eps} \end{figure} Figure 6:

Rates of SN Ic/SN II (blue long-dashed line), SN Ib/SN II (red short-dashed line), and SN Ibc/SN II (black solid line) for the reference case. The points are extracted from the data of Prieto et al. (2008): triangles (red) represent the SN Ib/SN II rate, and upside-down triangles (blue) the SN Ic/SN II rate. Each triangle corresponds to a sample of 11 core collapse SNe (color figure available online).

If we adopt $m^{\rm min}_{\rm He} = 0.4$ instead of 0.6 (see Fig. 7), we see that our models never produce type Ic SNe (except a small number at $Z = 0.4~ {Z}_\odot$). In that case, models would predict that most type Ibc supernovae would be type Ib, which is contrary to what is observed. This numerical experiment illustrates the strong dependence of the present results on the maximum quantity of helium in the ejecta that can be accommodated in a type Ic supernova event. In case it were confirmed that this value is as low as $0.4~{M}_\odot$, then present models would not provide a good fit to the observed data. Enhancing this quantity to values superior or equal to 0.55-0.6 suffices to provide a good fit.

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg7.eps} \end{figure} Figure 7:

Same as Fig. 6 but with $m^{\rm min}_{\rm He}=0.4~ {M}_\odot$ instead of 0.6 (see text).

How do the results change when one makes the hypothesis that when a BH is formed no SN event occurs? The situation for the variation with metallicity in the number ratio of type Ibc to type II SNe is shown in Fig. 5. In Fig. 8, we show what would be the predictions of the present models for the frequencies of type Ib and type Ic taken separately and adopting $m^{\rm min}_{\rm He}=0.6$. At low metallicity, no type Ibc SNe are expected. At $Z \sim 0.010 = 0.5~{Z}_\odot$, type Ib SNe begin to come out of the BH domain, making the SN Ib/SN II ratio begins to increase. For $Z \ge 0.02$, single-star models might still account for more than half of the type Ibc supernovae. However, at Z=0.02, all the SNe coming from single-star evolution would have type Ib, which is not consistent with the observations. At twice solar metallicity, all our models lead to a neutron star, and the ratios are identical to the one obtained in Fig. 6. These comparisons show that, considering that the hypothesis ``when a BH is formed no SN event occurs'' makes a big difference for metallicities below about 0.02. Above that metallicity, the results are only weakly affected.

If we take a lower maximal mass for neutron stars, this situation becomes more extreme: almost all the WR stars produce a BH, thus the rate of SN Ibc/SN II is zero or very low at every metallicity. This illustrates the high sensitivity of the above results on either the mass limit for forming BH and/or the possibility or impossibility to have a SN event when a BH is formed.

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg8.eps} \end{figure} Figure 8:

Same as Fig. 6, with the assumption that when a black hole is formed, no SN event occurs (color figure available online).

Probably, the hypothesis according to which no supernova event appears when a BH is formed is too restrictive. Until no clear understanding of the physics involved in core collapse supernova explosion is reached, it is difficult to draw firm conclusions at that point, especially when one considers the possible effects of rotation. In that context, the collapsar scenario for GRBs (Woosley 1993) needs the formation of a BH (Dessart et al. 2008), and this formation is at least accompanied in some cases by a type Ic supernova event. Also, the observation of the binary system GRO J1655-40 containing a BH (Israelian et al. 1999) suggests that a few stellar masses have been ejected when the black hole formed and therefore that an SN event occurred. This is suggested by the important chemical anomalies observed at the surface of the companion whose origin is attributed the companion being accreted a part of the SN ejecta. These examples show that there are at least some circumstances where the formation of a BH does not prevent an SN event to occur. Probably, reality is somewhere between the two extreme cases discussed before: ``all massive stars produce an SN, independently of the remnant type'', and ``only the stars leading to a neutron star produce an SN''.

As in Sect. 5, it is interesting here to infer how a decrease in the mass loss rate modifies the previous results. At first order, a decrease in the mass loss rates is equivalent to a shift in the curves presented in Figs. 5 to 8 to the right along the metallicity (horizontal) axis. In general, in that case, the theoretical ratios will be lower than the observed values.

8 Gamma ray burst progenitors

It is now admitted that there is a connection between so-called long-soft GRB events and type Ic supernovae (Woosley & Bloom 2006). As seen in the previous section, this kind of supernova can be produced by both WC or WO stars. In this way, WC and WO stars are natural candidates to be GRB progenitors. However, GRB are primarily found in metal-poor environments (Modjaz et al. 2008), while type Ic's SNe appear mainly at high metallicity. Many physical reasons have been invoked to explain why GRBs seem to occur only at low metallicities. Among them are the following (see the review by Woosley & Bloom 2006):

  • at low metallicity, stellar winds (even during the WR phases) are weaker, thus bringing away small quantities of angular momentum. Black holes are also more easily formed since higher final masses are obtained;
  • at low metallicity the transport of angular momentum between the core and the envelope is less efficient than at high metallicities, because of slower meridional currents in metal poor stars;
  • since the chemical mixing due to rotation is more efficient at low Z, homogeneous evolution is more easily obtained in metal-poor regions. Homogeneous evolution allows massive stars to produce a type Ic SN event without having to lose large amounts of mass. Indeed, a perfectly homogeneous evolution (actually never realized) would allow the formation of a pure CO core (and then lead to a type Ic SN event) at the end of the core He-burning phase without the need for the star to lose any mass!
  • the distribution of initial velocities at low metallicity might contain more fast rotators than at high metallicities (see Fig. 9 in Martayan et al. 2007).
At the moment, these arguments remain quite speculative.

What can be done presently is to compare the observed GRB frequency with the observed frequency of potential candidates. First, as mentioned above, type Ic supernovae at low metallicity do appear interesting candidates. From Fig. 9, one can see that the observed rate of type Ic supernovae from single star models is still above the estimated number ratio GRB/core collapse supernovae (CCSNe) even when one only considers low-metallicity, type Ic SNe. Here we supposed that the formation of a BH does not prevent an SN event. In case of course only rare circumstances would allow an SN to occur when a black hole is formed, then the situation might be very different (see below).

Other interesting candidates are the WO stars that primarily occur at low metallicity. Table 5 lists all the known WO stars. If we take the value $\log({\rm O} / {\rm H}) + 12 = 8.66$ (Grevesse et al. 2007) for the solar oxygen abundance, we see that 6 out of a total of 8 occur in regions with a metallicity under $Z/{Z}_\odot < 0.9$. These stars will explode as a type Ic SN. The frequency of type Ic's SNe with WO star progenitors is shown in Fig. 9. The expected rate is only marginally compatible with the observed GRB rate (assuming that the aperture angle of the bipolar jet is very small, typically around $1^{\rm o}$). This conclusion has already been obtained by Hirschi et al. (2005), so even restraining the progenitors of GRB to WO stars would still not match the observed frequency of GRBs.

 \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg9.eps} \end{figure} Figure 9:

Observed rate of type Ic SNe (triangles with error bars) and predicted rate of type Ic SNe whose progenitor is a WO star (dashed line) with respect to the total number of core collapse SNe (CCSNe). The gray rectangle represents the extension in metallicity (Modjaz et al. 2008) and rate (Podsiadlowski et al. 2004) of GRB events.

Table 5:   WO stars catalog.

From the above discussion, it can be deduced that exploding as a type Ic SN in metal-poor regions (or having a WO progenitor) is not a sufficient condition for a GRB. Physical characteristics shared by a subsample of the metal poor type Ic events exists that are needed to obtain a GRB event. Probably this physical characteristic is the high angular momentum in the core (Yoon et al. 2006; Woosley & Heger 2006).

If GRB would only occur for initially very fast-rotating stars, rotating so fast that these stars would follow a homogeneous evolution (Meynet & Maeder 2007; Yoon et al. 2006; Woosley & Heger 2006; Maeder 1987), can we expect to find any peculiar feature in the chemical composition of the ejecta testifying this previous homogeneous evolution? Or in other words, is there any difference in the chemical composition of the ejecta between a type Ic having ``normal WC or WO'' progenitors and those arising from a model that followed a homogeneous evolution during the MS phase? To that purpose we computed the evolution of a $60~{M}_\odot$ model followed up to the end of central Si-burning, with an initial metallicity of Z=0.002 and an initial rotation rate of $\Omega / \Omega_{\rm crit} = 0.75$. The model was computed including magnetic field, following Spruit (2002) and Maeder & Meynet (2005). This model follows a nearly homogeneous evolutionary track during the MS phase. The angular momentum content at the time of the presupernova is sufficient for producing a collapsar. In Tables 1 and 2, we present various characteristics of that model.

Table 1 shows that the model is in the WC phase at the end of its evolution. The mixing induced by the magnetic field produces larger cores than in the models without magnetic field. Thus, due to mass loss, the core is uncovered at an early stage of He-burning, leading to a WC star. The helium content of the ejecta ( $0.64~ {M}_\odot$) is also slightly greater than in our other models of similar mass and is also just above the value of $0.6~ {M}_\odot$ below which a type Ic is supposed to occur in our classification scheme. This implies that our homogenous progenitors lead to a type Ib SN. On the other hand, as mentioned above, the limit of 0.6 has not to be taken as a firm limit, and stars with higher He-content might also explode in a type Ic event. Moreover, the models having followed a homogeneous evolution will end their evolution with more angular momentum in their central regions. This may have a significant impact on the way the stars explode, as well as on the spectral features arising from the SN explosion!

It is also interesting to compare the chemical composition of the ejecta as they are obtained in the present models when they follow a homogeneous evolution or not. Compare for instance the 60 ${M}_\odot $ at Z=0.004 (heterogeneous case) with the 60 ${M}_\odot $ at Z=0.002 (homogeneous case). The differences in the masses of He, CNO elements, and Z are very small. Thus there is little chance from the observations of the composition of the ejecta in CNO elements to be able to distinguish between a normal and a homogeneous evolution.

A point that would be interesting to check in a later work is the following: since the star rotate very fast during a great part of its evolution, it will produce anisotropic stellar winds (see Fig. 10), typically bipolar winds and probably equatorial mass loss when the critical limit is reached (Maeder 1999). These anisotropic winds will shape the circumstellar environment of the star (van Marle et al. 2008, and see also Fig. 10) and it might be that some traces of the resulting particular morphology will still be present at the time of the SN event (for instance, equatorial mass loss probably gives rise to a slow equatorial expanding disk whose traces might still be present when the star explodes). In that case the circumstellar environment of GRB may be peculiar. This can in turn have an impact on some features in the spectrum. More generally, the circumstellar environment of stars that have been in a not too remote past very fast rotators may also be characterized by such features. It would of course be very interesting to find in the circumstellar environment of some stars, now slowly rotating, traces of their very rapid rotation in a previous phase of their evolution.

 \begin{figure} \par\includegraphics[angle=0,width=8.5cm]{11339F10.eps} \end{figure} Figure 10:

Young bi-modal shaped nebula blown from the highly aspherical, polar dominated winds around a $20~{M}_\odot$ star at Z=10-5. The star has evolved during $21~552~{\rm yr}$ since it has left the ZAMS. The maximal extension of the wind is $\sim $ $1.5\times 10^{18}~{\rm cm}$. This model was computed by the A-Maze code (Walder & Folini 2000) (color figure available online).

9 Conclusion

In Meynet & Maeder (2005) we showed that single-star models can account for the increase with the metallicity in the ratio (SN Ib+SN Ic)/SNII. In the present work we obtained the following additional results:

  • Single-star models indicate that the mass of helium ejected in type Ic supernovae is at least 0.3 ${M}_\odot $ (see also Eldridge & Tout 2004).
  • Type Ib supernovae may be produced by WNL, WNE, or less massive WC stars. Type Ic SNe are only produced by WC or WO stars;
  • If we consider that the ejecta of a type Ic supernova may contain up to $0.6~ {M}_\odot$ of He and that an SN event occurs even when a BH is formed, single-star models are able to reproduce the observational cumulative ratio of SN Ibc/SN II well along with both detailed ratios SN Ib/SN II and SN Ic/SN II. Because of the small number of observed events with known $\log({\rm O}/{\rm H})$, the actual estimations of the ratios SN Ib/SN II and of SN Ic/SN II have large uncertainties. Additional measurements will permit better constrained stellar models;
  • In the framework of the ``BH-no supernova event'' hypothesis, all type Ibc supernovae at low metallicity should come from close binary evolution. At solar and higher-than-solar metallicities, more than half of the type Ibc might still be explained by single-star models;
  • Our models predict too high a rate of GRB, even considering that only WO stars lead to this rare event. Probably GRB only originate in the most rapid rotators;
  • There is no way to distinguish on the basis of the chemical composition of the ejecta a type Ic supernova arising from a normal heterogeneous evolution from a type Ic supernova arising from a homogeneous evolution. The past fast rotation of the GRB progenitor may have left some imprint on the circumstellar matter distribution;
  • For $Z \ge 0.04$, all massive stars produce neutron stars. When the metallicity decreases, the range of initial mass producing neutron stars decreases. Below a metallicity of about 0.01, all stars more massive than about $30~{M}_\odot$ produce a black hole.
It does appear difficult on the basis of the variation in the ratio (SN Ib + SN Ic)/SN II with the metallicity to differentiate the relative importance of the single and binary channels. As can be seen from Fig. 5, scenarios with both a significant proportion of close binaries and with a modest one are able to reproduce the observational trend. To make progress it would be interesting to study the predictions of the close binary channel for the relative frequency of type Ib and type Ic supernovae and to compare them with the observations. Also it would be useful to have predictions similar to those presented here for what concerns the chemical composition of the ejecta.

One can also note that, in case the hypothesis ``no SN occurs when a BH forms'' is correct, the binary scenario is responsible for most, if not all, type Ic and type Ib events at low metallicities, while at high metallicities single stars would still account for more than one half of the type Ibc supernovae. In the framework of that hypothesis, the long-soft GRB (at least those associated to type Ic events) would originate in binary systems. Let us note however that some observations show that at least in some circumstances, a SN occurs when a BH is formed.

A promising way to study the respective importance of the single and close binary channel for the formation of type Ibc SN progenitors is to look for single-aged stellar populations with masses at the turn off between 8 and $25~{M}_\odot$. If such clusters very rarely show both RSG and WR stars, then this can be interpreted as these two types of stars coming from very different mass ranges just as would be expected in single-star scenarios. If on the contrary, the simultaneous occurrence of both RSG and WR stars is frequent, then it may indicate that there is an overlap of the initial mass ranges of progenitors of RSG and WR stars. In that case, close binary evolution is probably required (except maybe at high metallicity where the overlap also exists for single-star scenarios).

Acknowledgements
The authors would like to thank the referee Philippe Podsiadlowski, as well as John J. Eldridge and Stephen J. Smartt for their precious comments that improved this work.

References

Footnotes

... stage[*]
We checked using the solar metallicity models of Hirschi et al. (2004), which were pursued until the presupernova stage, that indeed this is the case.
...Prieto et al. (2008)[*]
In their recent paper, Prieto et al. (2008) analyze a sample of 77 core collapse SNe (and 38 type Ia SNe) in the redshift range 0.01 < z < 0.04. This sample is composed of 58 SNe II, 13 SNe Ic, 3 SNe Ib, and 3 SNe Ibc, all with a measurement of $\log({\rm O}/{\rm H}) + 12$ for the host galaxy.
...$Z = {Z}_\odot$[*]
In Meynet & Maeder (2005), we already made this comparison, although with a less sophisticated method to link the progenitor to the supernova type it produces.

All Tables

Table 1:   For all our models, various quantities are given here (see text for details).

Table 2:   Initial mass $M_{\rm ini}$, total released mass $M_{\rm rel}$, mass and fraction of ejected H ($M_{\rm H}$), He ( $M_{\rm He}$), C ($M_{\rm C}$), N ($M_{\rm N}$), O ($M_{\rm O}$) and of heavy elements ( $M_{\rm heavy}$), with all masses in ${M}_\odot $.

Table 3:   Mass limits between type II and type Ib SNe ( $M_{\rm lim}^{\rm II/Ib}$), and between type Ib SNe and Ic SNe ( $M_{\rm lim}^{\rm Ib/Ic}$) in ${M}_\odot $ (see text).

Table 4:   Mass ranges for various supernova progenitors, supernova types, and compact remnants at different metallicities (in ${M}_\odot $).

Table 5:   WO stars catalog.

All Figures

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg1.eps} \end{figure} Figure 1:

Ranges of masses of different types of supernova progenitors at different metallicities. For each area, the type of progenitor is indicated in the figure (see text). The shading on the right indicates the area where formation of a black hole (BH) is expected; elsewhere, the remnant is a neutron star (NS) (color figure available online).
In the text

  \begin{figure} \par\includegraphics[angle = 0, width=9cm]{11339Fg2.eps} \end{figure} Figure 2:

Same as Fig. 1 for various types of SNe. The SN type is given for $m^{\rm min}_{\rm H}=0$, $m^{\rm min}_{\rm He}=0.6$ and supposing that an SN occurs even when a BH is formed (see text).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg3.eps} \end{figure} Figure 3:

Correlation between the mass of nitrogen in the supernova ejecta and the mass of helium. The SN II region (red) and SN Ib region (blue) are represented, as well as the SN Ic region (on the left of the short dashed line). Masses are given in ${M}_\odot $ (color figure available online).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg4.eps} \end{figure} Figure 4:

For each metallicity, the left column represents the ratio type Ib SNe/type II SNe, and the right one the ratio type Ic/type II. The contribution of each kind of progenitor is also represented in each column. We used our reference case to built this figure (color figure available online).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg5.eps} \end{figure} Figure 5:

Rate of SN Ibc/SN II if all models produce an SN (solid line) or if models producing a BH do not explode in an SN (dashed line). Grey areas are the corresponding estimated errors from our models. Pentagons are observational data from Prieto et al. (2008), triangles are data from Prantzos & Boissier (2003), circle, and square, are measurements at solar metallicity from Cappellaro et al. (1999) and Smartt et al. (2009) respectively. The dotted line represents the binary models of Eldridge et al. (2008), and the dotted-dashed line the rate obtained with the binary models of Fryer et al. (2007).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg6.eps} \end{figure} Figure 6:

Rates of SN Ic/SN II (blue long-dashed line), SN Ib/SN II (red short-dashed line), and SN Ibc/SN II (black solid line) for the reference case. The points are extracted from the data of Prieto et al. (2008): triangles (red) represent the SN Ib/SN II rate, and upside-down triangles (blue) the SN Ic/SN II rate. Each triangle corresponds to a sample of 11 core collapse SNe (color figure available online).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg7.eps} \end{figure} Figure 7:

Same as Fig. 6 but with $m^{\rm min}_{\rm He}=0.4~ {M}_\odot$ instead of 0.6 (see text).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg8.eps} \end{figure} Figure 8:

Same as Fig. 6, with the assumption that when a black hole is formed, no SN event occurs (color figure available online).
In the text

  \begin{figure} \par\includegraphics[angle=0,width=9cm]{11339Fg9.eps} \end{figure} Figure 9:

Observed rate of type Ic SNe (triangles with error bars) and predicted rate of type Ic SNe whose progenitor is a WO star (dashed line) with respect to the total number of core collapse SNe (CCSNe). The gray rectangle represents the extension in metallicity (Modjaz et al. 2008) and rate (Podsiadlowski et al. 2004) of GRB events.
In the text

  \begin{figure} \par\includegraphics[angle=0,width=8.5cm]{11339F10.eps} \end{figure} Figure 10:

Young bi-modal shaped nebula blown from the highly aspherical, polar dominated winds around a $20~{M}_\odot$ star at Z=10-5. The star has evolved during $21~552~{\rm yr}$ since it has left the ZAMS. The maximal extension of the wind is $\sim $ $1.5\times 10^{18}~{\rm cm}$. This model was computed by the A-Maze code (Walder & Folini 2000) (color figure available online).
In the text

Copyright ESO 2009

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.