Aleksander Sądowski

Numerical simulations of super-critical black hole accretion flows in general relativity

A new general relativistic radiation magnetohydrodynamical code KORAL, is described, which employs the M1 scheme to close the radiation moment equations. The code has been successfully verified against a number of tests. Axisymmetric simulations of super-critical magnetized accretion on a non-rotating black hole (a=0.0) and a spinning black hole (a=0.9) are presented. The accretion rates in the two models are \dot M = 100-200 \dot M_Edd. These first general relativistic simulations of super-critical black hole accretion are potentially relevant to tidal disruption events and hyper-accreting supermassive black holes in the early universe. Both simulated models are optically and geometrically thick, and have funnels through which energy escapes in the form of relativistic gas, Poynting flux and radiative flux. The jet is significantly more powerful in the a=0.9 run. The net energy outflow rate in the two runs correspond to efficiencies of 5% (a=0) and 33% (a=0.9), as measured with respect to the mass accretion rate at the black hole. These efficiencies agree well with those measured in previous simulations of non-radiative geometrically thick disks. Furthermore, in the a=0.9 run, the outflow power appears to originate in the spinning black hole, suggesting that the associated physics is again similar in non-radiative and super-critical accretion flows. While the two simulations are efficient in terms of total energy outflow, both runs are radiatively inefficient. Their luminosities are only \sim 1-10 L_Edd, which corresponds to a radiative efficiency \sim 0.1%. Interestingly, most of the radiative luminosity emerges through the funnels, which subtend a very small solid angle. Therefore, measured in terms of a local radiative flux, the emitted radiation is highly super-Eddington.

GRMHD simulations of magnetized advection‐dominated accretion on a non‐spinning black hole: role of outflows

We present results from two long-duration general relativistic magneto-hydrodynamic (GRMHD) simulations of advection-dominated accretion around a non-spinning black hole. The first simulation was designed to avoid significant accumulation of magnetic flux around the black hole. This simulation was run for a time of 200 000 GM/c3 and achieved inflow equilibrium out to a radius ∼90 GM/c2. Even at this relatively large radius, the mass outflow rate is found to be only 60 per cent of the net mass inflow rate into the black hole. The second simulation was designed to achieve substantial magnetic flux accumulation around the black hole in a magnetically arrested disc. This simulation was run for a shorter time of 100 000 GM/c3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius ∼170 GM/c2. Here, becomes equal to at r ∼ 160 GM/c2. Since the mass outflow rates in the two simulations do not show robust convergence with time, it is likely that the true outflow rates are lower than our estimates. The effect of black hole spin on mass outflow remains to be explored. Neither simulation shows strong evidence for convection, though a complete analysis including the effect of magnetic fields is left for the future.

Global simulations of axisymmetric radiative black hole accretion discs in general relativity with a mean-field magnetic dynamo

We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency,

SLIM DISKS AROUND KERR BLACK HOLES REVISITED

\eta\approx0.04

Measuring black hole spin by the continuum-fitting method: effect of deviations from the Novikov–Thorne disc model

, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates (

Energy, momentum and mass outflows and feedback from thick accretion discs around rotating black holes

\dot{M} \lesssim 2\dot M_{\rm Edd}

Powerful radiative jets in supercritical accretion discs around non-spinning black holes

). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.

Relativistic slim disks with vertical structure

We investigate stationary slim accretion disks around Kerr black holes. We construct a new numerical method based on the relaxation technique. We systematically cover the whole parameter space relevant to stellar mass X-ray binaries. We also notice some non-monotonic features in the disk structure, overlooked in previous studies.

The Shakura-Sunyaev viscosity prescription with variable α (r)

The X-ray spectra of accretion discs of eight stellar mass black holes have been analysed to date using the thermal continuum-fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum-fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01, respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4–0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov–Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

Leaving the innermost stable circular orbit: the inner edge of a black-hole accretion disk at various luminosities

Using long-duration general relativistic magnetohydrodynamic simulations of radiatively inefficient accretion discs, the energy, momentum and mass outflow rates from such systems are estimated. Outflows occur via two fairly distinct modes: a relativistic jet and a sub-relativistic wind. The jet power depends strongly on the black hole spin and on the magnetic flux at the horizon. Unless these are very small, the energy output in the jet dominates over that in the wind. For a rapidly spinning black hole accreting in the magnetically arrested limit, it is confirmed that jet power exceeds the total rate of accretion of rest mass energy. However, because of strong collimation, the jet probably does not have a significant feedback effect on its immediate surroundings. The power in the wind is more modest and shows a weaker dependence on black hole spin and magnetic flux. Nevertheless, because the wind subtends a large solid angle, it is expected to provide efficient feedback on a wide range of scales inside the host galaxy. Empirical formulae are obtained for the energy and momentum outflow rates in the jet and the wind.