The Radial-Azimuthal Instability of Gas-Pressure-Dominated Accretion Disk with Advection

The radial-azimuthal instability of gas-pressure-dominated accretion disk with advection is examined in this paper. We find that the including of very little advection has significant effects on two acoustic modes, which are no longer complex conjugates of each other. They increase the instability of the O-mode and damp that of the I-mode. We also find that when the azimuthal perturbations are considered, the stability properties of disk are different from that in pure radial perturbation case. The increase of azimuthal wave number will stabilize the acoustic modes but make the viscous mode more unstable and does not change the thermal mode very much for optically thin disk. The I-mode is more stable. The O-mode, viscous mode and thermal mode tend to become more unstable with the increase of azimuthal perturbation wavenumber for optically thick disk. For a geometrically slim, advection-dominated disk, the increasing of azimuthal perturbations make thermal mode more unstable and acoustic mode more stable. This revised version was published online in July 2006 with corrections to the Cover Date.     

Local Stability Criterion for a Self-Gravitating Gaseous Disk of Spiral Galaxies

An improved linear stability theory of small-amplitude oscillations of a self-gravitating, infinitesimally thin gaseous disk of spiral galaxies has been developed by Bertin, Lau, Lin, Mark, Morozov, Polyachenko, and others in the approximation of moderately tightly wound gravity perturbations. In this regime, the generalized Lin–Shu type dispersion relation was also found by including higher order terms in the small parameter 1/kr for wavenumber k and radius r. It was shown that in the differentially rotating disks for nonaxisymmetric (spiral)perturbations Toomre's modified critical Q-parameter is larger than the standard one: the fact that the spiral perturbations in the nonuniformly rotating system are more unstable than the axisymmetric ones is taken into account in this modified local stability criterion. We use hydrodynamical simulations to test the validity of the modified local criterion. This revised version was published online in July 2006 with corrections to the Cover Date.     

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)     

含粘滞和磁场的薄吸积盘不稳定性

本文采用微扰方法导出色散方程，并在四种情况下详细讨论了薄吸积盘的不稳定性。结果表明：在纯粘滞和纯磁场盘中都存在脉动不稳定性。而且在吸积盘内同时考虑粘滞和磁场时，存在两种不稳定性，一种是脉动不稳定性，另一种是单调不稳定性。同时数值计算还表明，脉动不稳定性更可能存在于盘的内区，而单调不稳定性则只在盘的外区，对短波扰动才有意义。这些结果为解释ＢＬＬａｃ天体、Ｓｅｙｆｅｒｔ星系及类星体等活动星系核的光变现象进一步提供了理论依据。     

Non-axisymmetric instabilities in shocked adiabatic accretion flows

We investigate the linear stability of a shocked accretion flow on to a black hole in the adiabatic limit. Our linear analyses and numerical calculations show that, despite the post-shock deceleration, the shock is generally unstable to non-axisymmetric perturbations. The simulation results of Molteni, Tóth & Kuznetsov can be well explained by our linear eigenmodes. The mechanism of this instability is confirmed to be based on the cycle of acoustic waves between the corotation radius and the shock. We obtain an analytical formula to calculate the oscillation period from the physical parameters of the flow. We argue that the quasi-periodic oscillation should be a common phenomenon in accretion flows with angular momentum.     

Negative Energy And Angular Momentum Modes Of Thin Accretion Disks

This work derives the linearized equations of motion, the Lagrangian density, the Hamiltonian density, and the canonical angular momentum density for general perturbations [∝ exp (imφ) with m = 0, ± 1, ...] of a geometrically thin self-gravitating, homentropic fluid disk including the pressure. The theory is applied to “eccentric,” m = ± 1 perturbations of a geometrically thin Keplerian disk. We find m = 1 modes at low frequencies relative to the Keplerian frequency. Further, it is shown that these modes can have negative energy and negative angular momentum. The radial propagation of these low-frequency m = 1 modes can transport angular momentum away from the inner region of a disk and thus increase the rate of mass accretion. Depending on the radial boundary conditions there can be discrete low-frequency, negative-energy, m = 1 modes.     

