Non-axisymmetric instabilities in shocked adiabatic accretion flows

We investigate the linear stability of a shocked accretion flow on to a black hole in the adiabatic limit. Our linear analyses and numerical calculations show that, despite the post-shock deceleration, the shock is generally unstable to non-axisymmetric perturbations. The simulation results of Molteni, Tóth & Kuznetsov can be well explained by our linear eigenmodes. The mechanism of this instability is confirmed to be based on the cycle of acoustic waves between the corotation radius and the shock. We obtain an analytical formula to calculate the oscillation period from the physical parameters of the flow. We argue that the quasi-periodic oscillation should be a common phenomenon in accretion flows with angular momentum.     

Magnetic and vertical shear instabilities in accretion discs

The stability properties of magnetized discs rotating with angular velocity Ω = Ω( s ,  z ), dependent on both the radial and the vertical coordinates s and z , are considered. Such a rotation law is adequate for many astrophysical discs (e.g., galactic and protoplanetary discs, as well as accretion discs in binaries). In general, the angular velocity depends on height, even in thin accretion discs. A linear stability analysis is performed in the Boussinesq approximation, and the dispersion relation is obtained for short-wavelength perturbations. Any dependence of Ω on z can destabilize the flow. This concerns primarily small-scale perturbations for which the stabilizing effect of buoyancy is strongly suppressed due to the energy exchange with the surrounding plasma. For a weak magnetic field, instability of discs is mainly associated with vertical shear, whilst for an intermediate magnetic field the magnetic shear instability, first considered by Chandrasekhar and Velikhov, is more efficient. This instability is caused by the radial shear which is typically much stronger than the vertical shear. Therefore the growth time for the magnetic shear instability is much shorter than for the vertical shear instability. A relatively strong magnetic field can suppress both these instabilities. The vertical shear instability could be the source of turbulence in protoplanetary discs, where the conductivity is low.     

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.     

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.     

Interface modes and their instabilities in accretion disc boundary layers

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the innermost regions of accretion discs or boundary layers. Two simple models are considered in this paper. The magnetosphere–disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star–disc model involves a Keplerian disc terminated at the stellar atmosphere with high density and small density scaleheight. We find that the interface modes at the magnetosphere–disc boundary are generally unstable due to Rayleigh–Taylor and/or Kelvin–Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh–Taylor instability, and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin–Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh–Taylor or Kelvin–Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.     

Effects of the angular momentum of accreting matter on the flow structure in the subsonic settling and Bondi-Hoyle accretion regimes

Previously unexplored accretion regimes associated with the rotation of accreting matter, namely the perturbations of a quasi-spherical subsonic settling flow and Bondi-Hoyle accretion in the presence of axial rotation, are considered within the framework of ideal hydrodynamics. For subsonic settling accretion, the perturbations are shown to grow rapidly as the gravitating center is approached, so that the flow in the inner regions can no longer be considered quasi-spherical. For Bondi-Hoyle accretion, a vacuum cylindrical cavity is shown to be formed at large distances from the gravitating center near the flow axis, with the flow velocity outside this cavity being virtually independent of the distance to the rotation axis.     

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.     

Physics of accretion in the millisecond pulsar XTE J1751 − 305

We have undertaken an extensive study of X-ray data from the accreting millisecond pulsar XTE J1751 − 305 observed by RXTE and XMM–Newton during its 2002 outburst. In all aspects this source is similar to the prototypical millisecond pulsar SAX J1808.4 − 3658, except for the higher peak luminosity of 13 per cent of Eddington, and the optical depth of the hard X-ray source, which is larger by a factor ∼2. Its broad-band X-ray spectrum can be modelled by three components. We interpret the two soft components as thermal emission from a colder  ( kT ∼ 0.6 keV)  accretion disc and a hotter (∼1 keV) spot on the neutron star surface. We interpret the hard component as thermal Comptonization in plasma of temperature ∼40 keV and optical depth ∼1.5 in a slab geometry. The plasma is heated by the accretion shock as the material collimated by the magnetic field impacts on to the surface. The seed photons for Comptonization are provided by the hotspot, not by the disc. The Compton reflection is weak and the disc is probably truncated into an optically thin flow above the magnetospheric radius. Rotation of the emission region with the star creates an almost sinusoidal pulse profile with an rms amplitude of 3.3 per cent. The energy-dependent soft phase lags can be modelled by two pulsating components shifted in phase, which is naturally explained by a different character of emission of the optically thick spot and optically thin shock combined with the action of the Doppler boosting. The observed variability amplitude constrains the hotspot to lie within 3°–4° of the rotational pole. We estimate the inner radius of the optically thick accreting disc to be about 40 km. In that case, the absence of emission from the antipodal spot, which can be blocked by the accretion disc, gives the inclination of the system as ≳70°.     

