We present synthetic dust polarization maps of simulated molecular clouds (MCs) with the goal to systematically explore the origin of the relative orientation of the magnetic field ($\bf{B}$) with respect to the MC sub-structures identified in density ($n$; 3D) and column density ($N$; 2D). The polarization maps are generated with the radiative transfer code POLARIS, including self-consistently calculated efficiencies for radiative torque alignment. The MCs are formed in two sets of 3D MHD simulations: in (i) colliding flows (CF), and (ii) the SILCC-Zoom simulations. In 3D, for the CF simulations with an initial field strength below $\sim$5 $\mu$G, $\bf{B}$ is oriented parallel or randomly with respect to the $n$-structures. For CF runs with stronger initial fields and all SILCC-Zoom simulations, which have an initial field strength of 3 $\mu$G, a flip from parallel to perpendicular orientation occurs at high densities of $n_\text{trans}$ $\simeq$ 10$^2$ - 10$^3$ cm$^{-3}$. We suggest that this flip happens if the MC's mass-to-flux ratio, $\mu$, is close to or below the critical value of 1. This corresponds to a field strength around 3 - 5 $\mu$G. In 2D, we use the Projected Rayleigh Statistics (PRS) to study the orientation of $\bf{B}$. If present, the flip in orientation occurs at $N_\text{trans}$ $\simeq$ 10$^{21 - 21.5}$ cm$^{-2}$, similar to the observed transition value from sub- to supercritical magnetic fields in the ISM. However, projection effects can reduce the power of the PRS method: Depending on the MC or LOS, the projected maps of the SILCC-Zoom simulations do not always show the flip, although expected from the 3D morphology. Such projection effects can explain the variety of recently observed field configurations, in particular within a single MC. Finally, we do not find a correlation between the observed orientation of $\bf{B}$ and the $N$-PDF.

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.

Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM.

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($\Sigma_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $\alpha\sim-1.7$ ($\alpha\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $\Sigma_\mathrm{SFR}$. The cluster formation efficiency varies from $\Gamma\sim20\%$ in early merger phases to $\Gamma\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $\sigma\sim0.1-2$ Myr ($1\sigma$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.

We analyse the CO-dark molecular gas content of simulated molecular clouds from the SILCC-Zoom project. The simulations reach a resolution of 0.1 pc and include H$_2$ and CO formation, radiative stellar feedback and magnetic fields. CO-dark gas is found in regions with local visual extinctions $A_\text{V, 3D} \sim$ 0.2 - 1.5, number densities of 10 - 10$^{3}$ cm$^{-3}$ and gas temperatures of few 10 K - 100 K. CO-bright gas is found at number densities above 300 cm$^{-3}$ and temperatures below 50 K. The CO-dark gas fractions range from 40% to 95% and scale inversely with the amount of well-shielded gas ($A_\text{V, 3D}$ $\gtrsim$ 1.5), which is smaller in magnetised molecular clouds. We show that the density, chemical abundances and $A_\text{V, 3D}$ along a given line-of-sight cannot be properly determined from projected quantities. As an example, pixels with a projected visual extinction of $A_\text{V, 2D} \simeq$ 2.5 - 5 can be both, CO-bright or CO-dark, which can be attributed to the presence or absence of strong density enhancements along the line-of-sight. By producing synthetic CO(1-0) emission maps of the simulations with RADMC-3D, we show that about 15 - 65\% of the H$_2$ is in regions with intensities below the detection limit. Our clouds have $X_\text{CO}$-factors around 1.5 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ with a spread of up to a factor $\sim$ 4, implying a similar uncertainty in the derived total H$_2$ masses and even worse for individual pixels. Based on our results, we suggest a new approach to determine the H$_2$ mass, which relies on the availability of CO(1-0) emission and $A_\text{V, 2D}$ maps. It reduces the uncertainty of the clouds' overall H$_2$ mass to a factor of $\lesssim$ 1.8 and for individual pixels, i.e. on sub-pc scales, to a factor of $\lesssim$ 3.

We present a hydro-dynamical simulation at sub-parsec and few solar mass resolution of a merger between two gas-rich dwarf galaxies. Our simulation includes a detailed model for the multi-phase interstellar medium (ISM) and is able to follow the entire formation history of spatially resolved star clusters, including feedback from individual massive stars. Shortly after the merger we find a population of $\sim 900$ stellar clusters with masses above $10^{2.5}\; \rm{M_\odot}$ and a cluster mass function (CMF), which is well fitted with a power-law with a slope of $\alpha=-1.70\pm0.08$. We describe here in detail the formation of the three most massive clusters (M$_{*} \gtrsim 10^5$ M$_\odot$), which populate the high-mass end of the CMF. The simulated clusters form rapidly on a timescale of $6$-$8$ Myr in converging flows of dense gas. The embedded merger phase has extremely high star formation rate surface densities of $\Sigma_\mathrm{SFR}>10\; \mathrm{M}_\odot\; \mathrm{yr}^{-1}\; \mathrm{kpc}^{-2}$ and thermal gas pressures in excess of $P_{\rm th}\sim10^7 \; \mathrm{k}_{\rm B}\;(\rm K\;\mathrm{cm}^{-3})^{-1}$. The formation process is terminated by rapid gas expulsion driven by the first generation of supernovae, after which the cluster centers relax and both their structure and kinematics become indistinguishable from observed local globular clusters. The simulation presented here provides a general model for the formation of metal-poor globular clusters in chemically unevolved starbursting environments of low-mass dwarf galaxies, which are common at high redshifts.

Sources of X-rays such as active galactic nuclei and X-ray binaries are often variable by orders of magnitude in luminosity over timescales of years. During and after these flares the surrounding gas is out of chemical and thermal equilibrium. We introduce a new implementation of X-ray radiative transfer coupled to a time-dependent chemical network for use in 3D magnetohydrodynamical simulations. A static fractal molecular cloud is irradiated with X-rays of different intensity, and the chemical and thermal evolution of the cloud are studied. For a simulated $10^5$ M$_\odot$ fractal cloud an X-ray flux $<0.01$ erg cm$^{-2}$ s$^{-1}$ allows the cloud to remain molecular, whereas most of the CO and H$_2$ are destroyed for a flux of $>1$ erg cm$^{-2}$ s$^{-1}$. The effects of an X-ray flare, which suddenly increases the X-ray flux by $10^5 \times$ are then studied. A cloud exposed to a bright flare has 99% of its CO destroyed in 10-20 years, whereas it takes $>10^3$ years for 99% of the H$_2$ to be destroyed. CO is primarily destroyed by locally generated far-UV emission from collisions between non-thermal electrons and H$_2$; He$^+$ only becomes an important destruction agent when the CO abundance is already very small. After the flare is over, CO re-forms and approaches its equilibrium abundance after $10^3-10^5$ years. This implies that molecular clouds close to Sgr A$^*$ in the Galactic Centre may still be out of chemical equilibrium, and we predict the existence of clouds near flaring X-ray sources in which CO has been mostly destroyed but H is fully molecular.

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.

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

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

As part of the SILCC-ZOOM project we present our first sub-parsec resolution radiation-hydrodynamic simulations of two molecular clouds self-consistently forming from a turbulent, multi-phase ISM. The clouds have similar initial masses of few 10$^4$ M$_{\odot}$, escape velocities of ~5 km s$^{-1}$, and a similar initial energy budget. We follow the formation of star clusters with a sink based model and the impact of radiation from individual massive stars with the tree-based radiation transfer module TreeRay. Photo-ionizing radiation is coupled to a chemical network to follow gas heating, cooling and molecule formation and dissociation. For the first 3 Myr of cloud evolution we find that the overall star formation effciency is considerably reduced by a factor of ~4 to global cloud values of < 10 % as the mass accretion of sinks that host massive stars is terminated after <1 Myr. Despite the low effciency, star formation is triggered across the clouds. Therefore, a much larger region of the cloud is affected by radiation and the clouds begin to disperse. The time scale on which the clouds are dispersed sensitively depends on the cloud substructure and in particular on the amount of gas at high visual extinction. The damage of radiation done to the highly shielded cloud (MC1) is delayed. We also show that the radiation input can sustain the thermal and kinetic energy of the clouds at a constant level. Our results strongly support the importance of ionizing radiation from massive stars for explaining the low observed star formation effciency of molecular clouds

Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin of these features, we have developed a new episodic outflow model for the SPH code {\sc gandalf}, which mimics the accretion and ejection behaviour of FU Ori type stars. We apply this model to simulations of star formation, invoking two types of initial conditions: spherically symmetric cores in solid-body rotation with $\rho\propto r^{-2}$, and spherically symmetric turbulent cores with density proportional to the density of a Bonnor-Ebert sphere. For a wide range of model parameters, we find that episodic outflows lead to self-regulation of the ejected mass and momentum, and we achieve acceptable results, even with relatively low resolution. Using this model, we find that recently ejected outflow bullets produce a `Hubble wedge' in the position-velocity relation. However, once such a bullet hits the leading shock front, it decelerates and aligns with older bullets to form a `Hubble-law'. Bullets can be identified as bumps in the mass-velocity relation, which can be fit with a power-law, $dM/d\upsilon_{_{\rm RAD}}\propto\upsilon_{_{\rm RAD}}^{-1.5}$.

The formation of molecular hydrogen (H$_2$) and carbon monoxide (CO) is sensitive to the volume and column density distribution of the turbulent interstellar medium. In this paper, we study H$_2$ and CO formation in a large set of hydrodynamical simulations of periodic boxes with driven supersonic turbulence, as well as in colliding flows with the \textsc{Flash} code. The simulations include a non-equilibrium chemistry network, gas self-gravity, and diffuse radiative transfer. We investigate the spatial resolution required to obtain a converged H$_2$ and CO mass fraction and formation history. From the numerical tests we find that H$_2$ converges at a spatial resolution of $\lesssim0.2$~pc, while the required resolution for CO convergence is $\lesssim 0.04$~pc in gas with solar metallicity which is subject to a solar neighbourhood interstellar radiation field. We derive two critical conditions from our numerical results: the simulation has to at least resolve the densities at which (1) the molecule formation time in each cell in the computational domain is equal to the dissociation time, and (2) the formation time is equal to the the typical cell crossing time. For both H$_2$ and CO, the second criterion is more restrictive. The formulae we derive can be used to check whether molecule formation is converged in any given simulation.

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

The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{{\rm CII}}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{{\rm CII}}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{{\rm CII}}$ to change in environments with a lower or higher radiation field than simulated here.

Magnetic fields are ubiquitously observed in the interstellar medium (ISM) of present-day star-forming galaxies with dynamically relevant energy densities. Using three-dimensional magneto-hydrodynamic (MHD) simulations of the supernova (SN) driven ISM in the flux-freezing approximation (ideal MHD) we investigate the impact of the magnetic field on the chemical and dynamical evolution of the gas, fragmentation and the formation of molecular clouds. We follow the chemistry with a network of six species (H$^{+}$, H, H$_2$, C$^+$, CO, free electrons) including local shielding effects. We find that magnetic fields thicken the disc by a factor of a few to a scale height of $\sim100\,\mathrm{pc}$, delay the formation of dense (and molecular) gas by $\sim25\,\mathrm{Myr}$ and result in differently shaped gas structures. The magnetised gas fragments into fewer clumps, which are initially at subcritical mass-to-flux ratios, $M/\Phi\approx0.3(M/\Phi)_\mathrm{crit}$, and accrete gas preferentially parallel to the magnetic field lines until supercritial mass-to-flux ratios of up to order 10 are reached. The accretion rates onto molecular clouds scale with $\dot{M}\propto M^{1.5}$. The median of the inter-cloud velocity dispersion is $\sim2-5\,\mathrm{km\,s}^{-1}$ and lower than the internal velocity dispersion in the clouds ($\sim3-7\,\mathrm{km\,s}^{-1}$). However, individual cloud-cloud collisions occur at speeds of a few $10\,\mathrm{km\,s}^{-1}$.

Recent magneto-hydrodynamic simulations of the star-forming interstellar medium (ISM) with parsec scale resolution indicate that relativistic cosmic rays support the launching of galactic outflows on scales of a few kpc. If these fundamental constituents of the ISM are injected at the sites of supernova (SN) explosions, the outflows are smoother, colder, and denser than the highly structured, hot-phase driven outflows forming, e.g., by thermal SN energy injection alone. In this Letter we present computations of resonant Lyman-$\alpha$ (Ly$\alpha$) radiation transfer through snapshots of a suite of stratified disk simulations from the SILCC project. For a range of thermal, radiative, and kinetic feedback models only simulations including non-thermal cosmic rays produce Ly$\alpha$ spectra with enhanced red peaks and strong absorption at line center -- similar to observed systems. The absence of cosmic ray feedback leads to spectra incompatible with observations. We attribute this to the smoother neutral gas distribution of cosmic ray supported outflows within a few kpc from the disk midplane.

Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a turbulent medium. Synthetic C$^{18}$O observations of the simulated filaments reveal the existence of fibres very similar to the observed ones, i.e. elongated coherent features in the resulting PPV data-cubes. Analysis of the PPP data-cubes (i.e. 3D density fields) also reveals elongated coherent features, which we term sub-filaments. Unfortunately there is very poor correspondence between the fibres and the sub-filaments in the simulations. Both fibres and sub-filaments derive from inhomogeneities in the turbulent accretion flow onto the main filament. As a consequence, fibres are often affected by line-of-sight confusion. Similarly, sub-filaments are often affected by large velocity gradients, and even velocity discontinuities. These results suggest that extreme care should be taken when using velocity coherent features to constrain the underlying substructure within a filament.

We investigate the impact of cosmic rays (CRs) on galactic outflows from a multi-phase interstellar medium with solar neighbourhood conditions. The three-dimensional magneto-hydrodynamical simulations include CRs as a relativistic fluid in the advection-diffusion approximation. The thermal and chemical state of the ISM is computed with a non-equilibrium chemical network. We find that CRs (injected with 10 \% of the supernova energy) efficiently support the launching of outflows and strongly affect their phase structure. Outflows leaving the midplane are denser ($\rho \sim 10^{-26}\,\mathrm{g\,cm}^{-3}$), colder ($\sim 10^4\,\mathrm{K}$), and slower ($\sim 30\,\mathrm{km\,s}^{-1}$) if CRs are considered in addition to thermal SNe. The CR supported outflows are also smoother, in particular at larger heights ($> 1\,\mathrm{kpc}$ above the midplane) without the direct impact of SN explosions. Approximately $5\% - 25\%$ of the injected CR energy is lost via hadronic cooling. Smaller diffusion coefficients lead to slightly larger hadronic losses but allow for steeper CR pressure gradients, stronger outflows and larger accelerations. Up to a height of $z \sim1\,\mathrm{kpc}$ there are large volumes in approximate pressure equilibrium between thermal and CR component. At larger altitudes the CR pressure is $10-100$ times as large as the thermal counterpart. More than $\sim 1\,\mathrm{kpc}$ away from the midplane, CRs provide the dominant gas acceleration mechanism.

Photoionizing radiation and stellar winds from massive stars deposit energy and momentum into the interstellar medium (ISM). They might disperse the local ISM, change its turbulent multi-phase structure, and even regulate star formation. Ionizing radiation dominates the massive stars' energy output, but the relative effect of winds might change with stellar mass and the properties of the ambient ISM. We present simulations of the interaction of stellar winds and ionizing radiation of 12, 23, and 60 M$_{\odot}$ stars within a cold neutral (CNM, $n_{0}$ = 100 cm$^{-3}$), warm neutral (WNM, $n_{0}$ = 1, 10 cm$^{-3}$) or warm ionized (WIM, $n_{0}$ = 0.1 cm$^{-3}$) medium. The FLASH simulations adopt the novel tree-based radiation transfer algorithm TreeRay. With the On-the-Spot approximation and a temperature-dependent recombination coefficient, it is coupled to a chemical network with radiative heating and cooling. In the homogeneous CNM, the total momentum injection ranges from 1.6$\times$10$^{4}$ to 4$\times$10$^{5}$ M$_{\odot}$ km s$^{-1}$ and is always dominated by the expansion of the ionized H$_{\text{II}}$ region. In the WIM, stellar winds dominate (2$\times$10$^{2}$ to 5$\times$10$^{3}$ M$_{\odot}$ km s$^{-1}$), while the input from radiation is small ($\sim$ 10$^{2}$ M$_{\odot}$ km s$^{-1}$). The WNM ($n_{0}$ = 1 cm$^{-3}$) is a transition regime. Energetically, stellar winds couple more efficiently to the ISM ($\sim$ 0.1 percent of wind luminosity) than radiation ($<$ 0.001 percent of ionizing luminosity). For estimating the impact of massive stars, the strongly mass-dependent ratios of wind to ionizing luminosity and the properties of the ambient medium have to be considered.

We present high-resolution ($\sim$ 0.1 pc), hydrodynamical and magnetohydrodynamical simulations to investigate whether the observed level of molecular cloud (MC) turbulence can be generated and maintained by external supernova (SN) explosions. The MCs are formed self-consistently within their large-scale galactic environment following the non-equilibrium formation of H$_2$ and CO including (self-) shielding and important heating and cooling processes. The MCs inherit their initial level of turbulence from the diffuse ISM, where turbulence is injected by SN explosions. However, by systematically exploring the effect of individual SNe going off outside the clouds, we show that at later stages the importance of SN driven turbulence is decreased significantly. This holds for different MC masses as well as MCs with and without magnetic fields. The SN impact also decreases rapidly with larger distances. Nearby SNe ($d$ $\sim$ 25 pc) boost the turbulent velocity dispersions of the MC by up to 70 per cent (up to a few km s$^{-1}$). For $d$ $>$ 50 pc, however, their impact decreases fast with increasing $d$ and is almost negligible. For all probed distances the gain in velocity dispersion decays rapidly within a few 100 kyr. This is significantly shorter than the average timescale for a MC to be hit by a nearby SN under solar neighbourhood conditions ($\sim$ 2 Myr). Hence, at these conditions SNe are not able to sustain the observed level of MC turbulence. However, in environments with high gas surface densities and SN rates like the Central Molecular Zone, observed elevated MC dispersions could be triggered by external SNe.

We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this approach is a relatively low memory requirement, important in particular for the optical depth calculation, which needs to process information from many different directions. This feature also enables a general tree-based radiation transport algorithm that will be described in a subsequent paper, and delivers excellent scaling up to at least 1500 cores. Boundary conditions for gravity can be either isolated or periodic, and they can be specified in each direction independently, using a newly developed generalisation of the Ewald method. The gravity calculation can be accelerated with the {\em adaptive block update} technique by partially re-using the solution from the previous time-step. Comparison with the {\sc Flash} internal multi-grid gravity solver shows that tree based methods provide a competitive alternative, particularly for problems with isolated or mixed boundary conditions. We evaluate several multipole acceptance criteria (MACs) and identify a relatively simple APE MAC which provides high accuracy at low computational cost. The optical depth estimates are found to agree very well with those of the {\sc RADMC-3D} radiation transport code, with the tree solver being much faster. Our algorithm is available in the standard release of the {\sc FLASH} code in version 4.0 and later.

We present 3D "zoom-in" simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds - identified from the SILCC simulations - with a resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling, and the formation of H$_2$ and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to $\sim$ 10$^{-2}$ M$_\mathrm{sun}$ yr$^{-1}$ and final masses of $\sim$ 50 000 M$_\mathrm{sun}$. A spatial resolution of $\lesssim$ 0.1 pc is required for convergence with respect to the mass, velocity dispersion, and chemical abundances of the clouds, although these properties also depend on the cloud definition such as based on density thresholds, H$_2$ or CO mass fraction. To avoid grid artefacts, the progressive increase of resolution has to occur within the free-fall time of the densest structures (1 - 1.5 Myr) and $\gtrsim$ 200 time steps should be spent on each refinement level before the resolution is progressively increased further. This avoids the formation of spurious, large-scale, rotating clumps from unresolved turbulent flows. While CO is a good tracer for the evolution of dense gas with number densities $n \geq$ 300 cm$^{-3}$, H$_2$ is also found for $n \lesssim 30$ cm$^{-3}$ due to turbulent mixing and becomes dominant at column densities around 30 - 50 M$_\mathrm{sun}$ pc$^{-2}$. The CO-to-H$_2$ ratio steadily increases within the first 2 Myr whereas $X_\mathrm{CO}\simeq$ 1 - 4 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ is approximately constant since the CO(1-0) line quickly becomes optically thick.

We present synthetic continuum and $^{13}$CO and C$^{18}$O line emission observations of dense and cold filaments. The filaments are dynamically evolved using 3D-MHD simulations that include one of the largest on-the-fly chemical networks used to date, which models the detailed evolution of H$_2$ and CO. We investigate the reliability of observable properties, in particular filament mass and width, under different simulation conditions like magnetic field orientation and cosmic ray ionisation rate. We find that filament widths of $\sim$0.1 pc can be probed with both line and continuum emission observations with a high accuracy (deviations $\leq$ 20%). However, the width of more narrow filaments can be significantly overestimated by up to a factor of a few. Masses obtained via the dust emission are accurate within a few percent whereas the masses inferred from molecular line emission observations deviate from the actual mass by up to a factor of 10 and show large differences for different $J$ transitions. The inaccurate estimate of filament masses and widths of narrow filaments using molecular line observations can be attributed to (i) the non-isothermal state of the filaments, (ii) optical depth effects, and (iii) the subthermally excited state of CO, while inclination effects and opacity correction only influence the obtained masses and widths by less than 50%. Both, mass and width estimates, can be improved by using two isotopes to correct for the optical depth. Since gas and dust temperature generally differ (by up to 25 K), the filaments appear more gravitationally unstable if the (too low) dust temperature is used for the stability analysis.

We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible method to sample the stellar initial mass function allows us to follow the contribution to the ISRF, the metal output and the SN delay times of individual massive stars. We find that SNe play the dominant role in regulating the global star formation rate, shaping the multi-phase interstellar medium (ISM) and driving galactic outflows. Outflow rates (with mass-loading factors of a few) and hot gas fractions of the ISM increase with the number of SNe exploding in low-density environments where radiative energy losses are low. While PE heating alone can suppress star formation slightly more (a factor of a few) than SNe alone can do, it is unable to drive outflows and reproduce the multi-phase ISM that emerges naturally when SNe are included. These results are in conflict with recent results of Forbes et al. who concluded that PE heating is the dominant process suppressing star formation in dwarfs, about an order of magnitude more efficient than SNe. Potential origins for this discrepancy are discussed. In the absence of SNe and photoionization (mechanisms to disperse dense clouds), the impact of PE heating is highly overestimated owing to the (unrealistic) proximity of dense gas to the radiation sources. This leads to a substantial boost of the infrared continuum emission from the UV-irradiated dust and a far infrared line-to-continuum ratio too low compared to observations. Though sub-dominant in regulating star formation, the ISRF controls the abundance of molecular hydrogen via photodissociation.

We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3 $\mu$G). All magnetic runs saturate at mass weighted field strengths of $\sim$ 1 $-$ 3 $\mu$G but the ISM structure is notably different. With increasing initial field strengths (from $6\times 10^{-3}$ to 3 $\mu$G) the simulations develop an ISM with a more homogeneous density and temperature structure, with increasing mass (from 5% to 85%) and volume filling fractions (from 4% to 85%) of warm (300 K $<$ T $<$ 8000 K) gas, with decreasing volume filling fractions (VFF) from $\sim$ 35% to $\sim$ 12% of hot gas (T $> 10^5$ K) and with a decreasing H$_{2}$ mass fraction (from 70% to $<$ 1%). Meanwhile the mass fraction of gas in which the magnetic pressure dominates over the thermal pressure increases by a factor of 10, from 0.07 for an initial field of $6\times 10^{-3}$ $\mu$G to 0.7 for a 3 $\mu$G initial field. In all but the simulations with the highest initial field strength self-gravity promotes the formation of dense gas and H$_{2}$, but does not change any other trends. We conclude that magnetic fields have a significant impact on the multi-phase, chemical and thermal structure of the ISM and discuss potential implications and limitations of the model.

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionising radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed with a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionising radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionising radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to rho < 10^-25 g cm^-3, compared to 10^-23 g cm^-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Halpha emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 % to the total Halpha emission. Individual massive stars (M >= 30 M_sun) with short lifetimes are responsible for significant fluctuations in the Halpha luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by factors of ~10.

Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of H^+, H, H_2, C^+ and CO and take into account self-shielding by gas and dust using a tree-based radiation transfer method. Supernova explosions are injected either at random locations, at density peaks, or as a mixture of the two. For each simulation, we investigate how matter circulates between the ISM phases and find more sizeable transitions than considered in simple mass exchange schemes in the literature. The derived residence times in the ISM phases are characterised by broad distributions, in particular for the molecular, warm and hot medium. The most realistic simulations with random and mixed driving have median residence times in the molecular, cold, warm and hot phase around 17, 7, 44 and 1 Myr, respectively. The transition rates measured in the random driving run are in good agreement with observations of Ti gas-phase depletion in the warm and cold phases in a simple depletion model, although the depletion in the molecular phase is under-predicted. ISM phase definitions based on chemical abundance rather than temperature cuts are physically more meaningful, but lead to significantly different transition rates and residence times because there is no direct correspondence between the two definitions.

We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows.

We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code $\texttt{FLASH}$ (http://flash.uchicago.edu). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in $\texttt{FLASH}$.

Supernova (SN) blast waves inject energy and momentum into the interstellar medium (ISM), control its turbulent multiphase structure and the launching of galactic outflows. Accurate modelling of the blast wave evolution is therefore essential for ISM and galaxy formation simulations. We present an efficient method to compute the input of momentum, thermal energy, and the velocity distribution of the shock-accelerated gas for ambient media with uniform (and with stellar wind blown bubbles), power-law, and turbulent density distributions. Assuming solar metallicity cooling, the blast wave evolution is followed to the beginning of the momentum conserving snowplough phase. The model recovers previous results for uniform ambient media. The momentum injection in wind-blown bubbles depend on the swept-up mass and the efficiency of cooling, when the blast wave hits the wind shell. For power-law density distributions with $n(r) \sim$ $r^{-2}$ (for $n(r) > n_{_{\rm floor}}$) the amount of momentum injection is solely regulated by the background density $n_{_{\rm floor}}$ and compares to $n_{_{\rm uni}}$ = $n_{_{\rm floor}}$. However, in turbulent ambient media with log-normal density distributions the momentum input can increase by a factor of 2 (compared to the homogeneous case) for high Mach numbers. The average momentum boost can be approximated as $p_{_{\rm turb}}/\mathrm{p_{_{0}}}\ =23.07\, \left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-0.12} + 0.82 (\ln(1+b^{2}\mathcal{M}^{2}))^{1.49}\left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-1.6}$. The velocity distributions are broad as gas can be accelerated to high velocities in low-density channels. The model values agree with results from recent, computationally expensive, three-dimensional simulations of SN explosions in turbulent media.

We present synthetic ALMA observations of Keplerian, protostellar discs in the Class 0 stage studying the emission of molecular tracers like $^{13}$CO, C$^{18}$O, HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, and H$_2$CO. We model the emission of discs around low- and intermediate-mass protostars. We show that under optimal observing conditions ALMA is able to detect the discs already in the earliest stage of protostellar evolution, although the emission is often concentrated to the innermost 50 AU. Therefore, a resolution of a few 0.1" might be too low to detect Keplerian discs around Class 0 objects. We also demonstrate that under optimal conditions for edge-on discs Keplerian rotation signatures are recognisable, from which protostellar masses can be inferred. For this we here introduce a new approach, which allows us to determine protostellar masses with higher fidelity than before. Furthermore, we show that it is possible to reveal Keplerian rotation even for strongly inclined discs and that ALMA should be able to detect possible signs of fragmentation in face-on discs. In order to give some guidance for future ALMA observations, we investigate the influence of varying observing conditions and source distances. We show that it is possible to probe Keplerian rotation in inclined discs with an observing time of 2 h and a resolution of 0.1", even in the case of moderate weather conditions. Furthermore, we demonstrate that under optimal conditions, Keplerian discs around intermediate-mass protostars should be detectable up to kpc-distances.

Two-thirds of long duration gamma-ray bursts (GRBs) show soft X-ray absorption in excess of the Milky Way. The column densities of metals inferred from UV and optical spectra differ from those derived from soft X-ray spectra, at times by an order of magnitude, with the latter being higher. The origin of the soft X-ray absorption excess observed in GRB X-ray afterglow spectra remains a heavily debated issue, which has resulted in numerous investigations on the effect of hot material both internal and external to the GRB host galaxy on our X-ray afterglow observations. Nevertheless, all models proposed so far have either only been able to account for a subset of our observations (i.e. at z > 2), or they have required fairly extreme conditions to be present within the absorbing material. In this paper, we investigate the absorption of the GRB afterglow by a collisionally ionised and turbulent interstellar medium (ISM). We find that a dense (3 per cubic centimeters) collisionally ionised ISM could produce UV/optical and soft X-ray absorbing column densities that differ by a factor of 10, however the UV/optical and soft X-ray absorbing column densities for such sightlines and are 2-3 orders of magnitude lower in comparison to the GRB afterglow spectra. For those GRBs with a larger soft X-ray excess of up to an order of magnitude, the contribution in absorption from a turbulent ISM as considered here would ease the required conditions of additional absorbing components, such as the GRB circumburst medium and intergalactic medium.

We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $\Sigma_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.

We present simulations of star forming filaments incorporating on of the largest chemical network used to date on-the-fly in a 3D-MHD simulation. The network contains 37 chemical species and about 300 selected reaction rates. For this we use the newly developed package KROME (Grassi et al. 2014). We combine the KROME package with an algorithm which allows us to calculate the column density and attenuation of the interstellar radiation field necessary to properly model heating and ionisation rates. Our results demonstrate the feasibility of using such a complex chemical network in 3D-MHD simulations on modern supercomputers. We perform simulations with different strengths of the interstellar radiation field and the cosmic ray ionisation rate. We find that towards the centre of the filaments there is gradual conversion of hydrogen from H to H$_2$ as well as of C$^+$ over C to CO. Moreover, we find a decrease of the dust temperature towards the centre of the filaments in agreement with recent HERSCHEL observations.

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (\rho~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (\sigma_HI ~ 0.8\sigma_(300-8000K) ~ 10-20km/s, \sigma_H\alpha ~ 0.6\sigma_(8000-3e5K) ~ 20-30km/s), and strong outflows with mass loading factors of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (\rho ~ 1e-25-1e-24g/cc) and slow (v ~ 20-40km/s) but there is a high-velocity tail of up to v ~ 500km/s with \rho ~ 1e-28-1e-27g/cc.

We present simulations of star forming filaments incorporating - to our knowledge - the largest chemical network used to date on-the-fly in a 3D-MHD simulation. The network contains 37 chemical species and about 300 selected reaction rates. For this we use the newly developed package KROME (Grassi et al. 2014). Our results demonstrate the feasibility of using such a complex chemical network in 3D-MHD simulations on modern supercomputers. We perform simulations with different strengths of the interstellar radiation field and the cosmic ray ionisation rate and find chemical and physical results in accordance with observations and other recent numerical work.

The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

We present simulations of collapsing filaments studying the impact of turbulence and magnetic field morphologies on their evolution and star formation properties. We vary the mass per unit length of the filaments as well as the orientation of the magnetic field with respect to the major axis. We find that the filaments, which have no or a perpendicular magnetic field, typically reveal a smaller width than the universal width of 0.1 pc proposed by e.g. Arzoumanian et al. 2011. We show that this also holds in the presence of supersonic turbulence and that accretion driven turbulence is too weak to stabilize the filaments along their radial direction. On the other hand, we find that a magnetic field that is parallel to the major axis can stabilize the filament against radial collapse resulting in widths of 0.1 pc. Furthermore, depending on the filament mass and magnetic field configuration, gravitational collapse and fragmentation in filaments occurs either in an edge-on way, uniformly distributed across the entire length, or in a mixed way. In the presence of initially moderate density perturbations, a centralized collapse towards a common gravitational centre occurs. Our simulations can thus reproduce different modes of fragmentation observed recently in star forming filaments. Moreover, we find that turbulent motions influence the distance between individual fragments along the filament, which does not always match the results of a Jeans analysis.

Massive clumps within the swept-up shells of bubbles, like that surrounding the galactic HII region RCW 120, have been interpreted in terms of the Collect and Collapse (C&C) mechanism for triggered star formation. The cold, dusty clumps surrounding RCW 120 are arranged in an almost spherical shell and harbour many young stellar objects. By performing high-resolution, three-dimensional SPH simulations of HII regions expanding into fractal molecular clouds, we investigate whether the formation of massive clumps in dense, swept-up shells necessarily requires the C&C mechanism. In a second step, we use RADMC-3D to compute the synthetic dust continuum emission from our simulations, in order to compare them with observations of RCW 120 made with APEX-LABOCA at 870 micron. We show that a distribution of clumps similar to the one seen in RCW 120 can readily be explained by a non-uniform initial molecular cloud structure. Hence, a shell-like configuration of massive clumps does not imply that the C&C mechanism is at work. Rather, we find a hybrid form of triggering, which combines elements of C&C and Radiatively Driven Implosion (RDI). In addition, we investigate the reliability of deriving clump masses from their 870 micron emission. We find that for clumps with more than 100 M_sun the observational estimates are accurate to within a factor of two and that, even at these long wavelengths, it is important to account for the radiative heating from triggered, embedded protostars.

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

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $\Sigma_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H$_2$ and CO considering (self-) shielding, and supernova (SN) feedback. We explore SN explosions at different (fixed) rates in high-density regions (peak), in random locations (random), in a combination of both (mixed), or clustered in space and time (clustered). Only random or clustered models with self-gravity (which evolve similarly) are in agreement with observations. Molecular hydrogen forms in dense filaments and clumps and contributes 20% - 40% to the total mass, whereas most of the mass (55% - 75%) is in atomic hydrogen. The ionised gas contributes <10%. For high SN rates (0.5 dex above Kennicutt-Schmidt) as well as for peak and mixed driving the formation of H$_2$ is strongly suppressed. Also without self-gravity the H$_2$ fraction is significantly lower ($\sim$ 5%). Most of the volume is filled with hot gas ($\sim$90% within $\pm$2 kpc). Only for random or clustered driving, a vertically expanding warm component of atomic hydrogen indicates a fountain flow. Magnetic fields have little impact on the final disc structure. However, they affect dense gas ($n\gtrsim 10\;{\rm cm}^{-3}$) and delay H$_2$ formation. We highlight that individual chemical species, in particular atomic hydrogen, populate different ISM phases and cannot be accurately accounted for by simple temperature-/density-based phase cut-offs.

Lomax et al. have constructed an ensemble of 60 prestellar cores having masses, sizes, projected shapes, temperatures and non-thermal radial velocity dispersions that match, statistically, the cores in Ophiuchus; and have simulated the evolution of these cores using SPH. Each core has been evolved once with no radiative feedback from stars, once with continuous radiative feedback, and once with episodic radiative feedback. Here we analyse the multiplicity statistics from these simulations. With episodic radiative feedback, (i) the multiplicity frequency is ~60% higher than in the field; (ii) the multiplicity frequency and the mean semi-major axis both increase with primary mass; (iii) one third of multiple systems are hierarchical systems with more than two components; (iv) in these hierarchical systems the inner pairings typically have separations of a few au and mass ratios concentrated towards unity, whereas the outer pairings have separations of order 100 au and a flatter distribution of mass ratios. The binary statistics are compatible with observations of young embedded populations, and -- if wider orbits are disrupted preferentially by external perturbations -- with observations of mature field populations. With no radiative feedback, the results are similar to those from simulations with episodic feedback. With continuous radiative feedback, brown dwarfs are under-produced, the number of multiple systems is too low, and the statistical properties of multiple systems are at variance with observation. This suggests that star formation in Ophiuchus may only be representative of global star formation if accretion onto protostars, and hence radiative feedback, is episodic.

We investigate the early impact of single and binary supernova (SN) explosions on dense gas clouds with three-dimensional, high-resolution, hydrodynamic simulations. The effect of cloud structure, radiative cooling, and ionising radiation from the progenitor stars on the net input of kinetic energy, f_kin=E_kin/E_SN, thermal energy, f_therm=E_therm/E_SN, and gas momentum f_P=P/P_SN to the interstellar medium (ISM) is tested. For clouds with n=100 cm^{-3}, the momentum generating Sedov and pressure-driven snowplough phases are terminated early (~ 0.01 Myr) and radiative cooling limits the coupling to f_therm ~ 0.01, f_kin ~ 0.05, and f_P ~ 9, significantly lower than for the case without cooling. For pre-ionised clouds these numbers are only increased by ~ 50%, independent of the cloud structure. This only suffices to accelerate ~ 5% of the cloud to radial velocities >30km/s. A second SN might further enhance the coupling efficiencies if delayed past the Sedov phase of the first explosion. Such very low coupling efficiencies cast doubts on many galaxy-scale sub-resolution models for supernova feedback, most of which are validated a posteriori by qualitative agreement of galaxy properties with observations. Ionising radiation appears not to significantly enhance the immediate coupling of SNe to the surrounding gas as it drives the ISM into inert dense shells and cold clumps, a process which is unresolved in galaxy scale simulations. Our results support previous conclusions that supernovae might only drive a wind if a significant fraction explodes in low-density environments or if they are supported by processes other than ionising radiation.

We present the smoothed-particle hydrodynamics implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity gradients, a modified artificial viscosity switch with a modified strong limiter, and artificial conduction of thermal energy. With a series of idealized hydrodynamic tests we show that the pressure-entropy formulation is ideal for resolving fluid mixing at contact discontinuities but performs conspicuously worse at strong shocks due to the large entropy discontinuities. Including artificial conduction at shocks greatly improves the results. In simulations of Milky Way like disk galaxies a feedback-induced instability develops if too much artificial viscosity is introduced. Our modified artificial viscosity scheme prevents this instability and shows efficient shock capturing capability. We also investigate the star formation rate and the galactic outflow. The star formation rates vary slightly for different SPH schemes while the mass loading is sensitive to the SPH scheme and significantly reduced in our favored implementation. We compare the accretion behavior of the hot halo gas. The formation of cold blobs, an artifact of simple SPH implementations, can be eliminated efficiently with proper fluid mixing, either by conduction and/or by using a pressure-entropy formulation.

We present a novel implementation of cosmic rays (CR) in the magneto-hydrodynamic code FLASH. CRs are described as separate fluids with different energies. CR advection, energy dependent anisotropic diffusion with respect to the magnetic field and adiabatic losses to follow the evolution of spectra are taken into account. We present a first study of the transport and immediate (~150 kyr) dynamical impact of CRs on the turbulent magnetised interstellar medium around supernova remnants on scales up to 80 pc. CR diffusion quickly leads to an efficient acceleration of low-density gas (mainly perpendicular to the magnetic field) with accelerations up to two orders of magnitude above the thermal values. Peaked (at 1 GeV) CR injection spectra have a stronger impact on the dynamics than power-law spectra. For realistic magnetic field configurations low energy CRs (with smaller diffusion coefficients) distribute anisotropically with large spatial variations of a factor of ten and more. Adiabatic losses can change the local spectra perceptibly but do not have an integral effect on the dynamics at the spatial and temporal scales considered here. We discuss the potential global impact of CRs and anisotropic transport on the dynamical structure of the ISM and also detail the limitations of the model.

Most stars form in a clustered environment. Therefore, it is important to assess how this environment influences the evolution of protoplanetary discs around young stars. In turn, this affects their ability to produce planets and ultimately life. We present here for the first time 3D SPH/N-body simulations that include both the hydrodynamical evolution of the discs around their natal stars, as well as the dynamics of the stars themselves. The discs are viscously evolving, accreting mass onto the central star and spreading. We find penetrating encounters to be very destructive for the discs as in previous studies, although the frequency of such encounters is low. We also find, however, that encounter influence the disc radii more strongly than other disc properties such as the disc mass. The disc sizes are set by the competition between viscous spreading and the disruptive effect of encounters. As discs spread, encounters become more and more important. In the regime of rapid spreading encounters simply truncate the discs, stripping the outer portions. In the opposite regime, we find that the effect of many distant encounters is able to limit the disc size. Finally, we predict from our simulations that disc sizes are limited by encounters at stellar densities exceeding $\sim 2-3 \times 10^3 \ \mathrm{pc}^{-2}$.

We have simulated star formation in prestellar cores, using SPH and initial conditions informed by observations of the cores in Ophiuchus. Because the observations are limited to two spatial dimensions plus radial velocity, we cannot infer initial conditions for the collapse of a particular core. However, with a minimum of assumptions (isotropic turbulence with a power-law spectrum, a thermal mix of compressive and solenoidal modes, a critical Bonnor-Ebert density profile) we can generate initial conditions that match, in a statistical sense, the distributions of mass, projected size and aspect ratio, thermal and non-thermal one-dimensional velocity dispersion, observed in Ophiuchus. The time between core-core collisions in Ophiuchus is sufficiently long, that we can simulate single cores evolving is isolation, and therefore we are able to resolve masses well below the opacity limit. We generate an ensemble of 100 cores, and evolve them with no radiative feedback from the stars formed, then with continuous radiative feedback, and finally with episodic radiative feedback. With no feedback the simulations produce too many brown dwarfs, and with continuous feedback too few. With episodic radiative feedback, both the peak of the protostellar mass function (at ~ 0.2 M_sun) and the ratio of H-burning stars to brown dwarfs are consistent with observations. The mass of a star is not strongly related to the mass of the core in which it forms. Low-mass cores (M ~ 0.1 M_sun) tend to collapse into single objects, whereas high-mass cores (1 > M_sun) usually fragment into several objects with a broad mass range.

The spatial segregation between dSphs and dIrrs in the Local Group has long been regarded as evidence of an interaction with their host galaxies. In this paper, we assume that ram-pressure stripping is the dominant mechanism that removed gas from the dSphs and we use this to derive a lower bound on the density of the corona of the Milky Way at large distances (50-90 kpc) from the Galactic centre. At the same time, we derive an upper bound by demanding that the interstellar medium of the dSphs is in pressure equilibrium with the hot corona. We consider two dwarfs (Sextans and Carina) with well-determined orbits and star formation histories. Our approach introduces several novel features: we use the measured star formation histories of the dwarfs to derive the time at which they last lost their gas, and (via a modified version of the Kennicutt-Schmidt relation) their internal gas density at that time; we use a large suite of 2D hydrodynamical simulations to model the gas stripping; and we include supernova feedback tied to the gas content. Despite having very different orbits and star formation histories, we find results for the two dSphs that are in excellent agreement with one another. We derive an average particle density of the corona of the Milky Way at 50-90 kpc in the range 1.3-3.6 10^{-4} cm^{-3}. Including additional constraints from X-ray emission limits and pulsar dispersion measurements, we extrapolate Galactic coronal density profiles and we estimate the fraction of baryons that can exist within the virial radius of the Milky Way. For an isothermal corona (T=1.8 10^6 K) this is small, 10-20 % of the universal baryon fraction. Only a hot (T=3 10^6 K) and adiabatic corona can contain all of the Galaxy's missing baryons. Models for the Milky Way must explain why its corona is in a hot adiabatic thermal state or why a large fraction of its baryons lie beyond the virial radius.

Infrared shells and bubbles are ubiquitous in the Galaxy and can generally be associated with HII regions formed around young, massive stars. In this paper, we use high-resolution 3D SPH simulations to explore the effect of a single O7 star emitting photons at 10^49 1/s and located at the centre of a molecular cloud with mass 10^4 M_sun and radius 6.4 pc; the internal structure of the cloud is characterised by its fractal dimension, D (with 2.0 <= D <= 2.8), and the variance of its (log-normal) density distribution, sigma_0^2 (with 0.36 <= sigma_0^2 <= 1.42). Our study focuses on the morphology of the swept-up cold gas and the distribution and statistics of the resulting star formation. If the fractal dimension is low, the border of the HII region is dominated by extended shell-like structures, and these break up into a small number of massive high-density clumps which then spawn star clusters; star formation occurs relatively quickly, and delivers somewhat higher stellar masses. Conversely, if the fractal dimension is high, the border of the HII region is dominated by a large number of pillars and cometary globules, which contain compact dense clumps and tend to spawn single stars or individual multiple systems; star formation occurs later, the stellar masses are somewhat lower, and the stars are more widely distributed.

Observations indicate that the central portions of the Present-Day Prestellar Core Mass Function (CMF) and the Stellar Initial Mass Function (IMF) both have approximately log-normal shapes, but that the CMF is displaced to higher mass than the IMF by a factor F = 4+/-1. This has lead to suggestions that the shape of the IMF is directly inherited from the shape of the CMF - and therefore, by implication, that there is a self-similar mapping from the CMF onto the IMF. If we assume a self-similar mapping, it follows (i) that F = N0/eta, where eta is the mean fraction of a core's mass that ends up in stars, and N0 is the mean number of stars spawned by a single core; and (ii) that the stars spawned by a single core must have an approximately log-normal distribution of relative masses, with universal standard deviation sigma0. Observations can be expected to deliver ever more accurate estimates of F, but this still leaves a degeneracy between eta and N0; and sigma0 is also unconstrained by observation. Here we show that these parameters can be estimated by invoking binary statistics. Specifically, if (a) each core spawns one long-lived binary system, and (b) the probability that a star of mass M is part of this long-lived binary is proportional to M^alpha, current observations of the binary frequency as a function of primary mass and the distribution of mass ratios strongly favour eta = 1.0 +/- 0.3, N0 = 4.3 +/- 0.4, sigma0 = 0.3 +/- 0.03 and alpha = 0.9 +/- 0.6. The mapping from CMF to IMF is not necessarily self-similar - there are many possible motivations for a non self-similar mapping - but if it is not, then the shape of the IMF cannot be inherited from the CMF. Given the limited observational constraints currently available and the ability of a self-similar mapping to satisfy them, the possibility that the shape of the IMF is inherited from the CMF cannot be ruled out at this juncture.

Numerical simulations of star formation frequently rely on the implementation of sink particles, (a) to avoid expending computational resource on the detailed internal physics of individual collapsing protostars, (b) to derive mass functions, binary statistics and clustering kinematics (and hence to make comparisons with observation), and (c) to model radiative and mechanical feedback; sink particles are also used in other contexts, for example to represent accreting black holes in galactic nuclei. We present a new algorithm for creating and evolving sink particles in SPH simulations, which appears to represent a significant improvement over existing algorithms {\refrpt -- particularly in situations where sinks are introduced after the gas has become optically thick to its own cooling radiation and started to heat up by adiabatic compression}. (i) It avoids spurious creation of sinks. (ii) It regulates the accretion of matter onto a sink so as to mitigate non-physical perturbations in the vicinity of the sink. (iii) Sinks accrete matter, but the associated angular momentum is transferred back to the surrounding medium. With the new algorithm -- and modulo the need to invoke sufficient resolution to capture the physics preceding sink formation -- the properties of sinks formed in simulations are essentially independent of the user-defined parameters of sink creation, or the number of SPH particles used.

The role of feedback from massive stars is believed to be a key element in the evolution of molecular clouds. We use high-resolution 3D SPH simulations to explore the dynamical effects of a single O7 star located at the centre of a molecular cloud with mass 10^4M_sun and radius 6.4pc. The initial internal structure of the cloud is characterised by its fractal dimension, D=2.0 - 2.8, and its log-normal density PDF. (i) As regards star formation, in the short term ionising feedback is positive, in the sense that star formation occurs much more quickly in gas that is compressed by the high pressure of the ionised gas. However, in the long term ionising feedback is negative, in the sense that most of the cloud is dispersed with an outflow rate of up to ~0.01M_sun/yr, on a timescale comparable with the sound-crossing time for the ionised gas (~1-2Myr), and triggered star formation is therefore limited to a few percent of the cloud's mass. (ii) As regards the morphology of the ionisation fronts (IFs) bounding the HII region and the systematics of outflowing gas, we distinguish two regimes. For low D<=2.2, the initial cloud is dominated by large-scale structures, so the neutral gas tends to be swept up into a few extended coherent shells, and the ionised gas blows out through a few large holes between these shells; we term these HII regions "shell-dominated". Conversely, for high D>=2.6, the initial cloud is dominated by small-scale structures, and these are quickly overrun by the advancing IF, thereby producing neutral pillars whilst the ionised gas blows out through a large number of small holes between the pillars; we term these HII regions "pillar-dominated". (iii) As regards the injection of bulk kinetic energy, by ~1Myr, the expansion of the HII region has delivered a rms velocity of ~6km/s; this represents less than 0.1% of the total energy radiated by the O7 star.

The collapse of weakly turbulent prestellar cores is a critical stage in the process of star formation. Being highly non-linear and stochastic, the outcome of collapse can only be explored theoretically by performing large ensembles of numerical simulations. Standard practice is to quantify the initial turbulent velocity field in a core in terms of the amount of turbulent energy (or some equivalent) and the exponent in the power spectrum (n \equiv -d log Pk /d log k). In this paper, we present a numerical study of the influence of the details of the turbulent velocity field on the collapse of an isolated, weakly turbulent, low-mass prestellar core. We show that, as long as n > 3 (as is usually assumed), a more critical parameter than n is the maximum wavelength in the turbulent velocity field, {\lambda}_MAX. This is because {\lambda}_MAX carries most of the turbulent energy, and thereby influences both the amount and the spatial coherence of the angular momentum in the core. We show that the formation of dense filaments during collapse depends critically on {\lambda}_MAX, and we explain this finding using a force balance analysis. We also show that the core only has a high probability of fragmenting if {\lambda}_MAX > 0.5 R_CORE (where R_CORE is the core radius); that the dominant mode of fragmentation involves the formation and break-up of filaments; and that, although small protostellar disks (with radius R_DISK <= 20 AU) form routinely, more extended disks are rare. In turbulent, low-mass cores of the type we simulate here, the formation of large, fragmenting protostellar disks is suppressed by early fragmentation in the filaments.

We present simulations of initially stable isothermal clouds exposed to ionising radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN and an HEALPix-based photo-ionisation algorithm to simulate the propagation of the ionising radiation and the resulting dynamical evolution of the cloud. We find that the incident ionising flux, $\Phi_{_{\rm LyC}}$, is the critical parameter determining the cloud evolution. At moderate fluxes, a large fraction of the cloud mass is converted into stars. As the flux is increased, the fraction of the cloud mass that is converted into stars and the mean masses of the individual stars both decrease. Very high fluxes simply disperse the cloud. Newly-formed stars tend to be concentrated along the central axis of the cloud (i.e. the axis pointing in the direction of the incident flux). For given cloud parameters, the time, $t_{_\star}$, at which star formation starts is proportional to $\Phi_{_{\rm LyC}}^{-1/3}$. The pattern of star formation found in the simulations is similar to that observed in bright-rimmed clouds.

We study the influence of gas metallicity, turbulence, and non-equilibrium chemistry on the evolution of the two-phase interstellar medium (warm and cold atomic phases), and thereby constrain the initial conditions for star formation prevailing in turbulent gas. We perform high-resolution simulations in three dimensions, including a realistic non-equilibrium treatment of the ionization state of the gas, and examine both driven and decaying turbulence. This allows us to explore variations in the metallicity Z. In this paper, we study solar metallicity, Z=Z_sun, and low metallicity, Z=1d-3 Z_sun, gas. For driven, large-scale turbulence, we find that the influence of the metallicity on the amount of mass in the cold gas component is small. However, in decaying turbulent conditions this picture is much changed. While cold regions survive in the case of solar metallicity, they are quickly heated and dispersed in low-metallicity gas. This result suggests that star formation can be suppressed in environments of low metallicity, unless a strong turbulent driver is acting on time scales shorter than a few turbulent crossing times. Inter alia this finding could explain the overall inefficient star formation as well as the burst-like mode of star formation found in metal-poor, gas-rich systems like dwarf galaxies.

The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of dynamical scales involved, following this process with numerical simulations is challenging. We compare two-dimensional simulations of flow-driven cloud formation without self-gravity, using the Lagrangian Smoothed Particle Hydrodynamics (SPH) code VINE and the Eulerian grid code Proteus. Results are qualitatively similar for both methods, yet the variable spatial resolution of the SPH method leads to smaller fragments and thinner filaments, rendering the overall morphologies different. Thermal and hydro-dynamical instabilities lead to rapid cooling and fragmentation into cold clumps with temperatures below 300K. For clumps more massive than 1 Msun/pc, the clump mass function has an average slope of -0.8. The internal velocity dispersion of the clumps is nearly an order of magnitude smaller than their relative motion, rendering it subsonic with respect to the internal sound speed of the clumps, but supersonic as seen by an external observer. For the SPH simulations most of the cold gas resides at temperatures below 100K, while the grid-based models show an additional, substantial component between 100 and 300K. Independently of the numerical method our models confirm that converging flows of warm neutral gas fragment rapidly and form high-density, low-temperature clumps as possible seeds for star formation.

We study the structural evolution of turbulent molecular clouds under the influence of ionizing radiation emitted from a nearby massive star by performing a high resolution parameter study with the iVINE code. The temperature is taken to be 10K or 100K, the mean number density is either 100cm^3 or 300cm^3. Besides, the turbulence is varied between Mach 1.5 and Mach 12.5 and the main driving scale between 1pc and 8pc. We vary the ionizing flux by an order of magnitude. In our simulations the ionizing radiation enhances the initial turbulent density distribution and thus leads to the formation of pillar-like structures observed adjacent to HII regions in a natural way. Gravitational collapse occurs regularly at the tips of the structures. We find a clear correlation between the initial state of the turbulent cold cloud and the final morphology and physical properties of the structures formed. The most favorable regime for the formation of pillars is Mach 4-10. Structures and therefore stars only form if the initial density contrast between the high density unionized gas and the gas that is going to be ionized is lower than the temperature contrast between the hot and the cold gas. The density of the resulting pillars is determined by a pressure equilibrium between the hot and the cold gas. A thorough analysis of the simulations shows that the complex kinematical and geometrical structure of the formed elongated filaments reflects that of observed pillars to an impressive level of detail. In addition, we find that the observed line-of sight velocities allow for a distinct determination of different formation mechanisms. Comparing the current simulations to previous results and recent observations we conclude that e.g. the pillars of creation in M16 formed by the mechanism proposed here and not by the radiation driven implosion of pre-existing clumps.

We present simulations of initially stable isothermal clouds exposed to ionizing radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN (Hubber et al. 2010) and the HEALPix-based photoionization algorithm described in Bisbas et al. (2009). We find that the incident ionizing flux is the critical parameter determining the evolution: high fluxes simply disperse the cloud, whereas low fluxes trigger star formation. We find a clear connection between the intensity of the incident flux and the parameters of star formation.

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.

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)

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent prestellar cores, using SPH simulations. The initial conditions are slightly supercritical Bonnor-Ebert spheres, all with the same density profile, the same mass (M_O=6.1 Msun) and the same radius (R_O=17,000 AU), but having different initial turbulent velocity fields. Four hundred turbulent velocity fields have been generated, all scaled so that the mean Mach number is M=1. Then a subset of these, having a range of net angular momenta, j, has been evolved. The evolution of these turbulent cores is not strongly correlated with j. Instead it is moderated by the formation of filamentary structures due to converging turbulent flows. A high fraction (~ 82%) of the protostars forming from turbulent cores are attended by protostellar accretion discs, but only a very small fraction (~16%) of these discs is sufficiently cool and extended to develop non-linear gravitational instabilities and fragment.

Abridged: We use three-dimensional SPH simulations to investigate the collapse of low-mass prestellar cores and the formation and early evolution of protostellar discs. The initial conditions are slightly supercritical Bonnor-Ebert spheres in rigid rotation. The core mass and initial radius are held fixed at M_O=6.1 M_sun and R_O=17,000 AU, and the only parameter that we vary is the initial angular speed \Omega_O. Protostellar discs forming from cores with \Omega_O<1.35 10d-13 1/s have radii between 100 and 300 AU and are quite centrally concentrated; due to heating by gas infall onto the disc and accretion onto the central object, they are also quite warm, T>100 K, and therefore stable against gravitational fragmentation. In contrast, more rapidly rotating cores form discs which are less concentrated and cooler, and have radii between 400 and 1000 AU; as a consequence they are prone to gravitational fragmentation and the formation of multiple systems. We derive a criterion that predicts whether a rigidly rotating core having given M_O, R_O and \Omega_O will produce a protostellar disc which fragments whilst material is still infalling from the core envelope. We then apply this criterion to core samples for which M_O, R_O and \Omega_O have been estimated observationally. We conclude that the observed cores are stable against fragmentation at this stage, due to their low angular speeds and the heat delivered at the accretion shock where the infalling material hits the disc.

We present high resolution simulations on the impact of ionizing radiation of massive O-stars on the surrounding turbulent interstellar medium (ISM). The simulations are performed with the newly developed software iVINE which combines ionization with smoothed particle hydrodynamics (SPH) and gravitational forces. We show that radiation from hot stars penetrates the ISM, efficiently heats cold low density gas and amplifies over-densities seeded by the initial turbulence. The formation of observed pillar-like structures in star forming regions (e.g. in M16) can be explained by this scenario. At the tip of the pillars gravitational collapse can be induced, eventually leading to the formation of low mass stars. Detailed analysis of the evolution of the turbulent spectra shows that UV-radiation of O-stars indeed provides an excellent mechanism to sustain and even drive turbulence in the parental molecular cloud.

We present a three-dimensional, fully parallelized, efficient implementation of ionizing UV radiation for smoothed particle hydrodynamics (SPH) including self-gravity. Our method is based on the SPH/tree code VINE. We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from e.g. young massive stars on their turbulent parental molecular clouds. In this paper we describe the concept and the numerical implementation of the radiative transfer for a plain-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several $10^5$ years earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of proto-stellar envelopes by ionizing radiation.

We have studied the afterglow of the gamma-ray burst (GRB) of February 18, 2006. This is a nearby long GRB, with a very low peak energy, and is therefore classified as an X-ray Flash (XRF). XRF 060218 is clearly associated with a supernova -- dubbed SN 2006aj. We present early spectra for SN 2006aj as well as optical lightcurves reaching out to 50 days past explosion. Our optical lightcurves define the rise times, the lightcurve shapes and the absolute magnitudes in the U, V and R bands, and we compare these data with data for other relevant supernovae. SN 2006aj evolved quite fast, somewhat similarly to SN 2002ap, but not as fast as SN 1994I. Our spectra show the evolution of the supernova over the peak, when the U-band portion of the spectrum rapidly fades due to extensive line blanketing. We compare to similar spectra of very energetic Type Ic supernovae. Our first spectra are earlier than spectra for any other GRB-SN. The spectrum taken 12 days after burst in the rest frame is similar to somewhat later spectra of both SN 1998bw and SN 2003dh, implying a rapid early evolution. This is consistent with the fast lightcurve. From the narrow emission lines from the host galaxy we derive a redshift of z=0.0331+-0.0007. This makes XRF 060218 the second closest gamma-ray burst detected. The flux of these emission lines indicate a high-excitation state, and a modest metallicity and star formation rate of the host galaxy.

We present a new perspective on the crystallinity of dust in protoplanetary disks. The dominant crystallization by thermal annealing happens in the very early phases of disk formation and evolution. Both the disk properties and the level of crystallinity are thereby directly linked to the properties of the molecular cloud core from which the star+disk system was formed. We show that, under the assumption of single star formation, rapidly rotating clouds produce disks which, after the main infall phase (i.e. in the optically revealed class II phase), are rather massive and have a high accretion rate but low crystallinity. Slowly rotating clouds, on the other hand, produce less massive disks with lower accretion rate, but high levels of crystallinity. Cloud fragmentation and the formation of multiple stars complicates the problem and necessitates further study. The underlying physics of the model is insufficiently understood to provide the precise relationship between crystallinity, disk mass and accretion rate. But the fact that with `standard' input physics the model produces disks which, in comparison to observations, appear to have either too high levels of crystallinity or too high disk masses, demonstrates that the comparison of these models to observations can place strong contraints on the disk physics. The question to ask is not why some sources are so crystalline, but why some other sources have such a low level of crystallinity.

UVOIR bolometric light curves provide valuable insight into the nature of type Ia supernovae. We present an analysis of sixteen well-observed SNe Ia. Constraints are placed on several global parameters concerning the progenitor system, explosion mechanism and subsequent radiation transport. By fitting a radioactive decay energy deposition function to the quasi-exponential phase (50 to 100 days after maximum light), it is found that the ejected mass varies by at least a factor of two. This result suggests that a sub-Chandrasekhar mass model could be responsible for the progenitor system of some type Ia supernovae. We find that the range in the amount of synthesized (56)Ni indicates a significant variation in the burning mechanism. In order to explain a factor of ten range in the observed bolometric luminosity more detailed modeling of the explosion mechanism is required.

We discuss spectral energy distributions of a sample of Herbig Ae/Be stars in the context of a passive irradiated disk model. The data have been presented earlier by Meeus et al (2001), and preliminary interpretations of these data were given in that paper. While the spectra of Herbig Ae stars all show similarities, there is significant variation between the spectra, in particular in the shape of the mid-IR rise and in the presence or absence of a silicate feature. We explore the hypothesis that all these different spectra can be interpreted as pure disk spectra without additional components. Using the model of Dullemond, Dominik and Natta (2001) we deduce the disk parameters of a number of the sources, and find that for a large fraction of investigated sources, satisfactory fits can be obtained. The derived model parameters show that some group Ia sources can only be fit with radially increasing surface densities, indicating the presence of depleted inner disk regions. The steep-sloped SEDs of group IIa sources can be fit with very compact disks, probably representing disks with collapsed outer regions. The largest difficulties arise from sources that do not show significant silicate emission features. Our attempts to explain these objects with a pure geometric effect are only partially successful. It seems that these stars indeed require a strong depletion of small silicate grains.