Supersonic turbulence plays a pivotal role during the formation of molecular clouds and stars in galaxies. However, little is known about how the fraction of compressive and solenoidal modes in the velocity field evolves over time and how it depends on properties of the molecular cloud or the galactic environment. In this work, we carry out magnetohydrodynamical simulations of disc galaxies and study the time evolution of the turbulence driving parameter for an ensemble of clouds. We find that the time-averaged turbulence driving parameter is insensitive to the position of the cloud within the galaxy. The ensemble-averaged driving parameter is found to be rather compressive with $b\sim0.5-0.7$, indicating almost time-independent global star formation properties. However, each individual cloud shows a highly fluctuating driving parameter, which would strongly affect the cloud's star formation rate. We find that the mode of turbulence driving can rapidly change within only a few Myr, both from solenoidal to compressive and vice versa. We attribute these changes to cloud collisions and to tidal interactions with clouds or overdensities in the environment. Last, we find no significant differences in the average driving parameter between hydrodynamic and initially strongly magnetised galaxies. However, the magnetic field tends to reduce the overall fluctuation of the driving parameter. The average driving as well as its uncertainty are seen to be in agreement with recent constraints on the turbulence driving mode for solar neighbourhood clouds.

Complex turbulent motions of magnetized gas are ubiquitous in the interstellar medium. The source of this turbulence, however, is still poorly understood. Previous work suggests that compression caused by supernova shockwaves, gravity, or cloud collisions, may drive the turbulence to some extent. In this work, we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of contraction in turbulent, magnetized clouds from the warm neutral medium (WNM) of the ISM to the formation of cold dense molecular clouds, including radiative heating and cooling. We study different contraction rates and find that observed molecular cloud properties, such as the temperature, density, Mach number, and magnetic field strength, and their respective scaling relations, are best reproduced when the contraction rate equals the turbulent turnover rate. In contrast, if the contraction rate is significantly larger (smaller) than the turnover rate, the compression drives too much (too little) turbulence, producing unrealistic cloud properties. We find that the density probability distribution function evolves from a double log-normal representing the two-phase ISM, to a skewed, single log-normal in the dense, cold phase. For purely hydrodynamical simulations, we find that the effective driving parameter of contracting cloud turbulence is natural to mildly compressive (\mbox{$b\sim0.4$--$0.5$}), while for MHD turbulence, we find \mbox{$b\sim0.3$--$0.4$}, i.e., solenoidal to naturally mixed. Overall, the physical properties of the simulated clouds that contract at a rate equal to the turbulent turnover rate, indicate that large-scale contraction may explain the origin and evolution of turbulence in the ISM.

Magnetic fields are an elemental part of the interstellar medium in galaxies. However, their impact on gas dynamics and star formation in galaxies remains controversial. We use a suite of global magnetohydrodynamical simulations of isolated disc galaxies to study the influence of magnetic fields on the diffuse and dense gas in the discs. We find that the magnetic field acts in multiple ways. Stronger magnetised discs fragment earlier due to the shorter growth time of the Parker instability. Due to the Parker instability in the magnetised discs we also find cold ($T<50\,\mathrm{K}$) and dense ($n\sim10^3-10^4\,\mathrm{cm}^{-3}$) gas several hundred pc above/below the midplane without any form of stellar feedback. In addition, magnetic fields change the fragmentation pattern. While in the hydrodynamical case, the disc breaks up into ring-like structures, magnetised discs show the formation of filamentary entities that extent both in the azimuthal and radial direction. These kpc scale filaments become magnetically (super-)critical very quickly and allow for the rapid formation of massive giant molecular clouds. Our simulations suggest that major differences in the behaviour of star formation - due to a varying magnetisation - in galaxies could arise.

We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^\alpha$, with $\alpha$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models.

