zum Inhalt springen

Multi-band ray-tracing of ionizing photons

1 / 2
  • Fig. 1) Example SED for a star with surface temperature 50,000 K. Shown is the simple black body spectrum for T=50,000 K as well as a more realistic SEDs. The blue and red lines indicate how the ionizing part of the spectrum can be divided into two or more energy bands. Each energy band has its own mean photon energy and interacts with the gas slightly differently.
  • Fig. 2) The radial temperature profile of an HII region in a uniform medium simulated with the new multi-band ray-tracing method. The temperature rises towards larger radii, because the H and He ionization cross sections decrease with the photon energy, ie. photons with lower energies are absorbed preferably. The inner peak marks the edge of the He II region. The outer peak is at the edge of the HII region.

Author: Sebastian Vider

In the previous monthly highlight, an overview was provided on the subject of simulating the processes of the Interstellar Medium (ISM). In this segment, we will take a closer look at one of those processes, namely ionizing radiation, and how to simulate it. Ionizing radiation is defined as the radiation of photons with an energy greater than 13.6 eV. These photons possess the capacity to ionize hydrogen. Massive stars, due to their high surface temperatures, emit a substantial number of ionizing photons, leading to the formation of HII regions, which are regions of ionized hydrogen gas. Stars with surface temperatures above approximately 40,000 K emit photons with even higher energies, namely greater than 24.6 eV, and can therefore also ionize helium and create HeII regions. The temperature structure within the ionized regions is determined by a balance between heating and cooling processes. Heating is provided primarily by the ionizing photons, which transfer their excess energy, that is, the energy of a photon exceeding the 13.6 eV required to ionize H, to the gas. The predominant cooling processes are cooling due to metal lines and recombination cooling.

In order to model the radiative transport of the ionizing photons with the 3D hydrodynamics code Flash, we employ a ray-tracing algorithm. The standard version of the algorithm treats all ionizing photons equally and assumes that they all have the same mean energy, calculated from the stellar spectral energy distribution (SED) of the radiation source. The consequence is that the HII regions simulated with that method have a flat temperature profile. Furthermore, the method cannot produce the correct abundances of ionized elements, with the exception of hydrogen.

The new multi-band UV ray-tracing module allows the user to define up to 9 energy bands in the range between 13.6 eV to 100 eV (see Fig. 1). This makes it possible to accurately simulate the ionization states of several species such as He, O, C, Si, N, etc. and results in a more realistic temperature profile (see Fig. 2). This has a small impact on the dynamics of the region, but more importantly, it is now possible to generate much more accurate synthetic emission maps from ionized species such as [CII].