Influence of the Disk Wind on the Intrinsic Polarization of Young Stars

The behavior of the linear polarization parameters of UX Ori stars during their eclipses by circumstellar dust clouds is studied. A circumstellar disk with a disk wind creating a puffing in the dust sublimation zone is considered. We show that the disk puffing can strongly affect the degree of polarization and color index of the star during its eclipse. A strong wind can change the orientation of the plane of linear polarization. The scattered radiation from a thin disk is polarized perpendicularly to its plane, but the radiation from a disk with a strong wind can be polarized along the disk plane. A situation where the disk-scattered radiation is not polarized in a certain spectral band is possible owing to the disk puffing. There can be different orientations of the linear polarization of the disk radiation in different spectral bands.

     

Origin of the Cartwheel Galaxy: Disk Instability?

The linear theory and N-body simulations are used to present a new, alternative model of the galaxy A0035-324 (the “Cartwheel”), which is the most striking example of the relatively small class of ring galaxies. The model is based on the gravitational Jeans-type instability of both axisymmetric (radial) and nonaxisymmetric (spiral) small-amplitude gravity perturbations (e.g., those produced by spontaneous disturbances) of a dynamically cold subsystem (identified as the gaseous component) of an isolated disk galaxy. The simplified model of a galaxy is used in which stars (and a dark matter, if it exists at all) do not participate in the disk collective oscillations and just form a background charge. In the theory presented here, a case for both purely radial solutions and purely spiral solutions to the equations of motion of an infinitesimally thin gaseous disk is made, which is associated with both a radial density wave and a dominant spiral density wave which propagate outwards creating a rough ring and a number of spiral arms. Through three-dimensional numerical simulation of a collisionless set of many particles, I associate these gravitationally unstable axisymmetric waves and nonaxisymmetric waves with growing clumps of matter which take on the appearance of a ring and spokes of mass blobs.     

General relativistic analysis of structure and stability of charged fluid disks around compact objects

In this paper we give a detailed general relativistic formulation of the study of structure and stability of charged fluid disks around compact objects like black holes neglecting the self-gravitation of the disk itself. Having presented the general equations for equilibrium as well as for perturbations we solve explicitly the cases of rigidly and differentially rotating thin disks, with constant charge density and zero pressure, confined to the equatorial plane of the black hole. By using normal mode analysis we have analysed the stability of such disks under purely radial perturbations and find that the disks are generally stable. On leave of absence from Government College, Jagadalpur 494005     

Magneto-Acoustic Waves in Compressible Magnetically Twisted Flux Tubes

The oscillatory modes of a magnetically twisted compressible flux tube embedded in a compressible magnetic environment are investigated in cylindrical geometry. Solutions to the governing equations to linear wave perturbations are derived in terms of Whittaker’s functions. A general dispersion equation is obtained in terms of Kummer’s functions for the approximation of weak and uniform internal twist, which is a good initial working model for flux tubes in solar applications. The sausage, kink and fluting modes are examined by means of the derived exact dispersion equation. The solutions of this general dispersion equation are found numerically under plasma conditions representative of the solar photosphere and corona. Solutions for the phase speed of the allowed eigenmodes are obtained for a range of wavenumbers and varying magnetic twist. Our results generalise previous classical and widely applied studies of MHD waves and oscillations in magnetic loops without a magnetic twist. Potential applications to solar magneto-seismology are discussed.     

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.     

The equilibrium magnetized disks

Analytical solutions are obtained for gravitational and magnetic fields at given distributions of area matter density and surface currents.Numerical solutions ae obtained for a magnetized non-rotating disk, the equilibrium of which is governed by the balance of magnetic and gravitational forces. The models are shown to be unsable against short-wave perturbations. Pressure can make a small contribution (of the order of a few percent) towards stabilizng the disk against short-wave perturbations. Such disks may appear in some galactic nuclei.     

Linear kinetic theory of instabilities of a gravitatingstellar disk

