The stability of accretion tori - I. Long-wavelength modes of slender tori

Abstract

We elucidate the inviscid instabilities of an iesntropic torus found previously by Papaloizou and Pringle. The torus is a polytrope of index n and has a small ratio of minor radius a to orbital radius r0. In equilibrium it rotates on cylinders with angular velocity profile omega ~ r^-q. Linear modes are proportional to exp[i(m phi - w t)]. For small beta = m a / r0, we justify the use of height-averaged equations by appealing to approximate vertical hydrostatic equilibrium. The effective polytropic index for the resulting two-dimensional problem is N = n + 1/2; Thus an incompressible torus in three dimensions behaves compressibly in two. Height-averaged modes obey an ordinary differential equation which we solve numerically to obtain the growth rate as a function of q, n, and beta. The error made in predicting the growth rate of the actual three-dimensional system is small everywhere along the principal branch even for beta ~ .5, and is less than 1 percent for the fastest-growing mode. We analytically solve the artificial case n=-1/2 which is two-dimensionally incompressible, and show that is has all the qualitative features of the general case, except that it does not have a resonance at corotation. In the general case, with n > -1/2 and q < 2, the corotation resonance absorbs energy and angular momentum, so the growing and decaying modes do not occur in complex conjugate pairs. We solve a second special case, namely n=2-q=0, almost analytically in three dimensions, without height averaging; Papaloizou and Pringle asserted that this system is stable but we show that there is an unstable mode for small beta just as in the other systems. In fact this principal unstable branch, with corotation at the pressure maximum is qualitatively the same for all n and is essentially independent both of compressibility and the vortensity. Thus the modes are not sonic, nor are they similar to those of the Kelvin-Helmholtz instability. Instead they are composed of two edge waves, akin to surface waves in water although modified by shear and rotation, coupled across a forbidden region around corotation.