Induced Gravitational Collapse, Binary-Driven Hypernovae, Long Gramma-ray Bursts and Their Connection with Short Gamma-ray Bursts

TL;DR

This paper reviews the binary-driven hypernovae model linking long and short gamma-ray bursts through complex supernova and accretion processes, emphasizing the role of induced collapse and binary mergers in GRB phenomena.

Contribution

It provides a comprehensive overview of the BdHNe model, integrating analytic estimates and 3D simulations to explain GRB features and their connection to binary mergers.

Findings

Supernova explosions produce a new neutron star and trigger hypercritical accretion.

The formation of black holes from neutron stars explains various GRB phases.

Binary mergers produce gravitational waves linking long and short GRBs.

Abstract

Short and long Gamma-ray bursts (GRBs) originate in subclasses with specific energy release, spectra, duration, etc, and have binary progenitors. We review here the binary-driven hypernovae (BdHNe) subclass whose progenitor is a CO$_\textrm{core}$-neutron star (NS). The supernova (SN) explosion of the CO$_\textrm{core}$ produces at its center a new NS ($\nu$NS) and triggers a hypercritical accretion onto the NS. The NS can become a more massive NS or collapse into a black hole (BH). We summarize this topic from the first analytic estimates in 2012 to the most recent three-dimensional (3D) smoothed-particle-hydrodynamics (SPH) numerical simulations in 2018. Long GRBs are richer and more complex than previously thought. The SN and the accretion explain X-ray precursors. The NS accretion, its collapse and the BH formation produce asymmetries in the SN ejecta, implying a 3D GRB analysis. The newborn BH surrounded by the ejecta and the magnetic field inherited from the NS, are the \emph{inner engine} from which the electron-positron ($e^+e^-$) plasma and the high-energy emission initiate. The $e^+e^-$ impact on the ejecta converts the SN into a hypernova (HN). The plasma dynamics in the ejecta explains the ultrarelativistic prompt emission in the MeV domain and the mildly-relativistic flares of the early afterglow in the X-ray domain. The feedback of the $\nu$NS emission on the HN explains the X-ray late afterglow and its power-law regime. All the above is in contrast with GRB models attempting to explain all the GRB phases with the kinetic energy of anultrarelativistic jet, as traditionally proposed in the "collapsar-fireball" model. In addition, BdHNe in their different flavors lead to $\nu$NS-NS or $\nu$NS-BH binaries. These binaries merge by gravitational wave emission producing short GRBs, establishing a connection between long and short GRBs and their occurrence rates.

Tables

Table 1. Table 1: Summary of the Gamma-ray bursts (GRB) subclasses. This table is an extended version of the one presented in  Wang et al. ( 2019 ) with the addition of a column showing the local density rate, and it also updates the one in  Ruffini et al. ( 2016 , 2018a ) . We unify here all the GRB subclasses under two general names, BdHNe and BMs. Two new GRB subclasses are introduced; BdHN Type III and BM Type IV. In addition to the subclass name in “Class” column and “Type” column, as well as the previous names in “Previous Alias” column, we report the number of GRBs with known redshift identified in each subclass updated by the end of 2016 in “number” column (the value in a bracket indicates the lower limit). We recall as well the “in-state” representing the progenitors and the “out-state” representing the outcomes, as well as the the peak energy of the prompt emission, E p , i subscript 𝐸 p i E_{\rm p,i} , the isotropic Gamma-ray energy, E iso subscript 𝐸 iso E_{\rm iso}  defined in the 1 1 1  keV to 10 10 10  MeV energy range, the isotropic emission of ultra-high energy photons, E iso , Gev subscript 𝐸 iso Gev E_{\rm iso,Gev} , defined in the 0.1 0.1 0.1 – 100 100 100 GeV energy range, and the local observed rate ρ GRB subscript 𝜌 GRB \rho_{\rm GRB}   Ruffini et al. ( 2016 ) . We adopt as definition of kilonova a phenomenon more energetic than a nova (about 1000 times). A kilonova can be an infrared-optical counterpart of a NS-NS merger. In that case the transient is powered by the energy release from the decay of r-process heavy nuclei processed in the merger ejecta  Li and Paczyński ( 1998 ); Metzger et al. ( 2010 ); Tanvir et al. ( 2013 ); Berger et al. ( 2013 ) . FB-KN stands for fallback-powered kilonova  Rueda et al. ( 2018 , 2019 ) : a WD-WD merger can emit an infrared-optical transient, peaking at ∼ similar-to \sim 5 day post-merger, with the ejecta powered by accretion of fallback matter onto the newborn WD formed in the merger. The density rate of the GRB subclasses BdHN III (HN) and BM IV (FB-KN) have not yet been estimated.

Class	Type	Previous Alias	Number	In-State	Out-State	${𝑬}_{𝐩,𝐢}$ (MeV)	${𝑬}_{\mathrm{𝐢𝐬𝐨}}\text{(erg)}$	${𝑬}_{\mathrm{𝐢𝐬𝐨},\mathrm{𝐆𝐞𝐯}}$ (erg)	${𝝆}_{\mathrm{𝐆𝐑𝐁}}$ (Gpc-3 yr-1)
Binary-driven	I	BdHN	$329$	COcore-NS	$\nu$NS-BH	$\sim$0.2–$2$	$\sim$1052–${10}^{54}$	$\gtrsim$1052	${0.77}_{-0.08}^{+0.09}$
hypernova	II	XRF	$\left(30\right)$	COcore-NS	$\nu$NS-NS	$\sim$0.01–$0.2$	$\sim$1050–${10}^{52}$	$-$	${100}_{-34}^{+45}$
(BdHN)	III	HN	$\left(19\right)$	COcore-NS	$\nu$NS-NS	$\sim$0.01	$\sim$1048–${10}^{50}$	$-$	$-$
	IV	BH-SN	$5$	COcore-BH	$\nu$NS-BH	$\gtrsim$2	$>$1054	$\gtrsim$1053	$\lesssim$0.77${}_{-0.08}{}^{}{}_{}{}^{+0.09}$
	I	S-GRF	$18$	NS-NS	MNS	$\sim$0.2–$2$	$\sim$1049–${10}^{52}$	$-$	${3.6}_{-1.0}^{+1.4}$
Binary	II	S-GRB	$6$	NS-NS	BH	$\sim$2–$8$	$\sim$1052–${10}^{53}$	$\gtrsim$1052	$\left({1.9}_{-1.1}^{+1.8}\right)\times {10}^{-3}$
Merger	III	GRF	$\left(1\right)$	NS-WD	MNS	$\sim$0.2–$2$	$\sim$1049–${10}^{52}$	$-$	${1.02}_{-0.46}^{+0.71}$
(BM)	IV	FB-KN	$\left(1\right)$	WD-WD	NS/MWD	0.2	1051	$-$	$-$
	V	U-GRB	$\left(0\right)$	NS-BH	BH	$\gtrsim$2	$>$1052	$-$	$\approx$0.77${}_{-0.08}{}^{}{}_{}{}^{+0.09}$

Table 2. Table 2: Acronyms used in this work in alphabetical order.

Extended Wording	Acronym
Binary-driven hypernova	BdHN
Black hole	BH
Carbon-oxygen core	COcore
Gamma-ray burst	GRB
Gamma-ray flash	GRF
Induced gravitational collapse	IGC
Massive neutron star	MNS
Neutron star	NS
New neutron star created in the SN explosion	$\nu$NS
Short Gamma-ray burst	S-GRB
Short Gamma-ray flash	S-GRF
Supernova	SN
Ultrashort Gamma-ray burst	U-GRB
Ultra high-energy cosmic ray	UHECR
White dwarf	WD
X-ray flash	XRF

Table 3. Table 3: Critical NS mass in the non-rotating case and constants k 𝑘 k and p 𝑝 p needed to compute the NS critical mass in the non-rotating case given by Equation ( 6 ). The values are for the NL3, GM1 and TM1 EOS.

EOS	${𝑴}_{\mathrm{𝐜𝐫𝐢𝐭}}^{𝑱=\mathrm{𝟎}}$ $\left({𝑴}_{\odot }\right)$	$𝒑$	$𝒌$
NL3	$2.81$	$1.68$	$0.006$
GM1	$2.39$	$1.69$	$0.011$
TM1	$2.20$	$1.61$	$0.017$

Equations

$$\begin{array}[]{rclcrcl}\dot{M}_{B}(t)&=&\pi\rho_{\rm ej}R_{{\rm cap}}^{2}\sqrt{v_{{\rm rel}}^{2}+c_{\rm s,ej}^{2}},\qquad R_{{\rm cap}}(t)&=&\cfrac{2GM_{{\rm NS}}(t)}{v_{{\rm rel}}^{2}+c_{{\rm s,ej}}^{2}},\end{array}$$

$${\rho }_{\mathrm{ej}}(r,t)={\rho }_{\mathrm{ej}}^0(r/R_{\mathrm{star}}(t),t_0)\frac{M_{\mathrm{env}}(t)}{M_{\mathrm{env}}(0)}{\left(\frac{R_{\mathrm{star}}(0)}{R_{\mathrm{star}}(t)}\right)}^3,$$

$${\dot{M}}_{\mathrm{NS}}={\left(\frac{\partial M_{\mathrm{NS}}}{\partial M_b}\right)}_{J_{\mathrm{NS}}}{\dot{M}}_b+{\left(\frac{\partial M_{\mathrm{NS}}}{\partial J_{\mathrm{NS}}}\right)}_{M_b}{\dot{J}}_{\mathrm{NS}},{\dot{J}}_{\mathrm{NS}}=\xi l(r_{\mathrm{in}}){\dot{M}}_{\mathrm{B}},$$

$$\frac{M_b}{M_{\odot }}=\frac{M_{\mathrm{NS}}}{M_{\odot }}+\frac{13}{200}{\left(\frac{M_{\mathrm{NS}}}{M_{\odot }}\right)}^2\left(1-\frac{1}{137}{j}_{\mathrm{NS}}^{1.7}\right),$$

$$l_{\mathrm{lco}}=\frac{G{M}_{\mathrm{NS}}}{c}\left[2\sqrt{3}-0.37{\left(\frac{j_{\mathrm{NS}}}{M_{\mathrm{NS}}/M_{\odot }}\right)}^{0.85}\right].$$

$$M_{\mathrm{NS}}^{\mathrm{crit}}=M_{\mathrm{NS}}^{J=0}(1+k{j}_{\mathrm{NS}}^p),$$

$$S_{\mathrm{bubble}}\approx 16{\left(\frac{1.4M_{\odot }}{M_{\mathrm{NS}}}\right)}^{-7/8}{\left(\frac{M_{\odot }{\mathrm{s}}^{-1}}{{\dot{M}}_{\mathrm{B}}}\right)}^{1/4}{\left(\frac{10^6\mathrm{cm}}{r}\right)}^{3/8}{k}_B/\mathrm{nucleon},$$