Linear kinetic theory is developed to describe collective oscillations (and their instabilities) propagating in a rapidly rotating disk of stars, representing a highly flattened galaxy. The analysis is carried out for the special case of a self-gravitating, infinitesimally thin, and spatially inhomogeneous system, taking into account the effects both of thermal movements of stars and of gravitational encounters between stars and giant molecular clouds of an interstellar medium. The star–cloud encounters are described with the use of the Landau collision integral. The dynamics of gravity perturbations with rare interparticle encounters is considered. Such a disk is treated by employing the well elaborated mathematical formalisms from plasma perturbation theory using normal-mode analysis. In particular, the method of solving the Boltzmann equation is applied by integration along paths, neglecting the influence of star–cloud encounters on the distribution of stars in the zeroth-order approximation. We are especially interested in important kinetic effects due to wave–star resonances, which we have little knowledge about. The kinetic effects are introduced via a minor drift motion of stars which is computed from the equations of stellar motion in an unperturbed central force field of a galaxy. The dispersion laws for two main branches of disk's oscillations, that is the classical Jeans branch and an additional gradient branch, are deduced. The resonant Landau-type instabilities of hydrodynamically stable Jeans and gradient gravity perturbations is considered to be a long-term generating mechanism for propagating density waves, thereby leading to spiral-like and/or ring-like patterns in the flat galaxies. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bending instability of stellar disks: The stabilizing effect of a compact bulge

The saturation conditions for bending modes in inhomogeneous thin stellar disks that follow from an analysis of the dispersion relation are compared with those derived from N-body simulations. In the central regions of inhomogeneous disks, the reserve of disk strength against the growth of bending instability is smaller than that for a homogeneous layer. The spheroidal component (a dark halo, a bulge) is shown to have a stabilizing effect. The latter turns out to depend not only on the total mass of the spherical component, but also on the degree of mass concentration toward the center. We conclude that the presence of a compact (not necessarily massive) bulge in spiral galaxies may prove to be enough to suppress the bending perturbations that increase the disk thickness. This conclusion is corroborated by our N-body simulations in which we simulated the evolution of near-equilibrium, but unstable finite-thickness disks in the presence of spheroidal components. The final disk thickness at the same total mass of the spherical component (dark halo + bulge) was found to be much smaller than that in the simulations where a concentrated bulge is present.     

Detecting ionized bubbles in redshifted 21-cm maps

Linear transient phenomena induced by flow non-normality in thin self-gravitating astrophysical discs are studied using the shearing sheet approximation. The considered system includes two modes of perturbations: vortex and (spiral density) wave. It is shown that self-gravity considerably alters the vortex mode dynamics; its transient (swing) growth may be several orders of magnitude stronger than in the non-self-gravitating case and two to three times larger than the transient growth of the wave mode. Based on this finding, we comment on the role of vortex mode perturbations in a gravitoturbulent state. We also describe the linear coupling of the perturbation modes, caused by the differential character of disc rotation. The coupling is asymmetric: vortex mode perturbations are able to excite wave mode perturbations, but not vice versa. This asymmetric coupling lends additional significance to the vortex mode as a participant in spiral density waves and shock manifestations in astrophysical discs.     

Nonlocal density wave theory for gravitational instability of protoplanetary disks without sharp boundaries

The stability of a self‐gravitating infinitesimally thin gaseous disk rotating around a central mass is studied. Our global linear analysis concerns marginal stability, i.e. it yields the critical temperature for the onset of instability for any given ratio of the disk mass to the central mass. Both axisymmetric and low‐m nonaxisymmetric excitations are analysed. When the fractional disk mass increases, the symmetry character of the instability changes from rings (m = 0) to one‐armed trailing spirals (m = 1). The distribution of the surface density along the spiral arms is not uniform, but describes a sequence of maxima that might be identified with forming planets. The number of the mass concentrations decreases with increasing fractional disk mass. We also obtain solutions in the form of global nonaxisymmetric vortices, which are, however, never excited.     

