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, H 2 , 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 −1 ) and are denser (ρ ∼ 10 −24 − 10 −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.

We present simulations of the magnetized interstellar medium (ISM) in models of massive star forming (40 M ⊙ yr −1 ) disk galaxies with high gas surface densities (Σ gas ∼ 100 M ⊙ pc −2 ) similar to observed star forming high-redshift disks. We assume that type II supernovae deposit 10 per cent of their energy into the ISM as cosmic rays and neglect the additional deposition of thermal energy or momentum. With a typical Galactic diffusion coefficient for CRs (3 · 10 28 cm 2 s −1 ) we demonstrate that this process alone can trigger the local formation of a strong low density galactic wind maintaining vertically open field lines. Driven by the additional pressure gradient of the relativistic fluid the wind speed can exceed 10 3 km s −1 , much higher than the escape velocity of the galaxy. The global mass loading, i.e. the ratio of the gas mass leaving the galactic disk in a wind to the star formation rate becomes of order unity once the system has settled into an equilibrium. We conclude that relativistic particles accelerated in supernova remnants alone provide a natural and efficient mechanism to trigger winds similar to observed mass-loaded galactic winds in high-redshift galaxies. These winds also help explaining the low efficiencies for the conversion of gas into stars in galaxies as well as the early enrichment of the intergalactic medium with metals. This mechanism can be at least of similar importance than the traditionally considered momentum feedback from massive stars and thermal and kinetic feedback from supernova explosions.

We report the first results of the first global galactic-scale cosmic ray (CR)-MHD simulations of CR-driven dynamo. We investigate the dynamics of magnetized interstellar medium (ISM), which is dynamically coupled with CR gas. We assume that exploding stars deposit small-scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of CRs, accelerated in supernova shocks. We conduct numerical simulations with the aid of a new parallel MHD code PIERNIK. We find that the initial magnetization of galactic disks by exploding magnetized stars forms favorable conditions for the CR-driven dynamo. We demonstrate that dipolar magnetic fields supplied on small supernova remnant scales can be amplified exponentially by the CR-driven dynamo, to the present equipartition values, and transformed simultaneously to large galactic scales. The resulting magnetic field structure in an evolved galaxy appears spiral in the face-on view and reveals the so-called X-shaped structure in the edge-on view.

We investigate the impact of cosmic rays (CRs) on galactic outflows from a multiphase 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 (ρ ∼ 10 −26 g cm −3 ), colder (∼ 10 4 K), and slower (∼ 30 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 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 ∼ 1 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 ∼ 1 kpc away from the midplane, CRs provide the dominant gas acceleration mechanism.

We present the first numerical model of the magnetohydrodynamical cosmic-ray (CR) driven dynamo of the type proposed by Parker (1992). The driving force of the amplification process comes from CRs injected into the galactic disk in randomly distributed spherical regions representing supernova remnants. The underlying disk is differentially rotating. An explicit resistivity is responsible for the dissipation of the small-scale magnetic field component. We obtain amplification of the large-scale magnetic on a timescale 250 Myr.