X-rays from cusps of compact remnants near galactic centres

Compact remnants – stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low-mass stars, forming a high concentration cusp. Same physical region may also contain very high-density molecular clouds and accretion discs that are needed to fuel supermassive black hole (SMBH) activity. Here we estimate gas capture rates on to the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This 'diffuse' emission may be of importance for local moderately bright active galactic nuclei (AGNs), especially low-luminosity AGNs. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.     

Models for jet power in elliptical galaxies: a case for rapidly spinning black holes

The power of jets from black holes is expected to depend on both the spin of the black hole and the structure of the accretion disc in the region of the last stable orbit. We investigate these dependencies using two different physical models for the jet power: the classical Blandford–Znajek (BZ) model and a hybrid model developed by Meier. In the BZ case, the jets are powered by magnetic fields directly threading the spinning black hole while in the hybrid model, the jet energy is extracted from both the accretion disc as well as the black hole via magnetic fields anchored to the accretion flow inside and outside the hole's ergosphere. The hybrid model takes advantage of the strengths of both the Blandford–Payne and BZ mechanisms, while avoiding the more controversial features of the latter. We develop these models more fully to account for general relativistic effects and to focus on advection-dominated accretion flows (ADAFs) for which the jet power is expected to be a significant fraction of the accreted rest mass energy.
We apply the models to elliptical galaxies, in order to see if these models can explain the observed correlation between the Bondi accretion rates and the total jet powers. For typical values of the disc viscosity parameter  α∼ 0.04 –0.3  and mass accretion rates consistent with ADAF model expectations, we find that the observed correlation requires   j ≳ 0.9  ; that is, it implies that the black holes are rapidly spinning. Our results suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive jets powerful enough to heat the intracluster medium and quench cooling flows.     

Gravitational wave emission from a companion black hole in the presence of an accretion disc around a super-massive Kerr black hole

Gravitational wave signal characteristics from a binary black hole system in which the companion moves through the accretion disc of the primary are studied. We chose the primary to be a super-massive  ( M = 10⁸ M_⊙)  Kerr black hole and the companion to be a massive black hole  ( M = 10⁵ M_⊙)  to clearly demonstrate the effects. We show that the drag exerted on the companion by the disc is sufficient to reduce the coalescence time of the binary. The drag is primarily due to the fact that the accretion disc on a black hole deviates from a Keplerian disc and becomes sub-Keplerian due to inner boundary condition on the black hole horizon. We consider two types of accretion rates on to the companion. The companion is deeply immersed inside the disc and it can accrete at the Bondi rate which depends on the instantaneous density of the disc. However, an accretion disc can also form around the smaller black hole and it can accrete at its Eddington rate. Thus, this case is also studied and the results are compared. We find that the effect of the disc will be significant in reducing the coalescence time and one needs to incorporate this while interpreting gravitational wave signals emitted from such a binary system.     

原太阳吸积盘结构

计算了粘滞演化阶段原太阳吸积盘结构。采用稳态标准吸积盘模型来描述盘中湍动粘滞；忽略其径向能量传输，将垂直结构作为一维问题处理。假设盘作Keplerian较差旋转，处于流体力学平衡和局域热平衡，盘由粘滞耗散加热，能量通过对流和辐射向外传输。结果表明，对温度敏感的不透明度是决定盘结构的重要因素；原太阳吸积盘为冷的薄盘，盘中热对流不稳定性由外而内，由上而下地终结；行星的形成应首先开始于对流终结的区域。     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.     

Alignment and precession of a black hole with a warped accretion disc

We consider the shape of an accretion disc whose outer regions are misaligned with the spin axis of a central black hole and calculate the steady state form of the warped disc in the case where the viscosity and surface densities are power laws in the distance from the central black hole. We discuss the shape of the resulting disc in both the frame of the black hole and that of the outer disc. We note that some parts of the disc and also any companion star maybe shadowed from the central regions by the warp. We compute the torque on the black hole caused by the Lense–Thirring precession, and hence compute the alignment and precession time-scales. We generalize the case with viscosity and hence surface density independent of radius to more realistic density distributions for which the surface density is a decreasing function of radius. We find that the alignment time-scale does not change greatly but the precession time-scale is more sensitive. We also determine the effect on this time-scale if we truncate the disc. For a given truncation radius, the time-scales are less affected for more sharply falling density distributions.     

