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.     

A unified picture for the γ-ray and prompt optical emissions of GRB 990123

The prompt optical emission of GRB 990123 was uncorrelated to the γ-ray light curve and exhibited temporal properties similar to those of the steeply decaying, early X-ray emission observed by Swift at the end of many bursts. These facts suggest that the optical counterpart of GRB 990123 was the large-angle emission released during (the second pulse of) the burst. If the optical and γ-ray emissions of GRB 990123 have, indeed, the same origin then their properties require that (i) the optical counterpart was synchrotron emission and γ-rays arose from inverse-Compton scatterings (the 'synchrotron self-Compton model'), (ii) the peak energy of the optical-synchrotron component was at ∼20 eV and (iii) the burst emission was produced by a relativistic outflow moving at Lorentz factor  ≳450  and at a radius  ≳10¹⁵  cm, which is comparable to the outflow deceleration radius. Because the spectrum of GRB 990123 was optically thin above 2 keV, the magnetic field behind the shock must have decayed on a length-scale of  ≲1  per cent  of the thickness of the shocked gas, which corresponds to  10⁶–10⁷  plasma skin depths. Consistency of the optical counterpart decay rate and its spectral slope (or that of the burst, if they represent different spectral components) with the expectations for the large-angle burst emission represents the most direct test of the unifying picture proposed here for GRB 990123.     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the gamma-ray burst (GRB) prompt emission site from the explosion centre R , by determining the location of the electron's self-absorption frequency in the GRB prompt optical-to-X/γ-ray spectral energy distribution, assuming that the optical and the γ-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely the bulk Lorentz factor of the emitting source Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 300  . For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R -constraints for this non-detection sample and those for the detection sample.     

Radiation mechanisms and geometry of Cygnus X-1 in the soft state

We present X-ray/ γ -ray spectra of Cyg X-1 observed during the transition from the hard to the soft state and in the soft state by ASCA , RXTE and CGRO /OSSE in 1996 May and June. The spectra consist of a dominant soft component below ∼2 keV and a power-law-like continuum extending to at least ∼800 keV. We interpret them as emission from an optically thick, cold accretion disc and from an optically thin, non-thermal corona above the disc. A fraction f ≳0.5 of total available power is dissipated in the corona.
We model the soft component by multicolour blackbody disc emission taking into account the torque-free inner-boundary condition. If the disc extends down to the minimum stable orbit, the ASCA RXTE data yield the most probable black hole mass of M _X≈10 M_⊙ and an accretion rate,     , locating Cyg X-1 in the soft state in the upper part of the stable, gas-pressure-dominated, accretion-disc solution branch.
The spectrum of the corona is well modelled by repeated Compton scattering of seed photons from the disc off electrons with a hybrid, thermal/non-thermal distribution. The electron distribution can be characterized by a Maxwellian with an equilibrium temperature of kT _e∼30–50 keV, a Thomson optical depth of τ ∼0.3 and a quasi-power-law tail. The compactness of the corona is 2≲ℓ_h≲7, and a presence of a significant population of electron–positron pairs is ruled out.
We find strong signatures of Compton reflection from a cold and ionized medium, presumably an accretion disc, with an apparent reflector solid angle, Ω/2π∼0.5–0.7. The reflected continuum is accompanied by a broad iron K α line.     

Radiatively driven mass accretion on to galactic nuclei by circumnuclear starbursts

We examine the physical processes of radiatively driven mass accretion on to galactic nuclei, owing to intensive radiation from circumnuclear starbursts. The radiation from a starburst not only causes the inner gas disc to contract via radition flux force, but also extracts angular momentum owing to relativistic radiation drag, thereby inducing an avalanche of the surface layer of the disc. To analyse such a mechanism, the radiation–hydrodynamical equations are solved, including the effects of the radiation drag force as well as the radiation flux force. As a result, it is found that the mass accretion rate owing to the radiative avalanche is given by M ˙ ( r )= η ( L _*/ c ²)( r / R )² (Δ R / R )(1 −  e ^−τ) at radius r , where the efficiency η ranges from 0.2 up to 1, L _* and R are respectively the bolometric luminosity and the radius of the starburst ring, Δ R is the extent of the emission regions, and τ is the face-on optical depth of the disc. In an optically thick regime, the rate depends upon neither the optical depth nor the surface mass density distribution of the disc. The present radiatively driven mass accretion may provide a physical mechanism which enables mass accretion from 100-pc scales down to ∼ parsec scales, and it may eventually be linked to advection-dominated viscous accretion on to a massive black hole. The radiation–hydrodynamical and self-gravitational instabilities of the disc are briefly discussed. In particular, the radiative acceleration possibly builds up a dusty wall, which 'shades' the nucleus in edge-on views. This provides another version of the model for the formation of an obscuring torus.     

