Spectrally resolved cosmic rays

Cosmic rays (CRs) are an important component in the interstellar medium (ISM), which can heat and ionize the gas in dense star forming cores, alter the diffuse ISM and drive glactic outflows. The energy of CRs range from just above the thermal distribution in the non-relativistic regime up to the ultra-relativistic part with momenta up to 1020 eV. The details of the impact of CRs, they interaction with the gas, the cooling rates and transport properties strongly depend on the momentum of the CRs. Low energy CRs ionize the gas and can cool efficiently via Coulomb and ionization losses. GeV CRs contain most of the CR energy, which is comparable to the other energy components in the ISM and can have a direct dynamical impact. Ultra-relativistic CRs do not dominate in total energy, but provide important observational links via e.g. gamma rays.

In numerical simulations we would like to accurately incorporate CRs from the non-relativistic up to the ultra-relativistic limit. Ideally we would like to include a full spectrum in every computational, which requires to cover the large range in CR energies and the large dynamic range in the spectrum with efficient numerical methods. I have been working an efficient model that allows to compute the main proton processes including Coulomb & ionization losses, adiabatic processes and hadronic losses. The numerical model and an application to galaxies is published here:

 

Windtunnels

There is strong observational evidence that the gas flowing out of galaxies contains not only hot dilute gas but also cold molecular gas. This fact poses the question of whether the molecular can form within the outflo, i.e. above the galacitc disc, or whether it forms in the disc and is then transported to large altitudes above the disc via the hot gas or hot winds coming from star forming regions. We inverstigate this problem with idealised simulations, in which we expose a warm, dense cloud to a hot, low density wind. We follow how the wind can accelerate, destroy or compress the cloud. There are a number of paramters involed and we mainly focus on the density of the cloud and the magnetic field configuration in the wind. The movie below shows the comparison between three different wind configurations, namely a non-magnetic wind, a magnetic field that is orianted parallel to the direction of the wind and a field perpendicular to the wind

 

We find that magnetic fields perpendicular to the wind velocity are able to best accelerate the cloud, keep the cloud compact and thereby allow the formation of a coherent dense molecular phase. We find that it might be possible to form molecular gas inside the outflows of galaxies.

A second important fact is that self-gravity plays an important role, even if the initial conditions do not suggest that self-gravity matters. The continous compression of gas via the wind combined with radiative cooling of the gas gradually increases the density, decreases the Jeans mass and thus forms regions that are dominated by gravitational effects.

Here are more videos showing the different cloud densities and the chemical composition:

no. name length (sec.) mp4
1 low density, no self-gravity, different magnetic fields 20 mp4 (2M)
2 low density, self-gravity, different magnetic fields 20 mp4 (2M)
3 medium density, no self-gravity, different magnetic fields 20 mp4 (2.8M)
4 medium density, self-gravity, different magnetic fields 15 mp4 (1.9M)
5 high density, no self-gravity, different magnetic fields 10 mp4 (1.8M)
6 high density, self-gravity, different magnetic fields 20 mp4 (2.8M)

 

no. wind / gravity low density medium density high density
1 non-magnetic wind, no self-gravity mp4 mp4 mp4
2 non-magnetic wind, self-gravity mp4 mp4 mp4
3 wind with B parallel, no self-gravity mp4 mp4 mp4
4 wind with B parallel, self-gravity mp4 mp4 mp4
5 wind with B perpendicular, no self-gravity mp4 mp4 mp4
6 wind with B perpendicular, self-gravity mp4 mp4 mp4

 

SILCC with cosmic rays

Cosmic rays (CRs) are an important energy component in the interstellar medium. Being high-energy particles their coupling to the gas is not via direct particle-particle interactions as in the case for thermal gas because their cross section is very small. Instead CRs couple via the magnetic field. By gyrating around magnetic field lines and scattering off of magnetic irregularities, CRs transfer energy and momentum to the magntic field. Together with the approximation of ideal MHD, which freezes the magnetic field lines to the gas, CRs can effectively transfer energy and momentum to the gas and accelerate it. The interaction with the magnetic field allows them to move relative to the gas and effectively diffuse or stream through the ISM. By this mechanism the CR energy can be deposited far away from their production sites. The small effective cross sections also result in low cooling efficiencies for CRs with a momentum of a few GeV/c, which carry most of the CR energy in the ISM. As a result, CRs can prodive a long-lived energy reservoir with a vertical stratification in the galactic disk that can lauch galactic outflows.
We included CRs as a relativistic fluid into the SILCC setup in the diffusion-advection approximation. We investigate the impact of CRs in driving outflows in two papers:

We find that CRs are able to drive stable warm and smooth outflows with mass loading factors (ratio of outflow to star formation rate) of order unity. The simulation data (after acceptance of the papers) is publicly available at the SILCC data website, hosted at the Max-Planck Computing and Data Facility in Garching (MPCDF). An expample video for CR-driven outflows is shown in here:



or as mov: mov (7MB)

 

Movies from Girichidis et al. 2018a (column density)

 

large vertical extent

fps length small 590 px 800 px 1200 px
15 100 s mp4 mp4 mp4 mp4
25 60 s mp4 mp4 mp4 mp4
35 43 s mp4 mp4 mp4 mp4
50 30 s mp4 mp4 mp4 mp4

upper part

