Stability of power-law discs — II. The global spiral modes

This paper reports on the in-plane normal modes in the self-consistent and the cut-out power-law discs. Although the cut-out discs are remarkably stable to bisymmetric perturbations, they are very susceptible to one-armed modes. For this harmonic, there is no inner Lindblad resonance, thus removing a powerful stabilizing influence. A physical mechanism for the generation of the one-armed instabilities is put forward. Incoming trailing waves are reflected as leading waves at the inner cut-out, thus completing the feedback for the swing-amplifier. Growing three-armed and four-armed modes occur only at very low temperatures. However, neutral m  = 3 and m  = 4 modes are possible at higher temperatures for some discs. The rotation curve index β has a marked effect on stability. For all azimuthal wavenumbers, any unstable modes persist to higher temperatures and grow more vigorously if the rotation curve is rising (β < 0) than if the rotation curve is falling (β > 0). If the central regions or outer parts of the disc are carved out more abruptly, any instabilities become more virulent. The self-consistent power-law discs possess a number of unusual stability properties. There is no natural time-scale in the self-consistent disc. If a mode is admitted at some pattern speed and growth rate, then it must be present at all pattern speeds and growth rates. Our analysis — although falling short of a complete proof — suggests that such a two-dimensional continuum of non-axisymmetric modes does not occur and that the self-consistent power-law discs admit no global non-axisymmetric normal modes whatsoever. Without reflecting boundaries or cut-outs, there is no resonant cavity and no possibility of unstable growing modes. The self-consistent power-law discs certainly admit equi-angular spirals as neutral modes, together with a one-dimensional continuum of growing axisymmetric modes.     

Gravitationally unstable protoplanetary discs

The possibility that protoplanetary gaseous discs are dynamically unstable to axisymmetric and non-axisymmetric gravity perturbations (e.g. those produced by spontaneous disturbances) with characteristic scales larger than the vertical scale height is discussed analytically, using a local Wentzel–Kramers–Brillouin (WKB) approach. It is shown that such discs might be clumpy, and these gravitationally bound clumps may later collapse to become giant planets ('hot Jupiters'). The chief aim in this paper is to underscore a fact of vital importance for application in the planetary formation process: gravitationally unstable non-axisymmetric (spiral) perturbations can effectively transport both the angular momentum and the mass in a spatially inhomogeneous disc.     

Interface modes and their instabilities in accretion disc boundary layers

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the innermost regions of accretion discs or boundary layers. Two simple models are considered in this paper. The magnetosphere–disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star–disc model involves a Keplerian disc terminated at the stellar atmosphere with high density and small density scaleheight. We find that the interface modes at the magnetosphere–disc boundary are generally unstable due to Rayleigh–Taylor and/or Kelvin–Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh–Taylor instability, and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin–Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh–Taylor or Kelvin–Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.     

Reprocessed emission from warped accretion discs with application to X-ray iron line profiles

Fluorescent iron line profiles currently provide the best diagnostic for engine geometries of active galactic nuclei (AGN). Here we construct a method for calculating the relativistic iron line profile from an arbitrarily warped accretion disc, illuminated from above and below by hard X-ray sources. This substantially generalizes previous calculations of reprocessing by accretion discs by including non-axisymmetric effects. We include a relativistic treatment of shadowing by ray-tracing photon paths along Schwarzschild geodesics. We apply this method to two classes of warped discs, and generate a selection of resulting line profiles. New profile features include a time-varying line profile if the warp precesses about the disc, profile differences between 'twisted' and 'twist-free' warps and the possibility of steeper red and softer blue fall-offs than for flat discs. We discuss some qualitative implications of the line profiles in the context of Type I and II Seyfert AGN and other sources.     

Three-dimensional smoothed particle hydrodynamics simulations of radiation-driven warped accretion discs

We present three-dimensional smoothed particle hydrodynamics calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc, resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and, in general, the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios, using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios   q ∼ 0.1  a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (cf. X 1916−053). For less extreme mass ratios, the entire disc tilts out of the orbital plane (cf. Her X–1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction, we obtain both positive and negative superhumps simultaneously in the dissipation light curve (cf. V603 Aql).     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.     

Analytical potential–density pairs from complex-shifted Kuzmin–Toomre discs

The complex-shift method is applied to the Kuzmin–Toomre family of discs to generate a family of non-axisymmetric flat distributions of matter. These are then superposed to construct non-axisymmetric flat rings. We also consider triaxial potential–density pairs obtained from these non-axisymmetric flat systems by means of suitable transformations. The use of the imaginary part of complex-shifted potential–density pairs is also discussed.     

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.     

Tracing spiral density waves in M81

We use Spitzer IRAC 3.6 and 4.5 μm near-infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), optical B, V and I and Two-Micron All-Sky Survey K _s-band data to produce mass surface density maps of M81. The IRAC 3.6- and 4.5-μm data, whilst dominated by emission from old stellar populations, are corrected for small-scale contamination by young stars and polycyclic aromatic hydrocarbon emission. The I -band data are used to produce a mass surface density map by a   B − V   colour correction, following the method of Bell and de Jong. We fit a bulge and exponential disc to each mass map, and subtract these components to reveal the non-axisymmetric mass surface density. From the residual mass maps, we are able to extract the amplitude and phase of the density wave, using azimuthal profiles. The response of the gas is observed via dust emission in the 8-μm IRAC band, allowing a comparison between the phase of the stellar density wave and gas shock. The relationship between this angular offset and radius suggests that the spiral structure is reasonably long-lived and allows the position of corotation to be determined.     

