On the Rayleigh–Taylor instability of radio bubbles in galaxy clusters

We consider the Rayleigh–Taylor instability in the early evolution of the rarefied radio bubbles (cavities) observed in many cooling-flow clusters of galaxies. The top of a bubble becomes prone to the Rayleigh–Taylor instability as the bubble rises through the intracluster medium (ICM). We show that while the jet is powering the inflation, the deceleration of the bubble–ICM interface is able to reverse the Rayleigh–Taylor instability criterion. In addition, the inflation introduces a drag effect which increases substantially the instability growth time. The combined action of these two effects considerably delays the onset of the instability. Later on, when the magnitude of the deceleration drops or the jet fades, the Rayleigh–Taylor and the Kelvin–Helmholtz instabilities set in and eventually disrupt the bubble. We conclude that the initial deceleration and drag, albeit unable to prevent the disruption of a bubble, may significantly lengthen its lifetime, removing the need to invoke stabilizing magnetic fields.     

Accretion to magnetized stars through the Rayleigh–Taylor instability: global 3D simulations

We present results of 3D simulations of magnetohydrodynamics (MHD) instabilities at the accretion disc–magnetosphere boundary. The instability is Rayleigh–Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It manifests itself in the form of tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner disc dynamical time-scale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, Θ≲ 30°, between the star's rotation and magnetic axes, and is associated with higher accretion rates. The hotspots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems, as well as planet migration.     

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.     

Discovery of drifting high-frequency quasi-periodic oscillations in global simulations of magnetic boundary layers

We report on the numerical discovery of quasi-periodic oscillations (QPOs) associated with accretion through a non-axisymmetric magnetic boundary layer in the unstable regime, when two ordered equatorial streams form and rotate synchronously at approximately the angular velocity of the inner disc. The streams hit the star's surface producing hotspots. Rotation of the spots leads to high-frequency QPOs. We performed a number of simulation runs for different magnetospheric sizes from small to tiny, and observed a definite correlation between the inner disc radius and the QPO frequency: the frequency is higher when the magnetosphere is smaller. In the stable regime, a small magnetosphere forms and accretion through the usual funnel streams is observed, and the frequency of the star is expected to dominate the light curve. We performed exploratory investigations of the case in which the magnetosphere becomes negligibly small and the disc interacts with the star through an equatorial belt. We also performed investigation of somewhat larger magnetospheres where one or two ordered tongues may dominate over other chaotic tongues. In application to millisecond pulsars, we obtain QPO frequencies in the range of 350–990 Hz for one spot. The frequency associated with rotation of one spot may dominate if spots are not identical or antipodal. If the spots are similar and antipodal, then the frequencies are twice as high. We show that variation of the accretion rate leads to drift of the QPO peak.     

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.     

Linear and non-linear theory of a parametric instability of hydrodynamic warps in Keplerian discs

We consider the stability of warping modes in Keplerian discs. We find them to be parametrically unstable using two lines of attack, one based on three-mode couplings and the other on Floquet theory. We confirm the existence of the instability and investigate its non-linear development in three dimensions, via numerical experiment. The most rapidly growing non-axisymmetric disturbances are the most nearly axisymmetric (low- m ) ones. Finally, we offer a simple, somewhat speculative model for the interaction of the parametric instability with the warp. We apply this model to the masing disc in NGC 4258 and show that, provided the warp is not forced too strongly, parametric instability can fix the amplitude of the warp.     

MRI channel flows and their parasites

Local simulations of the magnetorotational instability (MRI) in accretion discs can exhibit recurrent coherent structures called channel flows. The formation and destruction of these structures may play a role in the development and saturation of MRI-induced turbulence, and consequently help us understand the time-dependent accretion behaviour of certain astrophysical objects. Previous investigations have revealed that channel solutions are attacked by various parasitic modes, foremost of which is an analogue of the Kelvin–Helmholtz instability. We revisit these instabilities and show how they relate to the classical instabilities of plasma physics, the kink and pinch modes. However, we argue that in most cases channels emerge from developed turbulence and are eventually destroyed by turbulent mixing, not by the parasites. The exceptions are clean isolated channels, which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures can achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.     

Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

Warp and eccentricity propagation in discs around black holes

We consider the inward propagation of warping and eccentric disturbances in discs around black holes under a wide variety of conditions. In our calculations, we use secular theories of warped and eccentric discs and assume the deformations to be stationary and propagating in a disc model similar to regions (a) and (b) of Shakura & Sunyaev discs. We find that the propagation of deformations to the innermost regions of the disc is facilitated for low viscous damping and high accretion rate. We relate our results to the possible excitation of trapped inertial modes, and to the observations of high-frequency quasi-periodic oscillations (QPOs) in black hole systems in the very high spectral state.     

Magnetic field dragging and the vertical structure of thin accretion discs

The presence of an imposed vertical magnetic field may drastically influence the structure of thin accretion discs. If the field is sufficiently strong, the rotation law can depart from the Keplerian one. We consider the structure of a disc for a given eddy magnetic diffusivity but neglect details of the energy transport. The magnetic field is assumed to be in balance with the internal energy of the accretion flow. The thickness of the disc as well as the turbulent magnetic Prandtl number and the viscosity, α , are the key parameters of our model. The calculations show that the radial velocity can reach the sound speed for a magnetic disc if the thickness is comparable to that of a non-magnetic one. This leads to a strong amplification of the accretion rate for a given surface density. The inclination angle of the magnetic field lines can exceed the critical value 30° (required to launch cold jets) even for a relatively small magnetic Prandtl number of order unity. The toroidal magnetic fields induced at the disc surface are smaller than predicted in previous studies.     

