We prescribe a method to study the effects of self-gravity of accretion disk around a black hole associated with long Gamma Ray Bursts (GRBs) in an evolving background Kerr metric. This is an extension to our previous work where we presented possible constraints for the final masses and spins of these astrophysical black holes. Incorporating the self-force of the accreting cloud around the black hole is a very important aspect due to the transient nature of the event, in which a huge amount of mass is accreted and changes the fundamental black hole parameters i.e. its mass and spin, during the process. Understanding of the GRBs engine is important because they are possible sources of high-energy particles and gravitational waves as most of the energy released from the dynamical evolution is in the form of gravitational radiation. Here, we describe the analytical framework we developed to employ in our numerical model. The numerical studies are planned for the future work.

We present the relativistic 2D and 3D GRMHD simulation of axisymmetric, inviscid, hydrodynamic accretion flows in a fixed Kerr black hole gravitational field. The flow is having low angular momentum with respect to Keplerian one. A relativistic fluid where its bulk velocity is comparable to the speed of light, flowing in the accretion disk very close to the horizon should be described by adiabtic index: 4=3 < g < 5=3 .The time dependent evolution of shock position and respective effect on mass accretion rate and oscillation frequency with varying adiabatic index has been studied. Here we present some of the results for adiabatic index = 1.4 in a 2D and 3D model.

We compute the evolution of a quasi-spherical, slowly rotating accretion flow around a black hole, whose mass and spin evolve adequately to the mass-energy transfer through the horizon. Our model is relevant for the central engine driving a long gamma ray burst, that originates from the collapse of a massive star. Our results show how much mass and spin a newly formed black hole should possess during collapsar to launch long GRB.

We study the role of adiabatic index in determining the critical points in the transonic low angular momentum accretion flow onto a black hole. We present the general relativistic 2D hydrodynamic simulations of axisymmetric, inviscid accretion flows in a fixed Kerr black hole gravitational field. A relativistic fluid where its bulk velocity is comparable to the speed of light, flowing in the accretion disk very close to the horizon can be described by an adiabatic index of 4/3 < {\gamma} < 5/3. The time dependent evolution of the shock position and respective effects on mass accretion rate and oscillation frequency with varying adiabatic index is discussed in the context of the observed microquasars.

We compute the evolution of a quasi-spherical, slowly rotating accretion flow around a black hole, whose mass and spin evolve adequately to the mass-energy transfer through the horizon. Our model is relevant for the central engine driving a long gamma ray burst, that originates from the collapse of a massive star. The computations of a GRB engine in a dynamically evolving spacetime metric are important specifically due to the transient nature of the event, in which a huge amount of mass is accreted and changes the fundamental black hole parameters, its mass and spin, during the process. We discuss the results in the context of angular momentum magnitude of the collapsing star. We also study the possible formation and evolution of shocks in the envelope, which may temporarily affect accretion. Our results are important for the limitations on the mass and spin range of black holes detected independently by electromagnetic observations of GRBs and gravitational waves. We speculate on the possible constraints for the final masses and spins of these astrophysical black holes. It is shown that the most massive BHs were rahter not formed in a powerful GRB explosion, if the cores of their progenitors were only weakly rotating.

We study the structure and evolution of the hyperaccreting disks and outflows in the gamma ray bursts central engines. The torus around a stellar mass black hole is composed of free nucleons, Helium, electron-positron pairs, and is cooled by neutrino emission. Accretion of matter powers the relativistic jets, responsible for the gamma ray prompt emission. The significant number density of neutrons in the disk and outflowing material will cause subsequent formation of heavier nuclei. We study the process of nucleosynthesis and its possible observational consequences. We also apply our scenario to the recent observation of the gravitational wave signal, detected on September 14th, 2015 by the two Advanced LIGO detectors, and related to an inspiral and merger of a binary black hole system. A gamma ray burst that could possibly be related with the GW150914 event was observed by the Fermi satellite. It had a duration of about 1 second and appeared about 0.4 seconds after the gravitational-wave signal. We propose that a collapsing massive star and a black hole in a close binary could lead to the event. The gamma ray burst was powered by a weak neutrino flux produced in the star remnant's matter. Low spin and kick velocity of the merged black hole are reproduced in our simulations. Coincident gravitational-wave emission originates from the merger of the collapsed core and the companion black hole.

