Context. Supersonic disordered flows accompany the formation and evolution of MCs. It has been argued that turbulence can support against gravitational collapse and form hierarchical sub-structures. Aims. We study the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Do the simulated flows agree with observations? Methods. We analyse three MCs that have formed self-consistently within kpc-scale numerical simulations of the ISM. The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using VSFs with/out density weighting or a density cutoff, and computed in one or three dimensions. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence consistent with Larson's size-velocity relationship. SN blast wave impacts produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Conclusions. We conclude that two different effects coincidentally reproduce Larson's size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs consistent with Larson's relationship. Later, contraction dominates, the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson's relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.

Theories suggest that filament fragmentation should occur on a characteristic fragmentation length-scale. This fragmentation length-scale can be related to filament properties, such as the width and the dynamical state of the filament. Here we present a study of a number of fragmentation analysis techniques applied to filaments, and their sensitivity to characteristic fragmentation length-scales. We test the sensitivity to both single-tier and two-tier fragmentation, i.e. when the fragmentation can be characterised with one or two fragmentation length-scales respectively. The nearest neighbour separation, minimum spanning tree separation and two-point correlation function are all able to robustly detect characteristic fragmentation length-scales. The Fourier power spectrum and the Nth nearest neighbour technique are both poor techniques, and require very little scatter in the core spacings for the characteristic length-scale to be successfully determined. We develop a null hypothesis test to compare the results of the nearest neighbour and minimum spanning tree separation distribution with randomly placed cores. We show that a larger number of cores is necessary to successfully reject the null hypothesis if the underlying fragmentation is two-tier, N>20. Once the null is rejected we show how one may decide if the observed fragmentation is best described by single-tier or two-tier fragmentation, using either Akaike's information criterion or the Bayes factor. The analysis techniques, null hypothesis tests, and model selection approaches are all included in a new open-source Python/C library called FragMent.

