Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around the first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called \textit{small-scale turbulent dynamo} and the \textit{large-scale mean-field dynamo}. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.

Explaining the initial mass function (IMF) of stars is a long-standing problem in astrophysics. The number of complex mechanisms involved in the process of star cluster formation, such as turbulence, magnetic fields and stellar feedback, make understanding and modeling the IMF a challenging task. In this paper, we aim to assert the importance of stellar heating feedback in the star cluster formation process and its effect on the shape of the IMF. We use an analytical sub-grid model to implement the radiative feedback in fully three-dimensional magnetohydrodynamical (MHD) simulations of star cluster formation, with the ultimate objective of obtaining numerical convergence on the IMF. We compare a set of MHD adaptive-mesh-refinement (AMR) simulations with three different implementations of the heating of the gas: 1) a polytropic equation of state (EOS), 2) a spherically symmetric stellar heating feedback, and 3) our newly developed polar heating model that takes into account the geometry of the accretion disc and the resulting shielding of stellar radiation by dust. For each of the three heating models, we analyse the distribution of stellar masses formed in ten molecular cloud simulations with different realizations of the turbulence to obtain a statistically representative IMF. We conclude that stellar heating feedback has a profound influence on the number of stars formed and plays a crucial role in controlling the IMF. We find that the simulations with the polar heating model achieve the best convergence on the observed IMF.

Magnetic fields play an important role for the formation of stars in both local and high-redshift galaxies. Recent studies of dynamo amplification in the first dark matter haloes suggest that significant magnetic fields were likely present during the formation of the first stars in the Universe at redshifts of 15 and above. In this work, we study how these magnetic fields potentially impact the initial mass function (IMF) of the first stars. We perform 200 high-resolution, three-dimensional (3D), magneto-hydrodynamic (MHD) simulations of the collapse of primordial clouds with different initial turbulent magnetic field strengths as predicted from turbulent dynamo theory in the early Universe, forming more than 1100 first stars in total. We detect a strong statistical signature of suppressed fragmentation in the presence of strong magnetic fields, leading to a dramatic reduction in the number of first stars with masses low enough that they might be expected to survive to the present day. Additionally, strong fields shift the transition point where stars go from being mostly single to mostly multiple to higher masses. However, irrespective of the field strength, individual simulations are highly chaotic, show different levels of fragmentation and clustering, and the outcome depends on the exact realisation of the turbulence in the primordial clouds. While these are still idealised simulations that do not start from cosmological initial conditions, our work shows that magnetic fields play a key role for the primordial IMF, potentially even more so than for the present-day IMF.

We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey, we find the vertical velocity dispersion ($\sigma_{v, z}$) to be positively correlated with measures of star-formation rate, stellar mass, HI gas mass, and rotational velocity. The greatest correlation is with star-formation rate surface density ($\Sigma_\text{SFR}$). Using the total sample, we find $\sigma_{v, z}$ increases slowly as a function of integrated star-formation rate in the range SFR $\in$ [$10^{-3}$, 1] M$_\odot$ yr$^{-1}$ from $17\pm3$ km s$^{-1}$ to $24\pm5$ km s$^{-1}$ followed by a steeper increase up to $\sigma_{v, z}$ $\sim 80$ km s$^{-1}$ for SFR $\gtrsim 1$ M$_\odot$ yr$^{-1}$. This is consistent with recent theoretical models that suggest a $\sigma_{v, z}$ floor driven by star-formation feedback processes with an upturn in $\sigma_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc.

Globular cluster progenitors may have been detected by \textit{HST}, and are predicted to be observable with \textit{JWST} and ground-based extremely-large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched (${\Delta}Y=+0.12$) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first $\sim$30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being $\sim$0.40 mag brighter at $\lambda \approx 2.0\, {\mu}m$, and to be 0.30 mag redder in the \textit{JWST}-NIRCam colour $(F070W-F200W)$, and to actually be fainter for $\lambda \lesssim 0.50 \, {\mu}m$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim 5\%$, though we show that the Lyman-$\alpha$ escape fraction could end up higher for helium-enriched stars.

