Issue 
A&A
Volume 649, May 2021



Article Number  A75  
Number of page(s)  11  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201937135  
Published online  27 May 2021 
The newborn black hole in GRB 191014C proves that it is alive
^{1}
ICRANet, Piazza della Repubblica 10, 65122 Pescara, Italy
email: rahim.moradi@icranet.org, jorge.rueda@icra.it, ruffini@icra.it
^{2}
ICRA, Dipartimento di Fisica, Sapienza Università di Roma, P.le Aldo Moro 5, 00185 Rome, Italy
^{3}
INAF, Osservatorio Astronomico d’Abruzzo, Via M. Maggini snc, 64100 Teramo, Italy
^{4}
ICRANetFerrara, Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
^{5}
Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
^{6}
INAF, Istituto de Astrofisica e Planetologia Spaziali, Via Fosso del Cavaliere 100, 00133 Rome, Italy
^{7}
INAF, Viale del Parco Mellini 84, 00136 Rome, Italy
Received:
18
November
2019
Accepted:
2
April
2021
A multidecade theoretical effort has been devoted to finding an efficient mechanism to use the rotational and electrodynamical extractable energy of a KerrNewman black hole (BH), to power the most energetic astrophysical sources such as gammaray bursts (GRBs) and active galactic nuclei. We show an efficient general relativistic electrodynamical process which occurs in the “inner engine” of a binary driven hypernova. The inner engine is composed of a rotating Kerr BH of mass M and dimensionless spin parameter α, a magnetic field of strength B_{0} aligned and parallel to the rotation axis, and a very lowdensity ionized plasma. Here, we show that the gravitomagnetic interaction between the BH and the magnetic field induces an electric field that accelerates electrons and protons from the environment to ultrarelativistic energies emitting synchrotron radiation. We show that in GRB 190114C the BH of mass M = 4.4 M_{⊙}, α = 0.4, and B_{0} ≈ 4 × 10^{10} G can lead to a highenergy (≳GeV) luminosity of 10^{51} erg s^{−1}. The inner engine parameters are determined by requiring (1) that the BH extractable energy explains the GeV and ultrahighenergy emission energetics, (2) that the emitted photons are not subjected to magneticpair production, and (3) that the synchrotron radiation timescale agrees with the observed highenergy timescale. We find for GRB 190114C a clear jetted emission of GeV energies with a semiaperture angle of approximately 60° with respect to the BH rotation axis.
Key words: black hole physics / magnetic fields / gammaray burst: individual: GRB 190114C / gammaray burst: general
© ESO 2021
1. Introduction
Rotating black holes (BHs) have traditionally been described by the Kerr (Kerr 1963) and the KerrNewman metrics (Newman et al. 1965) which assume three conditions: (i) they are in a matter vacuum, (ii) they are embedded in an asymptotically flat spacetime, and (iii) they fulfill global stationarity. Under these conditions, BHs are just a sink of energy, namely “dead BHs”. The discovery of the reversible and irreversible transformations in both these spacetimes (Christodoulou 1970; Christodoulou & Ruffini 1971) opened the conceptual possibility of extracting both rotational and electromagnetic energy from a KerrNewman BH. These results also led to the asymptotic massenergy formula relating the mass M of a BH to three independent parameters, the irreducible mass M_{irr}, the charge Q, and the angular momentum J, soon confirmed by Hawking (1972). The perspective that up to 50% of the massenergy of a KerrNewman BH could be extracted directed the attention to the alternative view of “alive BHs” whose extractable energy could be used as an astrophysical source (see, e.g., “Introducing the black hole” Ruffini & Wheeler 1971a and “On the energetics of black holes” by R. Ruffini in DeWitt & DeWitt 1973).
Since then an efficient process has been sought that is able to power the most energetic astrophysical sources, gammaray bursts (GRBs), and active galactic nuclei (AGNs), using the extractable energy from a BH. The theoretical framework has been constantly evolving (see, e.g., Tursunov & Dadhich 2019, for a review on this topic). As we show in this paper, the recent discovery of the birth of a BH in GRB 130427A (Ruffini et al. 2019a; Rueda & Ruffini 2020) demonstrates that the Kerr BH harbored in the inner engine of this source is indeed an enormous source of gigaelectron volt (GeV) energy. The main topic of this article is to reach a deeper understanding of the process of rotational energy extraction by further identifying the astrophysical setting, the boundary conditions, and the basic new physical laws that allow this process to become observable. Specifically, for the case of GRB 190114C, our goal is to infer the values of the independent physical component, the spectral distribution of the highenergy GeV emission, and the geometrical properties of the GeV and ultrahighenergy emissions.
Our approach is based on the binarydriven hypernova (BdHN) model of long GRBs (Rueda & Ruffini 2012; Fryer et al. 2014; Becerra et al. 2016; Ruffini et al. 2019a). The BdHN progenitor is a binary system composed of a carbonoxygen (CO) star and a neutron star (NS) companion. The collapse of the iron core of the evolved CO star forms a newborn NS (νNS) at its center and expels the stellar outermost layers, hence leading to a supernova (SN) explosion. The SN ejecta produces a hypercritical accretion process both onto the νNS and onto the NS companion. For very compact binaries (orbital period on the order of 5 min), the NS companion reaches the critical mass rapidly (a few seconds), undergoes gravitational collapse, and forms a rotating BH. We have called these long GRBs in which there is BH formation, BdHNe of type I (BdHN I). Their isotropic energy release is in the range 10^{53} − 10^{54} erg. Numerical simulations of the above process in one, two, and three dimensions have been presented in Fryer et al. (2014), Becerra et al. (2015, 2016, 2019), respectively. Only a fraction of BdHNe form BHs (380 BdHNe I have been identified; see Ruffini et al. 2021). In progenitors with longer binary periods, on the order of hours, no BHs are formed; the outcome is a binary NS with long GRBs in the range 10^{51} − 10^{53} erg (Wang et al. 2019). For even longer binary periods, on the order of days, even less energetic long GRBs are encountered, the BdHNe III, for example the case of GRB 060218 (Liang et al., in prep.).
We now return to GRB 190114C. It has already been shown that the collapse of the CO star, which triggers the complete GRB process in the presence of a binary NS companion, leads to a SN creating an additional NS (i.e., the νNS). The SN process is observed earlier at Xray (up to a few keV) and soft gammaray (up to a few MeV) wavelengths, and it has been referred to as “SNrise” (see, e.g., Wang et al. 2019). For GRB 190114C this occurs in the restframe interval t_{rf} ≲ 1.99 s. It carries an energy of E_{SNrise} = 2.82 × 10^{52} erg and is characterized by a blackbody plus cutoff powerlaw spectrum (Liang et al. 2019). The short duration of the SNrise finds a natural explanation in the BdHN model. In a BdHN I the companion NS is separated at only 10^{10} − 10^{11} cm (i.e., about 1 lightsecond) from the CO star, implying that only the first spike becomes observable before the expanding SN ejecta triggers the hypercritical accretion process onto the NS companion (see, e.g., Becerra et al. 2019; Wang et al. 2019).
The newborn BH is embedded in the magnetic field inherited from the NS (Rueda et al. 2020), and sits at the center of a cavity of very low density (see Ruffini et al. 2019b, for numerical simulations) of material from the SN ejecta. For GRB 190114C such a density has been estimated to be on the order of 10^{−14} g cm^{−3}. The cavity is carved during the accretion and subsequent gravitational collapse of the NS leading to the BH. The magnetic field remains anchored to material that did not participate in the BH formation (see Rueda et al. 2020 for a detailed discussion on the magnetic field around the newborn Kerr BH in a BdHN I).
The Kerr BH in the cavity is therefore not isolated and acts in conjunction with a test magnetic field of strength B_{0}, aligned with the BH rotation axis. An additional important feature is that there is no vacuum surrounding the BH. As we show in this article, a fully ionized, very lowdensity plasma is essential to allow the electrodynamical performance of the energy extraction process by the inner engine, which is necessarily nonstationary.
The operation procedure of the inner engine leads the mass and spin of the BH to decrease as functions of time, while the BH irreducible mass (M_{irr}) remains constant. The electrons accelerate to ultrahigh energies at the expense of the BH extractable energy^{1}
obtainable from the BH massenergy formula (Christodoulou 1970; Christodoulou & Ruffini 1971; Hawking 1971)
where J and M are respectively the angular momentum and the mass of the BH.
As we explain in this article, using the mathematical role of the PapapetrouWald solution (Papapetrou 1966; Wald 1974), a profound change of paradigm in relativistic astrophysics has been made possible by the inner engine (Ruffini et al. 2019a), namely the introduction of the effective charge given by the product of J and B_{0}:
This effective charge originates from the gravitomagnetic interaction of the Kerr BH with the surrounding magnetic field, left over by the collapse of the accreting NS to the BH still rooted in the surrounding material (see, e.g., Rueda et al. 2020). The existence of this effective charge finally explains the success of utilizing the concept of a KerrNewman BH as a temporary step to approach the analysis of quantum electrodynamical processes in the field of a rotating BH (see, e.g., Damour & Ruffini 1975).
We are now able to elaborate, with the use of quantum electrodynamics and general relativity, a novel and physically more complete treatment of the GRB highenergy engine in a globally neutral system, therefore satisfying Eq. (2). Starting from these general premises, the main focus of this article is the role of the newborn BH in giving origin to the GeV emission observed by FermiLAT in the context of the inner engine of a BdHN I. In Sect. 2 we mathematically describe the electromagnetic field surrounding the Kerr BH following the PapapetrouWald solution of the EinsteinMaxwell equations (Papapetrou 1966; Wald 1974). Section 3 summarizes the operation of the inner engine, including its energy budget and electric potential energy available for the acceleration of charged particles around the BH. The particle motion along the BH rotation axis and its relation to the inner engine contribution to ultrahighenergy cosmic rays (UHECRs) is presented in Sect. 4. In Sect. 5 we estimate the energy loss by synchrotron radiation for electrons moving outside the BH rotation axis. We obtain there the typical electron Lorentz factor, the corresponding pitch angles leading to highenergy (≳GeV) photons, and the radiation timescale. Section 6 presents an estimate of the energy and angular momentum extracted to the Kerr BH in the emission process, in the radiation timescale, implied by the BH massenergy formula. In Sect. 7 we present our method of inferring the inner engine parameters, namely the BH mass and spin, and the magnetic field strength from the three conditions required (the observed highenergy emission covered by the extractable energy of the BH, the observed highenergy luminosity equal to the synchrotron radiation value, the emitted highenergy photons able to freely escape from the system). In Sect. 8 we apply this framework to the case of GRB 190114C obtaining the corresponding inner engine parameters. Section 9 is dedicated to a comparison of our results with previous literature results. Finally, we outline the conclusions in Sect. 10.
2. Electric and magnetic fields around the BH
Following the considerations presented in Ruffini et al. (2015, 2019a) corresponding to GRB 130427A, we turn to a quantitative estimate of the inner engine via a solution of the EinsteinMaxwell equations of a Kerr BH embedded in a test, asymptotically aligned, uniform magnetic field (Papapetrou 1966; Wald 1974), hereafter the PapapetrouWald solution.
The BH rotation and the aligned magnetic field induce an electric field with the following radial and polar components:
The magnetic field components are
where , , , and and are respectively the geometric mass and specific angular momentum of the BH (see Table 1).
Units of the relevant physical quantities used in this article in the cgsGaussian and geometric system of units.
We here use the locally nonrotating observer, also known as zero angular momentum observer (ZAMO; see Bardeen 1970; Bardeen et al. 1972). Therefore, the electromagnetic field components ((4a), (4b) and (5a), (5b)) differ from those presented in Ruffini et al. (2019a), where the Carter’s observer was used.
For moderate dimensionless spin values, α ≲ 0.7, where , the electric and magnetic fields are accurately represented by the firstorder expansion (α ≪ 1 or ):
Thus, the electric field is mainly radial and inwardly directed. The electric field decreases with the square of the distance; therefore, it is maximum at the BH horizon, , and on the rotation axis θ = 0, so . The electric field vanishes for 3 cos θ_{±} − 1 = 0. Therefore, it is inwardly directed in the northern hemisphere for spherical polar angles (measured clockwise) −θ_{±} < θ < θ_{±}, where (see Eqs. (6a), (6b) and Fig. 1). Because of the equatorial symmetry it also points inward in the southern hemisphere for π − θ_{±} < θ < π + θ_{±}. In these regions, electrons are outwardly accelerated. In the remaining regions the electric field reverses sign, becoming outwardly directed (see Fig. 1). The value of θ_{±} is indeed accurately given by the slowrotation approximation; for instance, a numerical calculation shows that θ_{±} ≈ 54.74° −59.76° for α = 0.01−0.99. We show below in this article that, the electrons located in these northern and southern hemisphere cones of semiaperture angle of ≈60°, are outwardly accelerated with the appropriate pitch angles leading to GeV photons (see Sect. 5 and Fig. 2 for details). Clearly, being anisotropic, this “jetted” emission is not always visible. This feature has been crucial for the inference of the morphology of the BdHN I from the highenergy (GeV) data of long GRBs (Ruffini et al. 2021).
Fig. 1. Electric (blue lines) and magnetic (golden lines) field lines of the PapapetrouWald solution in the xzplane in Cartesian coordinates. The BH spin parameter is set to a/M = 0.3 and the magnetic field and the BH spin are aligned and parallel. The background is a densityplot of the electric field energy density which is decreasing from red to blue. The BH horizon is the black disk. Distances are in units of M and the fields in units of B_{0}. Outward electron acceleration occurs in the region limited by the dashed black lines, i.e., where the electric field is inwardly directed. In the northern hemisphere it covers spherical polar angles (measured clockwise) −θ_{±} < θ < θ_{±}, where θ_{±} ≈ 55°. By equatorial symmetry, in the southern hemisphere, it covers π − θ_{±} < θ < π + θ_{±}. 
It can be also seen that the magnetic field is everywhere nearly aligned with the BH rotation axis; at any distance we have B_{z} ≫ B_{x} ∼ B_{y}, and at distances it is perfectly aligned (i.e., B_{x} → 0, B_{y} → 0, and B_{z} → B_{0}). All these features can be seen in Fig. 1, which shows the electric and magnetic field lines given by the general expressions given by Eqs. (4a), (4b) and (5a), (5b).
3. Operation of the inner engine
The operation of the inner engine is based on three components naturally present in a BdHN I:

