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.