We take advantage of a set of molecular cloud simulations to demonstrate a possibility to uncover statistical properties of the gas density and velocity fields using reflected emission of a short (with duration much less than the cloud's light-crossing time) X-ray flare. Such situation is relevant for the Central Molecular Zone of our Galaxy where several clouds get illuminated by a $\sim110$ yr-old flare from the supermassive black hole Sgr A*. Due to shortness of the flare ($\Delta t\lesssim1.6$ yrs), only a thin slice ($\Delta z\lesssim0.5$ pc) of the molecular gas contributes to the X-ray reflection signal at any given moment, and its surface brightness effectively probes the local gas density. This allows reconstructing the density probability distribution function over a broad range of scales with virtually no influence of attenuation, chemo-dynamical biases and projection effects. Such measurement is key to understanding the structure and star-formation potential of the clouds evolving under extreme conditions in the CMZ. For cloud parameters similar to the currently brightest in X-ray reflection molecular complex Sgr A, the sensitivity level of the best available data is sufficient only for marginal distinction between solenoidal and compressive forcing of turbulence. Future-generation X-ray observatories with large effective area and high spectral resolution will dramatically improve on that by minimising systematic uncertainties due to contaminating signals. Furthermore, measurement of the iron fluorescent line centroid with sub-eV accuracy in combination with the data on molecular line emission will allow direct investigation of the gas velocity field.

Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galaxy, coupled with photoionisation models, to produce a suite of synthetic high resolution integral field spectroscopy (IFS) datacubes. We then degrade the datacubes, with a range of realistic models for spatial resolution (2 to 16 beams per galaxy scale length) and noise, to investigate and quantify how well the input metallicity gradient can be recovered as a function of resolution and signal-to-noise ratio (SNR) with the intention to compare with modern IFS surveys like MaNGA and SAMI. Given appropriate propagation of uncertainties and pruning of low SNR pixels, we show that a resolution of 3-4 telescope beams per galaxy scale length is sufficient to recover the gradient to ~10-20% uncertainty. The uncertainty escalates to ~60% for lower resolution. Inclusion of the low SNR pixels causes the uncertainty in the inferred gradient to deteriorate. Our results can potentially inform future IFS surveys regarding the resolution and SNR required to achieve a desired accuracy in metallicity gradient measurements.

Aims: Investigate how the strength of episodic accretion bursts depends on eccentricity. Methods: We investigate the binary trigger hypothesis in longer-period (>20yr) binaries by carrying out three-dimensional magnetohydrodynamical (MHD) simulations of the formation of low-mass binary stars down to final separations of ~10AU, including the effects of gas turbulence and magnetic fields. We run two simulations with an initial turbulent gas core of one solar mass each and two different initial turbulent Mach numbers, M = 0.1 and 0.2, for 6500yr after protostar formation. Results: We observe bursts of accretion at periastron during the early stages when the eccentricity of the binary system is still high. We find that this correlation between bursts of accretion and passing periastron breaks down at later stages, because of the gradual circularisation of the orbits. For eccentricities greater than e=0.2, we observe episodic accretion triggered near periastron. However, we do not find any strong correlation between the strength (the ratio of the burst accretion ratio to the quiescent accretion rate) of episodic accretion and eccentricity. We determine that accretion events are likely triggered by torques between the rotation of the circumstellar disc and the approaching binary stars. We compare our results with observational data of episodic accretion in short-period binaries and find good agreement between our simulations and the observations. Conclusions: We conclude that episodic accretion is a universal mechanism operating in eccentric young binary-star systems, independent of separation, and should be observable in long-period binaries as well as in short-period binaries, but that the strength will depend on the torques, and hence the separation at periastron.

The gas motions in the intracluster medium (ICM) are governed by stratified turbulence. Stratified turbulence is fundamentally different from Kolmogorov (isotropic, homogeneous) turbulence; kinetic energy not only cascades from large to small scales, but it is also converted into buoyancy potential energy. To understand the density and velocity fluctuations in the ICM, we conduct high-resolution ($1024^2\times 1536$ grid points) hydrodynamical simulations of subsonic turbulence (with rms Mach number $\mathcal{M}\approx 0.25$) and different levels of stratification, quantified by the Richardson number $\mathrm{Ri}$, from $\mathrm{Ri}=0$ (no stratification) to $\mathrm{Ri}=13$ (strong stratification). We quantify the density, pressure and velocity fields for varying stratification because observational studies often use surface brightness fluctuations to infer the turbulent gas velocities of the ICM. We find that the standard deviation of the logarithmic density fluctuations ($\sigma_s$), where $s=\ln(\rho/\left<\rho(z)\right>)$, increases with $\mathrm{Ri}$. For weakly stratified subsonic turbulence ($\mathrm{Ri}\lesssim10$, $\mathcal{M}<1$), we derive a new $\sigma_s$--$\mathcal{M}$--$\mathrm{Ri}$ relation, $\sigma_s^2=\ln(1+b^2\mathcal{M}^4+0.09\mathcal{M}^2\mathrm{Ri}H_P/H_S)$, where $b=1/3$--$1$ is the turbulence driving parameter, and $H_P$ and $H_S$ are the pressure and entropy scale heights respectively. We further find that the power spectrum of density fluctuations, $P(\rho_k/\left<\rho\right>)$, increases in magnitude with increasing $\mathrm{Ri}$, whereas the velocity power spectrum is invariant. Thus, the ratio between density and velocity power spectra strongly depends on $\mathrm{Ri}$. Pressure fluctuations, on the other hand, are independent of stratification and only depend on $\mathcal{M}$.

The physics regulating star formation (SF) in Hickson Compact Groups (HCG) has thus far been difficult to describe, due to their unique kinematic properties. In this study we expand upon previous works to devise a more physically meaningful SF relation able to better encompass the physics of these unique systems. We combine CO(1--0) data from the Combined Array from Research in Millimeter Astronomy (CARMA) to trace the column density of molecular gas $\Sigma_\mathrm{gas}$ and deep H$\alpha$ imaging taken on the Southern Astrophysical Research (SOAR) Telescope tracing $\Sigma_\mathrm{SFR}$ to investigate star formation efficiency across face-on HCG, NGC7674. We find a lack of universality in star formation, with two distinct sequences present in the $\Sigma_\mathrm{gas}-\Sigma_\mathrm{SFR}$ plane; one for inside and one for outside the nucleus. We devise a SF relation based on the multi-freefall nature of gas and the critical density, which itself is dependent on the virial parameter $\alpha_{\mathrm{vir}}$, the ratio of turbulent to gravitational energy. We find that our modified SF relation fits the data and describes the physics of this system well with the introduction of a virial parameter of about 5--10 across the galaxy. This $\alpha_{\mathrm{vir}}$ leads to an order-of-magnitude reduction in SFR compared to $\alpha_{\mathrm{vir}}\approx 1$ systems.

We investigate the turbulence driving mode of ionizing radiation from massive stars on the surrounding interstellar medium (ISM). We run hydrodynamical simulations of a turbulent cloud impinged by a plane-parallel ionization front. We find that the ionizing radiation forms pillars of neutral gas reminiscent of those seen in observations. We quantify the driving mode of the turbulence in the neutral gas by calculating the driving parameter $b$, which is characterised by the relation $\sigma_s^2 = \ln({1+b^2\mathcal{M}^2})$ between the variance of the logarithmic density contrast $\sigma_s^2$ (where $s = \ln({\rho/\rho_0})$ with the gas density $\rho$ and its average $\rho_0$), and the turbulent Mach number $\mathcal{M}$. Previous works have shown that $b\sim1/3$ indicates solenoidal (divergence-free) driving and $b\sim1$ indicates compressive (curl-free) driving, with $b\sim1$ producing up to ten times higher star formation rates than $b\sim1/3$. The time variation of $b$ in our study allows us to infer that ionizing radiation is inherently a compressive turbulence driving source, with a time-averaged $b\sim 0.76 \pm 0.08$. We also investigate the value of $b$ of the pillars, where star formation is expected to occur, and find that the pillars are characterised by a natural mixture of both solenoidal and compressive turbulent modes ($b\sim0.4$) when they form, and later evolve into a more compressive turbulent state with $b\sim0.5$--$0.6$. A virial parameter analysis of the pillar regions supports this conclusion. This indicates that ionizing radiation from massive stars may be able to trigger star formation by producing predominately compressive turbulent gas in the pillars.

The complexity of modern and upcoming computing architectures poses severe challenges for code developers and application specialists, and forces them to expose the highest possible degree of parallelism, in order to make the best use of the available hardware. The Intel$^{(R)}$ Xeon Phi$^{(TM)}$ of second generation (code-named Knights Landing, henceforth KNL) is the latest many-core system, which implements several interesting hardware features like for example a large number of cores per node (up to 72), the 512 bits-wide vector registers and the high-bandwidth memory. The unique features of KNL make this platform a powerful testbed for modern HPC applications. The performance of codes on KNL is therefore a useful proxy of their readiness for future architectures. In this work we describe the lessons learnt during the optimisation of the widely used codes for computational astrophysics P-Gadget-3, Flash and Echo. Moreover, we present results for the visualisation and analysis tools VisIt and yt. These examples show that modern architectures benefit from code optimisation at different levels, even more than traditional multi-core systems. However, the level of modernisation of typical community codes still needs improvements, for them to fully utilise resources of novel architectures.

Motivated by observations of outflowing galaxies, we investigate the combined impact of magnetic fields and radiative cooling on the evolution of cold clouds embedded in a hot wind. We perform a collection of three-dimensional adaptive mesh refinement, magnetohydrodynamical simulations that span two resolutions, and include fields that are aligned and transverse to the oncoming, super-Alfv\'enic material. Aligned fields have little impact on the overall lifetime of the clouds over the non-magnetized case, although they do increase the mixing between the wind and cloud material by a factor of $\approx 3.$ Transverse fields lead to magnetic draping, which isolates the clouds, but they also squeeze material in the direction perpendicular to the field lines, which leads to rapid mass loss. A resolution study suggests that the magnetized simulations have somewhat better convergence properties than non-magnetized simulations, and that a resolution of 64 zones per cloud radius is sufficient to accurately describe these interactions. We conclude that the combined effects of radiative cooling and magnetic fields are dependent on field orientation, but are unlikely to enhance cloud lifetimes beyond the effect of radiative cooling alone.

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.

Star formation is inefficient. Recent advances in numerical simulations and theoretical models of molecular clouds show that the combined effects of interstellar turbulence, magnetic fields and stellar feedback can explain the low efficiency of star formation. The star formation rate is highly sensitive to the driving mode of the turbulence. Solenoidal driving may be more important in the Central Molecular Zone, compared to more compressive driving agents in spiral-am clouds. Both theoretical and observational efforts are underway to determine the dominant driving mode of turbulence in different Galactic environments. New observations with ALMA, combined with other instruments such as CARMA, JCMT and the SMA begin to reveal the magnetic field structure of dense cores and protostellar disks, showing highly complex field geometries with ordered and turbulent field components. Such complex magnetic fields can give rise to a range of stellar masses and jet/outflow efficiencies in dense cores and protostellar accretion disks.

Stars form in highly-magnetised, supersonic turbulent molecular clouds. Many of the tools and models that we use to carry out star formation studies rely upon the assumption of cloud isotropy. However, structures like high-density filaments in the presence of magnetic fields, and magnetosonic striations introduce anisotropies into the cloud. In this study we use the two-dimensional (2D) power spectrum to perform a systematic analysis of the anisotropies in the column density for a range of Alfv\'en Mach numbers ($\mathcal{M}_A=0.1$--$10$) and turbulent Mach numbers ($\mathcal{M}=2$--$20$), with 20 high-resolution, three-dimensional (3D) turbulent magnetohydrodynamic simulations. We find that for cases with a strong magnetic guide field, corresponding to $\mathcal{M}_A<1$, and $\mathcal{M}\lesssim 4$, the anisotropy in the column density is dominated by thin striations aligned with the magnetic field, while for $\mathcal{M}\gtrsim 4$ the anisotropy is significantly changed by high-density filaments that form perpendicular to the magnetic guide field. Indeed, the strength of the magnetic field controls the degree of anisotropy and whether or not any anisotropy is present, but it is the turbulent motions controlled by $\mathcal{M}$ that determine which kind of anisotropy dominates the morphology of a cloud.

In this exploratory submission we present the visualization of the largest interstellar turbulence simulations ever performed, unravelling key astrophysical processes concerning the formation of stars and the relative role of magnetic fields. The simulations, including pure hydrodynamical (HD) and magneto-hydrodynamical (MHD) runs, up to a size of $10048^3$ grid elements, were produced on the supercomputers of the Leibniz Supercomputing Centre and visualized using the hybrid parallel (MPI+TBB) ray-tracing engine OSPRay associated with VisIt. Besides revealing features of turbulence with an unprecedented resolution, the visualizations brilliantly showcase the stretching-and-folding mechanisms through which astrophysical processes such as supernova explosions drive turbulence and amplify the magnetic field in the interstellar gas, and how the first structures, the seeds of newborn stars are shaped by this process.

The adiabatic index of H$_2\,$ ($\gamma_{\mathrm{H_2}}$) is non-constant at temperatures between $100-10^4\,\mathrm{K}$ due to the large energy spacing between its rotational and vibrational modes. For the formation of the first stars at redshifts 20 and above, this variation can be significant because primordial molecular clouds are in this temperature range due to the absence of efficient cooling by dust and metals. We study the possible importance of variations in $\gamma_{\mathrm{H_2}}$ for the primordial initial mass function by carrying out 80 3D gravito-hydrodynamic simulations of collapsing clouds with different random turbulent velocity fields, half using fixed $\gamma_{\rm H_2} = 7/5$ in the limit of classical diatomic gas (used in earlier works) and half using an accurate quantum mechanical treatment of $\gamma_{\mathrm{H_2}}$. We use the adaptive mesh refinement code FLASH with the primordial chemistry network from KROME for this study. The simulation suite produces almost 400 stars, with masses from $0.02 - 50$ M$_\odot$ (mean mass $\sim 10.5\,\mathrm{M_{\odot}}$ and mean multiplicity fraction $\sim 0.4$). While the results of individual simulations do differ when we change our treatment of $\gamma_{\mathrm{H_2}}$, we find no statistically significant differences in the overall mass or multiplicity distributions of the stars formed in the two sets of runs. We conclude that, at least prior to the onset of radiation feedback, approximating H$_2$ as a classical diatomic gas with $\gamma_{\rm H_2} = 7/5$ does not induce significant errors in simulations of the fragmentation of primordial gas. Nonetheless, we recommend using the accurate formulation of the H$_2\,$ adiabatic index in primordial star formation studies since it is not computationally more expensive and provides a better treatment of the thermodynamics.

Gaia DR2 provides an unprecedented sample of stars with full 6D phase-space measurements, creating the need for a self-consistent means of discovering and characterising the phase-space overdensities known as moving groups or associations. Here we present Chronostar, a new Bayesian analysis tool that meets this need. Chronostar uses the Expectation-Maximisation algorithm to remove the circular dependency between association membership lists and fits to their phase-space distributions, making it possible to discover unknown associations within a kinematic data set. It uses forward-modelling of orbits through the Galactic potential to overcome the problem of tracing backward stars whose kinematics have significant observational errors, thereby providing reliable ages. In tests using synthetic data sets with realistic measurement errors and complex initial distributions, Chronostar successfully recovers membership assignments and kinematic ages up to ~ 100 Myr. In tests on real stellar kinematic data in the phase-space vicinity of the ${\beta}$ Pictoris Moving Group, Chronostar successfully rediscovers the association without any human intervention, identifies 15 new likely members, corroborates 43 candidate members, and returns a kinematic age of 17.8+-1.2 Myr. In the process we also rediscover the Tucana-Horologium Moving Group, for which we obtain a kinematic age of 36.3+1.3-1.4 Myr.

Most gas in giant molecular clouds is relatively low-density and forms star inefficiently, converting only a small fraction of its mass to stars per dynamical time. However, star formation models generally predict the existence of a threshold density above which the process is efficient and most mass collapses to stars on a dynamical timescale. A number of authors have proposed observational techniques to search for a threshold density above which star formation is efficient, but it is unclear which of these techniques, if any, are reliable. In this paper we use detailed simulations of turbulent, magnetised star-forming clouds, including stellar radiation and outflow feedback, to investigate whether it is possible to recover star formation thresholds using current observational techniques. Using mock observations of the simulations at realistic resolutions, we show that plots of projected star formation efficiency per free-fall time $\epsilon_{\rm ff}$ can detect the presence of a threshold, but that the resolutions typical of current dust emission or absorption surveys are insufficient to determine its value. In contrast, proposed alternative diagnostics based on a change in the slope of the gas surface density versus star formation rate surface density (Kennicutt-Schmidt relation) or on the correlation between young stellar object counts and gas mass as a function of density are ineffective at detecting thresholds even when they are present. The signatures in these diagnostics sometimes taken as indicative of a threshold in observations, which we generally reproduce in our mock observations, do not prove to correspond to real physical features in the 3D gas distribution.

We examine the diagnostic power of rest-frame ultraviolet (UV) nebular emission lines, and compare them to more commonly used rest-frame optical emission lines, using the test case of a single star-forming knot of the bright lensed galaxy RCSGA 032727-132609 at redshift z~1.7. This galaxy has complete coverage of all the major rest-frame UV and optical emission lines from Magellan/MagE and Keck/NIRSPEC. Using the full suite of diagnostic lines, we infer the physical properties: nebular electron temperature (T_e), electron density (n_e), oxygen abundance (log(O/H)), ionisation parameter (log(q)) and interstellar medium (ISM) pressure (log(P/k)). We examine the effectiveness of the different UV, optical and joint UV-optical spectra in constraining the physical conditions. Using UV lines alone we can reliably estimate log(q), but the same is difficult for log(O/H). UV lines yield a higher (~1.5 dex) log(P/k) than the optical lines, as the former probes a further inner nebular region than the latter. For this comparison, we extend the existing Bayesian inference code IZI, adding to it the capability to infer ISM pressure simultaneously with metallicity and ionisation parameter. This work anticipates future rest-frame UV spectral datasets from the James Webb Space Telescope (JWST) at high redshift and from the Extremely Large Telescope (ELT) at moderate redshift.

Supersonic turbulence is a key player in controlling the structure and star formation potential of molecular clouds (MCs). The three-dimensional (3D) turbulent Mach number, $\mathcal{M}$, allows us to predict the rate of star formation. However, determining Mach numbers in observations is challenging because it requires accurate measurements of the velocity dispersion. Moreover, observations are limited to two-dimensional (2D) projections of the MCs and velocity information can usually only be obtained for the line-of-sight component. Here we present a new method that allows us to estimate $\mathcal{M}$ from the 2D column density, $\Sigma$, by analysing the fractal dimension, $\mathcal{D}$. We do this by computing $\mathcal{D}$ for six simulations, ranging between $1$ and $100$ in $\mathcal{M}$. From this data we are able to construct an empirical relation, $\log\mathcal{M}(\mathcal{D}) = \xi_1(\text{erfc}^{-1} [(\mathcal{D}-\mathcal{D}_{\text{min}})/\Omega] + \xi_2),$ where $\text{erfc}^{-1}$ is the inverse complimentary error function, $\mathcal{D}_{\text{min}} = 1.55 \pm 0.13$ is the minimum fractal dimension of $\Sigma$, $\Omega = 0.22 \pm 0.07$, $\xi_1 = 0.9 \pm 0.1$ and $\xi_2 = 0.2 \pm 0.2$. We test the accuracy of this new relation on column density maps from $Herschel$ observations of two quiescent subregions in the Polaris Flare MC, `saxophone' and `quiet'. We measure $\mathcal{M} \sim 10$ and $\mathcal{M} \sim 2$ for the subregions, respectively, which is similar to previous estimates based on measuring the velocity dispersion from molecular line data. These results show that this new empirical relation can provide useful estimates of the cloud kinematics, solely based upon the geometry from the column density of the cloud.

We probe the star formation properties of the gas in AzTEC-1 in the COSMOS field, one of the best resolved and brightest starburst galaxies at $z \approx 4.3$, forming stars at a rate > 1000 $\mathrm{M_{\odot}}\,\mathrm{yr^{-1}}$. Using recent ALMA observations, we study star formation in the galaxy nucleus and an off-center star-forming clump and measure a median star formation rate (SFR) surface density of $\Sigma^{\mathrm{nucleus}}_{\mathrm{SFR}} = 270\pm54$ and $\Sigma^{\mathrm{sfclump}}_{\mathrm{SFR}} = 170\pm38\,\mathrm{M_{\odot}}\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$, respectively. Following the analysis by Sharda et al. (2018), we estimate the molecular gas mass, freefall time and turbulent Mach number in these regions to predict $\Sigma_{\mathrm{SFR}}$ from three star formation relations in the literature. The Kennicutt-Schmidt (Kennicutt 1998, KS) relation, which is based on the gas surface density, underestimates the $\Sigma_{\mathrm{SFR}}$ in these regions by a factor 2-3. The $\Sigma_{\mathrm{SFR}}$ we calculate from the single-freefall model of Krumholz et al. 2012 (KDM) is consistent with the measured $\Sigma_{\mathrm{SFR}}$ in the nucleus and the star-forming clump within the uncertainties. The turbulence-regulated star formation relation by Salim et al. 2015 (SFK) agrees slightly better with the observations than the KDM relation. Our analysis reveals that an interplay between turbulence and gravity can help sustain high SFRs in high-redshift starbursts. It can also be extended to other high- and low-redshift galaxies thanks to the high angular resolution and sensitivity of ALMA observations.

Sharpless 2-27 (Sh2-27) is a nearby ${\rm H\,{\small II}}$ region excited by $\zeta$Oph. We present observations of polarized radio emission from 300 to 480$\,$MHz towards Sh2-27, made with the Parkes 64$\,$m Radio Telescope as part of the Global Magneto-Ionic Medium Survey. These observations have an angular resolution of $1.35^{\circ}$, and the data are uniquely sensitive to magneto-ionic structure on large angular scales. We demonstrate that background polarized emission towards Sh2-27 is totally depolarized in our observations, allowing us to investigate the foreground. We analyse the results of Faraday tomography, mapping the magnetised interstellar medium along the 165$\,$pc path to Sh2-27. The Faraday dispersion function in this direction has peaks at three Faraday depths. We consider both Faraday thick and thin models for this observation, finding that the thin model is preferred. We further model this as Faraday rotation of diffuse synchrotron emission in the Local Bubble and in two foreground neutral clouds. The Local Bubble extends for 80$\,$pc in this direction, and we find a Faraday depth of $-0.8 \pm 0.4\,$rad$\,$m$^{-2}$. This indicates a field directed away from the Sun with a strength of $-2.5\pm1.2\,\mu$G. The near and far neutral clouds are each about 30$\,$pc thick, and we find Faraday depths of $-6.6\pm0.6\,$rad$\,$m$^{-2}$ and $+13.7\pm0.8\,$rad$\,$m$^{-2}$, respectively. We estimate that the line-of-sight magnetic strengths in the near and far cloud are $B_{\parallel, \text{near}} \approx -15\,\mu\text{G}$ and $B_{\parallel, \text{far}} \approx +30\,\mu\text{G}$. Our results demonstrate that Faraday tomography can be used to investigate the magneto-ionic properties of foreground features in front of nearby ${\rm H\,{\small II}}$ regions.

Observations of interstellar gas clouds are typically limited to two-dimensional (2D) projections of the intrinsically three-dimensional (3D) structure of the clouds. In this study, we present a novel method for relating the 2D projected fractal dimension ($\mathcal{D}_{\text{p}}$) to the 3D fractal dimension ($\mathcal{D}_{\text{3D}}$) of turbulent clouds. We do this by computing the fractal dimension of clouds over two orders of magnitude in turbulent Mach number $(\mathcal{M} = 1-100)$, corresponding to seven orders of magnitude in spatial scales within the clouds. This provides us with the data to create a new empirical relation between $\mathcal{D}_{\text{p}}$ and $\mathcal{D}_{\text{3D}}$. The proposed relation is $\mathcal{D}_{\text{3D}}(\mathcal{D}_{\text{p}}) = \Omega_1 erfc ( \xi_1 erfc^{-1}[ (\mathcal{D}_{\text{p}} - \mathcal{D}_{\text{p,min}})/\Omega_2 ] + \xi_2 ) + \mathcal{D}_{\text{3D,min}}$, where the minimum 3D fractal dimension, $\mathcal{D}_{\text{3D,min}} = 2.06 \pm 0.35$, the minimum projected fractal dimension, $\mathcal{D}_{\text{p,min}} = 1.55 \pm 0.13$, $\Omega_1 = 0.47 \pm 0.18$, $\Omega_2 = 0.22 \pm 0.07$, $\xi_1 = 0.80 \pm 0.18$ and $\xi_2 = 0.26 \pm 0.19$. The minimum 3D fractal dimension, $\mathcal{D}_{\text{3D,min}} = 2.06 \pm 0.35$, indicates that in the high $\mathcal{M}$ limit the 3D clouds are dominated by planar shocks. The relation between $\mathcal{D}_{\text{p}}$ and $\mathcal{D}_{\text{3D}}$ of molecular clouds may be a useful tool for those who are seeking to understand the 3D structures of molecular clouds, purely based upon 2D projected data and shows promise for relating the physics of the turbulent clouds to the fractal dimension.

Most stars form in binaries and the evolution of their discs remains poorly understood. To shed light on this subject, we carry out 3D ideal MHD simulations with the AMR code FLASH of binary star formation for separations of $10-20\,\mathrm{AU}$. We run a simulation with no initial turbulence (NT), and two with turbulent Mach numbers of $\mathcal{M} = \sigma_v/c_s = 0.1$ and $0.2$ (T1 and T2) for $5000\,\mathrm{yr}$ after protostar formation. By the end of the simulations the circumbinary discs in NT and T1, if any, have radii of $\lesssim20\,\mathrm{AU}$ with masses $\lesssim0.02\,\mathrm{M}_\odot$, while T2 hosts a circumbinary disc with radius $\sim70-80\,\mathrm{AU}$ and mass $\sim0.12\,\mathrm{M}_\odot$. These circumbinary discs formed from the disruption of circumstellar discs and harden the binary orbit. Our simulated binaries launch large single outflows. We find that NT drives the most massive outflows, and also removes large quantities of linear and angular momentum. T2 produces the least efficient outflows concerning mass, momentum and angular momentum ($\sim$61 per cent, $\sim$71 per cent, $\sim$68 per cent of the respective quantities in NT). We conclude that while turbulence helps to build circumbinary discs which organise magnetic fields for efficient outflow launching, too much turbulence may also disrupt the ordered magnetic field structure required for magneto-centrifugal launching of jets and outflows. We also see evidence for episodic accretion during the binary star evolution. We conclude that the role of turbulence in building large circumbinary discs may explain some observed very old ($>10\,\mathrm{Myr}$) circumbinary discs. The longer lifetime of circumbinary discs may increase the likelihood of planet formation.

Recent observations suggest that dense gas clouds can survive even in hot galactic winds. Here we show that the inclusion of turbulent densities with different statistical properties has significant effects on the evolution of wind-swept clouds. We investigate how the initial standard deviation of the log-normal density field influences the dynamics of quasi-isothermal clouds embedded in supersonic winds. We compare uniform, fractal solenoidal, and fractal compressive cloud models in both 3D and 2D hydrodynamical simulations. We find that the processes of cloud disruption and dense gas entrainment are functions of the initial density distribution in the cloud. Fractal clouds accelerate, mix, and are disrupted earlier than uniform clouds. Within the fractal cloud sample, compressive clouds retain high-density nuclei, so they are more confined, less accelerated, and have lower velocity dispersions than their solenoidal counterparts. Compressive clouds are also less prone to Kelvin-Helmholtz and Rayleigh-Taylor instabilities, so they survive longer than solenoidal clouds. By comparing the cloud properties at the destruction time, we find that dense gas entrainment is more effective in uniform clouds than in either of the fractal clouds, and it is more effective in solenoidal than in compressive models. In contrast, mass loading into the wind is more efficient in compressive cloud models than in uniform or solenoidal models. Overall, wide density distributions lead to inefficient entrainment, but they facilitate mass loading and favour the survival of very dense gas in hot galactic winds.

This work presents a study of two Herbig Ae transitional discs, Oph IRS 48 and HD 169142; which both have reported rings in their dust density distributions. We use Keck-II/NIRC2 adaptive optics imaging observations in the L' filter (3.8 micron) to probe the regions of these discs inwards of ~20AU from the star. We introduce our method for investigating these transitional discs, which takes a forward modelling approach: making a model of the disc (using the Monte Carlo radiative transfer code RADMC), convolving it with point-spread functions of calibrator stars, and comparing the convolved models with the observational data. The disc surface density parameters are explored with a Monte Carlo Markov Chain technique. Our analysis recovers emission from both of the discs interior to the well known optically thick walls, modelled as a ring of emission at ~15AU in Oph IRS 48, and ~7AU for HD 169142, and identifies asymmetries in both discs. Given the brightness of the near-symmetric rings compared to the reported companion candidates, we suggest that the reported companion candidates can be interpreted as slightly asymmetric disc emission or illumination.

The atomic interstellar medium (ISM) is observed to be full of linear structures that are referred to as "fibers". Fibers exhibit similar properties to linear structures found in molecular clouds known as striations. Suggestive of a similar formation mechanism, both striations and fibers appear to be ordered, quasi-periodic, and well-aligned with the magnetic field. The prevailing formation mechanism for striations involves the excitation of fast magnetosonic waves. Based on this theoretical model, and through a combination of velocity centroids and column density maps, Tritsis et al. (2018) developed a method for estimating the plane-of-sky (POS) magnetic field from molecular cloud striations. We apply this method in two H\textsc{I} clouds with fibers along the same line-of-sight (LOS) towards the ultra-high-energy cosmic-ray (UHECR) hotspot, at the boundaries of Ursa Major. For the cloud located closer to Earth, where Zeeman observations from the literature were also available, we find general agreement in the distributions of the LOS and POS components of the magnetic field. We find relatively large values for the total magnetic field (ranging from $\sim$$\rm{10}$ to $\sim$$\rm{20} ~\rm{\mu G}$) and an average projection angle with respect to the LOS of $\sim$ 50$^\circ$. For the cloud located further away, we also find a large value for the POS component of the magnetic field of $15^{+8}_{-3}~\rm{\mu G}$. We discuss the potential of our new magnetic-field tomography method for large-scale application. We consider the implications of our findings for the accuracy of current reconstructions of the Galactic magnetic field and on the propagation of UHECR through the ISM.

Magnetic fields are known to play a crucial role in the star formation process, particularly in the formation of jets and outflows from protostellar discs. The magnetic field structure in star forming regions is not always uniform and ordered, often containing regions of magnetic turbulence. We present grid-based, magneto-hydrodynamical simulations of the collapse of a 1$\,\mathrm{M}_{\odot}$ cloud core, to investigate the influence of complex magnetic field structures on outflow formation, morphology and efficiency. We compare three cases: a uniform field, a partially turbulent field and a fully turbulent field, with the same magnetic energy in all three cases. We find that collimated jets are produced in the uniform-field case, driven by a magneto-centrifugal mechanism. Outflows also form in the partially turbulent case, although weaker and less collimated, with an asymmetric morphology. The outflows launched from the partially turbulent case carry the same amount of mass as the uniform-field case but at lower speeds, having only have 71$\%$ of the momentum of the uniform-field case. In the case of a fully turbulent field, we find no significant outflows at all. Moreover, the turbulent magnetic field initially reduces the accretion rate and later induces fragmentation of the disc, forming multiple protostars. We conclude that a uniform poloidal component of the magnetic field is necessary for the driving of jets.

Observations of bright and variable "reflected" X-ray emission from molecular clouds located within inner hundred parsec of our Galaxy have demonstrated that the central supermassive black hole, Sgr A*, experienced short and powerful flares in the past few hundred years. These flares offer a truly unique opportunity to determine 3D location of the illuminated clouds (with ~10 pc accuracy) and to reveal their internal structure (down to 0.1 pc scales). Short duration of the flare(s), combined with X-rays high penetration power and insensitivity of the reflection signal to thermo- and chemo-dynamical state of the gas, ensures that the provided diagnostics of the density and velocity fields is unbiased and almost free of the projection and opacity effects. Sharp and sensitive snapshots of molecular gas accessible with aid of future X-ray observatories featuring large collecting area and high angular (arcsec-level) and spectral (eV-level) resolution cryogenic bolometers will present invaluable information on properties of the supersonic turbulence inside the illuminated clouds, map their shear velocity field and allow cross-matching between X-ray data and velocity-resolved emission of various molecular species provided by ALMA and other ground-based facilities. This will highlight large and small-scale dynamics of the dense gas and help uncovering specifics of the ISM lifecycle and high-mass star formation under very extreme conditions of galactic centers. While the former is of particular importance for the SMBH feeding and triggering AGN feedback, the latter might be an excellent test case for star formation taking place in high-redshift galaxies.

We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H$\alpha$ gas velocity dispersion for our sample to be in the range $\bar{\sigma}_v \sim $[7, 30] km s$^{-1}$. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $\Delta \bar{\sigma}_v / \bar{\sigma}_v = - 0.29 \pm 0.18$ for the H$\alpha$ flux-weighted mean velocity dispersion.

Explaining the internal distribution and motions of stars and gas in galaxies is a key aspect in understanding their evolution. In previous work we identified five well resolved galaxies with atypical kinematics from a cosmological simulation; two had kinematically distinct cores (KDCs), and three had counter-rotating gas and stars (CRGD). In this paper, we show that i) the KDC galaxies have flattening of stellar [O/Fe] at large galacto-centric radii due to the minor mergers that gave rise to the KDCs, and ii) the CRGD galaxies have an abrupt transition in the gas metallicity maps, from high metallicity in the centre to very low metallicity further out. These galaxies are embedded in dark matter filaments where there is a ready supply of near-pristine gas to cause this effect. The non-linear increase in gas metallicity is also seen in the radial profiles, but when the metallicity gradients are measured, the difference is buried in the scatter of the relation. We also find that all five galaxies are fairly compact, with small effective radii given their stellar masses. This is because they have not experienced major mergers that kinematically heat the stars, and would have destroyed their unusual kinematics. In order to detect these signatures of minor mergers or accretion, the galaxy scaling relations or radial metallicity profiles are not enough, and it is necessary to obtain the 2D maps with integral field spectroscopy observations.