the Kerr metric that describes the gravitational field produced by the newborn rotating BH;

an asymptotically uniform magnetic field around the newborn BH fulfilling the PapapetrouWald solution (see Sect. 2);

a very lowdensity plasma around the newborn BH composed of ions and electrons of 10^{−14} g cm^{−3} (Ruffini et al. 2019b).
The inner engine operation follows these precise steps:

The magnetic field and the BH rotation induce an electric field as given by the PapapetrouWald solution (see Sect. 2). For an aligned and parallel magnetic field to the BH spin, the electric field is nearly radial and inwardly directed at, and about, the BH rotation axis within an angle θ_{±} (see Fig. 1).

The induced electric field accelerates electrons outwardly. The number of electrons that can be accelerated is set by the energy stored in the electric field Rueda & Ruffini (2020):
Here Ω_{+} = c^{2}∂M/∂J = c α/(2 r_{+}) is the socalled BH angular velocity; m_{n} the neutron mass; and , λ_{Pl} = ℏ/(m_{Pl}c), and are respectively the Planck energydensity, length, and mass. These expressions evidence the nature of the underlying physical process generating the electric field and the BH horizon: the electrodynamics of the PapapetrouWald solution (Ruffini et al. 2019a), the origin of its magnetic field from the binary NS companion (Rueda et al. 2020), and the smooth formation of the BH from the induced gravitational collapse process (Rueda & Ruffini 2012). Additional details on the above formulation are presented in Rueda & Ruffini (2020).

