Variability associated with α accretion disc theory for standard and advection-dominated discs

The α turbulent viscosity formalism for accretion discs must be interpreted as a mean field theory, modelling a steady state only on spatial or time-scales greater than those of the turbulence. The extent of the scale separation determines the relative precision error (RPE) of the predicted luminosity L _ν. Turbulence and the use of α implies that (1) field line stretching gives a magnetic pressure  α²/6 of the total pressure generally, and a one-to-one relation between α and the pressure ratio for thin discs, and (2) large turbulent scales in advection-dominated accretion flows (ADAFs) predict a lower L _ν precision than thin discs for a given observation duration and central mass. The allowed variability (or RPE) at frequency ν increases with the size of the contributing region. For X-ray binary ADAFs, the RPE ∼ 5 per cent at R  ≤ 1000 Schwarzchild radii ( R _s) for averages over  1000 s. However, current data for galaxies like NGC 4258 and M87 give RPEs in L _ν of 50–100 per cent even at R  ≤ 100  R _S. More data are required, but systematic deviations from ADAF predictions are more significant than random deviations, and may constrain properties of the turbulence, the accretion mode, the assumption of a steady state or the accretion rate.     

The circumnuclear material in the Galactic Centre: a clue to the accretion process

On the basis of 'sticky particle' calculations, it is argued that the gas features observed within 10 pc of the Galactic Centre — the circumnuclear disc (CND) and the ionized gas filaments — as well as the newly formed stars in the inner 1 pc can be understood in terms of tidal capture and disruption of gas clouds on low angular momentum orbits in a potential containing a point mass. The calculations demonstrate that a dissipative component forms a 'dispersion ring', an asymmetric elliptical torus precessing counter to the direction of rotation, and that this shape can be maintained for many orbital periods. For a range of plausible initial conditions, such a structure can explain the morphology and kinematics of the CND and of the most conspicuous ionized filament. While forming the dispersion ring, a small cloud with low specific angular momentum is drawn into a long filament which repeatedly collides with itself at high velocity. The compression in strong shocks is likely to lead to star formation even in the near tidal field of the point mass. This process may have general relevance to accretion on to massive black holes in normal and active galactic nuclei.     

Variable accretion and emission from the stellar winds in the Galactic Centre

We present numerical simulations of stellar wind dynamics in the central parsec of the Galactic Centre, studying in particular the accretion of gas on to Sgr A*, the supermassive black hole. Unlike our previous work, here we use state-of-the-art observational data on orbits and wind properties of individual wind-producing stars. Since wind velocities were revised upwards and non-zero eccentricities were considered, our new simulations show fewer clumps of cold gas and no conspicuous disc-like structure. The accretion rate is dominated by a few close 'slow-wind stars' ( v _w≤ 750 km s⁻¹), and is consistent with the Bondi estimate, but variable on time-scales of tens to hundreds of years. This variability is due to the stochastic infall of cold clumps of gas, as in earlier simulations, and to the eccentric orbits of stars. The present models fail to explain the high luminosity of Sgr A* a few hundred years ago implied by Integral observations, but we argue that the accretion of a cold clump with a small impact parameter could have caused it. Finally, we show the possibility of constraining the total mass-loss rate of the 'slow-wind stars' using near infrared observations of gas in the central few arcseconds.     

Magnetized accretion flows: effects of gas pressure

We study how axisymmetric magnetohydrodynamic (MHD) accretion flows depend on γ adiabatic index in the polytropic equation of state. This work is an extension of Mościbrodzka & Proga, where we investigated the γ dependence of two-dimensional Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various γ can give us insights into the problem of various modes of accretion observed in several types of accretion systems, such as black hole binaries (BHBs), active galactic nuclei (AGN) and gamma-ray bursts. We find that for  γ≳ 4/3  , the fast-rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for  γ= 5/3  , such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low γ, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHBs during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.     

Magnetized, mass-loaded, rotating accretion flows

