The dependence on system parameters of strange-mode instabilities in accretion discs

A systematic study of the dependence on disc parameters and input physics, such as opacity and the treatment of convection, of strange-mode instabilities in thin accretion discs, which have been discovered recently, is presented. The instabilities are found to exist for a wide range of parameters, are partly very robust, and their growth rates can reach the dynamical range. Even discs on galactic scales around massive black holes are affected by them. Two groups of instabilities can be distinguished, the first of which is related to the radiation-pressure-dominated part of the disc, and the second to helium/hydrogen ionization. By application of the NAR approximation, both of them can be shown to be of mechanical origin, and the classical κ -mechanism can be excluded as the instability mechanism. A heuristic model for strange-mode instabilities proposed in the context of stellar strange-mode instabilities in luminous stars seems to be applicable to the group associated with dominant radiation pressure.     

Magnetic field bending in accretion discs with dynamos

We consider the problem of poloidal magnetic field advection and bending of an initially vertical field owing to radial inflow in thin accretion discs. For a ratio of kinematic viscosity to magnetic diffusivity of order unity, significant bending of an externally applied vertical field cannot occur in a disc with no internal dynamo. However, we show that if poloidal field is generated by a dynamo operating near its critical state, then significant field bending may be possible. Our results are of particular relevance to wind launching from accretion discs.     

On the maximal value of the turbulent α ‐parameter in accretion discs