The maximum possible electron acceleration and energy is set by the electric potential energy difference from the horizon to infinity can be written as (Rueda & Ruffini 2020)
where a = J/M.

Along the polar axis radiation losses are absent (see below in Sect. 4), while at offaxis latitudes (see below in Sect. 5) the accelerated electrons emit synchrotron radiation. The radiation timescale τ_{rad} must fulfill
where L_{GeV} is the observed GeV luminosity.

After this, the energy ℰ has been used and emitted. The process restarts with a new angular momentum J = J_{0} − ΔJ, being ΔJ the angular momentum extracted to the Kerr BH by the event (see below Eqs. (19a)–(19c) in Sect. 6).

The above steps are repeated, with the same efficiency, if the density of plasma is sufficient, namely if the number of the particles is enough to cover the new value of the energy ℰ. Therefore, the inner engine evolves in a sequence of elementary processes, each emitting a welldefined, precise amount of energy.
For the sake of example, let us chose fiducial parameters B_{0} = 10^{11} G, M = 3 M_{⊙}, and α = 0.5. In this case the available energy is ℰ ≈ 3.39 × 10^{37} erg, and the maximum energy that an accelerated electron can gain is ΔΦ ≈ 1.06 × 10^{7} erg = 6.64 × 10^{18} eV.
4. Acceleration on the polar axis: Ultrahighenergy cosmic rays
Along the polar axis, θ = 0, the electric and magnetic fields only have the zcomponent and are thus parallel; see Eqs. (4a), (4b) and (5a), (5b), or Eqs. (6a) and (6b). Since the electron is accelerated by the electric field, this implies that the electron pitch angle, which is the angle between the electron’s injection velocity (into the magnetic field) and the magnetic field, is zero. Consequently, no radiation losses (by synchrotron emission) occur for motion along the BH rotation axis.
The electrons accelerate outward gaining the total electric potential energy, ΔΦ ∼ 10^{18} eV. Therefore, the maximum number of electrons that the inner engine can accelerate along the axis is
These ultrarelativistic electrons contribute to leptonic UHECRs. The timescale of this acceleration process along the polar axis is
This implies that the inner engine can accelerate electrons along the BH rotation axis at a rate
leading to a power
Since the electric and magnetic fields along the rotation axis (and nearly close to it) are parallel (see Fig. 1), the particles in that region are all accelerated (nearly) parallel to the BH rotation axis. Therefore, we do not expect the accelerated particles to have appreciable collisions able to reduce the above estimate of their maximum kinetic energy gain. Therefore, given by Eq. (13) is the maximum power available for UHECRs.
The extension of the considerations presented here to very massive BHs and AGN, the role of the accretion disk in these galactic configurations, and the possibility of accelerating protons to produce UHECRs by the BH have started to be addressed (Rueda & Ruffini 2020). We compare and contrast in Table 2 some of the inner engine physical properties applied to the case of GRB 190114C and to M 87*.
Some astrophysical properties of the inner engine for GRB 190114C and AGN, in the latter adopting as a proxy M 87* (Rueda & Ruffini 2020).
5. Acceleration at offaxis latitudes: Synchrotron radiation
For the present electric field, and assuming radial motion, the dynamics of the electrons in the electromagnetic field, for γ ≫ 1, is determined from (de Jager et al. 1996; Ruffini et al. 2019a)
where e is the elementary charge, γ is the electron Lorentz factor, ⟨χ⟩ is the injection angle between the direction of electron motion and the magnetic field (the pitch angle), and m_{e} is the electron mass. This equation is integrated assuming the electrons are injected near the horizon (where the electric field strength is αB_{0}/2), for selected values of the injection angle ⟨χ⟩, with an initial Lorentz factor γ = 1 at t = 0.
The synchrotron spectrum peaks roughly at the photon critical energy (see, e.g., Landau & Lifshitz 1975)
where in the last expression we introduced β = B_{0}/B_{c}, with G. Therefore, the synchrotron peak energy shifts from lower to higher energies (softtohard spectral evolution) as the electron accelerates. For example, the photon critical energy ϵ_{γ}, for γ ≳ 10^{3}, a magnetic field B_{0} = 10^{11} G (so β = 0.0023), and a pitch angle χ = 10° falls in the GeV regime.
During the acceleration, the Lorentz factor increases linearly with time up to an asymptotic maximum value (see Ruffini et al. 2019a, for details). This maximum value, set by the balance between the energy gain by acceleration in the electric field and energy loss by synchrotron radiation, is (Ruffini et al. 2019a)
which defines the maximum electron energy ϵ_{e} = γ_{max}m_{e}c^{2}. Associated with γ_{max}, by replacing Eq. (16) into (15) we obtain the maximum peak energy of the spectrum (Ruffini et al. 2019a)
and the synchrotron cooling timescale t = t_{c} for the above maximum photon critical energy is given by (Ruffini et al. 2019a)
For model parameters α = 0.5 and B_{0} = 10^{11} G, photons of energy 0.1−10 GeV (typical photon energy range detected by the FermiLAT) are emitted by electrons with pitch angles χ ≈ 0.23−23°, and electron energy ϵ_{e} = 1.98 × 10^{8} − 1.98 × 10^{10} eV, radiating on a timescale of t_{c} = 2.63 × 10^{−17} − 2.63 × 10^{−15} s. We show in Fig. 2 the pitch angle χ as a function of the maximum photon critical energy (spectrum peak energy) ϵ_{γ, max}, obtained from Eq. (17), in the energy range 0.1−10 GeV, and for three selected values of α. Figure 3 shows the contours of constant χ for electrons moving in the electromagnetic field of the PapapetrouWald solution shown in Fig. 1. In particular, we show pitch angles for which electrons emit photons of GeV energies (see also Fig. 2). It can be seen that this highenergy jetted emission occurs within an effective opening angle θ_{±} ≈ 60°. This anisotropic emission is essential to infer the BdHN I morphology from the GeV emission data of long GRBs (Ruffini et al. 2021).
Fig. 2. Pitch angle χ (in units of degrees) as a function of the photon critical energy ϵ_{γ} (in units of GeV) obtained from Eq. (17). The focus of the plot for photon critical energy is in the range 0.1−10 GeV. The solid, dashed, and dotted curves correspond to selected values of the BH spin parameter α = 0.1, 0.3, and 0.5, respectively. 
Fig. 3. Contours of constant pitch angle (χ) of electrons moving in the electromagnetic field of the PapapetrouWald solution shown in Fig. 1. For the present magnetic dominated case (E/B< 1), charged particles follow the magnetic field lines; therefore, (see, e.g., Kelner et al. 2015), where E_{} is the electric field component parallel to the magnetic field, and Ψ is the angle between B and −E (the minus sign is used because we are interested in the pitch angle of electrons). In the slowrotation regime (see Eqs. (6a) and (6b)), sin Ψ ≈ sin θ/(1−2Msin^{2}θ/r), so sin χ ≈ E/Bsin θ/(1−2Msin^{2}θ/r). The BH is indicated by the black disk. The background color map indicates the electric field energy density (lighter colors means more intense). 
6. Energy and angular momentum extraction, golden rule, and duration of the inner engine activity
At the end of every elementary process, all the energy ℰ has been emitted. The inner engine restarts its operation with the same magnetic field of B_{0}, but with a new slightly smaller angular momentum J = J_{0} − ΔJ, being ΔJ the angular momentum extracted by the process. From the BH massformula (2), keeping the irreducible mass constant (i.e., ΔM_{irr} = 0, and ΔMc^{2} = ℰ), we obtain a change in the BH angular momentum ΔJ in each event:
Here the last equation, a truly golden formula, relates the energy radiated (ΔE = ℰ) to the angular momentum extracted to the rotating BH (ΔJ).
For the fiducial parameters that we used above, M = 3 M_{⊙}, α = 0.5, and B_{0} = 10^{11} G, we have J ≈ 3.96 × 10^{49} g cm^{2} s^{−1}, M_{irr} ≈ 2.9 M_{⊙}, and ΔJ ≈ 1.0 × 10^{33} g cm^{2} s^{−1}, so a fractional change ΔJ/J ≈ 10^{−16}, implying that the activity can last for thousands of years or more, providing there is ionized plasma to feed the inner engine.
7. Inference of the BH mass, spin, and surrounding magnetic field
We require three physical and astrophysical conditions to obtain the three inner engine parameters, the BH mass and spin, M and α, and the strength of the magnetic field surrounding the BH, B_{0}. Following Ruffini et al. (2019a), who show that the use of solely the GeV emission data, after the ultrarelativistic prompt emission (UPE) phase (see Liang et al. 2019), is enough to determine the inner engine parameters (see Fig. 4). In particular, we show that this procedure serves to obtain a lower limit to the mass and spin of the BH. The most important point is that we obtain the value of the irreducible mass of the BH that is kept constant through the energy extraction process. This allows us to determine the time evolution of the BH mass and spin. This can be achieved by fulfilling the three following conditions.
Fig. 4. Red data points: restframe 0.1−100 GeV luminosity light curve of GRB 190114C during and after UPE phase obtained from FermiLAT. Green dashed line: best fit for powerlaw behavior of the luminosity following the UPE phase with slope of 1.2 ± 0.04 and amplitude of 7.75 ± 0.44 × 10^{52} erg s^{−1}. Black data points: restframe 0.3−10 keV luminosity expressed in the rest frame obtained from SwiftXRT. It follows a powerlaw behavior with an amplitude of A_{X} = (5.14 ± 2.03) × 10^{52} erg s^{−1} and a slope of α_{X} = 1.37 ± 0.05. 
7.1. Condition 1
First, we require that the rotational energy of the BH provides the energy budget for the observed GeV emission energetics,
which via Eqs. (1) and (2) leads to the following inequality between M, α, and E_{GeV}:
We recall that the maximum value of the efficiency parameter is η_{max} ≈ 0.293, which is attained for a maximally rotating BH, α_{max} = 1. It is also important to recall that, by keeping the BH irreducible mass constant in the energy extraction process, we are inferring a lower limit to the BH mass. An increasing M_{irr} with time implies a higher BH mass to explain the same GeV energetics.
7.2. Condition 2
We require that the GeV photons must be transparent to the magnetic e^{+}e^{−} pair creation process. The attenuation coefficient for this process is (see Daugherty & Harding 1983 and Sect. 5 in Ruffini et al. 2019a)
where ψ = β sin ⟨χ⟩ϵ_{γ}/(2m_{e}c^{2}). Substituting Eq. (17) into Eq. (22), becomes a function of ϵ_{γ} and the product αβ. For a given ϵ_{γ} and α, the lower the magnetic field, the larger the mean free path , as expected. When χ ≪ 1, the exponential term dominates, hence exponentially increases tending to become infinite. An orderofmagnitude estimate of the magnetic field can be obtained by requiring ψ ≪ 1,
which is independent of the photon peak energy. It should be noted that this constraint already restricts the magnetic field to be undercritical (β < 1), and as we shall see it is sufficient for explaining the GeV emission after the UPE phase. The constraint (23) is analogous to imposing a lower limit on . For instance, adopting a photon energy of 0.1 GeV, it can be checked that for α β = 1.298 × 10^{−2}, the mean free path is cm. Lower values of α β lead to much larger values of . It is very interesting that this value is comparable to G M_{⊙}/c^{2} ≈ 1.477 × 10^{5} cm. Therefore, requesting a value of α β lower than the abovementioned one, implies having a mean free path that is much larger than the BH horizon. Specifically, the highenergy photons are produced in the vicinity of the BH, but they can freely escape from the system. If we adopt as a fiducial value that 0.1 GeV photons have a sufficiently large mean free path (e.g., cm), we obtain (Ruffini et al. 2019a)
That we are in the exponentially increasing part of the mean free path is evident by the fact that, by requesting a mean free path which is 11 orders of magnitude larger than the one implied by (23), our upper limit to the magnetic field is decreased less than one order of magnitude. Therefore, our estimate of the magnetic field is not sensitive to the choice of the value of , providing it satisfies ≳10^{5} cm. This implies that the magnetic field strength of the inner engine is constrained to have a value, roughly speaking, in the range 10^{10} − 10^{11} G (see Fig. 5).
Fig. 5. Parameters of the inner engine of GRB 190114C. For this source we have L_{GeV, 1} = 1.47 × 10^{52} erg s^{−1} and E_{GeV} = 1.8 × 10^{53} erg. Upper panel: family of solutions of B_{0} as a function of α (blue curve), given by Eq. (27). We use here a photon energy ϵ_{γ} = 0.1 GeV (lower edge of the FermiLAT energy band). In the gray shaded region the mean free path is cm, while in the white shaded region it is cm. The curve separating the two regions is therefore given by the equality in (24). Lower panel: corresponding family of solutions M(α) (red curve), given by Eq. (21). 
7.3. Condition 3
The third condition (i.e., the closure equation) is obtained by requesting that the timescale of the synchrotron radiation, the cooling time t_{c} given by Eq. (18), be equal to the observed GeV emission timescale (Ruffini et al. 2019a)
where ℰ is the electrostatic energy available for the process (see Eqs. (7a) and (7b)). The subscript “1” refers to quantities evaluated at the beginning of the transparency of the GeV emission (i.e., at the end of the UPE phase) at t = t_{rf, UPE} (see Sect. 8). We refer to this as the first elementary impulsive event. Therefore, the third equation of the system is
where μ = M/M_{⊙} and M_{⊙} is the solar mass.
Therefore, having imposed these three conditions, we obtain the three inner engine parameters from the system of Eqs. (21), (24), and (26), as follows:

