Equations of general relativistic radiation hydrodynamics in Kerr space–time

Equations of fully general relativistic radiation hydrodynamics in Kerr space–time are derived. While the interactions between matter and radiation are introduced in the comoving frame, the derivatives used when describing the global evolutions of both the matter and the radiation are given in the Boyer–Lindquist frame (BLF) which is a frame fixed to the coordinate describing the central black hole. Around a rotating black hole, both the matter and the radiation are influenced by the frame-dragging effects due to the black hole's rotation. As a fixed frame, we use the locally non-rotating reference frame (LNRF) which is one of the orthonormal frame. While the special relativistic effects such as beaming effects are introduced by the Lorentz transformation connecting the comoving frame and the LNRF, the general relativistic effects such as frame dragging and gravitational redshift are introduced by the tetrads connecting the LNRF and the BLF.     

On the interpretation of the multicolour disc model for black hole candidates

We present a critical analysis of the usual interpretation of the multicolour disc model parameters for black hole candidates in terms of the inner radius and temperature of the accretion disc. Using a self-consistent model for the radiative transfer and the vertical temperature structure in a Shakura–Sunyaev disc, we simulate the observed disc spectra, taking into account Doppler blurring and gravitational redshift, and fit them with multicolour models. We show not only that such a model systematically underestimates the value of the inner-disc radius, but that when the accretion rate and/or the energy dissipated in the corona are allowed to change, the inner edge of the disc, as inferred from the multicolour model, appears to move even when it is in fact fixed at the innermost stable orbit.     

On some transonic aspects of general relativistic spherical accretion on to Schwarzschild black holes

The equations governing general relativistic, spherically symmetric, hydrodynamic accretion of polytropic fluid on to black holes are solved in the Schwarzschild metric to investigate some of the transonic properties of the flow. Only stationary solutions are discussed. For such accretion, it has been shown that real physical sonic points may form even for flow with   γ <4/3  or   γ >5/3  . The behaviour of some flow variables in the close vicinity of the event horizon is studied as a function of specific energy and the polytropic index of the flow.     

Effects of Kerr space–time on spectral features from X-ray illuminated accretion discs

We performed detailed calculations of the relativistic effects acting on both the reflection continuum and the iron line from accretion discs around rotating black holes. Fully relativistic transfer of both illuminating and reprocessed photons has been considered in Kerr space–time. We calculated overall spectra, line profiles and integral quantities, and present their dependences on the black hole angular momentum. We show that the observed EW of the lines is substantially enlarged when the black hole rotates rapidly and/or the source of illumination is near above the hole. Therefore, such calculations provide a way to distinguish between different models of the central source.     

For whom the disc tolls

