On the influence of cooling and heating processes on the Parker instability – II. Numerical simulations

We investigate the Parker instability (PI) influenced by thermal processes in a non-adiabatic, gravitationally stratified interstellar medium and discuss a model including the photoionization heating together with the supplemental heating mechanisms postulated by Reynolds, Haffner and Tufte. A cooling rate due to radiative losses is described by an approximation to the realistic cooling function of Dalgarno and McCray for ionized interstellar gas. An unperturbed initial state of the system simultaneously represents both a magnetohydrostatic and thermal equilibrium, and is thermally stable. We perform a set of 3D numerical magnetohydrodynamic simulations using the zeusmp code. We find that PI developing in the presence of non-adiabatic effects promotes a transition of gas in magnetic valleys to a thermally unstable regime. We find that the region of initially enhanced density due to PI starts to condense more as the result of thermal instability action. The density in this region rises above the classical isothermal limit of two times the equilibrium value at the mid-plane. The maximum density in an evolved system reaches 10–40 times the equilibrium value at the mid-plane, and the structures so formed attain oval shapes. These results lead to the conclusion that PI, operating in the presence of realistic cooling and heating processes, can trigger the formation of dense clouds, which may give rise to giant molecular complexes.