Steady states of the Parker instability

We study the linear properties, non-linear saturation, and a steady, strongly non-linear state of the Parker instability in galaxies. We consider magnetic buoyancy and its consequences with and without cosmic rays. Cosmic rays are described using the fluid approximation with anisotropic, non-Fickian diffusion. To avoid unphysical constraints on the instability (such as boundary conditions often used to specify an unstable background state), non-ideal magnetohydrodynamic equations are solved for deviations from a background state representing an unstable magnetohydrostatic equilibrium. We consider isothermal gas and neglect rotation. The linear evolution of the instability is in broad agreement with earlier analytical and numerical models; but we show that most of the simplifying assumptions of the earlier work do not hold, such that they provide only a qualitative rather than quantitative picture. In its non-linear stage the instability has significantly altered the background state from its initial state. Vertical distributions of both magnetic field and cosmic rays are much wider, the gas layer is thinner, and the energy densities of both magnetic field and cosmic rays are much reduced. The spatial structure of the non-linear state differs from that of any linear modes. A transient gas outflow is driven by the weakly non-linear instability as it approaches saturation.