Interfacial instabilities driven by self-gravity: first numerical results

The evolution of the interstellar medium (ISM) is driven by a variety of phenomena, including turbulence, shearing flows, magnetic fields and the thermal properties of the gas. Among the most important forces at work is self-gravity, which ultimately drives protostellar collapse. As part of an ongoing study of instabilities in the ISM, Hunter, Whitaker & Lovelace have discovered another process driven by self-gravity: the instability of an interface of discontinuous density. Theory predicts that this self-gravity driven interfacial instability persists in the static limit and in the absence of a constant background acceleration. Disturbances to a density interface are found to grow on a time-scale of the order of the free-fall time, even when the perturbation wavelength is much less than the Jeans length. Here we present the first numerical simulations of this instability. The theoretical growth rate is confirmed and the non-linear morphology displayed. The self-gravity interfacial instability is shown to be fundamentally different from the Rayleigh–Taylor instability, although both exhibit similar morphologies under the condition of a high density contrast, such as is commonly found in the ISM. Such instabilities are a possible mechanism by which observed features, such as the pillars of gas seen near the boundaries of interstellar clouds, are formed.     

Parker instability in a self-gravitating magnetized gas disk: II. Competition between gravitational and convective instabilities

Under influence of external gravity generated by Galactic all components excluding ISM, a magnetized gas disk may experience both Parker and convective instabilities. Growth rate of the convective instability increases with decreasing perturbation wavelength, and the convective motion makes sheet-like structures all over before the Parker instability forms structures of any meaningful size in the disk. Yet the Parker instability is thought to be an ideal route to form large-scale condensations in the Galaxy. In search of a means to curb convective activities in the Galactic ISM disk, the external gravity is replaced by self-gravity as a driving force of the Parker instability and the gravitational instability is invoked to reinforce the Parker instability. Perturbation of interchange mode is known to trigger convective instability in such disk and the one of undular mode to activate the Parker instability, while the gravitational instability can be triggered by both modes. Therefore, the resulting Jeans instability would help the Parker instability to overcome disrupting behavior of the convection. Dynamical properties of the disk can be characterized by ratio α of magnetic to gas pressure, adiabatic exponent γ, scale height H of the ISM, and disk thickness z_a. A linear stability analysis has been done to the disk, and the maximum growth rate of the Parker–Jeans instability is compared with that of the convective instability. The latter may or may not be higher than the former, depending on the disk parameters. The Parker–Jeans instability has chances to override the convective instability, when the disk is thicker than a certain value. In the disk thinner than the critical one, the Jeans instability can always suppress the convection. Since the growth rate of the convective instability is proportional to local gravitational acceleration, thereby in the general Galactic gravity, the convective instability works actively only in upper regions, we expect chaotic features to appear in regions of low density far from Galactic mid-plane.     

Generation of magnetosonic waves and formation of structures in galaxies

We study the generation of magnetosonic waves in galactic gaseous discs taking account of the magnetic field, differential rotation and self-gravity. The special case of perturbations is considered with the wavevector perpendicular to the magnetic field. The necessary condition of the amplification of seed perturbations is the presence of differential rotation and non-vanishing radial component of the magnetic field that can easily be satisfied in galactic discs. Differential rotation stretches the azimuthal field from the radial one and, therefore, we consider the generation of waves on the time-dependent background magnetic field. Basically, an amplification is rather efficient, and seed perturbations become non-linear already after several rotation periods for a wide range of wavelength. The generated magnetosonic waves can be either non-oscillatory or oscillatory depending on the parameters of gas. If perturbations are Jeans stable, then typically non-oscillatory waves are amplified. However, interplay between self-gravity, magnetic field and rotational shear can change qualitatively the classical Jeans instability, so that the latter becomes oscillatory and tends to be suppressed in galaxies.     

Scalar perturbations in two-temperature cosmological plasmas

