Interactions of a shock with a molecular cloud at various stages of its evolution due to thermal instability and gravity

Abstract

Using the adaptive mesh refinement code mg, we perform hydrodynamic simulations of the interaction of a shock with a molecular cloud evolving due to thermal instability (TI) and gravity. To explore the relative importance of these processes, three case studies are presented. The first follows the formation of a molecular cloud out of an initially quiescent atomic medium due to the effects of TI and gravity. The second case introduces a shock whilst the cloud is still in the warm atomic phase, and the third scenario introduces a shock once the molecular cloud has formed. The shocks accelerate the global collapse of the clouds with both experiencing local gravitational collapse prior to this. When the cloud is still atomic, the evolution is shock dominated and structures form due to dynamical instabilities within a radiatively cooled shell. While the transmitted shock can potentially trigger the TI, this is prevented as material is shocked multiple times on the order of a cloud-crushing time-scale. When the cloud is molecular, the post-shock flow is directed via the pre-existing structure through low-density regions in the inter-clump medium. The clumps are accelerated and deformed as the flow induces clump–clump collisions and mergers that collapse under gravity. For a limited period, both shocked cases show a mixture of Kolmogorov and Burgers turbulence-like velocity and logarithmic density power spectra, and strongly varying density spectra. The clouds presented in this work provide realistic conditions that will be used in future feedback studies.

Momentum and energy injection by a wind-blown bubble into an inhomogeneous interstellar medium

“…4 http://wonka.physics.ncsu.edu/pub/VH-1/ 10 −26 erg s −1 . The cooling curve, Λ(𝑇), is calculated assuming collisional ionization equilibrium and is constructed from 3 separate parts (for further details see Wareing, Pittard & Falle 2017a,b;Kupilas et al 2021;Pittard, Kupilas & Wareing 2022). We assume solar abundances with mass fractions 𝑋 H = 0.7381, 𝑋 He = 0.2485, and 𝑋 Z = 0.0134 (cf.…”

The Survival of Multiphase Dusty Clouds in Hot Winds

“…Recent absorption-line measurements by HST/COS detect massive reservoirs of cool 10 4 K gas in the circumgalactic medium (Werk et al 2013). The origins of such gas is difficult to ascertain as it may form via thermal instability (Field 1965;Sharma et al 2012;McCourt et al 2012;Voit et al 2015;Ji et al 2018;Butsky et al 2020;Falle et al 2020;Kupilas et al 2021) or be entrained in galactic outflows. While observational constraints of the atomic phase are emerging, the molecular component remains difficult to study.…”

How D-type H ii region expansion depends on numerical resolution

“…For the recombination cooling we use Λ rec = 6.1 × 10 −10 𝑘𝑇 0.11 erg cm 3 s −1 if 𝑇 100 K (Osterbrock 1989). The cooling curve, Λ(𝑇), is constructed from 3 separate parts (see also Wareing et al 2017a,b;Kupilas et al 2021). At low temperatures (𝑇 < 10 4 K) we use a fit to the data in Koyama & Inutsuka (2002), corrected by Vazquez-Semadeni et al (2007):…”