Star-forming gas clouds are strongly magnetized, and their ionization fractions are high enough to place them close to the regime of ideal magnetohydrodyamics on all but the smallest size scales. In this review we discuss the effects of magnetic fields on the star formation rate (SFR) in these clouds, and on the mass spectrum of the fragments that are the outcome of the star formation process, the stellar initial mass function (IMF). Current numerical results suggest that magnetic fields by themselves are minor players in setting either the SFR or the IMF, changing star formation rates and median stellar masses only by factors of $\sim 2-3$ compared to non-magnetized flows. However, the indirect effects of magnetic fields, via their interaction with star formation feedback in the form of jets, photoionization, radiative heating, and supernovae, could have significantly larger effects. We explore evidence for this possibility in current simulations, and suggest avenues for future exploration, both in simulations and observations.

We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}{\delta}$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with stellar mass ${\rm log(}M_*/M_{\odot})>10$. The HDSGs make up a significant fraction of non-passive cluster galaxies (15%; 17/115) and a smaller fraction (2.0%; 8/387) of the non-passive population in low-density environments. The majority (9/17) of cluster HDSGs show evidence for star formation at their centers, with the HDS regions found in the outer parts of the galaxy. Conversely, the $\rm{H}{\delta}$-strong signal is more evenly spread across the galaxy for the majority (6/8) of HDSGs in low-density environments, and is often associated with emission lines that are not due to star formation. We investigate the location of the HDSGs in the clusters, finding that they are exclusively within 0.6$R_{200}$ of the cluster centre, and have a significantly higher velocity dispersion relative to the cluster population. Comparing their distribution in projected-phase-space to those derived from cosmological simulations indicates that the cluster HDSGs are consistent with an infalling population that have entered the central 0.5$r_{200, 3D}$ cluster region within the last $\sim 1\,$Gyr. In the 8/9 cluster HDSGs with central star formation, the extent of star formation is consistent with that expected of outside-in quenching by ram-pressure stripping. Our results indicate that the cluster HDSGs are currently being quenched by ram-pressure stripping on their first passage through the cluster.

Most cosmological simulations of galaxy evolution include active galactic nucleus (AGN) feedback, typically seeding black holes with masses of $\geq 10^5\, h^{-1}\, \rm{M}_{\odot}$ when the dark matter halo exceeds a given threshold mass. Taylor & Kobayashi (2014) introduced a new model, which seeds black holes at $10^3\, h^{-1}\, \rm{M}_{\odot}$ based on gas properties alone, and motivated by the channel of black hole formation due to the collapse of the most massive first stars in the Universe. We compare the black hole mass when the dark matter halo mass is $10^{10}\, h^{-1}\, \rm{M}_{\odot}$ between the different seeding methods. We find that seeding based upon gas properties gives a distribution of black hole masses with $\langle \log M_{\rm{BH}} {/ \rm{M}_{\odot}} \rangle = (5.18 \pm 0.54)$ when dark matter halo mass is $10^{10}\, h^{-1}\, \rm{M}_{\odot}$, consistent with the {seeding criteria} used in other simulations. However, the evolution of individual galaxies can be strongly affected by the different seeding mechanisms. We also find that the mean value of the distribution of black hole masses at a given halo mass evolves over time, with higher masses at higher redshifts, indicative of downsizing. Our results can inform more physically motivated black hole and AGN feedback models in cosmological simulations and semi-analytic models.

How stars are born from clouds of gas is a rich physics problem whose solution will inform our understanding of not just stars but also planets, galaxies, and the universe itself. Star formation is stupendously inefficient. Take the Milky Way. Our galaxy contains about a billion solar masses of fresh gas available to form stars-and yet it produces only one solar mass of new stars a year. Accounting for that inefficiency is one of the biggest challenges of modern astrophysics. Why should we care about star formation? Because the process powers the evolution of galaxies and sets the initial conditions for planet formation and thus, ultimately, for life.

We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). The mean integrated specific star formation rate of star-forming galaxies in high-mass groups is lower than for galaxies in low-mass groups or that are ungrouped, with $\Delta \log(sSFR/\mathrm{yr^{-1}}) = 0.45 \pm 0.07$. This difference is seen at all galaxy stellar masses. In high-mass groups, star-forming galaxies more massive than $M_{*} \sim 10^{10} \, \mathrm{M_{\odot}}$ have centrally-concentrated star formation. These galaxies also lie below the star-formation main sequence, suggesting they may be undergoing outside-in quenching. Lower mass galaxies in high-mass groups do not show evidence of concentrated star formation. In groups less massive than $M_{G} = 10^{12.5} \, \mathrm{M_{\odot}}$ we do not observe these trends. In this regime we find a modest correlation between centrally-concentrated star formation and an enhancement in total star formation rate, consistent with triggered star formation in these galaxies.

Misalignment of gas and stellar rotation in galaxies can give clues to the origin and processing of accreted gas. Integral field spectroscopic observations of 1213 galaxies from the SAMI Galaxy Survey show that 11% of galaxies with fitted gas and stellar rotation are misaligned by more than 30 degrees in both field/group and cluster environments. Using SAMI morphological classifications and Sersic indices, the misalignment fraction is 45+/-6% in early-type galaxies, but only 5+/-1% in late-type galaxies. The distribution of position angle offsets is used to test the physical drivers of this difference. Slower dynamical settling time of the gas in elliptical stellar mass distributions accounts for a small increase in misalignment in early-type galaxies. However, gravitational dynamical settling time is insufficient to fully explain the observed differences between early- and late-type galaxies in the distributions of the gas/stellar position angle offsets. LTGs have primarily accreted gas close to aligned rather than settled from misaligned based on analysis of the skewed distribution of PA offsets compared to a dynamical settling model. Local environment density is less important in setting the misalignment fractions than morphology, suggesting that mergers are not the main source of accreted gas in these disks. Cluster environments are found to have gas misalignment driven primarily by cluster processes not by gas accretion.

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.

Recent studies of the diffuse parts of molecular clouds have revealed the presence of parallel, ordered low-density filaments termed striations. Flows along magnetic field lines, Kelvin-Helmholtz instabilities and hydromagnetic waves are amongst the various formation mechanisms proposed. Through a synergy of observational, numerical and theoretical analysis, previous studies singled out the hydromagnetic waves model as the only one that can account for the observed properties of striations. Based on the predictions of that model, we develop here a method for measuring the temporal evolution of striations through a combination of molecular and dust continuum observations. Our method allows us to not only probe temporal variations in molecular clouds but also estimate the strength of both the ordered and fluctuating components of the magnetic field projected on the plane-of-the-sky. We benchmark our new method against chemical and radiative transfer effects through two-dimensional magnetohydrodynamic simulations coupled with non-equilibrium chemical modelling and non-local thermodynamic equilibrium line radiative transfer. We find good agreement between theoretical predictions, simulations and observations of striations in the Taurus molecular cloud. We find a value of $\rm{27 \pm 7} ~\rm{\mu G}$ for the plane-of-sky magnetic field, in agreement with previous estimates via the Davis-Chandrasekhar-Fermi method, and a ratio of fluctuating to ordered component of the magnetic field of $\sim$ 10\%.

HCN(1$-$0) emission traces dense gas and correlates very strongly with star formation rates (SFRs) on scales from small Milky Way clouds to whole galaxies. The observed correlation offers strong constraints on the efficiency of star formation in dense gas, but quantitative interpretation of this constraint requires a mapping from HCN emission to gas mass and density. In this paper we provide the required calibration by postprocessing high-resolution simulations of dense, star-forming clouds to calculate their HCN emission ($L_{\rm HCN}$) and to determine how that emission is related to the underlying gas density distribution and star formation efficiency. We find that HCN emission traces gas with a luminosity-weighted mean number density of $0.8-1.7 \times 10^4\,{\rm cm}^{-3}$ and that HCN luminosity is related to mass of dense gas of $\gtrsim 10^4\,{\rm cm}^{-3}$ with a conversion factor of ${\alpha}_{\rm HCN} \approx 14\,\rm M_{\odot}/(K\,km\,s^{-1}\,{pc}^2)$. We also measure a new empirical relationship between the SFR per global mean freefall time (${\epsilon}_{\rm ff}$) and the SFR$-$HCN relationship, ${\rm SFR}/L_{\rm HCN} = 2.0 \times 10^{-7}\,({\epsilon}_{\rm ff}/0.01)^{1.1}\,\rm M_{\odot}\,{yr}^{-1}/(K\,km\,s^{-1}\,{pc}^2)$. The observed SFR$-$HCN correlation strongly constrains ${\epsilon}_{\rm ff} \approx 1\%$ with a factor of $\sim 3$ systematic uncertainty. The scatter in ${\epsilon}_{\rm ff}$ from cloud to cloud within the Milky Way is a factor of a few. We conclude that $L_{\rm HCN}$ is an effective tracer of dense gas and that the IR$-$HCN correlation is a significant diagnostic of the microphysics of star formation in dense gas.

We present the second major release of data from the SAMI Galaxy Survey. Data Release Two includes data for 1559 galaxies, about 50% of the full survey. Galaxies included have a redshift range 0.004 < z < 0.113 and a large stellar mass range 7.5 < log (M_star/M_sun) < 11.6. The core data for each galaxy consist of two primary spectral cubes covering the blue and red optical wavelength ranges. For each primary cube we also provide three spatially binned spectral cubes and a set of standardised aperture spectra. For each core data product we provide a set of value-added data products. This includes all emission line value-added products from Data Release One, expanded to the larger sample. In addition we include stellar kinematic and stellar population value-added products derived from absorption line measurements. The data are provided online through Australian Astronomical Optics' Data Central. We illustrate the potential of this release by presenting the distribution of ~350,000 stellar velocity dispersion measurements from individual spaxels as a function of R/R_e, divided in four galaxy mass bins. In the highest stellar mass bin (log (M_star/M_sun)>11), the velocity dispersion strongly increases towards the centre, whereas below log (M_star/M_sun)<10 we find no evidence for a clear increase in the central velocity dispersion. This suggests a transition mass around log (M_star/M_sun) ~10 for galaxies with or without a dispersion-dominated bulge.

We present observations of disc-bearing stars in Upper Scorpius (US) and Upper Centaurus-Lupus (UCL) with moderate resolution spectroscopy in order to determine the influence of multiplicity on disc persistence after ~5-20 Myr. Discs were identified using infra-red (IR) excess from the Wide-field Infra-red Survey Explorer (WISE) survey. Our survey consists of 55 US members and 28 UCL members, using spatial and kinematic information to assign a probability of membership. Spectra are gathered from the ANU 2.3m telescope using the Wide Field Spectrograph (WiFeS) to detect radial velocity variations that indicate the presence of a companion. We identify 2 double-lined spectroscopic binaries, both of which have strong IR excess. We find the binary fraction of disc-bearing stars in US and UCL for periods up to 20 years to be $0.06^{0.07}_{0.02}$ and $0.13^{0.06}_{0.03}$ respectively. Based on the multiplicity of field stars, we obtain an expected binary fraction of $0.12^{0.02}_{0.01}$. The determined binary fractions for disc-bearing stars does not vary significantly from the field, suggesting that the overall lifetime of discs may not differ between single and binary star systems.

We identify a counterpart to a Galactic supershell in diffuse radio polarisation, and use this to determine the magnetic fields associated with this object. GSH 006$-$15$+$7 has perturbed the polarised emission at 2.3$\,$GHz, as observed in the S-band Polarisation All Sky Survey (S-PASS), acting as a Faraday screen. We model the Faraday rotation over the shell, and produce a map of Faraday depth over the area across it. Such models require information about the polarised emission behind the screen, which we obtain from the Wilkinson Microwave Anisotropy Probe (WMAP), scaled from 23$\,$GHz to 2.3$\,$GHz, to estimate the synchrotron background behind GSH 006$-$15$+$7. Using the modelled Faraday thickness we determine the magnitude and the plane-of-the-sky structure of the line-of-sight magnetic field in the shell. We find a peak line-of-sight field strength of $|B_\parallel|_\text{peak} = 2.0\substack{+0.01 \\ -0.7}\,\mu$G. Our measurement probes weak magnetic fields in a low-density regime (number densities of $\sim0.4\,$cm$^{-3}$) of the ISM, thus providing crucial information about the magnetic fields in the partially-ionised phase.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

The evolution of self-gravitating clouds of isothermal gas forms the basis of many star formation theories. Therefore it is important to know under what conditions such a cloud will undergo homologous collapse into a single, massive object, or will fragment into a spectrum of smaller ones. And if it fragments, do initial conditions (e.g. Jeans mass, sonic mass) influence the mass function of the fragments, as predicted by many theories of star formation? In this paper we show that the relevant parameter separating homologous collapse from fragmentation is not the Mach number of the initial turbulence (as suspected by many), but the infall Mach number $\mathcal{M}_{\rm infall}\sim\sqrt{G M/(R c_s^2)}$, equivalent to the number of Jeans masses in the initial cloud $N_J$. We also show that fragmenting clouds produce a power-law mass function with slopes close to the expected -2 (i.e. equal mass in all logarithmic mass intervals). However, the low-mass cut-off of this mass function is entirely numerical; the initial properties of the cloud have no effect on it. In other words, if $\mathcal{M}_{\rm infall}\gg 1$, fragmentation proceeds without limit to masses much smaller than the initial Jeans mass.

Integral field spectroscopy surveys provide spatially resolved gas and stellar kinematics of galaxies. They have unveiled a range of atypical kinematic phenomena, which require detailed modelling to understand. We present results from a cosmological simulation that includes stellar and AGN feedback. We find that the distribution of angles between the gas and stellar angular momenta of galaxies is not affected by projection effects. We examine five galaxies ($\approx 6$ per cent of well resolved galaxies) that display atypical kinematics; two of the galaxies have kinematically distinct cores (KDC), while the other three have counter-rotating gas and stars. All five form the majority of their stars in the field, subsequently falling into cosmological filaments where the relative orientation of the stellar angular momentum and the bulk gas flow leads to the formation of a counter-rotating gas disc. The accreted gas exchanges angular momentum with pre-existing co-rotating gas causing it to fall to the centre of the galaxy. This triggers low-level AGN feedback, which reduces star formation. Later, two of the galaxies experience a minor merger (stellar mass ratio $\sim1/10$) with a galaxy on a retrograde orbit compared to the spin of the stellar component of the primary. This produces the KDCs, and is a different mechanism than suggested by other works. The role of minor mergers in the kinematic evolution of galaxies may have been under-appreciated in the past, and large, high-resolution cosmological simulations will be necessary to gain a better understanding in this area.

Using high-resolution (sub-kiloparsec scale) submillimeter data obtained by ALMA, we analyze the star formation rate (SFR), gas content and kinematics in SDP 81, a gravitationally-lensed star-forming galaxy at redshift 3. We estimate the SFR surface density ($\Sigma_{\mathrm{SFR}}$) in the brightest clump of this galaxy to be $357^{+135}_{-85}\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$, over an area of $0.07\pm0.02\,\mathrm{kpc}^2$. Using the intensity-weighted velocity of CO$\,$(5-4), we measure the turbulent velocity dispersion in the plane-of-the-sky and find $\sigma_{\mathrm{v,turb}} = 37\pm5\,\mathrm{km\,s}^{-1}$ for the star-forming clump, in good agreement with previous estimates along the line of sight. Our measurements of gas surface density, freefall time and turbulent Mach number reveal that the role of turbulence is vital to explaining the observed SFR in this clump. While the Kennicutt Schmidt (KS) relation predicts a SFR surface density of $\Sigma_{\mathrm{SFR,KS}} = 52\pm17\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$, the single-freefall model by Krumholz, Dekel and McKee (KDM) predicts $\Sigma_{\mathrm{SFR,KDM}} = 106\pm37\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$. In contrast, the multi-freefall (turbulence) model by Salim, Federrath and Kewley (SFK) gives $\Sigma_{\mathrm{SFR,SFK}} = 491^{+139}_{-194}\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$. Although the SFK relation overestimates the SFR in this clump (possibly due to the ignorance of magnetic field), it provides the best prediction among the available models. Finally, we compare the star formation and gas properties of this high-redshift galaxy to local star-forming regions and find that the SFK relation provides the best estimates of SFR in both local and high-redshift galaxies.