The dynamical evolution of substructure

The evolution of substructure embedded in non-dissipative dark haloes is studied through N -body simulations of isolated systems, both in and out of initial equilibrium, complementing cosmological simulations of the growth of structure. We determine by both analytic calculations and direct analysis of the N -body simulations the relative importance of various dynamical processes acting on the clumps, such as the removal of material by global tides, clump–clump heating, clump–clump merging and dynamical friction. The ratio of the internal clump velocity dispersion to that of the dark halo is an important parameter; as this ratio approaches a value of unity, heating by close encounters between clumps becomes less important, while the other dynamical processes continue to increase in importance. Our comparison between merging and disruption processes implies that spiral galaxies cannot be formed in a protosystem that contains a few large clumps, but can be formed through the accretion of many small clumps; elliptical galaxies form in a more clumpy environment than do spiral galaxies. Our results support the idea that the central cusp in the density profiles of dark haloes is the consequence of self-limiting merging of small, dense haloes. This implies that the collapse of a system of clumps/substructure is not sufficient to form a cD galaxy, with an extended envelope; plausibly, subsequent accretion of large galaxies is required. The post-collapse system is in general triaxial, with rounder systems resulting from fewer, but more massive, clumps. Persistent streams of material from disrupted clumps can be found in the outer regions of the final system, and at an overdensity of around 0.75, can cover 10 to 30 per cent of the sky.     

Density cusps: restrictions on non-axisymmetric models

Galactic nuclei are now generally thought to have density cusps in their centres, and the evidence is mounting that as a consequence they are unlikely to be triaxial. Self-consistent stellar dynamical models of non-axisymmetric cusps would be an interesting counter-argument to this conclusion. We consider 2D analogues of triaxial cusps: a sequence of non-axisymmetric, cuspy discs first described by Sridhar & Touma. Scale-free models with potential Φ ∝  r ^α are examined in detail. It is shown analytically for 0 < α ≲ 0.43 that self-consistent models with positive phase-space density do not exist. Numerical solutions of the combined Vlasov and Poisson equations suggest that the whole sequence of models with 0 < α < 1 are also unphysical. Together with existing work on cusps, we conclude on purely theoretical grounds that galactic nuclei are not expected to be triaxial.     

Coupled galactic density-wave modes in a composite system of thin stellar and gaseous discs

A rotating disc galaxy is modelled as a composite system consisting of thin stellar and gaseous discs, which are described by a two-fluid modal formalism. The composite disc system is assumed to retain axisymmetry in the background equilibrium. General density-wave perturbations in the two discs are coupled through the mutual gravitational interaction. We study the basic properties of open and tight spiral density-wave modes in such a composite disc system. Within the Lindblad resonances, perturbation enhancements of surface mass density in stellar and gaseous discs are in phase; this is also true during the initial growth phase of density-wave perturbations. Outside the Lindblad resonances, there exists a possible spiral density-wave branch for which perturbation enhancements of surface mass density in stellar and gaseous discs are out of phase. We discuss implications of these results on the critical parameters for global star formation in barred and normal spiral galaxies and on magnetohydrodynamic density waves within the Lindblad resonances.     

Instabilities of galactic discs in the presence of star formation

We discuss the stability of galactic discs in which the energy of interstellar clouds is gained in encounters with expanding supernova (SN) remnants and lost in inelastic collisions. Energy gain and loss processes introduce a phase difference between the pressure and density perturbations, making discs unstable on small scales for several recipes of star formation. This is in contrast to the standard stability analysis in which small-scale perturbations are stabilized by pressure. In the limit of small scales, the dispersion relation for the growth rate reduces to that of thermal instabilities in a fluid without gravity. If instabilities lead to star formation, then our results imply a secondary mode of star formation that operates on small scales and feeds on the existence of a primary mode on intermediate scales. This may be interpreted as positive feedback. Further, the standard stability criterion on intermediate scales is significantly modified.     

The excitation, propagation and dissipation of waves in accretion discs: the non-linear axisymmetric case

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing, and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic  ( γ =5/3)  . Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance  ℳ_r≳0.01  undergo shocks within a distance of order the resonance radius.     

Luminous hot accretion discs