The jet power extracted from a magnetized accretion disc

We consider the power of a relativistic jet accelerated by the magnetic field of an accretion disc. It is found that the power extracted from the disc is mainly determined by the field strength and configuration of the field far from the disc. Comparing it with the power extracted from a rotating black hole, we find that the jet power extracted from a disc can dominate over that from the rotating black hole. However, in some cases, the jet power extracted from a rapidly rotating hole can be more important than that from the disc, even if the poloidal field threading the hole is not significantly larger than that threading the inner edge of the disc. The results imply that the radio-loudness of quasars may be governed by its accretion rate, which might be regulated by the central black hole mass. It is proposed that the different disc field generation mechanisms might be tested against observations of radio-loud quasars if their black hole masses are available.     

Super-Eddington accretion discs around Kerr black holes

We calculate the structure of the accretion disc around a rapidly rotating black hole with a super-Eddington accretion rate. The luminosity and height of the disc are reduced by the advection effect. In the case of large viscosity parameter, α>0.03, the accretion flow deviates strongly from thermodynamic equilibrium and overheats in the central region. With increasing accretion rate, the flow temperature steeply increases, reaches maximum, and then falls off. The maximum is achieved in the advection-dominated regime of accretion. The maximum temperature in the disc around a massive black hole of M =10⁸ M⊙ with α=0.3 is of order 3×10⁸ K. The discs with large accretion rates can emit X-rays in quasars as well as in galactic black hole candidates.     

Hot accretion discs with thermal Comptonization and advection in luminous black hole sources

We solve for the structure of a hot accretion disc with unsaturated thermal Comptonization of soft photons and with advection, generalizing the classical model of Shapiro et al. The upper limit on the accretion rate due to advection constrains the luminosity to ≲ 0.15 y^3/5 α^7/5 of the Eddington limit, where y and α are the Compton and viscosity parameters, respectively. The characteristic electron temperature and Thomson optical depth of the inner flow at accretion rates within an order of magnitude of that upper limit are ∼ 10⁹ K and ∼ 1, respectively. The resulting spectra are then in close agreement with the X-ray and soft γ-ray spectra from black hole binaries in the hard state and Seyferts. At low accretion rates, bremsstrahlung becomes the dominant radiative process.     

Iron line profiles and self-shadowing from relativistic thick accretion discs

We present Fe Kα line profiles from and images of relativistic discs with finite thickness around a rotating black hole using a novel code. The line is thought to be produced by iron fluorescence of a relatively cold X-ray-illuminated material in the innermost parts of the accretion disc and provides an excellent diagnostic of accretion flows in the vicinity of black holes. Previous studies have concentrated on the case of a thin, Keplerian accretion disc. This disc must become thicker and sub-Keplerian with increasing accretion rates. These can affect the line profiles and in turn can influence the estimation of the accretion disc and black hole parameters from the observed line profiles. We here embark on, for the first time, a fully relativistic computation which offers key insights into the effects of geometrical thickness and the sub-Keplerian orbital velocity on the line profiles. We include all relativistic effects such as frame-dragging, Doppler boost, time dilation, gravitational redshift and light bending. We find that the separation and the relative height between the blue and red peaks of the line profile diminish as the thickness of the disc increases. This code is also well suited to produce accretion disc images. We calculate the redshift and flux images of the accretion disc and find that the observed image of the disc strongly depends on the inclination angle. The self-shadowing effect appears remarkable for a high inclination angle, and leads to the black hole shadow being in this case, completely hidden by the disc itself.     

A General Relativistic Calculation of Boundary Layer and Disc Luminosity for Accreting Non-magnetic Neutron Stars in Rapid Rotation

