Search

We present magnetohydrodynamic simulations of star formation in the multiphase interstellar medium to quantify the impact of non-ionising far-ultraviolet (FUV) radiation. This study is carried out within the framework of the \textsc{Silcc Project}. It incorporates the radiative transfer of ionising radiation and self-consistent modelling of variable FUV radiation from star clusters. Near young star clusters, the interstellar radiation field (ISRF) can reach values of $G_0 \approx 10^4$ (in Habing units), far exceeding the canonical solar neighbourhood value of $G_0 = 1.7$. However, our findings suggest that FUV radiation has minimal impact on the integrated star formation rate compared to other feedback mechanisms such as ionising radiation, stellar winds, and supernovae. Only a slight decrease in star formation burstiness, related to increased photoelectric heating efficiency by the variable FUV radiation field, is detectable. Dust near star-forming regions can be heated up to 60 K via the photoelectric (PE) effect, showing a broad temperature distribution. PE heating rates for variable FUV radiation models show higher peak intensities but lower average heating rates than static ISRF models. Simulations of solar neighbourhood conditions without stellar winds or ionising radiation but with self-consistent ISRF and supernovae show high star formation rates $\sim10^{-1}\,\mathrm{M_\odot\,yr^{-1}\,kpc^{-2}}$, contradicting expectations. Our chemical analysis reveals increased cold neutral medium volume-filling factors (VFF) outside the vicinity of stellar clusters with a variable ISRF. Simultaneously, the thermally unstable gas is reduced, and a sharper separation of warm and cold gas phases is observed. The variable FUV field also promotes a diffuse molecular gas phase with VFF of $\sim5-10$~per cent., Comment: submitted to MNRAS

We present an analysis of the cold gas phase in a low metallicity starburst generated in a high-resolution hydrodynamical simulation of a gas-rich dwarf galaxy merger as part of the GRIFFIN project. The simulations resolve (4 M$_\odot$ gas phase mass resolution, $\sim$ 0.1 pc spatial resolution) the multi-phase interstellar medium with a non-equilibrium chemical heating/cooling network at temperatures below $10^4$ K. Massive stars are sampled individually and interact with the ISM through the formation of HII regions and supernova explosions. In the extended starburst phase, the ISM is dominated by cold ($T_\mathrm{gas} < 300$ K) filamentary clouds with self-similar internal structures. The clouds have masses of $10^{2.6}$ - $10^{5.6}$ M$_\odot$ with a power law mass function, $dN/dM \propto M^\alpha$ with $\alpha = -1.78 (\pm 0.08)$. They also follow the Larson relations, in good agreement with observations. We trace the lifecycle of the cold clouds and find that they follow an exponential lifetime distribution and an e-folding time of $\sim$ 3.5 Myr. Clouds with peak masses below $10^4$ M$_\odot$ follow a power law relation with their average lifetime $\tau_\mathrm{life} \propto M^{0.3}_\mathrm{max}$ which flattens out for higher cloud masses at $ < 10$ Myr. A similar relation exists between cloud size at peak mass and lifetime. This simulation of the evolution of a realistic galactic cold cloud population supports the rapid formation and disruption of star-forming clouds by stellar radiation and supernovae on a timescale less than 10 Myr., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society Main Journal

Supernovae (SNe) inject $\sim 10^{51}$ erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the FLASH code to calculate the cooling radiation from shock-heated gas using collisional excitation data from MAPPINGS V. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ($[$O III$]$, $[$N II$]$, $[$S II$]$, H${\alpha}$, H${\beta}$) are usually best observable. Our shock detection scheme shows that [S II] and [N II] emissions arise from the thin shell surrounding the SNR, while [O III], H$\rm \alpha$, and H$\rm \beta$ originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR's complex structure and its projection onto the plane of the sky because the escaping line luminosity can be reduced by 10 -- 80\% due to absorption along the line-of-sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity., Comment: 18 pages, 13 figures, accepted for publication in MNRAS

We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities $\Sigma_\mathrm{gas} =$ 10, 30, 50, and 100 $\mathrm{M_\odot\,pc^{-2}}$. The SILCC project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionising UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between $\Sigma_\mathrm{gas}$ and the star formation rate surface density $\Sigma_\mathrm{SFR}$. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher $\Sigma_\mathrm{gas}$ and the mass loading factors of the CR-supported outflows are of order unity independent of $\Sigma_\mathrm{SFR}$. Similar to observations, vertical velocity dispersions of the warm ionised medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma_\mathrm{z} \propto \Sigma_\mathrm{SFR}^a$, with $a \sim 0.20$. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor $\sim 2.2$ higher than that of the CNM, in agreement with local observations. For $\Sigma_\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$ the WIM motions become supersonic., Comment: 20 pages, 12 figures, accepted by MNRAS

We present simulations of the multi-phase interstellar medium (ISM) at solar neighbourhood conditions including thermal and non-thermal ISM processes, star cluster formation, and feedback from massive stars: stellar winds, hydrogen ionising radiation computed with the novel TreeRay radiative transfer method, supernovae (SN), and the injection of cosmic rays (CR). N-body dynamics is computed with a 4th-order Hermite integrator. We systematically investigate the impact of stellar feedback on the self-gravitating ISM with magnetic fields, CR advection and diffusion and non-equilibrium chemical evolution. SN-only feedback results in strongly clustered star formation with very high star cluster masses, a bi-modal distribution of the ambient SN densities, and low volume-filling factors (VFF) of warm gas, typically inconsistent with local conditions. Early radiative feedback prevents an initial starburst, reduces star cluster masses and outflow rates. Furthermore, star formation rate surface densities of $\Sigma_{\dot{M}_\star} = 1.4-5.9 \times 10^{-3}$ $\mathrm{M}_\odot\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$, VFF$_\mathrm{warm} = 60-80$ per cent as well as thermal, kinetic, magnetic, and cosmic ray energy densities of the model including all feedback mechanisms agree well with observational constraints. On the short, 100 Myr, timescales investigated here, CRs only have a moderate impact on star formation and the multi-phase gas structure and result in cooler outflows, if present. Our models indicate that at low gas surface densities SN-only feedback only captures some characteristics of the star-forming ISM and outflows/inflows relevant for regulating star formation. Instead, star formation is regulated on star cluster scales by radiation and winds from massive stars in clusters, whose peak masses agree with solar neighbourhood estimates., Comment: 23 pages, 16 figures, accepted by MNRAS, corrected an incorrect reference