fps length small 590 px 800 px 1200 px
15 100 s mp4 mp4 mp4 mp4
25 60 s mp4 mp4 mp4 mp4
35 43 s mp4 mp4 mp4 mp4
50 30 s mp4 mp4 mp4 mp4

 

lower part

fps length small 590 px 800 px 1200 px
15 100 s mp4 mp4 mp4 mp4
25 60 s mp4 mp4 mp4 mp4
35 43 s mp4 mp4 mp4 mp4
50 30 s mp4 mp4 mp4 mp4

 

Movies from Girichidis et al. 2016 (density, density, CR energy)

 

name time (sec.) avi (large) mov (large) mp4 (large)
25fps 55 avi, 1200 (69M) mov, 1200 (26M) mp4, 1200 (23M)
35fps 39 avi, 1200 (50M) mov, 1200 (21M) mp4, 1200 (19M)
50fps 28 avi, 1200 (35M) mov, 1200 (18M) mp4, 1200 (16M)

 

Movies from Girichidis et al. 2016 (density, density)

 

name time (sec.) mov (large) mp4 (large)
800px, 25fps 50 mov, 0800 (6M) mp4, 0800 (7M)
1024px, 25fps 50 mov, 1024 (8M) mp4, 1024 (9M)
1200px, 25fps 50 mov, 1200 (11M) mp4, 1200 (12M)
1280px, 25fps 50 mov, 1280 (11M) mp4, 1280 (12M)

Magnetic fields in the interstellar medium

Magnetic fields are everywhere in the interstellar medium and are locally strong enough to modify the dynamics. Diffuse warm gas in the ISM is very likely to be supported by magnetic pressure. In dense, star forming regions the magnetic fields can be very strong such that the gas flow is channeld by the field structure. Using three-dimensional magneto-hydrodynamical simulations we study the impact of magnetic fields on the ISM, in particular the transition from diffuse warm atomic to dense colde molecular gas. In addition, we investigate the oriantation of the magnetic field compared to the orientation of the gas structures and the gas flows.

I published two papers related to the magnetic fields in the ISM (the first is linked the second will be there soon):

The following video shows the total column density of a fraction of the galactic disc egde-on (top) and face-on (bottom) for three different magnetic fields; no field (left-hand panels), a 3μG field (middle panel), and a 6μG field (right-hand panel).

other video formats

name mov (15 fps) mp4 (15 fps) mov (25 fps) mp4 (25 fps)
column density (800 px), ref mov mp4 mov mp4
column density (1200 px), ref mov mp4 mov mp4
column density (1280 px), ref mov mp4 mov mp4

 

There are two main differences between the simulations. The first is that the runs including magnetic fields have smaller density contrasts. The dense gas is embedded in larger cloud complexes rather than in locally isolated clouds. The second is that dense gas forms later because the additional magnetic pressure supports the gas against fast compression and the resulting faster cooling. This can be nicely seen in the formation of molecular gas here:

and some other video formats

name mov (15 fps) mp4 (15 fps) mov (25 fps) mp4 (25 fps)
H2 column density (800 px) mov mp4 mov mp4
H2 column density (1200 px) mov mp4 mov mp4
H2 column density (1280 px) mov mp4 mov mp4

My research focusses on the theoretical and numerical studies of the multi-phase interstallar medium (ISM). The computations range from the onset and composition of galactic outflows on kpc scale over the chemical evolution in the disc and the formation of molecular clouds down to cold and self-gravitating regions on sub-parsec scales. Besides a chemical evolution and megnetic fields I focus on the dynamcial impact of cosmic rays in the ISM. My projects can be split into different groups.

Cosmic rays in the interstellar medium

Cosmic rays (CRs) are charged high energy particles that have comparable energy densities in the ISM as the magnetic and kinetic one. Due to their relaltivistic nature CRs have relatively complicated transport and coupling properties. I use numerical simulations with two fluid components to investigate how CRs change the structures in the interstellar medium and drive outflows. The video shows MHD-CR simulations of the supernova-driven ISM with different energy injection types. The left panel shows a simulation with purely thermal energy injection for the SNe. The middle panel uses only CR energy injection and the right panel the combined effect of thermal and CR energy injection.



other video formats avi (49MB), mov (2.5MB)

SILCC: SImulating the LifeCycle of molecular Clouds

Within the SILCC project we numerically investigate the formation process of molecular clouds, their chemcial and magnetic structure as well as their destruction in the SN-driven interstellar medium. We use hydrodynmical simulations of stratified boxes for our studies including a chemical network that follows the abundances of ionised, atomic and molecular hydrogen. The video shows the cuts through the centre of the box for the density and the temperature (left two panels) as well as projections of the density (total, ionised hydrogen, atomic hydrogen, molecular hydrogen and CO)



other video formats avi (44MB), mov (8MB)

Turbulent self-gravitating gas

The coldest and densest condensations in the ISM are star-forming regions which mark the border line between molecular clouds that are supported by thermal pressure and turbulent motions and strongly self-gravitating clumps.

Statistical descriptions of turbulence

The complexity of turbulent motions requires simplifications in the models to understand the origin of statistical properties like the energy transfer between different scales and the composition of turbulent modes. I therefore also study simulations of isothermal driven turbulence in periodic boxes. The video below shows the projected density structures of compressively driven turbulence on the left and solinoidally driven turbulence on the right for a supersonic flow with Mach number 8.



other video formats avi (13MB), mov (1.4MB)