In order to determine the causes of kinematic asymmetry in the H$\alpha$ gas in the SAMI Galaxy Survey sample, we investigate the comparative influences of environment and intrinsic properties of galaxies on perturbation. We use spatially resolved H$\alpha$ velocity fields from the SAMI Galaxy Survey to quantify kinematic asymmetry ($\overline{v_{asym}}$) in nearby galaxies and environmental and stellar mass data from the GAMA survey. {We find that local environment, measured as distance to nearest neighbour, is inversely correlated with kinematic asymmetry for galaxies with $\mathrm{\log(M_*/M_\odot)}>10.0$, but there is no significant correlation for galaxies with $\mathrm{\log(M_*/M_\odot)}<10.0$. Moreover, low mass galaxies ($\mathrm{\log(M_*/M_\odot)}<9.0$) have greater kinematic asymmetry at all separations, suggesting a different physical source of asymmetry is important in low mass galaxies.} We propose that secular effects derived from gas fraction and gas mass may be the primary causes of asymmetry in low mass galaxies. High gas fraction is linked to high $\frac{\sigma_{m}}{V}$ (where $\sigma_m$ is H$\alpha$ velocity dispersion and $V$ the rotation velocity), which is strongly correlated with $\overline{v_{asym}}$, and galaxies with $\log(M_*/M_\odot)<9.0$ have offset $\overline{\frac{\sigma_{m}}{V}}$ from the rest of the sample. Further, asymmetry as a fraction of dispersion decreases for galaxies with $\log(M_*/M_\odot)<9.0$. Gas mass and asymmetry are also inversely correlated in our sample. We propose that low gas masses in dwarf galaxies may lead to asymmetric distribution of gas clouds, leading to increased relative turbulence.

Interstellar gas clouds are often both highly magnetized and supersonically turbulent, with velocity dispersions set by a competition between driving and dissipation. This balance has been studied extensively in the context of gases with constant mean density. However, many astrophysical systems are contracting under the influence of external pressure or gravity, and the balance between driving and dissipation in a contracting, magnetized medium has yet to be studied. In this paper we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of compression in a turbulent, magnetized medium that resembles the physical conditions inside molecular clouds. We find that in some circumstances the combination of compression and magnetic fields leads to a rate of turbulent dissipation far less than that observed in non-magnetized gas, or in non-compressing magnetized gas. As a result, a compressing, magnetized gas reaches an equilibrium velocity dispersion much greater than would be expected for either the hydrodynamic or the non-compressing case. We use the simulation results to construct an analytic model that gives an effective equation of state for a coarse-grained parcel of the gas, in the form of an ideal equation of state with a polytropic index that depends on the dissipation and energy transfer rates between the magnetic and turbulent components. We argue that the reduced dissipation rate and larger equilibrium velocity dispersion produced by compressing, magnetized turbulence has important implications for the driving and maintenance of turbulence in molecular clouds, and for the rates of chemical and radiative processes that are sensitive to shocks and dissipation.

