We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (sphNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the sphNG or GIZMO runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the sphNG and GIZMO runs and secondary spiral arms are more pronounced. By resolving the Jeans' length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and sphNG/GIZMO. Although more similar, sphNG displays different density distributions and vertical mass profiles to all modes of GIZMO (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations.

We present a method to simultaneously model the dust far-infrared spectral energy distribution (SED) and the total infrared $-$ carbon monoxide (CO) integrated intensity $(S_{\rm IR}-I_{\rm CO})$ relationship. The modelling employs a hierarchical Bayesian (HB) technique to estimate the dust surface density, temperature ($T_{\rm eff}$), and spectral index at each pixel from the observed far-infrared (FIR) maps. Additionally, given the corresponding CO map, the method simultaneously estimates the slope and intercept between the FIR and CO intensities, which are global properties of the observed source. The model accounts for correlated and uncorrelated uncertainties, such as those present in Herschel observations. Using synthetic datasets, we demonstrate the accuracy of the HB method, and contrast the results with common non-hierarchical fitting methods. As an initial application, we model the dust and gas on 100 pc scales in the Magellanic Clouds from Herschel FIR and NANTEN CO observations. The slopes of the $\log S_{\rm IR}-\log I_{\rm CO}$ relationship are similar in both galaxies, falling in the range 1.1$-$1.7. However, in the SMC the intercept is nearly 3 times higher, which can be explained by its lower metallicity than the LMC, resulting in a larger $S_{\rm IR}$ per unit $I_{\rm CO}$. The HB modelling evidences an increase in $T_{\rm eff}$ in regions with the highest $I_{\rm CO}$ in the LMC. This may be due to enhanced dust heating in the densest molecular regions from young stars. Such simultaneous dust and gas modelling may reveal variations in the properties of the ISM and its association with other galactic characteristics, such as star formation rates and/or metallicities.

We model driven, compressible, isothermal, turbulence with Mach numbers ranging from the subsonic ($\mathcal{M} \approx 0.65$) to the highly supersonic regime ($\mathcal{M}\approx 16 $). The forcing scheme consists both solenoidal (transverse) and compressive (longitudinal) modes in equal parts. We find a relation $\sigma_{s}^2 = \mathrm{b}\log{(1+\mathrm{b}^2\mathcal{M}^2)}$ between the Mach number and the standard deviation of the logarithmic density with $\mathrm{b} = 0.457 \pm 0.007$. The density spectra follow $\mathcal{D}(k,\,\mathcal{M}) \propto k^{\zeta(\mathcal{M})}$ with scaling exponents depending on the Mach number. We find $\zeta(\mathcal{M}) = \alpha \mathcal{M}^{\beta}$ with a coefficient $\alpha$ that varies slightly with resolution, whereas $\beta$ changes systematically. We extrapolate to the limit of infinite resolution and find $\alpha = -1.91 \pm 0.01,\, \beta =-0.30\pm 0.03$. The dependence of the scaling exponent on the Mach number implies a fractal dimension $D=2+0.96 \mathcal{M}^{-0.30}$. We determine how the scaling parameters depend on the wavenumber and find that the density spectra are slightly curved. This curvature gets more pronounced with increasing Mach number. We propose a physically motivated fitting formula $\mathcal{D}(k) = \mathcal{D}_0 k^{\zeta k^{\eta}}$ by using simple scaling arguments. The fit reproduces the spectral behaviour down to scales $k\approx 80$. The density spectrum follows a single power-law $\eta = -0.005 \pm 0.01$ in the low Mach number regime and the strongest curvature $\eta = -0.04 \pm 0.02$ for the highest Mach number. These values of $\eta$ represent a lower limit, as the curvature increases with resolution.

Turbulence is a dominant feature operating in gaseous flows across nearly all scales in astrophysical environments. Accordingly, accurately estimating the statistical properties of such flows is necessary for developing a comprehensive understanding of turbulence. We develop and employ a hierarchical Bayesian fitting method to estimate the parameters describing the scaling relationships of the velocity power spectra of supersonic turbulence. We demonstrate the accuracy and other advantages of this technique compared with ordinary linear regression methods. Using synthetic power spectra, we show that the Bayesian method provides accurate parameter and error estimates. Commonly used normal linear regression methods can provide estimates that fail to recover the underlying slopes, up to 70% of the instances, even when considering the 2std.dev. uncertainties. Additionally, we apply the Bayesian methods to analyse the statistical properties of compressible turbulence in 3D numerical simulations. We model driven, isothermal, turbulence with rms Mach numbers in the highly supersonic regime M~15. We study the influence of purely solenoidal (divergence-free) and purely compressive (curl-free) forcing on the scaling exponent of the power spectrum. In simulations with solenoidal forcing and 1024^3 resolution, our results indicate that there is no extended inertial range with a constant scaling exponent. The bottleneck effect results in a curved power spectrum at all wave numbers and is more pronounced in the transversal modes compared with the longitudinal modes. Therefore, this effect is stronger in stationary turbulent flows driven by solenoidal forcing compared to the compressive one. The longitudinal spectrum driven with compressive forcing is the only spectrum with constant scaling exponent z=-1.94 +- 0.01, corresponding to slightly shallower slopes than the Burger's prediction.

We compute structure functions and Fourier spectra of 2D centroid velocity (CV) maps in order to study the gas dynamics of typical molecular clouds (MCs) in numerical simulations. We account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of the gas and use a 3D radiative transfer tool to model the CO line emission in a post-processing step. We perform simulations using three different initial mean number densities of n_0 = 30, 100 and 300 cm^{-3} to span a range of typical values for dense gas clouds in the solar neighbourhood. We compute slopes of the centroid velocity increment structure functions (CVISF) and of Fourier spectra for different chemical components: the total density, H2 number density, 12CO number density as well as the integrated intensity of 12CO (J=1-0) and 13CO (J=1-0). We show that optical depth effects can significantly affect the slopes derived for the CVISF, which also leads to different scaling properties for the Fourier spectra. The slopes of CVISF and Fourier spectra for H2 are significantly steeper than those for the different CO tracers, independent of the density and the numerical resolution. This is due to the larger space-filling factor of H2 as it is better able to self-shield in diffuse regions, leading to a larger fractal co-dimension compared to CO.

The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple appoximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form $\mathrm{P}_v(\rho)\propto\rho^{-1.54}$ for the (volume-weighted) PDF and $\mathrm{P}_m(\rho)\propto\rho^{-0.54}$ for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates. Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured star formation rates are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall.

We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population. Assuming standard conversion factors to estimate Sigma_SFR and Sigma_mol from the observations, we find that the KS parameters vary between galaxies, indicating that no universal relationship holds for all galaxies. The KS slope of the whole population is 0.76, with the 2sigma range extending from 0.58 to 0.94. These results imply that the molecular gas depletion time is not constant, but varies from galaxy to galaxy, and increases with the molecular gas surface density. Therefore, other galactic properties besides just Sigma_mol affect Sigma_SFR, such as the gas fraction or stellar mass. The non-universality of the KS relationship indicates that a comprehensive theory of star formation must take into account additional physical processes that may vary from galaxy to galaxy.

The probability density function (PDF) of the gas density in subsonic and supersonic, isothermal, driven turbulence is analyzed with a systematic set of hydrodynamical grid simulations with resolutions up to 1024^3 cells. We performed a series of numerical experiments with root mean square (r.m.s.) Mach number M ranging from the nearly incompressible, subsonic (M=0.1) to the highly compressible, supersonic (M=15) regime. We study the influence of two extreme cases for the driving mechanism by applying a purely solenoidal (divergence-free) and a purely compressive (curl-free) forcing field to drive the turbulence. We find that our measurements fit the linear relation between the r.m.s. Mach number and the standard deviation of the density distribution in a wide range of Mach numbers, where the proportionality constant depends on the type of the forcing. In addition, we propose a new linear relation between the standard deviation of the density distribution and the standard deviation of the velocity in compressible modes, i.e. the compressible component of the r.m.s. Mach number. In this relation the influence of the forcing is significantly reduced, suggesting a linear relation between the standard deviation of the density distribution and the standard deviation of the velocity in compressible modes, independent of the forcing, ranging from the subsonic to the supersonic regime.

We present a systematic study of the influence of different forcing types on the statistical properties of supersonic, isothermal turbulence in both the Lagrangian and Eulerian frameworks. We analyse a series of high-resolution, hydrodynamical grid simulations with Lagrangian tracer particles and examine the effects of solenoidal (divergence-free) and compressive (curl-free) forcing on structure functions, their scaling exponents, and the probability density functions of the gas density and velocity increments. Compressively driven simulations show a significantly larger density contrast, a more intermittent behaviour, and larger fractal dimension of the most dissipative structures at the same root mean square Mach number. We show that the absolute values of Lagrangian and Eulerian structure functions of all orders in the integral range are only a function of the root mean square Mach number, but independent of the forcing. With the assumption of a Gaussian distribution for the probability density function of the velocity increments on large scales, we derive a model that describes this behaviour.