The microquasar IGR J17091-3624, which is the recently discovered analogue of the well known source GRS 1915+105, exhibits quasi-periodic outbursts, with a period of 5-70 seconds, and regular amplitudes, referred to as "heartbeat state". We argue that these states are plausibly explained by accretion disk instability, driven by the dominant radiation pressure. Using our GLobal Accretion DIsk Simulation hydrodynamical code, we model these outbursts quantitatively. We also find a correlation between the presence of massive outflows launched from the accretion disk and the stabilization of its oscillations. We verify the theoretical predictions with the available timing and spectral observations. Furthermore, we postulate that the underlying non-linear differential equations that govern the evolution of an accretion disk are responsible for the variability pattern of several other microquasars, including XTE J1550-564, GX 339-4, and GRO J1655-40. This is based on the signatures of deterministic chaos in the observed lightcurves of these sources, which we found using the recurrence analysis method. We discuss these results in the frame of the accretion disk instability model.

Data from the Fermi Gamma-ray Burst Monitor satellite observatory suggested that the recently discovered gravitational wave source, a pair of two coalescing black holes, was related to a gamma-ray burst. The observed high-energy electromagnetic radiation (above 50 keV) originated from a weak transient source and lasted for about 1 second. Its localization is consistent with the direction to GW150914. We speculate about the possible scenario for the formation of a gamma-ray burst accompanied by the gravitational-wave signal. Our model invokes a tight binary system consisting of a massive star and a black hole which leads to the triggering of a collapse of the star's nucleus, the formation of a second black hole, and finally to the binary black hole merger. For the most-likely configuration of the binary spin vectors with respect to the orbital angular momentum in the GW150914 event, the recoil speed (kick velocity) acquired by the final black hole through gravitational wave emission is of the order of a few hundred km/s and this might be sufficient to get it closer to the envelope of surrounding material and capture a small fraction of matter from the remnant of the host star. The gamma-ray burst is produced by the accretion of this remnant matter onto the final black hole. The moderate spin of the final black hole suggests that the gamma-ray burst jet is powered by weak neutrino emission rather than the Blandford-Znajek mechanism, and hence explains the low power available for the observed GRB signal.

Hardly any of the observed black hole accretion disks in X-Ray binaries and active galaxies shows constant flux. When the local stochastic variations of the disk occur at specific regions where a resonant behaviour takes place, there appear the Quasi-Periodic Oscillations (QPOs). If the global structure of the flow and its non-linear hydrodynamics affects the fluctuations, the variability is chaotic in the sense of deterministic chaos. Our aim is to solve a problem of the stochastic versus deterministic nature of the black hole binaries vari- ability. We use both observational and analytic methods. We use the recurrence analysis and we study the occurence of long diagonal lines in the recurrence plot of observed data series and compare it to the sur- rogate series. We analyze here the data of two X-Ray binaries - XTE J1550-564, and GX 339-4 observed by Rossi X-ray Timing Explorer. In these sources, the non-linear variability is expected because of the global conditions (such as the mean accretion rate) leading to the pos- sible instability of an accretion disk. The thermal-viscous instability and fluctuations around the fixed-point solution occurs at high accretion rate, when the radiation pressure gives dominant contribution to the stress tensor.

The dynamical system studied in previous papers of this series, namely a bound time-like geodesic motion in the exact static and axially symmetric space-time of an (originally) Schwarzschild black hole surrounded by a thin disc or ring, is considered to test whether the often employed "pseudo-Newtonian" approach (resorting to Newtonian dynamics in gravitational potentials modified to mimic the black-hole field) can reproduce phase-space properties observed in the relativistic treatment. By plotting Poincar\'e surfaces of section and using two recurrence methods for similar situations as in the relativistic case, we find similar tendencies in the evolution of the phase portrait with parameters (mainly with mass of the disc/ring and with energy of the orbiters), namely those characteristic to weakly non-integrable systems. More specifically, this is true for the Paczy\'nski--Wiita and a newly suggested logarithmic potential, whereas the Nowak--Wagoner potential leads to a different picture. The potentials and the exact relativistic system clearly differ in delimitation of the phase-space domain accessible to a given set of particles, though this mainly affects the chaotic sea whereas not so much the occurrence and succession of discrete dynamical features (resonances). In the pseudo-Newtonian systems, the particular dynamical features generally occur for slightly smaller values of the perturbation parameters than in the relativistic system, so one may say that the pseudo-Newtonian systems are slightly more prone to instability. We also add remarks on numerics (a different code is used than in previous papers), on the resemblance of dependence of the dynamics on perturbing mass and on orbital energy, on the difference between the Newtonian and relativistic Bach--Weyl rings, and on the relation between Poincar\'e sections and orbital shapes within the meridional plane.