We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{\alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{\beta}, [N II]/H{\alpha}, [S II]/H{\alpha}, and [O I]/H{\alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{\sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions.

We analyse the velocity dispersion properties of 472 z~0.9 star-forming galaxies observed as part of the KMOS Redshift One Spectroscopic Survey (KROSS). The majority of this sample is rotationally dominated (83 +/- 5% with v_C/sigma_0 > 1) but also dynamically hot and highly turbulent. After correcting for beam smearing effects, the median intrinsic velocity dispersion for the final sample is sigma_0 = 43.2 +/- 0.8 km/s with a rotational velocity to dispersion ratio of v_C/sigma_0 = 2.6 +/- 0.1. To explore the relationship between velocity dispersion, stellar mass, star formation rate and redshift we combine KROSS with data from the SAMI survey (z~0.05) and an intermediate redshift MUSE sample (z~0.5). While there is, at most, a weak trend between velocity dispersion and stellar mass, at fixed mass there is a strong increase with redshift. At all redshifts, galaxies appear to follow the same weak trend of increasing velocity dispersion with star formation rate. Our results are consistent with an evolution of galaxy dynamics driven by disks that are more gas rich, and increasingly gravitationally unstable, as a function of increasing redshift. Finally, we test two analytic models that predict turbulence is driven by either gravitational instabilities or stellar feedback. Both provide an adequate description of the data, and further observations are required to rule out either model.

[abridged] We determine the relationship between turbulence energy and gas density variance for 15 molecular clouds in the Solar neighbourhood. We use the linewidths of the CO molecule as the probe of the turbulence energy (sonic Mach number, $\cal{M}_\mathrm{s}$) and three-dimensional models to reconstruct the density probability distribution function ($\rho$-PDF) of the clouds, derived using near-infrared extinction and Herschel dust emission data, as the probe of the density variance ($\sigma_\mathrm{s}$). We find no significant correlation between $\cal{M}_\mathrm{s}$ and $\sigma_\mathrm{s}$ among the studied clouds, however, we also cannot rule out a weak correlation. In the context of turbulence-dominated gas, the range of the $\cal{M}_\mathrm{s}$ and $\sigma_\mathrm{s}$ values corresponds with the model predictions. The data cannot constrain whether or not the turbulence driving parameter, $b$, and/or thermal-to-magnetic pressure ratio, $\beta$, vary among the sample clouds. Most clouds are not in agreement with field strengths stronger than given by $\beta \lesssim 0.05$. A model with $b^2 \beta / (\beta+1) = 0.30 \pm 0.06$ provides an adequate fit to the cloud sample as a whole. When considering the average behaviour of the sample, we can rule out three regimes: (i) strong compression combined with a weak magnetic field ($b \gtrsim 0.7$ and $\beta \gtrsim 3$), (ii) weak compression ($b \lesssim 0.35$), and (iii) strong magnetic field ($\beta \lesssim 0.1$). Including independent magnetic field strength estimates to the analysis, the data rule out solenoidal driving ($b < 0.4$) for the majority of the Solar neighbourhood clouds. However, most clouds have $b$ parameters larger than unity, which indicates a discrepancy with the turbulence-dominated picture; we discuss the possible reasons for this.

We report an integral-field spectroscopic (IFS) observation of a gravitationally lensed spiral galaxy A1689B11 at redshift $z=2.54$. It is the most ancient spiral galaxy discovered to date and the second kinematically confirmed spiral at $z\gtrsim2$. Thanks to gravitational lensing, this is also by far the deepest IFS observation with the highest spatial resolution ($\sim$ 400 pc) on a spiral galaxy at a cosmic time when the Hubble sequence is about to emerge. After correcting for a lensing magnification of 7.2 $\pm$ 0.8, this primitive spiral disk has an intrinsic star formation rate of 22 $\pm$ 2 $M_{\odot}$ yr$^{-1}$, a stellar mass of 10$^{9.8 \pm 0.3}$$M_{\odot}$ and a half-light radius of $r_{1/2}=2.6 \pm 0.7$ kpc, typical of a main-sequence star-forming (SF) galaxy at $z\sim2$. However, the H\alpha\ kinematics show a surprisingly tranquil velocity field with an ordered rotation ($V_{\rm c}$ = 200 $\pm$ 12 km/s) and uniformly small velocity dispersions ($V_{\rm \sigma, mean}$ = 23 $\pm$ 4 km/s and $V_{\rm \sigma, outer-disk}$ = 15 $\pm$ 2 km/s). The low gas velocity dispersion is similar to local spiral galaxies and is consistent with the classic density wave theory where spiral arms form in dynamically cold and thin disks. We speculate that A1689B11 belongs to a population of rare spiral galaxies at $z\gtrsim2$ that mark the formation epoch of thin disks. Future observations with JWST will greatly increase the sample of these rare galaxies and unveil the earliest onset of spiral arms.

We investigate star formation occurring in idealised giant molecular clouds, comparing structures that evolve in isolation versus those undergoing a collision. Two different collision speeds are investigated and the impact of photoionising radiation from the stars is determined. We find that a colliding system leads to more massive star formation both with and without the addition of feedback, raising overall star formation efficiencies (SFE) by a factor of 10 and steepening the high-mass end of the stellar mass function. This rise in SFE is due to increased turbulent compression during the cloud collision. While feedback can both promote and hinder star formation in the isolated system, it increases the SFE by approximately 1.5 times in the colliding case when the thermal speed of the resulting HII regions matches the shock propagation speed in the collision.

We present an analysis of the global stellar populations of galaxies in the SAMI Galaxy Survey. Our sample consists of 1319 galaxies spanning four orders of magnitude in stellar mass and includes all morphologies and environments. We derive luminosity-weighted, single stellar population equivalent stellar ages, metallicities and alpha enhancements from spectra integrated within one effective radius apertures. Variations in galaxy size explain the majority of the scatter in the age--mass and metallicity--mass relations. Stellar populations vary systematically in the plane of galaxy size and stellar mass, such that galaxies with high stellar surface mass density are older, more metal-rich and alpha-enhanced than less dense galaxies. Galaxies with high surface mass densities have a very narrow range of metallicities, however, at fixed mass, the spread in metallicity increases substantially with increasing galaxy size (decreasing density). We identify residual correlations with morphology and environment. At fixed mass and size, galaxies with late-type morphologies, small bulges and low Sersic n are younger than early-type, high n, high bulge-to-total galaxies. Age and metallicity both show small residual correlations with environment; at fixed mass and size, galaxies in denser environments or more massive halos are older and somewhat more metal rich than those in less dense environments. We connect these trends to evolutionary tracks within the size--mass plane.

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central.

Protostellar discs play an important role in star formation, acting as the primary mass reservoir for accretion onto young stars and regulating the extent to which angular momentum and gas is released back into stellar nurseries through the launching of powerful disc winds. In this study, we explore how disc structure relates to the properties of the wind-launching region, mapping out the regions of protostellar discs where wind launching could be viable. We combine a series of 1.5D semi-analytic, steady-state, vertical disc-wind solutions into a radially extended 1+1.5D model, incorporating all three diffusion mechanisms (Ohm, Hall and ambipolar). We observe that the majority of mass outflow via disc winds occurs over a radial width of a fraction of an astronomical unit, with outflow rates attenuating rapidly on either side. We also find that the mass accretion rate, magnetic field strength and surface density profile each have significant effects on both the location of the wind-launching region and the ejection/accretion ratio $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. Increasing either the accretion rate or the magnetic field strength corresponds to a shift of the wind-launching region to smaller radii and a decrease in $\dot{M}_{\rm out}/\dot{M}_{\rm in}$, while increasing the surface density corresponds to launching regions at larger radii with increased $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. Finally, we discover a class of disc winds containing an ineffective launching configuration at intermediate radii, leading to two radially separated regions of wind launching and diminished $\dot{M}_{\rm out}/\dot{M}_{\rm in}$. We find that the wind locations and ejection/accretion ratio are consistent with current observational and theoretical estimates.

We present a set of numerical experiments designed to systematically investigate how turbulence and magnetic fields influence the morphology, energetics, and dynamics of filaments produced in wind-cloud interactions. We cover 3D magnetohydrodynamic systems of supersonic winds impacting clouds with turbulent density, velocity, and magnetic fields. We find that log-normal density distributions aid shock propagation through clouds, increasing their velocity dispersion and producing filaments with expanded cross sections and highly-magnetised knots and sub-filaments. In self-consistently turbulent scenarios the ratio of filament to initial cloud magnetic energy densities is ~1. The effect of Gaussian velocity fields is bound to the turbulence Mach number: Supersonic velocities trigger a rapid cloud expansion; subsonic velocities only have a minor impact. The role of turbulent magnetic fields depends on their tension and is similar to the effect of radiative losses: the stronger the magnetic field or the softer the gas equation of state, the greater the magnetic shielding at wind-filament interfaces and the suppression of Kelvin-Helmholtz instabilities. Overall, we show that including turbulence and magnetic fields is crucial to understanding cold gas entrainment in multi-phase winds. While cloud porosity and supersonic turbulence enhance the acceleration of clouds, magnetic shielding protects them from ablation and causes Rayleigh-Taylor-driven sub-filamentation. Wind-swept clouds in turbulent models reach distances ~15-20 times their core radius and acquire bulk speeds ~0.3-0.4 of the wind speed in one cloud-crushing time, which are three times larger than in non-turbulent models. In all simulations the ratio of turbulent magnetic to kinetic energy densities asymptotes at ~0.1-0.4, and convergence of all relevant dynamical properties requires at least 64 cells per cloud radius.

We investigate the energy sources of random turbulent motions of ionised gas from H$\alpha$ emission in eight local star-forming galaxies from the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. These galaxies satisfy strict pure star-forming selection criteria to avoid contamination from active galactic nuclei (AGN) or strong shocks/outflows. Using the relatively high spatial and spectral resolution of SAMI, we find that -- on sub-kpc scales our galaxies display a flat distribution of ionised gas velocity dispersion as a function of star formation rate (SFR) surface density. A major fraction of our SAMI galaxies shows higher velocity dispersion than predictions by feedback-driven models, especially at the low SFR surface density end. Our results suggest that additional sources beyond star formation feedback contribute to driving random motions of the interstellar medium (ISM) in star-forming galaxies. We speculate that gravity, galactic shear, and/or magnetorotational instability (MRI) may be additional driving sources of turbulence in these galaxies.

We carry out magnetohydrodynamical simulations with FLASH of the formation of a single, a tight binary ($a\sim$2.5 AU) and a wide binary star ($a\sim$45 AU). We study the outflows and jets from these systems to understand the contributions the circumstellar and circumbinary discs have on the efficiency and morphology of the outflow. In the single star and tight binary case we obtain a single pair of jets launched from the system, while in the wide binary case two pairs of jets are observed. This implies that in the tight binary case the contribution of the circumbinary disc on the outflow is greater than that in the wide binary case. We also find that the single star case is the most efficient at transporting mass, linear and angular momentum from the system, while the wide binary case is less efficient ($\sim$50$\%, \sim$33$\%, \sim$42$\%$ of the respective quantities in the single star case). The tight binary's efficiency falls between the other two cases ($\sim$71$\%, \sim$66$\%, \sim$87$\%$ of the respective quantities in the single star case). By studying the magnetic field structure we deduce that the outflows in the single star and tight binary star case are magnetocentrifugally driven, whereas in the wide binary star case the outflows are driven by a magnetic pressure gradient.

Star formation in galaxies relies on the availability of cold, dense gas, which, in turn, relies on factors internal and external to the galaxies. In order to provide a simple model for how star formation is regulated by various physical processes in galaxies, we analyse data at redshift $z=0$ from a hydrodynamical cosmological simulation that includes prescriptions for star formation and stellar evolution, active galactic nuclei (AGN), and their associated feedback processes. This model can determine the star formation rate (SFR) as a function of galaxy stellar mass, gas mass, black hole mass, and environment. We find that gas mass is the most important quantity controlling star formation in low-mass galaxies, and star-forming galaxies in dense environments have higher SFR than their counterparts in the field. In high-mass galaxies, we find that black holes more massive than $\sim10^{7.5}$ M$_\odot$ can be triggered to quench star formation in their host; this mass scale is emergent in our simulations. Furthermore, this black hole mass corresponds to a galaxy bulge mass $\sim2\times10^{10}$ M$_\odot$, consistent with the mass at which galaxies start to become dominated by early types ($\sim3\times10^{10}$ M$_\odot$, as previously shown in observations by Kauffmann et al.). Finally, we demonstrate that our model can reproduce well the SFR measured from observations of galaxies in the GAMA and ALFALFA surveys.

The outskirts of galaxy clusters are characterised by the interplay of gas accretion and dynamical evolution involving turbulence, shocks, magnetic fields and diffuse radio emission. The density and velocity structure of the gas in the outskirts provide an effective pressure support and affect all processes listed above. Therefore it is important to resolve and properly model the turbulent flow in these mildly overdense and relatively large cluster regions; this is a challenging task for hydrodynamical codes. In this work, grid-based simulations of a galaxy cluster are presented. The simulations are performed using adaptive mesh refinement (AMR) based on the regional variability of vorticity, and they include a subgrid scale model (SGS) for unresolved turbulence. The implemented AMR strategy is more effective in resolving the turbulent flow in the cluster outskirts than any previously used criterion based on overdensity. We study a cluster undergoing a major merger, which drives turbulence in the medium. The merger dominates the cluster energy budget out to a few virial radii from the centre. In these regions the shocked intra-cluster medium is resolved and the SGS turbulence is modelled, and compared with diagnostics on larger length scale. The volume-filling factor of the flow with large vorticity is about 60% at low redshift in the cluster outskirts, and thus smaller than in the cluster core. In the framework of modelling radio relics, this point suggests that upstream flow inhomogeneities might affect pre-existing cosmic-ray population and magnetic fields, and the resulting radio emission.

We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the turbulence within each cloud is characterised by the ratio of the density contrast (sigma_rho) to the average Mach number (Mach) within the cloud, b = sigma_rho/Mach. As shown in previous works, b ~ 1/3 indicates solenoidal (divergence-free) driving and b ~ 1 indicates compressive (curl-free) driving. We find that the average b value of all the clouds formed in the simulations has a lower limit of b > 0.2. Importantly, we find that b has a broad distribution, covering values from purely solenoidal to purely compressive driving. Tracking the evolution of individual clouds reveals that the b value for each cloud does not vary significantly over their lifetime. Finally, we perform a resolution study with minimum cell sizes of 8, 4, 2 and 1 pc and find that the average b value increases with increasing resolution. Therefore, we conclude that our measured b values are strictly lower limits and that a resolution better than 1 pc is required for convergence. However, regardless of the resolution, we find that b varies by factors of a few in all cases, which means that the effective driving mode alters significantly from cloud to cloud.

Stars form in cold molecular clouds. However, molecular gas is difficult to observe because the most abundant molecule (H2) lacks a permanent dipole moment. Rotational transitions of CO are often used as a tracer of H2, but CO is much less abundant and the conversion from CO intensity to H2 mass is often highly uncertain. Here we present a new method for estimating the column density of cold molecular gas (Sigma_gas) using optical spectroscopy. We utilise the spatially resolved H-alpha maps of flux and velocity dispersion from the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey. We derive maps of Sigma_gas by inverting the multi-freefall star formation relation, which connects the star formation rate surface density (Sigma_SFR) with Sigma_gas and the turbulent Mach number (Mach). Based on the measured range of Sigma_SFR = 0.005-1.5 M_sol/yr/kpc^2 and Mach = 18-130, we predict Sigma_gas = 7-200 M_sol/pc^2 in the star-forming regions of our sample of 260 SAMI galaxies. These values are close to previously measured Sigma_gas obtained directly with unresolved CO observations of similar galaxies at low redshift. We classify each galaxy in our sample as 'Star-forming' (219) or 'Composite/AGN/Shock' (41), and find that in Composite/AGN/Shock galaxies the average Sigma_SFR, Mach, and Sigma_gas are enhanced by factors of 2.0, 1.6, and 1.3, respectively, compared to Star-forming galaxies. We compare our predictions of Sigma_gas with those obtained by inverting the Kennicutt-Schmidt relation and find that our new method is a factor of two more accurate in predicting Sigma_gas, with an average deviation of 32% from the actual Sigma_gas.

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.

Previous studies have shown that star formation depends on the driving of molecular cloud turbulence, and differences in the driving can produce an order of magnitude difference in the star formation rate. The turbulent driving is characterised by the parameter $\zeta$, with $\zeta=0$ for compressive, curl-free driving (e.g. accretion or supernova explosions), and $\zeta=1$ for solenoidal, divergence-free driving (e.g. Galactic shear). Here we develop a new method to measure $\zeta$ from observations of synchrotron emission from molecular clouds. We calculate statistics of mock synchrotron intensity images produced from magnetohydrodynamic simulations of molecular clouds, in which the driving was controlled to produce different values of $\zeta$. We find that the mean and standard deviation of the log-normalised synchrotron intensity are sensitive to $\zeta$, for values of $\zeta$ between $0$ (curl-free driving) and $0.5$ (naturally-mixed driving). We quantify the dependence of zeta on the direction of the magnetic field relative to the line of sight. We provide best-fit formulae for $\zeta$ in terms of the log-normalised mean and standard deviation of synchrotron intensity, with which $\zeta$ can be determined for molecular clouds that have similar Alfv\'enic Mach number to our simulations. These formulae are independent of the sonic Mach number. Signal-to-noise ratios larger than $5$, and angular resolutions smaller than $5\%$ of the cloud diameter, are required to apply these formulae. Although there are no firm detections of synchrotron emission from molecular clouds, by combining Green Bank Telescope and Very Large Array observations it should be possible to detect synchrotron emission from molecular clouds, thereby constraining the value of $\zeta$.

Magnetic fields play an important role in astrophysical accretion discs, and in the interstellar and intergalactic medium. They drive jets, suppress fragmentation in star-forming clouds and can have a significant impact on the accretion rate of stars. However, the exact amplification mechanisms of cosmic magnetic fields remain relatively poorly understood. Here I start by reviewing recent advances in the numerical and theoretical modelling of the 'turbulent dynamo', which may explain the origin of galactic and inter-galactic magnetic fields. While dynamo action was previously investigated in great detail for incompressible plasmas, I here place particular emphasis on highly compressible astrophysical plasmas, which are characterised by strong density fluctuations and shocks, such as the interstellar medium. I find that dynamo action works not only in subsonic plasmas, but also in highly supersonic, compressible plasmas, as well as for low and high magnetic Prandtl numbers. I further present new numerical simulations from which I determine the growth of the turbulent (un-ordered) magnetic field component ($B_\mathrm{turb}$) in the presence of weak and strong guide fields ($B_0$). I vary $B_0$ over 5 orders of magnitude and find that the dependence of $B_\mathrm{turb}$ on $B_0$ is relatively weak, and can be explained with a simple theoretical model in which the turbulence provides the energy to amplify $B_\mathrm{turb}$. Finally, I discuss some important implications of magnetic fields for the structure of accretion discs, the launching of jets, and the star formation rate of interstellar clouds.

Star formation in the Galactic disc is primarily controlled by gravity, turbulence, and magnetic fields. It is not clear that this also applies to star formation near the Galactic Centre. Here we determine the turbulence and star formation in the CMZ cloud G0.253+0.016. Using maps of 3mm dust emission and HNCO intensity-weighted velocity obtained with ALMA, we measure the volume-density variance $\sigma_{\rho/\rho_0} = 1.3 \pm 0.5$ and turbulent Mach number $\mathcal{M} = 11 \pm 3$. Combining these with turbulence simulations to constrain the plasma $\beta = 0.34 \pm 0.35$, we reconstruct the turbulence driving parameter $b = 0.22 \pm 0.12$ in G0.253+0.016. This low value of $b$ indicates solenoidal (divergence-free) driving of the turbulence in G0.253+0.016. By contrast, typical clouds in the Milky Way disc and spiral arms have a significant compressive (curl-free) driving component ($b > 0.4$). We speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this may reduce the star formation rate by a factor of 7 compared to nearby clouds.

Star formation is primarily controlled by the interplay between gravity, turbulence, and magnetic fields. However, the turbulence and magnetic fields in molecular clouds near the Galactic Center may differ substantially from spiral-arm clouds. Here we determine the physical parameters of the central molecular zone (CMZ) cloud G0.253+0.016, its turbulence, magnetic field and filamentary structure. Using column-density maps based on dust-continuum emission observations with ALMA+Herschel, we identify filaments and show that at least one dense core is located along them. We measure the filament width W_fil=0.17$\pm$0.08pc and the sonic scale {\lambda}_sonic=0.15$\pm$0.11pc of the turbulence, and find W_fil~{\lambda}_sonic. A strong velocity gradient is seen in the HNCO intensity-weighted velocity maps obtained with ALMA+Mopra, which is likely caused by large-scale shearing of G0.253+0.016, producing a wide double-peaked velocity PDF. After subtracting the gradient to isolate the turbulent motions, we find a nearly Gaussian velocity PDF typical for turbulence. We measure the total and turbulent velocity dispersion, 8.8$\pm$0.2km/s and 3.9$\pm$0.1km/s, respectively. Using magnetohydrodynamical simulations, we find that G0.253+0.016's turbulent magnetic field B_turb=130$\pm$50$\mu$G is only ~1/10 of the ordered field component. Combining these measurements, we reconstruct the dominant turbulence driving mode in G0.253+0.016 and find a driving parameter b=0.22$\pm$0.12, indicating solenoidal (divergence-free) driving. We compare this to spiral-arm clouds, which typically have a significant compressive (curl-free) driving component (b>0.4). Motivated by previous reports of strong shearing motions in the CMZ, we speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this reduces the star formation rate (SFR) by a factor of 6.9 compared to typical nearby clouds.

In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks --- fast, intermediate and slow --- distinctly compress and heat up the molecular gas, and so provide an important probe of the physical conditions within a turbulent cloud. Here we introduce the publicly available algorithm, SHOCKFIND, to extract and characterise the mixture of shock families in MHD turbulence. The algorithm is applied to a 3-dimensional simulation of a magnetised turbulent molecular cloud, and we find that both fast and slow MHD shocks are present in the simulation. We give the first prediction of the mixture of turbulence-driven MHD shock families in this molecular cloud, and present their distinct distributions of sonic and Alfvenic Mach numbers. Using subgrid one-dimensional models of MHD shocks we estimate that ~0.03 % of the volume of a typical molecular cloud in the Milky Way will be shock heated above 50 K, at any time during the lifetime of the cloud. We discuss the impact of this shock heating on the dynamical evolution of molecular clouds.

Understanding turbulence is critical for a wide range of terrestrial and astrophysical applications. Here we present first results of the world's highest-resolution simulation of turbulence ever done. The current simulation has a grid resolution of 10048^3 points and was performed on 65536 compute cores on SuperMUC at the Leibniz Supercomputing Centre (LRZ). We present a scaling test of our modified version of the FLASH code, which updates the hydrodynamical equations in less than 3 micro seconds per cell per time step. A first look at the column density structure of the 10048^3 simulation is presented and a detailed analysis is provided in a forthcoming paper.

Star formation is inefficient. Only a few percent of the available gas in molecular clouds forms stars, leading to the observed low star formation rate (SFR). The same holds when averaged over many molecular clouds, such that the SFR of whole galaxies is again surprisingly low. Indeed, considering the low temperatures, molecular clouds should be highly gravitationally unstable and collapse on their global mean freefall timescale. And yet, they are observed to live about 10-100 times longer, i.e., the SFR per freefall time is only a few percent. Thus, other physical mechanisms must provide support against quick global collapse. Magnetic fields, turbulence and stellar feedback have been proposed as stabilising agents controlling star formation, but it is still unclear which of these processes is the most important and what their relative contributions are. Here I present high-resolution, adaptive-mesh-refinement simulations of star cluster formation that include turbulence, magnetic fields, and protostellar jet/outflow feedback. These simulations produce nearly realistic star formation rates consistent with observations, but only if turbulence, magnetic fields and feedback are included simultaneously.

Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve $\nabla \cdot {\bf B}=0$ to machine precision, albeit at significant computational expense.

We perform a comparison between the smoothed particle magnetohydrodynamics (SPMHD) code, Phantom, and the Eulerian grid-based code, Flash, on the small-scale turbulent dynamo in driven, Mach 10 turbulence. We show, for the first time, that the exponential growth and saturation of an initially weak magnetic field via the small-scale dynamo can be successfully reproduced with SPMHD. The two codes agree on the behaviour of the magnetic energy spectra, the saturation level of magnetic energy, and the distribution of magnetic field strengths during the growth and saturation phases. The main difference is that the dynamo growth rate, and its dependence on resolution, differs between the codes, caused by differences in the numerical dissipation and shock capturing schemes leading to differences in the effective Prandtl number in Phantom and Flash.

Probability distribution functions (PDFs) of column densities are an established tool to characterize the evolutionary state of interstellar clouds. Using simulations, we show to what degree their determination is affected by noise, line-of-sight contamination, field selection, and the incomplete sampling in interferometric measurements. We solve the integrals that describe the convolution of a cloud PDF with contaminating sources and study the impact of missing information on the measured column density PDF. The effect of observational noise can be easily estimated and corrected for if the root mean square (rms) of the noise is known. For $\sigma_{noise}$ values below 40% of the typical cloud column density, $N_{peak}$, this involves almost no degradation of the accuracy of the PDF parameters. For higher noise levels and narrow cloud PDFs the width of the PDF becomes increasingly uncertain. A contamination by turbulent foreground or background clouds can be removed as a constant shield if the PDF of the contamination peaks at a lower column or is narrower than that of the observed cloud. Uncertainties in the definition of the cloud boundaries mainly affect the low-column density part of the PDF and the mean density. As long as more than 50% of a cloud are covered, the impact on the PDF parameters is negligible. In contrast, the incomplete sampling of the uv plane in interferometric observations leads to uncorrectable distortions of the PDF of the produced maps. An extension of ALMA's capabilities would allow us to recover the high-column density tail of the PDF but we found no way to measure the intermediate and low column density part of the underlying cloud PDF in interferometric observations.

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

Filaments are ubiquitous in the universe. Recent observations have revealed that stars and star clusters form preferentially along dense filaments. Understanding the formation and properties of filaments is therefore a crucial step in understanding star formation. Here we perform three-dimensional high-resolution magnetohydrodynamical simulations that follow the evolution of molecular clouds and the formation of filaments and stars. We apply a filament detection algorithm and compare simulations with different combinations of physical ingredients: gravity, turbulence, magnetic fields and jet/outflow feedback. We find that gravity-only simulations produce significantly narrower filament profiles than observed, while simulations that include turbulence produce realistic filament properties. For these turbulence simulations, we find a remarkably universal filament width of 0.10 +/- 0.02 pc, which is independent of the star formation history of the clouds. We derive a theoretical model that provides a physical explanation for this characteristic filament width, based on the sonic scale (lambda_sonic) of molecular cloud turbulence. Our derivation provides lambda_sonic as a function of the cloud diameter L, the velocity dispersion sigma_v, the gas sound speed c_s, and the ratio of thermal to magnetic pressure, plasma beta. For typical cloud conditions in the Milky Way spiral arms, we find lambda_sonic = 0.04-0.16 pc, in excellent agreement with the filament width of 0.05-0.15 pc from observations. Consistent with the theoretical model assumptions, we find that the velocity dispersion inside the filaments is subsonic and supersonic outside. We further explain the observed p=2 scaling of the filament density profile, rho ~ r^(-p) with the collision of two planar shocks forming a filament at their intersection.

Filamentary structures are ubiquitous in the interstellar medium, yet their formation, internal structure, and longevity have not been studied in detail. We report the results from a comprehensive numerical study that investigates the characteristics, formation, and evolution of filaments arising from magnetohydrodynamic interactions between supersonic winds and dense clouds. Here we improve on previous simulations by utilising sharper density contrasts and higher numerical resolutions. By following multiple density tracers, we find that material in the envelopes of the clouds is removed and deposited downstream to form filamentary tails, while the cores of the clouds serve as footpoints and late-stage outer layers of these tails. Aspect ratios >12, subsonic velocity dispersions ~0.1-0.3 of the wind sound speed, and magnetic field amplifications ~100 are found to be characteristic of these filaments. We also report the effects of different magnetic field strengths and orientations. The magnetic field strength regulates vorticity production: sinuous filamentary towers arise in non-magnetic environments, while strong magnetic fields inhibit small-scale perturbations at boundary layers making tails less turbulent. Magnetic field components aligned with the direction of the flow favour the formation of pressure-confined flux ropes inside the tails, whilst transverse components tend to form current sheets. Softening the equation of state to nearly isothermal leads to suppression of dynamical instabilities and further collimation of the tail. Towards the final stages of the evolution, we find that small cloudlets and distorted filaments survive the break-up of the clouds and become entrained in the winds, reaching velocities ~0.1 of the wind speed.

The origin of strong magnetic fields in the Universe can be explained by amplifying weak seed fields via turbulent motions on small spatial scales and subsequently transporting the magnetic energy to larger scales. This process is known as the turbulent dynamo and depends on the properties of turbulence, i.e. on the hydrodynamical Reynolds number and the compressibility of the gas, and on the magnetic diffusivity. While we know the growth rate the magnetic energy in the linear regime, the saturation level, i.e. the ratio of magnetic energy to turbulent kinetic energy that can be reached, is not known from analytical calculations. In this paper we present the first scale-dependent saturation model based on an effective turbulent resistivity which is determined by the turnover timescale of turbulent eddies and the magnetic energy density. The magnetic resistivity increases compared to the Spitzer value and the effective scale on which the magnetic energy spectrum is at its maximum moves to larger spatial scales. This process ends when the peak reaches a characteristic wavenumber k* which is determined by the critical magnetic Reynolds number. The saturation level of the dynamo also depends on the type of turbulence and differs for the limits of large and small magnetic Prandtl numbers Pm. With our model we find saturation levels between 43.8% and 1.3% for Pm>>1 and between 2.43% and 0.135% for Pm<<1, where the higher values refer to incompressible turbulence and the lower ones to highly compressible turbulence.

Star formation is inefficient. Only a few percent of the available gas in molecular clouds forms stars, leading to the observed low star formation rate (SFR). The same holds when averaged over many molecular clouds, such that the SFR of whole galaxies is again surprisingly low. Indeed, considering the low temperatures, molecular clouds should be highly gravitationally unstable and collapse on their global mean freefall timescale. And yet, they are observed to live about 10-100 times longer, i.e., the SFR per freefall time (SFR_ff) is only a few percent. Thus, other physical mechanisms must counteract the quick global collapse. Turbulence, magnetic fields and stellar feedback have been proposed as regulating agents, but it is still unclear which of these processes is the most important and what their relative contributions are. Here we run high-resolution simulations including gravity, turbulence, magnetic fields, and jet/outflow feedback. We confirm that clouds collapse on a mean freefall time, if only gravity is considered, producing stars at an unrealistic rate. In contrast, if turbulence, magnetic fields, and feedback are included step-by-step, the SFR is reduced by a factor of 2-3 with each additional physical ingredient. When they all act in concert, we find a constant SFR_ff = 0.04, currently the closest match to observations, but still about a factor of 2-4 higher than the average. A detailed comparison with other simulations and with observations leads us to conclude that only models with turbulence producing large virial parameters, and including magnetic fields and feedback can produce realistic SFRs.

Whilst the star formation rate (SFR) of molecular clouds and galaxies is key in understanding galaxy evolution, the physical processes which determine the SFR remain unclear. This uncertainty about the underlying physics has resulted in various different star formation laws, all having substantial intrinsic scatter. Extending upon previous works that define the column density of star formation (Sigma_SFR) by the gas column density (Sigma_gas), we develop a new universal star formation (SF) law based on the multi-freefall prescription of gas. This new SF law relies predominantly on the probability density function (PDF) and on the sonic Mach number of the turbulence in the star-forming clouds. By doing so we derive a relation where the star formation rate (SFR) correlates with the molecular gas mass per multi-freefall time, whereas previous models had used the average, single-freefall time. We define a new quantity called maximum (multi-freefall) gas consumption rate (MGCR) and show that the actual SFR is only about 0.4% of this maximum possible SFR, confirming the observed low efficiency of star formation. We show that placing observations in this new framework (Sigma_SFR vs. MGCR) yields a significantly improved correlation with 3-4 times reduced scatter compared to previous SF laws and a goodness-of-fit parameter R^2=0.97. By inverting our new relationship, we provide sonic Mach number predictions for kpc-scale observations of Local Group galaxies as well as unresolved observations of local and high-redshift disk and starburst galaxies that do not have independent, reliable estimates for the turbulent cloud Mach number.

The density variance - Mach number relation of the turbulent interstellar medium is relevant for theoretical models of the star formation rate, efficiency, and the initial mass function of stars. Here we use high-resolution hydrodynamical simulations with grid resolutions of up to 1024^3 cells to model compressible turbulence in a regime similar to the observed interstellar medium. We use Fyris Alpha, a shock-capturing code employing a high-order Godunov scheme to track large density variations induced by shocks. We investigate the robustness of the standard relation between the logarithmic density variance (sigma_s^2) and the sonic Mach number (M) of isothermal interstellar turbulence, in the non-isothermal regime. Specifically, we test ideal gases with diatomic molecular (gamma = 7/5) and monatomic (gamma = 5/3) adiabatic indices. A periodic cube of gas is stirred with purely solenoidal forcing at low wavenumbers, leading to a fully-developed turbulent medium. We find that as the gas heats in adiabatic compressions, it evolves along the relationship in the density variance - Mach number plane, but deviates significantly from the standard expression for isothermal gases. Our main result is a new density variance - Mach number relation that takes the adiabatic index into account: sigma_s^2 = ln {1+b^2*M^[(5*gamma+1)/3]} and provides good fits for b*M <= 1. A theoretical model based on the Rankine-Hugoniot shock jump conditions is derived, sigma_s^2 = ln {1+(gamma+1)*b^2*M^2/[(gamma-1)*b^2*M^2+2]}, and provides good fits also for b*M > 1. We conclude that this new relation for adiabatic turbulence may introduce important corrections to the standard relation, if the gas is not isothermal.

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment -- magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures -- models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The subgrid-scale energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ($\beta_\mathrm{p} = 0.25, 1, 2.5, 5, 25$) and sonic Mach numbers ($M_s = 2, 2.5, 4$). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and $\alpha-\beta-\gamma$ closures. We find only mediocre performance of the kinetic eddy-viscosity and $\alpha-\beta-\gamma$ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlation and constant coefficient values over time and and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

The interstellar medium of galaxies is governed by supersonic turbulence, which likely controls the star formation rate (SFR) and the initial mass function (IMF). Interstellar turbulence is non-universal, with a wide range of Mach numbers, magnetic fields strengths, and driving mechanisms. Although some of these parameters were explored, most previous works assumed that the gas is isothermal. However, we know that cold molecular clouds form out of the warm atomic medium, with the gas passing through chemical and thermodynamic phases that are not isothermal. Here we determine the role of temperature variations by modelling non-isothermal turbulence with a polytropic equation of state (EOS), where pressure and temperature are functions of gas density, P~rho^Gamma, T~rho^(Gamma-1). We use grid resolutions of 2048^3 cells and compare polytropic exponents Gamma=0.7 (soft EOS), Gamma=1 (isothermal EOS), and Gamma=5/3 (stiff EOS). We find a complex network of non-isothermal filaments with more small-scale fragmentation occurring for Gamma<1, while Gamma>1 smoothes out density contrasts. The density probability distribution function (PDF) is significantly affected by temperature variations, with a power-law tail developing at low densities for Gamma>1. In contrast, the PDF becomes closer to a lognormal distribution for Gamma<=1. We derive and test a new density variance - Mach number relation that takes Gamma into account. This new relation is relevant for theoretical models of the SFR and IMF, because it determines the dense gas mass fraction of a cloud, from which stars form. We derive the SFR as a function of Gamma and find that it decreases by a factor of ~5 from Gamma=0.7 to Gamma=5/3.

Column-density maps of molecular clouds are one of the most important observables in the context of molecular cloud- and star-formation (SF) studies. With the Herschel satellite it is now possible to determine the column density from dust emission. We use observations and simulations to demonstrate how LOS contamination affects the column density probability distribution function (PDF). We apply a first-order approximation (removing a constant level) to the molecular clouds of Auriga, Maddalena, Carina and NGC3603. In perfect agreement with the simulations, we find that the PDFs become broader, the peak shifts to lower column densities, and the power-law tail of the PDF flattens after correction. All PDFs have a lognormal part for low column densities with a peak at Av~2, a deviation point (DP) from the lognormal at Av(DP)~4-5, and a power-law tail for higher column densities. Assuming a density distribution rho~r^-alpha, the slopes of the power-law tails correspond to alpha(PDF)=1.8, 1.75, and 2.5 for Auriga, Carina, and NGC3603 (alpha~1.5-2 is consistent gravitational collapse). We find that low-mass and high-mass SF clouds display differences in the overall column density structure. Massive clouds assemble more gas in smaller cloud volumes than low-mass SF ones. However, for both cloud types, the transition of the PDF from lognormal shape into power-law tail is found at the same column density (at Av~4-5 mag). Low-mass and high-mass SF clouds then have the same low column density distribution, most likely dominated by supersonic turbulence. At higher column densities, collapse and external pressure can form the power-law tail. The relative importance of the two processes can vary between clouds and thus lead to the observed differences in PDF and column density structure.

The turbulent dynamo may explain the origin of cosmic magnetism. While the exponential amplification of magnetic fields has been studied for incompressible gases, little is known about dynamo action in highly-compressible, supersonic plasmas, such as the interstellar medium of galaxies and the early Universe. Here we perform the first quantitative comparison of theoretical models of the dynamo growth rate and saturation level with three-dimensional magnetohydrodynamical simulations of supersonic turbulence with grid resolutions of up to 1024^3 cells. We obtain numerical convergence and find that dynamo action occurs for both low and high magnetic Prandtl numbers Pm = nu/eta = 0.1-10 (the ratio of viscous to magnetic dissipation), which had so far only been seen for Pm >= 1 in supersonic turbulence. We measure the critical magnetic Reynolds number, Rm_crit = 129 (+43, -31), showing that the compressible dynamo is almost as efficient as in incompressible gas. Considering the physical conditions of the present and early Universe, we conclude that magnetic fields need to be taken into account during structure formation from the early to the present cosmic ages, because they suppress gas fragmentation and drive powerful jets and outflows, both greatly affecting the initial mass function of stars.

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ~ 1000 times lower resolution than would be required without SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about 1/4 of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ~ 1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ~ 3 and may thus be essential for understanding the characteristic mass of the stellar initial mass function.

We introduce a new method for observationally estimating the fraction of momentum density (${\rho}{\mathbf{v}}$) power contained in solenoidal modes (for which $\nabla \cdot {\rho}{\mathbf{v}} = 0$) in molecular clouds. The method is successfully tested with numerical simulations of supersonic turbulence that produce the full range of possible solenoidal/compressible fractions. At present the method assumes statistical isotropy, and does not account for anisotropies caused by (e.g.) magnetic fields. We also introduce a framework for statistically describing density--velocity correlations in turbulent clouds.

The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M_Sun, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H_2 and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of 2-3.

The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (rho-PDF), which determines the star formation rates predicted with analytical models. This function has remained unconstrained by observations. We have developed an approach to quantify rho-PDFs and establish their relation to star formation. The rho-PDFs instigate a density threshold of star formation and allow us to quantify the star formation efficiency above it. The rho-PDFs provide new constraints for star formation theories and correctly predict several key properties of the star-forming interstellar medium.

We use Principal Component Analysis (PCA) to study the gas dynamics in numerical simulations of typical MCs. Our simulations account for the non-isothermal nature of the gas and include a simplified treatment of the time-dependent gas chemistry. We model the CO line emission in a post-processing step using a 3D radiative transfer code. We consider mean number densities n_0 = 30, 100, 300 cm^{-3} that span the range of values typical for MCs in the solar neighbourhood and investigate the slope \alpha_{PCA} of the pseudo structure function computed by PCA for several components: the total density, H2 density, 12CO density, 12CO J = 1 -> 0 intensity and 13CO J = 1 -> 0 intensity. We estimate power-law indices \alpha_{PCA} for different chemical species that range from 0.5 to 0.9, in good agreement with observations, and demonstrate that optical depth effects can influence the PCA. We show that when the PCA succeeds, the combination of chemical inhomogeneity and radiative transfer effects can influence the observed PCA slopes by as much as ~ +/- 0.1. The method can fail if the CO distribution is very intermittent, e.g. in low-density clouds where CO is confined to small fragments.

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

Observations of external galaxies and of local star-forming clouds in the Milky Way have suggested a variety of star formation laws, i.e., simple direct relations between the column density of star formation (Sigma_SFR: the amount of gas forming stars per unit area and time) and the column density of available gas (Sigma_gas). Extending previous studies, we show that these different, sometimes contradictory relations for Milky Way clouds, nearby galaxies, and high-redshift discs and starbursts can be combined in one universal star formation law in which Sigma_SFR is about 1% of the local gas collapse rate, Sigma_gas/t_ff, but a significant scatter remains in this relation. Using computer simulations and theoretical models, we find that the observed scatter may be primarily controlled by physical variations in the Mach number of the turbulence and by differences in the star formation efficiency. Secondary variations can be induced by changes in the virial parameter, turbulent driving and magnetic field. The predictions of our models are testable with observations that constrain both the Mach number and the star formation efficiency in Milky Way clouds, external disc and starburst galaxies at low and high redshift. We also find that reduced telescope resolution does not strongly affect such measurements when Sigma_SFR is plotted against Sigma_gas/t_ff.

We presents results from Smoothed Particle Magnetohydrodynamics simulations of collapsing molecular cloud cores, and dynamo amplification of the magnetic field in the presence of Mach 10 magnetised turbulence. Our star formation simulations have produced, for the first time ever, highly collimated magnetised protostellar jets from the first hydrostatic core phase. Up to 40% of the initial core mass may be ejected through this outflow. The primary difficulty in performing these simulations is maintaining the divergence free constraint of the magnetic field, and to address this issue, we have developed a new divergence cleaning method which has allowed us to stably follow the evolution of these protostellar jets for long periods. The simulations performed of supersonic MHD turbulence are able to exponentially amplify magnetic energy by up to 10 orders of magnitude via turbulent dynamo. To reduce numerical dissipation, a new shock detection algorithm is utilised which is able to track magnetic shocks throughout a large range of magnetic field strengths.

Molecular clouds are supersonically turbulent. This turbulence governs the initial mass function and the star formation rate. In order to understand the details of star formation, it is therefore essential to understand the properties of turbulence, in particular the probability distribution of density in turbulent clouds. We present formaldehyde volume density measurements of a non-star-forming cloud along the line of sight towards W49A. We use these measurements in conjunction with total mass estimates from 13CO to infer the shape of the density probability distribution function. This method is complementary to measurements of turbulence via the column density distribution and should be applicable to any molecular cloud with detected CO. We show that turbulence in this cloud is probably compressively driven, with a compressive-to-total Mach number ratio $b = \mathcal{M}_C/\mathcal{M}>0.4$. We measure the standard deviation of the density distribution, constraining it to the range $1.5 < \sigma_s < 1.9$ assuming that the density is lognormally distributed. This measurement represents an essential input into star formation laws. The method of averaging over different excitation conditions to produce a model of emission from a turbulent cloud is generally applicable to optically thin line observations.

Compressible turbulence shapes the structure of the interstellar medium of our Galaxy and likely plays an important role also during structure formation in the early Universe. The density PDF and the power spectrum of such compressible, supersonic turbulence are the key ingredients for theories of star formation. However, both the PDF and the spectrum are still a matter of debate, because theoretical predictions are limited and simulations of supersonic turbulence require enormous resolutions to capture the inertial-range scaling. To advance our limited knowledge of compressible turbulence, we here present and analyse the world's largest simulations of supersonic turbulence. We compare hydrodynamic models with numerical resolutions of 256^3-4096^3 mesh points and with two distinct driving mechanisms, solenoidal (divergence-free) driving and compressive (curl-free) driving. We find convergence of the density PDF, with compressive driving exhibiting a much wider and more intermittent density distribution than solenoidal driving. Analysing the power spectrum of the turbulence, we find a pure velocity scaling close to Burgers turbulence with P(v) k^(-2) for both driving modes in our hydrodynamical simulations with Mach = 17. The spectrum of the density-weighted velocity rho^(1/3)v, however, does not provide the previously suggested universal scaling for supersonic turbulence. We find that the power spectrum P(rho^(1/3)v) scales with wavenumber as k^(-1.74) for solenoidal driving, close to incompressible Kolmogorov turbulence, k^(-5/3), but is significantly steeper with k^(-2.10) for compressive driving. We show that this is consistent with a recent theoretical model for compressible turbulence that predicts P(rho^(1/3)v) k^(-19/9) in the presence of a strong div(v) component as is produced by compressive driving and remains remarkably constant throughout the supersonic turbulent cascade.

We examine the physical parameters that affect the accumulation of gas in molecular clouds to high column densities where the formation of stars takes place. In particular, we analyze the dense gas mass fraction (DGMF) in a set of self-gravitating, isothermal, magnetohydrodynamic turbulence simulations including sink particles to model star formation. We find that the simulations predict close to exponential DGMFs over the column density range N(H2) = 3-25 x 10^{21} cm^{-2} that can be easily probed via, e.g., dust extinction measurements. The exponential slopes correlate with the type of turbulence driving and also with the star formation efficiency. They are almost uncorrelated with the sonic Mach number and magnetic-field strength. The slopes at early stages of cloud evolution are steeper than at the later stages. A comparison of these predictions with observations shows that only simulations with relatively non-compressive driving (b ~< 0.4) agree with the DGMFs of nearby molecular clouds. Massive infrared dark clouds can show DGMFs that are in agreement with more compressive driving. The DGMFs of molecular clouds can be significantly affected by how compressive the turbulence is on average. Variations in the level of compression can cause scatter to the DGMF slopes, and some variation is indeed necessary to explain the spread of the observed DGMF slopes. The observed DGMF slopes can also be affected by the clouds' star formation activities and statistical cloud-to-cloud variations.

A key parameter to the description of all star formation processes is the density structure of the gas. In this letter, we make use of probability distribution functions (PDFs) of Herschel column density maps of Orion B, Aquila, and Polaris, obtained with the Herschel Gould Belt survey (HGBS). We aim to understand which physical processes influence the PDF shape, and with which signatures. The PDFs of Orion B (Aquila) show a lognormal distribution for low column densities until Av 3 (6), and a power-law tail for high column densities, consistent with a rho r^-2 profile for the equivalent spherical density distribution. The PDF of Orion B is broadened by external compression due to the nearby OB stellar aggregates. The PDF of a quiescent subregion of the non-star-forming Polaris cloud is nearly lognormal, indicating that supersonic turbulence governs the density distribution. But we also observe a deviation from the lognormal shape at Av>1 for a subregion in Polaris that includes a prominent filament. We conclude that (i) the point where the PDF deviates from the lognormal form does not trace a universal Av-threshold for star formation, (ii) statistical density fluctuations, intermittency and magnetic fields can cause excess from the lognormal PDF at an early cloud formation stage, (iii) core formation and/or global collapse of filaments and a non-isothermal gas distribution lead to a power-law tail, and (iv) external compression broadens the column density PDF, consistent with numerical simulations.

(Abridged) The small-scale dynamo may play a substantial role in magnetizing the Universe under a large range of conditions, including subsonic turbulence at low Mach numbers, highly supersonic turbulence at high Mach numbers and a large range of magnetic Prandtl numbers Pm, i.e. the ratio of kinetic viscosity to magnetic resistivity. Low Mach numbers may in particular lead to the well-known, incompressible Kolmogorov turbulence, while for high Mach numbers, we are in the highly compressible regime, thus close to Burgers turbulence. In this study, we explore whether in this large range of conditions, a universal behavior can be expected. Our starting point are previous investigations in the kinematic regime. Here, analytic studies based on the Kazantsev model have shown that the behavior of the dynamo depends significantly on Pm and the type of turbulence, and numerical simulations indicate a strong dependence of the growth rate on the Mach number of the flow. Once the magnetic field saturates on the current amplification scale, backreactions occur and the growth is shifted to the next-larger scale. We employ a Fokker-Planck model to calculate the magnetic field amplification during the non-linear regime, and find a resulting power-law growth that depends on the type of turbulence invoked. For Kolmogorov turbulence, we confirm previous results suggesting a linear growth of magnetic energy. For more general turbulent spectra, where the turbulent velocity v_t scales with the characteristic length scale as u_\ell\propto \ell^{\vartheta}, we find that the magnetic energy grows as (t/T_{ed})^{2\vartheta/(1-\vartheta)}, with t the time-coordinate and T_{ed} the eddy-turnover time on the forcing scale of turbulence. For Burgers turbulence, \vartheta=1/2, a quadratic rather than linear growth may thus be expected, and a larger timescale until saturation is reached.

We study the star formation efficiency (SFE) in simulations and observations of turbulent, magnetized, molecular clouds. We find that the probability density functions (PDFs) of the density and the column density in our simulations with solenoidal, mixed, and compressive forcing of the turbulence, sonic Mach numbers of 3-50, and magnetic fields in the super- to the trans-Alfvenic regime, all develop power-law tails of flattening slope with increasing SFE. The high-density tails of the PDFs are consistent with equivalent radial density profiles, rho ~ r^(-kappa) with kappa ~ 1.5-2.5, in agreement with observations. Studying velocity-size scalings, we find that all the simulations are consistent with the observed v ~ l^(1/2) scaling of supersonic turbulence, and seem to approach Kolmogorov turbulence with v ~ l^(1/3) below the sonic scale. The velocity-size scaling is, however, largely independent of the SFE. In contrast, the density-size and column density-size scalings are highly sensitive to star formation. We find that the power-law slope alpha of the density power spectrum, P(rho,k) ~ k^alpha, or equivalently the Delta-variance spectrum of column density, DV(Sigma,l) ~ l^(-alpha), switches sign from alpha < 0 for SFE ~ 0 to alpha > 0 when star formation proceeds (SFE > 0). We provide a relation to compute the SFE from a measurement of alpha. Studying the literature, we find values ranging from alpha = -1.6 to +1.6 in observations covering scales from the large-scale atomic medium, over cold molecular clouds, down to dense star-forming cores. From those alpha values, we infer SFEs and find good agreement with independent measurements based on young stellar object (YSO) counts, where available. Our SFE-alpha relation provides an independent estimate of the SFE based on the column density map of a cloud alone, without requiring a priori knowledge of star-formation activity or YSO counts.

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}.

We derive and compare six theoretical models for the star formation rate (SFR) - the Krumholz & McKee (KM), Padoan & Nordlund (PN), and Hennebelle & Chabrier (HC) models, and three multi-freefall versions of these, suggested by HC - all based on integrals over the log-normal distribution of turbulent gas. We extend all theories to include magnetic fields, and show that the SFR depends on four basic parameters: (1) virial parameter alpha_vir; (2) sonic Mach number M; (3) turbulent forcing parameter b, which is a measure for the fraction of energy driven in compressive modes; and (4) plasma beta=2(M_A/M)^2 with the Alfven Mach number M_A. We compare all six theories with MHD simulations, covering cloud masses of 300 to 4x10^6 solar masses and Mach numbers M = 3 to 50 and M_A = 1 to infinity, with solenoidal (b=1/3), mixed (b=0.4) and compressive turbulent (b=1) forcings. We find that the SFR increases by a factor of four between M=5 and 50 for compressive forcing and alpha_vir~1. Comparing forcing parameters, we see that the SFR is more than 10x higher with compressive than solenoidal forcing for Mach 10 simulations. The SFR and fragmentation are both reduced by a factor of two in strongly magnetized, trans-Alfvenic turbulence compared to hydrodynamic turbulence. All simulations are fit simultaneously by the multi-freefall KM and multi-freefall PN theories within a factor of two over two orders of magnitude in SFR. The simulated SFRs cover the range and correlation of SFR column density with gas column density observed in Galactic clouds, and agree well for star formation efficiencies SFE = 1% to 10% and local efficiencies epsilon = 0.3 to 0.7 due to feedback. We conclude that the SFR is primarily controlled by interstellar turbulence, with a secondary effect coming from magnetic fields.

The first galaxies form due to gravitational collapse of primordial halos. During this collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev theory, which describes the small-scale dynamo analytically, to study magnetic field amplification for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic Reynolds number, Re^(1/2). In the case of highly compressible turbulence (Burgers turbulence) the growth rate increases proportional to Re^(1/3). With a detailed chemical network we are able to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the growth rate of the small-scale dynamo quantitatively and predict the evolution of the small-scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy-timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant field structures can thus be expected already during the formation of the first objects in the Universe.

Stars form by the gravitational collapse of interstellar gas. The thermodynamic response of the gas can be characterized by an effective equation of state. It determines how gas heats up or cools as it gets compressed, and hence plays a key role in regulating the process of stellar birth on virtually all scales, ranging from individual star clusters up to the galaxy as a whole. We present a systematic study of the impact of thermodynamics on gravitational collapse in the context of high-redshift star formation, but argue that our findings are also relevant for present-day star formation in molecular clouds. We consider a polytropic equation of state, P = k rho^Gamma, with both sub-isothermal exponents Gamma < 1 and super-isothermal exponents Gamma > 1. We find significant differences between these two cases. For Gamma > 1, pressure gradients slow down the contraction and lead to the formation of a virialized, turbulent core. Weak magnetic fields are strongly tangled and efficiently amplified via the small-scale turbulent dynamo on timescales corresponding to the eddy-turnover time at the viscous scale. For Gamma < 1, on the other hand, pressure support is not sufficient for the formation of such a core. Gravitational contraction proceeds much more rapidly and the flow develops very strong shocks, creating a network of intersecting sheets and extended filaments. The resulting magnetic field lines are very coherent and exhibit a considerable degree of order. Nevertheless, even under these conditions we still find exponential growth of the magnetic energy density in the kinematic regime.

We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.

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 study the amplification of magnetic fields during the formation of primordial halos. The turbulence generated by gravitational infall motions during the formation of the first stars and galaxies can amplify magnetic fields very efficiently and on short timescales up to dynamically significant values. Using the Kazantsev theory, which describes the so-called small-scale dynamo - a magnetohydrodynamical process converting kinetic energy from turbulence into magnetic energy - we can then calculate the growth rate of the small-scale magnetic field. Our calculations are based on a detailed chemical network and we include non-ideal magnetohydrodynamical effects such as ambipolar diffusion and Ohmic dissipation. We follow the evolution of the magnetic field up to larger scales until saturation occurs on the Jeans scale. Assuming a weak magnetic seed field generated by the Biermann battery process, both Burgers and Kolmogorov turbulence lead to saturation within a rather small density range. Such fields are likely to become relevant after the formation of a protostellar disk and, thus, could influence the formation of the first stars and galaxies in the Universe.

It is widely accepted that supersonic, magnetised turbulence plays a fundamental role for star formation in molecular clouds. It produces the initial dense gas seeds out of which new stars can form. However, the exact relation between gas compression, turbulent Mach number, and magnetic field strength is still poorly understood. Here, we introduce and test an analytical prediction for the relation between the density variance and the root-mean-square Mach number in supersonic, isothermal, magnetised turbulent flows. We approximate the density and velocity structure of the interstellar medium as a superposition of shock waves. We obtain the density contrast considering the momentum equation for a single magnetised shock and extrapolate this result to the entire cloud. Depending on the field geometry, we then make three different assumptions based on observational and theoretical constraints: B independent of density, B proportional to the root square of the density and B proportional to the density. We test the analytically derived density variance--Mach number relation with numerical simulations, and find that for B proportional to the root square of the density, the variance in the logarithmic density contrast, $\sigma_{\ln \rho/\rho_0}^2=\ln[1+b^2\mathscr{M}^2\beta_0/(\beta_0+1)]$, fits very well to simulated data with turbulent forcing parameter b=0.4, when the gas is super-Alfv\'enic. However, this result breaks down when the turbulence becomes trans-Alfv\'enic or sub-Alfv\'enic, because in this regime the turbulence becomes highly anisotropic. Our density variance--Mach number relations simplify to the purely hydrodynamic relation as the ratio of thermal to magnetic pressure $\beta_0$ approaches infinite.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

Here we summarize our recent results of high-resolution computer simulations on the turbulent amplification of weak magnetic seed fields showing that such fields will be exponentially amplified also during the gravitational collapse reminiscent to the situation during primordial star formation. The exponential magnetic field amplification is driven by the turbulent small-scale dynamo that can be only observed in computer simulations if the turbulent motions in the central core are sufficiently resolved. We find that the Jeans length, which determines the central core region, has to be resolved by at least 30 grid cells to capture the dynamo activity. We conclude from our studies that strong magnetic fields will be unavoidably created already during the formation of the first stars in the Universe, potentially influencing their evolution and mass distribution.

The small-scale dynamo is a process by which turbulent kinetic energy is converted into magnetic energy, and thus is expected to depend crucially on the nature of turbulence. In this work, we present a model for the small-scale dynamo that takes into account the slope of the turbulent velocity spectrum v(l) ~ l^theta, where l and v(l) are the size of a turbulent fluctuation and the typical velocity on that scale. The time evolution of the fluctuation component of the magnetic field, i.e., the small-scale field, is described by the Kazantsev equation. We solve this linear differential equation for its eigenvalues with the quantum-mechanical WKB-approximation. The validity of this method is estimated as a function of the magnetic Prandtl number Pm. We calculate the minimal magnetic Reynolds number for dynamo action, Rm_crit, using our model of the turbulent velocity correlation function. For Kolmogorov turbulence (theta=1/3), we find that the critical magnetic Reynolds number is approximately 110 and for Burgers turbulence (theta=1/2) approximately 2700. Furthermore, we derive that the growth rate of the small-scale magnetic field for a general type of turbulence is Gamma ~ Re^((1-theta)/(1+theta)) in the limit of infinite magnetic Prandtl numbers. For decreasing magnetic Prandtl number (down to Pm approximately larger than 10), the growth rate of the small-scale dynamo decreases. The details of this drop depend on the WKB-approximation, which becomes invalid for a magnetic Prandtl number of about unity.

We study the influence of initial conditions on the magnetic field amplification during the collapse of a magnetised gas cloud. We focus on the dependence of the growth and saturation level of the dynamo generated field on the turbulent properties of the collapsing cloud. In particular, we explore the effect of varying the initial strength and injection scale of turbulence and the initial uniform rotation of the collapsing magnetised cloud. In order to follow the evolution of the magnetic field in both the kinematic and the nonlinear regime, we choose an initial field strength of $\simeq 1\,\mkG$ with the magnetic to kinetic energy ratio, $E_{\rm m}/E_{\rm k} \sim 10^{-4}$. Both gravitational compression and the small-scale dynamo initially amplify the magnetic field. Further into the evolution, the dynamo-generated magnetic field saturates but the total magnetic field continues to grow because of compression. The saturation of the small-scale dynamo is marked by a change in the slope of $B/\rho^{2/3}$ and by a shift in the peak of the magnetic energy spectrum from small scales to larger scales. For the range of initial Mach numbers explored in this study, the dynamo growth rate increases as the Mach number increases from $v_{\rm rms}/c_{\rm s}\sim 0.2$ to 0.4 and then starts decreasing from $v_{\rm rms}/c_{\rm s}\sim 1.0$. We obtain saturation values of $E_{\rm m}/E_{\rm k} = 0.2 - 0.3$ for these runs. Simulations with different initial injection scales of turbulence also show saturation at similar levels. For runs with different initial rotation of the cloud, the magnetic energy saturates at $E_{\rm m}/E_{\rm k}\sim 0.2 - 0.4$ of the equipartition value. (Abridged)

We present results from high-resolution three-dimensional simulations of the turbulent interstellar medium that study the influence of the nature of the turbulence on the formation of molecular hydrogen. We have examined both solenoidal (divergence-free) and compressive (curl-free) turbulent driving, and show that compressive driving leads to faster H2 formation, owing to the higher peak densities produced in the gas. The difference in the H2 formation rate can be as much as an order of magnitude at early times, but declines at later times as the highest density regions become fully molecular and stop contributing to the total H2 formation rate. We have also used our results to test a simple prescription suggested by Gnedin et al. (2009) for modeling the influence of unresolved density fluctuations on the H2 formation rate in large-scale simulations of the ISM. We find that this approach works well when the H2 fraction is small, but breaks down once the highest density gas becomes fully molecular.

We study the mass-to-flux ratio (M/\Phi) of clumps and cores in simulations of supersonic, magnetohydrodynamical turbulence for different initial magnetic field strengths. We investigate whether the (M/\Phi)-ratio of core and envelope, R = (M/\Phi)_{core}/(M/\Phi)_{envelope} can be used to distinguish between theories of ambipolar diffusion and turbulence-regulated star formation. We analyse R for different Lines-of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average and median values of |R| for different times and initial magnetic field strengths are typically greater, but close to unity, 2) the average and median values of |R| saturate at average values of |R| ~ 1 for smaller magnetic fields, 3) values of |R| < 1 for small magnetic fields in the envelope are caused by field reversals when turbulence twists the field lines such that field components in different directions average out. Finally, we propose two mechanisms for generating values |R| ~< 1 for the weak and strong magnetic field limit in the context of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo leads to a significantly increased flux in the core and we find |R| ~< 1. Second, in the strong field limit, field reversals in the envelope also lead to values |R| ~< 1. These reversals are less likely to occur in the core region where the velocity field is more coherent and the internal velocity dispersion is typically subsonic.