Possible quasi-periodic oscillations from unstable accretion: 3D magnetohydrodynamic simulations

We investigate the photometric variability of magnetized stars, particularly neutron stars, accreting through a magnetic Rayleigh–Taylor-type instability at the disc–magnetosphere interface, and compare it with the variability during stable accretion, with the goal of looking for possible quasi-periodic oscillations (QPOs). The light curves during stable accretion show periodicity at the star's frequency and sometimes twice that, due to the presence of two funnel streams that produce antipodal hotspots near the magnetic poles. On the other hand, light curves during unstable accretion through tongues penetrating the magnetosphere are more chaotic due to the stochastic behaviour of the tongues, and produce noisier power spectra. However, the power spectra do show some signs of quasi-periodic variability. Most importantly, the rotation frequency of the tongues and the resulting hotspots are close to the inner-disc orbital frequency, except in the most strongly unstable cases. There is therefore a high probability of observing QPOs at that frequency in longer simulations. In addition, the light curves in the unstable regime show periodicity at the star's rotation frequency in many of the cases investigated here, again except in the most strongly unstable cases which lack funnel flows and the resulting antipodal hotspots. The noisier power spectra result in the fractional rms amplitudes of the Fourier peaks being smaller.
We also study in detail the effect of the misalignment angle between the rotation and magnetic axes of the star on the variability, and find that at misalignment angles  ≳25°  the star's period always appears in the light curves.     

The radial–azimuthal instability of a hot two-temperature accretion disc with advection

The radial–azimuthal instability of a hot two-temperature accretion disc with advection is examined in this paper. We find that the inclusion of very little advection has significant effects on two acoustic modes for a geometrically thin, cooling-dominated two-temperature disc, but has no effect on acoustic modes for a geometrically slim, cooling-dominated two-temperature disc. We also find that, when azimuthal perturbations are considered, the stability properties of the disc are different from those in the pure radial perturbation case. An increase of the azimuthal wavenumber will stabilize the acoustic modes but make the viscous and thermal modes more unstable for a geometrically thin, cooling-dominated two-temperature disc. It makes the thermal mode more unstable and the acoustic mode more stable, but only affects the instability of the viscous mode for short-wavelength perturbations for a geometrically slim, cooling-dominated two-temperature disc. For a geometrically slim, advection-dominated two-temperature disc, the increase of the azimuthal perturbation makes the I- and O-modes more stable and the thermal mode more unstable, but has no effect on the viscous mode.     

Iron Kα line profiles and the inner boundary condition of accretion flows

Recent X-ray observations have shown evidence for exceptionally broad and skewed iron Kα emission lines from several accreting black hole systems. The lines are assumed to be due to fluorescence of the accretion disc illuminated by a surrounding corona and require a steep emissivity profile increasing into the innermost radius. This appears to question both standard accretion disc theory and the zero-torque assumption for the inner boundary condition, both of which predict a much less extreme profile. Instead it argues that a torque may be present due to magnetic coupling with matter in the plunging region or even to the spinning black hole itself. Discussion so far has centred on the torque acting on the disc. However, the crucial determinant of the iron line profile is the radial variation of the power radiated in the corona. Here we study the effects of different inner boundary conditions on the coronal emissivity and on the profiles of the observable Fe Kα lines. We argue that in the extreme case where a prominent highly redshifted component of the iron line is detected, requiring a steep emissivity profile in the innermost part and a flatter one outside, energy from the gas plunging into the black hole is being fed directly to the corona.     

Behaviour of dissipative accretion flows around black holes

We investigate the behaviour of dissipative accreting matter close to a black hole, as this provides important observational features of galactic and extragalactic black hole candidates. We find a complete set of global solutions in the presence of viscosity and synchrotron cooling. We show that advective accretion flow can have a standing shock wave and the dynamics of the shock is controlled by the dissipation parameters (both viscosity and cooling). We study the effective region of the parameter space for standing as well as oscillating shock. We find that the shock front always moves towards the black hole as the dissipation parameters are increased. However, viscosity and cooling have opposite effects in deciding the solution topologies. We obtain two critical cooling parameters that separate the nature of the accretion solution.     

