We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly anisotropic radiation transport as can be expected in the vicinity of an accreting planet in a protoplanetary disk. To reduce the computational cost caused by the timescale disparity between radiation and matter fields, we integrate their evolution equations separately in an operator-split way, using substepping to evolve the radiation equations. We further increase the code's efficiency by adopting the reduced speed of light approximation (RSLA). Our integration scheme for the evolution equations of radiation fields relies on implicit-explicit schemes, in which radiation-matter interaction terms are integrated implicitly while fluxes are integrated via Godunov-type solvers. The module is suitable for general astrophysical computations in 1, 2, and 3 dimensions in Cartesian, spherical and cylindrical coordinates, and can be implemented on rotating frames. We demonstrate the algorithm performance on different numerical benchmarks, paying particular attention to the applicability of the RSLA for computations of physical processes in protoplanetary disks. We show 2D simulations of vertical convection in disks and 3D simulations of gas accretion by planetary cores, which are the first of their kind to be solved with a two-moment approach.

We estimate the neutrino emission from the decay chain of the $\pi$-meson and $\mu$-lepton, produced by proton-proton inelastic scattering in energetic ($E_{\rm iso}\gtrsim 10^{52}$ erg) long gamma-ray bursts (GRBs), within the type I binary-driven hypernova (BdHN) model. The BdHN I progenitor is binary system composed of a carbon-oxygen star (CO$_{\rm core}$) and a neutron star (NS) companion. The CO$_{\rm core}$ explosion as supernova (SN) triggers a massive accretion process onto the NS. For short orbital periods of few minutes, the NS reaches the critical mass, hence forming a black hole (BH). Recent numerical simulations of the above scenario show that the SN ejecta becomes highly asymmetric, creating a \textit{cavity} around the newborn BH site, due to the NS accretion and gravitational collapse. Therefore, the electron-positron ($e^{\pm}$) plasma created in the BH formation, during its isotropic and expanding self-acceleration, engulfs different amounts of ejecta baryons along different directions, leading to a direction-dependent Lorentz factor. The protons engulfed inside the high-density ($\sim10^{23}$~particle/cm$^3$) ejecta reach energies in the range $1.24\lesssim E_p\lesssim 6.14$ GeV and interact with the unshocked protons in the ejecta. The protons engulfed from the low density region around the BH reach energies $\sim 1$ TeV and interact with the low-density ($\sim1$~particle/cm$^3$) protons of the interstellar medium (ISM). The above interactions give rise, respectively, to neutrino energies $E_{\nu}\leq 2$ GeV and $10\leq E_{\nu}\leq 10^3$ GeV, and for both cases we calculate the spectra and luminosity.

Within the binary-driven hypernova I (BdHN I) scenario, the gamma-ray burst GRB190114C originates in a binary system composed of a massive carbon-oxygen core (CO$_{core}$), and a binary neutron star (NS) companion. As the CO$_{core}$ undergoes a supernova explosion with the creation of a new neutron star ($\nu$NS), hypercritical accretion occurs onto the companion binary neutron star until it exceeds the critical mass for gravitational collapse. The formation of a black hole (BH) captures $10^{57}$ baryons by enclosing them within its horizon, and thus a cavity of approximately $10^{11}$ cm is formed around it with initial density $10^{-7}$ g/cm$^3$. A further depletion of baryons in the cavity originates from the expansion of the electron-positron-photon ($e^{+}e^{-}\gamma$) plasma formed at the collapse, reaching a density of $10^{-14}$ g/cm$^3$ by the end of the interaction. It is demonstrated here using an analytical model complemented by a hydrodynamical numerical simulation that part of the $e^{+}e^{-}\gamma$ plasma is reflected off the walls of the cavity. The consequent outflow and its observed properties are shown to coincide with the featureless emission occurring in a time interval of duration $t_{rf}$, measured in the rest frame of the source, between $11$ and $20$ s of the GBM observation. Moreover, similar features of the GRB light curve were previously observed in GRB 090926A and GRB 130427A, all belonging to the BdHN I class. This interpretation supports the general conceptual framework presented in Ruffini et al. (2019) and guarantees that a low baryon density is reached in the cavity, a necessary condition for the operation of the "inner engine" of the GRB presented in an accompanying article (Ruffini & Moradi 2019).

We motivate how the most recent progress in the understanding the nature of the GeV radiation in most energetic gamma-ray bursts (GRBs), the binary-driven hypernovae (BdHNe), has led to the solution of a forty years unsolved problem in relativistic astrophysics: how to extract the rotational energy from a Kerr black hole for powering synchrotron emission and ultra high-energy cosmic rays. The "inner engine" is identified in the proper use of a classical solution introduced by Wald in 1974 duly extended to the most extreme conditions found around the newborn black hole in a BdHN. The energy extraction process occurs in a sequence impulsive processes each accelerating protons to $10^{21}$ eV in a timescale of $10^{-6}$ s and in presence of an external magnetic field of $10^{14}$ G. Specific example is given for a black hole of initial angular momentum $J=0.3\,M^2$ and mass $M\approx 3\,M_\odot$ leading to the GeV radiation of $10^{49}$ erg$\cdot$s$^{-1}$. The process can energetically continue for thousands of years.