We adopt the equality in Eq. (21), which implies that we will obtain a lower limit to the BH and spin;

We replace it into the equality of Eq. (24), which implies that we are adopting the upper limit to the magnetic field strength (for a given α);

We obtain the following expression for β as a function of α and of the observables E_{GeV} and L_{GeV} (Ruffini et al. 2019a):
where we have substituted Eq. (17) into Eq. (18) to express t_{c} as a function of the peak photon energy ϵ_{γ}, instead of the pitch angle χ.
Therefore, the BH horizon r_{+} is a function of μ and α, but in view of Eq. (21) it becomes a function of E_{GeV} and α. Given the observational quantities E_{GeV} (integrated after the UPE phase) and the luminosity L_{GeV, 1} (at the end of the UPE phase), Eq. (27) gives a family of solutions of β as a function of α. The solution of this equation together with Eq. (24) gives the values of β and α. With the knowledge of α and E_{GeV}, we obtain μ from Eq. (21).
8. Application to GRB 190114C
We now turn to apply the above procedure to GRB 190114C. For the observational properties of this source we follow the results of Liang et al. (2019). The UPE phase has been shown in that paper to end at the restframe time t_{rf, UPE} = 3.99 s, so we limit our analysis to longer times. Therefore, for GRB 190114C, we have t_{rad, 1} = t_{rf, UPE} = 3.99 s. The 0.1−100 GeV luminosity observed by FermiLAT, at t ≥ t_{rf, UPE}, is well fitted by a powerlaw function analogous to the case of GRB 130427A (Ruffini et al. 2019a) (see our Fig. 4):
The associated integrated isotropic energy observed by FermiLAT, from t_{1} = t_{rf, UPE} = 3.99 s to ∼10^{4} s, is E_{GeV} = (1.8 ± 1.3) × 10^{53} erg, and the luminosity at t = t_{rad, 1} is L_{GeV, 1} = 1.47 × 10^{52} erg s^{−1} (see Liang et al. 2019).
For the above numbers, and assuming the minimum energy budget requirement, specifically assuming the equality in Eq. (20), the inner engine parameters are (see Fig. 5) magnetic field B_{0} ≈ 3.9 × 10^{10} G, and spin and BH mass α = 0.41 and M = 4.45 M_{⊙}, respectively. The corresponding BH irreducible mass is M_{irr} = 4.35 M_{⊙}. For the above spin value, Eq. (17) leads to a pitch angle for the emission of 0.1 GeV photons, θ ≈ π/9.
The inequality (20) implies that the above mass and spin values of the BH must be considered as lower limits. As we show in Sect. 6, the inner engine can be longlasting so it can continue to emit and so will eventually radiate more than the observed E_{GeV} that we have used limiting ourselves to the first 10^{4} s of emission. However, in view of the powerlaw behavior of the GeV luminosity, most of the energy is emitted in this early evolution so the BH parameters do not change significantly if we consider the extrapolation of the energy budget. For example, let us assume that the powerlaw luminosity (28) extends for 1000 yr. This would increase the total GeV energy radiated by 25%, and recalculating all the parameters we obtain α = 0.46, M = 4.50 M_{⊙}, and M_{irr} = 4.38 M_{⊙}.
9. Comparison with previous literature
9.1. Long GRBs in the traditional model
We first recall some key features of the traditional model of long GRBs. To this end, we are facilitated by the book by Bing Zhang (Zhang 2018), which includes an extensive number of references. The traditional GRB model roots can be found in the papers by Rees & Meszaros (1992), Mészáros & Rees (1997), and Woosley (1993). The model proposed by Rees & Meszaros (1992) is based on a single system: GRBs are explained by a single BH from which an ultrarelativistic blastwave originates and whose expansion follows the Blandford–McKee selfsimilar solution (Blandford & McKee 1976). Woosley (1993) connected the GRB origin to a Kerr BH emitting an ultrarelativistic jet that originates from matter accretion onto the BH. The BH was proposed to form from the direct collapse of a massive star, called a failed SN or a “collapsar”, leading to a BH in the mass range 5−10 M_{⊙}.
In these models the afterglows are explained via the kinetic energy of the ultrarelativistic blastwave, which can reach very high bulk Lorentz factors, Γ ∼ 1000, to be released when it interacts with the circumburst medium (Waxman & Piran 1994; Sari & Piran 1995; Wijers et al. 1997; Sari 1997; Sari et al. 1998). The observed spectrum is proposed to be produced by synchrotron and synchrotron selfCompton (SSC) emission from accelerated electrons during the deceleration of the ultrarelativistic blastwave radiating at distances 10^{16} − 10^{18} cm. As pointed out by Zhang (2018), these models based on a ultrarelativistic blastwave radiating at distances 10^{16} − 10^{18} cm have been applied to jointly explain in the jetted emission several observations:

the Xray afterglow, the steep and shallow decay of the luminosity in the NousekZhang phase (Nousek et al. 2006; Zhang et al. 2006), the Xray flares and the gammaray flares;

the emission in the MeV and in the keV observed by the Neils Gehrels Swift Observatory, as well as the emission in the optical and in the radio, as well as the emission in the TeV recently observed by MAGIC (Mirzoyan et al. 2019; MAGIC Collaboration 2019a,b);

the highenergy (GeV) emission observed in some long GRBs by FermiLAT.
Within the traditional model, all the above emissions are explained using solely the kinetic energy of the ultrarelativistic blastwave with Γ ∼ 10^{3}, and radiating at 10^{16} − 10^{18} cm. It becomes clear that, requiring to the single kinetic energy of an ultrarelativistic blastwave to account for the entire energetics of all the observed radiation, at all wavelengths, from the prompt to the afterglow, results in an extreme request to the energy reservoir of the GRB engine.
Within the traditional collapsarfireball model, the presence of a mildly relativistic expanding component has been introduced in RamirezRuiz et al. (2002), called a cocoon, which moves sideways to the jet. However, as is clearly elucidated in Nakar & Piran (2017; see also references therein), the emission in the Xrays from this cocoon is too low with respect to be observed Xray afterglow of long GRBs, unless the cocoon Lorentz factor becomes Γ > 10. The possibility of a mildly relativistic component in the traditional model is interesting for its implications for the nature of the lowenergy sources such as GRB 060218 (see, e.g., Nakar 2015). These sources have energies 10^{49} − 10^{51} erg, which is the range of energies of BdHNe II and III. However, this is beyond the scope of this article, which is dedicated to BdHNe I, which are characterized by the energy range 10^{52} − 10^{54} erg. In conclusion, an explanation of the Xray afterglow in the traditional model needs ultrarelativistic values of the Lorentz factor (see also Zhang 2018, for a review on the subject).
9.2. Long GRBs in the BdHN model
As we note in Sect. 1, the BdHNe have a binary progenitor composed of a CO star and a companion NS. The GRB is composed of independent physical process identified by a timeresolved spectral analysis. Some key results are the following:

In the analysis of the data of the XRT detector on board the Neils Gehrels Swift satellite of the gammaray flare, the Xray flares, the flareplateau, and the early afterglow phases (the NousekZhang phase), after the ultrarelativistic prompt radiation phase, showed that the emitter in these phases is in mildly relativistic expansion with Γ ≲ 5 (see Ruffini et al. 2018a, for details). A similar upper limit Γ ≲ 3 was obtained in the case of GRB 151027A (Ruffini et al. 2018b), and for GRB 130427A the corresponding upper limit on the bulk Lorentz factor is Γ ≲ 2 (Ruffini et al. 2018c). Therefore, these stringent upper limits on Γ exclude any ultrarelativistic motion following the UPE phase, contrary to the prediction of traditional GRB models based on the ultrarelativistic blastwave.

The highenergy GeV emission follows from the action of the inner engine presented in this work, powered by the BH rotational energy extraction process. In the case of GRB 190114C studied in this work, this corresponds to t_{rf} ≳ 3.99 s (see Fig. 4 and Liang et al. 2019). It is characterized by an afterglow in the GeV radiation which, when expressed in the rest frame, follows a powerlaw luminosity (see Eq. (28) and Fig. 4), and it carries an energy of E_{GeV} = (1.8 ± 1.3) × 10^{53} erg.

