Magnetohydrodynamic turbulence in warped accretion discs

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.Comment: to appear in the proceedings of the 20th Texas Symposium on Relativistic Astrophysics, J. C. Wheeler & H. Martel (eds.

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $\alpha_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.Comment: 12 pages, 18 figures, MNRAS accepte

Possible potentials responsible for stable circular relativistic orbits

Bertrand's theorem in classical mechanics of the central force fields attracts us because of its predictive power. It categorically proves that there can only be two types of forces which can produce stable, circular orbits. In the present article an attempt has been made to generalize Bertrand's theorem to the central force problem of relativistic systems. The stability criterion for potentials which can produce stable, circular orbits in the relativistic central force problem has been deduced and a general solution of it is presented in the article. It is seen that the inverse square law passes the relativistic test but the kind of force required for simple harmonic motion does not. Special relativistic effects do not allow stable, circular orbits in presence of a force which is proportional to the negative of the displacement of the particle from the potential center.Comment: 11 pages, Latex fil

Hypercritical Advection Dominated Accretion Flow

In this note we study the accretion disc that arises in hypercritical accretion of $\dot M\sim 10^8 M_{\rm Edd}$ onto a neutron star while it is in common envelope evolution with a massive companion. In order to raise the temperature high enough that the disc might cool by neutrino emission, Chevalier found a small value of the $\alpha$-parameter, where the kinematic coefficient of shear viscosity is $\nu=\alpha c_s H$, with $c_s$ the velocity of sound and $H$ the disc height; namely, $\alpha\sim 10^{-6}$ was necessary for gas pressure to dominate. He also considered results with higher values of $\alpha$, pointing out that radiation pressure would then predominate. With these larger $\alpha$'s, the temperatures of the accreting material are much lower, \lsim 0.35 MeV. The result is that neutrino cooling during the flow is negligible, satisfying very well the advection dominating conditions. The low temperature of the accreting material means that it cannot get rid of its energy rapidly by neutrino emission, so it piles up, pushing its way through the accretion disc. An accretion shock is formed, far beyond the neutron star, at a radius \gsim 10^8 cm, much as in the earlier spherically symmetric calculation, but in rotation. Two-dimensional numerical simulation shows that an accretion disc is reformed inside of the accretion shock, allowing matter to accrete onto the neutron star with pressure high enough so that neutrinos can carry off the energy.Comment: 6 pages, ApJ, submitte

Radio emission from the high-mass X-ray binary BP Cru: first detection

BP Cru is a well known high-mass X-ray binary composed of a late B hypergiant (Wray 977) and a neutron star, also observed as the X-ray pulsar GX 301-2. No information about emission from BP Cru in other bands than X-rays and optical has been reported to date in the literature, though massive X-ray binaries containing black holes can have radio emission from a jet. In order to assess the presence of a radio jet, we searched for radio emission towards BP Cru using the Australia Compact Array Telescope during a survey for radio emission from Be/X-ray transients. We probed the 41.5d orbit of BP Cru with the Australia Telescope Compact Array not only close to periastron but also close to apastron. BP Cru was clearly detected in our data on 4, possibly 6, of 12 occasions at 4.8 and 8.6 GHz. Our data suggest that the spectral index of the radio emission is modulated either by the X-ray flux or the orbital phase of the system. We propose that the radio emission of BP Cru probably arises from two components: a persistent component, coming from the mass donor Wray 977, and a periodic component connected to the accretion onto the neutron star, possibly coming from a (weak and short lived) jet.Comment: 2 figures, accepted for publication in A+A letter

An Incoherent α−Ω\alpha-\Omega Dynamo in Accretion Disks

We use the mean-field dynamo equations to show that an incoherent alpha effect in mirror-symmetric turbulence in a shearing flow can generate a large scale, coherent magnetic field. We illustrate this effect with simulations of a few simple systems. In accretion disks, this process can lead to axisymmetric magnetic domains whose radial and vertical dimensions will be comparable to the disk height. This process may be responsible for observations of dynamo activity seen in simulations of dynamo-generated turbulence involving, for example, the Balbus-Hawley instability. In this case the magnetic field strength will saturate at $\sim (h/r)^2$ times the ambient pressure in real accretion disks. The resultant dimensionless viscosity will be of the same order. In numerical simulations the azimuthal extent of the simulated annulus should be substituted for $r$. We compare the predictions of this model to numerical simulations previously reported by Brandenburg et al. (1995). In a radiation pressure dominated environment this estimate for viscosity should be reduced by a factor of $(P_{gas}/P_{radiation})^6$ due to magnetic buoyancy.Comment: 23 pages, uses aaste

Oscillations of tori in the pseudo-Newtonian potential

Context. The high-frequency quasi-periodic oscillations (HF QPOs) in neutron star and stellar-mass black hole X-ray binaries may be the result of a resonance between the radial and vertical epicyclic oscillations in strong gravity. Aims. In this paper we investigate the resonant coupling between the epicyclic modes in a torus in a strong gravitational field. Methods. We perform numerical simulations of axisymmetric constant angular momentum tori in the pseudo-Newtonian potential. The epicyclic motion is excited by adding a constant radial velocity to the torus. Results. We verify that slender tori perform epicyclic motions at the frequencies of free particles, but the epicyclic frequencies decrease as the tori grow thicker. More importantly, and in contrast to previous numerical studies, we do not find a coupling between the radial and vertical epicyclic motions. The appearance of other modes than the radial epicyclic motion in our simulations is rather due to small numerical deviations from exact equilibrium in the initial state of our torus. Conclusions. We find that there is no pressure coupling between the two axisymmetric epicyclic modes as long as the torus is symmetric with respect to the equatorial plane. However we also find that there are other modes in the disc that may be more attractive for explaining the HF QPOs.Comment: 8 pages, 9 figure

Propagation of non-linear circularly polarised Alfven waves in a homogeneous medium

We study the evolution of non-linear circularly polarised Alfven waves by solving numerically the time-dependent equations of magnetohydrodynamics (MHD) in one dimension. We examine the behaviour of the waves and find that different physical mechanisms are relevant in different ranges of beta. In a low beta plasma the wave may undergo a parametric decay. This is because the wave excites a density enhancement that travels slower than the wave itself and thus interacts with the wave. When beta is greater or equal to one the density enhancement does not interact with the wave and no decay takes place, instead the Alfven wave is reflected against the density enhancement. The reflection zone propagates with the speed 1/n of the Alfven speed. Because of that the magnetic flux is conserved which results in an amplification of the oscillating magnetic field by a factor 1/n. We find that n depends on beta, and that in particular it is less or equal to one for values of beta ~ 1 and is greater or equal to one for beta >> 1. We discuss the relevance of these mechanisms to the acceleration of the solar wind, and the triggering of MHD turbulence in the polar wind region. In particular these simulations can explain the presence of inward propagating Alfven waves in the solar corona

An alpha theory of time-dependent warped accretion discs

The non-linear fluid dynamics of a warped accretion disc was investigated in an earlier paper by developing a theory of fully non-linear bending waves in a thin, viscous disc. That analysis is here extended to take proper account of thermal and radiative effects by solving an energy equation that includes viscous dissipation and radiative transport. The problem is reduced to simple one-dimensional evolutionary equations for mass and angular momentum, expressed in physical units and suitable for direct application. This result constitutes a logical generalization of the alpha theory of Shakura & Sunyaev to the case of a time-dependent warped accretion disc. The local thermal-viscous stability of such a disc is also investigated.Comment: 16 pages, 3 figures, to be published in MNRA