We study the properties of density perturbations of a two-component plasma with a temperature difference on a homogeneous and isotropic background. For this purpose, we extend the general relativistic gauge-invariant and covariant (GIC) perturbation theory to include a multifluid with a particular equation of state (ideal gas) and imperfect fluid terms due to the relative energy flux between the two species. We derive closed sets of GIC vector and subsequently scalar evolution equations. We then investigate solutions in different regimes of interest. In particular, we study long-wavelength and arbitrary-wavelength Langmuir and ion-acoustic perturbations. The harmonic oscillations are superposed on a Jeans-type instability. We find a generalized Jeans criterion for collapse in a two-temperature plasma, which states that the species with the largest sound velocity determines the Jeans wavelength. Furthermore, we find that within the limit for gravitational collapse, initial perturbations in either the total density or charge density lead to a growth in the initial temperature difference. These results are relevant for the basic understanding of the evolution of inhomogeneities in cosmological models.     

Thermal instability in a weakly ionized plasma

We revisit the problem of clump formation due to thermal instabilities in a weakly ionized plasma with the help of a linear perturbation analysis, as discussed by Nejad-Asghar & Ghanbari. In the absence of a magnetic field and ambipolar diffusion the characteristic equation reduces to the thermal instability described by Field. We derive the critical wavelengths, which separate the spatial ranges of stability and instability. Contrary to the original analysis of Nejad-Asghar & Ghanbari, perturbations with a wavelength larger than the critical wavelength destabilize the cloud. Moreover, the instability regime of isentropic perturbations is drastically reduced. Isobaric modes with real values of the critical wavelength appear only if the density dependence of the cooling rate is more pronounced than the temperature dependence. Isentropic modes arise only if the power of the density in the cooling rate is smaller than 1/2, which is not fulfilled for CO cooling. We find that ambipolar diffusion is not a dominating heating process in molecular gas.     

The Formation of Knots and Filaments in Shocks

We present a survey of different kinds of instabilities in the context of radiative colliding flows which greatly contribute to structure formation. In particular, this includes analytical results for different kinds of thin shell instabilities (DI, NDI, NTSI). New numerical results for the non-linear evolution of such instabilities in two dimensions, and their coupling with the thermal cooling instability are presented. The astrophysical implications are briefly outlined, in particular the formation of knots and filaments. This revised version was published online in July 2006 with corrections to the Cover Date.     

Self-gravitational instability of rotating viscous Hall plasma with arbitrary radiative heat-loss functions and electron inertia

The effects of arbitrary radiative heat-loss functions and Hall current on the self-gravitational instability of a homogeneous, viscous, rotating plasma has been investigated incorporating the effects of finite electrical resistivity, finite electron inertia and thermal conductivity. A general dispersion relation is obtained using the normal mode analysis with the help of relevant linearized perturbation equations of the problem, and a modified Jeans criterion of instability is obtained. The conditions of modified Jeans instabilities and stabilities are discussed in the different cases of our interest. We find that the presence of arbitrary radiative heat-loss functions and thermal conductivity modifies the fundamental Jeans criterion of gravitational instability into a radiative instability criterion. The Hall parameter affects only the longitudinal mode of propagation and it has no effect on the transverse mode of propagation. For longitudinal propagation, it is found that the condition of radiative instability is independent of the magnetic field, Hall parameter, finite electron inertia, finite electrical resistivity, viscosity and rotation; but for the transverse mode of propagation it depends on the finite electrical resistivity, the strength of the magnetic field, and it is independent of rotation, electron inertia and viscosity. From the curves we find that the presence of thermal conductivity, finite electrical resistivity and density-dependent heat-loss function has a destabilizing influence, while viscosity and magnetic field have a stabilizing effect on the growth rate of an instability. The effect of arbitrary heat-loss functions is also studied on the growth rate of a radiative instability.     

Saturated magnetic field amplification at supernova shocks

Cosmic ray streaming instabilities at supernova shocks are discussed in the quasi-linear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the non-resonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the non-resonant instability is actually driven by a compensating current in the background plasma. The non-resonant instability can develop into a non-linear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic ray pressures is predominantly determined by non-resonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.     

Evolution of the cosmological density distribution function