In parallel, the Xray afterglow emission observed by the Swift satellite originates from the synchrotron radiation produced in the expanding SN ejecta, threaded by the magnetic field of the νNS, and aided by the injection of particles and the pulsarlike radiation from the νNS into the SN ejecta (Ruffini et al. 2018c; Wang et al. 2019; Rueda et al. 2020). These processes are mainly powered by the rotational energy of the νNS and have led to a significant progress in understanding the origin of the Xray afterglow emission (see, e.g., the case of GRB 130427A in Ruffini et al. 2018c, and GRB 160509A, GRB 160625B, GRB 180728A, and GRB 190114C in Rueda et al. 2020). In these analyses the spin of the νNS and the strength and structure of its magnetic field have been inferred. In the case of GRB 190114C, the luminosity expressed in the rest frame follows a powerlaw behavior L_{X} = A_{X}t^{−αX}, where A_{X} = (5.14 ± 2.03) × 10^{52} erg s^{−1} and α_{X} = 1.37 ± 0.05 and carries an energy E_{X} = 2.11 × 10^{52} erg; see Fig. 4 (Ruffini et al. 2021; see also Liang et al. 2019). This interpretation of the Xray afterglow in the BdHN model conforms with the observational upper limits on the Γ factor of the Xray afterglow emitter summarized in point 1 above (see Ruffini et al. 2018a, for details).
In this way, being the total energetics divided into the different components of the system and their associated different physical phenomena, the energetic request to each emission episode in the BdHN becomes affordable.
9.3. Process of BH energy extraction
Having indicated the main differences between the traditional GRB model and the BdHN model regarding the Xray and the GeV afterglow emissions, we focus now on the mechanism of the highenergy (GeV) emission, which is intimately related to the physics of the GRB central engine.
There is a vast literature devoted to magnetic fields around BHs and how they may act in a mechanism that could extract the rotational energy of a Kerr BH. An early attempt in the absence of a charge by a matterdominated magnetized plasma accreting in a disk around a preexisting Kerr BH was presented in Ruffini & Wilson (1975). The effective potential describing the circular orbit of massive particles around a Kerr BH was adopted (see Ruffini & Wheeler 1971b, in problem 2 of Sect. 104 in Landau & Lifshitz 1975). The infinite conductivity condition, F_{μν}u^{ν} = 0, where F_{μν} is the electromagnetic field tensor and u^{ν} the plasma fourvelocity, was used there leading to E ⋅ B = 0. Under these conditions, the acceleration of particles and processes of energy extraction were not possible.
This work was further developed by Blandford & Znajek (1977); in order to overcome the condition E ⋅ B = 0 in the magnetosphere, they adopted the concepts of gaps and spontaneous e^{+}e^{−} pair creation, closely following the seminal ideas of pulsar theory by Sturrock (1971) and Ruderman & Sutherland (1975). They imposed a forcefree condition, F_{μν}J^{ν} = 0, where J^{ν} is the current density, as well as gaps outside the BH horizon. The aim was to produce an ultrarelativistic matterdominated plasma whose bulk kinetic energy could be used to explain the energetics of a jet at large distances from the BH.
There is also another direction in the literature following the work of MacDonald & Thorne (1982). It extends the work of Blandford & Znajek (1977) and looks at the problem of matterdominated accretion in presence of a magnetic field anchored to a rotating surrounding disk. Specifically, they proposed an analogy of a rotating BH immersed in a magnetic field with a rotating conductive sphere and/or with the analogy of such a BH and the surrounding magnetosphere as an electric circuit. Independent of the analogies, the underlying physical system remains the same as that proposed by Blandford & Znajek (1977).
The present model is mainly motivated by fitting the GeV emission of GRBs. There is no matterdominated disk accretion. There is instead a very lowdensity ionized plasma fulfilling an acceleration electrodynamical process around a newly born BH. We use the PapapetrouWald solution in which the electromagnetic field is naturally characterized by regions where E ⋅ B ≠ 0 (see Sect. 2, Fig. 1, and Wald 1974). This feature naturally allows the acceleration of particles without the need of introducing any gaps. There is no ultrarelativistic matterdominated plasma outflow. The accelerated charged particles emit synchrotronradiation photons that carry off energy and angular momentum close to the BH. The BH in our scenario is not preexisting: it is smoothly formed by the hypercritical accretion onto the binary companion NS. The magnetic field, characterizing the PapapetrouWald solution, is amplified during the process of gravitational collapse of the binary companion NS (Rueda et al. 2020). There is no room in this model for the gravitational stable circular orbits around the Kerr BH. The particles are accelerated by an ultrarelativistic electrodynamical process.
Our description is also different with respect to recent GRB literature. For instance, in Metzger et al. (2011), Beniamini et al. (2017), and references therein, the presence of a magnetized wind, outflow, or jet is powered by a central engine. In these works the engine is represented by a NS endowed with an ultrahigh magnetic field, a magnetar, that loses its rotational energy via magneticdipole braking, in complete analogy to pulsars. The magnetar powers the outflows that produce the GRB emission at large radii on the order of 10^{15} cm. These models focus on the explanation of the (MeV) GRB prompt and the (Xray) afterglow emission using the rotational energy of a magnetar, so they do not look either to the physics of BHs, or to the GeV emission that are the topics of the present article.
The understanding of the complex nature of a BdHN requires the knowledge of different episodes, which in some cases are strictly independent, and their description can occur independently of each other.
For instance, the existence of hyperenergetic SN, the SNrise, radiates off 10^{52} erg in the case of GRB 190114C (Liang et al. 2019). In parallel, the interaction of the SN ejecta with the magnetic field of the νNS and its pulsarlike emission, explain the observed Xray afterglow (Ruffini et al. 2018c; Wang et al. 2019; Rueda et al. 2020). This emission is produced at distances 10^{12} − 10^{16} cm from the binary progenitor.
In the present work we address the most energetic GRB component, the GeV emission originating close to the horizon, at distances of 10^{6} cm, starting in the case of GRB 190114C at a restframe time of 3.99 s after the trigger.
After the clarification of these concepts, we will be ready to describe the optically thick subMeV emission in the time interval 1.99−3.99 s, which comprises the 55% of the energy of GRB 190114C, overcoming the compactness problem using our classic approach of the fireshell model (see Ruffini et al. 1999, 2000; Bianco et al. 2001; Moradi et al., in prep.).
10. Conclusions
The inner engine theory applied in this work to GRB 190114C represents an authentic full paradigm shift from the traditional model of long GRBs based on the emission of an ultrarelativistic blastwave, somehow powered by a Kerr BH. It seems too expensive for nature to accelerate matter in bulk, against the gravitational pull of the BH, to a large distance of ∼10^{16} − 10^{17} cm and with Γ ∼ 10^{3} to guarantee the transparency of highenergy radiation. For instance, the explanation of the GRB 190114C highenergy emission needs an ultrarelativistic blastwave with a kinetic energy on the order of 10^{55} erg (see, e.g., MAGIC Collaboration 2019a,b). It is clear that such energy cannot be powered by extracting the rotational energy of a Kerr BH of a few M_{⊙}, which will be a few 10^{53} erg (see Eq. (21)).
We have shown that the inner engine can nicely explain the GeV emission by accelerating electrons in the vicinity of the Kerr BH, which radiate their kinetic energy gain via synchrotron emission. The number of particles needed by the inner engine to explain the observed highenergy emission is relatively low. Let us adopt the derived inner engine parameter for GRB 190114C: M = 4.4 M_{⊙}, α = 0.4, and B_{0} = 4 × 10^{10} G. For instance, from Eq. (17) we obtain that for this α a photon peak energy of 10 GeV is obtained for an electron pitch angle χ ≈ 0.2° (see also Fig. 2). Using Eq. (16), this implies an electron Lorentz factor γ ≈ 6.76 × 10^{4}, which corresponds to an electron energy ϵ_{e} = γm_{e} c^{2} ≈ 5.53 × 10^{−2} erg = 3.45 × 10^{10} eV. Therefore, the number of such electrons needed to power the GeV emission of total energy E_{GeV} = 1.8 × 10^{53} erg ≈ 0.1 M_{⊙}c^{2}, is N_{e} = E_{GeV}/ϵ_{e} = 3.25 × 10^{54}, which for ionized matter implies a mass of m_{p}N_{e} ≈ 2.73 × 10^{−3} M_{⊙}, where m_{p} is the proton mass.
Therefore, the inner engine uses a more efficient electrodynamical process that produces observable highenergy emission in the vicinity of the BH. In fact the acceleration is not based on a bulkexpanding motion. Every single electron is accelerated from its initial velocity up to an asymptotic value defined by the maximum electric potential energy available for their acceleration, which depends only on the external magnetic field strength and the BH spin parameter; see Eqs. (8a) and (8b). These accelerated electrons radiate mainly at high energies in the GeV domain. The radiation of the inner engine (e.g., at keV to MeV energies) is negligible (with respect to the observed values). The observed radiation in the keV to MeV energy domains is explained by a different mechanism in a BdHN I; see (Rueda et al. 2020). The observed luminosity of GeV allows us to estimate the mass and spin of the BH.
We have determined the parameters of the inner engine of GRB 190114C using only the GeV emission data after the UPE phase. We asked the system to satisfy three physical conditions. First, that the GeV energetics is paid by the extractable energy of the BH (see Eq. (21)); second that the system is transparent to GeV photons produced by the synchrotron radiation of the accelerated electrons (see Eq. (24)); and third that the synchrotron radiation timescale explains the observed GeV emission timescale (see Eq. (26)) with the aid of Eq. (18). In order to be fulfilled, this last constraint implies that the GeV emission is emitted from electrons being accelerated with the appropriate pitch angles (see Figs. 2 and 3). These pitch angles occur within a cone of approximately 60° from the BH rotation axis (see Fig. 3), which is a key result for the interpretation of the morphology of the BdHN I (Ruffini et al. 2021).
From this procedure, we have obtained the inner engine parameters of GRB 190114C: B_{0} ≈ 3.9 × 10^{10} G, α ≈ 0.41, and M = 4.45 M_{⊙}. The corresponding irreducible mass of the BH is M_{irr} = 4.35 M_{⊙}. It is worth recalling that both M_{irr} and B_{0} are kept constant and this should be all over the evolution. The corresponding BH parameters for GRB 130427A are dimensionless spin α = 0.47, mass M = 2.3 M_{⊙}, and irreducible mass M_{irr} = 2.2 M_{⊙} (Ruffini et al. 2019a). The above are the first two BH masses derived directly from the GRB observations, and in both cases they are above the theoretical values of the NS critical mass enforcing the validity of the BdHN I model: the BH are formed by smooth hypercritical accretion of the HN ejecta on the NS binary companion.
Since here we only used the GeV emission data, the BH parameters that we have obtained, namely mass and spin, have to be considered as lower limits. Thus, it is clear that even a slightly higher mass (or spin) of the BH can guarantee even larger and longer emission of the inner engine.
Our analysis paves the way to additional research; the data from the different energy bands (e.g., the higher energy bands; MAGIC Collaboration 2019a,b) might provide additional information on the energy distribution of the electrons injected by the electric field into the magnetic field, and on the pitch angle distribution for the synchrotron emission. Figure 3 shows, for the electromagnetic field configuration of the PapapetrouWald solution (see Fig. 1), the contours of constant pitch angle and constant electric energy density.
Before concluding, it is worth recalling some crucial aspects of the inner engine here applied to the case of GRB 190114C. The nature of the emission results from considering the physical process leading to the electric and magnetic fields and the BH formation (see Sect. 3 and Rueda et al. 2020). This is fundamental to show that the emission process leading to the observed luminosity is not continuous, but discrete. The timescale of the emission in GRBs is too short to be probed directly by current observational facilities. Direct evidence of the value and discreteness might come instead from the observation of large BHs of 10^{8} − 10^{10} M_{⊙} in AGN. For instance, in the case of M 87*, for fiducial parameters M = 6 × 10^{9} M_{⊙}, α = 0.1, and B_{0} = 10 G, the inner engine theory predicts a highenergy (GeV) emission with a luminosity of a few 10^{43} erg s^{−1}, with a timescale of up to tenths of seconds (see Table 2). Emission at higher energies (e.g., in the TeV band), would be characterized by a lower luminosity and a longer timescale. The timescale for UHECR emission is instead approximately half a day (see Table 2 and Rueda & Ruffini 2020).
We can therefore conclude, in the light of the results of this article and the previous articles in this series, that all BdHN I are powered by three independent sources of energy. The BdHN I is triggered by the SN explosion originating from the collapse of the CO_{core} generating a νNS. The accretion of the SN onto the νNS (see Sect. 9.2 and Ruffini et al. 2018c; Wang et al. 2019; Rueda et al. 2020), gives origin to the Xray afterglow observed by Swift. The hypercritical accretion of the SN onto the binary companion NS gives origin to the BH as soon as the NS reaches the critical mass. This smooth accretion process is alternative to the direct gravitational collapse of a massive star. This happens in GRB 190114C at t_{rf} = 1.99 s. The further accretion of the SN ejecta onto the newly born BH generates the prompt gammaray radiation observed in GRB 190114C between 1.99 s and 3.99 s (Moradi et al., in prep.). The further accretion of the SN ejecta onto the newly born BH leads to a process of energy extraction from the inner engine that generates the jetted highenergy (≳GeV) emission. This radiation, as is shown in this article using the PapapetrouWald solution (see Sect. 2), is emitted close to the BH horizon and within an angle of nearly 60° from the BH rotation axis (see Sect. 5 and Fig. 3).
We use cgsGaussian units throughout, unless otherwise specified. Careted symbols stand for quantities in geometric units; for example, denotes geometric mass. See Table 1 for details on the units and conversion factors between the cgsGaussian and geometric systems of units.
Acknowledgments
We thank the Referee for the reiterated suggestions which have certainly improved the presentation of our results. An extended correspondence with the Referee has been addressed since our initial submission in order to broader the context of our work with respect to the previous literature and which have materialized in Sect. 9.
References
 Bardeen, J. M. 1970, ApJ, 162, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347 [NASA ADS] [CrossRef] [Google Scholar]
 Becerra, L., Cipolletta, F., Fryer, C. L., Rueda, J. A., & Ruffini, R. 2015, ApJ, 812, 100 [CrossRef] [Google Scholar]
 Becerra, L., Bianco, C. L., Fryer, C. L., Rueda, J. A., & Ruffini, R. 2016, ApJ, 833, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Becerra, L., Ellinger, C. L., Fryer, C. L., Rueda, J. A., & Ruffini, R. 2019, ApJ, 871, 14 [Google Scholar]
 Beniamini, P., Giannios, D., & Metzger, B. D. 2017, MNRAS, 472, 3058 [NASA ADS] [CrossRef] [Google Scholar]
 Bianco, C. L., Ruffini, R., & Xue, S.S. 2001, A&A, 368, 377 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Blandford, R. D., & McKee, C. F. 1976, Phys. Fluids, 19, 1130 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 [Google Scholar]
 Christodoulou, D. 1970, Phys. Rev. Lett., 25, 1596 [NASA ADS] [CrossRef] [Google Scholar]
 Christodoulou, D., & Ruffini, R. 1971, Phys. Rev. D, 4, 3552 [NASA ADS] [CrossRef] [Google Scholar]
 Damour, T., & Ruffini, R. 1975, Phys. Rev. Lett., 35, 463 [Google Scholar]
 Daugherty, J. K., & Harding, A. K. 1983, ApJ, 273, 761 [Google Scholar]
 de Jager, O. C., Harding, A. K., Michelson, P. F., et al. 1996, ApJ, 457, 253 [NASA ADS] [CrossRef] [Google Scholar]
 DeWitt, C., & DeWitt, B. S. 1973, in Proceedings, Ecole d’Eté de Physique Théorique: Les Astres Occlus (New York, NY: Gordon and Breach), Les Houches Summer School, 23 [Google Scholar]
 Fryer, C. L., Rueda, J. A., & Ruffini, R. 2014, ApJ, 793, L36 [NASA ADS] [CrossRef] [Google Scholar]
 Hawking, S. W. 1971, Phys. Rev. Lett., 26, 1344 [NASA ADS] [CrossRef] [Google Scholar]
 Hawking, S. W. 1972, Commun. Math. Phys., 25, 152 [Google Scholar]
 Kelner, S. R., Prosekin, A. Y., & Aharonian, F. A. 2015, AJ, 149, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Kerr, R. P. 1963, Phys. Rev. Lett., 11, 237 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Landau, L. D., & Lifshitz, E. M. 1975, The Classical Theory of Fields (Oxford: Pergamon Press) [Google Scholar]
 Liang, L., Ruffini, R., Rueda, J. A., et al. 2019, ArXiv eprints [arXiv:1910.12615] [Google Scholar]
 MacDonald, D., & Thorne, K. S. 1982, MNRAS, 198, 345 [Google Scholar]
 MAGIC Collaboration (Acciari, V. A., et al.) 2019a, Nature, 575, 455 [Google Scholar]
 MAGIC Collaboration (Acciari, V. A., et al.) 2019b, Nature, 575, 459 [Google Scholar]
 Mészáros, P., & Rees, M. J. 1997, ApJ, 482, L29 [Google Scholar]
 Metzger, B. D., Giannios, D., Thompson, T. A., Bucciantini, N., & Quataert, E. 2011, MNRAS, 413, 2031 [NASA ADS] [CrossRef] [Google Scholar]
 Mirzoyan, R., Noda, K., Moretti, E., et al. 2019, GCN Circ., 23701, 1 [Google Scholar]
 Nakar, E. 2015, ApJ, 807, 172 [Google Scholar]
 Nakar, E., & Piran, T. 2017, ApJ, 834, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Newman, E. T., Couch, E., Chinnapared, K., et al. 1965, J. Math. Phys., 6, 918 [CrossRef] [Google Scholar]
 Nousek, J. A., Kouveliotou, C., Grupe, D., et al. 2006, ApJ, 642, 389 [NASA ADS] [CrossRef] [Google Scholar]
 Papapetrou, A. 1966, Annales de L’Institut Henri Poincare Section (A) Physique Theorique, 4, 83 [Google Scholar]
 RamirezRuiz, E., Celotti, A., & Rees, M. J. 2002, MNRAS, 337, 1349 [NASA ADS] [CrossRef] [Google Scholar]
 Rees, M. J., & Meszaros, P. 1992, MNRAS, 258, 41P [Google Scholar]
 Ruderman, M. A., & Sutherland, P. G. 1975, ApJ, 196, 51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rueda, J. A., & Ruffini, R. 2012, ApJ, 758, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Rueda, J. A., & Ruffini, R. 2020, Eur. Phys. J. C, 80, 300 [EDP Sciences] [Google Scholar]
 Rueda, J. A., Ruffini, R., Karlica, M., Moradi, R., & Wang, Y. 2020, ApJ, 893, 148 [Google Scholar]
 Ruffini, R., & Wheeler, J. A. 1971a, Phys. Today, 24, 30 [CrossRef] [Google Scholar]
 Ruffini, R., & Wheeler, J. A. 1971b, ESRO, 52, 45 [Google Scholar]
 Ruffini, R., & Wilson, J. R. 1975, Phys. Rev. D, 12, 2959 [NASA ADS] [CrossRef] [Google Scholar]
 Ruffini, R., Salmonson, J. D., Wilson, J. R., & Xue, S.S. 1999, A&AS, 138, 511 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruffini, R., Salmonson, J. D., Wilson, J. R., & Xue, S.S. 2000, A&A, 359, 855 [Google Scholar]
 Ruffini, R., Wang, Y., Enderli, M., et al. 2015, ApJ, 798, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Ruffini, R., Wang, Y., Aimuratov, Y., et al. 2018a, ApJ, 852, 53 [Google Scholar]
 Ruffini, R., Becerra, L., Bianco, C. L., et al. 2018b, ApJ, 869, 151 [Google Scholar]
 Ruffini, R., Karlica, M., Sahakyan, N., et al. 2018c, ApJ, 869, 101 [Google Scholar]
 Ruffini, R., Moradi, R., Rueda, J. A., et al. 2019a, ApJ, 886, 82 [Google Scholar]
 Ruffini, R., Melon Fuksman, J. D., & Vereshchagin, G. V. 2019b, ApJ, 883, 191 [Google Scholar]
 Ruffini, R., Moradi, R., Rueda, J. A., et al. 2021, MNRAS, in press [arXiv:2103.09142] [Google Scholar]
 Sari, R. 1997, ApJ, 489, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Sari, R., & Piran, T. 1995, ApJ, 455, L143 [NASA ADS] [CrossRef] [Google Scholar]
 Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Sturrock, P. A. 1971, ApJ, 164, 529 [NASA ADS] [CrossRef] [Google Scholar]
 Tursunov, A., & Dadhich, N. 2019, Universe, 5, 125 [Google Scholar]
 Wald, R. M. 1974, Phys. Rev. D, 10, 1680 [Google Scholar]
 Wang, Y., Rueda, J. A., Ruffini, R., et al. 2019, ApJ, 874, 39 [Google Scholar]
 Waxman, E., & Piran, T. 1994, ApJ, 433, L85 [Google Scholar]
 Wijers, R. A. M. J., Rees, M. J., & Meszaros, P. 1997, MNRAS, 288, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Woosley, S. E. 1993, ApJ, 405, 273 [Google Scholar]
 Zhang, B. 2018, The Physics of Gammaray Bursts (Cambridge: Cambridge University Press) [Google Scholar]
 Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Units of the relevant physical quantities used in this article in the cgsGaussian and geometric system of units.