Both observational and theoretical research over the past decade has demonstrated that the probability distribution function (PDF) of the gas density in turbulent molecular clouds is a key ingredient for understanding star formation. It has recently been argued that the PDF of molecular clouds is a pure power-law distribution. It has been claimed that the log-normal part is ruled out when using only the part of the PDF up/down to which it is complete, that is where the column density contours are still closed. By using the results from high-resolution magnetohydrodynamical simulations of molecular cloud formation and evolution, we find that the column density PDF is indeed composed of a log-normal and, if including self-gravity, a power-law part. We show that insufficient sampling of a molecular cloud results in closed contours that cut off the log-normal part. In contrast, systematically increasing the field of view and sampling the entire cloud yields a completeness limit at the lower column densities, which also recovers the log-normal part. This demonstrates that the field of view must be sufficiently large for the PDF to be complete down to its log-normal part, which has important implications for predictions of star-formation activity based on the PDF.

Observations show that galaxies and their interstellar media are pervaded by strong magnetic fields with energies in the diffuse component being at least comparable to the thermal and even as large or larger than the turbulent energy. Such strong magnetic fields prevent the formation of stars because patches of the interstellar medium are magnetically subcritical. Here we present the results from global numerical simulations of strongly magnetised and self-gravitating galactic discs, which show that the buoyancy of the magnetic field due to the Parker instability leads at first to the formation of giant filamentary regions. These filamentary structures become gravitationally unstable and fragment into $\sim10^5 M_{\odot}$ clouds that attract kpc long, coherent filamentary flows that build them into GMCs. Our results thus provide a solution to the long-standing problem of how the transition from sub- to supercritical regions in the interstellar medium proceeds.

Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments

Astrochemistry plays a central role during the process of star formation, both in the primordial regime as well as in the present-day Universe. We revisit here the chemistry in both regimes, focusing first on the chemistry under close to primordial conditions, as observed in the so-called Caffau star SDSS J102915+172927, and subsequently discuss deuteration processes in present-day star-forming cores. In models of the high-redshift Universe, the chemistry is particularly relevant to determine the cooling, while it also serves as an important diagnostic in the case of present-day star formation.

Molecular clouds are the principle stellar nurseries of our universe, keeping them in the focus of both observational and theoretical studies. From observations, some of the key properties of molecular clouds are well known but many questions regarding their evolution and star formation activity remain open. While numerical simulations feature a large number and complexity of involved physical processes, this plenty of effects may hide the fundamentals that determine the evolution of molecular clouds and enable the formation of stars. Purely analytical models, on the other hand, tend to suffer from rough approximations or a lack of completeness, limiting their predictive power. In this paper, we present a model that incorporates central concepts of astrophysics as well as reliable results from recent simulations of molecular clouds and their evolutionary paths. Based on that, we construct a self-consistent semi-analytical framework that describes the formation, evolution and star formation activity of molecular clouds, including a number of feedback effects to account for the complex processes inside those objects. The final equation system is solved numerically but at much lower computational expense than, e.g., hydrodynamical descriptions of comparable systems. The model presented in this paper agrees well with a broad range of observational results, showing that molecular cloud evolution can be understood as an interplay between accretion, global collapse, star formation and stellar feedback.

Molecular clouds are to a great extent influenced by turbulent motions in the gas. Numerical and observational studies indicate that the star formation rate and efficiency crucially depend on the mixture of solenoidal and compressive modes in the turbulent acceleration field, which can be quantified by the turbulent driving parameter b. For purely solenoidal (divergence-free) driving previous studies showed that b=1/3 and for entirely compressive (curl-free) driving b=1. In this study, we determine the evolution of the turbulent driving parameter b in magnetohydrodynamical simulations of molecular cloud formation and evolution. The clouds form due to the convergence of two flows of warm neutral gas. We explore different scenarios by varying the magnitude of the initial turbulent perturbations in the flows. We show that the driving mode of the turbulence within the cloud strongly fluctuates with time and exhibits no clear correlation with typical cloud properties, such as the cloud mass and the (Alfven) Mach number. We specifically find that $b$ strongly varies from b=0.3 to b=0.8 on timescales t<5 Myr, where the timescale and range of variation can change from cloud to cloud. This rapid change of b from solenoidal to compressive driving is primarily associated with global contraction of the cloud and subsequent onset of star formation. We conclude that the effective turbulence driving parameter should be treated as a free parameter that can vary from solenoidal to compressive in both time and space.