$$F_{\mathrm{bubble}}\approx 2\times 10^{40}{\left(\frac{M_{\mathrm{NS}}}{1.4M_{\odot }}\right)}^{-7/2}\left(\frac{{\dot{M}}_{\mathrm{B}}}{M_{\odot }{\mathrm{s}}^{-1}}\right){\left(\frac{R_{\mathrm{NS}}}{10^6\mathrm{cm}}\right)}^{3/2}{\left(\frac{r_0}{r}\right)}^{8/3}{\mathrm{erg}\mathrm{s}}^{-1}{\mathrm{cm}}^{-2}.$$

$$T_{\mathrm{acc}}\approx {\left(\frac{3P_{\mathrm{shock}}}{4\sigma /c}\right)}^{1/4}={\left(\frac{7}{8}\frac{{\dot{M}}_{\mathrm{acc}}{v}_{\mathrm{acc}}c}{4\pi R_{\mathrm{NS}}^2\sigma }\right)}^{1/4},$$

$${ϵ}_{e^{-}{e}^{+}}\approx 8.69\times 10^{30}{\left(\frac{k_{B}T}{1\mathrm{MeV}}\right)}^9\mathrm{MeV}{\mathrm{cm}}^{-3}{\mathrm{s}}^{-1},$$

$$\Delta r_{\nu }=\frac{{ϵ}_{e^{-}{e}^{+}}}{\nabla {ϵ}_{e^{-}{e}^{+}}}=\frac{\Delta r_{\mathrm{ER}}}{9}\approx 0.08R_{\mathrm{NS}}.$$

$${\dot{M}}_{\mathrm{eff}}\approx \Delta M_{\nu }\frac{L_{\nu }}{E_g},$$

$$L_{\nu }\approx 4\pi R_{\mathrm{NS}}^2\Delta r_{\nu }{ϵ}_{e^{-}{e}^{+}},$$

$$L_{\mathrm{acc}}=({\dot{M}}_b-{\dot{M}}_{\mathrm{NS}})c^2={\dot{M}}_b{c}^2\left[1-{\left(\frac{\partial M_{\mathrm{NS}}}{\partial J_{\mathrm{NS}}}\right)}_{M_b}l-{\left(\frac{\partial M_{\mathrm{NS}}}{\partial M_b}\right)}_{J_{\mathrm{NS}}}\right].$$

$$\Delta E_{\mathrm{acc}}\equiv \int L_{\mathrm{acc}}dt,$$

$$E_{r_{+}}\approx \frac{1}{2}\alpha B_0=6.5\times 10^{15}\cdot \alpha \left(\frac{B_0}{B_c}\right)\frac{\mathrm{V}}{\mathrm{cm}},$$

$$\alpha (B_0/B_c)\ge 2,B_0/B_c\ge 2,$$

$$\frac{{\beta }^5}{4[\ln (1+\beta )-(1-\beta )\beta {]}^2}{\left(\frac{1+\beta }{1-\beta }\right)}^{1/2}=\frac{D_L(z)}{1+z}\frac{1}{t_2-t_1}\left(\sqrt{\frac{F_{\mathrm{bb},\mathrm{obs}}(t_2)}{\sigma T_{\mathrm{obs}}^4(t_2)}}-\sqrt{\frac{F_{\mathrm{bb},\mathrm{obs}}(t_1)}{\sigma T_{\mathrm{obs}}^4(t_1)}}\right),$$

$$\Delta \varphi =-{\int }_{\infty }^{r_{+}}Edr=E_{r_{+}}{r}_{+}=9.7\times 10^{20}\cdot \alpha \left(\frac{B_0}{B_c}\right)\left(\frac{M}{M_{\odot }}\right)(1+\sqrt{1-{\alpha }^2})\frac{\mathrm{V}}{e},$$

$$\mathcal{E}=\frac{1}{2}{E}_{r_{+}}^2{r}_{+}^3\approx 7.5\times 10^{41}\cdot {\alpha }^2{\left(\frac{B_0}{B_c}\right)}^2{\left(\frac{M}{M_{\odot }}\right)}^3(1+\sqrt{1-{\alpha }^2}{)}^3\mathrm{erg},$$

$$N_p=\frac{\mathcal{E}}{{ϵ}_p}\approx 4.8\times 10^{32}\alpha \left(\frac{B_0}{B_c}\right){\left(\frac{M}{M_{\odot }}\right)}^2(1+\sqrt{1-{\alpha }^2}{)}^2.$$

$$\Delta t_{\mathrm{el}}=\frac{\Delta \varphi }{E_{r_{+}}c}=\frac{r_{+}}{c}\approx 4.9\times 10^{-6}\left(\frac{M}{M_{\odot }}\right)(1+\sqrt{1-{\alpha }^2})\mathrm{s}.$$

$$\frac{d\mathcal{E}}{dt}\approx \frac{\mathcal{E}}{\Delta t_{\mathrm{el}}}=1.5\times 10^{47}\cdot {\alpha }^2{\left(\frac{B_0}{B_c}\right)}^2{\left(\frac{M}{M_{\odot }}\right)}^2(1+\sqrt{1-{\alpha }^2}{)}^2\mathrm{erg}\cdot {\mathrm{s}}^{-1}.$$

$$\frac{a}{a_0}=\frac{M_0-\Delta M}{M_0-2a_0\Delta M/r},$$

Full text

\pubvolume

5 \issuenum5 \articlenumber110

\historyReceived: 25 February 2019; Accepted: 5 May 2019; Published: 9 May 2019 \updatesyes

\TitleInduced Gravitational Collapse, Binary-Driven Hypernovae, Long Gramma-ray Bursts and Their Connection with Short Gamma-ray Bursts \AuthorJ. A. Rueda *1,2,3,**, R. Ruffini 1,2,4 and Y. Wang 1,2

\corresCorrespondence: jorge.rueda@icra.it

1 Introduction

1.1 The Quest for the Binary Nature of GRB Progenitors

We first recall that GRBs have been traditionally classified by a phenomenological division based on the duration of the time-interval in which the 90% of the total isotropic energy in Gamma-rays is emitted, the $T_{90}$. Long GRBs are those with $T_{90}>2$ s and short GRBs the sources with $T_{90}<2$ s Mazets et al. (1981); Klebesadel (1992); Dezalay et al. (1992); Kouveliotou et al. (1993); Tavani (1998).

In the case of short bursts, rapid consensus was reached in the scientific community that they could be the product of mergers of NS-NS and/or NS-BH binaries (see e.g., the pioneering works Goodman (1986); Paczynski (1986); Eichler et al. (1989); Narayan et al. (1991)). We shall return on this issue below by entering into the description of their properties and also to introduce additional mergers of compact-star object binaries leading to short bursts.

For long bursts, possibly the most compelling evidence of the necessity of a binary progenitor comes from the systematic and spectroscopic analysis of the GRBs associated with SNe, the so-called GRB-SNe, started with the pioneering discovery of the spatial and temporal concomitance of GRB 980425 Pian et al. (2000) and SN 1998bw Galama et al. (1998). Soon after, many associations of other nearby GRBs with type Ib/c SNe were evidenced (see e.g., Della Valle (2011); Cano et al. (2017)).

There are models in the literature attempting an explanation of both the SN and the GRB within the same astrophysical system. For instance, GRBs have been assumed to originate from a violent SN from the collapse of a massive and fast rotating star, a “collapsar” Woosley and Bloom (2006). A very high rotating rate of the star is needed to produce a collimated, jet emission. This traditional picture adopts for the GRB dynamics the “fireball” model based on the existence of a single ultrarelativistic collimated jet Blandford and McKee (1976); Shemi and Piran (1990); Meszaros et al. (1993); Piran et al. (1993); Mao and Yi (1994). There is a vast literature devoted to this “traditional" approach and we refer the reader to it for additional details (see, e.g., Piran (1999, 2004); Mészáros (2002, 2006); Berger (2014); Kumar and Zhang (2015), and references therein).

Nevertheless, it is worth to mention here some of the most important drawbacks of the aforementioned“traditional” approach and which has motivated the introduction of an alternative model, based on a binary progenitor, for the explanation of long GRBs:

• •

SNe Ic as the ones associated with GRBs lack hydrogen and helium in their spectra. It has been recognized that they most likely originate in helium stars, COcore, or Wolf-Rayet stars, that have lost their outermost layers (see e.g., Smith et al. (2011), and references therein). The pre-SN star, very likely, does not follow a single-star evolution but it belongs to a tight binary with a compact star companion (e.g., a NS). The compact star strips off the pre-SN star outermost layers via binary interactions such as mass-transfer and tidal effects (see e.g., Nomoto and Hashimoto (1988); Iwamoto et al. (1994); Fryer et al. (2007); Yoon et al. (2010); Smith et al. (2011)).

• •

Denoting the beaming angle by $\theta_{j}$, to an observed isotropic energy $E_{\rm iso}$ it would correspond to a reduced intrinsic source energy released $E_{s}=f_{b}E_{\rm iso}<E_{\rm iso}$, where $f_{b}=(1-\cos\theta_{j})\sim\theta^{2}_{j}/2<1$. Extremely small beaming factors $f_{b}\sim 1/500$ (i.e., $\theta_{j}\sim 1^{\circ}$) are inferred to reduce the observed energetics of $E_{\rm iso}\sim 10^{54}$ erg to the expected energy release by such a scenario $\sim$1051 erg Frail et al. (2001). However, the existence of such extremely narrow beaming angles have never been observationally corroborated Covino et al. (2006); Sato et al. (2007); Burrows et al. (2009).

• •

An additional drawback of this scenario is that it implies a dense and strong wind-like circumburst medium (CBM) in contrast with the one observed in most GRBs (see e.g., Izzo et al. (2012)). Indeed, the average CBM density inferred from GRB afterglows is of the order of 1 baryon per cubic centimeter Ruffini (2011). The baryonic matter component in the GRB process is represented by the so-called baryon load Ruffini et al. (2000). The GRB $e^{+}e^{-}$ plasma should engulf a limited amount of baryons in order to be able to expand at ultrarelativistic velocities with Lorentz factors $\Gamma\gtrsim 100$ as requested by the observed non-thermal component in the prompt Gamma-ray emission spectrum Shemi and Piran (1990); Piran et al. (1993); Meszaros et al. (1993). The amount of baryonic mass $M_{B}$ is thus limited by the prompt emission to a maximum value of the baryon-load parameter, $B=M_{B}c^{2}/E_{e^{+}e-}\lesssim 10^{-2}$, where $E_{e^{+}e-}$ is the total energy of the $e^{+}e-$ plasma Ruffini et al. (2000).

• •