We investigate the formation of protostellar clusters during the collapse of dense molecular cloud cores with a focus on the evolution of potential and kinetic energy, the degree of substructure, and the early phase of mass segregation. Our study is based on a series of hydrodynamic simulations of dense cores, where we vary the initial density profile and the initial turbulent velocity. In the three-dimensional adaptive mesh refinement simulations, we follow the dynamical formation of filaments and protostars until a star formation efficiency of 20%. Despite the different initial configurations, the global ensemble of all protostars in a setup shows a similar energy evolution and forms sub-virial clusters with an energy ratio $E_\mathrm{kin}/|E_\mathrm{pot}|\sim0.2$. Concentrating on the innermost central region, the clusters show a roughly virialised energy balance. However, the region of virial balance only covers the innermost $\sim10-30%$ of all the protostars. In all simulations with multiple protostars, the total kinetic energy of the protostars is higher than the kinetic energy of the gas cloud, although the protostars only contain 20% of the total mass. The clusters vary significantly in size, mass, and number of protostars, and show different degrees of substructure and mass segregation. Flat density profiles and compressive turbulent modes produce more subclusters then centrally concentrated profiles and solenoidal turbulence. We find that dynamical relaxation and hence dynamical mass segregation is very efficient in all cases from the very beginning of the nascent cluster, i.e., during a phase when protostars are constantly forming and accreting.

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.

