Accretion of a massive magnetized torus on a rotating black hole

We present numerical simulations of the axisymmetric accretion of a massive magnetized plasma torus on a rotating black hole. We use a realistic equation of state, which takes into account neutrino cooling and energy loss due to nucleus dissociations. The calculation are performed in the ideal relativistic MHD approximation using an upwind conservative scheme that is based on a linear Riemann solver and the constrained transport method to evolve the magnetic field. The gravitational attraction of the black hole is introduced via the Kerr metric in the Kerr–Schild coordinates. We simulate various magnetic field configurations and torus models, both optically thick and thin for neutrinos.We have found an effect of alternation of the magnetic field orientation in the ultrarelativistic jet formed as a result of the collapse. The calculations show evidence for heating of the wind surrounding the collapsar by the shock waves generated at the jet–wind border. It is shown that the neutrino cooling does not significantly change either the structure of the accretion flow or the total energy release of the system. The angular momentum of the accreting matter defines the time scale of the accretion. Due to the absence of the magnetic dynamo in our calculations, the initial strength and topology of the magnetic field determines the magnetization of the black hole, jet formation properties and the total energy yield. We estimate the total energy of accretion which transformed to jets as 1.3 × 1052 ergs which was sufficient to explain hypernova explosions like GRB 980425 or GRB 030329.