Person:

Sadowski, Aleksander

We explore the variability properties of long, high cadence GRMHD simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both disk- and jet-dominated simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A∗ and the size of its image at 1.3 mm. We find that the disk-dominated models produce short timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, jet-dominated models generate only slow variability, at lower flux levels. Neither set of models show any X-ray flares, which most likely indicate that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The disk-dominated models show strong, short-lived mm/IR flares, with short (. 1 hr) time lags between the mm and IR wavelengths, that arise from strong-field gravitational lensing of magnetic flux tubes near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

Slim-disk models describe advective accretion flows at high luminosities, while reducing to the standard thin disk form in the low luminosity limit. We have developed a new spectral model, slimbb, within the framework of XSPEC, which describes fully relativistic slim-disk accretion and includes photon ray-tracing that starts from the disk photosphere, rather than the equatorial plane. We demonstrate the features of this model by applying it to RXTE spectra of the persistent black-hole X-ray binary LMC X-3. LMC X-3 has the virtues of exhibiting large intensity variations while maintaining itself in soft spectral states which are well described using accretion-disk models, making it an ideal candidate to test the aptness of slimbb. Our results demonstrate consistency between the low-luminosity (thin-disk) and high luminosity (slim-disk) regimes. The results also illustrate that advection alone does not solve the problem of the origin of the surprisingly soft high-luminosity spectra in LMC X-3. We show that X-ray continuum-fitting in the high accretion rate regime can powerfully test black-hole accretion disk models.

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.

We present results from two long-duration 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,000GM/c^3 and achieved inflow equilibrium out to a radius \sim90GM/c^2. Even at this relatively large radius, the mass outflow rate \dot{M}{out} is found to be only 60% of the net mass inflow rate \dot{M}{BH} 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,000GM/c^3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius \sim170GM/c^2. Here, \dot{M}{out} becomes equal to \dot{M}{BH} at r\sim 160GM/c^2. 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.

We present a set of four three-dimensional, general relativistic, radiation magnetohydrodynamical simulations of black hole accretion at supercritical mass accretion rates, dot{M} > dot{M}_Edd. We use these simulations to study how disc properties are modified when we vary the black hole mass, the black hole spin, or the mass accretion rate. In the case of a non-rotating black hole, we find that the total efficiency is of the order of 3 per cent dot{M} c^2, approximately a factor of 2 less than the efficiency of a standard thin accretion disc. The radiation flux in the funnel along the axis is highly super-Eddington, but only a small fraction of the energy released by accretion escapes in this region. The bulk of the 3 per cent dot{M} c^2 of energy emerges farther out in the disc, either in the form of photospheric emission or as a wind. In the case of a black hole with a spin parameter of 0.7, we find a larger efficiency of about 8 per cent dot{M} c^2. By comparing the relative importance of advective and diffusive radiation transport, we show that photon trapping is effective near the equatorial plane. However, near the disc surface, vertical transport of radiation by diffusion dominates. We compare the properties of our fiducial three-dimensional run with those of an equivalent two-dimensional axisymmetric model with a mean-field dynamo. The latter simulation runs nearly 100 times faster than the three-dimensional simulation, and gives very similar results for time-averaged properties of the accretion flow, but does not reproduce the time-variability.

We study energy flows in geometrically thick accretion discs, both optically thick and thin, using general relativistic, three-dimensional simulations of black hole accretion flows. We find that for non-rotating black holes the efficiency of the total feedback from thick accretion discs is 3% - roughly half of the thin disc efficiency. This amount of energy is ultimately distributed between outflow and radiation, the latter scaling weakly with the accretion rate for supercritical accretion rates, and returned to the interstellar medium. Accretion on to rotating black holes is more efficient because of the additional extraction of rotational energy. However, the jet component is collimated and likely to interact only weakly with the environment, whereas the outflow and radiation components cover a wide solid angle.

In general relativity static gaseous atmospheres may be in hydrostatic balance in the absence of a supporting stellar surface, provided that the luminosity is close to the Eddington value. We construct analytic models of optically thin, spherically symmetric shells supported by the radiation pressure of a luminous central body in the Schwarzschild metric. Opacity is assumed to be dominated by Thomson scattering. The inner parts of the atmospheres, where the luminosity locally has supercritical values, are characterized by a density and pressure inversion. The atmospheres are convectively and Rayleigh-Taylor stable, and there is no outflow of gas.

We construct models of static, spherically symmetric shells supported by the radiation flux of a luminous neutron star in the Schwarzschild metric. The atmospheres are disconnected from the star and levitate above its surface. Gas pressure and density inversion appear in the inner region of these atmospheres, which is a purely relativistic phenomenon. We account for the scattering opacity dependence on temperature green by using the Klein-Nishina formula. The relativistic M1 closure scheme for the radiation tensor provides a general relativity-consistent treatment of the photon flux and radiation tensor anisotropy. In this way, we are able to address atmospheres of both large and moderate/low optical depths with the same set of equations. We discuss properties of the levitating atmospheres and find that they may indeed be optically thick, with the distance between star surface and the photosphere expanding as luminosity increases. These results may be relevant for the photosphereric radius expansion X-ray bursts.

We describe a set of simulations of super-critical accretion onto a non-rotating supermassive BH. The accretion flow is radiation pressure dominated and takes the form of a geometrically thick disk with twin low-density funnels around the rotation axis. For accretion rates & 10M˙ Edd, there is sufficient gas in the funnel to make this region optically thick. Radiation from the disk first flows into the funnel, after which it accelerates the optically thick funnel gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted velocity ≈ 0.3c. Much of the radiative luminosity is converted into kinetic energy by the time the escaping gas becomes optically thin. These jets are not powered by black hole rotation or magnetic driving, but purely by radiation. Their characteristic beaming angle is ∼ 0.2 radians. For an observer viewing down the axis, the isotropic equivalent luminosity of total energy is as much as 1048 erg s−1 for a 107M BH accreting at 103 Eddington. Therefore, energetically, the simulated jets are consistent with observations of the most powerful tidal disruption events, e.g., Swift J1644. The jet velocity is, however, too low to match the Lorentz factor γ > 2 inferred in J1644. There is no such conflict in the case of other tidal disruption events. Since favorably oriented observers see isotropic equivalent luminosities that are highly superEddington, the simulated models can explain observations of ultra-luminous X-ray sources, at least in terms of luminosity and energetics, without requiring intermediate mass black holes. The spectrum remains to be worked out. Finally, since the simulated jets are baryon-loaded and have mildly relativistic velocities, they match well the jets observed in SS433. The latter are, however, more collimated than the simulated jets. This suggests that, even if magnetic fields are not important for acceleration, they may perhaps still play a role in confining the jet.