Stability of magnetic flux tubes in close binary stars

Photometric and Doppler imaging observations of active binaries indicate the existence of starspots at preferred longitudes (position angles with respect to the companion star). We investigate the stability of magnetic flux tubes in the convection zone of close, fast‐rotating binary stars and explore whether the observed preferred longitudes could be caused by tidal forces and the deformation of the active star. We assume a synchronized binary system with spin axes perpendicular to the orbital plane and a rotation period of a few days. The tidal force and the deviation from spherical structure are considered in lowest‐order perturbation theory. The magnetic field is in the form of toroidal magnetic flux rings, which are stored in mechanical equilibrium within the stably stratified overshoot region beneath the convection zone until the field has grown sufficiently strong for the undulatory instability to initiate the formation of rising loops. Frequencies and geometry of stable as well as growth rates of unstable eigenmodes are determined by linear stability analysis. Particular consideration is given to the question whether the effects of tidal forces and perturbations of the stellar structure can force a rising flux loop to enter the convection zone at specific longitudes.     

Evolutionary processes in protoplanetary accretion disks: the propagation of axisymmetric shock waves

In this paper we investigate both the global and the local hydrodynamics of axisymmetric accretion disks around young stellar objects under the simultaneous action of viscosity, self-gravity and pressure forces. For simplicity, we take for the global model a polytropic equation of state, make the infinitely thin disk approximation and characterize the surface density and temperature profiles in the disk as power laws in the radial distance r from the protostar. We solve the problem of the general density profile of a Keplerian disk showing that self-gravity could not be an important factor for the fast formation of the rocky cores of giant gaseous planets in our solar system. Under the hypothesis that the unperturbed rotation curve of the disk is nearly Keplerian throughout the radial extent, we can estimate with our polytropic model a lower limit for the resulting masses M_d(r) of stable disks up to 100 AU. These masses are in the range of the so-called minimum mass solar nebular (_d/M_s ≈ 0.01–0.02).By adopting a simplified viscosity model, where the height-integrated turbulent dynamical viscosity ν is a function of the surface density σ like η ∝ σ^Γ, we derive in the local shearing sheet model linearized evolution equations for small density perturbations describing both a diffusion process and the propagation of acoustic density waves. We solve a special initial value problem and calculate the appropriate Green's function. The analytical solutions so obtained describe in the case Γ < 0 the successive formation of quasi-stationary ring-shaped density structures in a disk with a definite mode of maximum instability, whereas in the case Γ > Γ_c the density wave equation describes the propagation of an “overstable” ring-shaped acoustic density wavelet to the outer ranges of the accretion disk. Whereas the group velocity of the wave packet is subsonic, the phase velocities of individual wave crests in the wave packet are supersonic. The mode of maximum instability, the growth rate and the number of growing waves in the wavelet are controlled by Γ and α. Our present knowledge concerning turbulent viscosity in protoplanetary disks is not sufficient to decide whether or not the case Γ > Γ_c is realized.The suggested structuring processes in the linear theory should initiate in the non-linear regime the formation of narrow ring-shaped density shock waves moving through the protoplanetary disk. These non-linear waves could produce extremely spatially and temporally heterogeneous temperature regions in the disk. We speculate that ring-shaped density waves, excited by inner boundary conditions and which have dominated the disk's evolution at early times, are responsible both for the fast growth of dust to planetesimals and at least for the rapid accretion of the rocky cores of giant gaseous planets in the protoplanetary accretion disk (shock wave trigger hypothesis). We derive provisional scaling rules for planetary systems regarding the spacing of orbits as a function of the mass ratio of the protoplanetary disk to the protostar. However, further analytical work and linear as well as nonlinear numerical simulations of density waves excited by inner boundary conditions are needed to consolidate the results and speculations of our linear wave mechanics in the future.     

On the Hall Instability in Protostellar Disks

Astronomy Letters - Small perturbations of a protostellar disk with vertical and azimuthal magnetic field components are considered in terms of Hall magnetohydrodynamics. The dispersion equation...