We present a semi-analytical investigation of a simple one-dimensional, steady-state model for a mass-loaded, rotating, magnetized, hydrodynamical flow. Our approach is analogous to one used in early studies of magnetized winds. The model represents the infall towards a central point mass of the gas generated in a cluster of stars surrounding it, as is likely to occur in some active nuclei and starburst galaxies. We describe the properties of the different classes of infall solutions. We find that the flow becomes faster than the fast-mode speed, and hence decoupled from the centre, only for a limited range of parameter values, and when magnetic stresses are ineffective. Such flow is slowed as it approaches a centrifugal barrier, implying the existence of an accretion disc. When the flow does not become super-fast and the magnetic torque is insufficient, no steady solution extending inward to the centre exists. Finally, with a larger magnetic torque, solutions representing steady sub-Alfvénic flows are found, which can resemble spherical hydrodynamical infall. Such solutions, if applicable, would imply that rotation is not important and that any accretion disc formed would be of very limited size.     

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.     

Constraining the number of compact remnants near Sgr A*

Due to dynamical friction stellar mass black holes and neutron stars are expected to form high-density cusps in the inner parsec of our Galaxy. These compact remnants, expected to number around 20 000, may be accreting cold dense gas present there, and give rise to potentially observable X-ray emission. Here we build a simple but detailed time-dependent model of such emission. The possibility that these accretion flows are radiatively inefficient is taken into account and brings in some uncertainty in the conclusions. Despite this uncertainty, we find that at least several X-ray sources of this nature should be detectable with Chandra at any one time. Turning this issue around, we also ask a question of what current observational constraints might be telling us about the total number of compact remnants. In our 'best guess' model, a cusp of ∼40 000 remnants overpredicts the number of discrete sources and the total X-ray luminosity of the inner parsec, and is hence ruled out. In the most radiatively inefficient scenario that we consider, the radiative efficiency is set to be as small as  ɛ= 10⁻⁵  . In this rather unlikely scenario, a cusp of ∼40 000 black holes would be allowed by the data, but several individual sources should still be visible. Future observations of the distribution and orbits of the cold ionized gas in the inner parsec of our Galaxy will put tighter constraints on the cusp of compact remnants.     

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.     

Stochastic wobble of accretion discs and jets from turbulent rocket torques

Models of accretion discs and their associated outflows often incorporate assumptions of axisymmetry and symmetry across the disc plane. However, for turbulent discs these symmetries only apply to averaged quantities and do not apply locally. The local asymmetries can induce local imbalances in outflow power across the disc mid-plane, which can in turn induce local tilting torques. Here we calculate the effect of the resulting stochastic torques on disc annuli that are a consequence of standard mean field accretion disc models. The torques induce a random walk of the vector perpendicular to the plane of each averaged annulus. This random walk is characterized by a radially dependent diffusion coefficient which we calculate for small angle tilt. We use the coefficient to calculate a radially dependent time-scale for annular tilt and associated jet wobble. The wobble time depends on the square of the wander angle so the age of a given system determines the maximum wobble angle. We apply this to examples of blazars, young stellar objects and binary engines of pre-planetary nebulae and microquasars. It is noteworthy that for an averaging time   t _w∼ 3 d  , we estimate a wobble angle for jets in SS 433 of  θ∼ 0.8°  , not inconsistent with observational data. In general the non-periodic nature of the stochastic wobble could distinguish it from faster periodic jet precession.     

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.     

On the precession of accretion discs in X-ray binaries

In this Letter, recent results on the nodal precession of accretion discs in close binaries are applied to the discs in some X-ray binary systems. The ratio between the tidally forced precession period and the binary orbital period is given, as well as the condition required for the rigid precession of gaseous Keplerian discs. Hence the minimum precessional period that may be supported by a fluid Keplerian disc is determined. It is concluded that near-rigid body precession of tilted accretion discs can occur and generally reproduce observationally inferred precession periods, for reasonable system parameters. In particular, long periods in SS 433, Her X-1, LMC X-4 and SMC X-1 can be fitted by the tidal model. It is also found that the precession period that has been tentatively put forward for Cyg X-2 cannot be accommodated by a tidally precessing disc model for any realistic choice of system parameters.