1. Introduction
Gamma-ray bursts (GRBs) have been known for nearly 50 yr (Klebesadel, Strong, & Olson Reference Piran1973). Their properties and origins have been extensively studied (see Piran Reference Piran2004; Woosley & Bloom Reference Woosley and Bloom2006; Kumar & Zhang Reference Kumar and Zhang2015, for reviews). Since the Burst and Transient Source Experiment (BATSE) era (1990–2000), GRBs have been divided into two groups according to their duration and the spectral hardness ratio (Kouveliotou et al. Reference Kouveliotou1993). The short hard GRBs are believed to result from a black hole–neutron star (BH–NS) merger or NS–NS merger, while the long soft GRBs result from the collapsing of massive stars (Piran Reference Piran2004). It is clear that at least part of long GRBs are collapsars as they are connected with type Ic supernovae (SNe), as was discovered in the connection of GRB 980425 and SN 1998bw (Galama et al. Reference Galama1998), and then spectrally confirmed in the connection of GRB 030329 and SN 2003dh (Stanek et al. Reference Stanek2003). Hereafter, we use GRBs to denote long GRBs.
Type Ic supernova (SN Ic) is one of the final fates of a massive star (Filippenko Reference Filippenko1997; Heger et al. Reference Heger, Fryer, Woosley, Langer and Hartmann2003). Since the discovery of the connection between the long GRB and the SN Ic, the SN Ic can be divided into several sub-classes, which are normal SN Ic, GRB-associated SN Ic, X-ray flash (XRF)-associated SN Ic,Footnote a broad-line SN Ic (SN Ic-BL), and super-luminous supernova Ic (SLSN Ic).
A unified picture is always a fundamental goal for the researchers, as it reveals a simple essence for the diverse phenomena, such as the unified models for different types of active galactic nuclei (Antonucci Reference Antonucci1993). de Rújula (Reference de Rújula2005) even suggested a unified model for almost all the high- energy phenomena by their ‘cannon ball’ model (Dar & de Rújula Reference Dar and de Rújula2004). Yamazaki, Ioka, & Nakamura (Reference Yamazaki, Ioka and Nakamura2004) provided a unified picture for different types of GRBs by counting the numbers of sub-shells of the jets: short GRBs with a few sub-shells, long GRBs with a lot of sub-shells, and XRFs (and X-ray-rich GRBs) being observed off-axis. However, one should notice that the essential difference between the long GRBs and short GRBs is the progenitor. For the different types of GRBs based on their spectra, a unified picture has been suggested for XRFs, X-ray-rich GRBs, and GRBs by observing at different viewing angle of the same top-hat jet (Lamb, Donaghy, & Graziani Reference Lamb, Donaghy and Graziani2005), or intrinsically non-uniform jet with either a power law distribution (Rossi, Lazzati, & Rees Reference Rossi, Lazzati and Rees2002) or a Gaussian distribution (Zhang et al. Reference Zhang, Dai, Lloyd-Ronning and Mészáros2004a).
Different studies have been made to reveal the physics of the different types of SNe Ic, as well as their connection to the GRBs. Just after the first detection of the association of GRB and SN (GRB 980425/SN 1998bw (Galama et al. Reference Galama1998)), Cen (Reference Cen1998) suggested a unified scenario in which the explosion of SN also launches a pair of jets with opposite directions, which produces the GRB, and the asymmetry of the explosion provides a kick for the high-speed pulsar.
Through the numerical simulation of a relativistic jet propagating through the stellar envelope, Zhang, Woosley, & MacFadyen Zhang et al. (Reference Zhang, Woosley and MacFadyen2003) suggested that the XRFs and the GRBs may be a consequence of different viewing angles, and the GRB-associated SN may be powered by the non-relativistic part of the jets. Lazzati et al. (Reference Lazzati, Morsony, Blackwell and Begelman2012) performed a hydrodynamic simulation by keeping the total energy budget as constant. They found that a longer-lasting jet can overcome the envelope, while the shorter one may just power the SN component to be more energetic, and if it is short enough, the SN component cannot even be distinguished from a normal SN. Most recently, López-Cámara, Lazzati, & Morsony (Reference López-Cámara, Lazzati and Morsony2016) showed the GRB jets propagation through the progenitor envelope, depends largely on the parameters choosen.
Lei, Wang, & Ma (Reference Lei, Wang and Ma2005) suggested that the GRB jets are launched by the Blandford–Znajek process (Blandford & Znajek Reference Blandford and Znajek1977), while the SN component is launched by the spin energy transferred from the central BH to the accreting disk in a magnetic coupling process (van Putten Reference van Putten1999; van Putten, Della Valle, & Levinson Reference van Putten, Della Valle and Levinson2011).
Nomoto et al. (Reference Nomoto2006) tried to unify the GRBs and SNe by the mass of 56Ni, which depends on the rate of the energy deposited to form the 56Ni. They outlined the consequence of the deposition rate dropping as GRB connected with a hypernova, a less luminous SN, and without an SN connection. When the deposition rate is smaller than a critical value 3 × 1051 erg s−1, the process of nucleosynthesis changes dramatically and no considerable 56Ni is produced. After the discovery of GRBs 060505 and 060614, which are GRBs with no SN component observations at a deep flux limit (Della Valle et al. Reference Della Valle2006; Fynbo et al. Reference Fynbo2006), in the same scenario, with numerical simulation, Tominaga et al. (Reference Tominaga2007) showed that this kind of GRBs may have very little 56Ni been synthesised.
Kelly, Kirshner, & Pahre (Reference Kelly, Kirshner and Pahre2008) found long GRBs and SNe Ic have similar locations in host galaxies. Both GRBs and type Ic SNe are located in similar star-forming region.
After the discovery of SLSNe (e.g. Gal-Yam Reference Gal-Yam2012), people found the connection between them and GRBs. 111209A was observed connected with a possible SLSN SN 2011kl (Greiner et al. Reference Greiner2015; Kann et al. Reference Kann2016). Nakauchi et al. (Reference Nakauchi, Kashiyama, Suwa and Nakamura2013) suggested blue supergiant model for the connection of super-long GRBs and SLSNe. van den Heuvel & Portegies Zwart (2013) proposed that SLSNe and long GRBs are both from young dense star clusters of their host galaxies. Yu & Li (Reference Yu and Li2017) found the flares of them share similar empirical correlation between the luminosity and timescale. Margalit et al. (Reference Margalit, Metzger, Thompson, Nicholl and Sukhbold2018) suggested a misaligned magnetar model for magnetar thermalisation and jet formation to connect SLSNe and GRBs. Contini (Reference Contini2018) considered the spectral similarity among the host galaxies of the GRBs, SLSNe, star bursts, and active galactic nuclei.
Fan et al. (Reference Fan, Zhang, Xu, Liang and Zhang2011) suggested an alternated unified picture based on the central engine that: an energetic SN associated with a normal GRB (e.g. GRB 030329/SN 2003dh (Stanek et al. Reference Stanek2003)) comes from a powerful central engine plus immediate jets launching, an energetic SN associated with a low-luminosity GRB (e.g. GRB 980425/SN 1998bw (Galama et al. Reference Galama1998)) comes from a powerful central engine but with delayed jet launching from the envelope, and a less-energetic SN associated with a low-luminosity GRB (e.g. XRF 060218/SN 2006aj (Campana et al. Reference Campana2006)) originates from an essentially less powerful central engine (slowly rotating magnetar). On the other hand, Nakar (Reference Nakar2015) tried to unify the low-luminosity GRBs and normal long GRBs, by suggesting that they are both originated from a massive core collapsar, while the progenitor of the low-luminosity GRBs has an extra extended envelope. This envelope prevents the jets penetrating and it results in a low-luminosity GRB. They used this scenario to explain the association of the low-luminosity GRB 060218/SN 2006aj.
Recently, Wang et al. (Reference Wang, Wang, Dai and Wu2015) suggested that the energy from the 56Ni and the rotational energy of the newly born NS may power the luminous SNe Ic, with different fractions for different SNe Ic. Gilkis, Soker, & Papish (Reference Gilkis, Soker and Papish2016) showed that the jet feedback provides a variety of energy from the core-collapsing star, which can power both the super-energetic SN and the GRB jets. Piran et al. (Reference Piran, Nakar, Mazzali and Pian2017) showed the observational evidences for a possible connection between a normal core-collapsing SN and GRB by the choked jets.
In the diagram of kinetic energy versus velocity of the ejecta, Margutti et al. (Reference Margutti2014) showed a continuous sequence between normal SNe Ic and GRBs/SNe, and the total energy for both SNe (Margutti et al. Reference Margutti2014) and GRBs (Bloom, Frail, & Kulkarni Reference Bloom, Frail and Kulkarni2003) covers a similar relatively narrow range (∼1051−1053 erg).Footnote b This indicates some inner connections in GRBs/SNe, as also suggested by Gilkis et al. (Reference Gilkis, Soker and Papish2016), Piran et al. (Reference Piran, Nakar, Mazzali and Pian2017), and Yu et al. (Reference Yu, Zhu, Li, Lü and Zou2017). Motivated by the efforts on the unified scenario and the total energy of GRBs and SNe Ic are in the same order of magnitude, we suggest an inner connection for these different classes of SNe Ic, based on the behaviour of the core collapsing launched the pair of jets. A schematic diagram is shown in Figure 1. If the jets with energy ∼1051−1053 erg are not launched by the central engine, it appears as a normal SN Ic with no broad emission lines. The remaining classes occur when the launched jets plough into the envelope of the progenitor star. If the jets successfully passes through the envelope, and one points to the observer, it appears as a GRB, such as GRB 980425, while if it does not point to the observer, we can only see an SN Ic-BL, and an off-axis afterglow is expected, such as the possible candidate SN 2002ap. If the jets cannot successfully pass through the envelope, but most of the energy is transferred into the relativistic shock breakout, we can see it as an XRF, such as GRB 060218, while if most of the energy is deeply dissipated into the envelope, a major part of the energy is transferred into optical photons and appears as an SLSN.Footnote c We describe the scenario and its evidences in Section 2 and discuss it in Section 3.
A diagram for the different core-collapsing SNe Ic without considering the effect of the line of sight. If the jet successfully passes through the envelope, a GRB is associated with an SN Ic-BL. If the jet is blocked and the energy is released mainly as a relativistic shock breakout, an XRF is associated with the SN Ic-BL. If the jet is totally blocked, and most of the energy goes to the SN, it appears as an SLSN. While if there is no jet launched, it appears as a normal SN Ic.
2. Scenario and evidences
As suggested by Bromberg et al. (Reference Bromberg, Nakar, Piran and Sari2012), the collapsar-originated GRBs may also represent as short bursts, as the GRB jets may break out of the envelope at different times.Footnote d It is also possible that the jets may not even overcome the envelope. However, the energy must come out, though it cannot come out as γ-rays. A very promising possibility is that the energy of the GRB jets is converted into optical photons and results in an SLSN. Numerical simulations also show that the relativistic jet may or may not successfully break the envelope (Zhang et al. Reference Zhang, Woosley and MacFadyen2003; Woosley & Bloom Reference Woosley and Bloom2006). Based on the fact that the core collapsar may or may not launch a pair of jets, and the jets may or may not penetrate the envelope, we propose the unified scenario as shown in Figure 2.
The full scenario is the following. The progenitor of SN Ic may or may not have a pair of jets. In the case of no jets, it becomes a normal SN Ic (e.g. SN 1994I). If it launches a pair of jets, there are two situations depending on whether the jets overcome the envelope of the progenitor or not. If the jets cannot pass through the envelope, the total energy of the jets is injected into the envelope. It becomes an SLSN (e.g. SN 2007bi). On the other hand, if some energy escapes as a shock breakout, one can see the XRF-connected SN (e.g. SN 2006aj—GRB 060218; SN 2008D—XRF 080109). If the jet passes through the envelope successfully and is pointing to the Earth, we see a GRB (e.g. SN 1998bw—GRB 980425). On the other hand, if the jet is not pointing to the Earth, we only see an SN Ic-BL (e.g. SN 2002ap). Energy injection of the jets can create relative motion of particles in the envelope to produce broad-line feature in the SN spectrum. This is a common result for SN Ic–XRF, normal SN Ic-BL, and SN Ic–GRB.
In the following subsections, we show evidences or indications for supporting this scenario.
2.1. Luminosity distributions
One piece of evidence for this scenario is the continuity among different types of SNe Ic in the luminosity space, which can be seen in the luminosity-distance diagram, as shown in Figure 3. The peak V-band absolute magnitude versus the luminosity distance is plotted for normal SNe Ic, SNe Ic-BL, SNe Ic/XRF, SNe Ic/GRB, and SLSNe. On the bottom right of the panel, one can see a big empty region below the data of SNe Ic/GRB. This is a selection effect. For the normal SNe, they are not luminous enough to trigger the survey at far distances, while for the GRB-connected SNe, the position is determined by the GRB satellites and follow-up afterglow monitoring, and then more powerful telescopes with longer observation time are used, so the much dimmer SN component can be identified. Therefore, in this blank region, there should be more SNe Ic hidden. On the upper/lower left of the panel, the blank regions do mean there are no objects.
Because of the observational selection effect, one cannot derive the individual properties of the SNe population by directly counting the observed events. An alternative estimation could be performed based on the edge of the blank area on the upper left and bottom left of the panel, which shows roughly 
${L_{{\rm{p}},V}} \propto D_{{\rm{L}},e}^{0.8}$, for 
${M_{{\rm{p}},V}} \lt - 18.6$ (
${L_{{\rm{p}},V}} \gt 8.7 \times {10^{42}}{\kern 1pt} {\rm{erg}}{\kern 1pt} {{\rm{s}}^{ - 1}}$), and 
${L_{{\rm{p}},V}} \propto D_{{\rm{L}},e}^{ - 1.2}$, for 
${M_{{\rm{p}},V}} \ge - 18.6$ (
${L_{{\rm{p}},V}} \lt 8.7 \times {10^{42}}{\kern 1pt} {\rm{erg}}{\kern 1pt} {{\rm{s}}^{ - 1}}$), where L p,V is the peak V-band luminosity directly from the absolute magnitude, and D L,e is the luminosity distance for the edge. Assuming a uniform space distribution of the SNe, the D L,e means the average distance between two SNe with the same luminosity, that is, 
${D_{{\rm{L}},e}} \sim {[{1 \over {n(L)}}]^{1/3}}$, where n(L) is the number density at a certain luminosity. Therefore, we can get a rough luminosity function for all the SNe Ic (Zou Reference Zou2017):
\begin{equation} n(L) \propto\left\{ {\matrix{{L_{{\rm{p}},V}^{2.5},} \hfill& {{L_{{\rm{p}},V}} \lt 8.7 \times {{10}^{42}}{\kern 1pt} {\rm{e}}rg{\kern 1pt} {s^{ - 1}},} \hfill\\{L_{{\rm{p}},V}^{ - 3.75},} \hfill& {{L_{{\rm{p}},V}} \ge8.7 \times {{10}^{42}}{\kern 1pt} {\rm{erg}}{\kern 1pt} {{\rm{s}}^{ - 1}},} \hfill\cr} } \right. \end{equation}
which is roughly similar to the luminosity function of SNe Ia (Yasuda & Fukugita Reference Yasuda and Fukugita2010) while different from that of GRBs (Yu et al. Reference Yu, Wang, Dai and Cheng2015). In comparing with luminosity-distance diagram for all kinds of SNe shown in Figure 2 of Richardson, Branch, & Baron (Reference Richardson, Branch and Baron2006), one can see no edge on either the upper- or lower-left area. It is simply because different types of SNe obey different luminosity functions. The edge as clearly shown in our Figure 3 suggests that different types of SNe Ic are under the same population in some sense.
Schematic diagram for the scenario of different SNe Ic. Depending on whether there is a pair of jets, whether the jet successfully passes through the envelope, whether there is a strong shock breakout, and whether the jet is pointing to us, the observed SN Ic can be a normal SN Ic, an SLSN Ic, an SN Ic connected to an XRF, a broad-line SN Ic, or an SN Ic connected to a GRB, respectively.
Luminosity-distance diagram (Richardson et al. 2014) (V-band peak absolute magnitude of the SN component vs. the luminosity distance in unit of Mpc) for all the SNe Ic, including normal SNe Ic (yellow triangles), SNe Ic-BL (green rectangles), SNe Ic associated with XRFs (marked as XRF, pink rectangles with crosses), SNe Ic associated with GRBs (marked as GRB, red dots), and SLSN Ic (SLSN-Ic, blue crosses). The two black solid lines show the edge of the data points, with slopes ∼2 and ∼−3, respectively. The intersection locates at M p,V ∼−18.6. Data are mainly from the Asiago Supernova Catalogue (online updating data at http://cdsarc.u-strasbg.fr/viz-bin/Cat?B/sn) Barbon et al. Reference Barbon, Buondí, Cappellaro and Turatto1999, and a few other individuals are taken from literature which can be found in the references of the text.
Concentrating on the peak magnitude in Figure 3, one can see, roughly, that the luminosity is increasing in the sequence: normal SNe Ic, SNe Ic-BL,Footnote e and SLSN. They are overlapping, while a Kolmogorov–Smirnov (K–S) test also shows that the SNe Ic-BL is intrinsically brighter than the normal SNe Ic (
${M_{{\rm{p,V,Ic}}}} = - 18.0 \pm 0.5$, 
${M_{{\rm{p,V,Ic - BL}}}} = - 18.3 \pm 0.8$ for V-band, and 
${M_{{\rm{p,R,Ic}}}} = - 18.3 \pm 0.6$, 
${M_{{\rm{p,V,Ic - BL}}}} = - 19.0 \pm 1.1$ for R-band) (Drout et al. Reference Drout2011). For the normal SNe Ic, there is no jet and consequently no extra energy. Therefore, the luminosity is the lowest. For the SNe Ic-BL, they are either associated with on-axis GRB/XRF or off-axis GRB/XRF (no GRB/XRF observed). Some extra energy from the jets goes into the progenitor envelope. Therefore, the luminosity of the SN component is brighter. For the SLSNe, most energy of the jets is deposited into the envelope. Therefore, the luminosity is the highest. Notice the intrinsic variety of the SNe and the jets; the difference of the SNe Ic is not strictly departed.
From the point of view of energy, SN-XRF, SN-GRB, and SLSN all have similar amounts of total energy. Normal SN Ic and SN Ic-BL should have smaller total energy because of no jet component or jet component is pointing to other direction. As the peak luminosity of the SLSN reaches more than 1044 erg/s, and peaks at around 100 d, the total radiation energy reaches to ∼1051 erg (Gal-Yam Reference Gal-Yam2012). The energy of the GRBs is also of the order of 1051 erg.Footnote f The inferred kinetic energy of a normal SN Ic is ∼1051 erg, and it is ∼1052 erg for an SN Ic-BL (Drout et al. Reference Drout2011). These figures are consistent with the inner connection, that roughly speaking, if the jet energy is less than ∼1051 erg, there is no observational effect and it appears as a normal SN Ic, while if the jet energy is larger than this value, it appears as an SLSN or a GRB, etc. Notice again, they should also have some diversity because of the diversity of the intrinsic energy reservoirs, which could be any of the GRB central engines, such as a BH accretion disk system or a magnetar.
2.2. SN Ic-GRB
As the GRB comes from the relativistic jet, for the GRB-connected SN Ic, it has to produce a pair of jets. One directed to the observer, and the other in the opposite direction. The spectra of many events have been confirmed, such as GRB 980425/SN 1998bw (Galama et al. Reference Galama1998), GRB 030329/SN 2003dh (Stanek et al. Reference Stanek2003), GRB 031203/SN 2003lw (Malesani et al. Reference Malesani2004; Thomsen et al. Reference Thomsen2004), GRB 091127/SN 2009nz (Cobb et al. Reference Cobb2010), GRB 101219B/SN 2010ma (Sparre et al. Reference Sparre2011), GRB 120422A/SN 2012bz (Melandri et al. Reference Melandri2012), GRB 130427A/SN 2013cq (Xu et al. Reference Xu2013), and GRB 130702A/SN 2013dx (D’Elia et al. Reference D’Elia2015). Greiner et al. (Reference Greiner2015) suggested that GRB 111209A/SN 2011kl is an SLSN though its peak magnitude is –20 and it is dimmer than the usual SLSNe. It may come from a sub-group as a ‘blue supergiant’ (Nakauchi et al. Reference Nakauchi, Kashiyama, Suwa and Nakamura2013). There might be a link between SLSNe and GRB-connected SNe. By accumulating more and more samples, one can tell what fractions of the long GRBs are associated with SNe Ic, and what fractions of the SNe Ic are associated with GRBs (Guetta & Della Valle Reference Heger, Fryer, Woosley, Langer and Hartmann2007). Soderberg, Frail, & Wieringa (Reference Soderberg, Frail and Wieringa2004) suggested that less than 6% of the SNe Ic produces relativistic jets by searching for radio emission from misaligned jets. This may reveal the nature that which properties are crucial for producing a pair of relativistic jets.
Figure 4 shows the K–S test for the peak magnitude of GRB-associated SNe and normal SNe Ic. Bloom et al. (Reference Bloom, Frail and Kulkarni2003) (Figure 10 in Piran Reference Piran2004) plotted the comparison between the GRB-associated SNe and the local normal SNe Ib/c. The distributions were quite different. However, the data were very limited. Similar cumulative plots have been shown in Hjorth & Bloom (Reference Hjorth, Bloom, Kouveliotou, Wijers and Woosley2012) and Hjorth (Reference Hjorth2013), and the relatively brighter peak luminosities for GRB-associated SNe were suggested to be the result of a bias against the faint systems. Cano et al. (Reference Cano2011a) found the distributions are similar by neglecting those samples not having host galaxy extinction information. Here we assume the sample is complete, as some of the GRB afterglows have been observed in a deep flux, such as GRBs 060605 and 060614 (Fynbo et al. Reference Fynbo2006). We increased the peak magnitude of GRB-associated SNe by 1 magnitude, to look for any similarity between these two samples. It shows that these two samples may share a similar distribution, which implies a similar origin, for example, a similar range of the massive star. A possible explanation might be that the SN component is dominated by the mass of the progenitor, while the jet is dominated by the rotation of the progenitor, and the rotation is not related to the mass. Therefore, the SN with or without a jet shows a similar distribution. Notice the magnitude of GRB-associated SNe is increased by 1, which means they are relatively brighter than the normal SNe Ic. The reason might be that, for the SNe with jets, part of the jet energy may deposit into the envelope, finally ending up as the optical luminosity. However, this constant off-set is likely observational bias. The normal SN Ic in the distance of GRB-associated SNe will not be found by the routine survey. Once they are associated with GRB, an extensive search will be carried so apparent dimmer but intrinsic more powerful SNe can be detected. Furthermore, the beaming effect of GRBs makes them difficult to detect in the smaller nearby volume but more easily detected at larger cosmological distances. As a conclusion, the similar distribution of the magnitude of the normal SNe Ic and the modified magnitude (1 magnitude added) of the SNe Ic with GRB connection is interesting, but the reason might be complex and the similarity might be a coincidence.
K–S test for the peak magnitude of GRB-associated SNe (red) and normal SNe Ic (cyan), of which the GRB-associated SNe are increased by 1 magnitude. The red dashed line indicates the maximum difference. Data are taken from the same source as in Figure 3.
2.3. SN Ic-XRF
It is possible that the jet cannot fully penetrate the envelope, either because the envelope is too thick or the energy of the jet is relatively low. However, the jet is still energetic enough to power a relativistic shock breakout (see Nakar & Sari (Reference Nakar and Sari2010, Reference Nakar and Sari2012) for an extensive modelling). This breakout may radiate thermal X-rays, which can be detected as an XRF or an X-ray-rich GRB (or presents as a low-luminosity GRB (Barniol Duran et al. Reference Barniol Duran, Nakar, Piran and Sari2015). The light curve of the XRF is smooth compared with a normal GRB, and the spectrum is thermal rather than non-thermal. Several XRF/SN systems have been observed, such as GRB 060218/SN 2006aj (Campana et al. Reference Campana2006; Pian et al. Reference Pian2006; Soderberg et al. Reference Soderberg, Nakar, Berger and Kulkarni2006b) and GRB 100316D/SN 2010bh (Cano et al. Reference Cano2011b; Bufano et al. Reference Bufano2012). This phenomenon has attracted greatest interests and is successfully explained by a shock breakout model with some tuning or modifications (Campana et al. Reference Campana2006; Waxman, Mészáros, & Campana Reference Waxman, Mészáros and Campana2007; Li Reference Li2007; Nakar Reference Nakar2015) (see Irwin & Chevalier (Reference Irwin and Chevalier2016) however for an alternative jet model).
The off-axis case has probably been observed already, such as SN 2007gr (Paragi et al. Reference Paragi2010) and SN 2009bb (Soderberg et al. Reference Soderberg, Frail and Wieringa2010b). For SN 2007gr, a mildly relativistic jet directed far away from the observer was suggested (Paragi et al. Reference Paragi2010; Xu, Nagataki, & Huang (Reference Xu, Nagataki and Huang2011).Footnote g Very recently, numerical simulation on the radio afterglow of SN 2009bb also showed the consistency to the off-axis scenario (De Colle, Kumar, & Aguilera-Dena Reference De Colle, Kumar and Aguilera-Dena2018).
SN 2008D (Mazzali, Iwamoto, & Nomoto Reference Mazzali2008) is an SN Ib with a shock breakout (Chevalier & Fransson Reference Chevalier and Fransson2008), which is connected to XRF 080109.Footnote h This may be a link among SNe Ib, SNe Ic, and GRBs. For SNe Ib, as the existence of helium envelopes, the jets, if they exist, penetrate more hardly. It can only appear as an XRF rather than a normal GRB.
2.4. SN Ic-BL
Besides powering the relativistic jets, the central engine may also provide energy to the envelope of the progenitor. This extra energy may result in the ejecta of the SN having a higher speed compared to normal SN without a pair of jets. This is why all the SNe Ic connected with GRBs/XRFs are broad-line SNe. GRBs are highly beaming events. That means a great part of GRBs are not observed because of the direction of the jet, while the percentage depends on the jet angle. However, the beaming effect for the SN component is negligible. Therefore, an off-axis GRB jets cannot be observed as a GRB, but the SN can still present as a broad-line SN. A similar effect applies to XRF-connected SN Ic. However, as there is no relativistic jet propagating into the surrounding medium, one cannot expect an orphan afterglow.
Several SNe-BL have been observed, such as SN 2002ap (Kinugasa et al. Reference Kinugasa2002), SN 2003jd (Mazzali et al. Reference Mazzali2005), SN 2010ah (Mazzali et al. Reference Mazzali2013), SN 2012ap (Margutti et al. Reference Margutti2014), and PTF10qts (Walker et al. Reference Walker2014). People have paid high attention to the probability that these events are also off-beamed GRBs, such as Kawabata et al. (Reference Kawabata2002), Totani (Reference Totani2003, for SN 2002ap), and other SNe-BL (Soderberg et al. Reference Soderberg, Nakar, Berger and Kulkarni2006a). However, there is no strong evidence for the existence of relativistic jets, though Mazzali et al. (Reference Mazzali2005) suggested SN 2003jd is an SN with off-axis GRB jets from the double-peaked emission line. Recently, Corsi et al. (Reference Corsi2016) searched the radio emission from the SN Ic-BL trying to find a signal from the off-axis jet, setting a probability ≤45% of being associated with a GRB. This scenario can be confirmed if an off-axis emission of a relativistic jet associated with an SN Ic-BL is detected, either in the optical band (Nakar, Piran, & Granot Reference Nakar, Piran and Granot2002; Zou et al. Reference Zou, Wu and Dai2007) or in the radio band (Levinson et al. Reference Levinson, Ofek, Waxman and Gal-Yam2002). The reason for the missing optical and radio signals might be due to the high beaming effect of the jet from a GRB, or even the absence of a jet from an XRF. For example, based on the similarity between SN 2003jd and SN 2006aj, Valenti et al. (Reference Valenti2008) suggested SN 2003jd may also be connected with an XRF.
The birth rate of the SN Ic-BL should be related to the SNe Ic connected with GRBs/XRFs, as these are just geometrical effects. Lacking data, here we apply a very rough estimation similar to that used in equation (1), that is, we take the distance D of the nearest object as the average distance between two objects and get the rate 
$\propto {D^{ - 3}}$. As shown in Figure 3, the nearest SN Ic-BL is 9.2 Mpc for SN 2002ap (Barbon et al. Reference Barbon, Buondí, Cappellaro and Turatto1999), and the nearest SN Ic-GRB is 35 Mpc for SN 1998bw (Galama et al. Reference Galama1998). (The nearest SN Ic-XRF is 139 Mpc for SN 206aj (Campana et al. Reference Campana2006), much farther than the nearest SN Ic-GRB. Hence, the birth rate of SN Ic-XRF is negligible comparing with the birth rate of SN Ic-GRB.) A comparison of the birth rate of the SNe Ic-GRBs and the SNe Ic-BL is ∼(9.2 Mpc)3/(35 Mpc)3 = 0.018. Considering the GRBs are beamed because of the relativistic jet, while SNe Ic-BL can be seen in any direction, this rate indicates that the solid angle of the jet (with half opening angle 
${\theta _j}$): 
${{\theta _j^2} \over 2} \sim 0.018$. This number is consistent with the estimation of the solid angle being 10−3−10−2 from the jet opening angle (Kumar & Zhang Reference Kumar and Zhang2015).
2.5. SLSN Ic
If the envelope of the progenitor star is too massive, the jets can never overcome the envelope. Then all the energy of the jet will be deposited into the envelope. The mass–radius (M−R) relation of a Wolf-Rayet (WR) star is given by Schaerer & Maeder (Reference Schaerer and Maeder1992),
\begin{equation} \log {R \over {{{\rm{R}}_ \odot}}} = - 0.6629 + 0.5840\log {M \over {{{\rm{M}}_ \odot}}}, \end{equation}
where R⊙ and M⊙ are the radius and mass of the Sun, respectively. Taking SLSN 2007bi as an example, the inferred mass of the WR star of SN 2007bi is ≃43M⊙ (Moriya et al. Reference Moriya, Tominaga, Tanaka, Maeda and Nomoto2010). Consequently, the radius is ∼2R⊙. Notice the shocked envelope propagates in mildly relativistic velocity (Zhang et al. Reference Zhang, Woosley and MacFadyen2003; Zhang, Woosley, & Heger Reference Zhang, Woosley and Heger2004b). By taking 0.1c as the velocity, it takes the jet ∼50 s to pass through the envelope. Roughly we can say that only if the central engine last longer than ∼50 s, it can support the jet to pass through the envelope. The typical timescale of the long GRBs is ∼26 s (Kouveliotou et al. Reference Kouveliotou1993), which can marginally fail to support the jet in overcoming the envelope. Therefore, most of the jet energy is deposited into the envelope. On the other hand, the WR star mass of the SLSN is about 10 times larger than that of the GRB progenitor star, which also makes it reasonable that the jets cannot pass through the envelope. In Figure 3, it appears that SLSN Ic have a larger intrinsic power than those GRB-connected SN Ic. However, as we have emphasised in our picture that GRB-connected SN Ic are those with a jet penetrating the envelope, whereas SLSN Ic are those where the jet energy is deposited in the envelope. If we take into account the GRB energy in the GRB-connected SN Ic, the combined energy is actually similar to that of SLSN Ic.
The mass of the initial jet from the central engine is the same as the jet from a normal GRB, of which the total energy is 1051−1052 erg, the initial Lorentz factor of the jet is ∼100, and consequently the mass is about 0.01 M⊙ (see MacFadyen, Woosley, & Heger Reference MacFadyen, Woosley and Heger2001, for example). With the propagation of the jet inside the envelope, more materials are accumulated onto the jet and it can be slowed down to about 0.1c. Considering the total energy is fully deposited into the kinetic energy of the envelope, the accumulated mass of the pair of jets is consequently about 0.1−1 M⊙. Considering the extreme case of SN 2007bi, the kinetic energy was suggested being 3.6 × 1052 erg (Moriya et al. Reference Moriya, Tominaga, Tanaka, Maeda and Nomoto2010). The corresponding mass of the pair of jets is roughly 4 M⊙ if the velocity is taken as 0.1c. Taking the full mass of the envelope as 40 M⊙ as suggested in Moriya et al. (Reference Moriya, Tominaga, Tanaka, Maeda and Nomoto2010), which indicates the opening solid angle of a jet is roughly 
${\pi \over 5}$. That is, to acquire a mass of 2 M⊙ along its motion through the star, the opening angle (from symmetry axis to the edge) of each jet should be 26°. However, if the kinetic energy is deposited into a wider angle in the envelope, the accumulated mass is heavier and the final average velocity is smaller. Notice that this opening angle is not the same as the opening angle of the relativistic jet just launched from the central engine. See Figure 1 in Nakar & Piran (Reference Nakar and Piran2017) for a schematic view.Footnote i Each relativistic jet has much smaller opening angle. With the propagation in the envelope, all the kinetic energy of the relativistic jet is deposited into the envelope. The outflow becomes non-relativistic and the opening angle becomes much larger.
The energy from the jets may enhance the luminosity of the SN dramatically, either by nucleosynthesis in terms of storing the energy in 56Ni (e.g. Arnett Reference Arnett1982; Chatzopoulos et al. Reference Chatzopoulos, Wheeler, Vinko, Horvath and Nagy2013) or by being stored as thermal energy of the envelope (e.g. Couch, Wheeler, & Milosavljević Reference Couch, Wheeler and Milosavljević2009) or by transferring into kinetic energy of the envelope and then interacting with the circumstellar medium (CSM) (e.g. Moriya et al. Reference Moriya, Tominaga, Blinnikov, Baklanov and Sorokina2011; Chatzopoulos et al. Reference Chatzopoulos, Wheeler, Vinko, Horvath and Nagy2013). Together with the fact that the massive star itself may produce a brighter SN, the final result may present as an SLSN. However, the observations make it difficult to distinguish between the Ib and Ic. The SLSNe are divided into three classes, that is, radioactively powered SLSN-R, hydrogen-poor SLSN-I, and hydrogen-rich SLSN-II (Gal-Yam Reference Gal-Yam2012), where SLSN-R is also type Ic. Recently, there have been several SLSN-R detected, such as SN 1999as (Gal-Yam Reference Gal-Yam2012), SN 2007bi (Gal-Yam et al. Reference Gal-Yam2009), PTF12dam (Nicholl et al. Reference Nicholl2013), PS1-11ap (McCrum et al. Reference McCrum2014), and SN 2015bn (Nicholl et al. Reference Nicholl2016). In Figure 3, we only take SLSN-R into account. There are three leading models: magnetar spin down, 56Ni decay, and ejecta-CSM interaction (Chatzopoulos et al. Reference Chatzopoulos, Wheeler, Vinko, Horvath and Nagy2013). For the last two models, the central engine could be the energy injection from the GRB jets.
Recently, a very luminous SLSN ASASSN-15lh has been discovered with a peak bolometric luminosity (2.2 ± 0.2) × 1045 erg s−1 (Dong et al. Reference Dong2016). The magnetar model can hardly provide so much energy, and a quark star model has to be invoked (Dai et al. Reference Dai, Wang, Wang, Wang and Yu2016). On the other hand, there is no such energy limit for a jet driving an SLSN. Very recently, Chen et al. (Reference Chen2017) found the metallicity of the progenitors of the SLSN is relatively lower. This is consistent with the low metallicity stars having more extensive envelopes.
The total radiated energy of SLSNe is of the order of 1051 erg (Gal-Yam Reference Gal-Yam2012). If they are powered by the relativistic jets, the energy should be comparable with GRBs, for which the radiated energy is truly also of the same order of 1051 erg (Frail et al. Reference Frail2001). Here we assume the unknown radiation efficiency from jet energy to γ-rays of GRBs is of the same order as that of the jet energy to optical photons of SNe. Considering the range of the jet energy, which varies from 1.1 × 1050 to 3.3 × 1053 erg (Shivvers & Berger Reference Shivvers and Berger2011), we predict that, with accumulating data of SLSNe, there should be a fraction with energy up to 1053 erg.
2.6. Normal SN Ic
A type Ic SN is both H and He poor in the spectrum, which indicates the progenitor is a very massive star (Smartt Reference Smartt2009). The peak absolute magnitude is about –18, and the kinetic energy is about 1051 erg (Drout et al. Reference Drout2011). There is no evidence that a normal SN Ic has any jet. Berger et al. (Reference Berger, Kulkarni, Frail and Soderberg2003) carried out a radio search of 33 SNe Ib and Ic observed from 1999 to 2002, and none of them showed clear jet emission. It is also possible that there is a jet as in the SLSN case, and the total energy of the jet is much less than 1051 erg. Consequently, the deposited energy has no observational effect comparable to a normal SN Ic.
Smartt (Reference Smartt2009) showed a picture where a fast-rotating WR star may produce an SN Ibc-BL, while a slow-rotating WR star will produce a normal SN. This is consistent with our picture, as a rotating WR star may also launch a pair of jets.
3. Conclusions and discussions
Relating two or three sub-classes of SNe Ic has been suggested widely in the literatures. In this paper we suggested an inner connection for all different sub-classes of SNe Ic. It is mainly based on the fate of the central engine launched jets. Roughly speaking, the more energy of the jets that is deposited into the envelope of the progenitor, the stronger the luminosity of the SN component. In addition to collecting previous precise evidence support of our full picture, we also find new evidences in support of this inner connection. Figure 2 schematically illustrates a full picture of the inner connection. We also propose a luminosity-distance diagram and use it to indicate the existence of a single luminosity function for all different types of SNe Ic, which obeys some inner connection. We estimate that SLSN Ic and GRB-associated SN Ic have similar intrinsic power if the energy of the GRB is included. This evidence supports our picture that the difference between SLSN and GRB-associated SNe is due to the absorption of the jet energy. As there is a small part of the GRBs which are much more energetic, we also predicted that much more luminous SLSNe will be observed with accumulating data.
This full picture also shows a hint why the GRBs only associate with type Ic SNe. The type II core-collapsing SNe may also launch a pair of jets. The critical parameters to launch jets might be the mass of the initial main sequence star, its magnetic field, spin, etc. But the most important factor is that for the type II SNe, they should have a more extensive envelope due to the hydrogen shell. This prevents the jets penetrating the envelope, and a relativistic jet cannot form. Consequently, it cannot produce a GRB. Maybe XRF 080109 is an example that the jets are launched from a type II SN, but have not penetrated. The other hint is the SLSNe with type II SNe. If the jets cannot pass through the envelope, but the energy is high enough, they may still power an SLSN, which is type II.
More detailed models should consider the extent of the energy budget for individual events. Failed jets may produce SLSN. While if it is not so energetic, it may produce a normal SN Ic-BL. Therefore, the type SN Ic-BL may not only be associated with off-beamed GRBs/XRFs. If the energy of the jets is even weaker, say ≪1051 erg, the jets may not even penetrate a normal SN Ic progenitor. Because of the jet energy is negligible, there is no obvious observational evidence, and they are indistinguishable from an SN Ic with no jets launched. The inner engine could be sufficiently energetic to penetrate an extensive envelope. It is possible that part of the energy of the jets penetrate into the envelope, while the jets still pass through it. Consequently, the connected SN is brighter than a normal SN, but dimmer than an SLSN, as seen in SN 2011kl (connection with GRB 111209A) (Greiner et al. Reference Greiner2015; Kann et al. Reference Kann2016). The progenitor of GRB 111209A is also suggested as a blue supergiant by Gendre et al. (Reference Gendre2013). These might be minorities.
We also notice there are some long GRBs, which have been observed in a deep flux limit but the SN component did not appear, such as GRB 060605, GRB 060614 (Fynbo et al. Reference Fynbo2006), and GRB 100418A (Niino et al. Reference Niino2012). This is either because the associated SN essentially has low luminosity (Valenti et al. Reference Valenti2009) or because the GRB belongs to other sub-classes, such as, a merger-originated GRB (Zhang et al. Reference Zhang2007), an intermediate-mass BH-powered GRB (Gao, Lu, & Zhang Reference Gao, Lu and Zhang2010), or a macronova-associated GRB (Yang et al. Reference Yang2015). According to these explanations, they are very likely not correlated with a core-collapsing SN. Future gravitational wave detection [such as LISA (Bender et al. Reference Bender1998) and TianQin (Luo et al. Reference Luo2016)] should be a key discriminator for the single or double origins of this kind of long GRBs. An alternative possibility is that, the progenitor star collapses directly into a BH rather than an NS first. With no NS stage, the SN cannot form, while only jet-launched GRB appears.
It is unclear if the total energy of the jets and isotropic ejecta are the same order of magnitude for all SNe Ic due to the diversity of SNe’s energy reservoirs. These two should be compatible with either BH powering or NS powering central engines only if they can power a pair of relativistic jets. Therefore, it is compatible for the magnetar models. Whether the jets can be produced may depend on the spin and the magnetic field.
In this paper we have proposed a simple full picture of inner connection for the normal SNe Ic, SLSNe Ic, XRF-connected SNe, SN Ic-BL, and GRB-connected SNe. We want to remark that we have selected the interpretations of some observed cases favouring our scenario. In fact there are some other possible interpretations, for example, Pignata et al. (Reference Pignata2011) argued that the type SN Ic-BL, SN 2009bb is not an off-beamed GRB. However, all the different interpretations are far from decisive and are not always consistent with each other. A selection of some interpretations to provide a logically clear scenario is still helpful to understand the death of the massive stars.
Author ORCIDs.
Y.-C. Zou, http://orcid.org/0000-0002-5400-3261
Acknowledgements
Y. C. Z. thanks the helpful discussions with Tsvi Piran, Zigao Dai, Yongfeng Huang, Xiangyu Wang, Zhuo Li, Xuefeng Wu, Yunwei Yu, Fayin Wang, Weihua Lei, Jumpei Takata, Dong Xu, Shanqin Wang, and Wei Xie. Figures 2 and 3 are plotted by using Google Presentation and R, respectively. We are grateful for the anonymous reviewer for the insightful and detailed comments, and for Prof. Kevin Mackeown for a critical reading of the whole manuscript. K. S. C. is supported by a GRF grant of Hong Kong Government under 17310916. Y. C. Z. is supported by the National Basic Research Program of China (973 Program, Grant No. 2014CB845800) and by the National Natural Science Foundation of China (Grant Nos. U1738132 and U1231101).