GRBs and SNe have markedly different energetics. SNe emit energies in the range $10^{49}$–$10^{51}$ erg, while GRBs emit in the range $10^{49}$–$10^{54}$ erg. Thus, the origin of GRB energetics point to the gravitational collapse to a stellar-mass BH. The SN origin points to evolutionary stages of a massive star leading to a NS or to a complete disrupting explosion, but not to a BH. The direct formation of a BH in a core-collapse SN is currently ruled out by the observed masses of pre-SN progenitors, $\lesssim$18 $M_{\odot}$ Smartt (2015). It is theoretically known that massive stars with such a relatively low mass do not lead to a direct collapse to a BH (see Smartt (2009, 2015) for details) .

• •

It was recently shown in Ruffini et al. (2018) that the observed thermal emission in the X-ray flares present in the early (rest-frame time $t\sim 10^{2}$ s) afterglow implies an emitter of size $\sim$1012 cm expanding at mildly-relativistic velocity, e.g., $\Gamma\lesssim 4$. This is clearly in contrast with the “collapsar-fireball” scenario in which there is an ultrarelativistic emitter (the jet) with $\Gamma\sim 10^{2}$–$10^{3}$ extending from the prompt emission all the way to the afterglow.

Therefore, it seems most unlikely that the GRB and the SN can originate from the same single-star progenitor. Following this order of ideas, it was introduced for the explanation of the spatial and temporal coincidence of the two phenomena the concept of induced gravitational collapse (IGC) Ruffini et al. (2001, 2008). Two scenarios for the GRB-SN connection have been addressed: Ruffini et al. Ruffini et al. (2001) considered that the GRB was the trigger of the SN. However, for this scenario to happen it was shown that the companion star had to be in a very fine-tuned phase of its stellar evolution Ruffini et al. (2001). Ruffini et al. Ruffini et al. (2008) proposed an alternative scenario in a compact binary: the explosion of a Ib/c SN triggering an accretion process onto a NS companion. The NS, reaching the critical mass value, gravitationally collapses leading to the formation of a BH. The formation of the BH consequently leads to the emission of the GRB. Much more about this binary scenario has been discovered since its initial proposal; its theoretical studies and the search for its observational verification have led to the formulation of a much rich phenomenology which will be the main subject of this article.

Therefore, both short and long GRBs appear to be produced by binary systems, well in line with the expectation that most massive stars belong to binary systems (see, e.g., Sana et al. (2012); Smith (2014), and references therein). The increasing amount and quality of the multiwavelength data of GRBs have revealed the richness of the GRB phenomenon which, in a few seconds, spans different regimes from X-ray precursors to the Gamma-rays of the prompt emission, to the optical and X-rays of the early and late afterglow, to the optical emission of the associated SNe and, last but not least, the presence or absence of high-energy GeV emission. This, in addition to the multiyear effort of reaching a comprehensive theoretical interpretation of such regimes, have lead to the conclusion that GRBs separate into subclasses, each with specific energy release, spectra, duration, among other properties and, indeed, all with binary progenitors Ruffini et al. (2016); Rueda et al. (2017); Ruffini et al. (2018a, b); Wang et al. (2019).

1.2 GRB Subclasses

Up to 2017 we had introduced seven GRB subclasses summarized in Table 1. In addition, we have recently introduced in Rueda et al. (2018, 2019) the possibility of a further GRB subclass produced by WD-WD mergers. We now give a brief description of all the GRB subclasses identified. In Wang et al. (2019) we have renominated the GRB subclasses introduced in Ruffini et al. (2016) and in Rueda et al. (2018, 2019), and inserted them into two groups: binary-driven hypernovae (BdHNe) and compact-object binary mergers. Below we report both the old and the new names to facilitate the reader when consulting our works prior to Wang et al. (2019).

1. i.

X-ray flashes (XRFs). These systems have COcore-NS binary progenitors in which the NS companion does not reach the critical mass for gravitational collapse Becerra et al. (2016, 2015). In the SN explosion, the binary might or might not be disrupted depending on the mass loss and/or the kick imparted Postnov and Yungelson (2014). Thus XRFs lead either to two NSs ejected by the disruption, or to binaries composed of a newly-formed $\sim$1.4–1.5 $M_{\odot}$ NS (hereafter $\nu$NS) born at the center of the SN, and a massive NS (MNS) which accreted matter from the SN ejecta. Some observational properties are: Gamma-ray isotropic energy $E_{\rm iso}\lesssim 10^{52}$ erg, rest-frame spectral peak energy $E_{p,i}\lesssim 200$ keV and a local observed rate of $\rho_{\rm XRF}=100^{+45}_{-34}$ Gpc*-3* yr*-1* Ruffini et al. (2016). We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019), this class has been divided into BdHN type II, the sources with $10^{50}\lesssim E_{\rm iso}\lesssim 10^{52}$ erg, and BdHN type III, the sources with $10^{48}\lesssim E_{\rm iso}\lesssim 10^{50}$ erg. 2. ii.

Binary-driven hypernovae (BdHNe). Originate in compact COcore-NS binaries where the accretion onto the NS becomes high enough to bring it to the point of gravitational collapse, hence forming a BH. We showed that most of these binaries survive to the SN explosion owing to the short orbital periods ($P\sim 5$ min) for which the mass loss cannot be considered as instantaneous, allowing the binary to keep bound even if more than half of the total binary mass is lost Fryer et al. (2015). Therefore, BdHNe produce $\nu$NS-BH binaries. Some observational properties are: $E_{\rm iso}\gtrsim 10^{52}$ erg, $E_{p,i}\gtrsim 200$ keV and a local observed rate of $\rho_{\rm BdHN}=0.77^{+0.09}_{-0.08}$ Gpc*-3* yr*-1* Ruffini et al. (2016). We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019) this class has been renominated as BdHN type I. 3. iii.

BH-SN. These systems originate in COcore (or Helium or Wolf-Rayet star)-BH binaries, hence the hypercritical accretion of the SN explosion of the COcore occurs onto a BH previously formed in the evolution path of the binary. They might be the late evolutionary stages of X-ray binaries such as Cyg X-1 Giacconi and Ruffini (1978); Belczynski et al. (2011), or microquasars Mirabel and Rodríguez (1998). Alternatively, they can form following the evolutionary scenario XI in Fryer et al. (1999). If the binary survives to the SN explosion BH-SNe produce $\nu$NS-BH, or BH-BH binaries when the central remnant of the SN explosion collapses directly to a BH (see, although, Smartt (2009, 2015)). Some observational properties are: $E_{\rm iso}\gtrsim 10^{54}$ erg, $E_{p,i}\gtrsim 2$ MeV and an upper limit to their rate is $\rho_{\rm BH-SN}\lesssim\rho_{\rm BdHN}=0.77^{+0.09}_{-0.08}$ Gpc*-3* yr*-1*, namely the estimated observed rate of BdHNe type I which by definition covers systems with the above $E_{\rm iso}$ and $E_{p,i}$ range Ruffini et al. (2016). We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019) this class has been renominated as BdHN type IV.

We proceed with the short bursts which are amply thought to originate from compact-object binary mergers (BMs). First, we discuss the traditionally proposed BMs namely NS-NS and/or NS-BH mergers Goodman (1986); Paczynski (1986); Eichler et al. (1989); Narayan et al. (1991); Meszaros and Rees (1997); Rosswog et al. (2003); Lee et al. (2004); Berger (2014). These BMs can be separated into three subclasses Fryer et al. (2015); Ruffini et al. (2015, 2016):

1. iv.

Short Gamma-ray flashes (S-GRFs). They are produced by NS-NS mergers leading to a MNS, namely when the merged core does not reach the critical mass of a NS. Some observational properties are: $E_{\rm iso}\lesssim 10^{52}$ erg, $E_{p,i}\lesssim 2$ MeV and a local observed rate of $\rho_{\rm S-GRF}=3.6^{+1.4}_{-1.0}$ Gpc*-3* yr*-1* Ruffini et al. (2016). We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019) this class has been renominated as BM type I. 2. v.

Authentic short GRBs (S-GRBs). They are produced by NS-NS mergers leading to a BH, namely when the merged core reaches the critical mass of a NS, hence it forms a BH as a central remnant Ruffini et al. (2016, 2015); Muccino et al. (2013). Some observational properties are: $E_{\rm iso}\gtrsim 10^{52}$ erg, $E_{p,i}\gtrsim 2$ MeV and a local observed rate of $\rho_{\rm S-GRB}=\left(1.9^{+1.8}_{-1.1}\right)\times 10^{-3}$ Gpc*-3* yr*-1* Ruffini et al. (2016). We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019) this class has been renominated as BM type II. 3. vi.

Ultra-short GRBs (U-GRBs). This is a theoretical GRB subclass subjected for observational verification. U-GRBs are expected to be produced by $\nu$NS-BH mergers whose binary progenitors can be the outcome of BdHNe type I (see II above) or of BdHNe type IV (BH-SN; see III above). The following observational properties are expected: $E_{\rm iso}\gtrsim 10^{52}$ erg, $E_{p,i}\gtrsim 2$ MeV and a local observed rate similar to the one of BdHNe type I since we have shown that most of them are expected to remain bound Fryer et al. (2015), i.e., $\rho_{\rm U-GRB}\approx\rho_{\rm BdHN}=0.77^{+0.09}_{-0.08}$ Gpc*-3* yr*-1* Ruffini et al. (2016). We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019) this class has been renominated as BM type V.

Besides the existence of the above three subclasses of long bursts and three subclasses of short bursts in which the presence of NSs plays a fundamental role, there are two subclasses of bursts in which there is at least a WD component.

1. vii.

Gamma-ray flashes (GRFs). These sources show an extended and softer emission, i.e., they have hybrid properties between long and short bursts and have no associated SNe Della Valle et al. (2006). It has been proposed that they are produced by NS-WD mergers Ruffini et al. (2016). These binaries are expected to be very numerous Cadelano et al. (2015) and a variety of evolutionary scenarios for their formation have been proposed Bildsten and Cutler (1992); Fryer et al. (1999); Tauris et al. (2000); Lazarus et al. (2014). GRFs form a MNS and not a BH Ruffini et al. (2016). Some observational properties are: $10^{51}\lesssim E_{\rm iso}\lesssim 10^{52}$ erg, $0.2\lesssim E_{p,i}\lesssim 2$ MeV and a local observed rate of $\rho_{\rm GRF}=1.02^{+0.71}_{-0.46}$ Gpc*-3* yr*-1* Ruffini et al. (2016). It is worth noting that this rate is low with respect to the one expected from the current number of known NS-WD in the Galaxy Cadelano et al. (2015). From the GRB observations only one NS-WD merger has been identified (GRB 060614 Caito et al. (2009)). This implies that most NS-WD mergers are probably under the threshold of current X and Gamma-ray instruments. We refer the reader to Table 1 and Ruffini et al. (2016, 2018a) for further details on this class. In Wang et al. (2019) this class has been renominated as BM type III. 2. viii.