Report on the Nordita Workdays on Quasi‐Peridic Oscillations (QPOs). (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Critical properties of spherically symmetric black hole accretion in Schwarzschild geometry

The stationary, spherically symmetric, polytropic and inviscid accretion flow in the Schwarzschild metric has been set-up as an autonomous first-order dynamical system, and it has been studied completely analytically. Of the three possible critical points in the flow, the one that is physically realistic behaves like the saddle point of the standard Bondi accretion problem. One of the two remaining critical points exhibits the strange mathematical behaviour of being either a saddle point or a centre-type point, depending on the values of the flow parameters. The third critical point is always unphysical and behaves like a centre-type point. The treatment has been extended to pseudo-Schwarzschild flows for comparison with the general relativistic analysis.     

Models for the X‐ray spectra and variability of luminous accreting black holes

The X‐ray spectra of luminous Seyfert 1 galaxies often appear to be reflection dominated. In a number of Narrow Line Seyfert 1 (NLS1) galaxies and galactic black holes in the very high state, the variability of the continuum and of the iron line are decoupled, the reflected component being often much less variable than the continuum. These properties have been interpreted as effects of gravitational light bending. In this framework, we present detailed Monte‐Carlo simulations of the reflection continuum in the Kerr metric. These calculations confirm that the spectra and variability behaviour of these sources can be reproduced by the light bending model. As an alternative to the light bending model, we show that similar observational properties are expected from radiation pressure dominated discs subject to violent clumping instabilities and, as a result, have a highly inhomogeneous two‐phase structure. In this model, most of the observed spectral and variability features originate from the complex geometrical structure of the inner regions of near‐Eddington accretion flows and are therefore a signature of accretion physics rather than general relativity. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-axisymmetric relativistic Bondi–Hoyle accretion on to a Schwarzschild black hole

We present the results of an exhaustive numerical study of fully relativistic non-axisymmetric Bondi–Hoyle accretion on to a moving Schwarzschild black hole. We have solved the equations of general relativistic hydrodynamics with a high-resolution shock-capturing numerical scheme based on a linearized Riemann solver. The numerical code was previously used to study axisymmetric flow configurations past a Schwarzschild black hole. We have analysed and discussed the flow morphology for a sample of asymptotically high Mach number models. The results of this work reveal that initially asymptotic uniform flows always accrete on to the hole in a stationary way, which closely resembles the previous axisymmetric patterns. This is in contrast with some Newtonian numerical studies where violent flip-flop instabilities were found. As discussed in the text, the reason can be found in the initial conditions used in the relativistic regime, as they cannot exactly duplicate the previous Newtonian setups where the instability appeared. The dependence of the final solution on the inner boundary condition as well as on the grid resolution has also been studied. Finally, we have computed the accretion rates of mass and linear and angular momentum.     

Spectral variability in transonic discs around black holes

Transonic discs with accretion rates relevant to intrinsically bright Galactic X-ray sources ( L ≈10³⁸–10³⁹ erg s⁻¹) exhibit a time-dependent cyclic behaviour due to the onset of a thermal instability driven by radiation pressure. In this paper we calculate radiation spectra emitted from thermally unstable discs to provide detailed theoretical predictions for observationally relevant quantities. The emergent spectrum has been obtained by solving self-consistently the vertical structure and radiative transfer in the disc atmosphere. We focus on four particular stages of the disc evolution, the maximal evacuation stage and three intermediate stages during the replenishment phase. The disc is found to undergo rather dramatic spectral changes during the evolution, emitting mainly in the 1–10 keV band during outburst and in the 0.1–1 keV band off-outburst. Local spectra, although different in shape from a blackbody at the disc effective temperature, may be characterized in terms of a hardening factor f . We have found that f is more or less constant, both in radius and in time, with a typical value ∼ 1.65.     

Extreme gravitational lensing near rotating black holes

We describe a new approach to calculating photon trajectories and gravitational lensing effects in the strong gravitational field of the Kerr black hole. These techniques are applied to explore both the imaging and spectral properties of photons emitted from an accretion disc, which perform multiple orbits of the central mass before escaping to infinity. Viewed at large inclinations, these higher-order photons contribute ∼20 per cent of the total luminosity of the system for a Schwarzschild hole, while for an extreme Kerr black hole this fraction rises to ∼60 per cent. In more realistic models, these photons will be reabsorbed by the disc at large distances from the hole, but this returning radiation could provide a physical mechanism to resolve the discrepancy between the predicted and observed optical/ultraviolet colours in active galactic nuclei. Conversely, at low inclinations, higher-order images reintercept the disc plane close to the black hole, so need not be absorbed by the disc if this is within the plunging region. These photons form a bright ring carrying approximately 10 per cent of the total disc luminosity for a Schwarzschild black hole. The spatial separation between the inner edge of the disc and the ring is similar to the size of the event horizon. This is resolvable for supermassive black holes with proposed X-ray interferometery missions such as the Microarcsecond X-ray Imaging Mission (MAXIM), and so has the potential to provide an observational test of strong field gravity.     

Self-consistent computation of γ-ray spectra due to proton–proton interactions in black hole systems

In the inner regions of an accretion disc around a black hole, relativistic protons can interact with ambient matter to produce electrons, positrons and γ-rays. The resultant steady-state electron and positron particle distributions are self-consistently computed taking into account Coulomb and Compton cooling,  e⁻ e⁺  pair production (due to γ–γ annihilation) and pair annihilation. While earlier works used the diffusion approximation to obtain the particle distributions, here we solve a more general integro-differential equation that correctly takes into account the large change in particle energy that occurs when the leptons Compton scatter off hard X-rays. Thus this formalism can also be applied to the hard state of black hole systems, where the dominant ambient photons are hard X-rays. The corresponding photon energy spectrum is calculated and compared with broad-band data of black hole binaries in different spectral states. The results indicate that the γ-ray spectra  ( E > 0.8 MeV)  of both the soft and hard spectral states and the entire hard X-ray/γ-ray spectrum of the ultrasoft state could be due to p–p interactions. These results are consistent with the hypothesis that there always exists in these systems a γ-ray spectral component due to p–p interactions that can contribute between 0.5 and 10 per cent of the total bolometric luminosity. The model predicts that GLAST would be able to detect black hole binaries and provide evidence for the presence of non-thermal protons, which in turn would give insight into the energy dissipation process and jet formation in these systems.