Magnetic shear-flow instability in thin accretion discs

The possibility that the magnetic shear-flow instability (also known as the 'Balbus–Hawley' instability) might give rise to turbulence in a thin accretion disc is investigated through numerical simulations. The study is linear and the fluid disc is supposed to be incompressible and differentially rotating with a simple velocity profile with Ω∝ R ^{− q} . The simplicity of the model is counterbalanced by the fact that the study is fully global in all three spatial directions with boundaries on each side; finite diffusivities are also allowed. The investigation is also carried out for several values of the azimuthal wavenumber of the perturbations in order to analyse whether non-axisymmetric modes might be preferred, which may produce, in a non-linear extension of the study, a self-sustained magnetic field.
  We find the final pattern steady, with similar kinetic and magnetic energies and the angular momentum always transported outwards. Despite the differential rotation, there are only small differences for the eigenvalues for various non-axisymmetric eigensolutions. Axisymmetric instabilities are by no means preferred; in fact for Prandtl numbers between 0.1 and 1, the azimuthal wavenumbers m =0,1,2(10¹⁶ g s^-1). All three quantities appear to be equally readily excited. The equatorial symmetry is quadrupolar for the magnetic field and dipolar for the flow field system. The maximal magnetic field strength required to cause the instability is almost independent of the magnetic Prandtl number. With typical white dwarf values, a magnetic amplitude of 10⁵ G is estimated.     

Magnetospheric accretion and pre-main-sequence stellar masses

We present a method of determining lower limits on the masses of pre-main-sequence (PMS) stars and so constraining the PMS evolutionary tracks. This method uses the redshifted absorption feature observed in some emission-line profiles of T Tauri stars, indicative of infall. The maximum velocity of the accreting material measures the potential energy at the stellar surface, which, combined with an observational determination of the stellar radius, yields the stellar mass. This estimate is a lower limit owing to uncertainties in the geometry and projection effects. Using available data, we show that the computed lower limits can be larger than the masses derived from PMS evolutionary tracks for M   0.5 M. Our analysis also supports the notion that accretion streams do not impact near the stellar poles but probably hit the stellar surface at moderate latitudes.     

Spread of matter over a neutron-star surface during disk accretion: Deceleration of rapid rotation

The accretion of hot slowly rotating gas onto a supermassive black hole is considered. The important case where the velocities of turbulent pulsations at the Bondi radius r _B are low, compared to the speed of sound c _s, is studied. Turbulence is probably responsible for the appearance of random average rotation. Although the angular momentum at r _B is low, it gives rise to the centrifugal barrier at a depth r _c = l ² /GM _BH ≪ r _B, that hinders supersonic accretion. The numerical solution of the problem of hot gas accretion with finite angular momentum is found taking into account electron thermal conductivity and bremsstrahlung energy losses of two temperature plasma for density and temperature near Bondi radius similar to observed in M87 galaxy. The saturation of the Spitzer thermal conductivity was also taken into account. The parameters of the saturated electron thermal conductivity were chosen similar to the parameters used in the numerical simulations of interaction of the strong laser beam radiation with plasma targets. These parameters are confirmed in the experiments. It is shown that joint action of electron thermal conductivity and free-free radiation leads to the effective cooling of accreting plasma and formation of the subsonic settling of accreting gas above the zone of a centrifugal barrier. A toroidal condensation and a hollow funnel that separates the torus from the black hole emerge near the barrier. The barrier divides the flow into two regions: (1) the settling zone with slow subKeplerian rotation and (2) the zone with rapid supersonic nearly Keplerian rotation. Existence of the centrifugal barrier leads to significant decrease of the accretion rate Ṁ in comparison with the critical Bondi solution for γ = 5/3 for the same values of density and temperature of the hot gas near Bondi radius. Shear instabilities in the torus and related friction cause the gas to spread slowly along spirals in the equatorial plane in two directions.As a result, outer (r > r _c) and inner (r < r _c) disks are formed. The gas enters the immediate neighborhood of the black hole or the zone of the internal ADAF flow along the accretion disk (r < r _c). Since the angular momentum is conserved, the outer disk removes outward an excess of angular momentum along with part of the matter falling into the torus. It is possible, that such outer Keplerian disk was observed by Hubble Space Telescope around the nucleus of the M87 galaxy in the optical emission lines. We discuss shortly the characteristic times during which the accretion of the gas with developed turbulence should lead to the changes in the orientation of the torus, accretion disk and, possibly, of the jet.     