We calculate the disc and boundary layer luminosities for accreting rapidly rotating neutron stars with low magnetic fields in a fully general relativistic manner. Rotation increases the disc luminosity and decreases the boundary layer luminosity. A rapid rotation of the neutron star substantially modifies these quantities as compared with the static limit. For a neutron star rotating close to the centrifugal mass shed limit, the total luminosity has contribution only from the extended disc. For such maximal rotation rates, we find that well before the maximum stable gravitational mass configuration is reached, there exists a limiting central density, for which particles in the innermost stable orbit will be more tightly bound than those at the surface of the neutron star. We also calculate the angular velocity profiles of particles in Keplerian orbits around the rapidly rotating neutron star. The results are illustrated for a representative set of equation of state models of neutron star matter.     

Magnetohydrodynamical non-radiative accretion flows in two dimensions

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided that the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow with previous hydrodynamical simulations. The net mass accretion rate is small compared with the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared with the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.     

Evolution of protoplanetary discs driven by the MRI, self-gravity and hydrodynamical turbulence

We study the viscous evolution of protoplanetary discs driven by the combined action of magnetohydrodynamic turbulence, resulting from the magneto-rotational instability (MRI), self-gravity torques, parametrized in terms of an effective viscosity and an additional viscous agent of unspecified origin. The distribution of torques driving the evolution of the disc is calculated by analysing where in the disc the MRI develops and, to incorporate the effect of self-gravity, calculating the Toomre parameter. We find that, generally, discs rapidly evolve towards a configuration where the intermediate regions, from a fraction of an au to a few au, are stable against the MRI due to their low-ionization degree. As an additional source of viscosity is assumed to operate in those regions, subsequent evolution of the disc is eruptive. Brief episodes of high mass accretion ensue as the criterion for the development of the MRI is met in the low-ionization region. The radial distribution of mass and temperature in the disc differs considerably from disc models with constant α parameter or layered accretion models, with potentially important consequences on the process of planet formation.     

The effect of radiation stress on gas motions in AGN accretion discs

In this paper we look at one of the effects of irradiation on a warped accretion disc in the context of active galactic nuclei (AGN). A warp will catch a substantial amount of the radiation emitted by the central object. We consider the fluid motions that may arise inside a warped disc when the surface is subject to a radiation stress, and also the net mass flows that result. We find that, to first order, we have a balance of the viscous and Coriolis-type forces. The radial radiation stress causes outward motion of the surface layer, but only the azimuthal Poynting–Robertson drag leads to an increase in the net accretion rate. We investigate the distribution of the velocity perturbations and find them to be significant in determining the local structure of the disc.
An unexpected result is that the picture changes significantly when we take into account the periodic illumination of the warped disc. A type of resonance at the local Keplerian rotation frequency causes a flow that penetrates the whole thickness of the disc; these flows are faster than the flows due to unchanging illumination. They totally dominate the induced flows in terms of sheer mass, but significant impact on disc structure still occurs only near the surface, where velocity perturbations typically go up to some kilometres per second.     

Application of a relativistic accretion disc model to X-ray spectra of LMC X-1 and GRO J1655-40

We present a general relativistic accretion disc model and its application to the soft-state X-ray spectra of black hole binaries. The model assumes a flat, optically thick disc around a rotating Kerr black hole. The disc locally radiates away the dissipated energy as a blackbody. Special and general relativistic effects influencing photons emitted by the disc are taken into account. The emerging spectrum, as seen by a distant observer, is parametrized by the black hole mass and spin, the accretion rate, the disc inclination angle and the inner disc radius.
We fit the ASCA soft-state X-ray spectra of LMC X-1 and GRO J1655-40 by this model. We find that, having additional limits on the black hole mass and inclination angle from optical/UV observations, we can constrain the black hole spin from X-ray data. In LMC X-1 the constraint is weak, and we can only rule out the maximally rotating black hole. In GRO J1655-40 we can limit the spin much better, and we find 0.68 a 0.88 . Accretion discs in both sources are radiation-pressure dominated. We do not find Compton reflection features in the spectra of any of these objects.     

A massive thick disc around a massive black hole and its runaway instability

