Shatter or not: role of temperature and metallicity in the evolution of thermal instability
ABSTRACT We test how metallicity variation (a background gradient and fluctuations) affects the physics of local thermal instability using analytical calculations and idealized, high-resolution 1D hydrodynamic simulations. Although the cooling function (Λ[T, Z]) and the cooling time (tcool) depend on gas temperature and metallicity, we find that the growth rate of thermal instability is explicitly dependent only on the derivative of the cooling function relative to temperature (∂ln Λ/∂ln T) and not on the metallicity derivative (∂ln Λ/∂ln Z). For most of 104 K ≲ T ≲ 107 K, both the isobaric and isochoric modes (occurring at scales smaller and larger than the sonic length covered in a cooling time [cstcool], respectively) grow linearly, and at higher temperatures (≳107 K) the isochoric modes are stable. We show that even the non-linear evolution depends on whether the isochoric modes are linearly stable or unstable. For the stable isochoric modes, we observe the growth of small-scale isobaric modes but this is distinct from the non-linear fragmentation of a dense cooling region. For unstable isochoric perturbations we do not observe large density perturbations at small scales. While very small clouds (∼min[cstcool]) form in the transient state of non-linear evolution of the stable isochoric thermal instability, most of them merge eventually.