We present a new calculation for the evolution of the one-point probability distribution function (PDF) of the cosmological density field based on an exact statistical treatment. Using the Chapman–Kolmogorov equation and second-order Eulerian perturbation theory we propagate the initial density distribution into the non-linear regime. Our calculations yield the moment generating function, allowing a straightforward derivation of the skewness of the PDF to second order. We find a new dependence on the initial perturbation spectrum. We compare our results with other approximations to the one-point PDF, and with N -body simulations. We find that our distribution accurately models the evolution of the one-point PDF of dark matter.     

Numerical simulations of rotating cooling flows in galaxy cluster environments

We have used 2D numerical simulations to study the evolution of galaxy cluster cooling flows undergoing a rotational perturbation. We show that such rotations in the intracluster medium may arise from cluster/subcluster mergers. Our galaxy cluster initial conditions involve spherically symmetric, steady-state cooling flows with varying mass-dropout strengths. The rotational perturbation serves to break the symmetry for each of the initial cooling flows, resulting in the formation of thin, gaseous disc-like structure extending radially out to ∼10 kpc. Disc-like structure formed for low mass-dropout strength simulations appears to contain cooling condensations whereas disc-like structure in higher mass-dropout strength simulations appears smooth. This is due to the influence of mass-dropout on the degree of cooling, which serves to reduce the strength of thermal instabilities by the removal of 'cold' gas from the flow. Morphological comparisons of the disc-like structure formed in our simulations are made to structure observed in the X-ray emitting gas of A4059. Comparisons of the gas dynamics within the disc-like structure are also made to the solid-body rotation profile observed from emission-line gas within the central galaxy of Hydra A. The influence of grid effects on the simulations is also discussed.     

Global lopsided instability in a purely stellar galactic disc

It is shown that pure exponential discs in spiral galaxies are capable of supporting slowly varying discrete global lopsided modes, which can explain the observed features of lopsidedness in the stellar discs. Using linearized fluid dynamical equations with the softened self-gravity and pressure of the perturbation as the collective effect, we derive self-consistently a quadratic eigenvalue equation for the lopsided perturbation in the galactic disc. On solving this, we find that the ground-state mode shows the observed characteristics of the lopsidedness in a galactic disc, namely the fractional Fourier amplitude A ₁, increases smoothly with the radius. These lopsided patterns precess in the disc with a very slow pattern speed with no preferred sense of precession. We show that the lopsided modes in the stellar disc are long-lived because of a substantial reduction (approximately a factor of 10 compared to the local free precession rate) in the differential precession. The numerical solution of the equations shows that the ground-state lopsided modes are either very slowly precessing stationary normal mode oscillations of the disc or growing modes with a slow growth rate depending on the relative importance of the collective effect of the self-gravity. N -body simulations are performed to test the spontaneous growth of lopsidedness in a pure stellar disc. Both approaches are then compared and interpreted in terms of long-lived global   m = 1  instabilities, with almost zero pattern speed.     

Density wave formation in differentially rotating disk galaxies: Hydrodynamic simulation of the linear regime

Most rapidly and differentially rotating disk galaxies, in which the sound speed (thermal velocity dispersion) is smaller than the orbital velocity, display graceful spiral patterns. Yet, over almost 240 yr after their discovery in M51 by Charles Messier, we still do not fully understand how they originate. In this first paper of a series, the dynamical behavior of a rotating galactic disk is examined numerically by a high-order Godunov hydrodynamic code. The code is implemented to simulate a two-dimensional flow driven by an internal Jeans gravitational instability in a nonresonant wave–“fluid” interaction in an infinitesimally thin disk composed of stars or gas clouds. A goal of this work is to explore the local and linear regimes of density wave formation, employed by Lin, Shu, Yuan and many others in connection with the problem of spiral pattern of rotationally supported galaxies, by means of computer-generated models and to compare those numerical results with the generalized fluid-dynamical wave theory. The focus is on a statistical analysis of time-evolution of density wave structures seen in the simulations. The leading role of collective processes in the formation of both the circular and spiral density waves (“heavy sound”) is emphasized. The main new result is that the disk evolution in the initial, quasilinear stage of the instability in our global simulations is fairly well described using the local approximation of the generalized wave theory. Certain applications of the simulation to actual gas-rich spiral galaxies are also explored.     

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.     