In this short paper we show that making turbulence two‐rather than three‐dimensional may increase the effective turbulent viscosity by about 40 %. Dimensionless hydrodynamical viscosity parameters up to α_max = 0.25 M_t² may be obtained in this approach, which are in better agreement with the observational data on non‐stationary accretion than the values obtained in numerical simulations. However, the α ‐parameter values known from observations are still several times higher (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Accretion disc viscosity: how big is alpha?

We consider observational and theoretical estimates of the accretion disc viscosity parameter α. We find that in thin, fully ionized discs, the best observational evidence suggests a typical range α∼ 0.1–0.4, whereas the relevant numerical simulations tend to derive estimates for α which are an order of magnitude smaller. We discuss possible reasons for this apparent discrepancy.     

Disk-Corona Model in Active Galactic Nuclei： an Observational Test

We compiled a sample of 98 radio-quiet active galactic nuclei observed by ASCA, Chandra, XMM-Newton, INTEGRAL and Swift with the aim of testing the formation of hot corona and the magnetic shear stress operating in a disk-corona system. We found a strong correlation between the hard X-ray luminosity, bolometric luminosity L Bol and Eddington luminosity LEdd, in the sense that the fraction f of hard X-ray to the bolometric luminosity is inversely proportional to the Eddington ratio. This correlation favors the shear stress tensor being of the form of trφ∝ Pgas, with which the disk-corona structure is stable.     

The photoevaporation of discs around young stars in massive clusters

We present models in which the photoevaporation of discs around young stars by an external ultraviolet source (as computed by Adams et al.) is coupled with the internal viscous evolution of the discs. These models are applied to the case of the Orion Nebula Cluster (ONC), where the presence of a strong ultraviolet field from the central OB stars, together with a detailed census of circumstellar discs and photoevaporative flows, is well established. In particular we investigate the constraints that are placed on the initial disc properties in the ONC by the twin requirement that most stars possess a disc on a scale of a few astronomical unit (au), but that only a minority (<20 per cent) are resolved by Hubble Space Telescope ( HST ) at a scale of 50 au. We find that these requirements place very weak constraints on the initial radius distribution of circumstellar discs: the resulting size distribution readily forgets the initial radius distribution, owing to the strong positive dependence of the photoevaporation rate on disc radius. Instead, the scarcity of large discs reflects the relative scarcity of initially massive discs (with mass  >0.1 M_⊙  ). The ubiquity of discs on a small scale, on the other hand, mainly constrains the time-span over which the discs have been exposed to the ultraviolet field (<2 Myr). We argue that the discs that are resolved by HST represent a population of discs in which self-gravity was important at the time that the dominant central OB star switched on, but that, according to our models, self-gravity is unlikely to be important in these discs at the present time. We discuss the implications of our results for the so-called proplyd lifetime problem.     

Non-linear bending waves in Keplerian accretion discs

The non-linear dynamics of a warped accretion disc is investigated in the important case of a thin Keplerian disc with negligible viscosity and self-gravity. A one-dimensional evolutionary equation is formally derived that describes the primary non-linear and dispersive effects on propagating bending waves other than parametric instabilities. It has the form of a derivative non-linear Schrödinger (DNLS) equation with coefficients that are obtained explicitly for a particular model of a disc. The properties of this equation are analysed in some detail and illustrative numerical solutions are presented. The non-linear and dispersive effects both depend on the compressibility of the gas through its adiabatic index Γ. In the physically realistic case Γ < 3, non-linearity does not lead to the steepening of bending waves but instead enhances their linear dispersion. In the opposite case Γ > 3, non-linearity leads to wave steepening and solitary waves are supported. The effects of a small effective viscosity, which may suppress parametric instabilities, are also considered. This analysis may provide a useful point of comparison between theory and numerical simulations of warped accretion discs.     

Evolution of initially localized perturbations in stratified ionized discs

A detailed solution of an initial value problem of a vertically localized initial perturbation in rotating magnetized vertically stratified disc is presented. The appropriate linearized magnetohydrodynamics equations are solved by employing the Wentzel–Kramers–Brillouin (WKB) approximation and the results are verified numerically. The eigenfrequencies as well as eigenfunctions are explicitly obtained. It is demonstrated that the initial perturbation remains confined within the disc. It is further shown that thin enough discs are stable but as their thickness grows increasing number of unstable modes participate in the solution of the initial value problem. However, it is demonstrated that due to the localization of the initial perturbation, the growth time of the instability is significantly longer than the calculated inverse growth rate of the individual unstable eigenfunctions.     

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.     

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.     

X-ray irradiation in XTE J1817−330 and the inner radius of the truncated disc in the hard state

The key aspect of the very successful truncated disc model for the low/hard X-ray spectral state in black hole binaries is that the geometrically thin disc recedes from the last stable orbit at the transition to this state. This has recently been challenged by direct observations of the low/hard state disc from CCD data. We reanalyse the Swift and RXTE campaign covering the 2006 outburst of XTE J1817−330, and show that these data actually strongly support the truncated disc model as the transition spectra unambiguously show that the disc begins to recede as the source leaves the disc-dominated soft state. The disc radius inferred for the proper low/hard state is less clear-cut, but we show that the effect of irradiation from the energetically dominant hot plasma leads to an underestimate of the disc radius by a factor of 2–3 in this state. This may also produce the soft excess reported in some hard-state spectra. The inferred radius becomes still larger when the potential difference in stress at the inner boundary, increased colour temperature correction from incomplete thermalization of the irradiation, and loss of observable disc photons from Comptonization in the hot plasma is taken into account. We conclude that the inner disc radius in XTE J1817−330 in the low/hard spectral state is at least six to eight times that seen in the disc-dominated high/soft state, and that recession of the inner disc is the trigger for the soft-hard-state transition, as predicted by the truncated disc models.     

Temporal behaviour of global perturbations in compressible axisymmetric flows with free boundaries

The dynamics of small global perturbations in the form of a linear combination of a finite number of non‐axisymmetric eigenmodes is studied in the two‐dimensional approximation. The background flow is assumed to be an axisymmetric perfect fluid with adiabatic index γ = 5/3 rotating with a power law angular velocity distribution Γ ∝ r^–q , 1.5 < q < 2.0, confined by free boundaries in the radial direction. The substantial transient growth of acoustic energy of optimized perturbations is discovered. An optimal energy growth G is calculated numerically for a variety of parameters. Its value depends essentially on the perturbation azimuthal wavenumber m and increases for higher values of m. The closer the rotation profile to the Keplerian law, the larger growth factors can be obtained but over a longer time. The highest acoustic energy increase found numerically is of order ∼10² over ∼6 typical Keplerian periods. Slow neutral eigenmodes with corotation radius beyond the outer boundary mostly contribute to the transient growth. The revealed linear temporal behaviour of perturbations may play an important role in angular momentum transfer in toroidal flows near compact relativistic objects (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.