High-mass stars are expected to form from dense prestellar cores. Their precise formation conditions are widely discussed, including their virial condition, which results in slow collapse for super-virial cores with strong support by turbulence or magnetic fields, or fast collapse for sub-virial sources. To disentangle their formation processes, measurements of the deuterium fractions are frequently employed to approximately estimate the ages of these cores and to obtain constraints on their dynamical evolution. We here present 3D magneto-hydrodynamical simulations including for the first time an accurate non-equilibrium chemical network with 21 gas-phase species plus dust grains and 213 reactions. With this network we model the deuteration process in fully depleted prestellar cores in great detail and determine its response to variations in the initial conditions. We explore the dependence on the initial gas column density, the turbulent Mach number, the mass-to-magnetic flux ratio and the distribution of the magnetic field, as well as the initial ortho-to-para ratio of H2. We find excellent agreement with recent observations of deuterium fractions in quiescent sources. Our results show that deuteration is rather efficient, even when assuming a conservative ortho-to-para ratio of 3 and highly sub-virial initial conditions, leading to large deuterium fractions already within roughly a free-fall time. We discuss the implications of our results and give an outlook to relevant future investigations.

We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer timescales, supernova feedback increases the amount of gas with moderate temperatures (T~300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4.

Molecular clouds, the birthplaces of stars in galaxies, form dynamically from the diffuse atomic gas of the interstellar medium (ISM). The ISM is also threaded by magnetic fields which have a large impact on its dynamics. In particular, star forming regions must be magnetically supercrit- ical in order to accomodate gas clumps which can collapse under their own weight. Based on a parameter study of three dimensional magneto-hydrodyamical (MHD) simulations, we show that the long-standing problem of how such supercritical regions are generated is still an open issue.

We use magnetohydrodynamical simulations of converging flows to investigate the process of molecular cloud formation and evolution out of the magnetised ISM. Here, we investigate whether the observed subcritical HI clouds can become supercritical and hence allow the formation of stars within them. To do so, we vary the turbulent Mach number of the flows, as well as the initial magnetic field strength. We show that dense cores are able to build up under all conditions, but that star formation in these cores is either heavily delayed or completely suppressed if the initial field strength is B>3 microGauss. To probe the effect of magnetic diffusion, we introduce a tilting angle between the flows and the uniform background magnetic field, which mimics non--ideal MHD effects. Even with highly inclined flows, the formed cores are devoid of star formation, because no magnetically supercritical regions are build up. Hence we conclude, that the problem of how supercritical cloud cores are generated still persists.

We present an analysis of probability distribution functions (pdfs) of column density in different zones of the star-forming region Perseus and its diffuse environment based on the map of dust opacity at 353 GHz available from the Planck archive. The pdf shape can be fitted by a combination of a lognormal function and an extended power-law tail at high densities, in zones centred at the molecular cloud Perseus. A linear combination of several lognormals fits very well the pdf in rings surrounding the cloud or in zones of its diffuse neighbourhood. The slope of the mean density scaling law $\langle\rho\rangle_L \propto L^\alpha$ is steep ($\alpha=-1.93$) in the former case and rather shallow ($\alpha=-0.77\pm0.11$) in the rings delineated around the cloud. We interpret these findings as signatures of two distinct physical regimes: i) a gravoturbulent one which is characterized by nearly linear scaling of mass and practical lack of velocity scaling; and ii) a predominantly turbulent one which is best described by steep velocity scaling and by invariant for compressible turbulence $\langle\rho\rangle_L u_L^3/L$, describing a scale-independent flux of the kinetic energy per unit volume through turbulent cascade. The gravoturbulent spatial domain can be identified with the molecular cloud Perseus while a relatively sharp transition to predominantly turbulent regime occurs in its vicinity.