Fallback kilonovae (FB-KNe). This is a recently introduced GRB subclass having as progenitors WD-WD mergers Rueda et al. (2018, 2019). The WD-WD mergers of interest are those that do not produce type Ia SNe but that lead to a massive ($M\sim 1~{}M_{\odot}$), fast rotating ($P\sim 1$–$10$ s), highly-magnetized ($B\sim 10^{9}$–$10^{10}$ G) WD. Some observational properties are: $E_{\rm iso}\lesssim 10^{51}$ erg, $E_{p,i}\lesssim 2$ MeV and a local observed rate $\rho_{\rm FB-KN}=(3.7$–$6.7)\times 10^{5}$ Gpc*-3* yr*-1* Rueda et al. (2018, 2019); Maoz and Hallakoun (2017); Maoz et al. (2018). The coined name FB-KN is due to the fact that they are expected to produce an infrared-optical transient by the cooling of the ejecta expelled in the dynamical phase of the merger and heated up by fallback accretion onto the newly-formed massive WD.

The density rates for all GRB subclasses have been estimated assuming no beaming Ruffini et al. (2016, 2018a); Rueda et al. (2018, 2019). The GRB density rates have been analyzed in Ruffini et al. (2016) following the method suggested in Sun et al. (2015).

1.3 The Specific Case of BdHNe

We review in this article the specific case of BdHNe type I and II. As we have mentioned, the progenitor system is an exploding COcore as a type Ic SN in presence of a NS companion Ruffini et al. (2016); Wang et al. (2019). Figure 1 shows a comprehensive summary of the binary path leading to this variety of compact binaries that are progenitors of the above subclasses of long GRBs and that, at the same time, have an intimate connection with the short GRBs.

We emphasize on the theoretical framework concerning the COcore-NS binaries which have been extensively studied by our group in a series of publications Rueda and Ruffini (2012); Izzo et al. (2012); Fryer et al. (2014, 2015); Becerra et al. (2015, 2016). The COcore explodes as SN producing an accretion process onto the NS. For sufficiently compact binaries, e.g., orbital periods of the order of few minutes, the accretion is highly super-Eddington (hypercritical) leading to the possibility of the IGC of the NS once it reaches the critical mass, and forms a BH (see Figure 2).

If the binary is not disrupted by the explosion, BdHNe produces new binaries composed of a new NS ($\nu$NS) formed at the center of the SN, and a more massive NS or a BH companion (see Figure 1).

In the case of BH formation, the rotation of the BH together with the presence of the magnetic field inherited from the NS and the surrounding matter conform to what we have called the inner engine of the high-energy emission Ruffini et al. (2018a, b, 2019a, 2019b). The electromagnetic field of the engine is mathematically described by the Wald solution Wald (1974). The above ingredients induce an electric field around the BH which under the BdHN conditions is initially overcritical, creating electron-positron ($e^{+}e^{-}$) pair plasma which self-accelerates to ultrarelativistic velocities and whose transparency explains to the GRB prompt emission in Gamma-rays. The electric field is also able to accelerate protons which along the rotation axis lead to ultra high-energy cosmic rays (UHECRs) of up to $10^{21}$ eV. In the other directions the acceleration process lead to proton-synchrotron radiation which explains the GeV emission Ruffini et al. (2018a, b). The interaction/feedback of the GRB into the SN makes it become the hypernova (HN) Ruffini et al. (2018); Becerra et al. (2019) observed in the optical, powered by nickel decay, a few days after the GRB trigger. The SN shock breakout and the hypercritical accretion can be observed as X-ray precursors Becerra et al. (2016). The $e^{+}e^{-}$ feedback onto the SN ejecta also produces gamma- and X-ray flares observed in the early afterglow Ruffini et al. (2018). The synchrotron emission by relativistic electrons from the $\nu$NS in the expanding magnetized HN ejecta and the $\nu$NS pulsar emission explain the early and late X-ray afterglow Ruffini et al. (2018).

The article is organized as follows. In Section 2 we summarized following a chronological order the (1D, 2D and 3D) numerical simulations of BdHNe up to the year 2016, mentioning their salient features. A detailed explanation of the main ingredients of the calculations (equations of motion, accretion modeling, NS evolution equations, critical mass, accretion-zone hydrodynamics, neutrino emission and accretion energy release) can be found in Section 3. The most recent 3D smoothed-particle-hydrodynamics (SPH) numerical simulations of 2018 are presented in Section 4. Section 5 is devoted to the consequences on these simulations on the analysis and interpretation of the GRB multiwavelength data. In Section 6 we present an analysis of the binary gravitational binding of BdHNe progenitors, so it is shown that most BdHNe type I are expected to be NS-BH binaries. The cosmological evolutionary scenario leading to the formation of BdHN, their occurrence rate and connection with short GRBs is presented in Section 7.

We show in Table 2 a summary of acronyms used in this work.

2 A Chronological Summary of the IGC Simulations: 2012–2016

2.1 First Analytic Estimates

The IGC scenario was formulated in 2012 Rueda and Ruffini (2012) presenting a comprehensive astrophysical picture supporting this idea as well as a possible evolutionary scenario leading to the progenitor COcore-NS binaries. It was also there presented an analytic formula for the accretion rate onto the NS companion on the basis of the following simplified assumptions: (1) a uniform density profile of the pre-SN COcore; (2) the ejecta was evolved following an homologous expansion; (3) the mass of the NS (assumed to be initially $1.4~{}M_{\odot}$) and the COcore (in the range $4$–$8~{}M_{\odot}$) were assumed nearly constant. So, it was shown that the accretion rate onto the NS is highly super-Eddington, namely it is hypercritical, reaching values of up to $0.1~{}M_{\odot}$ s*-1* for compact binaries with orbital periods of the order of a few minutes. This estimate implied that the hypercritical accretion could induce the gravitational collapse of the NS which, in a few seconds, would reach the critical mass with consequence formation of a BH. A first test of this IGC first model in real data was soon presented in the case of GRB 090618 Izzo et al. (2012).

2.2 First Numerical Simulations: 1D Approximation

The first numerical simulations were implemented in 2014 in Fryer et al. (2014) via a 1D code including (see Figure 3): (1) the modeling of the SN via the 1D core-collapse SN code of Los Alamos Fryer et al. (1999); (2) the microphysics experienced by the inflow within the accretion region including the neutrino ($\nu$) emission and hydrodynamics processes such as shock formation; (3) with the above it was followed by the evolution of the material reaching the Bondi-Hoyle capture region and the subsequent in-fall up to the NS surface. Hypercritical accretion rates in the range $10^{-3}$–$10^{-1}~{}M_{\odot}$ s*-1* were inferred, confirming the first analytic estimates and the IGC of the NS companion for binary component masses similar to the previous ones and for orbital periods of the order of $5$ min.

The above simulations were relevant in determining that the fate of the system is mainly determined by the binary period ($P$); the SN ejecta velocity ($v_{\rm ej}$) and the NS initial mass. $P$ and $v_{\rm ej}$ enter explicitly in the Bondi-Hoyle accretion rate formula through the capture radius expression, and implicitly via the ejecta density since they influence the decompression state of the SN material at the NS position.

2.3 2D Simulations including Angular Momentum Transfer

Soon after, in 2015, we implemented in Becerra et al. (2015) a series of improvements to the above calculations by relaxing some of the aforementioned assumptions (see Figure 4). We adopted for the ejecta a density profile following a power-law with the radial distance and evolved it with an homologous expansion. The angular momentum transport, not included in the previous estimates, was included. With this addition it was possible to estimate the spin-up of the NS companion by the transfer of angular momentum from the in-falling matter which was shown to circularize around the NS before being accreted. General relativistic effects were also introduced, when calculating the evolution of the structure parameters (mass, radius, spin, etc) of the accreting NS, in the NS gravitational binding energy, and in the angular momentum transfer by the circularized particles being accreted from the innermost circular orbit.

One of the most important results of Becerra et al. (2015) was that, taking into account that the longer the orbital $P$ the lower the accretion rate, it was there computed the maximum orbital period ($P_{\rm max}$) for which the NS reaches the critical mass for gravitational collapse, so for BH formation. The dependence of $P_{\rm max}$ on the initial mass of the NS was also there explored. The orbital period $P_{\rm max}$ was then presented as the separatrix of two families of long GRBs associated with these binaries: at the time we called them Family-1, the systems in which the NS does not reach the critical mass, and Family-2 the ones in which it reaches the critical mass and forms a BH. It can be seen that the Family-1 and Family-2 long GRBs evolve subsequently into the concepts of XRFs and BdHNe, respectively.

2.4 First 3D Simulations

A great step toward the most recent simulations was achieved in 2016 in Becerra et al. (2016) where an SPH-like simulation was implemented in which the SN ejecta was emulated by “point-like” particles. The mass and number of the particles populating each layer were assigned, for self-consistency, according to the power-law density profile. The initial velocity of the particles of each layer was set, in agreement with the chosen power-law density profile, following a radial velocity distribution; i.e., $v\propto r$.

The evolution of the SN particles was followed by Newtonian equations of motion in the gravitational field of the NS companion, also taking into account the orbital motion which was included under the assumption that the NS performs a circular orbit around the COcore center that acts as the common center-of-mass, namely assuming that the mass of the pre-SN core is much larger than the NS mass.

The accretion rate onto the NS was computed, as in Becerra et al. (2015), using the Bondi-Hoyle accretion formula and, every particle reaching the Bondi-Hoyle surface, was removed from the system. The maximum orbital period $P_{\textrm{max}}$ in which the NS collapses by accretion could be further explored including the dependence on the mass of the pre-SN CO${}_{\textrm{core}}$, in addition to the dependence on the NS mass.

A detailed study of the hydrodynamics and the neutrino emission in the accretion region on top the NS surface was performed. Concerning the neutrino emission, several $\nu$ and antineutrino ($\bar{\nu}$) production processes were considered and showed that electron-positron annihilation ($e^{+}e^{-}\to\nu\bar{\nu}$) overcomes by orders of magnitude any other mechanism of neutrino emission in the range of accretion rates $10^{-8}$–$10^{-2}~{}M_{\odot}$ s*-1*, relevant for XRFs and BdHNe. The neutrino luminosity can reach values of up to $10^{52}$ erg s*-1* and the neutrino mean energy of $20$ MeV for the above upper value of the accretion rate. For the reader interested in the neutrino emission, we refer to Becerra et al. (2018) for a detailed analysis of the neutrino production in XRFs and BdHNe including flavor oscillations experienced by the neutrinos before abandoning the system.

Concerning the hydrodynamics, the evolution of the temperature and density of outflows occurring during the accretion process owing to convective instabilities was estimated. It was there shown the interesting result that the temperature of this outflow and its evolution can explain the early (i.e., precursors) X-ray emission that has been observed in some BdHNe and in XRFs, exemplified there analyzing the early X-ray emission observed in GRB 090618, a BdHN I, and in GRB 060218, a BdHN II (an XRF).