Some astrophysical properties of the inner engine for GRB 190114C and AGN, in the latter adopting as a proxy M 87* (Rueda & Ruffini 2020).
All Figures
Fig. 1. Electric (blue lines) and magnetic (golden lines) field lines of the PapapetrouWald solution in the xzplane in Cartesian coordinates. The BH spin parameter is set to a/M = 0.3 and the magnetic field and the BH spin are aligned and parallel. The background is a densityplot of the electric field energy density which is decreasing from red to blue. The BH horizon is the black disk. Distances are in units of M and the fields in units of B_{0}. Outward electron acceleration occurs in the region limited by the dashed black lines, i.e., where the electric field is inwardly directed. In the northern hemisphere it covers spherical polar angles (measured clockwise) −θ_{±} < θ < θ_{±}, where θ_{±} ≈ 55°. By equatorial symmetry, in the southern hemisphere, it covers π − θ_{±} < θ < π + θ_{±}. 

In the text 
Fig. 2. Pitch angle χ (in units of degrees) as a function of the photon critical energy ϵ_{γ} (in units of GeV) obtained from Eq. (17). The focus of the plot for photon critical energy is in the range 0.1−10 GeV. The solid, dashed, and dotted curves correspond to selected values of the BH spin parameter α = 0.1, 0.3, and 0.5, respectively. 