We find a new two-temperature hot branch of equilibrium solutions for stationary accretion discs around black holes. In units of Eddington accretion rate defined as 10 L _Edd c ², the accretion rates to which these solutions correspond are within the range ̇ ₁≲ ̇ ≲1, where ̇ ₁ is the critical rate of advection-dominated accretion flow (ADAF). In these solutions, the energy loss rate of the ions by Coulomb energy transfer between the ions and electrons is larger than the viscously heating rate and it is the advective heating together with the viscous dissipation that balances the Coulomb cooling of ions. When ̇ ₁≲ ̇ ≲ ̇ ₂, where ̇ ₂∼5 ̇ ₁<1, the accretion flow remains hot throughout the disc. When ̇ ₂≲ ̇ ≲1, Coulomb interaction will cool the inner region of the disc within a certain radius ( r _tr∼several tens of Schwarzschild radii or larger depending on the accretion rate and the outer boundary condition) and the disc will collapse on to the equatorial plane and form an optically thick cold annulus. Compared with ADAF, these hot solutions are much more luminous because of the high accretion rate and efficiency; therefore, we call them luminous hot accretion discs.     

The paradox of the scale-free discs

Scale-free discs have no preferred length or time-scale. The question has been raised whether such discs have a continuum of unstable linear modes or perhaps no unstable modes at all. We resolve this paradox by analysing the particular case of a gaseous, isentropic disc with a completely flat rotation curve (the Mestel disc) exactly . The heart of the matter is this: what are the correct boundary conditions to impose at the origin or central cusp? We argue that the linear stability problem is ill-posed and that similar ambiguities may afflict general disc models with power-law central cusps. From any finite radius, waves reach the origin after finite time but with logarithmically divergent phase. Instabilities exist, but their pattern speeds depend upon an undetermined phase with which waves are reflected from the origin. For any definite choice of this phase, there is an infinite but discrete set of growing modes. The ratio of growth rate to pattern speed is independent of the central phase. This ratio is derived in closed form for non-self-gravitating normal modes and is shown to agree with approximate results obtained from the shearing sheet in the short-wavelength limit. This provides the first exact, analytically solved stability analysis for a differentially rotating disc. For self-gravitating normal modes, the ratio of growth rate to pattern is found numerically by solving recurrence relations in Mellin-transform space.     

Magnetic fields and the dynamics of spiral galaxies

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D magnetohydrodynamics simulations of galactic discs subject to a spiral potential using cold gas, warm gas and a two-phase mixture of both. Recent hydrodynamic simulations have demonstrated the formation of interarm spurs as well as spiral arm molecular clouds, provided the interstellar medium model includes a cold H  i phase. We find that the main effect of adding a magnetic field to these calculations is to inhibit the formation of structure in the disc. However, provided a cold phase is included, spurs and spiral arm clumps are still present if β≳ 0.1 in the cold gas. A caveat to the two-phase calculations though is that by assuming a uniform initial distribution, β≳ 10 in the warm gas, emphasizing that models with more consistent initial conditions and thermodynamics are required. Our simulations with only warm gas do not show such structure, irrespective of the magnetic field strength.
Furthermore, we find that the introduction of a cold H  i phase naturally produces the observed degree of disorder in the magnetic field, which is again absent from simulations using only warm gas. Whilst the global magnetic field follows the large-scale gas flow, the magnetic field also contains a substantial random component that is produced by the velocity dispersion induced in the cold gas during the passage through a spiral shock. Without any cold gas, the magnetic field in the warm phase remains relatively well ordered apart from becoming compressed in the spiral shocks. Our results provide a natural explanation for the observed high proportions of disordered magnetic field in spiral galaxies and we thus predict that the relative strengths of the random and ordered components of the magnetic field observed in spiral galaxies will depend on the dynamics of spiral shocks.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Dust-penetrated arm classes: insights from rising and falling rotation curves

In the last decade, near-infrared imaging has highlighted the decoupling of gaseous and old stellar discs: the morphologies of optical (Population I) tracers compared to the old stellar disc morphology, can be radically different. Galaxies which appear multi-armed and even flocculent in the optical may show significant grand-design spirals in the near-infrared. Furthermore, the optically determined Hubble classification scheme does not provide a sound way of classifying dust-penetrated stellar discs: spiral arm pitch angles (when measured in the near-infrared) do not correlate with Hubble type. The dust-penetrated classification scheme of Block & Puerari provides an alternative classification based on near-infrared morphology, which is thus more closely linked to the dominant stellar mass component. Here we present near-infrared K -band images of 14 galaxies, on which we have performed a Fourier analysis of the spiral structure in order to determine their near-infrared pitch angles and dust-penetrated arm classes. We have also used the rotation curve data of Mathewson et al. to calculate the rates of shear in the stellar discs of these galaxies. We find a correlation between near-infrared pitch angle and rate of shear: galaxies with wide open arms (the γ class) are found to have rising rotation curves, while those with falling rotation curves belong to the tightly wound α bin. The major determinant of near-infrared spiral arm pitch angle is the distribution of matter within the galaxy concerned. The correlation reported in this study provides the physical basis underpinning spiral arm classes in the dust-penetrated regime and underscores earlier spectroscopic findings by Burstein and Rubin that Hubble type and mass distributions are unrelated.