Following Fermi and NOT observations, Ruffini et al. (2019b) soon identified GRB 190114C as BdHN I at z=0.424, it has been observed since, with unprecedented accuracy, [...] all the way to the successful optical observation of our predicted supernova (SN). This GRB is a twin of GRB 130427A. Here we take advantage of the GBM data and identify in it three different Episodes. Episode 1 represents the precursor which includes the SN breakout and the creation of the new neutron star ($\nu$NS), the hypercritical accretion of the SN ejecta onto the NS binary companion, exceeding the NS critical mass at $t_{rf}=1.9$s. Episode 2 starting at $t_{rf}=1.9$s includes three major events: the formation of the BH, the onset of the GeV emission and the onset of the ultra-relativistic prompt emission (UPE), which extends all the way up to $t_{rf}=3.99$s. Episode 3 which occurs at times following $t_{rf}=3.99$s reveals the presence of a cavity carved out in the SN ejecta by the BH formation. We perform an in depth time-resolved spectral analysis on the entire UPE with the corresponding determination of the spectra best fit by a cut-off power-law and a black body (CPL+BB) model, and then we repeat the spectral analysis in 5 successive time iterations in increasingly shorter time bins: we find a similarity in the spectra in each stage of the iteration revealing clearly a self-similar structure. We find a power-law dependence of the BB temperature with index $-1.56\pm0.38$, a dependence with index $-1.20\pm0.26$ for the gamma-ray luminosity confirming a similar dependence with index $-1.20\pm0.36$ which we find as well in the GeV luminosity, both expressed in the rest-frame. We thus discover in the realm of relativistic astrophysics the existence of a self-similar physical process and power-law dependencies, extensively described in the micro-physical world by the classical works of Heisenberg-Landau-Wilson.

We present a numerical implementation for the solution of the relativistic radiation hydrodynamics and magnetohydrodynamics equations, designed as an independent module within the freely available code PLUTO. The radiation transfer equations are solved under the grey approximation and imposing the M1 closure, which allows the radiation transport to be handled in both the free-streaming and diffusion limits. Equations are integrated following an implicit-explicit scheme, where radiation-matter interaction terms are integrated implicitly, whereas transport and all of the remaining source terms are solved explicitly by means of the same Godunov-type solvers included in PLUTO. Among these, we introduce a new Harten-van Leer-contact (HLLC) solver for optically thin radiation transport. The code is suitable for multidimensional computations in Cartesian, spherical and cylindrical coordinates, using either a single processor or parallel architectures. Adaptive grid computations are also made possible, by means of the CHOMBO library. The algorithm performance is demonstrated through a series of numerical benchmarks by investigating various different configurations with a particular emphasis on the behavior of the solutions in the free-streaming and diffusion limits.

We analyze GRB 151027A within the binary-driven hypernova (BdHN) approach, with progenitor a carbon-oxygen core on the verge of a supernova (SN) explosion and a binary companion neutron star (NS). The hypercritical accretion of the SN ejecta onto the NS leads to its gravitational collapse into a black hole (BH), to the emission of the GRB and to a copious $e^+e^-$ plasma. The impact of this $e^+e^-$ plasma on the SN ejecta explains {the} early soft X-ray flare observed in long GRBs. We here apply this approach to the UPE and to the hard X-ray flares. We use GRB 151027A as a prototype. From the time-integrated and the time-resolved analysis we identify a double component in the UPE and confirm its ultra-relativistic nature. We confirm the mildly-relativistic nature of the soft X-ray flare, of the hard X-ray flare and of the ETE. We show that the ETE identifies the transition from a SN to the HN. We then address the theoretical justification of these observations by integrating the hydrodynamical propagation equations of the $e^+ e^-$ into the SN ejecta, the latter independently obtained from 3D smoothed-particle-hydrodynamics simulations. We conclude that the UPE, the hard X-ray flare and the soft X-ray flare do not form a causally connected sequence: Within our model they are the manifestation of \textbf{the same} physical process of the BH formation as seen through different viewing angles, implied by the morphology and the $\sim 300$~s rotation period of the HN ejecta.

We analyze the early X-ray flares in the GRB "flare-plateau-afterglow" (FPA) phase observed by Swift-XRT. The FPA occurs only in one of the seven GRB subclasses: the binary-driven hypernovae (BdHNe). This subclass consists of long GRBs with a carbon-oxygen core and a neutron star (NS) binary companion as progenitors. The hypercritical accretion of the supernova (SN) ejecta onto the NS can lead to the gravitational collapse of the NS into a black hole. Consequently, one can observe a GRB emission with isotropic energy $E_{iso}\gtrsim10^{52}$~erg, as well as the associated GeV emission and the FPA phase. Previous work had shown that gamma-ray spikes in the prompt emission occur at $\sim 10^{15}$--$10^{17}$~cm with Lorentz gamma factor $\Gamma\sim10^{2}$--$10^{3}$. Using a novel data analysis we show that the time of occurrence, duration, luminosity and total energy of the X-ray flares correlate with $E_{iso}$. A crucial feature is the observation of thermal emission in the X-ray flares that we show occurs at radii $\sim10^{12}$~cm with $\Gamma\lesssim 4$. These model independent observations cannot be explained by the "fireball" model, which postulates synchrotron and inverse Compton radiation from a single ultra relativistic jetted emission extending from the prompt to the late afterglow and GeV emission phases. We show that in BdHNe a collision between the GRB and the SN ejecta occurs at $\simeq10^{10}$~cm reaching transparency at $\sim10^{12}$~cm with $\Gamma\lesssim4$. The agreement between the thermal emission observations and these theoretically derived values validates our model and opens the possibility of testing each BdHN episode with the corresponding Lorentz gamma factor.