In the text 
Fig. 3. Contours of constant pitch angle (χ) of electrons moving in the electromagnetic field of the PapapetrouWald solution shown in Fig. 1. For the present magnetic dominated case (E/B< 1), charged particles follow the magnetic field lines; therefore, (see, e.g., Kelner et al. 2015), where E_{} is the electric field component parallel to the magnetic field, and Ψ is the angle between B and −E (the minus sign is used because we are interested in the pitch angle of electrons). In the slowrotation regime (see Eqs. (6a) and (6b)), sin Ψ ≈ sin θ/(1−2Msin^{2}θ/r), so sin χ ≈ E/Bsin θ/(1−2Msin^{2}θ/r). The BH is indicated by the black disk. The background color map indicates the electric field energy density (lighter colors means more intense). 

In the text 
Fig. 4. Red data points: restframe 0.1−100 GeV luminosity light curve of GRB 190114C during and after UPE phase obtained from FermiLAT. Green dashed line: best fit for powerlaw behavior of the luminosity following the UPE phase with slope of 1.2 ± 0.04 and amplitude of 7.75 ± 0.44 × 10^{52} erg s^{−1}. Black data points: restframe 0.3−10 keV luminosity expressed in the rest frame obtained from SwiftXRT. It follows a powerlaw behavior with an amplitude of A_{X} = (5.14 ± 2.03) × 10^{52} erg s^{−1} and a slope of α_{X} = 1.37 ± 0.05. 

In the text 
Fig. 5. Parameters of the inner engine of GRB 190114C. For this source we have L_{GeV, 1} = 1.47 × 10^{52} erg s^{−1} and E_{GeV} = 1.8 × 10^{53} erg. Upper panel: family of solutions of B_{0} as a function of α (blue curve), given by Eq. (27). We use here a photon energy ϵ_{γ} = 0.1 GeV (lower edge of the FermiLAT energy band). In the gray shaded region the mean free path is cm, while in the white shaded region it is cm. The curve separating the two regions is therefore given by the equality in (24). Lower panel: corresponding family of solutions M(α) (red curve), given by Eq. (21). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.