Escaping the maze: a statistical sub-grid model for cloud-scale density structures in the interstellar medium

Abstract

The interstellar medium (ISM) is a turbulent, highly structured multi-phase medium. State-of-the-art cosmological simulations of the formation of galactic discs usually lack the resolution to accurately resolve those multi-phase structures. However, small-scale density structures play an important role in the life cycle of the ISM, and determine the fraction of cold, dense gas, the amount of star formation and the amount of radiation and momentum leakage from cloud-embedded sources. Here, we derive a $statisticalmodel$ to calculate the unresolved small-scale ISM density structure from coarse-grained, volume-averaged quantities such as the $gasclumpingfactor$, $\mathcal{C}$, and mean density ${⟨\rho ⟩}_V$. Assuming that the large-scale ISM density is statistically isotropic, we derive a relation between the three-dimensional clumping factor, ${\mathcal{C}}_{\rho }$, and the clumping factor of the $4\pi$ column density distribution on the cloud surface, ${\mathcal{C}}_{\Sigma }$, and find ${\mathcal{C}}_{\Sigma }={\mathcal{C}}_{\rho }^{2/3}$. Applying our model to calculate the covering fraction, i.e., the $4\pi$ sky distribution of optically thick sight-lines around sources inside interstellar gas clouds, we demonstrate that small-scale density structures lead to significant differences at fixed physical ISM density. Our model predicts that gas clumping increases the covering fraction by up to 30 per cent at low ISM densities compared to a uniform medium. On the other hand, at larger ISM densities, gas clumping suppresses the covering fraction and leads to increased scatter such that covering fractions can span a range from 20 to 100 per cent at fixed ISM density. All data and example code is publicly available at GitHub.

16 pages with 8 figures, 14 pages main text with 7 figures, 1 page references, 1 page appendix with 1 figure, accepted by MNRAS on April 1st 2022

References in corpus (24)

• The NumPy array: a structure for efficient numerical computation
• Analytical theory for the initial mass function: CO clumps and prestellar cores
• The Auriga Project: the properties and formation mechanisms of disc galaxies across cosmic time
• Momentum Injection by Supernovae in the Interstellar Medium
• Physical properties of molecular clouds for the entire Milky Way disk
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• Molecular Hydrogen and Global Star Formation Relations in Galaxies
• Numerical simulations of compressively driven interstellar turbulence: I. Isothermal gas
• Environmental Dependence of the Kennicutt-Schmidt Relation in Galaxies
• On the interplay between star formation and feedback in galaxy formation simulations
• The Physics of Star Cluster Formation and Evolution
• A Detailed Study of Feedback from a Massive Star
• Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies
• SILCC VI -- Multi-phase ISM structure, stellar clustering, and outflows with supernovae, stellar winds, ionising radiation and cosmic rays
• Physical Processes in Star Formation
• The shapes of column density PDFs - The importance of the last closed contour
• Density probability distribution functions of diffuse gas in the Milky Way
• The natural emergence of the correlation between H2 and star formation rate surface densities in galaxy simulations
• NIHAO -- XXV. Convergence in the cusp-core transformation of cold dark matter haloes at high star formation thresholds
• The density structure of supersonic self-gravitating turbulence
• First extragalactic measurement of the turbulence driving parameter: ALMA observations of the star-forming region N159E in the Large Magellanic Cloud
• Probability distribution functions of gas in M31 and M51
• Star formation in evolving molecular clouds
• Reconstructing three-dimensional densities from two-dimensional observations of molecular gas