Analysing the atolls: X-ray spectral transitions of accreting neutron stars

We systematically analyse all the available X-ray spectra of disc accreting neutron stars (atolls and millisecond pulsars) from the RXTE data base. We show that while all these have similar spectral evolution as a function of mass accretion rate, there are also subtle differences. There are two different types of hard/soft transition, those where the spectrum softens at all energies, leading to a diagonal track on a colour–colour diagram, and those where only the higher energy spectrum softens, giving a vertical track. The luminosity at which the transition occurs is correlated with this spectral behaviour, with the vertical transition at   L / L _Edd∼ 0.02  while the diagonal one is at ∼0.1. Superimposed on this is the well-known hysteresis effect, but we show that classic, large-scale hysteresis occurs only in the outbursting sources, indicating that its origin is in the dramatic rate of change of mass accretion rate during the disc instability. We show that the long-term mass accretion rate correlates with the transition behaviour, and speculate that this is due to the magnetic field being able to emerge from the neutron star surface for low average mass accretion rates. While this is not strong enough to collimate the flow except in the millisecond pulsars, its presence may affect the inner accretion flow by changing the properties of the jet.     

Perturbative analysis of two-temperature radiative shocks with multiple cooling processes

The structure of the hot downstream region below a radiative accretion shock, such as that of an accreting compact object, may oscillate because of a global thermal instability. The oscillatory behaviour depends on the functional forms of the cooling processes, the energy exchanges of electrons and ions in the shock-heated matter, and the boundary conditions. We analyse the stability of a shock with unequal electron and ion temperatures, where the cooling consists of thermal bremsstrahlung radiation which promotes instability, plus a competing process which tends to stabilize the shock. The effect of transverse perturbations is considered also. As an illustration, we study the special case in which the stabilizing cooling process is of order 3/20 in density and 5/2 in temperature, which is an approximation for the effects of cyclotron cooling in magnetic cataclysmic variables. We vary the efficiency of the second cooling process, the strength of the electron–ion exchange and the ratio of electron and ion pressures at the shock, to examine particular effects on the stability properties and frequencies of oscillation modes.     

Photon bubbles in accretion discs

We show that radiation-dominated accretion discs are likely to suffer from a 'photon bubble' instability similar to that described by Arons in the context of accretion on to neutron star polar caps. The instability requires a magnetic field for its existence. In an asymptotic regime appropriate to accretion discs, we find that the overstable modes obey the remarkably simple dispersion relation
ο²=−i gkF ( B , k ).
Here g is the vertical gravitational acceleration, B is the magnetic field, and F is a geometric factor of order unity that depends on the relative orientation of the magnetic field and the wavevector. In the non-linear outcome it seems likely that the instability will enhance vertical energy transport and thereby change the structure of the innermost parts of relativistic accretion discs.     

The stability of accretion discs with inflow driven purely by magnetic winds

A model is presented for an accretion disc in which the inflow is driven purely by the angular momentum removed in a centrifugally accelerated magnetic wind. Turbulent discs around compact stars are considered, with the required magnetic field being generated in the disc by a simple dynamo. The turbulent magnetic Prandtl number, N _p, measures the ratio of turbulent viscosity to turbulent magnetic diffusivity. Formally, the hypothetical limit   N _p→ 0  corresponds to the magnetic wind torque dominating the viscous torque, but in practice the inflow is magnetically controlled for   N _p≲ 0.1  .
The suggestion by previous authors that purely magnetic wind-driven discs may be unstable is investigated. A detailed steady solution is found which allows perturbations to the thermal balance and vertical equilibrium to be calculated, and hence the effect of perturbations to the magnetic diffusivity, η, to be assessed. For a standard parametrized form of η, the wind-driven angular momentum balance is found to be linearly unstable. An increase in the inflow rate leads to increased bending of the poloidal magnetic field and an enhanced wind mass loss rate. This increases the angular momentum loss rate which drives further inflow. There is a resultant increase in η, due to the temperature perturbation, but this does not relieve field bending sufficiently to prevent the instability.