We investigate the equilibrium charge distribution of dust grains in the interstellar medium (ISM). Our treatment accounts for collisional charging by electrons and ions, photoelectric charging due to a background interstellar radiation field, the collection of suprathermal cosmic ray electrons and photoelectric emission due to a cosmic ray induced ultraviolet radiation field within dense molecular clouds. We find that the charge equilibrium assumption is valid throughout the multi-phase ISM conditions investigated here, and should remain valid for simulations with resolutions down to AU scales. The charge distribution of dust grains is size, composition, and ISM environment dependent: local radiation field strength, $G$, temperature, $T$, and electron number density, $n_{\mathrm{e}}$. The charge distribution is tightly correlated with the `charging parameter', $G\sqrt{T}/n_{\mathrm{e}}$. In the molecular medium, both carbonaceous and silicate grains have predominantly negative or neutral charges with narrow distributions. In the cold neutral medium, carbonaceous and silicate grains vary from negative and narrow distributions, to predominantly positive and wide distributions depending on the magnitude of the charging parameter. In the warm neutral medium, grains of all sizes are positively charged with wide distributions. We derive revised parametric expressions that can be used to recover the charge distribution function of carbonaceous and silicate grains from 3.5 {\AA} to 0.25 $\mu$m as a function of the size, composition and ambient ISM parameters. Finally, we find that the parametric equations can be used in environments other than Solar neighborhood conditions, recovering the charge distribution function of dust grains in photon dominated regions.

In order to understand the origin of observed molecular cloud properties, it is critical to understand how clouds interact with their environments during their formation, growth, and collapse. It has been suggested that accretion-driven turbulence can maintain clouds in a highly turbulent state, preventing runaway collapse, and explaining the observed non-thermal velocity dispersions. We present 3D, AMR, MHD simulations of a kiloparsec-scale, stratified, supernova-driven, self-gravitating, interstellar medium, including diffuse heating and radiative cooling. These simulations model the formation and evolution of a molecular cloud population in the turbulent interstellar medium. We use zoom-in techniques to focus on the dynamics of the mass accretion and its history for individual molecular clouds. We find that mass accretion onto molecular clouds proceeds as a combination of turbulent and near free-fall accretion of a gravitationally bound envelope. Nearby supernova explosions have a dual role, compressing the envelope, boosting accreted mass, but also disrupting parts of the envelope and eroding mass from the cloud's surface. It appears that the inflow rate of kinetic energy onto clouds from supernova explosions is insufficient to explain the net rate of charge of the cloud kinetic energy. In the absence of self-consistent star formation, conversion of gravitational potential into kinetic energy during contraction seems to be the main driver of non-thermal motions within clouds. We conclude that although clouds interact strongly with their environments, bound clouds are always in a state of gravitational contraction, close to runaway, and their properties are a natural result of this collapse.

We present synthetic dust polarisation maps of 3D magneto-hydrodynamical simulations of molecular clouds before the onset of stellar feedback. The clouds are modelled within the SILCC-Zoom project and are embedded in their galactic environment. The radiative transfer is carried out with POLARIS for wavelengths from 70 $\mu$m to 3 mm at a resolution of 0.12 pc, and includes self-consistently calculated alignment efficiencies for radiative torque alignment. We explore the reason of the observed depolarisation in the center of molecular clouds: We find that dust grains remain well aligned even at high densities ($n$ $>$ 10$^3$ cm$^{-3}$) and visual extinctions ($A_\text{V}$ $>$ 1). The depolarisation is rather caused by strong variations of the magnetic field direction along the LOS due to turbulent motions. The observed magnetic field structure thus resembles best the mass-weighted, line-of-sight averaged field structure. Furthermore, it differs by only a few 1$^\circ$ for different wavelengths and is little affected by the spatial resolution of the synthetic observations. Noise effects can be reduced by convolving the image. Doing so, for $\lambda$ $\gtrsim$ 160 $\mu$m the observed magnetic field traces reliably the underlying field in regions with intensities $I$ $\gtrsim$ 2 times the noise level and column densities above 1 M$_\text{sun}$ pc$^{-2}$. Here, typical deviations are $\lesssim$ 10$^\circ$. The observed structure is less reliable in regions with low polarisation degrees and possibly in regions with large column density gradients. Finally, we show that a simplified and widely used method without self-consistent dust alignment efficiencies can provide a good representation of the observable polarisation structure with deviations below 5$^\circ$.

We investigate the impact of the far ultraviolet (FUV) heating rate on the stability of the three-phase interstellar medium using three-dimensional simulations of a $1$ kpc$^2$, vertically-extended domain. The FUV heating rate sets the range of thermal pressures across which the cold ($\sim10^2$ K) and warm ($\sim10^4$ K) neutral media (CNM and WNM) can coexist in equilibrium. Even absent a variable star formation rate regulating the FUV heating rate, the gas physics keeps the pressure in the two-phase regime: because radiative heating and cooling processes happen on shorter timescales than sound wave propagation, turbulent compressions tend to keep the interstellar medium within the CNM-WNM pressure regime over a wide range of heating rates. The thermal pressure is set primarily by the heating rate with little influence from the hydrostatics. The vertical velocity dispersion adjusts as needed to provide hydrostatic support given the thermal pressure: when the turbulent pressure $\langle\rho\rangle\sigma_z^2$ is calculated over scales $\gtrsim500$ pc, the thermal plus turbulent pressure approximately equals the weight of the gas. The warm gas volume filling fraction is $0.2<f_w<0.8$ over a factor of less than three in heating rate, with $f_w$ near unity at higher heating rates and near zero at lower heating rates. We suggest that cosmological simulations that do not resolve the CNM should maintain an interstellar thermal pressure within the two-phase regime.

A large fraction of the gas in the Galaxy is cold, dense, and molecular. If all this gas collapsed under the influence of gravity and formed stars in a local free-fall time, the star formation rate in the Galaxy would exceed that observed by more than an order of magnitude. Other star-forming galaxies behave similarly. Yet observations and simulations both suggest that the molecular gas is indeed gravitationally collapsing, albeit hierarchically. Prompt stellar feedback offers a potential solution to the low observed star formation rate if it quickly disrupts star-forming clouds during gravitational collapse. However, this requires that molecular clouds must be short-lived objects, raising the question of how so much gas can be observed in the molecular phase. This can occur only if molecular clouds form as quickly as they are destroyed, maintaining a global equilibrium fraction of dense gas. We therefore examine cloud formation timescales. We first demonstrate that supernova and superbubble sweeping cannot produce dense gas at the rate required to match the cloud destruction rate. On the other hand, Toomre gravitational instability can reach the required production rate. We thus argue that, although dense, star-forming gas may last only around a single global free-fall time, the dense gas in star-forming galaxies can globally exist in a state of dynamic equilibrium between formation by gravitational instability, and disruption by stellar feedback. At redshift z >~ 2, the Toomre instability timescale decreases, resulting in a prediction of higher molecular gas fractions at early times, in agreement with observations.

Molecular cloud observations show that clouds have non-thermal velocity dispersions that scale with the cloud size as $\sigma\propto R^{1/2}$ at constant surface density, and for varying surface density scale with both the cloud`s size and surface density, $\sigma^2 \propto R \Sigma$. The energy source driving these chaotic motions remains poorly understood. We describe the velocity dispersions observed in a cloud population formed in a kiloparsec-scale numerical simulation of a magnetized, supernova-driven, self-gravitating, interstellar medium, including diffuse heating and radiative cooling. We compare the relationships between velocity dispersion, size, and surface density measured in the simulated cloud population to those found in observations of Galactic molecular clouds. We find that external supernova explosions can not drive turbulent motions of the observed magnitudes within dense clouds. On the other hand, self-gravity also induces non-thermal motions as gravitationally bound clouds begin to collapse in our model, and by doing so their internal velocity dispersions recover the observed relations. Energy conservation suggests that the observed behavior is consistent with the kinetic energy being proportional to the gravitational energy. However, the clouds in our model show no sign of reaching a stable equilibrium state at any time, even for strongly magnetized clouds. We conclude that gravitationally bound molecular clouds are always in a state of gravitational collapse and their properties are a natural result of this chaotic collapse. In order to agree with observed star formation efficiencies, this process must be terminated by the early destruction of the clouds, presumably from internal stellar feedback.

We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_\alpha$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.