A most important result of these simulations was the possibility of having a first glance of the morphology acquired by the SN ejecta: the matter density, initially spherically symmetric, becomes highly asymmetric due to the accretion process and the action of the gravitational field of the NS companion (see Figure 5).

3 The Hypercritical Accretion Process

We now give details of the accretion process within the IGC scenario following Fryer et al. (2014); Becerra et al. (2015); Fryer et al. (2015); Becerra et al. (2016). There are two main physical conditions for which hypercritical (i.e., highly super-Eddington) accretion onto the NS occurs in XRFs and BdHNe. The first is that the photons are trapped within the inflowing material and the second is that the shocked atmosphere on top of the NS becomes sufficiently hot ($T\sim 10^{10}$ K) and dense ($\rho\gtrsim 10^{6}$ g cm*-3*) to produce a very efficient neutrino-antineutrino ($\nu\bar{\nu}$) cooling emission. In this way the neutrinos become mainly responsible for releasing the energy gained by accretion, allowing hypercritical accretion to continue.

3.1 Accretion Rate and NS Evolution

The first numerical simulations of the IGC were performed in Fryer et al. (2014), including: (1) realistic SN explosions of the COcore; (2) the hydrodynamics within the accretion region; (3) the simulated evolution of the SN ejecta up to their accretion onto the NS. Becerra et al. Becerra et al. (2015) then estimated the amount of angular momentum carried by the SN ejecta and how much is transferred to the NS companion by accretion. They showed that the SN ejecta can circularize for a short time and form a disc-like structure surrounding the NS before being accreted. The evolution of the NS central density and rotation angular velocity (the NS is spun up by accretion) was computed from full numerical solutions of the axisymmetric Einstein equations. The unstable limits of the NS are set by the mass-shedding (or Keplerian) limit and the critical point of gravitational collapse given by the secular axisymmetric instability, (see, e.g., Becerra et al. (2015) for details).

The accretion rate of the SN ejecta onto the NS is given by:

[TABLE]

where $G$ is the gravitational constant, $\rho_{\rm ej}$ and $c_{{\rm s,ej}}$ are the density and sound speed of the ejecta, $R_{\rm cap}$ and $M_{\rm NS}$ are the NS gravitational capture radius (Bondi-Hoyle radius) and gravitational mass, and $v_{{\rm rel}}$ the ejecta velocity relative to the NS: $\vec{v}_{{\rm rel}}=\vec{v}_{{\rm orb}}-\vec{v}_{{\rm ej}}$; $|\vec{v}_{\rm orb}|=\sqrt{G(M_{\rm core}+M_{\rm NS})/a}$, and $\vec{v}_{\rm ej}$ is the velocity of the supernova ejecta (see Figure 2).

Numerical simulations of the SN explosions suggest the adopted homologous expansion of the SN, i.e., $v_{\rm ej}(r,t)=nr/t$, where $r$ is the position of each layer from the SN center and $n$ is the expansion parameter. The density evolves as

[TABLE]

where $M_{\rm env}(t)$ the mass of the COcore envelope, $R_{\rm star}(t)$ is the radius of the outermost layer, and $\rho_{\rm ej}^{0}$ is the pre-SN COcore density profile; $\rho_{\rm ej}(r,t_{0})=\rho_{\rm core}(R_{\rm core}/r)^{m}$, where $\rho_{\rm core}$, $R_{\rm core}$ and $m$ are the profile parameters obtained from numerical simulations. Typical parameters of the COcore mass are (3.5–9.5) $M_{\odot}$ corresponding to $(15$–$30)~{}M_{\odot}$ zero-age-main-sequence (ZAMS) progenitors (see Fryer et al. (2014); Becerra et al. (2015) for details). The binary period is limited from below by the request of having no Roche lobe overflow by the COcore before the SN explosion Fryer et al. (2014). For instance, for a COcore of 9.5 $M_{\odot}$ forming a binary system with a 2 $M_{\odot}$ NS, the minimum orbital period allowed by this condition is $P_{\rm min}\approx 5$ min. For these typical binary and pre-SN parameters, Equation (1) gives accretion rates $10^{-4}$–$10^{-2}M_{\odot}$ s*-1*.

We adopt an initially non-rotating NS companion so its exterior spacetime at time $t=0$ is described by the Schwarzschild metric. The SN ejecta approach the NS with specific angular momentum, $l_{\rm acc}=\dot{L}_{\rm cap}/\dot{M}_{B}$, circularizing at a radius $r_{\rm circ}\geq r_{\rm lco}$ if $l_{\rm acc}\geq l_{\rm lso}$ with $r_{\rm lco}$ the radius of the last circular orbit (LCO). For a non-rotating NS $r_{\rm lco}=6GM_{\rm NS}/c^{2}$ and $l_{\rm lco}=2\sqrt{3}GM_{\rm NS}/c$. For typical parameters, $r_{\rm circ}/r_{\rm lco}\sim 10$–$10^{3}$.

The accretion onto the NS proceeds from the radius $r_{\rm in}$. The NS mass and angular angular momentum evolve as Becerra et al. (2015); Cipolletta et al. (2017):

[TABLE]

where $M_{b}$ is the NS baryonic mass, $l(r_{\rm in})$ is the specific angular momentum of the accreted material at $r_{\rm in}$, which corresponds to the angular momentum of the LCO, and $\xi\leq 1$ is a parameter that measures the efficiency of angular momentum transfer. In this picture we have $\dot{M}_{b}=\dot{M}_{B}$.

For the integration of Equations (1) and (3) we have to supply the values of the two partial derivatives in Equation (3). They are obtained from the relation of the NS gravitational mass, $M_{\rm NS}$, with $M_{b}$ and $J_{\rm NS}$, namely from the knowledge of the NS binding energy. For this we use the general relativistic calculations of rotating NSs presented in Cipolletta et al. (2015). They show that, independent on the nuclear EOS, the following analytical formula represents the numerical results with sufficient accuracy (error $<2\%$):

[TABLE]

where $j_{\rm NS}\equiv cJ_{\rm NS}/(GM_{\odot}^{2})$.

In the accretion process, the NS gains angular momentum and therefore spins up. To evaluate the amount of angular momentum transferred to the NS at any time we include the dependence of the LCO specific angular momentum as a function of $M_{\rm NS}$ and $J_{\rm NS}$. For corotating orbits, the following relation is valid for the NL3, TM1 and GM1 EOS Cipolletta et al. (2017); Becerra et al. (2015):

[TABLE]

The NS continues the accretion until it reaches an instability limit or up to when all the SN ejecta overcomes the NS Bondi-Hoyle region. We take into account the two main instability limits for rotating NSs: the mass-shedding or Keplerian limit and the secular axisymmetric instability limit. The latter defines critical NS mass. For the aforementioned nuclear EOS, the critical mass can be approximately written as Cipolletta et al. (2015):

[TABLE]

where $k$ and $p$ are EOS-dependent parameters (see Table 3). These formulas fit the numerical results with a maximum error of 0.45%.

Additional details and improvements of the hypercritical accretion process leading to XRFs and BdHNe were presented in Becerra et al. (2016). Specifically:

The density profile included finite size/thickness effects and additional COcore progenitors, leading to different SN ejecta masses being considered. 2. 2.

In Becerra et al. (2015) the maximum orbital period, $P_{\rm max}$, over which the accretion onto NS companion is not sufficient to bring it to the critical mass, was inferred. Thus, binaries with $P>P_{\rm max}$ lead to XRFs while the ones with $P\lesssim P_{\rm max}$ lead to BdHNe. Becerra et al. Becerra et al. (2016) extended the determination of $P_{\rm max}$ for all the possible initial values of the NS mass. They also examined the outcomes for different values of the angular momentum transfer efficiency parameter. 3. 3.

The expected luminosity during the process of hypercritical accretion for a wide range of binary periods covering both XRFs and BdHNe was estimated. 4. 4.

It was shown that the presence of the NS companion originates asymmetries in the SN ejecta (see, e.g., Figure 6 in Becerra et al. (2016)). The signatures of such asymmetries in the X-ray emission was there shown in the specific example of XRF 060218.

3.2 Hydrodynamics in the Accretion Region

The accretion rate onto the NS can be as high as $\sim$10*-2*–$10^{-1}~{}M_{\odot}$ s*-1*. For such accretion rates:

The magnetic pressure is much smaller than the random pressure of the infalling material, therefore the magnetic-field effects on the accretion process are negligible Fryer et al. (1996); Rueda and Ruffini (2012). 2. 2.

The photons are trapped within the infalling matter, hence the Eddington limit does not apply and hypercritical accretion occurs. The trapping radius is defined as Chevalier (1989): $r_{\rm trapping}={\rm min}\{\dot{M}_{B}\kappa/(4\pi c),R_{\rm cap}\}$, where $\kappa$ is the opacity. Fryer et al. (2014) estimated a Rosseland mean opacity of $\approx$5 $\times$ 103 cm2 g*-1* for the COcores. This, together with our typical accretion rates, lead to $\dot{M}_{B}\kappa/(4\pi c)\sim 10^{13}$–$10^{19}$ cm. This radius is much bigger than the Bondi-Hoyle radius. 3. 3.

The above condition, and the temperature-density values reached on top of the NS surface, lead to an efficient neutrino cooling which radiates away the gain of gravitational energy of the infalling material Zel’dovich et al. (1972); Ruffini and Wilson (1973); Fryer et al. (1996); Rueda and Ruffini (2012); Fryer et al. (2014).

The accretion shock moves outward as the material piles onto the NS. Since the post-shock entropy is inversely proportional to the shock radius position, the NS atmosphere is unstable with respect to Rayleigh-Taylor convection at the beginning of the accretion process. Such instabilities might drive high-velocity outflows from the accreting NS Fryer et al. (2006); Fryer (2009). The entropy at the base of the atmosphere is Fryer et al. (1996):

[TABLE]

where $k_{B}$ is the Boltzmann constant. The material expands and cools down adiabatically, i.e., $T^{3}/\rho$ = constant. In the case of a spherically symmetric expansion, $\rho\propto 1/r^{3}$ and $k_{B}T_{\rm bubble}=195\,S_{\rm bubble}^{-1}\left(10^{6}\,{\rm cm}/r\right)$ MeV. In the more likely case that the material expand laterally we have Fryer (2009): $\rho\propto 1/r^{2}$, i.e., $T_{\rm bubble}=T_{0}(S_{\rm bubble})\left(r_{0}/r\right)^{2/3}$, where $T_{0}(S_{\rm bubble})$ is obtained from the above equation at $r=r_{0}\approx R_{\rm NS}$. This implies a bolometric blackbody flux at the source from the rising bubbles:

[TABLE]

The above thermal emission has been shown Fryer et al. (2014) to be a plausible explanation of the early X-ray (precursor) emission observed in some GRBs. The X-ray precursor observed in GRB 090618 Izzo et al. (2012, 2012) is explained adopting an accretion rate of $10^{-2}~{}M_{\odot}$ s*-1*, the bubble temperature drops from $50$ keV to $15$ keV while expanding from $r\approx 10^{9}$ cm to $6\times 10^{9}$ cm (see Figure 6). More recently, the X-ray precursor has been observed in GRB 180728A and it is well explained by a bubble of $\sim$7 keV at \sim$$10^{10} cm and an accretion rate of $10^{-3}~{}M_{\odot}$ s*-1* (see Wang et al. (2019) for details).

3.3 Neutrino Emission and Effective Accretion Rate

For the accretion rate conditions characteristic of our models at peak $\sim$10*-4*–$10^{-2}~{}M_{\odot}$ s*-1*, pair annihilation dominates the neutrino emission and electron neutrinos remove the bulk of the energy Becerra et al. (2016). The $e^{+}e^{-}$ pairs producing the neutrinos are thermalized at the matter temperature. This temperature is approximately given by:

[TABLE]

where $P_{\rm shock}$ is the pressure of the shock developed on the accretion zone above the NS surface, $\dot{M}_{\rm acc}$ is the accretion rate, $v_{\rm acc}$ is the velocity of the infalling material, $\sigma$ is the Stefan-Boltzmann constant and $c$ the speed of light. It can be checked that, for the accretion rates of interest, the system develops temperatures and densities $T\gtrsim 10^{10}$ K and $\rho\gtrsim 10^{6}$ g cm*-3*; respectively (see Figure 7).

Under these conditions of density and temperature the neutrino emissivity of the $e^{+}e^{-}$ annihilation process can be estimated by the simple formula Yakovlev et al. (2001):

[TABLE]

where $k_{B}$ is the Boltzmann constant.

The accretion zone is characterized by a temperature gradient with a typical scale height $\Delta r_{\rm ER}=T/\nabla T\approx 0.7~{}R_{\rm NS}$. Owing to the aforementioned strong dependence of the neutrino emission on temperature, most of the neutrinos are emitted from a spherical shell around the NS of thickness

[TABLE]

Equations (9) and (10) imply the neutrino emissivity satisfies $\epsilon_{e^{-}\!e^{+}}\propto\dot{M}^{9/4}_{\rm acc}$ as we had anticipated. These conditions lead to the neutrinos to be efficient in balancing the gravitational potential energy gain, allowing the hypercritical accretion rates. The effective accretion onto the NS can be estimated as:

[TABLE]

where $\Delta M_{\nu}$, $L_{\nu}$ are, respectively, the mass and neutrino luminosity in the emission region, and $E_{g}=(1/2)GM_{\rm NS}\Delta M_{\nu}/(R_{\nu}+\Delta r_{\nu})$ is half the gravitational potential energy gained by the material falling from infinity to the $R_{\rm NS}+\Delta r_{\nu}$. The neutrino luminosity is

[TABLE]

with $\epsilon_{e^{-}\!e^{+}}$ being the neutrino emissivity in Equation (10). For $M_{\rm NS}=2~{}M_{\odot}$ and temperatures $1$–10 MeV, the Equations (12) and (13) result $\dot{M}_{\rm eff}\approx 10^{-10}$–$10^{-1}~{}M_{\odot}$ s*-1* and $L_{\nu}\approx 10^{48}$–$10^{57}$ MeV s*-1*.

Therefore, the neutrino emission can reach luminosities of up to $10^{57}$ MeV s*-1*, mean neutrino energies $20$–$30$ MeV, and neutrino densities $10^{31}$ cm*-3*. Along their path from the vicinity of the NS surface outward, such neutrinos experience flavor transformations dictated by the neutrino to electron density ratio. We have determined in Becerra et al. (2018) the neutrino and electron on the accretion zone and use them to compute the neutrino flavor evolution. For normal and inverted neutrino-mass hierarchies and within the two-flavor formalism ($\nu_{e}\nu_{x}$), we estimated the final electronic and non-electronic neutrino content after two oscillation processes: (1) neutrino collective effects due to neutrino self-interactions where the neutrino density dominates and, (2) the Mikheyev-Smirnov-Wolfenstein (MSW) effect, where the electron density dominates. We find that the final neutrino content is composed by $\sim$55% ($\sim$62%) of electronic neutrinos, i.e., $\nu_{e}+\bar{\nu}_{e}$, for the normal (inverted) neutrino-mass hierarchy (see Figure 8). This is a first step toward the characterization of a novel source of astrophysical MeV-neutrinos in addition to core-collapse SNe. We refer the reader to Becerra et al. (2018) for additional details of the flavor-oscillations as well as the final neutrino spectra after such a process.

3.4 Accretion Luminosity

The energy release in a time-interval $dt$, when an amount of mass $dM_{b}$ with angular momentum $l\dot{M}_{b}$ is accreted, is Becerra et al. (2016):

[TABLE]

This is the amount of gravitational energy gained by the matter by infalling to the NS surface that is not spent in NS gravitational binding energy. The total energy release in the time interval from $t$ to $t+dt$,

[TABLE]

is given by the NS binding energy difference between its initial and final state. The typical luminosity is $L_{\rm acc}\approx\Delta E_{\rm acc}/\Delta t_{\rm acc}$, where $\Delta t_{\rm acc}$ is the duration of the accretion process.

The value of $\Delta t_{\rm acc}$ is approximately given by the flow time of the slowest layers of the SN ejecta to the NS companion position. If we denote the velocity of these layers by $v_{\rm inner}$, we have $\Delta t_{\rm acc}\sim a/v_{\rm inner}$, where $a$ is the binary separation. For $a\sim 10^{11}$ cm and $v_{\rm inner}\sim 10^{8}$ cm s*-1*, $\Delta t_{\rm acc}\sim 10^{3}$ s. For shorter separations, e.g., $a\sim 10^{10}$ cm ($P\sim 5$ min), $\Delta t_{\rm acc}\sim 10^{2}$ s. For a binary with $P=5$ min, the NS accretes $\approx$1 $M_{\odot}$ in $\Delta t_{\rm acc}\approx 100$ s. From Equation (4) one obtains that the binding energy difference of a $2~{}M_{\odot}$ and a $3~{}M_{\odot}$ NS, is $\Delta E_{\rm acc}\approx 13/200(3^{2}-2^{2})~{}M_{\odot}c^{2}\approx 0.32~{}M_{\odot}c^{2}$. This leads to $L_{\rm acc}\approx 3\times 10^{-3}~{}M_{\odot}c^{2}\approx 0.1\,\dot{M_{b}}c^{2}$. The accretion power can be as high as $L_{\rm acc}\sim 0.1\dot{M_{b}}c^{2}\sim 10^{47}$–$10^{51}$ erg s*-1* for accretion rates in the range $\dot{M_{b}}\sim 10^{-6}$–$10^{-2}~{}M_{\odot}$ s*-1*.

4 New 3D SPH Simulations

We have recently presented in Becerra et al. (2019) new, 3D hydrodynamic simulations of the IGC scenario by adapting the SPH code developed at Los Alamos, SNSPH Fryer et al. (2006), which has been tested and applied in a variety of astrophysical situations Fryer and Warren (2002); Young et al. (2006); Diehl et al. (2008); Batta et al. (2017).

The time $t=0$ of the simulation is set as the time at which the SN shock breaks out the COcore external radius. We calculate the accretion rate both onto the NS companion and onto the $\nu$NS (via fallback), and calculate the evolution of other binary parameters such as the orbital separation, eccentricity, etc. Figure 9 shows an example of simulation for a binary system composed of a COcore of mass $\approx$6.85 $M_{\odot}$, the end stage of a ZAMS progenitor star of $M_{\rm zams}=25~{}M_{\odot}$, and a $2~{}M_{\odot}$ NS companion. The initial orbital period is $\approx$5 min.

The accretion rate onto both stars was estimated from the flux of SPH particles falling, per unit time, into the Bondi-Hoyle accretion region of the NS (see Figure 10). It is confirmed that the accretion onto the NS companion occurs from a disk-like structure formed by the particles that circularize before being accreted; see vortexes in the upper panel of Figure 9 and the disk structure is clearly seen in the lower panel.

Several binary parameters were explored thanks to the new code. We performed simulations changing the COcore mass, the NS companion mass, the orbital period, the SN explosion energy (so the SN kinetic energy or velocity). We also explored intrinsically asymmetric SN explosion. We checked if the $\nu$NS and/or the NS companion reach the mass-shedding (Keplerian) limit or the secular axisymmetric instability, i.e., the critical mass. The NS can also become just a more massive, fast rotating, stable NS when the accretion is moderate. All this was done for various NS nuclear equations of state (NL3, TM1 and GM1).

We followed the orbital evolution up to the instant when most of the ejecta has abandoned the system to determine if the system remains bound or becomes unbound by the explosion. We thus assessed the COcore-NS parameters leading to the formation of $\nu$NS-NS (from XRFs) or $\nu$NS-BH (from BdHNe) binaries. The first proof that BdHNe remain bound leading to $\nu$NS-BH binaries was presented in Fryer et al. (2015) (see next section).

5 Consequences on GRB Data Analysis and Interpretation

In a few seconds a BdHN shows different physical processes that lead to a specific sequence of observables at different times and at different wavelengths. Starting with the at-times-observed X-ray precursors, to the Gamma-ray prompt emission, to the GeV emission, to the early and late X-ray afterglow in which, respectively, are observed flares and a distinct power-law luminosity.

5.1 X-ray Precursor

X-ray precursors can comprise the presence of both the SN shock breakout as well as the hypercritical accretion onto the NS companion until it reaches the critical mass. These processes have been identified in Izzo et al. (2012); Fryer et al. (2014); Becerra et al. (2016); Wang et al. (2019).

The conversion of the SN shockwave kinetic energy (see Arnett (1996) for details on the SN physics) into electromagnetic energy imply that about $10^{50}$ erg can be emitted.

Once it reached the NS companion, the ejecta induced a hypercritical accretion onto the NS at a rate $\sim$10${}^{-3}M_{\odot}$ s*-1* for an assumed orbital separation of few $10^{10}$ cm. As we have recalled (see Figure 6 in Section 3), the accretion process triggers the expansion of thermal convective bubbles on top of the NS owing to the Rayleigh-Taylor instability Izzo et al. (2012); Fryer et al. (2014); Becerra et al. (2016); Wang et al. (2019).

It is of special interest to refer the reader to the results presented in Wang et al. (2019) on GRB 180728A, a BdHN II. It has been identified in the precursor of this GRB, for the first time, the presence of both the emergence of the SN shockwave as well as the hypercritical accretion process. From this the binary parameters have been extracted and further confirmed by the analysis of the prompt and the afterglow emission.