Jet and accretion power in the most powerful Fermi blazars

Among the blazars detected by the Fermi satellite, we have selected the 23 blazars that in the 3 months of survey had an average γ-ray luminosity above 10⁴⁸ erg s⁻¹. For 17 out of the 23 sources we found and analysed X-ray and optical–ultraviolet data taken by the Swift satellite. With these data, implemented by archival and not simultaneous data, we construct the spectral energy distributions, and interpreted them with a simple one-zone, leptonic, synchrotron and inverse Compton model. When possible, we also compare different high-energy states of single sources, like 0528+134 and 3C 454.3, for which multiple good sets of multiwavelength data are available. In our powerful blazars the high energy emission always dominates the electromagnetic output, and the relatively low level of the synchrotron radiation often does not hide the accretion disc emission. We can then constrain the black hole mass and the disc luminosity. Both are large (i.e. masses equal or greater than  10⁹ M _⊙  and disc luminosities above 10 per cent of Eddington). By modelling the non-thermal continuum we derive the power that the jet carries in the form of bulk motion of particles and fields. On average, the jet power is found to be slightly larger than the disc luminosity, and proportional to the mass accretion rate.     

Ultraluminous X-ray source 1E 0953.8+6918 (M81 X-9): an intermediate-mass black hole candidate and its environs

We present a ROSAT and ASCA study of the Einstein source X-9 and its relation to a shock-heated shell-like optical nebula in a tidal arm of the M81 group of interacting galaxies. Our ASCA observation of the source shows a flat and featureless X-ray spectrum well described by a multicolour disc blackbody model. The source most likely represents an optically thick accretion disc around an intermediate-mass black hole  ( M ∼10² M_⊙)  in its high/soft state, similar to other variable ultraluminous X-ray sources observed in nearby disc galaxies. Using constraints derived from both the innermost stable orbit around a black hole and the Eddington luminosity, we find that the black hole is fast-rotating and that its mass is between ∼80 M_⊙–1.5×10² M_⊙. The inferred bolometric luminosity of the accretion disc is ∼(1.1×10⁴⁰ erg s⁻¹)/(cos  i ). Furthermore, we find that the optical nebula is very energetic and may contain large amounts of hot gas, accounting for a soft X-ray component as indicated by archival ROSAT PSPC data. The nebula is apparently associated with X-9; the latter may be powering the former and/or they could be formed in the same event (e.g. a hypernova). Such a connection, if confirmed, could have strong implications for understanding both the birth of intermediate-mass black holes and the formation of energetic interstellar structures.     

A unified treatment of the gamma-ray burst 021211 and its afterglow

The gamma-ray burst (GRB) 021211 had a simple light curve, containing only one peak and the expected Poisson fluctuations. Such a burst may be attributed to an external shock, offering the best chance for a unified understanding of the gamma-ray burst and afterglow emissions. We analyse the properties of the prompt (burst) and delayed (afterglow) emissions of GRB 021211 within the fireball model. Consistency between the optical emission during the first 11 min (which, presumably, comes from the reverse shock heating of the ejecta) and the later afterglow emission (arising from the forward shock) requires that, at the onset of deceleration (∼2 s), the energy density in the magnetic field in the ejecta, expressed as a fraction of the equipartition value  (ɛ _B )  , is larger than in the forward shock at 11 min by a factor of approximately 10³. We find that synchrotron radiation from the forward shock can account for the gamma-ray emission of GRB 021211; to explain the observed GRB peak flux requires that, at 2 s,  ɛ _B   in the forward shock is larger by a factor 100 than at 11 min. These results suggest that the magnetic field in the reverse shock and early forward shock is a frozen-in field originating in the explosion and that most of the energy in the explosion was initially stored in the magnetic field. We can rule out the possibility that the ejecta from the burst for GRB 021211 contained more than 10 electron–positron pairs per proton.     

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.     

Evidence for energy injection and a fine-tuned central engine at optical wavelengths in GRB 070419A

We present a comprehensive multiwavelength temporal and spectral analysis of the 'fast rise exponential decay' GRB 070419A. The early-time emission in the γ-ray and X-ray bands can be explained by a central engine active for at least 250 s, while at late times the X-ray light curve displays a simple power-law decay. In contrast, the observed behaviour in the optical band is complex (from 10² up to 10⁶ s). We investigate the light-curve behaviour in the context of the standard forward/reverse shock model; associating the peak in the optical light curve at ∼450 s with the fireball deceleration time results in a Lorenz factor  Γ≈ 350  at this time. In contrast, the shallow optical decay between 450 and 1500 s remains problematic, requiring a reverse shock component whose typical frequency is above the optical band at the optical peak time for it to be explained within the standard model. This predicts an increasing flux density for the forward shock component until   t ∼ 4 × 10⁶ s  , inconsistent with the observed decay of the optical emission from   t ∼ 10⁴ s  . A highly magnetized fireball is also ruled out due to unrealistic microphysic parameters and predicted light-curve behaviour that is not observed. We conclude that a long-lived central engine with a finely tuned energy injection rate and a sudden cessation of the injection is required to create the observed light curves, consistent with the same conditions that are invoked to explain the plateau phase of canonical X-ray light curves of γ-ray bursts.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

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.     