On an excitation mechanism for trapped inertial waves in discs around black holes

According to one model, high-frequency quasi-periodic oscillations (QPOs) can be identified with inertial waves, trapped in the inner regions of accretion discs around black holes due to relativistic effects. In order to be detected, their amplitudes need to reach large enough values via some excitation mechanism. We work out in detail a non-linear coupling mechanism suggested by Kato, in which a global warping or eccentricity of the disc has a fundamental role. These large-scale deformations combine with trapped modes to generate 'intermediate' waves of negative energy that are damped as they approach either their corotation resonance or the inner edge of the disc, resulting in amplification of the trapped waves. We determine the growth rates of the inertial modes, as well as their dependence on the spin of the black hole and the properties of the disc. Our results indicate that this coupling mechanism can provide an efficient excitation of trapped inertial waves, provided the global deformations reach the inner part of the disc with non-negligible amplitude.     

Relativistic and Newtonian diskoseismology

The normal mode oscillations of thin accretion disks around black holes and other compact objects are analyzed and contrasted with those in stars. For black holes, the most robust modes are gravitationally trapped near the radius at which the radial epicyclic frequency is maximum. Their eigenfrequencies depend mainly on the mass and angular momentum of the black hole. The fundamental g-mode has recently been seen in numerical simulations of black hole accretion disks. For stars such as white dwarfs, the modes are trapped near the inner boundary (magnetospheric or stellar) of the accretion disk. Their eigenfrequencies are approximately multiples of the (Keplerian) angular velocity of the inner edge of the disk. The relevance of these modes to the high frequency quasi-periodic oscillations observed in the power spectra of accreting binaries will be discussed. In contrast to most stellar oscillations, most of these modes are unstable in the presence of viscosity (if the turbulent viscosity induced by the magnetorotational instability acts hydrodynamically).     

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.     

Corotational damping of discoseismic c modes in black hole accretion discs

Discoseismic c modes in accretion discs have been invoked to explain low-frequency variabilities observed in black hole X-ray binaries. These modes are trapped in the innermost region of the disc and have frequencies much lower than the rotation frequency at the disc inner radius. We show that because the trapped waves can tunnel through the evanescent barrier to the corotational wave zone, the c modes are damped due to wave absorption at the corotation resonance. We calculate the corotational damping rates of various c modes using the Wentzel-Kramers-Brillouin (WKB) approximation. The damping rate varies widely depending on the mode frequency, the black hole spin parameter and the disc sound speed, and is generally much less than 10 per cent of the mode frequency. A sufficiently strong excitation mechanism is needed to overcome this corotational damping and make the mode observable.     

The role of Kelvin–Helmholtz instability in dusty and partially ionized outflows

We investigate the linear theory of Kelvin–Helmholtz instability at the interface between a partially ionized dusty outflow and the ambient material analytically. We model the interaction as a multifluid system in a planar geometry. The unstable modes are independent from the charge polarity of the dust particles. Although our results show a stabilizing effect for charged dust particles, the growth time-scale of the growing modes gradually becomes independent of the mass or charge of the dust particles when the magnetic-field strength increases. We show that growth time-scale decreases with increasing the magnetic field. Also, as the mass of the dust particles increases, the growth time-scale of the unstable mode increases.     

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.     

Dynamical instability of laminar axisymmetric flows of ideal incompressible fluid

The instability of axisymmetric flows of ideal incompressible fluid with respect to infinitesimal perturbations with the nonconservation of angular momentum is investigated by numerically integrating the differential equations of hydrodynamics. The problem has been solved for two types of rotation profiles of an unperturbed flow: with zero and nonzero pressure gradients at the flow boundaries. Both rigid and free boundary conditions have been considered. The stability of axisymmetric flows with free boundaries is of great importance in disk accretion problems. Our calculations have revealed a crucial role of the flow pattern near the boundaries in the instability of the entire main flow. When the pressure gradient at the boundaries is zero, there is such a limiting scale of perturbations in azimuthal coordinate that longer-wavelength perturbations grow, while growing shorter-wavelength perturbations do not exit. In addition, for a fixed radial flow extent, there exists a nonzero minimum amplitude of the deviation of the angular velocity from the Keplerian one at which the instability vanishes. For a nonzero pressure gradient at the boundaries, the flow is unstable with respect to perturbations of any scale and at any small deviation of the angular velocity from the Keplerian one.     

The effective temperature distribution of steady-state, mass-losing accretion discs

Mass loss appears to be a common phenomenon among astrophysical accretion disc systems. An outflow emanating from an accretion disc can act as a sink for mass, angular momentum and energy, and can therefore alter the dissipation rates and effective temperatures across the disc. Here, the radial distributions of dissipation rate and effective temperature across a Keplerian, steady-state, mass-losing accretion disc are derived, using a simple, parametric approach that is sufficiently general to be applicable to many types of dynamical disc–wind models.
Effective temperature distributions for mass-losing accretion discs in cataclysmic variables are shown explicitly, with parameters chosen to describe both radiation-driven and centrifugally driven outflows. For realistic wind mass-loss rates of a few per cent, only centrifugally driven outflows – particularly those in which mass loss is concentrated in the inner disc – are likely to alter the effective temperature distribution of the disc significantly. Accretion discs that drive such outflows could produce spectra and eclipse light curves that are noticeably different from those produced by standard, conservative discs.     

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.