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⁻¹.     

A model of magnetically induced disc–corona for black hole binaries

We propose a model of magnetic connection (MC) of a black hole with its surrounding accretion disc based on large-scale magnetic field. The MC gives rise to transport of energy and angular momentum between the black hole and the disc, and the closed field lines pipe the hot matter evaporated from the disc, and shape it in the corona above the disc to form a magnetically induced disc–corona system, in which the corona has the same configuration as the large-scale magnetic field. We numerically solve the dynamic equations in the context of the Kerr metric, in which the large-scale magnetic field is determined by dynamo process and equipartition between magnetic pressure and gas pressure. Thus we can obtain a global solution rather than assuming the distribution of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) The rotational energy of a fast-spinning black hole can be extracted, enhancing the dissipation in the accretion disc, (2) the closed field lines provide a natural channel for corona matter escaping from disc and finally falling into black hole and (3) the scope of the corona can be bounded by the conservation of magnetic flux. We simulate the high-energy spectra of this system by using Monte Carlo method, and find that the relative hardness of the spectra decreases as accretion rate or black hole spin a _* increases. We fit the typical X-ray spectra of three black hole binaries  (GRO J1655−40, XTE 1118+480 and GX 339−4)  in the low/hard or very high state.     

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.     

Black hole masses from power density spectra: determinations and consequences

We analyse the scaling of the X-ray power density spectra with the mass of the black hole in the examples of Cyg X-1 and the Seyfert 1 galaxy NGC 5548. We show that the high-frequency tail of the power density spectrum can be successfully used for the determination of the black hole mass. We determine the masses of the black holes in six broad-line Seyfert 1 galaxies, five narrow-line Seyfert 1 galaxies and two quasi-stellar objects (QSOs) using the available power density spectra. The proposed scaling is clearly appropriate for other Seyfert galaxies and QSOs. In all but one of the normal Seyferts, the resulting luminosity to Eddington luminosity ratio is smaller than 0.15, with the source MCG -6-15-30 being an exception. The applicability of the same scaling to a narrow-line Seyfert 1 is less clear and there may be a systematic shift between the power spectra of NLS1 and S1 galaxies of the same mass, leading to underestimation of the black hole mass. However, both the method based on variability and the method based on spectral fitting show that those galaxies have relatively low masses and a high luminosity to Eddington luminosity ratio, supporting the view of those objects as analogues of galactic sources in their high, soft or very high state, based on the overall spectral shape. The bulge masses of their host galaxies are similar to that of normal Seyfert galaxies, so they do not follow the black hole mass–bulge mass relation for Seyfert galaxies, being evolutionarily less advanced, as suggested by Mathur. The bulge mass–black hole mass relation in our sample is consistent with being linear, with the black hole to bulge ratio ∼0.03 per cent, similar to Wandel and Laor for low-mass objects, but significantly shifted from the relation of Magorrian et al. and McLure & Dunlop.     

On the effect of coronal outflow on spectra formation in galactic black hole systems

We present the results of both analytical and numerical calculations of the amplitude of the reflection component in X-ray spectra of galactic black hole systems. We take into account the anisotropy of Compton scattering and the systematic relativistic bulk motion of the hot plasma. In the case of the single scattering approximation, the reflection from the disc surface is significantly enhanced owing to the anisotropy of Compton scattering. On the other hand, the calculations of multiple scattering obtained using the Monte Carlo method show that the anisotropy effect is much weaker in that case. Therefore, the enhanced back-scattered flux may affect the observed spectra only if the disc surface is highly ionized, which reduces the absorption in the energy band corresponding to the first Compton scattering.     

Corotational damping of discoseismic c modes in black hole accretion discs

Discoseismic c modes in accretion discs have been invoked to explain low-frequency variabilities observed in black hole X-ray binaries. These modes are trapped in the innermost region of the disc and have frequencies much lower than the rotation frequency at the disc inner radius. We show that because the trapped waves can tunnel through the evanescent barrier to the corotational wave zone, the c modes are damped due to wave absorption at the corotation resonance. We calculate the corotational damping rates of various c modes using the Wentzel-Kramers-Brillouin (WKB) approximation. The damping rate varies widely depending on the mode frequency, the black hole spin parameter and the disc sound speed, and is generally much less than 10 per cent of the mode frequency. A sufficiently strong excitation mechanism is needed to overcome this corotational damping and make the mode observable.     

The structure of black hole magnetospheres I. Schwarzschild black holes

We introduce a multipolar scheme for describing the structure of stationary, axisymmetric, force-free black hole magnetospheres in the '3+1' formalism. We focus here on Schwarzschild spacetime, giving a complete classification of the separable solutions of the stream equation. We show a transparent term-by-term analogy of our solutions with the familiar multipoles of flat-space electrodynamics. We discuss electrodynamic processes around disc-fed black holes in which our solutions find natural applications: (i) 'interior' solutions in studies of the BlandfordZnajek process of extracting the rotational energy of holes, and of the formation of relativistic jets in active galactic nuclei and 'microquasars'; (ii) 'exterior' solutions in studies of accretion disc dynamos, disc-driven winds and jets. On the strength of existing numerical studies, we argue that the poloidal field structures found here are also expected to hold with good accuracy for rotating black holes, except for the cases of the maximum possible rotation rates. We show that the closed-loop exterior solutions found here are not in contradiction with the MacdonaldThorne theorem, as these solutions, which diverge logarithmically on the horizon of the hole , only apply to those regions that exclude .     

Mass and spin co-evolution during the alignment of a black hole in a warped accretion disc

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with a misaligned accretion disc to study BH spin precession and alignment jointly with BH mass M _BH and spin parameter a evolution, under the assumption that the disc is continually fed, in its outer region, by matter with angular momentum fixed on a given direction     . We develop an iterative scheme based on the adiabatic approximation to study the BH–disc co-evolution: in this approach, the accretion disc transits through a sequence of quasi-steady warped states (Bardeen–Petterson effect) and interacts with the BH until the spin   J _BH  aligns with     . For a BH aligning with a corotating disc, the fractional increase in mass is typically less than a few per cent, while the spin modulus can increase up to a few tens of per cent. The alignment time-scale     is of  ∼10⁵–10⁶ yr  for a maximally rotating BH accreting at the Eddington rate. BH–disc alignment from an initially counter-rotating disc tends to be more efficient compared to the specular corotating case due to the asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and opposite angular momentum when crossing the innermost stable orbit, so that the spin modulus decreases faster and so the relative inclination angle.