Hydrodynamical non-radiative accretion flows in two dimensions

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows that cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function ν . We study the flows generated as the amplitude and form of ν are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by and a large fraction of the mass participating in the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long-time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. The two-dimensional structure of the time-averaged flow is marginally stable to the Høiland criterion, indicating that convection is efficient. Near the equatorial plane, the radial profiles of the time-averaged variables are power laws with an index that depends on the radial scaling of the shear stress. A stress in which ν ∝ r ^1/2 recovers the widely studied self-similar solution corresponding to an ' α -disc'. We find that, regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.     

On the absence of winds in advection-dominated accretion flows

We show that recently published assertions that advection-dominated accretion flows (ADAFs) require the presence of strong winds are unfounded because they assume that low radiative efficiency in flows accreting at low rates on to black holes implies vanishing radial energy and angular momentum fluxes through the flow (which is also formulated in terms of the 'Bernoulli function' being positive). This, however, is a property only of self-similar solutions which are an inadequate representation of global accretion flows. We recall general properties of accretion flows on to black holes and show that such, necessarily transonic, flows may have either d positive or negative Bernoulli function depending on the flow viscosity. Flows with low viscosities ( α ≲0.1 in the α -viscosity model) have a negative Bernoulli function. Without exception, all 2D and 1D numerical models of low-viscosity flows constructed to date experience no significant outflows. At high viscosities the presence of outflows depends on the assumed viscosity, on the equation of state and on the outer boundary condition. The positive sign of the Bernoulli function invoked in this context is irrelevant to the presence of outflows. As an illustration, we recall 2D numerical models with moderate viscosity that have positive values of the Bernoulli function and experience no outflows. ADAFs, therefore, do not differ from this point of view from thin Keplerian discs: they may have, but they do not have to have, strong winds.     

Delayed X-ray emission from fallback in compact-object mergers

When double neutron star or neutron star–black hole binaries merge, the final remnant may comprise a central solar-mass black hole surrounded by a  ∼0.01–0.1 M_⊙  torus. The subsequent evolution of this disc may be responsible for short γ-ray bursts (SGRBs). A comparable amount of mass is ejected into eccentric orbits and will eventually fallback to the merger site after ∼0.01 s. In this paper, we investigate analytically the fate of the fallback matter, which may provide a luminous signal long after the disc is exhausted. We find that matter in the eccentric tail returns at a super-Eddington rate and eventually (≳0.1 s) is unable to cool via neutrino emission and accrete all the way to the black hole. Therefore, contrary to previous claims, our analysis suggests that fallback matter is not an efficient source of late-time accretion power and unlikely to cause the late-flaring activity observed in SGRB afterglows. The fallback matter rather forms a radiation-driven wind or a bound atmosphere. In both the cases, the emitting plasma is very opaque and photons are released with a degraded energy in the X-ray band. We therefore suggest that compact binary mergers could be followed by an 'X-ray renaissance', as late as several days to weeks after the merger. This might be observed by the next generation of X-ray detectors.     

Fractal features in accretion discs

Fractal concepts have been introduced in the accretion disc as a new feature. Due to the fractal nature of the flow, its continuity condition undergoes modifications. The conserved stationary fractal flow admits only saddle points and centre-type points in its phase portrait. Completely analytical solutions of the equilibrium point conditions indicate that the fractal properties enable the flow to behave like an effective continuum of lesser density, and facilitate the generation of transonicity. However, strongly fractal flows inhibit multitransonicity from developing. The mass accretion rate exhibits a fractal scaling behaviour, and the entire fractal accretion disc is stable under linearized dynamic perturbations.     

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.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

Unipolar outflows and global meridional circulations in rotating accretion flows

Using two-dimensional simulations of non-radiative viscous rotating black hole accretion flows, we show that the flows with α ∼0.1–0.3 self-organize to form stationary unipolar or bipolar outflows accompanied by global meridional circulations. The required energy comes, with efficiency ∼0.001–0.01, from the matter directly accreted on to the black hole. Observational implications are discussed.     

Parallel tracks in infrared versus X-ray emission in black hole X-ray transient outbursts: a hysteresis effect?

We report the discovery of a new hysteresis effect in black hole X-ray binary state transitions, that of the near-infrared (NIR) flux (which most likely originates in the jets) versus X-ray flux. We find, looking at existing data sets, that the IR emission of black hole X-ray transients appears to be weaker in the low/hard state rise of an outburst than the low/hard state decline of an outburst at a given X-ray luminosity. We discuss how this effect may be caused by a shift in the radiative efficiency of the inflowing or outflowing matter, or variations in the disc viscosity or the spectrum/power of the jet. In addition we show that there is a correlation (in slope but not in normalization) between IR and X-ray luminosities on the rise and decline, for all three low-mass black hole X-ray binaries with well-sampled IR and X-ray coverage:   L _NIR∝ L ^0.5–0.7_X  . In the high/soft state this slope is much shallower;   L _NIR∝ L ^0.1–0.2_X  , and we find that the NIR emission in this state is most likely dominated by the viscously heated (as opposed to X-ray heated) accretion disc in all three sources.     

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.