We investigate the impact of different initial conditions for the initial density profile and the initial turbulence on the formation process of protostellar clusters. We study the collapse of dense molecula r cloud cores with three-dimensional adaptive mesh refinement simulations. We focus our discussion on the distribution of the gas among the protostellar objects in the turbulent dynamical cores. Despite the large variations in the initial configurations and the resulting gas and cluster morphology we find that all stellar clusters follow a very similar gas accretion behaviour. Once secondary protostars begin to form, the central region of a cluster is efficiently shielded from further accretion. Hence, objects located close to the centre are starved of material, as indicated by a strong decrease of the central acc retion rate. This Fragmentation Induced Starvation occurs not only in rotationally supported discs and filaments, but also in more spherically symmetric clusters with complex chaotic motions.

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

We discuss the amplification of magnetic fields by the small-scale dynamo, a process that could efficiently produce strong magnetic fields in the first galaxies. In addition, we derive constraints on the primordial field strength from the epoch of reionization.

We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity.

Turbulence plays a major role in the formation and evolution of molecular clouds. The problem is that turbulent velocities are convolved with the density of an observed region. To correct for this convolution, we investigate the relation between the turbulence spectrum of model clouds, and the statistics of their synthetic observations obtained from Principal Component Analysis (PCA). We apply PCA to spectral maps generated from simulated density and velocity fields, obtained from hydrodynamic simulations of supersonic turbulence, and from fractional Brownian motion fields with varying velocity, density spectra, and density dispersion. We examine the dependence of the slope of the PCA structure function, alpha_PCA, on intermittency, on the turbulence velocity (beta_v) and density (beta_n) spectral indexes, and on density dispersion. We find that PCA is insensitive to beta_n and to the log-density dispersion sigma_s, provided sigma_s < 2. For sigma_s>2, alpha_PCA increases with sigma_s due to the intermittent sampling of the velocity field by the density field. The PCA calibration also depends on intermittency. We derive a PCA calibration based on fBms with sigma_s<2 and apply it to 367 CO spectral maps of molecular clouds in the Galactic Ring Survey. The average slope of the PCA structure function, <alpha_PCA>=0.62\pm0.2, is consistent with the hydrodynamic simulations and leads to a turbulence velocity exponent <beta_v>=2.06\pm0.6 for a non-intermittent, low density dispersion flow. Accounting for intermittency and density dispersion, the coincidence between the PCA slope of the GRS clouds and the hydrodynamic simulations suggests beta_v~1.9, consistent with both Burgers and compressible intermittent turbulence.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme's robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma beta. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma beta. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512^3 grid cells we identify a narrow inertial range with the expected power-law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate twice as high from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A_V^tail = 2-4 mag, or equivalently, at \Sigma^tail = 40-80 Ms pc^{-2}. The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n = 10^3 cm^{-3}. This is 5-10 times larger than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A_V^tail may be linked to a transition between this population and the surrounding medium. The star formation rates in the clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A_V^tail, dropping abruptly below that threshold. These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A_V^tail.

We present a detailed parameter study of collapsing turbulent cloud cores, varying the initial density profile and the initial turbulent velocity field. We systematically investigate the influence of different initial conditions on the star formation process, mainly focusing on the fragmentation, the number of formed stars, and the resulting mass distributions. Our study compares four different density profiles (uniform, Bonnor-Ebert type, $\rho\propto r^{-1.5}$, and $\rho\propto r^{-2}$), combined with six different supersonic turbulent velocity fields (compressive, mixed, and solenoidal, initialised with two different random seeds each) in three-dimensional simulations using the adaptive-mesh refinement, hydrodynamics code FLASH. The simulations show that density profiles with flat cores produce hundreds of low-mass stars, either distributed throughout the entire cloud or found in subclusters, depending on the initial turbulence. Concentrated density profiles always lead to the formation of one high-mass star in the centre of the cloud and, if at all, low-mass stars surrounding the central one. In uniform and Bonnor-Ebert type density distributions, compressive initial turbulence leads to local collapse about 25% earlier than solenoidal turbulence. However, central collapse in the steep power-law profiles is too fast for the turbulence to have any significant influence. We conclude that (I) the initial density profile and turbulence mainly determine the cloud evolution and the formation of clusters, (II) the initial mass function (IMF) is not universal for all setups, and (III) that massive stars are much less likely to form in flat density distributions. The IMFs obtained in the uniform and Bonnor-Ebert type density profiles are more consistent with the observed IMF, but shifted to lower masses.

Recent FERMI observations provide a lower limit of 10^{-15} G for the magnetic field strength in the intergalactic medium (IGM). This is consistent with theoretical expectations based on the Biermann battery effect, which predicts such IGM fields already at redshifts z~10. During gravitational collapse, such magnetic fields can be amplified by compression and by turbulence, giving rise to the small-scale dynamo. On scales below the Jeans length, the eddy turnover timescale is much shorter than the free-fall timescale, so that saturation can be reached during collapse. This scenario has been tested and confirmed with magneto-hydrodynamical simulations following the collapse of a turbulent, weakly magnetized cloud. Based on a spectral analysis, we confirm that turbulence is injected on the Jeans scale. For the power spectrum of the magnetic field, we obtain the Kazantsev slope which is characteristic for the small-scale dynamo. A calculation of the critical length scales for ambipolar diffusion and Ohmic dissipation shows that these scales are always small enough to allow significant amplification of the magnetic field by small-scale eddies. We discuss potential implications for the protostellar accretion disk, with particular focus on the magneto-rotational instability, which may change the morphology of the disk and reduce the accretion rate by a factor of a few.

Cosmic structure formation is characterized by the complex interplay between gravity, turbulence, and magnetic fields. The processes by which gravitational energy is converted into turbulent and magnetic energies, however, remain poorly understood. Here, we show with high-resolution, adaptive-mesh simulations that MHD turbulence is efficiently driven by extracting energy from the gravitational potential during the collapse of a dense gas cloud. Compressible motions generated during the contraction are converted into solenoidal, turbulent motions, leading to a natural energy ratio of E_sol/E_tot of approximately 2/3. We find that the energy injection scale of gravity-driven turbulence is close to the local Jeans scale. If small seeds of the magnetic field are present, they are amplified exponentially fast via the small-scale dynamo process. The magnetic field grows most efficiently on the smallest scales, for which the stretching, twisting, and folding of field lines, and the turbulent vortices are sufficiently resolved. We find that this scale corresponds to about 30 grid cells in the simulations. We thus suggest a new minimum resolution criterion of 30 cells per Jeans length in (magneto)hydrodynamical simulations of self-gravitating gas, in order to resolve turbulence on the Jeans scale, and to capture minimum dynamo amplification of the magnetic field. Due to numerical diffusion, however, any existing simulation today can at best provide lower limits on the physical growth rates. We conclude that a small, initial magnetic field can grow to dynamically important strength on time scales significantly shorter than the free-fall time of the cloud.

