Tilted accretion discs in cataclysmic variables: tidal instabilities and superhumps

We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three-dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhumps might arise from a tilted disc, but we find no evidence that the companion gravitational tilt instability can grow fast enough in a fluid disc to create a measurable inclination. The origin of negative superhumps in the light curves of cataclysmic variables remains a puzzle.     

Effects of local thermodynamics and of stellar mass ratio on accretion disc stability in close binaries

Inflow kinematics at the inner Lagrangian point L1, gas compressibility, and physical turbulent viscosity play a fundamental role on accretion disc dynamics and structure in a close binary (CB). Physical viscosity supports the accretion disc development inside the primary gravitational potential well, developing the gas radial transport, converting mechanical energy into heat. The Stellar‐Mass‐Ratio (SMR) between the compact primary and the secondary star (M₁/M₂) is also effective, not only in the location of the inner Lagrangian point, but also in the angular kinematics of the mass transfer and in the geometry ofthe gravitational potential wells. In this work we pay attention in particular to the role ofthe SMR, evaluating boundaries, separating theoretical domains in compressibility‐viscosity graphs where physical conditions allow a well‐bound disc development, as a function ofmass transfer kinematic conditions. In such domains, the lower is the gas compressibility (the higher the polytropic index γ), the higher is the physical viscosity (α) requested. In this work, we show how the boundaries of such domains vary as a function of M₁/M₂. Conclusions as far as dwarf novae outbursts are concerned, induced by mass transfer rate variations, are also reported. The smaller M₁/M₂, the shorter the duration of the active‐to‐quiet and vice‐versa transitional phases. Time‐scales are of the order of outburst duration of SU Uma, OY Car, Z Cha and SS Cyg‐like objects. Moreover, conclusions as far as active‐quiet‐active phenomena in a CB, according to viscous‐thermal instabilities, in accordance to such domains, are also reported (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On accretion flow penetration of magnetospheres

We address the problem of plasma penetration of astrophysical magnetospheres, an important issue in a wide variety of contexts, ranging from accretion in cataclysmic variables to flows in protostellar systems. We point out that under well-defined conditions, penetration can occur without any turbulent mixing (driven, for example, by Rayleigh–Taylor or Kelvin–Helmholtz instabilities) caused by charge polarization effects, if the inflowing plasma is bounded in the direction transverse to both the flow velocity and the magnetic field. Depolarization effects limit the penetration depth, which nevertheless can, under specific circumstances, be comparable to the size of the magnetosphere. We discuss the effect of ambient medium on plasma propagation across the stellar magnetic field and determine the criteria for deep magnetosphere penetration. We show that, under conditions appropriate to magnetized white dwarfs in AM Her type cataclysmic variables, charge polarization effects can lead to deep penetration of the magnetosphere.     

Hydrodynamical activity in thin accretion disks

An asymptotic treatment of thin accretion disks, introduced by Kluźniak and Kita [Kluźniak, W., Kita, D., 2000. Three-dimensional structure of an alpha accretion disk. Available from: <arXiv:astro-ph/0006266v1> (KK)] for a steady-state disk flow, is extended to a time-dependent problem. Transient growth of axisymmetric disturbances is analytically shown to occur on the global disk scale. The implications of this result on the theory of hydrodynamical thin accretion disks, as well as future prospects, are discussed.     

Towards a new standard model for black hole accretion

We briefly review recent developments in black hole accretion disk theory, emphasizing the vital role played by magnetohydrodynamic (MHD) stresses in transporting angular momentum. The apparent universality of accretion-related outflow phenomena is a strong indicator that large-scale MHD torques facilitate vertical transport of angular momentum. This leads to an enhanced overall rate of angular momentum transport and allows accretion of matter to proceed at an interesting rate. Furthermore, we argue that when vertical transport is important, the radial structure of the accretion disk is modified at small radii and this affects the disk emission spectrum. We present a simple model demonstrating how energetic, magnetically-driven outflows modify the emergent disk emission spectrum with respect to that predicted by standard accretion disk theory. A comparison of the predicted spectra against observations of quasar spectral energy distributions suggests that mass accretion rates inferred using the standard disk model may be severely underestimated.     

Linkage between accretion disks and blazars

The magnetic field in an accretion disk is estimated assuming that all of the angular momentum within prescribed accretion disk radii is removed by a jet. The magnetic field estimated at the base of the jet is extrapolated to the blazar emission region using a model for a relativistic axisymmetric jet combined with some simplifying assumptions based on the relativistic nature of the flow. The extrapolated magnetic field is compared with estimates based upon the synchrotron and inverse Compton emission from three blazars, MKN 501, MKN 421 and PKS 2155-304. The magnetic fields evaluated from pure synchrotron self-Compton models are inconsistent with the magnetic fields extrapolated in this way. However, in two cases inverse Compton models in which a substantial part of the soft photon field is generated locally agree well, mainly because these models imply magnetic field strengths consistent with an important Poynting Flux component. This comparison is based on estimating the mass accretion rate from the jet energy flux. Further comparisons along these lines will be facilitated by independent estimates of the mass accretion rate in blazars and by more detailed models for jet propagation near the black hole.     

General relativistic MHD simulations of black hole accretion disks and jets

Accretion disks orbiting black holes power high-energy systems such as X-ray binaries and Active Galactic Nuclei. Observations are providing increasingly detailed quantitative information about such systems. This data has been interpreted using standard toy-models that rely on simplifying assumptions such as regular flow geometry and a parameterized stress. Global numerical simulations offer a way to investigate the basic physical dynamics of accretion flows without these assumptions and, in principle, lead to a genuinely predictive theory. In recent years we have developed a fully three-dimensional general relativistic magnetohydrodynamic simulation code that evolves time-dependent inflows into Kerr black holes. Although the resulting global simulations of black hole accretion are still somewhat simplified, they have brought to light a number of interesting results. These include the formation of electro-magnetically dominated jets powered by the black hole’s rotation, and the presence of strong stresses in the plunging region of the accretion flow. The observational consequences of these features are gradually being examined. Increasing computer power and increasingly sophisticated algorithms promise a bright future for the computational approach to black hole accretion.     

Relativistic reflection X-ray spectra of accretion disks

We have calculated the relativistic reflection component of the X-ray spectra of accretion disks in active galactic nuclei (AGN). Our calculations have shown that the spectra can be significantly modified by the motion of the accretion flow, and the gravity and rotation of the central black hole. The absorption edges in the spectra suffer severe en- ergy shifts and smearing, and the degree of distortion depends on the system parameters, in particular, the inner radius of the accretion disk and the disk viewing inclination angles. The effects are significant. Fluorescent X-ray emission lines from the inner accretion disk could be a powerful diagnostic of space-time distortion and dynamical relativistic effects near the event horizons of accreting black holes. However, improper treatment of the re- flection component in fitting the X-ray continuum could give rise to spurious line-like features. These features mimic the true fluorescent emission lines and may mask their relativistic signatures. Fully relativistic models for reflection continua together with the emission lines are needed in order to extract black-hole parameters from the AGN X-ray spectra.     

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary disks. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface. Magnetically-driven mixing has implications for disk chemistry and evolution of the grain population, and the effective viscous response of the disk determines whether planets migrate inwards or outwards. However, the weak ionisation of protoplanetary disks means that magnetic fields may not be able to effectively couple to the matter. I examine the magnetic diffusivity in a minimum solar nebula model and present calculations of the ionisation equilibrium and magnetic diffusivity as a function of height from the disk midplane at radii of 1 and 5 AU. Dust grains tend to suppress magnetic coupling by soaking up electrons and ions from the gas phase and reducing the conductivity of the gas by many orders of magnitude. However, once grains have grown to a few microns in size their effect starts to wane and magnetic fields can begin to couple to the gas even at the disk midplane. Because ions are generally decoupled from the magnetic field by neutral collisions while electrons are not, the Hall effect tends to dominate the diffusion of the magnetic field when it is able to partially couple to the gas, except at the disk surfaces where the low density of neutrals permits the ions to remain attached to the field lines. For a standard population of 0.1 μm grains the active surface layers have a combined column Σ_active≈2 g cm⁻² at 1 AU; by the time grains have aggregated to 3 μm, Σ_active≈80 g cm⁻². Ionisation in the active layers is dominated by stellar X-rays. In the absence of grains, X-rays maintain magnetic coupling to 10% of the disk material at 1 AU (i.e. Σ_active≈150 g cm⁻²). At 5 AU the Σ_active≈Σ_total once grains have aggregated to 1 μm in size.     

Angular momentum transport by thermal emission in black hole accretion disks

We calculate the amount of angular momentum that thermal photons carry out of a viscous black hole accretion disk, due to the strong Doppler shift imparted to them by the high orbital velocity of the radiating disk material. While thermal emission can not drive accretion on its own, we show that along with disk heating it does nonetheless result in a loss of specific angular momentum, thereby contributing to an otherwise viscosity‐driven accretion flow. In particular, we show that the fraction of the angular momentum that is lost to thermal emission at a radius r in a standard, multi‐color disk is ∼0.4r_s/r, where r_s is the Schwarzschild radius of the black hole. We briefly highlight the key similarties between this effect and the closely related Poynting‐Robertson effect (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dissipative-acoustic instability in accretion disks at a nonlinear stage

Nonstationary hydrodynamic models of a viscous accretion disk around a central compact object were constructed. Two different numerical methods (TVD and SPH) are used to study the dynamics of dissipatively unstable acoustic perturbations at the nonlinear stage in terms of the standard α-disk model. The standard disk accretion in the Shakura-Sunyaev model is unstable against acoustic waves for various parameters of the system. If the α parameter, which specifies the level of turbulent viscosity, exceeds α?0.03, then a complex nonstationary system of small-scale weak shock waves is formed. The growth rate of the perturbations is higher in the central disk region. For α?0.2, the relative shock amplitude can exceed 50% of the equilibrium disk parameters. The reflection of waves from the disk boundaries and their nonlinear interaction are important factors that can produce unsteady accretion. The luminosity of such a disk undergoes quasi-periodic oscillations at a level of several percent (?5%) of the equilibrium level.     

Current‐driven instabilities in weakly ionized disks

Cool weakly ionized gaseous rotating disk form the basis for many models in astrophysics objects. Instabilities against perturbations in such disks play an important role in the theory of the formation of stars and planets. Traditionally, axisymmetric magnetohydrodynamic (MHD) and recently Hall‐MHD instabilities have been thoroughly studied as providers of an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, the Hall instability against axisymmetric perturbations in incompressible rotating fluid in external poloidal and toroidal magnetic field is considered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic activity in accretion disc boundary layers

We use three-dimensional magnetohydrodynamic simulations to study the structure of the boundary layer between an accretion disc and a non-rotating, unmagnetized star. Under the assumption that cooling is efficient, we obtain a narrow but highly variable transition region in which the radial velocity is only a small fraction of the sound speed. A large fraction of the energy dissipation occurs in high-density gas adjacent to the hydrostatic stellar envelope, and may therefore be reprocessed and largely hidden from view of the observer. As suggested by Pringle , the magnetic field energy in the boundary layer is strongly amplified by shear, and exceeds that in the disc by an order of magnitude. These fields may play a role in generating the magnetic activity, X-ray emission and outflows in disc systems where the accretion rate is high enough to overwhelm the stellar magnetosphere.     

The evolution and efficiency of energy release of magnetized black-hole accretion disks

On the basis of the equivalent circuit model, we investigate two different mechanisms of extracting energy of rotation and angular momentum from a black hole by magnetic field, namely, the Blandford-Znajek (BZ) process and the magnetic coupling (MC) process. The contributions to the efficiency of energy release via pure accretion process, BZ process and MC process are compared in detail by studying the evolutionary characteristics of the spin parameter of the black hole at the center of the magnetized accretion disk. It is shown that the MC process is an important new mechanism of extracting energy from the rotating black hole and its efficiency of energy release is almost as high as that of the BZ process. The efficiency of energy release via pure accretion process is higher than those of BZ process and magnetic coupling process. However, when the rotation of a black hole approaches that of an extreme Kerr black hole,the efficiency of energy release is mainly due to the contributions of BZ process and MC process.     

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.     

Magnetic fields of accretion disks and outflows in prograde and retrograde black holes

We show that for the accretion disk with equipartition between magnetic and radiative pressures, prograde black holes generate outflowing energy in jets more efficiently than retrograde black holes do. Both viscous radiative and irradiative disks provide more efficient outflow jets in the case of a prograde black hole than in the case of a retrograde black hole. Our results confirm the conclusion of Tchekhovskoy & McKinney (2012) that, for the same absolute value of the spin, prograde black holes with geometrically thick accretion disks generate outflows several times more efficiently than retrograde black holes do. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)