Small-scale inviscid accretion discs around black holes

Gas falling quasi-spherically on to a black hole forms an inner accretion disc if its specific angular momentum l exceeds l ∗∼ r _g c , where r _g is the Schwarzschild radius. The standard disc model assumes l ≫ l ∗. We argue that, in many black hole sources, accretion flows have angular momenta just above the threshold for disc formation, l ≳ l ∗, and assess the accretion mechanism in this regime. In a range l ∗< l < l _cr, a small-scale disc forms in which gas spirals fast into the black hole without any help from horizontal viscous stresses. Such an 'inviscid' disc, however, interacts inelastically with the feeding infall. The disc–infall interaction determines the dynamics and luminosity of the accretion flow. The inviscid disc radius can be as large as 14 r _g, and the energy release peaks at 2 r _g. The disc emits a Comptonized X-ray spectrum with a break at ∼100 keV. This accretion regime is likely to take place in wind-fed X-ray binaries and is also possible in active galactic nuclei.     

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.     

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.     

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.     

Outbursts of irradiated accretion discs

I solve analytically the viscous evolution of an irradiated accretion disc, as seen during outbursts of soft X-ray transients. The solutions predict steep power-law X-ray decays L _X ∼ (1 + t/t_visc)⁻⁴, changing to L _X ∼ (1 − t/t'_visc)⁴ at late times, where t _visc, t '_visc are viscous time-scales. These forms closely resemble the approximate exponential and linear decays inferred by King and Ritter in these two regimes. The decays are much steeper than for unirradiated discs because the viscosity is a function of the central accretion rate rather than of local conditions in the disc.     

Effects of non-thermal tails in Maxwellian electron distributions on synchrotron and Compton processes

We investigate how the presence of a non-thermal tail beyond a Maxwellian electron distribution affects the synchrotron process as well as Comptonization in plasmas with parameters typical for accretion flows on to black holes. We find that the presence of the tail can significantly increase the net (after accounting for self-absorption) cyclo-synchrotron emission of the plasma, which then provides seed photons for Compton upscattering. Thus, the luminosity in the thermally Comptonized spectrum is enhanced as well. The importance of these effects increases with both increasing Eddington ratio and black hole mass. The enhancement of the Comptonized synchrotron luminosity can be as large as ∼10³ and ∼10⁵ for stellar and supermassive black holes, respectively, when the energy content in the non-thermal tail is 1 per cent.
The presence of the tail only weakly hardens the thermal Comptonization spectrum but it leads to the formation of a high-energy tail beyond the thermal cut-off, which two effects are independent of the nature of the seed photons. Since observations of high-energy tails in Comptonization spectra can constrain the non-thermal tails in the electron distribution and thus the Comptonized synchrotron luminosity, they provide upper limits on the strength of magnetic fields in accretion flows. In particular, the measurement of an MeV tail in the hard state of Cyg X-1 by McConnell et al. implies the magnetic field strength in this source to be at most an order of magnitude below equipartition.     

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.     

Spectral variability of ultraluminous X-ray sources

We study spectral variability of 11 ultraluminous X-ray sources (ULX) using archived XMM–Newton and Chandra observations. We use three models to describe the observed spectra: a power law, a multicolour disc (MCD) and a combination of these two models. We find that seven ULXs show a correlation between the luminosity L _X and the photon index Γ. Furthermore, four out of these seven ULXs also show spectral pivoting in the observed energy band. We also find that two ULXs show an   L _X–Γ  anticorrelation. The spectra of four ULXs in the sample can be adequately fitted with a MCD model. We compare these sources to known black hole binaries (BHB) and find that they follow similar paths in their luminosity–temperature diagrams. Finally, we show that the 'soft excess' reported for many of these ULXs at ∼0.2 keV seems to roughly follow a trend   L _soft∝ T ^−3.5  when modelled with a power law plus a 'cool' MCD model. This is contrary to the   L ∝ T ⁴  relation that is expected from theory and what is seen for many accreting BHBs. The observed trend could instead arise from disc emission beamed by an outflowing wind around a  ∼10 M_⊙  black hole.