Search

We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σ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-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density Σ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 Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized 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 ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼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. [ABSTRACT FROM AUTHOR]

We analyse observational signatures of magnetic fields for simulations of a Milky Way-like disc with supernova-driven interstellar turbulence and self-consistent chemical processes. In particular, we post-process two simulations data sets of the SILCC Project for two initial amplitudes of the magnetic field, B 0 = 3 and 6 μG, to study the evolution of Faraday rotation measures (RM) and synchrotron luminosity. For calculating the RM, three different models of the electron density n e are considered. A constant electron density, and two estimations based on the density of ionized species and the fraction of the total gas, respectively. Our results show that the RM profiles are extremely sensitive to the n e models, which assesses the importance of accurate electron distribution observations/estimations for the magnetic fields to be probed using Faraday RMs. As a second observable of the magnetic field, we estimate the synchrotron luminosity in the simulations using a semi-analytical cosmic ray model. We find that the synchrotron luminosity decreases over time, which is connected to the decay of magnetic energy in the simulations. The ratios between the magnetic, the cosmic ray, and the thermal energy density indicate that the assumption of equipartition does not hold for most regions of the ISM. In particular, for the ratio of the cosmic ray to the magnetic field energy the assumption of equipartition could lead to a wrong interpretation of the observed synchrotron emission. [ABSTRACT FROM AUTHOR]

Molecular outflows contributing to the matter cycle of star-forming galaxies are now observed in small and large systems at low and high redshift. Their physical origin is still unclear. In most theoretical studies, only warm ionized/neutral and hot gas outflowing from the interstellar medium is generated by star formation. We investigate an in situ H2 formation scenario in the outflow using high-resolution simulations, including non-equilibrium chemistry and self-gravity, of turbulent, warm, and atomic clouds with densities 0.1, 0.5, and |$1\, \mathrm{cm}^{-3}$| exposed to a magnetized hot wind. For cloud densities |$\gtrsim 0.5\, \mathrm{cm}^{-3}$|⁠ , a magnetized wind triggers H2 formation before cloud dispersal. Up to 3 per cent of the initial cloud mass can become molecular on |$\sim \! 10\, \mathrm{Myr}$| time-scales. The effect is stronger for winds with perpendicular B -fields and intermediate density clouds (⁠|$n_\mathrm{c}\sim 0.5\, \mathrm{cm}^{-3}$|⁠). Here, H2 formation can be boosted by up to one order of magnitude compared to isolated cooling clouds independent of self-gravity. Self-gravity preserves the densest clouds well past their |$\sim \! 15\, \mathrm{Myr}$| cloud crushing time-scales. This model could provide a plausible in situ origin for the observed molecular gas. All simulations form warm ionized gas, which represents an important observable phase. The amount of warm ionized gas is almost independent of the cloud density but solely depends on the magnetic field configuration in the wind. For low-density clouds (⁠|$0.1\, \mathrm{cm}^{-3}$|⁠), up to 60 per cent of the initially atomic cloud mass can become warm and ionized. [ABSTRACT FROM AUTHOR]

We investigate the impact of cosmic rays (CRs) on galactic outflows from a multiphase interstellar medium with solar neighbourhood conditions. The three-dimensional magnetohydrodynamical simulations include CRs as a relativistic fluid in the advection-diffusion approximation. The thermal and chemical state of the interstellar medium is computed with a non-equilibrium chemical network. We find that CRs [injected with 10  per cent of the supernova (SN) 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 per cent − 25 per cent 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 \sim 1\, \mathrm{kpc,}$$ there are large volumes in approximate pressure equilibrium between the thermal and the 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. [ABSTRACT FROM AUTHOR]

The large vertical scale heights of the diffuse ionized gas (DIG) in disc galaxies are challenging to model, as hydrodynamical models including only thermal feedback seem to be unable to support gas at these heights. In this paper, we use a three-dimensional Monte Carlo radiation transfer code to post-process disc simulations of the Simulating the Life-Cycle of Molecular Clouds project that include feedback by cosmic rays. We show that the more extended discs in simulations including cosmic ray feedback naturally lead to larger scale heights for the DIG which are more in line with observed scale heights. We also show that including a fiducial cosmic ray heating term in our model can help to increase the temperature as a function of disc scale height, but fails to reproduce observed DIG nitrogen and sulphur forbidden line intensities. We show that, to reproduce these line emissions, we require a heating mechanism that affects gas over a larger density range than is achieved by cosmic ray heating, which can be achieved by fine tuning the total luminosity of ionizing sources to get an appropriate ionizing spectrum as a function of scale height. This result sheds a new light on the relation between forbidden line emissions and temperature profiles for realistic DIG gas distributions. [ABSTRACT FROM AUTHOR]

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionizing radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed by a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionizing 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. Ionizing radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to ρ < 10-25 g cm-3, compared to 10-23g 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 Hα emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 per cent to the total Hα emission. Individual massive stars (M ≥ 30M☉) with short lifetimes are responsible for significant fluctuations in the Hα luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by a factor of ~10. [ABSTRACT FROM AUTHOR]

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more selfconsistent 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 (sHI 0.8s300-8000K ~ 10-20 km s-1, sHa 0.6s8000-3 105 K ~ 20-30 km s-1), and strong outflows with mass loading factors (ratio of outflow to star formation rate) 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 (~ 10-25-10-24 g cm-3) and slow (v ~ 20-40 km s-1) but there is a high-velocity tail of up to v ~ 500 km s-1 with ~ 10-28-10-27 g cm-3. [ABSTRACT FROM AUTHOR]