We aim to better understand how the spatial structure of molecular clouds is governed by turbulence. For that, we study the large-scale spatial distribution of low density molecular gas and search for characteristic length scales. We employ a 35 square degrees 13CO 1-0 molecular line survey of Cygnus X and visual extinction (A_V) maps of 17 Galactic clouds to analyse the spatial structure using the Delta-variance method. This sample contains a large variety of different molecular cloud types with different star forming activity. The Delta-variance spectra obtained from the A_V maps show differences between low-mass star-forming (SF) clouds and massive giant molecular clouds (GMC) in terms of shape of the spectrum and its power-law exponent beta. Low-mass SF clouds have a double-peak structure with characteristic size scales around 1 pc (though with a large scatter around this value) and 4 pc. GMCs show no characteristic scale in the A_V-maps, which can partly be ascribed to a distance effect due to a larger line-of-sight (LOS) confusion. The Delta-variance for Cygnus, determined from the 13CO survey, shows characteristic scales at 4 pc and 40 pc, either reflecting the filament structure and large-scale turbulence forcing or - for the 4 pc scale - the scale below which the 13CO 1-0 line becomes optically thick. Though there are different processes that can introduce characteristic scales, i.e. geometry, decaying turbulence the transition scale from supersonic to subsonic turbulence (the sonic scale), line-of-sight effects and energy injection due to expanding supernova shells, outflows, HII-regions, and the relative contribution of these effects strongly varies from cloud to cloud, it is remarkable that the resulting turbulent structure of molecular clouds shows similar characteristics.

We examine the relation between the density variance and the mean-square Mach number in supersonic, isothermal turbulence, assumed in several recent analytic models of the star formation process. From a series of calculations of supersonic, hydrodynamic turbulence driven using purely solenoidal Fourier modes, we find that the `standard' relationship between the variance in the log of density and the Mach number squared, i.e., sigma^2_(ln rho/rhobar)=ln (1+b^2 M^2), with b = 1/3 is a good fit to the numerical results in the supersonic regime up to at least Mach 20, similar to previous determinations at lower Mach numbers. While direct measurements of the variance in linear density are found to be severely underestimated by finite resolution effects, it is possible to infer the linear density variance via the assumption of log-normality in the Probability Distribution Function. The inferred relationship with Mach number, consistent with sigma_(rho/rhobar) ~ b M with b=1/3, is, however, significantly shallower than observational determinations of the relationship in the Taurus Molecular Cloud and IC5146 (both consistent with b~ 0.5), implying that additional physics such as gravity is important in these clouds and/or that turbulent driving in the ISM contains a significant compressive component. Magnetic fields are not found to change this picture significantly, in general reducing the measured variances and thus worsening the discrepancy with observations.

We formulate and implement the Euler equations with SGS dynamics and provide numerical tests of an SGS turbulence energy model that predicts the turbulent pressure of unresolved velocity fluctuations and the rate of dissipation for highly compressible turbulence. We test closures for the turbulence energy cascade by filtering data from high-resolution simulations of forced isothermal and adiabatic turbulence. Optimal properties and an excellent correlation are found for a linear combination of the eddy-viscosity closure that is employed in LES of weakly compressible turbulence and a term that is non-linear in the Jacobian matrix of the velocity. Using this mixed closure, the SGS turbulence energy model is validated in LES of turbulence with stochastic forcing. It is found that the SGS model satisfies several important requirements: 1. The mean SGS turbulence energy follows a power law for varying grid scale. 2. The root mean square (RMS) Mach number of the unresolved velocity fluctuations is proportional to the RMS Mach number of the resolved turbulence, independent of the forcing. 3. The rate of dissipation and the turbulence energy flux are constant. Moreover, we discuss difficulties with direct estimates of the turbulent pressure and the dissipation rate on the basis of resolved flow quantities that have recently been proposed. In combination with the energy injection by stellar feedback and other unresolved processes, the proposed SGS model is applicable to a variety of problems in computational astrophysics. Computing the SGS turbulence energy, the treatment of star formation and stellar feedback in galaxy simulations can be improved. Further, we expect that the turbulent pressure on the grid scale affects the stability of gas against gravitational collapse.

We implemented sink particles in the Adaptive Mesh Refinement (AMR) code FLASH to model the gravitational collapse and accretion in turbulent molecular clouds and cores. Sink particles are frequently used to measure properties of star formation in numerical simulations, such as the star formation rate and efficiency, and the mass distribution of stars. We show that only using a density threshold for sink particle creation is insufficient in case of supersonic flows, because the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional physical collapse indicators have to be considered. We apply our AMR sink particle module to the formation of a star cluster, and compare it to a Smoothed Particle Hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas and sink particle properties between the AMR and SPH code.

Cosmological hydrodynamical simulations of primordial star formation suggest that the gas within the first star-forming halos is turbulent. This has strong implications on the subsequent evolution, in particular on the generation of magnetic fields. Using high-resolution numerical simulations, we show that in the presence of turbulence, weak seed magnetic fields are exponentially amplified by the small-scale dynamo during the formation of the first stars. We conclude that strong magnetic fields are generated during the birth of the first stars in the universe, potentially modifying the mass distribution of these stars and influencing the subsequent cosmic evolution. We find that the presence of the small-scale turbulent dynamo can only be identified in numerical simulations in which the turbulent motions in the central core are resolved with at least 32 grid cells.

The formation of massive stars is a highly complex process in which it is not clear whether the star-forming gas is in global gravitational collapse or in an equilibrium state, supported by turbulence. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, the filament containing the well-known sources DR21 and DR21(OH), we expect to find observational signatures that allow to discriminate between the two views. We use molecular line data from our 13CO 1-0, CS 2-1, and N2H+ 1-0 survey of the Cygnus X region obtained with the FCRAO and high-angular resolution observations of CO, CS, HCO+, N2H+, and H2CO, obtained with the IRAM 30m telescope. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of highest column-density, i.e. dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to rotation. Infall signatures in optically thick line profiles of HCO+ and 12CO are observed along and across the whole DR21 filament. From modelling the observed spectra, we obtain a typical infall speed of 0.6 km/s and mass accretion rates of the order of a few 10^-3 Msun/yr for the two main clumps constituting the filament. These massive (4900 and 3300 Msun) clumps are both gravitationally contracting. All observed kinematic features in the DR21 filament can be explained if it is formed by the convergence of flows at large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity at larger scales or from other processes can not be settled with the present study. The observed velocity field and velocity dispersion are consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent flows.

We compare the statistics of driven, supersonic turbulence at high Mach number using FLASH a widely used Eulerian grid-based code and PHANTOM, a Lagrangian SPH code at resolutions of up to 512^3 in both grid cells and SPH particles. We find excellent agreement between codes on the basic statistical properties: a slope of k^-1.95 in the velocity power spectrum for hydrodynamic, Mach 10 turbulence, evidence in both codes for a Kolmogorov-like slope of k^-5/3 in the variable rho^1/3 v as suggested by Kritsuk et al. and a log-normal PDF with a width that scales with Mach number and proportionality constant b=0.33-0.5 in the density variance-Mach number relation. The measured structure function slopes are not converged in either code at 512^3 elements. We find that, for measuring volumetric statistics such as the power spectrum slope and structure function scaling, SPH and grid codes give roughly comparable results when the number of SPH particles is approximately equal to the number of grid cells. In particular, to accurately measure the power spectrum slope in the inertial range, in the absence of subgrid models, requires at least 512^3 computational elements in either code. On the other hand the SPH code was found to be better at resolving dense structures, giving max. densities at a resolution of 128^3 particles that were similar those resolved in the grid code at 512^3 cells, reflected also in the high density tail of the PDF. We find SPH to be more dissipative at comparable numbers of computational elements in statistics of the velocity field, but correspondingly less dissipative than the grid code in the statistics of density weighted quantities such as rho^1/3 v. For SPH simulations of high Mach number turbulence we find it important to use sufficient non-linear beta-viscosity to prevent particle interpenetration in shocks (we require beta=4 instead of the default beta=2).

We study two limiting cases of turbulence forcing in numerical experiments: solenoidal (divergence-free) forcing, and compressive (curl-free) forcing, and compare our results to observations reported in the literature. We solve the equations of hydrodynamics on grids with up to 1024^3 cells for purely solenoidal and purely compressive forcing. Eleven lower-resolution models with mixtures of both forcings are also analysed. We find velocity dispersion--size relations consistent with observations and independent numerical simulations, irrespective of the type of forcing. However, compressive forcing yields stronger turbulent compression at the same RMS Mach number than solenoidal forcing, resulting in a three times larger standard deviation of volumetric and column density probability distributions (PDFs). We conclude that the strong dependence of the density PDF on the type of forcing must be taken into account in any theory using the PDF to predict properties of star formation. We supply a quantitative description of this dependence. We find that different observed regions show evidence of different mixtures of compressive and solenoidal forcing, with more compressive forcing occurring primarily in swept-up shells.

We introduce a method for calculating the probability density function (PDF) of a turbulent density field in three dimensions using only information contained in the projected two-dimensional column density field. We test the method by applying it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds. To a good approximation, the PDF of log(normalised column density) is a compressed, shifted version of the PDF of log(normalised density). The degree of compression can be determined observationally from the column density power spectrum, under the assumption of statistical isotropy of the turbulence.

We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu (1977), and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a sole density threshold for sink particle creation is insufficient in supersonic flows, if the density threshold is below the opacity limit. In that case, the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.

In this paper we discuss recent applications of the Smoothed Particle Hydrodynamics (SPH) method to the simulation of supersonic turbulence in the interstellar medium, as well as giving an update on recent algorithmic developments in solving the equations of magnetohydrodynamics (MHD) in SPH. Using high resolution calculations (up to 134 million particles), we find excellent agreement with grid-based results on a range of measures including the power spectrum slope in both the velocity field and the density-weighted velocity rho^(1/3) v, the latter showing a Kolmogorov-like k^-5/3 scaling as proposed by Kritsuk et al. (2007). We also find good agreement on the statistics of the Probability Distribution Function (PDF) and structure functions, independently confirming the scaling found by Schmidt, Federrath & Klessen (2008). On Smoothed Particle Magnetohydrodynamics (SPMHD) we have recently wasted a great deal of time and effort investigating the vector potential as an alternative to the Euler potentials formulation, in the end concluding that using the vector potential has even more severe problems than the standard (B-field based) SPMHD approach.

We investigate the influence of the turbulence forcing on the mass distributions of gravitationally unstable cores by postprocessing data from simulations of non-selfgravitating isothermal supersonic turbulence with varying resolution. In one set of simulations solenoidal forcing is applied, while the second set uses purely compressive forcing to excite turbulent motions. From the resulting density field, we compute the mass distribution of gravitationally unstable cores by means of a clump-finding algorithm. Using the time-averaged probability density functions of the mass density, semi-analytic mass distributions are calculated from analytical theories. We apply stability criteria that are based on the Bonnor-Ebert mass resulting from the thermal pressure and from the sum of thermal and turbulent pressure. Although there are uncertainties in the application of the clump-finding algorithm, we find systematic differences in the mass distributions obtained from solenoidal and compressive forcing. Compressive forcing produces a shallower slope in the high-mass power-law regime compared to solenoidal forcing. The mass distributions also depend on the Jeans length resulting from the choice of the mass in the computational box, which is freely scalable for non-selfgravitating isothermal turbulence. Provided that all cores are numerically resolved and most cores are small compared to the length scale of the forcing, the normalised core mass distributions are found to be close to the semi-analytic models. Especially for the high-mass tails, the Hennebelle-Chabrier theory implies that the additional support due to turbulent pressure is important.

We introduce and test an expression for calculating the variance of a physical field in three dimensions using only information contained in the two-dimensional projection of the field. The method is general but assumes statistical isotropy. To test the method we apply it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds, and demonstrate that it can recover the 3D normalised density variance with ~10% accuracy if the assumption of isotropy is valid. We show that the assumption of isotropy breaks down at low sonic Mach number if the turbulence is sub-Alfvenic. Theoretical predictions suggest that the 3D density variance should increase proportionally to the square of the Mach number of the turbulence. Application of our method will allow this prediction to be tested observationally and therefore constrain a large body of analytic models of star formation that rely on it.

We analyze the statistics and star formation rate obtained in high-resolution numerical experiments of forced supersonic turbulence, and compare with observations. We concentrate on a systematic comparison of solenoidal (divergence-free) and compressive (curl-free) forcing, which are two limiting cases of turbulence driving. Our results show that for the same RMS Mach number, compressive forcing produces a three times larger standard deviation of the density probability distribution. When self-gravity is included in the models, the star formation rate is more than one order of magnitude higher for compressive forcing than for solenoidal forcing.

We present results from high-resolution three-dimensional simulations of turbulent interstellar gas that self-consistently follow its coupled thermal, chemical and dynamical evolution, with a particular focus on the formation and destruction of H2 and CO. We quantify the formation timescales for H2 and CO in physical conditions corresponding to those found in nearby giant molecular clouds, and show that both species form rapidly, with chemical timescales that are comparable to the dynamical timescale of the gas. We also investigate the spatial distributions of H2 and CO, and how they relate to the underlying gas distribution. We show that H2 is a good tracer of the gas distribution, but that the relationship between CO abundance and gas density is more complex. The CO abundance is not well-correlated with either the gas number density n or the visual extinction A_V: both have a large influence on the CO abundance, but the inhomogeneous nature of the density field produced by the turbulence means that n and A_V are only poorly correlated. There is a large scatter in A_V, and hence CO abundance, for gas with any particular density, and similarly a large scatter in density and CO abundance for gas with any particular visual extinction. This will have important consequences for the interpretation of the CO emission observed from real molecular clouds. Finally, we also examine the temperature structure of the simulated gas. We show that the molecular gas is not isothermal. Most of it has a temperature in the range of 10--20 K, but there is also a significant fraction of warmer gas, located in low-extinction regions where photoelectric heating remains effective.

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

In a systematic study, we compare the density statistics in high resolution numerical experiments of supersonic isothermal turbulence, driven by the usually adopted solenoidal (divergence-free) forcing and by compressive (curl-free) forcing. We find that for the same rms Mach number, compressive forcing produces much stronger density enhancements and larger voids compared to solenoidal forcing. Consequently, the Fourier spectra of density fluctuations are significantly steeper. This result is confirmed using the Delta-variance analysis, which yields power-law exponents beta~3.4 for compressive forcing and beta~2.8 for solenoidal forcing. We obtain fractal dimension estimates from the density spectra and Delta-variance scaling, and by using the box counting, mass size and perimeter area methods applied to the volumetric data, projections and slices of our turbulent density fields. Our results suggest that compressive forcing yields fractal dimensions significantly smaller compared to solenoidal forcing. However, the actual values depend sensitively on the adopted method, with the most reliable estimates based on the Delta-variance, or equivalently, on Fourier spectra. Using these methods, we obtain D~2.3 for compressive and D~2.6 for solenoidal forcing, which is within the range of fractal dimension estimates inferred from observations (D~2.0-2.7). The velocity dispersion to size relations for both solenoidal and compressive forcing obtained from velocity spectra follow a power law with exponents in the range 0.4-0.5, in good agreement with previous studies.

The probability density function (PDF) of the gas density in turbulent supersonic flows is investigated with high-resolution numerical simulations. In a systematic study, we compare the density statistics of compressible turbulence driven by the usually adopted solenoidal forcing (divergence-free) and by compressive forcing (curl-free). Our results are in agreement with studies using solenoidal forcing. However, compressive forcing yields a significantly broader density distribution with standard deviation ~3 times larger at the same rms Mach number. The standard deviation-Mach number relation used in analytical models of star formation is reviewed and a modification of the existing expression is proposed, which takes into account the ratio of solenoidal and compressive modes of the turbulence forcing.

The statistical properties of turbulence are considered to be universal at sufficiently small length scales, i. e., independent of boundary conditions and large-scale forces acting on the fluid. Analyzing data from numerical simulations of supersonic turbulent flow driven by external forcing, we demonstrate that this is not generally true for the two-point velocity statistics of compressible turbulence. However, a reformulation of the refined similarity hypothesis in terms of the mass-weighted velocity rho^(1/3)v yields scaling laws that are almost insensitive to the forcing. The results imply that the most intermittent dissipative structures are shocks closely following the scaling of Burgers turbulence.

We performed numerical simulations of supersonic isothermal turbulence driven by mostly compressive large-scale forcing, using both a static grid and adaptive mesh refinement with an effective resolution N=768^3. After a transient phase dominated by shocks, turbulence evolves into a steady state with an RMS Mach number about 2.5, in which cloud-like structures of over-dense gas are surrounded by highly rarefied gas. The index of the turbulence energy spectrum function beta = 2.0 in the shock-dominated phase. As the flow approaches statistical equilibrium, the spectrum flattens, with beta = 1.9. For the scaling exponent of the root mean square velocity fluctuation, we obtain gamma = 0.43 from the velocity structure functions of second order. These results are well within the range of observed scaling properties for the velocity dispersion in molecular clouds. Calculating structure functions of order p=1,...,5, we find for all scaling exponents significant deviations from the Kolmogorov-Burgers model proposed by Boldyrev. Our results are very well described by a general log-Poisson model with a higher degree of intermittency, which implies an influence of the forcing on the scaling properties. Contrary to previous numerical results for isothermal turbulence, we obtain a skewed probability density function of the mass density fluctuations that is not consistent with log-normal statistics and entails a substantially higher fraction of mass in the density peaks than implied by the Padoan-Nordlund relation between the variance of the density fluctuations and the Mach number. In conclusion, it seems necessary to account for the production mechanism of turbulence in the ISM.

We use 3-dimensional numerical simulations of self-gravitating compressible turbulent gas in combination with Lagrangian tracer particles to investigate the mixing process of molecular hydrogen (H2) in interstellar clouds. Tracer particles are used to represent shock-compressed dense gas, which is associated with H2. We deposit tracer particles in regions of density contrast in excess of ten times the mean density. Following their trajectories and using probability distribution functions, we find an upper limit for the mixing timescale of H2, which is of order 0.3 Myr. This is significantly smaller than the lifetime of molecular clouds, which demonstrates the importance of the turbulent mixing of H2 as a preliminary stage to star formation.