We investigate the runaway instability of configurations consisting of a massive dense but non-self-gravitating thick disc/torus surrounding a massive black hole (MBH). We limit our model parameters to values that result in a self-consistent thick disc around an MBH. We identify, analytically, the index of the angular momentum distribution that will form a thick disc as an initial equilibrium state, and obtain the mass ratio of the disc to the central black hole for which the disc is dominated by the radiation pressure. In our theoretical framework we find that a self-consistent thick disc with constant angular momentum leads to a runaway instability on a dynamical time-scale. However, even a slight increase of the specific angular momentum outwards has a strong stabilizing effect on the accretion process. Finally, we discuss our results and present possible applications to high-energy emission.     

Quasi-stars: accreting black holes inside massive envelopes

We study the structure and evolution of 'quasi-stars', accreting black holes embedded within massive hydrostatic gaseous envelopes. These configurations may model the early growth of supermassive black hole seeds. The accretion rate on to the black hole adjusts so that the luminosity carried by the convective envelope equals the Eddington limit for the total mass,   M _*+ M _BH≈ M _*  . This greatly exceeds the Eddington limit for the black hole mass alone, leading to rapid growth of the black hole. We use analytic models and numerical stellar structure calculations to study the structure and evolution of quasi-stars. We show that the photospheric temperature of the envelope scales as   T _ph∝ M ^−2/5_BH M ^7/20_*  , and decreases with time while the black hole mass increases. Once   T _ph < 10⁴ K  , the photospheric opacity drops precipitously and T _ph hits a limiting value, analogous to the Hayashi track for red giants and protostars, below which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) opacities, this limiting temperature is approximately 4000 K. After a quasi-star reaches this limiting temperature, it is rapidly dispersed by radiation pressure. We find that black hole seeds with masses between 10³ and  10⁴ M_⊙  could form via this mechanism in less than a few Myr.     

Black hole mergers: can gas discs solve the 'final parsec' problem?

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework.
The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high-photon-energy variability on the relatively long orbital time-scale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary.
We also find that the disc can drive the binary to merger on a reasonable time-scale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e.   a ₀≲ 0.05  pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.     

On equilibrium rotation of the central object of an accretion disc

The hydrodynamic interaction of an accretion disc with its central object is reanalysed within the framework of the slim-disc approximation. Arguments are presented against an interpretation of the total angular momentum flux as an eigenvalue of the system. A simple intuitive consideration is provided, which shows that the central object may be in a state of stationary rotation even if the disc imposes the constraint of a finite angular momentum flux into it. It is argued that equilibrium rotation is characterized by vanishing viscous torque rather than by zero total angular momentum flux. As a consequence, the central object can be in a state of stationary rotation below the break-up limit, although its angular momentum increases. Despite accretion, even for positive total angular momentum flux and subcritical rotation, central objects are spun down within a considerable range of their parameters. The results are illustrated by application to FU Orionis systems.     

Radio and X-ray properties of relativistic beaming models for ultraluminous X-ray sources

We calculate the broad-band radio–X-ray spectra predicted by microblazar and microquasar models for ultraluminous X-ray sources (ULXs), exploring the possibility that their dominant power-law component is produced by a relativistic jet, even at near-Eddington mass accretion rates. We do this by first constructing a generalized disc–jet theoretical framework in which some fraction of the total accretion power, P _a, is efficiently removed from the accretion disc by a magnetic torque responsible for jet formation. Thus, for different black hole masses, mass accretion rates and magnetic coupling strength, we self-consistently calculate the relative importance of the modified disc spectrum, as well as the overall jet emission due to synchrotron and Compton processes. In general, transferring accretion power to a jet makes the disc fainter and cooler than a standard disc at the same mass accretion rate; this may explain why the soft spectral component appears less prominent than the dominant power-law component in most bright ULXs. We show that the apparent X-ray luminosity and spectrum predicted by the microquasar model are consistent with the observed properties of most ULXs. We predict that the radio synchrotron jet emission is too faint to be detected at the typical threshold of radio surveys to date. This is consistent with the high rate of non-detections over detections in radio counterpart searches. Conversely, we conclude that the observed radio emission found associated with a few ULXs cannot be due to beamed synchrotron emission from a relativistic jet.