The formation of primordial stars triggered by the evolution of large-scale density perturbations

This paper summarizes the basic idea of coupling of density perturbation modes in Newtonian cosmology and the formation of primordial stars. The way of deriving the second-order differential equation governing the growth of two coupled adiabatic density perturbation modes in the matter era is briefly discussed. A Jeans criterion for the growth of short wavelength perturbations is given. A mechanism is proposed for the cosmological origin of Population III stars.     

A possible mechanism governing the production and evolution of Ellerman bombs

A hydrodynamic-magnetofluid hybrid analysis of lower chromospheric shear flows in the Sun may explain the occurrence and time development of Ellerman bombs. The analysis assumes that the erupting material forming the bomb is driven initially by the Kelvin-Helmholtz fluid instability applied to the interface between two atmospheric fluid layers, characterized by a steep density change across the boundary and driven by flow fields around sunspots. The ensuing instability eventually evolves into a magnetofluid phenomenon by virtue of the trapping and bending of the interfacial magnetic field, giving rise to a dense globule of material entering, and persisting in, the upper layers and due to Ohmic dissipation having a significantly enhanced temperature compared with ambient material.     

Global m= 1 instabilities and lopsidedness in disc galaxies

Lopsidedness is common in spiral galaxies. Often, there is no obvious external cause, such as an interaction with a nearby galaxy, for such features. Alternatively, the lopsidedness may have an internal cause, such as a dynamical instability. In order to explore this idea, we have developed a computer code that searches for self-consistent perturbations in razor-thin disc galaxies and performed a thorough mode-analysis of a suite of dynamical models for disc galaxies embedded in an inert dark matter halo with varying amounts of rotation and radial anisotropy.
Models with two equal-mass counter-rotating discs and fully rotating models both show growing lopsided modes. For the counter-rotating models, this is the well-known counter-rotating instability, becoming weaker as the net rotation increases. The m = 1 mode of the maximally rotating models, on the other hand, becomes stronger with increasing net rotation. This rotating m = 1 mode is reminiscent of the eccentricity instability in near-Keplerian discs.
To unravel the physical origin of these two different m = 1 instabilities, we studied the individual stellar orbits in the perturbed potential and found that the presence of the perturbation gives rise to a very rich orbital behaviour. In the linear regime, both instabilities are supported by aligned loop orbits. In the non-linear regime, other orbit families exist that can help support the modes. In terms of density waves, the counter-rotating m = 1 mode is due to a purely growing Jeans-type instability. The rotating m = 1 mode, on the other hand, grows as a result of the swing amplifier working inside the resonance cavity that extends from the disc centre out to the radius where non-rotating waves are stabilized by the model's outwardly rising Q profile.     

Nonlinear thin shell instabilities in molecular clouds

Observations of molecular clouds point to the existence of supersonic, turbulent flows. Therefore, any theory which attempts to describe molecular cloud evolution and star formation must include a consideration of the dynamics of colliding flows. Previous studies have considered the collision of supersonic streams or clouds. The resultant instabilities provide a mechanism which may give rise to observable cloud morphologies and enhance the star formation rate. One such instability is the nonlinear thin shell instability (NTSI) of a shock-bounded slab. This process is driven by ram pressure and efficient cooling. In this study, I use numerical simulations to examine the head-on collision of supersonic gas streams in a cold, molecular gas. A dense slab forms in the collision midplane and is prone to a number of instabilities, including the NTSI. The thermodynamic processes involved are found to have a controlling influence upon the instability and fragmentation of the slab. Although some minimal amount of cooling is needed to drive the instability, too rapid a cooling rate gives rise to smaller wavelength instabilities which wipe out the NTSI. The growth rate of the NTSI in a gas undergoing molecular cooling corresponds to a timescale of order 10¹² s, in general agreement with the theoretical value for an isothermal gas. The NTSI may provide a viable mechanism for the instigation of rapid star formation.     

Time-dependent models of the structure and stability of self-gravitating protoplanetary discs

Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disc is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (  α∼ 10⁻³  or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several astronomical units at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the     relation for very young systems may allow a test of the model.