5.2 GRB Prompt Emission

A BdHN I leaves as a remnant a $\nu$NS-BH binary surrounded by the asymmetric SN ejecta (see Figure 5 and Becerra et al. (2016, 2019)). The asymmetric ejecta includes a “cavity” of $\sim$1011 cm of very low-density matter around the newborn BH. The hydrodynamics inside such a low-density cavity have been recently studied by numerical simulations in Ruffini et al. (2019b).

The asymmetric character acquired by the SN ejecta implies that the $e^{+}e-$ plasma, expanding from the BH site in all directions with equal initial conditions, experiences a different dynamics along different directions. The reason for this is that the $e^{+}e-$ plasma engulfs different amounts of baryonic mass (see Figure 11). This leads to observable signatures as a function of the viewing angle.

The newborn Kerr BH, surrounded by ejecta and immersed in a test magnetic field (likely the one left by the magnetized, collapsed NS), represents what we have called the inner engine of the high-energy emission Ruffini et al. (2018a, b, 2019a, 2019b). The rotating BH, of mass $M$ and angular momentum $J$, in the presence of the magnetic field $B_{0}$, induces an electromagnetic field described by the Wald solution Wald (1974).

The induced electric field at the BH horizon $r_{+}=M(1+\sqrt{1-\alpha^{2}})$ is Ruffini et al. (2018a, b)

[TABLE]

where $\alpha=J/M^{2}$ is the dimensionless angular momentum of the BH and $B_{c}=m_{e}^{2}c^{3}/(e\hbar)\approx 4.4\times 10^{13}$ G. This field acquires values over the critical one, $E_{c}=m_{e}^{2}c^{3}/(e\hbar)$ if the following conditions are verified:

[TABLE]

where the second condition comes from the constraint that a rotating BH must satisfy: $\alpha\leq 1$. The above huge value of the electric field (16) guarantees the production of the $e^{+}e^{-}$ pair plasma around the newborn BH via the quantum electrodynamics (QED) process of vacuum polarization Ruffini et al. (2010).

In the direction pointing from the COcore to the accreting NS outwards and lying on the orbital plane, the aforementioned cavity represents a region of low baryonic contamination Ruffini et al. (2019b, a). The $e^{+}e-$ plasma can then self-accelerate to Lorentz factors $\Gamma\sim 10^{2}$–$10^{3}$ reaching transparency and impacting on the CBM filaments as described in Preparata et al. (1998); Ruffini et al. (1999, 2000). At transparency, MeV-photons are emitted which are observed in the ultrarelativistic prompt emission. This picture has been successfully applied and verified on plenty of GRBs, e.g., GRBs 050904, 080319B, 090227, 090618 and 101023 Penacchioni et al. (2012); Izzo et al. (2012); Patricelli et al. (2012); Muccino et al. (2013).

5.3 Early X-ray Afterglow: Flares

It was recently addressed in Ruffini et al. (2018) the role of X-ray flares as a powerful tool to differentiate the BdHN model from the “collapsar-fireball” model Woosley (1993).

First, it is known that the GRB prompt emission shows Gamma-ray spikes occurring at $10^{15}$–$10^{17}$ cm from the source and have Lorentz factor $\Gamma\sim 10^{2}$–$10^{3}$.

Second, the thermal emission observed in the X-ray flares of the early (rest-frame time $t\sim 10^{2}$ s) afterglow of BdHNe, implies occurrence radii $\sim$1012 cm expanding at mildly-relativistic velocity, e.g., $\Gamma\lesssim 4$ Ruffini et al. (2018) (see below). The latter observational fact evidences that the X-ray afterglow is powered by a mildly-relativistic emitter. These model-independent observations contrast with the assumption of an ultrarelativistic expansion starting from the GRB prompt emission and extending to the afterglow. Such a “traditional” approach to GRBs has been adopted in a vast number of articles over decades as it is summarized in review articles (see, e.g., Piran (1999, 2004); Mészáros (2002, 2006); Berger (2014); Kumar and Zhang (2015).

In the other directions, the GRB $e^{+}e-$ plasma impacts the SN ejecta at approximately $10^{10}$ cm, evolves carrying a large amount of baryons reaching transparency at radii $10^{12}$ cm with a mildly $\Gamma\lesssim 4$. The theoretical description and the consequent numerical simulation have been addressed in Ruffini et al. (2018).

Such a mildly-relativistic photospheric emission is experimentally demonstrated by the thermal radiation observed in the early X-ray afterglow and in the X-ray flares Starling et al. (2012); Wang et al. (2018). For instance, in the early hundreds of seconds, GRB 090618 is found to have a velocity of $\beta\sim 0.8$ Ruffini et al. (2014); Page et al. (2011), GRB 081008 has a velocity $\beta\sim 0.9$ Ruffini et al. (2018), and GRB 130427A has a velocity of $\beta\sim 0.9$ as well Ruffini et al. (2015); Wang et al. (2015); Ruffini et al. (2018). We emphasize that the mildly-relativistic photo-sphere velocity is derived from the data in a model-independent way, summarising from Ruffini et al. (2018):

[TABLE]

The left-hand side is a function of velocity $\beta$, the right-hand side is only from observables, $D_{L}(z)$ is the luminosity distance and $z$ the cosmological redshift. From the observed thermal flux $F_{\mathrm{bb,obs}}$ and temperature $T_{\mathrm{obs}}$ in two times $t_{1}$ and $t_{2}$, the velocity $\beta$ is obtained. This model-independent equation has been derived in a fully relativistic way so it remains valid in the Newtonian non-relativistic regime.

An additional, and very important prediction of this scenario, is that the injection of energy and momentum from the GRB plasma into the ejecta transforms the SN into an HN (see Ruffini et al. (2018) for the specific case of GRB 151027A).

5.4 Late X-ray Afterglow

We have shown in Ruffini et al. (2018) that the synchrotron emission by relativistic electrons from the $\nu$NS, injected into the expanding magnetized HN ejecta, together with the $\nu$NS pulsar emission that extracts its rotational energy, power the X-ray afterglow. This includes the early part and the late power-law behavior. An exceptional by-product of this analysis is that it gives a glimpse on the $\nu$NS magnetic field strength and structure (dipole+quadrupole).

Based on the above model Ruffini et al. (2018), GRB 130427A (a BdHN I) and GRB 180728A (a BdHN II) have been analyzed in Wang et al. (2019). The explanation of the afterglow data of GRB 130427A led to an initial $1$ ms rotation period for the $\nu$NS. For GRB 180728A, a slower spin of $2.5$ ms was there obtained. A simple analysis showed how this result is in agreement with the BdHN I and II nature of these GRBs. First, we recall that compact binary systems have likely synchronized components with the orbital period. Second, we can infer the orbital period from the analysis of the X-ray precursor and the prompt emission (see Wang et al. (2019) for the procedure). Then, we can infer the COcore rotation period too. Finally, assuming angular momentum conservation in the core-collapse SN process, we can estimate the rotation period of the $\nu$NS formed at the SN center. This method led to a binary separation remarkably in agreement with the one inferred from the precursor and the prompt emission, demonstrating the self-consistency of this scenario Wang et al. (2019).

5.5 High-Energy GeV Emission

We turn back again to the already introduced inner engine. The joint action of rotation and magnetic field induces an electric potential Ruffini et al. (2018a, b)

[TABLE]

capable to accelerate protons to ultrarelativistic velocities and energies up to $\epsilon_{p}=e\Delta\phi\approx 10^{21}$ eV.

Along the rotation axis, there are no radiation losses and so the inner engine leads to UHECRs. In the off-polar directions, the protons radiate synchrotron photons, e.g., at GeV and TeV energies.

In Ruffini et al. (2018a) it has been estimated that the available electrostatic energy to accelerate protons is

[TABLE]

so the number of protons that the inner engine can accelerate is

[TABLE]

The timescale of the first elementary process is dictated by the acceleration time, i.e.,:

[TABLE]

so the emission power of the inner engine is approximately:

[TABLE]

The timescale of the subsequent processes depends crucially on the time required to rebuild the electric field. It has been shown that this condition implies an essential role of the density profile of the ionic matter surrounding the BH and its evolution with time Ruffini et al. (2018b, 2019a).

For a BH mass of the order of the NS critical mass, say $M\sim 3~{}M_{\odot}$, a BH spin parameter $\alpha\sim 0.3$, and a strength of the magnetic field $B_{0}\sim 10^{14}$ G, the above numbers are in agreement with the observed GeV emission data. See, for instance, in Ruffini et al. (2018b) and Ruffini et al. (2019a), respectively, the details of the analysis of GRB 130427A and GRB 190114C. We refer to Ruffini et al. (2018a, b) for details on the synchrotron emission of the accelerated protons in the above magnetic field.

5.6 Additional Considerations

The strong dependence of $P_{\rm max}$ on the initial mass of the NS companion opens the interesting possibility of producing XRFs and BdHNe from binaries with similar short (e.g., $P\sim$ few minutes) orbital periods and COcore properties: while a system with a massive (e.g., $\gtrsim$2 $M_{\odot}$) NS companion would lead to a BdHN, a system with a lighter (e.g., $\lesssim$1.4 $M_{\odot}$) NS companion would lead to an XRF. This predicts systems with a similar initial SN, leading to a similar $\nu$NS, but with different GRB prompt and afterglow emission. Given that the GRB energetics are different, the final SN kinetic energy should also be different being that it is larger for the BdHNe. This has been clearly shown by specific examples in Wang et al. (2019).

There are also additional novel features unveiled by the new 3D SPH simulations which can be observable in GRB light-curves and spectra, e.g.,:

1. (1)

the hypercritical accretion occurs not only on the NS companion but also on the $\nu$NS and with a comparable rate. 2. (2)

This implies that BdHNe might be also be able to form, in special cases, BH-BH binaries. Since the system remains bound the binary will quickly merge by emitting gravitational waves. Clearly, no electromagnetic emission is expected from these mergers. However, the typically large cosmological distances of BdHNe would make it extremely difficult to detect their gravitational waves e.g., by LIGO/Virgo. 3. (3)

Relatively weak SN explosions produce a long-lived hypercritical accretion process leading and enhance, at late times, the accretion onto the $\nu$NS. The revival of the accretion process at late times is a unique feature of our binary and does not occur for single SNe, namely in the absence of the NS companion. This feature increases the probability of detection of weak SNe by X-ray detectors via the accretion phase in an XRF/BdHN. 4. (4)

For asymmetric SN explosions the accretion rate shows a quasi-periodic behavior that might be detected by X-rays instruments, possibly allowing a test of the binary nature and the identification of the orbital period of the progenitor.

6 Post-Explosion Orbits and Formation of NS-BH Binaries

The SN explosion leaves as a central remnant a $\nu$NS and the induced gravitational collapse of the NS companion leads to BH formation. Therefore, BdHNe potentially leads to $\nu$NS-BH binaries, providing the binary keeps bound. This question was analyzed via numerical simulations in Fryer et al. (2015).

Typical binaries become unbound during an SN explosion because of mass loss and the momentum imparted (kick) to the $\nu$NS by the explosion. A classical astrophysical result shows that, assuming the explosion as instantaneous (sudden mass loss approximation), disruption occurs if half of the binary mass is lost. For this reason the fraction of massive binaries that can produce double compact-object binaries is usually found to be very low (e.g., $\sim$0.001–1%) Fryer et al. (1999); Dominik et al. (2012); Postnov and Yungelson (2014).

Assuming instantaneous mass loss, the post-explosion semi-major axis is Hills (1983):

[TABLE]

where $a_{0}$ and $a$ are the initial and final semi-major axes respectively, $M_{0}$ is the (initial) binary mass, $\Delta M$ is the change of mass (in this case the amount of mass loss), and $r$ is the orbital separation before the explosion. For circular orbits, the system is unbound if it loses half of its mass. For the very tight BdHNe, however, additional effects have to be taken into account to determine the fate of the binary.

The shock front in an SN moves at roughly $10^{4}$ km s*-1*, but the denser, lower-velocity ejecta, can move at velocities as low as $10^{2}$–$10^{3}$ km s*-1* Fryer et al. (2014). This implies that the SN ejecta overcomes an NS companion in a time 10–1000 s. For wide binaries this time is a small fraction of the orbital period and the “instantaneous” mass-loss assumption is perfectly valid. BdHNe have instead orbital periods as short as 100–1000 s, hence the instantaneous mass-loss approximation breaks down.

We recall the specific examples studied in Fryer et al. (2015): close binaries in an initial circular orbit of radius $7\times 10^{9}$ cm, COcore radii of $(1$–$4)\times 10^{9}$ cm with a 2.0 $M_{\odot}$ NS companion. The COcore leaves a central 1.5 $M_{\odot}$ NS, ejecting the rest of the core. The NS leads to a BH with a mass equal to the NS critical mass. For these parameters it was there obtained that even if 70% of the mass is lost the binary remains bound, providing the explosion time is of the order of the orbital period ($P=180$ s) with semi-major axes of less than $10^{11}$ cm (see Figure 12).

The tight $\nu$NS-BH binaries produced by BdHNe will, in due time, merge owing to the emission of gravitational waves. For the above typical parameters the merger time is of the order of $10^{4}$ year, or even less (see Figure 12). We expect little baryonic contamination around such merger site since this region has been cleaned-up by the BdHN. These conditions lead to a new family of sources which we have called ultrashort GRBs, U-GRBs.

7 BdHN Formation, Occurrence Rate and Connection with Short GRBs

7.1 An Evolutionary Scenario

The X-ray binary and SN communities have introduced a new evolutionary scenario for the formation of compact-object binaries (NS-NS or NS-BH). After the collapse of the primary star forming a NS, the binary undergoes mass-transfer episodes finally leading to the ejection of both the hydrogen and helium shells of the secondary star. These processes lead naturally to a binary composed of a COcore and an NS companion (see Figure 1). In the X-ray binary and SN communities these systems are called “ultra-stripped” binaries Tauris et al. (2015). These systems are expected to comprise $0.1$–$1\%$ of the total SNe Tauris et al. (2013).

The existence of ultra-stripped binaries supports our scenario from the stellar evolution side. In the above studies most of the binaries have orbital periods in the range $3\times 10^{3}$–$3\times 10^{5}$ s which are longer with respect to the short periods expected in the BdHN scenario. Clearly, XRF and BdHN progenitors should be only a small subset that result from the binaries with initial orbital separation and component masses leading to COcore-NS binaries with short orbital periods, e.g., $100$–$1000$ s for the occurrence of BdHNe. This requires fine-tuning both of the COcore mass and the binary orbit. From an astrophysical point of view the IGC scenario is characterized by the BH formation induced by the hypercritical accretion onto the NS companion and the associated GRB emission. Indeed, GRBs are a rare phenomenon and the number of systems approaching the conditions for their occurrence must be low (see Fryer et al. (2015) for details).

7.2 Occurrence Rate

If we assume that XRFs and BdHNe can be final stages of ultra-stripped binaries, then the percentage of the ultra-stripped population leading to these long GRBs must be very small. The observed occurrence rate of XRFs and BdHNe has been estimated to be $\sim$100 Gpc*-3* yr*-1* and $\sim$1 Gpc*-3* yr*-1*, respectively Ruffini et al. (2016), namely the $0.5\%$ and $0.005\%$ of the Ibc SNe rate, $2\times 10^{4}$ Gpc*-3* yr*-1* Guetta and Della Valle (2007). It has been estimated that $(0.1$–$1\%)$ of the SN Ibc could originate from ultra-stripped binaries Tauris et al. (2013), which would lead to an approximate density rate of $(20$–$200)$ Gpc*-3* yr*-1*. This would imply that a small fraction ($\lesssim$5%) of the ultra-stripped population would be needed to explain the BdHNe while, roughly speaking, almost the whole population would be needed to explain the XRFs (see Table 1). These numbers, while waiting for a confirmation by further population synthesis analyses, would suggest that most SNe originated from ultra-stripped binaries should be accompanied by an XRF. It is interesting that the above estimates are consistent with traditional estimates that only $\sim$0.001–1% of massive binaries lead to double compact-object binaries Fryer et al. (1999); Dominik et al. (2012); Postnov and Yungelson (2014).

7.3 Connection with Short GRBs

It is then clear that XRFs and BdHNe lead to $\nu$NS-NS and $\nu$NS-BH binaries. In due time, the emission of gravitational waves shrink their orbit leading to mergers potentially detectable as short GRBs. This implies a connection between the rate of long and short GRBs. It is clear from the derived rates (see Table 1 and Ruffini et al. (2016, 2018a)) that the short GRB population is dominated by the low-luminosity class of short Gamma-ray flashes (S-GRFs), double NS mergers that do not lead to BH formation. It can be seen that it is sufficient $\lesssim$4% of XRFs to explain the S-GRFs population, which would be consistent with the fact that many XRF progenitor binaries will get disrupted by the SN explosion. Therefore, by now, the observed rates of the GRB subclasses are consistent with the interesting possibility of a connection between the progenitors of the long and the ones of the short GRBs.

8 Conclusions

It is by now clear that short and long Gamma-ray bursts subclassify into eight different families and have as progenitors binary systems of a variety of flavors (see Table 1). We have focused in this work on the specific class of BdHNe of two types: type I and type II BdHNe, what in our old classification Ruffini et al. (2016) we called BdHNe and XRFs, respectively.

We have devoted this article mostly to the theoretical aspects of the induced gravitational collapse scenario and its evolution into BdHN as a complete model of long GRBs. We have also discussed, although briefly, the observable features of the model and how they compare with the observational data, providing to the reader the appropriate references for deepening this important aspect.

BdHNe I and II have as a common progenitor a CO${}_{\textrm{core}}$-NS binary. The CO${}_{\textrm{core}}$ explodes as type Ic SN, forming at its center a new NS, which we denote $\nu$NS, and produces onto the NS companion a hypercritical accretion process accompanied by an intense neutrino emission. The intensity of the accretion process and the neutrino emission depends mainly on the binary period, being more intense for tighter binaries. The NS companion in such an accretion process can reach or not the critical mass for gravitational collapse, i.e., to form a BH. The former binaries leading to a BH by accretion are the BdHNe I, while the ones in which the NS companion becomes just a more massive NS, are the BdHNe II (the old XRFs) (see Table 1).

We have reviewed the results of the numerical simulations performed of the above physical process starting from the 1D ones all the way to the latest 3D SPH ones. The simulation of this binary process has opened our eyes to a new reality: long GRBS are much richer and more complex systems than every one of us thought before, with the 3D morphology of the SN ejecta, that becomes asymmetric by the accretion process, playing a fundamental role in the GRB analysis.

We have recalled the relevance of each of the following processes in a BdHN:

1. (1)

the SN explosion; 2. (2)

the hypercritical accretion onto the NS companion; 3. (3)

the NS collapse with consequent BH formation; 4. (4)

the initiation of the inner engine; 5. (5)

the $e^{+}e^{-}$ plasma production; 6. (6)

the $e^{+}e^{-}$ plasma feedback onto the SN which converts the SN into a HN; 7. (7)

the formation of the cavity around the newborn BH; 8. (8)

the transparency of the $e^{+}e^{-}$ plasma along different directions; 9. (9)

the HN emission powered by the $\nu$NS; 10. (10)

the action of the inner engine in accelerating protons leading to UHECRs and to the high-energy emission.

The aforementioned involved physical processes in a BdHN have specific signatures observable (and indeed observed) in the long GRB multiwavelength lightcurves and spectra. We have recalled for each process its energetics, spectrum, and associated Lorentz factor: from the mildly-relativistic X-ray precursor, to the ultrarelativistic prompt Gamma-ray emission, to the mildly-relativistic X-ray flares of the early afterglow, to the mildly-relativistic late afterglow and to the high-energy GeV emission.

All of the above is clearly in contrast with a simple GRB model attempting to explain the entire GRB process with the kinetic energy of an ultrarelativistic jet extending through all of the above GRB phases, as in the traditional collapsar-fireball model.

If the binaries keep bound during the explosion, BdHNe I lead to $\nu$NS-BH binaries and BdHNe II lead to $\nu$NS-NS binaries. In due time, via gravitational wave emission, such binaries merge producing short GRBs. This unveiled clear interconnection between long and short GRBs and their occurrence rates needs to be accounted for in the cosmological evolution of binaries within population synthesis models for the formation of compact-object binaries.

We have taken the opportunity to include a brief summary of very recent developments published during the peer-review process of this article. These results cover the explanation of the observed GeV emission in BdHNe Ruffini et al. (2018a, b, 2019a, 2019b). One of the most relevant aspects of this topic is that it requests the solution of one of the fundamental problems in relativistic astrophysics: how to extract the rotational energy from a BH. This implies the role of a magnetic field around the newborn BH and the presence of surrounding matter as predicted in a BdHN. We have called this part of the system the inner engine of the high-energy emission. The BH rotation and surrounding magnetic field, for appropriate values induces an electric field via the Wald’s mechanism Wald (1974). Such an electric field is of paramount importance in accelerating surrounding protons to ultrarelativistic velocities leading to the high-energy emission via proton-synchrotron radiation. The details of this exciting new topic are beyond the scope of the present article but we encourage the reader to go through the above references for complementary details.

\funding

This research received no external funding.

\conflictsofinterest

The authors declare no conflict of interest.

\reftitle

References